(57)【特許請求の範囲】
【請求項１】
複数の受信アンテナで受信した複数の受信信号から、複数の送信アンテナで送信した複数の送信信号を分離して推定する最尤復号化方法であって、
前記受信信号に対応するチャネルインパルス応答から、チャネル行列を生成する第１の手順と、
前記生成されたチャネル行列を複数の小行列に分割し、前記分割された小行列に対応する前記受信信号の部分を特定し、前記分割された小行列の逆行列を用いて前記分割された小行列を変換する第２の手順と、
前記変換された小行列をＱＲ分解することによって三角行列を取得し、前記取得した三角行列を用いて前記各部分の受信信号を変換する第３の手順と、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ決定する第４の手順と、を含むことを特徴とする最尤復号化方法。
【請求項２】
前記第４の手順では、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ以上選択し、
前記各部分の決定された組み合わせ候補間のユークリッド距離が最も短くなるように、前記各部分の組み合わせ候補を一つ決定することを特徴とする請求項１に記載の最尤復号化方法。
【請求項３】
前記第２の手順で用いられる逆行列は、前記分割された小行列の対角項成分に雑音成分を付加して逆行列を求めることによって生成されることを特徴とする請求項１に記載の最尤復号化方法。
【請求項４】
前記第２の手順では、前記生成されたチャネル行列の行及び列の公約数で、前記生成されたチャネル行列の行及び列を分けて、複数の小行列に分割することを特徴とする請求項１に記載の最尤復号化方法。
【請求項５】
前記第２の手順では、前記生成されたチャネル行列の分割によって得られる複数の小行列の数が最も多くなるように、前記生成されたチャネル行列を分割することを特徴とする請求項１に記載の最尤復号化方法。
【請求項６】
複数の受信アンテナで受信した複数の受信信号から、複数の送信アンテナで送信した複数の送信信号を分離して推定する最尤復号装置であって、
前記受信信号に対応するチャネルインパルス応答から、チャネル行列を生成するチャネル行列生成部と、
前記生成されたチャネル行列を複数の小行列に分割し、前記分割された小行列に対応する前記受信信号の部分を特定し、前記分割された小行列の逆行列を用いて前記分割された小行列を変換する行列分割処理部と、
前記変換された小行列をＱＲ分解することによって三角行列を取得するＱＲ分解処理部と、
前記取得した三角行列を用いて前記各部分の受信信号を変換する信号変換部と、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ決定するＭＬＤ処理部と、を備えることを特徴とする最尤復号装置。
【請求項７】
前記最尤復号装置は、前記生成されたチャネル行列に基づいて、チャネル行列の分割数を決定し、前記決定された分割数に基づいて前記チャネル行列を分割するよう前記行列分割処理部及び前記ＱＲ分解処理部に指示する分割制御部を備えることを特徴とする請求項６に記載の最尤復号装置。
【請求項８】
前記ＭＬＤ処理部は、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ以上選択し、
前記各部分の決定された組み合わせ候補間のユークリッド距離が最も短くなるように、前記各部分の組み合わせ候補を一つ決定することを特徴とする請求項６に記載の最尤復号装置。
【請求項９】
前記行列分割処理部で用いられる逆行列は、前記分割された小行列の対角項成分に雑音成分を付加して逆行列を求めることによって生成されることを特徴とする請求項６に記載の最尤復号装置。
【請求項１０】
前記行列分割処理部は、前記生成されたチャネル行列の行及び列の公約数で、前記生成されたチャネル行列の行及び列を分けて、複数の小行列に分割することを特徴とする請求項６に記載の最尤復号装置。
【請求項１１】
前記行列分割処理部は、前記生成されたチャネル行列の分割によって得られる複数の小行列の数が最も多くなるように、前記生成されたチャネル行列を分割することを特徴とする請求項６に記載の最尤復号装置。
【請求項１２】
複数の送信アンテナによって送信された複数の送信信号を受信する複数の受信アンテナと、
前記受信アンテナで受信された複数の受信信号を増幅するＲＦ受信回路部と、
前記増幅された受信信号を復号化するベースバンド信号処理部と、を備える受信機であって、
前記ベースバンド信号処理部は、
前記振幅が増幅された受信信号から複数の送信アンテナで送信した複数の送信信号を分離して推定する最尤復号装置と、
前記分離された複数の送信信号の対数尤度を計算する計算部と、
前記計算された対数尤度に基づいて前記分離された複数の送信信号を復号化する復号化部と、を備え、
前記最尤復号装置は、
前記受信信号に対応するチャネルインパルス応答から、チャネル行列を生成するチャネル行列生成部と、
前記生成されたチャネル行列を複数の小行列に分割し、前記分割された小行列に対応する前記受信信号の部分を特定し、前記分割された小行列の逆行列を用いて前記分割された小行列を変換する行列分割処理部と、
前記変換された小行列をＱＲ分解することによって三角行列を取得するＱＲ分解処理部と、
前記取得した三角行列を用いて前記各部分の受信信号を変換する信号変換部と、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ決定するＭＬＤ処理部と、を備えることを特徴とする受信機。
【請求項１３】
前記最尤復号装置は、前記生成されたチャネル行列に基づいて、チャネル行列の分割数を決定し、前記決定された分割数に基づいて前記チャネル行列を分割するよう前記行列分割処理部及び前記ＱＲ分解処理部に指示する分割制御部を備えることを特徴とする請求項１２に記載の受信機。
【請求項１４】
前記ＭＬＤ処理部は、
前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ以上選択し、
前記各部分の決定された組み合わせ候補間のユークリッド距離が最も短くなるように、前記各部分の組み合わせ候補を一つ決定することを特徴とする請求項１２に記載の受信機。
【請求項１５】
前記行列分割処理部で用いられる逆行列は、前記分割された小行列の対角項成分に雑音成分を付加して逆行列を求めることによって生成されることを特徴とする請求項１２に記載の受信機。
【請求項１６】
前記行列分割処理部は、前記生成されたチャネル行列の行及び列の公約数で、前記生成されたチャネル行列の行及び列を分けて、複数の小行列に分割することを特徴とする請求項１２に記載の受信機。
【請求項１７】
前記行列分割処理部は、前記生成されたチャネル行列の分割によって得られる複数の小行列の数が最も多くなるように、前記生成されたチャネル行列を分割することを特徴とする請求項１２に記載の受信機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の送受信アンテナを用いて高速かつ安定な通信を行う無線通信方法及び無線通信装置に関し、特に、受信信号の最尤復号の処理に関する。
【背景技術】
【０００２】
複数の送信受信アンテナを用いて無線通信を行うＭＩＭＯ（Ｍｕｌｔｉｐｌｅ  Ｉｎｐｕｔ  Ｍｕｌｔｉｐｌｅ  Ｏｕｔｐｕｔ）方式は、高い伝送速度を達成できる。しかし、高い伝送速度を達成するために、受信信号に混在する複数の送信信号の干渉を正確に分離し、検出する処理が必要である。最尤復号（ＭＬＤ：Ｍａｘｉｍｕｍ  Ｌｉｋｅｌｉｈｏｏｄ  Ｄｅｔｅｃｔｉｏｎ）は、最も優れた特性が得られる検出方法であるが、処理が複雑になるため、膨大な演算量を必要とする。そこで、実用的な近似方法としてＱＲＭ−ＭＬＤ処理が提案されている。
【０００３】
非特許文献１には、逆行列演算（ゼロフォーシング）の処理を予め実行し、一つの送信信号の候補数を一定数（例えば、４）に限定してから全ての送信信号をＭＬＤ処理によって推定する技術が記載されている。
【０００４】
また、特許文献１には、非特許文献１に記載されている逆行列演算の処理を実行後、各送信信号のＳＮ比に応じて各送信信号の候補数を設定し、全ての送信信号をＭＬＤ処理によって推定する技術が記載されている。
【０００５】
ここで、図１０及び図１１を用いて従来のＱＲＭ−ＭＬＤ処理について説明する。
【０００６】
図１０は、従来技術のＱＲＭ−ＭＬＤ処理を実行するＱＲＭ−ＭＬＤ処理部の構成図である。
【０００７】
ＱＲＭ−ＭＬＤ処理部は、各チャネル推定部７１、チャネル行列生成部７２、ＱＲ分解処理部７３、信号変換部７４、及びＭＬＤ処理部７５を備える。
【０００８】
各チャネル推定部７１は、既知のパイロット信号を用いて、各伝搬路のチャネルインパルス応答を推定する。
【０００９】
チャネル行列生成部７２は、各チャネル推定部７１で推定されたチャネルインパルス応答を行列要素とするチャネル行列Ｈを生成する。
【００１０】
ＱＲ分解処理部７３は、チャネル行列生成部７２で生成されたチャネル行列をＱＲ分解する。例えば、送信アンテナ及び受信アンテナの数が４本ずつの場合、送信信号Ｔと受信信号Ｒとの関係はＲ＝ＨＴであり、次式で表される。
【００１１】
【数１】
0005074148.tif 000002
【００１２】
そして、チャネル行列のＱＲ分解は、Ｈ＝ＱＨ’であり、次式で表される。
【００１３】
【数２】
0005074148.tif 000003
【００１４】
ＱＲ分解は、一意の行列の変換であり、右辺の一つ目の行列Ｑは、ユニタリー行列（複素共役転置との行列積が単位行列に等しい）である。また、二つ目の行列Ｈ’は、上三角行列である。
【００１５】
次に、行列Ｑの複素共役転置行列をＱ*と表し、式（１）の両辺に左側からＱ*を乗算して左辺をＱ*Ｒ＝Ｚとおく場合、右辺はＱ*ＨＴ＝Ｑ*（ＱＨ’）Ｔ＝Ｈ’Ｔとなり、次式で表される。
【００１６】
【数３】
0005074148.tif 000004
【００１７】
信号変換部７４は、ＱＲ分解によって得られるユニタリー行列の複素共役転置行列を受信信号に乗算し、受信信号を新たな信号に変換する。例えば、式（３）で表されるように、受信信号行列Ｒに複素共役転置行列Ｑ*を乗算し、信号行列Ｚに変換する。
【００１８】
ＭＬＤ処理部７５は、ＭＬＤ処理によって送信信号を推定する。
【００１９】
次に、ＱＲ分解後のＭＬＤ処理の詳細を説明する。まず、式（３）においてｔ4に着目すると、ｚ4＝ｈ44’ｔ4が成立する。ここで、変復調にＱＰＳＫ方式を用いることを考えると、送信信号のシンボル候補は、多値数に応じて４種類となる。各シンボル候補について、ｈ44’ｔ4を計算し、ｚ4との２乗ユークリッド距離を算出する。算出されたユークリッド距離のうち、最も短いシンボルの候補が、正しい送信信号である確率が高いと推定される。次に、ｔ3に着目すると、ｚ3＝ｈ33’ｔ3＋ｈ34’ｔ4が成立するため、ｔ3及びｔ4のシンボル候補の組み合わせ（４×４＝１６種類）について、ｈ33’ｔ3＋ｈ34’ｔ4を計算し、ｚ3との２乗ユークリッド距離を算出する。ｚ3との２乗ユークリッド距離及びｚ4との２乗ユークリッド距離を合わせて、１６種類のシンボル候補のユークリッド距離を算出する。算出されたユークリッド距離のうち、最も短いシンボルの候補が正しい信号であると推定される。同様の処理をｔ1まで続けるのがＭＬＤ処理である。なお、ｔ1では、４の４乗（＝２５６）種類のシンボル候補について距離計算が必要となり、一般に、多値数Ｃのシンボルで送信アンテナ数Ｎを用いて送信する場合、ＣのＮ乗という膨大な演算量が必要となる。そこで、演算量を削減するため、Ｍアルゴリズムが適用される。
【００２０】
図１１は、従来技術のＭアルゴリズムの処理の動作図である。
【００２１】
まず、送信信号ｔ4の候補として、Ｃ1〜Ｃ4の４種類の信号レプリカを生成する。信号レプリカとは、受信装置で仮に設定される信号である。具体的には、推定されたチャネルインパルス応答に基づいて受信信号として想定された信号である。
【００２２】
次に、４種類の信号レプリカについて、送信信号ｔ3の候補を４種類ずつ信号レプリカとして生成し、１６種類の送信信号ｔ3の候補を設定する。そして、設定された各送信信号の候補と変換信号Ｚとの２乗ユークリッド距離を計算し、計算された２乗ユークリッド距離が近い順に（ｔ3、ｔ4）の組み合わせを絞り込む。例えば、図１１に示されるようにＭ＝３の場合、（ｔ3、ｔ4）の組み合わせを３種類の候補に絞り込む。
【００２３】
次に、ｔ3で絞り込まれた３種類の送信信号の候補について送信信号ｔ2の信号レプリカを生成し、２乗ユークリッド距離を計算する。
【００２４】
最後に、送信信号ｔ1についても同様に、ｔ2で絞り込まれた３種類の送信信号の候補についてＭアルゴリズムの処理を実行し、最終的に２乗ユークリッド距離が最短となる送信信号（ｔ1、ｔ2、ｔ3、ｔ4）の組み合わせを決定する。つまり、途中の処理の段階で組み合わせの候補を絞り込むことによって、最適解を見逃す可能性は生じるが、指数関数的な演算量の増大を抑えることができる。
【特許文献１】特開２００６−２１１１３１号公報
【非特許文献１】古田浩之、池田哲臣、“ミリ波モバイルカメラ用ＭＩＭＯ−ＯＦＤＭ伝送装置の試作とレイリーフェージング環境下での伝送性能の評価"、電子情報通信学会無線通信システム研究会技術報告、ＲＣＳ２００５−１４１、２００６年１月、ｐｐ．１０１−１０６
【発明の開示】
【発明が解決しようとする課題】
【００２５】
前述した従来技術では、各送信信号の複数の候補を個別に推定し、それから全ての送信信号の複数の候補を組み合わせてＭＬＤ処理を行っていたため、個別に送信信号の複数の候補を推定する場合、他の送信信号の影響が反映されず、推定解の品質が劣化する可能性が高かった。
【課題を解決するための手段】
【００２６】
本発明の代表的な一例を示せば以下の通りである。すなわち、複数の受信アンテナで受信した複数の受信信号から、複数の送信アンテナで送信した複数の送信信号を分離して推定する最尤復号化方法であって、前記受信信号に対応するチャネルインパルス応答から、チャネル行列を生成し、前記生成されたチャネル行列を複数の小行列に分割し、前記分割された小行列に対応する前記受信信号の部分を特定し、前記分割された小行列の逆行列を用いて前記分割された小行列を変換し、前記変換された小行列をＱＲ分解することによって三角行列を取得し、前記取得した三角行列を用いて前記各部分の受信信号を変換し、前記変換された受信信号に対応する前記送信信号の各部分の組み合わせ候補を一つ決定する。
【発明の効果】
【００２７】
本発明の一形態によると、各チャネル行列を小行列に分割してからＭＬＤ処理を実行することによって演算量を削減することができる。また、各小行列において複数の送信信号のＭＬＤ処理を実行して、複数の送信信号の候補を推定することによって送信信号間の影響が反映され、推定解の品質を向上することができる。
【発明を実施するための最良の形態】
【００２８】
以下、本発明の実施の形態について図面を参照しながら説明する。
【００２９】
＜実施の形態１＞
まず、本発明の第１の実施の形態のＭＩＭＯ通信システムについて説明する。
【００３０】
図１Ａは、本発明の第１の実施の形態のＭＩＭＯ通信システムの構成図である。
【００３１】
ＭＩＭＯ通信システムは、送信機１、受信機２、送信アンテナ３（３Ａ〜３Ｄ）、及び受信アンテナ４（４Ａ〜４Ｄ）を備える。以下、３Ａ〜３Ｄを集約して３として説明することもある。また、４Ａ〜４Ｄを集約して４として説明することもある
送信機１は、送信データ（ｔ1〜ｔ4）を送信信号として、異なる４本の送信アンテナ３によって、同一の周波数で同時に送信する。なお、送信信号は、各伝搬路のチャネルインパルス応答（ｈ11〜ｈ44）の影響を受ける。
【００３２】
受信機２は、チャネルインパルス応答（ｈ11〜ｈ44）の影響を受けた信号を異なる４本の受信アンテナ４によって受信する。また、受信信号には複数の送信データが混在しているため、混在している送信信号を受信信号から分離して推定する信号分離処理を実行し、推定信号（ｒ1’〜ｒ4’）、すなわち、送信信号の推定解を求める。
【００３３】
なお、図１Ａには、送信アンテナ３及び受信アンテナ４がそれぞれ４本ずつ示されているが、異なるアンテナをいくつ備えてもよい。
【００３４】
図１Ｂは、本発明の第１の実施の形態の受信機２の構成図である。
【００３５】
受信機２は、ＲＦ受信回路部５及びベースバンド信号処理部６を備える。
【００３６】
ＲＦ受信回路部５は、受信アンテナ数と同数の受信回路７（７Ａ〜７Ｄ）を備える。以下、７Ａ〜７Ｄを集約して７として説明することもある。受信回路７は、受信アンテナ４から受信した信号の振幅を増幅し、増幅した受信信号をベースバンド信号処理部に入力する。
【００３７】
ベースバンド信号処理部６は、ＱＲＭ−ＭＬＤ処理部８、対数尤度計算部９、及び復号化処理部１０を備える。
【００３８】
ＱＲＭ−ＭＬＤ処理部８は、受信信号に混在する送信信号の干渉を分離する。なお、ＱＲＭ−ＭＬＤ処理部８の詳細は、後述する図４で説明する。
【００３９】
対数尤度計算部９は、分離された信号の対数尤度を計算する。
【００４０】
復号化処理部１０は、計算された対数尤度に基づいて、復号化処理（例えば、ターボ復号化処理）を実行する。
【００４１】
本発明の第１の実施の形態のＱＲＭ−ＭＬＤ処理部８は、全体の演算量を削減するため、チャネル行列を複数の小行列に分割して変形し、変形された小行列毎にＱＲ分解する。さらに、ＱＲ分解から得られる三角行列を用いてＭＬＤ処理を実行する。そして、最後に各々のＭＬＤ処理の結果を組合せて推定解を得る。
【００４２】
以上の処理を分割ＱＲＭ−ＭＬＤ処理と呼び、信号分割数２の場合のＱＲＭ−ＭＬＤ処理部８の構成図を図２に示す。なお、信号分割数は、チャネル行列の分割に対応して送受信信号が分割される数である。
【００４３】
図２は、本発明の第１の実施の形態のＱＲＭ−ＭＬＤ処理部８の構成図である。
【００４４】
ＱＲＭ−ＭＬＤ処理部８は、各チャネル推定部２１、チャネル行列生成部２２、行列分割処理部２３、ＱＲ分解処理部２４（２４Ａ、２４Ｂ）、信号変換部２５（２５Ａ、２５Ｂ）、ＭＬＤ処理部２６（２６Ａ、２６Ｂ）及び分割制御部２８を備える。以下、２４Ａ、２４Ｂを集約して２４として説明することもある。また、２５Ａ、２５Ｂを集約して２５として説明することもある。また、２６Ａ、２６Ｂを集約して２６として説明することもある。
【００４５】
各チャネル推定部２１は、既知のパイロット信号を用いて、各伝搬路のチャネルインパルス応答を推定する。
【００４６】
チャネル行列生成部２２は、各チャネル推定部２１で推定されたチャネルインパルス応答を行列要素とするチャネル行列を生成する。
【００４７】
行列分割処理部２３は、チャネル行列生成部２２で生成されたチャネル行列を小行列に分割して変形する。
【００４８】
例えば、行列分割処理部２３は、式（１）の右辺に現れる４行４列のチャネル行列を、２行２列の四つの小行列式（４）〜（７）に分割する。
【００４９】
【数４】
0005074148.tif 000005
【００５０】
【数５】
0005074148.tif 000006
【００５１】
【数６】
0005074148.tif 000007
【００５２】
【数７】
0005074148.tif 000008
【００５３】
これらの四つの小行列を用いると、式（１）で表された４種類ずつの送受信信号は、信号分割数と同じ ２種類ずつに分割されて、式（８）及び式（９）に表される式で記述される。
【００５４】
【数８】
0005074148.tif 000009
【００５５】
【数９】
0005074148.tif 000010
【００５６】
さらに、式（８）及び式（９）を用いて送信信号（ｔ1，ｔ2）に関して変形すると、式（１０）が得られる。また、送信信号（ｔ3，ｔ4）に関して変形すると、式（１１）が得られる。
【００５７】
【数１０】
0005074148.tif 000011
【００５８】
【数１１】
0005074148.tif 000012
【００５９】
ここで、式（１０）及び式（１１）のＨ-1は、行列Ｈの逆行列を示す。
【００６０】
以上のように、送信信号（ｔ1，ｔ2，ｔ3，ｔ4）は、式（１０）と式（１１）の二つに分割される。
【００６１】
なお、チャネル行列が小行列に分割される場合には、チャネル行列の行及び列が、割り切れる数（チャネル行列の行及び列の数の公約数）でそれぞれ分けられるように分割される。そうすることによって、分割された小行列のサイズが等しくなり、並列演算処理を実施する場合に処理時間を最短にできる。また、チャネル行列の分割によって得られる小行列のサイズが最も小さく（すなわち小行列数が最も多く）なるように、チャネル行列が分割されるとよい。後述する分割されたチャネル行列に対するＱＲ分解処理では、例えば行列サイズがｎ行ｎ列においてｎの３乗程度の演算量を要する。例えばサイズ（ｎ／２）行（ｎ／２）列の４行列に分割した場合、逐次演算処理であっても、行列数は４倍となるが各行列の演算は（１／２）＾３＝１／８倍となり、全体の演算量は減少する。つまり小行列のサイズを小さくする程、演算量が削減できる。
【００６２】
ＱＲ分解処理部２４は、行列分割処理部２３で分割された各チャネル行列をＱＲ分解する。ここで、チャネル行列は、式（１０）及び式（１１）の右辺の括弧内に表される行列とする。
【００６３】
信号変換部２５は、ＱＲ分解によって得られるユニタリー行列の複素共役転置行列を受信信号に乗算し、受信信号を新たな信号に変換する。
【００６４】
ＭＬＤ処理部２６は、分割された送信信号の組合せ（ｔ1，ｔ2）及び（ｔ3，ｔ4）を、各々のＭＬＤ処理部（２６Ａ、２６Ｂ）で推定する。
【００６５】
分割制御部２８は、チャネル行列生成部２２で生成されたチャネル行列に基づいて、チャネル行列の分割を制御する。例えば、分割制御部２８は、送受信アンテナの数に依存してチャネル行列の分割数（信号分割数）を決定し、決定された分割数でチャネル行列を分割するよう行列分割処理部２３に指示する。さらに、分割制御部２８は、決定された分割数に応じて、行列分割処理部２３で分割された各小行列をＱＲ分解するよう複数（動作するＱＲ分解処理部２４と同じ数）のＱＲ分解処理部２４に指示する。
【００６６】
図３は、本発明の第１の実施の形態の分割ＱＲＭ−ＭＬＤ処理のフローチャートである。
【００６７】
分割ＱＲＭ−ＭＬＤ処理は、受信信号がＱＲＭ−ＭＬＤ処理部８に入力されたときに開始される。
【００６８】
まず、各チャネル推定部２１は、各チャネル推定部２１で送信アンテナ３と受信アンテナ４との間の伝搬路のチャネルインパルス応答を推定する（３０１）。具体的に、各チャネル推定部２１は、既知のパイロット信号を用いて伝搬路のチャネルインパルス応答を推定する。
【００６９】
次に、チャネル行列生成部２２は、ステップ３０１で推定されたチャネルインパルス応答を用いてチャネル行列を生成する（３０２）。
【００７０】
次に、行列分割処理部２３は、ステップ３０２で生成されたチャネル行列から分割チャネル行列を生成する（３０３）。分割チャネル行列を生成する具体的な処理については、図４を用いて後述する。なお、送受信アンテナの数が異なる場合には、生成されるチャネル行列が非正方型となるため、式（１０）及び式（１１）に表れる逆行列を計算するのは一般に困難となる。それを回避するための方法については後述する。
【００７１】
次に、ＱＲ分解処理部２４は、ステップ３０３で生成された各分割チャネル行列にＱＲ分解を実行する。そして、信号変換部２５は、ＱＲ分解された各分割チャネル行列に基づいて受信信号を変換する（３０４）。具体的には、ＱＲ分解処理部２４は、ステップ３０３で生成された各分割チャネル行列にＱＲ分解を実行することによって、ユニタリー行列の複素共役転置行列及び三角行列を得る。そして、信号変換部２５は、得られたユニタリー行列の複素共役転置行列を受信信号に乗算し、受信信号を変換する。
【００７２】
次に、ＭＬＤ処理部２６は、ステップ３０４で信号変換された受信信号にＭＬＤ処理を実行することによって送信信号を推定する（３０５）。
【００７３】
具体的に、ＭＬＤ処理部２６は、ステップ３０４で得られた変換信号にＭＬＤ処理を実行し、送信信号の候補を推定する。送信信号（ｔ1，ｔ2）の候補を推定する場合には、例えば、ＭＬＤ処理において設定される送信信号の信号レプリカと、ステップ３０４で変換された受信信号とのユークリッド距離が最も短い組み合わせ候補を一つ決定する。同様に、送信信号（ｔ3，ｔ4）についても組み合わせ候補を一つ決定する。そして、送信信号（ｔ1，ｔ2，ｔ3，ｔ4）が推定値として決定される。
【００７４】
次に、送受信アンテナの数が異なり、ステップ３０２で生成されるチャネル行列が非正方型となる場合の逆行列の計算方法を説明する。
【００７５】
非正方型のチャネル行列に対応するため、正方型のチャネル行列を生成してから逆行列を計算する一般化逆行列を、式（１２）に基づいて計算する。ここで、Ｈ+は行列Ｈの一般化逆行列（正方型でのＨ-1に対応する）、σnは雑音成分、Ｉは単位行列を表す。式（１２）の右辺において、σn2Ｉの加算処理を一般化逆行列演算に追加した事によって、行列の対角成分が０に近い場合に逆行列計算が困難になるという問題を回避する効果を有する。
【００７６】
【数１２】
0005074148.tif 000013
【００７７】
図４は、本発明の第１の実施の形態の分割チャネル行列作成処理のフローチャートである。
【００７８】
分割チャネル行列作成処理は、図３のステップ３０２でチャネル行列が生成された後に実行される。
【００７９】
まず、行列分割処理部２３は、信号分割数を設定する（４０１）。具体的に、行列分割処理部２３は、分割制御部２８で決定される信号分割数の情報を受信し、受信した情報を信号分割数として設定する。
【００８０】
次に、行列分割処理部２３は、チャネル行列生成部２２で生成されたチャネル行列を、ステップ４０１で設定された信号分割数に対応した小行列に分割する（４０２）。例えば、信号分割数が２の場合には、行及び列がそれぞれ二つに分割されて４個の小行列に分割される。また、信号分割数が３の場合には９個の小行列に分割される。
【００８１】
次に、行列分割処理部２３は、ステップ４０２で分割された各小行列の逆行列、又は、式（１２）に基づいた一般化逆行列を求める（４０３）。ここで、式（１２）右辺の逆行列計算は、例えば、代表的なＬＵ分解による数値計算で求めることができる。
【００８２】
次に、行列分割処理部２３は、ステップ４０２で分割された小行列とステップ４０３で求めた（一般化）逆行列とを用いて、式（１０）及び式（１１）の右辺の括弧内の計算を行い、分割チャネル行列を生成する（４０４）。そして、処理は終了する。
【００８３】
図５は、本発明の第１の実施の形態のＭＬＤ処理部２６の動作例を示す。
【００８４】
図１１に示した従来のＭアルゴリズムの動作例と比較すると、第１の実施形態のＭＬＤ処理では、送信信号（ｔ1，ｔ2）及び（ｔ3，ｔ4）が独立して決定される。独立して送信信号（ｔ1，ｔ2）及び（ｔ3，ｔ4）が決定されることによって、ツリーの深さが短縮（図５に示す例では半減）されるので、全体として演算量の削減が可能となる。
【００８５】
＜実施の形態２＞
本発明の第２の実施の形態では、第１の実施の形態において各小行列においてＭＬＤ処理を実行して複数の送信信号の推定解を残し、最後に解の品質を向上するため、残された複数の推定解を組み合せた統合ＭＬＤ処理を行う。
【００８６】
図６は、本発明の第２の実施の形態のＱＲＭ−ＭＬＤ処理部８の構成図である。
【００８７】
第１の実施の形態と異なるのは、ＱＲＭ−ＭＬＤ処理部８が、ＭＬＤ処理部２６の後に統合ＭＬＤ処理部を備える点である。
【００８８】
統合ＭＬＤ処理部２７は、各ＭＬＤ処理部２６で実行された送信信号の推定結果を組み合わせて統合ＭＬＤ処理を実行する。
【００８９】
図７は、本発明の第２の実施の形態の分割ＱＲＭ−ＭＬＤ処理のフローチャートである。
【００９０】
ステップ３０１からステップ３０４は、図３のステップ３０１からステップ３０４と対応するため説明を省略する。
【００９１】
次に、ＭＬＤ処理部２６は、送信信号（ｔ1，ｔ2）が含まれる小行列にＭＬＤ処理を実行し、送信信号（ｔ1，ｔ2）の複数の候補を選択する（７０５）。
【００９２】
次に、ＭＬＤ処理部２６は、送信信号（ｔ3，ｔ4）が含まれる別の小行列にＭＬＤ処理を実行し、同様に送信信号（ｔ3，ｔ4）の複数の候補を選択する（７０６）。
【００９３】
次に、統合ＭＬＤ処理部２７は、ステップ３０２で生成されるチャネル行列をＱＲ分解し、ＱＲ分解された結果に基づいて受信信号を信号変換する（７０７）。
【００９４】
次に、統合ＭＬＤ処理部２７は、ステップ７０７で信号変換された受信信号について、ステップ７０５及びステップ７０６で選択された送信信号（ｔ1，ｔ2）及び（ｔ3，ｔ4）の複数の組み合わせ候補を用いて統合ＭＬＤ処理を行う。そして、統合ＭＬＤ処理の結果に基づいて送信信号（ｔ1，ｔ2，ｔ3，ｔ4）を決定する。
【００９５】
具体的に、統合ＭＬＤ処理部２７は、ステップ７０５で選択された送信信号（ｔ1，ｔ2）の複数の組み合わせ候補と、ステップ７０６で選択された送信信号（ｔ3，ｔ4）の複数の組み合わせ候補とのユークリッド距離を算出し、全体のユークリッド距離が最も短くなる送信信号（ｔ1，ｔ2，ｔ3，ｔ4）の候補を推定値として決定する。
【００９６】
図８は、本発明の第２の実施の形態の統合ＭＬＤ処理部２７の動作例を示す。
【００９７】
送信信号（ｔ3，ｔ4）は、予め推定された候補を適用するのでＭＬＤ処理がスキップされる。また、送信信号（ｔ1，ｔ2）も同様に候補が推定されているので、送信信号（ｔ1，ｔ2）と（ｔ3，ｔ4）とを組み合わせるＭＬＤ処理のみが実行される。したがって、全ての送信信号について信号レプリカを作成して新たにＭＬＤ処理を実行するのに比較して、演算量が少なくて済む。また、第１の実施の形態と異なり最後に統合ＭＬＤ処理を実行することによって、全ての送信信号の影響が考慮されるため、送信信号の推定解の品質を向上することができる。
【００９８】
図９は、本発明の第１の実施の形態及び第２の実施の形態の演算量を従来技術と比較した評価図である。
【００９９】
変復調方式に１６ＱＡＭを使用し、送信信号の候補を選択する数を１６とし、送受信アンテナ数をそれぞれ４×４（送信アンテナが４本、受信アンテナが４本）及び６×６（送信アンテナが６本、受信アンテナが６本）とした場合の複素乗算数を比較する。
【０１００】
複素乗算数を比較するために考慮する処理は、ＱＲ分解、各送信信号の信号レプリカ作成及び２乗ユークリッド距離の計算である。送信アンテナ数をＮ、変復調多値数をＣ、選択候補数をＳ、信号分割数をＧとすると、複素乗算数は、ＱＲ分解ではＮの３乗程度の演算量が必要である。また、信号レプリカ作成では、Ｃ＋ＳＣ（Ｎ／Ｇ−１）（Ｎ／Ｇ＋２）／２の演算量が必要である。また、２乗ユークリッド距離の計算では、Ｃ＋ＳＣ（Ｎ／Ｇ−１）の演算量が必要である。また、分割行列作成計算では、｛（１０／Ｇ）Ｎ＾３−１３ＧＮ｝／３の演算量が必要である。
【０１０１】
従来型に着目すると、送受信アンテナ数が４×４の場合に、演算量が３２５３回であるのに対し、送受信アンテナ数が６×６の場合に、演算量が６８６２回であり、送受信アンテナ数の２乗程度で演算量が増加している。
【０１０２】
一方、本発明の第１の実施の形態で信号分割数を２として小行列数を４とした場合には、各送受信アンテナ数で演算量が１７５２回及び４２０６回であり、従来型よりも演算量が減少（３９％〜４５％減）していることが分かる。さらに、送受信アンテナ数６×６で信号分割数を３として小行列数を９とした場合には、演算量が２７０６回であり、一層の削減効果（６１％減）が得られている。
【０１０３】
また、本発明の第２の実施の形態では、統合ＭＬＤ処理が追加される分、第１の実施の形態よりも演算量が増加する。しかし、従来型と比較すると信号分割数が２の場合には、送受信アンテナ数が４×４で演算量が２０４０回（４２％減）であり、送受信アンテナ数が６×６で演算量が４５２０回（３４％減）であり、信号分割数が３の場合には、送受信アンテナ数が６×６で演算量が３２８２回（５２％減）であり、演算量が減少している事が確認できる。
【図面の簡単な説明】
【０１０４】
【図１Ａ】本発明の第１の実施の形態のＭＩＭＯ通信システムの構成図である。
【図１Ｂ】本発明の第１の実施の形態の受信機２の構成図である。
【図２】本発明の第１の実施の形態のＱＲＭ−ＭＬＤ処理部８の構成図である。
【図３】本発明の第１の実施の形態の分割ＱＲＭ−ＭＬＤ処理のフローチャートである。
【図４】本発明の第１の実施の形態の分割チャネル行列作成処理のフローチャートである。
【図５】本発明の第１の実施の形態のＭＬＤ処理部２６の動作例を示す。
【図６】本発明の第２の実施の形態のＱＲＭ−ＭＬＤ処理部８の構成図である。
【図７】本発明の第２の実施の形態の分割ＱＲＭ−ＭＬＤ処理のフローチャートである。
【図８】本発明の第２の実施の形態の統合ＭＬＤ処理部２７の動作例を示す。
【図９】本発明の第１の実施の形態及び第２の実施の形態の演算量を従来技術と比較した評価図である。
【図１０】従来技術のＱＲＭ−ＭＬＤ処理を実行するＱＲＭ−ＭＬＤ処理部の構成図である。
【図１１】従来技術のＭアルゴリズムの処理の動作図である。
【符号の説明】
【０１０５】
７１  各チャネル推定部
７２  チャネル行列生成部
７３  ＱＲ分解処理部
７４  信号変換部
７５  ＭＬＤ処理部
【図１Ａ】
0005074148.tif 000014
【図１Ｂ】
0005074148.tif 000015
【図２】
0005074148.tif 000016
【図３】
0005074148.tif 000017
【図４】
0005074148.tif 000018
【図５】
0005074148.tif 000019
【図６】
0005074148.tif 000020
【図７】
0005074148.tif 000021
【図８】
0005074148.tif 000022
【図９】
0005074148.tif 000023
【図１０】
0005074148.tif 000024
【図１１】
0005074148.tif 000025