(57)【特許請求の範囲】
【請求項１】
MIMO-OFDM通信によりデータを送受信する空間多重OFDMシステムのMIMO-OFDM通信方法において、
送信側において、送信データに時空間符号化（STBC）処理あるいは周波数空間符号化（SFBC）処理を施して複数の送信アンテナにマッピングするステップ、
前記符号化処理によりマッピングされた送信データにそれぞれOFDM変調を施して送信アンテナより送信するステップ、
を有し、受信側において、
送信装置の各送信アンテナから各受信アンテナまでのチャネル応答を要素とするチャネル応答マトリクスを推定し、該チャネル応答マトリクスに所定の変換を施すステップ、
該変換により得られたチャネルマトリクスを用いて、受信信号にZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムに従ったMIMOデコード処理を施して送信データストリームを出力するステップ、
前記チャネルマトリクスを用いて前記送信データストリームの信号対雑音比を重み係数として計算するステップ、
パスメトリックに前記重み係数を乗算することにより前記デコード処理により得られた送信データストリームに重み付け復号処理を施すステップ、
を有し、更に、受信側において、
復号処理に得られた復号データに誤りが存在するかチェックするステップ、
誤りが存在すれば送信側に再送を要求するステップ、を有し
また、送信側において、再送に際して、前記送信データストリームと該送信データストリームを送信する送信アンテナの組み合わせを変更して再送するステップ、を有し
前記送信データストリームを送信する送信アンテナの組み合わせは、各送信データストリームと該送信データストリームを送信するアンテナのすべての組み合わせの中から、スループットが最大となる組み合わせ、あるいは、最小特異値が最大となる組み合わせ、あるいはチャネルの条件数が最小となる組み合わせを選択して送信することを特徴とするMIMO-OFDM通信方法。ただし、最小特異値は、送信データストリームと該送信データストリームを送信するアンテナの組み合わせに応じたチャネルマトリックスを特異値分解して得られた対角マトリックスの対角要素の最小値、チャネルの条件数は、前記特異値分解して得られた対角マトリックスの対角要素の最大値（最大特異値）と最小特異値との比である。
【請求項２】
MIMO-OFDM通信によりデータを送受信する空間多重OFDMシステムのMIMO-OFDM通信装置において、
送信データに時空間符号化（STBC）処理あるいは周波数空間符号化（SFBC）処理を施して複数の送信アンテナにマッピングする符号化部、
前記符号化処理によりマッピングされた送信データにそれぞれOFDM変調を施して送信アンテナより送信するOFDM送信部、
送信装置の送信アンテナから受信アンテナまでのチャネル応答を要素とするチャネル応答マトリクスを推定し、該チャネル応答マトリクスに所定の変換を施すチャネル取得部、
該チャネル取得部により取得されたチャネルマトリクスを用いて、前記複数の受信アンテナで受信した受信信号にZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムに従ったMIMOデコード処理を施して送信データストリームを出力するMIMO処理部、
前記チャネルマトリクスを用いて前記送信データストリームの信号対雑音比を重み係数として計算する重み係数計算部、
パスメトリックに前記重み係数を乗算することにより前記デコード処理により得られた送信データストリームに重み付け復号処理を施す復号部、
を有し、更に、
復号処理に得られた復号データに誤りが存在するかチェックする誤り検出部、
誤りが存在すれば送信側に再送を要求する再送要求部、
再送に際して、前記送信データストリームと該送信データストリームを送信する送信アンテナの組み合わせを変更するアンテナ変更制御部、
を備え、前記アンテナ変更制御部は、
各送信データストリームと該送信データストリームを送信するアンテナのすべての組み合わせの中から、スループットが最大となる組み合わせ、あるいは、最小特異値が最大となる組み合わせ、あるいはチャネルの条件数が最小となる組み合わせを選択して送信することを特徴とするMIMO-OFDM通信装置。ただし、最小特異値は、送信データストリームと該送信デーストリームタを送信するアンテナの組み合わせに応じたチャネルマトリックスを特異値分解して得られた対角マトリックスの対角要素の最小値、チャネルの条件数は、前記特異値分解して得られた対角マトリックスの対角要素の最大値（最大特異値）と最小特異値との比である。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はMIMO−OFDM通信方法およびMIMO-OFDM通信装置に係わり、特に、MIMO−OFDM通信によりデータを送受する次世代無線ＬＡＮなどのディジタル無線通信システム及び空間多重OFDMシステムにおいてシステムの性能改善と複雑度の減少を可能にしたMIMO−OFDM通信方法およびMIMO-OFDM通信装置に関する。
【背景技術】
【０００２】
今日の無線通信システムでは、複数の送信アンテナから異なるデータストリームを並列に送信することにより、送信アンテナ数に比例して伝送容量を増大させる空間多重伝送技術が注目されている。異なる送信アンテナは、お互いに無相関になるように配置され、各アンテナから送信されるデータストリームは、それぞれ独立のフェージング伝搬路を通り、受信アンテナで受信される。
ここで、お互いに無相関になるように配置された複数の受信アンテナを利用して、多入力多出力（Multiple-Input Multiple-Output, MIMO）システムとすることにより、自由度の高いチャネル相関行列が生成でき、空間多重されたデータストリームを分離する際のSNRを向上することができる。
【０００３】
図２１はMIMOシステムの構成説明図であり、TRは送信局、RVは受信局である。送信アンテナATT1〜ATTM の数Mと同じ数のデータストリームS1〜SMが、それぞれの送信部TX1〜TXMでデータ変調・オーバーサンプリング・D/A変換・直交変調・周波数アップコンバート・帯域制限フィルタリングなどの処理を経て、各送信アンテナATT1〜ATTMから送信される。各アンテナATT1〜ATTMから送信された信号は、独立のフェージングチャネルhmn（ｍ＝１〜M，ｎ＝１〜N）を通り、空間で多重された後、N本の受信アンテナATR1〜ATＲNで受信される。各受信アンテナで受信された信号は、受信部RX1〜RXNでフィルタリング・周波数ダウンコンバート・直交検波・A/D変換処理を経て、ｙ1〜ｙNの受信データストリームが生成される。各受信データストリームは、M個の送信データストリームが多重された形になっているため、データ処理装置DPUにおいて全ての受信データストリームに対して信号処理を行うことにより、送信データストリームが分離・再生される。
【０００４】
受信信号より送信データストリームS1〜SMを分離する信号処理のアルゴリズムには、チャネル相関行列の逆行列を用いるZF(Zero-Forcing)アルゴリズムやMMSE（Minimum Mean Square Error）アルゴリズムといった線形アルゴリズムと、BLAST（Bell Laboratories Layered Space-Time）アルゴリズムに代表される非線形アルゴリズムがある。また、MLD（Maximum Likelihood Decoding）アルゴリズムなどの相関行列の逆行列演算を使用しない方法も知られている。
今、送信データストリームをM次元の複素行列Sで、受信データストリームをN次元の複素行列Yで表すと、次式の関係がある。
【数１】
0005397428.tif 000002
ここで、HはN×Mの複素チャネル行列(h11〜hNM)、VはN次元のAWGNである。
ZFアルゴリズムでは、次式により送信データストリームを推定する。
【数２】
0005397428.tif 000003
ここで、*は、行列の複素共役転置を表し、H*Hはチャネル相関行列と呼ばれる。
MMSEアルゴリズムでは、次式により送信データストリームを推定する。
【数３】
0005397428.tif 000004
ここで、ρは、受信アンテナ当りのSNRに相当する。MMSEではSNRを精度よく推定する必要が生じるが、ZFにおける雑音強調の影響を低減することができるため、一般にZFより特性が優れている。
【０００５】
MLDアルゴリズムでは、次式により送信データストリームを推定する。
【数４】
0005397428.tif 000005
ここで、Qは変調データの信号点配置の数で、QPSKでQ=4, 16QAMでQ=16, 64QAMでQ=64となる。このように、MLDでは、多値変調の演算量が膨大となり、かつ演算量は送信アンテナ数に対して指数関数的に増大してしまう。
MIMO−OFDM送受信は、図２１における送信局TRの送信部、受信局RVの受信部をOFDM（Orthogonal Frequency Division Multiplexing：直交周波数分割多重）送信部、OFDM受信部で構成したものである。
【０００６】
図２２は従来のMIMO−OFDMシステムの構成図である。送信側において、FEC エンコーダ１は畳み込み符号化により入力データストリームに誤り検出訂正用の符号化処理を施し、パンクチュア部２は符号化されたデータ列に対して符号化率に応じてパンクチュア処理を施し、空間インタリーブ部（Spatial Interleaver）３はパンクチュアしたデータビット列を複数個、図ではNs個のデータストリームに分ける。ついで、周波数インタリーブ部(Frequency Interleaver)４1〜４Nsは入力された直列データ（サブキャリア信号成分）の位置を入れ替える周波数インタリーブを行なう。コンステレーションマッパー(Constellation Mapper）５1〜５Nsは、適応制御部(図示せず)より指示されるデータ変調方式(BPSK, QPSK, 16QAM, 64QAMなど)に従ってサブキャリア数の信号成分をコンステレーションマッピングする。ついで、プリコーダ６はNs組の各並列データに所定のプリコーディングマトリックスを乗算し、Ns個のデータストリームをNt個のデータストリームに変換する。なお、プリコーディング処理は送信ダイバーシチ効果を得るため及びチャネル間干渉を軽減するために行なわれる。
IFFT(Inverse Fast Fourier Transform)部７1〜７Ntはそれぞれ直列入力するデータをサブキャリア数の並列データに変換し、しかる後、該並列データにIFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。ガードインターバル挿入部８1〜８Ntは、IFFT 部から入力するOFDM信号にガードインターバルを挿入し、送信部(TX)９1〜９Ntはガードインターバルが挿入されたOFDM信号をDA変換し、ついで、OFDM信号の周波数をベースバンドから無線帯域に変換し、高周波増幅してアンテナ１０1〜１０Ntより送信する。
【０００７】
送信アンテナ１０1〜１０Ntから送信された信号は、フェージングチャネル（伝搬路）を経て、受信装置の受信アンテナ１１1〜１１Nrにより受信され、受信回路（Rx）１２1〜１２Nrはアンテナにより受信されたＲＦ信号をベースバンド信号に変換し、該ベースバンド信号をディジタルにAD変換して出力する。シンボル切出し部１３1〜１３NrはガードインターバルGIを削除すると共にFFTタイミングでOFDMシンボルを切出してＦＦＴ部１４1〜１４Nrに入力する。FFT部１４1〜１４Nrは切り出されたOFDMシンボル毎にFFT処理を行ない、周波数領域のサブキャリアサンプルに変換する。
チャネル推定部１５は、サブキャリア毎にチャネル推定を行ない、チャネル状態情報CSI（チャネルH(k))を送信側とMIMO処理部１６に入力する。なお、送信側のプリコーダ６は通知されたチャネルに基づいてプリコーディングマトリックスを決定し、図示しない適応制御部は適応制御によりデータ変調方式を決定し、各コンステレーションマッパー５1〜５Nsに入力する。
MIMO処理部１６は推定されたチャネルを用いてサブキャリア毎に、所定の信号アルゴリズムに従って受信信号よりNs個の送信データストリームを分離出力する。
以後、コンステレーションデマッパー（Constellation Demapper）１７1〜１７Ns、周波数デインタリーブ部(Frequency Deinterleaver)１８1〜１８Ns、空間デインタリーブ部（Spatial Deinterleaver）１９、デパンクチュア部２０、復号部２１は送信側の処理と逆順にDemapping処理、周波数デインタリーブ処理、空間デインタリーブ処理、復号処理を行なって送信データを復調、復号して出力する。
【０００８】
送信側のプリコーダ６は、空間多重プリコーディング方式、固有モード転送プリコーディング方式、Limited feedback preceding方式（限定フィードバックプリコーディング方式）のいずれかに従ったプリコーディングマトリックスF(k)をNs個のデータストリームに乗算してプリコーディングを行なう。
プリコーダ６はチャネルH(k)が不明であればプリコーディングマトリクスF(k)としてidentity matrixを用い、これを送信データストリームに乗算してプリコーディングを実行する(空間多重プリコーディング方式)。また、送信側と受信側におけるチャネル状態情報(CSI)が既知であれば、プリコーダ６はチャネルマトリックスH(k)を特異値分解(SVD)し、分解されたマトリックスをプリコーディングマトリックスF(k)とし、該マトリックスF(k)を送信データストリームに乗算してプリコーディングを実行する(固有モード転送プリコーディング方式)。また、チャネルに適応したプリコーディングマトリクスF(k)が受信側からフィードバックされてくれば、プリコーダ６は該プリコーディングマトリクスF(k)を送信データストリームに乗算してプリコーディングを実行する(限定フィードバックプリコーディング方式)。
固有モード転送プリコーディング方式や限定フィードバックプリコーディング方式では、送信側がチャネル状態情報(CSI)或いは部分CSIを利用する必要がある。更に、固有モード転送プリコーディング方式の場合、送信側は適応変調制御とパワー制御を行うことが要求されている。
【０００９】
図２３は従来のMIMO−OFDM通信装置の別の構成図であり、図２１のMIMO−OFDM通信装置と異なる点は、プリコーダ６の替わりにSTBC部６′を設けている点及びSTBC部６′が送信ダイバーシチ（Diversity）を得るために、Ns個のデータストリームのそれぞれを時空間符号化（STBC：Space Time Block Coding)して、Nt本のアンテナにマッピングしている点である。すなわち、STBC部６′のSTBCエンコーダ６1′〜６Ns′は、コンステレーションマッパー５1〜５Nsからそれぞれ入力するデータにSTBCエンコード処理を施して出力する。図２４はSTBCエンコード処理の説明図であり、STBCエンコーダ６i′は、周期Tの連続する２シンボルデータ［ｘ0，ｘ1］を２系列のシンボルデータ列に変換する。第１のデータ列は［ｘ0，−ｘ1*］であり、第２のデータ列は［ｘ1，ｘ0*］である。*は複素共役を意味する。この２系列のデータがそれぞれOFDM処理されて２つの送信アンテナに入力する。図２４では１つの系列のOFDM送信構成のみを示しているが、他の系列も同様の構成を備えている。
図２２の従来のMIMO−OFDMシステム(空間多重プリコーディング方式、固有モード転送プリコーディング方式、限定フィードバックプリコーディング方式)では、高性能化するために、MLDアルゴリズムを使用する。また、図２３の従来の混成STBC空間多重OFDMシステムあるいは混成SFBC空間多重OFDMシステムにおいても高性能化するために、MLDアルゴリズムを使用する。SFBC はSpace Frequency Block Codingの略である。
しかし、MLDアルゴリズムは演算数が非常に多く、実現するために複雑度が膨大になる問題がある。そこで、複雑度を減らすため、ZFアルゴリズム或いはMMSEアルゴリズムを使用し、性能を改善するために繰り返しデコードアルゴリズム(VBLASTアルゴリズム)などを利用する。しかし、複雑度を減らすためにZFアルゴリズム或いはMMSEアルゴリズムを使用すれば誤り率性能が劣化する問題点が発生する。また、性能を改善するためにVBLASTアルゴリズムを使用すると、複雑度と処理遅延が増える問題点が発生する。また、従来のMIMO-OFDM通信では受信状態に基づいて適応変調制御やパワー制御を採用する必要があるが、これらの制御はシステムの複雑度をさらに大きくし、また、コストや消費電力などの点で好ましくない。以上から、システムの複雑度を減少し、かつ、性能を改善するMIMO-OFDM通信方法が要望されている。
従来技術としてMIMO−OFDMシステムにおいて、チャネルパラメータ推定を正確に行なうものがある（特許文献１参照）。また、別の従来技術としてMIMO-OFDMシステムにおいて、ZFとMLDを組み合わせて演算処理を簡易化する復号器がある（特許文献２参照）。さらに、別の従来技術としてMIMO-OFDMシステムにおいて、V-BLASTを用いた場合の演算処理を簡易化するものがある（特許文献３参照）。しかし、これらの従来技術はZFアルゴリズム或いはMMSEアルゴリズムによりMIMO復号の複雑度を減らし、かつ性能の向上を図るものではない。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２００２−４４０５１
【特許文献２】特表２００６−５０９３９６
【特許文献３】特開２００３−２４４１０３
【発明の概要】
【発明が解決しようとする課題】
【００１１】
以上より、本発明の目的は、受信側のMIMOデコード処理にZFアルゴリズム或いはMMSEアルゴリズムを用いても誤り率性能を改善することである。
本発明の別の目的は、簡単な信号処理でMIMO−OFDM通信(空間多重、固有モード転送、限定フィードバックプリコーディング方式)における誤り率性能を改善することである。
本発明の別の目的は、送信側で適応変調制御やパワー制御を行わなくても誤り率性能を改善することである。
本発明の別の目的は、簡単な信号処理で混成STBC空間多重OFDMシステムあるいは混成SFBC空間多重OFDMシステムにおける誤り率性能を改善することである。
本発明の別の目的は、再送性能を向上することである。
【課題を解決するための手段】
【００１２】
本願発明は、MIMO-OFDM通信によりデータを送受信する空間多重OFDMシステムのMIMO-OFDM通信方法及びMIMO-OFDM通信装置である。
・MIMO-OFDM通信方法
本発明のMIMO-OFDM通信方法は、送信側において、送信データに時空間符号化（STBC）処理あるいは周波数空間符号化（SFBC）処理を施して複数の送信アンテナにマッピングするステップ、前記符号化処理によりマッピングされた送信データにそれぞれOFDM変調を施して送信アンテナより送信するステップ、を有し、受信側において、送信装置の各送信アンテナから各受信アンテナまでのチャネル応答を要素とするチャネル応答マトリクスを推定し、該チャネル応答マトリクスに所定の変換を施すステップ、該変換により得られたチャネルマトリクスを用いて、受信信号にZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムに従ったMIMOデコード処理を施して送信データストリームを出力するステップ、前記チャネルマトリクスを用いて前記送信データストリームの信号対雑音比を重み係数として計算するステップ、パスメトリックに前記重み係数を乗算することにより前記デコード処理により得られた送信データストリームに重み付け復号処理を施すステップ、を有し、更に、受信側において、復号処理に得られた復号データに誤りが存在するかチェックするステップ、誤りが存在すれば送信側に再送を要求するステップ、を有し、また、送信側において、再送に際して、前記送信データストリームと該送信データストリームを送信する送信アンテナの組み合わせを変更して再送するステップ、を有し前記送信データストリームを送信する送信アンテナの組み合わせは、各送信データストリームと該送信データストリームを送信するアンテナのすべての組み合わせの中から、スループットが最大となる組み合わせ、あるいは、最小特異値が最大となる組み合わせ、あるいはチャネルの条件数が最小となる組み合わせを選択して送信する。ただし、最小特異値は、送信データストリームと該送信データストリームを送信するアンテナの組み合わせに応じたチャネルマトリックスを特異値分解して得られた対角マトリックスの対角要素の最小値、チャネルの条件数は、前記特異値分解して得られた対角マトリックスの対角要素の最大値（最大特異値）と最小特異値との比である。
【００１３】
・MIMO-OFDM通信装置
本発明のMIMO-OFDM通信装置は、送信データに時空間符号化（STBC）処理あるいは周波数空間符号化（SFBC）処理を施して複数の送信アンテナにマッピングする符号化部、前記符号化処理によりマッピングされた送信データにそれぞれOFDM変調を施して送信アンテナより送信するOFDM送信部、送信装置の送信アンテナから受信アンテナまでのチャネル応答を要素とするチャネル応答マトリクスを推定し、該チャネル応答マトリクスに所定の変換を施すチャネル取得部、該チャネル取得部により取得されたチャネルマトリクスを用いて、前記複数の受信アンテナで受信した受信信号にZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムに従ったMIMOデコード処理を施して送信データストリームを出力するMIMO処理部、前記チャネルマトリクスを用いて前記送信データストリームの信号対雑音比を重み係数として計算する重み係数計算部、パスメトリックに前記重み係数を乗算することにより前記デコード処理により得られた送信データストリームに重み付け復号処理を施す復号部、を有し、更に、復号処理に得られた復号データに誤りが存在するかチェックする誤り検出部、誤りが存在すれば送信側に再送を要求する再送要求部、再送に際して、前記送信データストリームと該送信データストリームを送信する送信アンテナの組み合わせを変更するアンテナ変更制御部、を備え、前記アンテナ変更制御部は、各送信データストリームと該送信データストリームを送信するアンテナのすべての組み合わせの中から、スループットが最大となる組み合わせ、あるいは、最小特異値が最大となる組み合わせ、あるいはチャネルの条件数が最小となる組み合わせを選択して送信する。ただし、最小特異値は、送信データストリームと該送信デーストリームタを送信するアンテナの組み合わせに応じたチャネルマトリックスを特異値分解して得られた対角マトリックスの対角要素の最小値、チャネルの条件数は、前記特異値分解して得られた対角マトリックスの対角要素の最大値（最大特異値）と最小特異値との比である。
【発明の効果】
【００１４】
本発明によれば、複数の受信アンテナで受信した受信信号にZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムに従ったMIMOデコード処理を施して送信データストリームを出力し、送信側のプリコーディング方式に応じた重み係数を計算し、パスメトリックに該重み係数を乗算して前記デコード処理により得られた送信データに重み付け復号処理を施すようにしたから、複雑度を抑制しつつ、誤り率性能を改善することができる。
本発明によれば、送信側のプリコーディング方式が限定フィードバックプリコーディング方式あるいは空間多重プリコーディング方式であれば、データストリームの信号対雑音比SINRを前記重み係数とするから、SINRが大きいときの復号時のパスメトリックの重み係数を大きくできるため復号における貢献度を大きくでき、結果的に誤り率が減少でき復号性能を向上できる。また、送信側と受信側のチャネルが既知であれば、前記送信側のプリコーディング方式として固有モード転送プリコーディング方式を採用でき、復号性能を向上できる。したがって、本発明によれば、MIMO−OFDM通信(空間多重、固有モード転送、限定フィードバックプリコーディング方式)における誤り率性能を改善することができる。
本発明によれば、チャネルの平均条件数を計算し、該平均条件数が閾値以下であればプリコーディング方式として、空間多重プリコーディング方式を採用し、平均条件数が閾値より大きければプリコーディング方式として、限定フィードバックプリコーディング方式あるいは固有モード転送プリコーディング方式を採用するように構成したから、簡単にプ
リコーディング方式を切り替える事ができる。
本発明によれば、信号対雑音比SINRをパスメトリックの重み係数として復号するため性能を向上でき、結果的に適応変調制御やパワー制御を行う必要がなくなる。
また本発明によれば、パンクチュア部あるいは周波数インタリーブあるいはコンステレーションマッピング部の後段に送信データストリームを複数の送信データストリームに分ける空間インタリーブ部を設けたから、送信ダイバーシチ効果を受けることができる。
本発明では、送信側において送信データをSTBC符号化し、あるいはSFBC符号化し、受信側においてチャネル応答マトリクスを変換してチャンネルマトリックスを求め、該チャネルマトリクスのSINRをViterbiのパスMetricに重み付けるようにしたから、簡単な信号処理で、誤り率性能の改善が可能となった。
また、本発明では、受信エラーが検出された場合、送信データを再送し、この再送時に送信アンテナを並べ替えて再送信を行うから、再送により正しくデータを受信できる確率を向上することができる。
更に、本発明によれば変調方式により、ZFデコード或いはMMSEデコードを選択するため、変調方式に応じた適切なデコード方式を採用できるため、誤り率性能が向上する。
また、本発明によれば、受信側の簡単な信号処理により、無線システムの複雑度、コスト、消費電力などの低減と共に性能改善が可能となった。
【図面の簡単な説明】
【００１５】
【図１】第１実施例のMIMO−OFDM送受信システムの構成図である。
【図２】s=1、Ns=3の空間インタリーブの説明図である。
【図３】受信機におけるプリコーディング方式の切替制御構成図である。
【図４】SNRに対するエラーレートを示すシミュレーション結果例である。
【図５】第２実施例のMIMO−OFDM送受信システムの構成図である。
【図６】第３実施例のMIMO−OFDM送受信システムの構成図である。
【図７】第４実施例のMIMO−OFDM送受信システムの構成図である。
【図８】第５実施例のMIMO−OFDM送受信システムの構成図である。
【図９】空間インタリーブ処理の説明図である。
【図１０】第６実施例のMIMO−OFDM送受信システムの構成図である。
【図１１】第７実施例のMIMO−OFDM送受信システムの構成図である。
【図１２】混成STBC空間多重OFDMシステムにおけるMIMO−OFDM通信装置の構成図(第８実施例)である。
【図１３】２データストリームを４アンテナにマッピングするSTBCエンコード処理説明図である。
【図１４】２データストリームを３アンテナにマッピングするSTBCエンコード処理説明図である。
【図１５】STBC−IFFT処理説明図である。
【図１６】混成SＦBC空間多重OFDMシステムにおけるMIMO−OFDM通信装置の構成図(第９実施例)である。
【図１７】２データストリームを４アンテナにマッピングするSFBCエンコード処理説明図である。
【図１８】２データストリームを３アンテナにマッピングするSFBCエンコード処理説明図である。
【図１９】SFBC−IFFT処理説明図である。
【図２０】再送機能を持たせたMIMO−OFDM通信装置の構成図(第１０実施例)である。
【図２１】MIMOシステムの構成説明図である。
【図２２】従来のMIMO−OFDMシステムの構成図である。
【図２３】従来の別のMIMO−OFDM通信装置の構成図である。
【図２４】STBCエンコード処理の説明図である。
【発明を実施するための形態】
【００１６】
（Ａ）本発明の概要
MIMO−OFDM通信システムの受信処理に適用して好適な新規な重み付きビタビ（Viterbi）デコーディング技術を提案する。MIMO−OFDM通信システムの受信装置は、MIMO復号チャネルのSINR(Signal to Interference Noise Ratio)を計算し、計算されたSINRを信号の信頼度とみなし、該SINRをViterbiデコーディングのパスメトリック（path metric）の重み係数として用いる。この方法によれば、適応変調制御やパワー制御を行わなくても、簡単な信号処理で、すなわち、低複雑度のZFデコードアルゴリズムあるいはMMSEデコードアルゴリズムによる簡単な信号処理で受信データの誤り率を減少できる。
【００１７】
（Ｂ）第１実施例
（B-1）全体の動作
図１は第１実施例のMIMO−OFDM通信システムの構成図である。なお、送信側はNt 個の送信アンテナを有し、受信側はNr 個の受信アンテナを有するものとし、また、送信データストリーム数をNs であるとし、Ns≦min[Nr, Nt]である。
送信側において、FEC エンコーダ５１は畳み込み符号化により入力データストリームに誤り検出訂正用の符号化処理を施し、パンクチュア部５２は符号化されたデータ列に対して符号化率に応じてパンクチュア処理を施し、空間インタリーブ部（Spatial Interleaver）５３はパンクチュアしたデータビット列を複数個、図ではNs個のデータストリームに分ける。具体的に、空間インタリーブ部５３はsビットの連続ブロックを周期的に各データストリームに送る。ここでsは以下の三つの中で一つ選択すればよい。
0005397428.tif 000006
ただし、NBPSCはサブキャリア当たりのCodedビットの数である。例えば、３)の場合、データ変調がBPSK変調であればs=1、QPSK変調であればs=2、１６QAM変調であればs=4である。図２はs=1、Ns=3の空間インタリーブの説明図である。
周波数インタリーブ部(Frequency Interleaver)５４1〜５４Nsは入力された直列のデータを並列データに変換してサブキャリア数の信号成分を作成し、各サブキャリア信号成分の位置を入れ替える周波数インタリーブを行なう。コンステレーションマッパー(Constellation Mapper）５５1〜５５Nsは所定のデータ変調方式、たとえば、QPSK変調に従ってサブキャリア数の各信号成分をコンステレーションマッピングする。
【００１８】
ついで、プリコーダ５６はNs個の各データストリームに所定のプリコーディングマトリックスを乗算し、Ns個のデータストリームをNt個のデータストリームに変換する。すなわち、プリコーダ５６はプリコーディングマトリクスF(k)を使って、以下の式で、Ns 個送信データストリームをNt 個送信アンテナにマッピングする。
【数５】
0005397428.tif 000007
ただし、kはサブキャリア番号、s(k)はNs×１の送信データベクトルで、x(k)はNt×１の送信アンテナデータベクトルである。
IFFT部５７1〜５７Ntはそれぞれ並列入力するサブキャリアサンプルにIFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。ガードインターバル挿入部５８1〜５８Ntは、IFFT 部から入力するOFDM信号にガードインターバルを挿入し、送信部(TX)５９1〜５９Ntはガードインターバルが挿入されたOFDM信号をDA変換し、ついで、OFDM信号の周波数をベースバンドから無線帯域に変換し、高周波増幅してアンテナ６０1〜６０Ntより送信する。
送信アンテナ６０1〜６０Ntから送信された信号は、フェージングチャネル（伝搬路）を経て、受信装置の受信アンテナ６１1〜６１Nrにより受信され、受信回路（Rx）６２1〜６２Nrはアンテナにより受信されたＲＦ信号をベースバンド信号に変換し、該ベースバンド信号をディジタルにAD変換して出力する。シンボル切出し部６３1〜６３NrはガードインターバルGIを削除すると共にFFTタイミングでOFDMシンボルを切出してFFT部６４1〜６４Nr に入力する。FFT部６４1〜６４Nr は切り出されたOFDMシンボル毎にFFT処理を行ない、周波数領域のサブキャリアサンプルに変換する。
チャネル推定回路６５は周知の方法によりサブキャリア毎にチャネル推定を行ない、コードワード決定部８１はコードブックより最適のプリコーディングマトリクスF(k)に対応するコードワードを求め、そのインデックス番号を送信側のプリコーダ５６に通知する。これにより、プリコーダ５６は受信側と同一のコードブックを用いてプリコーディングマトリクスF(k)を求め、（5)式の演算を行なう。また、チャネル推定部６５は推定した各サブキャリアのチャネルマトリックスH(k)をMIMO処理部６６及びSINR計算部７２に入力する。
【００１９】
FFT部６４1〜６４Nr から出力する周波数領域の信号y(k)は以下の式で表現できる。
【数６】
0005397428.tif 000008
ただし、H(k)はk番目のサブキャリアのチャネルマトリックス、n(k)はNr×1の雑音ベクトルを表し、y(k)はNr×1の受信信号ベクトルを表わし、各受信アンテナの信号は送信信号の重なりとなっている。y(k)とＨ(k)とＦ(k)は受信機において測定により既知である。送信信号
0005397428.tif 000009
を得るために、MIMO処理部６６はMLDアルゴリズムを使えるが、複雑度が膨大となる。そこで、複雑度の小さい次式で示される線形デコーダを採用する。
【数７】
0005397428.tif 000010
ただし、Ｇ(k)はNs×Nrのマトリックスである。ZF(Zero Forcing)線形デコーダを使う場合、Ｇ(k)は次式
【数８】
0005397428.tif 000011
で表現される。(.)+はマトリックスpseudo-inverse を表す。
一方、MMSE (Minimum Mean Square Error) 線形デコーダを使う場合、Ｇ(k)は次式
【数９】
0005397428.tif 000012
(.)-1はマトリックス反転、( )′は複素共役転置を表す。また、N0はアンテナ受信雑音電力、εsはk番目のサブキャリアのトータルの送信電力、INsはNs×Nsの単位行列である。
【００２０】
SINR計算部７２は後述するようにチャネルマトリックスH(k)、プリコーディングマトリクスF(k)を用いてSINRi(k)を計算して復号部７１に入力する。
コンステレーションデマッパー（Constellation Demapper）６７1〜６７Ns、周波数デインタリーブ部(Frequency Deinterleaver)６８1〜６８Ns、空間デインタリーブ部（Spatial Deinterleaver）６９、デパンクチュア部７０、復号部７１は、(7)式により分離した
0005397428.tif 000013
の各データストリームに対して、送信側の処理と逆順にDemapping処理、周波数デインタリーブ処理、空間デインタリーブ処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。なお、復号部７１は入力されたSINRを信号の信頼度とみなし、該SINRをViterbiデコーディングのpath metricに重み付け乗算し、重み付の誤り訂正デコーダ処理を実行する。
【００２１】
（B-2）プリコーディング方式
プリコーディングマトリクスF(k)を決定する方式には、(1)限定フィードバックプリコーディング方式と、(2)固有モード転送プリコーディング方式と、(3)空間多重プリコーディング方式とがある。チャネルH(k)が未知であれば、プリコーダ５６は空間多重プリコーディング方式を採用し、プリコーディングマトリクスF(k)として単位行列(identity matrix)を用いて(5)式の演算を行い、空間多重転送を行なう。また、送信側と受信側のチャネルH(k)が既知であれば、プリコーダ５６は固有モード転送プリコーディング方式を採用し、チャネルマトリックスH(k)を特異値分解(SVD)し、分解された右特異マトリックスをプリコーディングマトリックスF(k)として(5)式の演算を行なう。また、チャネルH(k)に基づいてプリコーディングマトリクスF(k)を特定するインデックスが決定されて受信側からフィードバックされてくれば、限定フィードバックプリコーディング方式を採用し、プリコーダ５６は該インデックスに応じたプリコーディングマトリクスF(k)を求めて(5)式の演算を行なう。以下、各システムについて詳述する。
【００２２】
（ａ） 限定フィードバックプリコーディング方式
限定フィードバックプリコーディング方式の原理を説明する。事前にデータストリーム数Ns、送信アンテナ数Ntとフィードバックビット数(L)とに基づいてコードブックを作成しておき、送受信機が該コードブックを有しているものとする。
受信機のコードワード決定部８１は、チャネルH(k)を使って、コードブックの中から最適なコードワードを選択し、Limited feedback channel(限定フィードバックチャネル)を利用して、コードワードのインデックスを送信機にフィードバックする。送信機はフィードバックされたインデックスを使って、コードブックから、プリコーディングマトリックスF(k)を選択し、該マトリックスF(k)を送信データストリームに乗算してプリコーディング処理を行う。コードブックは予め決めたユニタリーマトリックス（unitary matrix）の集合であり、コードワードはコードブックの中の１つのユニタリーマトリックスを定義する。ユニタリーマトリックスはUであれば、次式
U′U=I
を満足する。なお、限定フィードバックプリコーディング方式は、次の論文に記述されている。
IEEE TRANSACTION ON INFORMATION THEPRY, VOL.51, NO.8, AUGUST ”David J. Love and Robert W. Heath. Jr, Limited Feedback Unitary Precoding for Spatial Multiplexing Systems" 2005Page 2967-2976
【００２３】
以下は、限定フィードバックプリコーディング方式の詳細動作手順である。
(1)まず、受信機のコードワード決定部８１はコードワードを選択する。このコードワードの選択は以下のように行なう。すなわち、チャネル推定部６５において、チャネル推定を行い、チャネルH(k)を求める。ついで、コードワード決定部８１は求められたチャネルH(k)を使って、下記1)〜3)の基準を利用して、最適なコードワードを選択する。
（1-1) MMSE(Minimum mean square error )選択基準
受信機のMSE (Mean square error )は、Ｈ(k)とＦ(k)を既知として、以下の式で表現できる。
【数１０】
0005397428.tif 000014
コードワードはMSEのTrace(tr)或いはDeterminant(det)を最小化するように、以下の式に基づいて
【数１１】
0005397428.tif 000015
コードブックの中から選択される。なお、Trace(tr)はマトリックスの対角要素の加算である。
上式に替えて、コードワードは簡略化した以下の式を使って求めることができる。
【数１２】
0005397428.tif 000016
【００２４】
（1-2） 容量(Capacity )選択基準
コードワード決定部８１は、容量(Capacity )、換言すれば所定時間あたりの通信情報量を最大化するように、次式
【数１３】
0005397428.tif 000017
によりコードブックの中からコードワードを選択することができる。または、簡略化した以下の式を使うこともできる。
【数１４】
0005397428.tif 000018
注目すべき点は、原理上、(12b)、(14)式から最大容量選択基準と最小Det(MSE)基準とが同等であることである。
【００２５】
（1-3） 特異値(Singular value)選択基準
コードワード決定部８１は、最小特異値(Minimum singular value)を最大化するように、次式
【数１５】
0005397428.tif 000019
によりコードブックの中からコードワードを選択する。(15)式は、コードブックの中のコードワードの各候補Fj(k)についてマトリックスH(k)Fj(k)を求め、該マトリックス[H(k)Fj(k)]を特異値分解して最小の特異値λminを求め、全最小特異値のうち最大の最小特異値を与える候補Fj(k)を意味している。
以上では、特異値分解によりコードワードを決定する場合について説明したが、固有値分解（Eigen Value Decomposition）によっても同様にコードワードを決定することができる。
(2)コードワードが求まれば、コードワード決定部８１は該コードワードのインデックスを送信機のプリコーダ５６にフィードバックする。すなわち、受信機は選択された最適なコードワードのインデックス（Ｌビット）をlimited feedbackチャネルで送信機に転送する。
【００２６】
(3)送信機のプリコーダ５６はインデックスを受信すれば、内蔵のコードブックを参照して受信したインデックスに応じたマトリックスを選択する。そして、プリコーダ５６はこの選択したマトリックスをMIMO プリコーディングマトリックスとして使用する。
(4)受信機のMIMO処理部６６は以上のコードワード選択処理と並行して、送信信号
0005397428.tif 000020
をZF(Zero Forcing)線形デコーダあるいはMMSE (Minimum Mean Square Error) 線形デコーダを用いて分離する。
【００２７】
(5)また、SINR計算部７２は、ビタビ（Viterbi）復号処理におけるパスメトリック（Path metric）の重みとして用いるSINRを計算する。
(5-1) ZFデコーダを使用する場合のSINRの計算
SINR計算部７２は、送信信号分離にZFデコーダ(ZFデコードアルゴリズム)を使用する場合には、各サブキャリアのSINRを以下の式で計算して復号部７１に入力する。
【数１６】
0005397428.tif 000021
ただし、iはデータストリームのインデックスである。また[・]i,iはマトリックスのi行i列の要素である。
復号部７１は計算されたSINRをビタビ（Viterbi）のパスメトリック（Path metric）に次式で示すように重み付けする。
【数１７】
0005397428.tif 000022
ただし、iはデータストリームのインデックス、nはデマップしたビットのインデックスである。例えば、16QAM変調の場合、b1〜b4の4ビットにデマップされるから、n=1であればb1を意味し、n=2であればb2を意味する。また、
0005397428.tif 000023
はソフトDecisionの値であり、
0005397428.tif 000024
は参考値である。
以上のようにビタビのパスメトリックをSINRで重み付けすることにより、ZFデコードアルゴリズムによる雑音増強を低減して、性能を改善することができる。
【００２８】
(5-2) MMSEデコーダを使用する場合のSINRの計算
以上はZFデコードアルゴリズムを用いた場合であるが、SINR計算部７２は、送信信号分離にMMSEデコーダ(MMSEデコードアルゴリズム)を使用する場合には、各サブキャリアのSINRを以下の式で計算して復号部７１に入力する。
【数１８】
0005397428.tif 000025
復号部７１は計算されたSINRをビタビ（Viterbi）のパスメトリック（Path metric）に次式で示すように重み付けする。
【数１９】
0005397428.tif 000026
以上のようにビタビのパスメトリックをSINRで重み付けすることにより、MMSEデコードアルゴリズムによる雑音増強を低減して、性能を改善することができる。
なお、MMSEデコーダされた信号に(16)式で計算した
0005397428.tif 000027
を重みとして使用することできる。逆に、ZFデコーダされた信号に(18)式で計算した
0005397428.tif 000028
を重みとして使用することできる。
また、複雑度を減らすために、SINRがある閾値以下の場合、SINR=0にし、SINRがある閾値以上の場合、SINRを1としてもよい。
【００２９】
(5-3)ビタビパスメトリックに重み付けした効果
(7)式に(6)式を代入すると以下の式になる。
【数２０】
0005397428.tif 000029
線形デコーダ ZFとMMSEでは、上式の第２項G(k)n(k)が雑音増強となる。G(k)n(k)が大きい場合は、デコードした信号は雑音が大きく増強され、線形デコーダの性能が劣化する。特に、周波数選択性MIMOチャネルの場合は、この問題が大きくなる。本発明ではビタビパスメトリックにSINRを重み付けすることにより、これらの雑音増強の低減が可能になる。すなわち、本発明では各データストリーム、各サブキャリアのSINRを計算し、このSINRを復調した信号の信頼度として採用し、式(17),(19)のように、重み付きビタビ復号を行う。これにより、SINRが小さい場合(雑音増強が大きい)、ビタビのパスメトリックの貢献度を小さくし、逆に、SINRが大きい場合(雑音増強が小さい)、ビタビのパスメトリックの貢献度を大きくし、結果的に、性能改善を図ることができる。
【００３０】
(ｂ)固有モード転送プリコーディング方式
固有モード転送プリコーディング方式は、送信機と受信機のチャネル状態情報CSI(チャネルH（k）)が既知である場合に採用できる。固有モード転送プリコーディング方式においては、チャネルマトリックスH（k）を特異値分解（Singular Value Decomposition: SVD）し、この分解により得られるマトリックスをプリコーディングマトリックスとして用いて、高性能のシステムを構築する。
すなわち、サブキャリア kに対し、チャネルマトリックスＨ(k)は、Nr×Ntの行列である。Ｈ(k)の特異値分解は次式で表現できる。
【数２１】
0005397428.tif 000030
(21a)式におけるＵ(k)は左特異マトリックスで、Ｖ(k)は右特異マトリックスと呼ばれる。Ｄ(k)は対角マトリックスであり、(21b)式に示す対角要素を有している。λ1(k),….,λm(k)は固有モードの利得に比例するもので、大きい順に順序付けられており、第１固有モードλ1(k)が最大、第ｍ固有モードλm(k)が最小の利得を有する。すなわち、λ1(k)> λ2(k) >…. >λm(k)である。ｍはチャネルＨ(k)のrankである。最適なプリコーディングマトリックスは右特異マトリックスＶ(k)の左側のNs列である。
【００３１】
また、上記の分析から、理想的な場合、すなわち、送信機と受信機のチャネル状態情報CSI(チャネルH（k）)が既知で、かつ完全に一致する場合、固有値λi(k) を次式で示すようにビタビ復号におけるパスメトリック（Path metric）の重みとする。すなわち、復号部７１は上記固有値λi(k)をパスメトリックの重みとして、該パスメトリックに乗算して復号することにより、ZFデコードアルゴリズムやMMSEデコードアルゴリズムによる雑音増強を低減して、性能改善ができる。
【数２２】
0005397428.tif 000031
なお、固有モード転送プリコーディング方式を実現するにはチャネル量子化誤差や時変のチャネルなどの問題が存在し、実現困難な場合が多い。かかる場合には、限定フィードバックプリコーディング方式によりSINRを求めてパスメトリックの重みとするのが現実的である。
【００３２】
(ｃ) 空間多重プリコーディング方式
空間多重プリコーディング方式は、送信側でH(k)が未知の場合にプリコーディングマトリクスF(k)として単位行列（identity matrix）を用いる方式である。すなわち、空間多重プリコーディング方式の場合、プリコーディングマトリックスＦ(k)はidentity matrixとなる。このため、SINR計算部７２は次式によりSNIRを簡単に計算することができる。
【数２３】
0005397428.tif 000032
ただし、(23a)式はZFデコードアルゴリズムを採用する場合、(23b)式はMMSEデコードアルゴリズムを用いる場合である。
復号部７１は上記計算されたSINRをパスメトリックの重みとして、該パスメトリックに乗算して復号することにより、すなわち、(17)式あるいは(19)式に従って復号することによりZFデコードアルゴリズムやMMSEデコードアルゴリズムによる雑音増強を低減して、性能を改善することができる。
【００３３】
（ｅ）システムの切替制御
・第1の切替制御
チャネルの条件と通信装置の能力によりプリコーディング方式を、限定フィードバックプリコーディング方式、固有モード転送プリコーディング方式、空間多重プリコーディング方式に切り換えてOFDM−MIMO通信をすることが可能である。
具体的に、 (24)式を用いてサブキャリア毎にチャネルの条件数を計算し、その平均値をチャネル平均条件数として求め、該チャネル平均条件数が良い場合、すなわち、チャネル平均条件数がある閾値以下の場合に空間多重プリコーディング方式を使用し、チャネル平均条件数が悪い場合、すなわち、チャネルの平均条件数が閾値以上の場合には、通信装置の能力により、固有モード転送プリコーディング方式か限定フィードバックプリコーディング方式を使用する。
【００３４】
チャネルの条件数は次式
【数２４】
0005397428.tif 000033
で計算できる。上式において、分子はチャネルマトリックスH(k)を特異値分解して得られる対角マトリックスの対角要素の最大値(最大特異値)であり、分母は対角要素の最小値(最小特異値)である。チャネル条件数が小さいほど各伝搬路のチャネル(伝搬特性)は均一であり、受信状態が良好であるといえる。なお、信頼性の通信を確保するため、場合により、データストリーム数を減らすことも可能である。
【００３５】
以上では、特異値分解によりチャネル条件数を算出する場合について説明したが、固有値分解（Eigen Value Decomposition）によっても同様にチャネル条件数を算出することができる。
また、以上では、限定フィードバックプリコーディング方式、固有モード転送プリコーディング方式、空間多重プリコーディング方式を切り換える場合であるが、チャネル平均条件数に基づいて空間多重プリコーディング方式と限定フィードバックプリコーディング方式の切替を行ない、固有モード転送プリコーディング方式を使用しないようにもできる。
また、プリコーダは、チャネルH(k)が未知の初期時に空間多重プリコーディング方式を使用し、チャネルH(k)に基づいてプリコーディングマトリックスが求まるようになったら限定フィードバックプリコーディング方式を使用するように切替制御することもできる。
図３は受信機におけるプリコーディング方式の切替制御構成図であり、２００aは受信部、２００bは送信部である。受信部２００aには、チャネル推定部６５、コードワード決定部８１、コードブック８２に加えて、チャネル条件数演算部８３、プリコーディング方式決定部８４が設けられている。チャネル条件数演算部８３は(24)式に基づいてチャネル平均条件数を演算し、プリコーディング方式決定部はチャネル平均条件数に基づいて空間多重プリコーディング方式を使用するか、限定フィードバックプリコーディング方式を使用するか決定して、プリコーディング指示データPCを送信部２００bを介して送信機のプリコーダ５６に送信する。プリコーダ５６は指示されたプリコーディング方式に基づいてプリコーディングマトリックスを決定して(5)式の演算を行なう。
【００３６】
・第2の切替制御
MIMO処理部６６はデコードアルゴリズムとしてZFデコードアルゴリズム或いはMMSEデコードアルゴリズムを使用できるが、変調方式(Modulation)により、ZFデコードアルゴリズムとMMSEデコードアルゴリズムを選択的に使用することが可能である。例えば、BPSK，QPSKの場合にはMMSEデコードアルゴリズムを使い、16QAM、64QAMの場合にはZFデコードアルゴリズムを使う。これは以下の理由による。
すなわち、MMSEデコードアルゴリズムは雑音に強く、ZFデコードアルゴリズムは雑音に対してMMSEデコードアルゴリズムほど強くない。BPSK，QPSKは雑音の強いときに採用する変調方式であり、従ってBPSK，QPSK変調方式に際してはMMSEデコードアルゴリズムを使用する。また、16QAM、64QAMは受信品質が良好なときに採用する変調方式であり、従って16QAM、64QAM 変調方式に際してはZFデコードアルゴリズムを使用する。
【００３７】
(Ｂ-３) シミュレーション結果
図４は信号対雑音比SNRに対するエラーレートを示すシミュレーション結果例であり、縦軸はパケットエラーレートPER、横軸は信号対雑音比（SNR）を示している。（Ａ）はZFデコーダのシミュレーション結果、（Ｂ）はMMSEデコーダのシミュレーション結果で、変調方式／符号化率が16QAM-1/2の場合である。また、空間多重プリコーディング方式の場合、送信アンテナ数は２、受信アンテナ数は２である。固有モード転送プリコーディング方式と限定フィードバックプリコーディング方式（LFP）の場合、送信データストリーム数、送信アンテナ数、受信アンテナ数はそれぞれ２、４、２である。
SP１，SP2は空間多重プリコーディング方式の特性であり、SP1は重み無し、SP2は重み付けした場合である。LFP1，LPF2は限定フィードバックプリコーディング方式の特性であり、LFP1は重み無し、LFP2は重み付けした場合である。MP1,MP2は固有モード転送プリコーディング方式の特性であり、MP1は重み無し、MP2は重み付けした場合である。
図８（Ａ），（ｂ）から明らかなように、ZFデコードアルゴリズム,MMSEデコードアルゴリズムに基づいてMIMO処理する場合、何れの場合も、重み係数付けを行なうとすべてのプリコーディング方式において性能改善が確認できる。
また、３つのプリコーディング方式のうち、固有モード転送プリコーディング方式が最も性能が良く、次に限定フィードバックプリコーディング方式の性能が良く、空間多重プリコーディング方式の性能は３番目である。
【００３８】
（B-４）MISO,SIMO,SISOシステムへの適用
以上では、MIMO（Multiple Input Multiple Output）の場合について説明したが、本発明はMISO (Multiple Input Single Output) およびSIMO（Single Input Multiple Output）にも適用できる。
MISO或いはSIMOの場合、SINRは次式
【数２５】
0005397428.tif 000034
により計算できる。したがって、復号部７１はViterbiデコーディングにおける Path metricを次式
【数２６】
0005397428.tif 000035
で表現して復号処理を行なう。なお、上式において、iは送信アンテナ(MISOの場合)或いは受信アンテナ(SIMOの場合)のインデックスである。
また、本発明はSISO (Single Input Single Output)にも適用できる。SISOの場合、SINRは次式
【数２７】
0005397428.tif 000036
により計算できる。したがって、復号部７１はViterbiデコーディングにおける Path metricを次式
【数２８】
0005397428.tif 000037
で表現して復号処理を行なう。
【００３９】
（Ｃ）第２実施例
図５は第２実施例のMIMO−OFDM送受信システムの構成図であり、図１の第１実施例と同一部分には同一符号を付している。第１実施例と異なる点は、空間インタリーブする前に周波数インタリーブする点である。
送信側は、FEC エンコーダ５１、パンクチュア部５２、周波数インタリーブ部５４をそれぞれ1つ備えている。空間インタリーブ部５３は空間ダイバーシチを利用するために、周波数インタリーブされたデータビットをNs個のデータストリームに分けてコンステレーションマッパー５５1〜５５Nsに入力する。コンステレーションマッパー５５1〜５５Nsはコンステレーションマッピングし、プリコーダ５６は所定の方式に従ってプリコーディングする。以降の手順は第1実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、送信側の処理と逆順にコンステレーションデマッピング処理、空間デインタリーブ処理、周波数デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
第２実施例では、周波数インタリーブ部５４の回路規模を小さくでき、しかも、空間ダイバーシチを利用できる。
【００４０】
（Ｄ）第３実施例
図６は第３実施例のMIMO−OFDM送受信システムの構成図であり、図１の第１実施例と同一部分には同一符号を付している。第１実施例と異なる点は、空間インタリーブする前に周波数インタリーブ、コンステレーションマッピングする点である。
送信側は、FEC エンコーダ５１、パンクチュア部５２、周波数インタリーブ部５４、コンステレーションマッパー５５をそれぞれ1つ備えている。空間インタリーブ部５３は空間ダイバーシチを利用するために、コンステレーションマッピングされたデータ列をNs個のデータストリームに分けてプリコーダ５６に入力し、プリコーダ５６は所定の方式に従ってプリコーディングする。以降の手順は第1実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、送信側の処理と逆順に空間デインタリーブ処理、コンステレーションデマッピング処理、周波数デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
第３実施例では、周波数インタリーブ部５４、コンステレーションマッパー５５の回路規模を小さくでき、しかも、空間ダイバーシチを利用できる。
【００４１】
（Ｅ）第４実施例
図７は第４実施例のMIMO−OFDM送受信システムの構成図であり、図１の第１実施例と同一部分には同一符号を付している。第１実施例と異なる点は、Ns個の送信データストリームが入力し、各送信データストリームに対応して、FEC エンコーダ５１1〜５１Ns、パンクチュア部５２1〜５２Ns、周波数インタリーブ部５４1〜５４Ns、コンステレーションマッパー５５1〜５５Nsが設けられ、空間インタリーブ部が存在しない点である。
プリコーダ５６はコンステレーションマッピングされた信号に対して所定の方式に従ってプリコーディングする。以降の手順は第１実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、分離されたデータストリート毎に送信側の処理と逆順にコンステレーションデマッピング処理、周波数デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
第４実施例によれば、各データストリームを並列処理できるため、回路の動作周波数を低減することが可能である。ただし、空間ダイバーシチは利用できないため、特に固有モード転送プリコーディング方式の場合に性能劣化が若干ある。
【００４２】
（Ｆ）第５実施例
図８は第５実施例のMIMO−OFDM送受信システムの構成図であり、図７の第４実施例と同一部分には同一符号を付している。第４実施例と異なる点は、パンクチュア部５２1〜５２Nsと周波数インタリーブ部５４1〜５４Nsの間に空間インタリーブ部５３を設けた点である。
空間インタリーブ部５３は空間ダイバーシチを利用するために、各データストリームをパンクツアリングした後に各データストリームを並べ替える。そして、並べ替えた各データストリーム毎に周波数インタリーブ、コンステレーションマッピングを施してプリコーダ５６に入力する。プリコーダ５６はコンステレーションマッピングされた信号に対して所定の方式に従ってプリコーディング処理を施す。以降の手順は第１実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、分離されたデータストリート毎に送信側の処理と逆順にコンステレーションデマッピング処理、周波数デインタリーブ処理、空間デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
図９は空間インタリーブ処理の説明図であり、入力データストリーム、出力データストリームがそれぞれ３の例である。空間インタリーブ部５３は、まず、第１入力データストリームの最初の３つのデータを３つの出力データストリームに１個づつ分配し、次に、第２データストリームの最初の３つのデータを３つの出力データストリームに分配し、次に
、第３データストリームの最初の３つのデータを３つの出力データストリームに分配し、以後、同様の並び替えを行なって空間インタリーブ処理を行なう。
第５実施例によれば、回路規模は第４実施例より少し増えるが、各データストリームが並列処理できるため回路動作周波数を低減することが可能であり、また、空間ダイバーシチを利用することができる。
【００４３】
（Ｇ）第６実施例
図１０は第６実施例のMIMO−OFDM送受信システムの構成図であり、図７の第４実施例と同一部分には同一符号を付している。第４実施例と異なる点は、周波数インタリーブ部５４1〜５４Nsとコンステレーションマッパー５５1〜５５Nsとの間に空間インタリーブ部５３を設けた点である。
空間インタリーブ部５３は空間ダイバーシチを利用するために、各データストリームに周波数インタリーブ処理を施した後に各データストリームを並べ替える。そして、並べ替えた各データストリーム毎にコンステレーションマッピングを施してプリコーダ５６に入力する。プリコーダ５６はコンステレーションマッピングされた信号に対して所定の方式に従ってプリコーディング処理を施す。以降の手順は第１実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、分離されたデータストリート毎に送信側の処理と逆順にコンステレーションデマッピング処理、空間デインタリーブ処理、周波数デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
第６実施例によれば、回路規模は第４実施例より少し増えるが、各データストリームが並列処理できるため回路動作周波数を低減することが可能であり、また、空間ダイバーシチを利用することができる。
【００４４】
（Ｈ）第７実施例
図１１は第７実施例のMIMO−OFDM送受信システムの構成図であり、図７の第４実施例と同一部分には同一符号を付している。第４実施例と異なる点は、コンステレーションマッパー５５1〜５５Nsとプリコーダ５６の間に空間インタリーブ部５３を設けた点である。
空間インタリーブ部５３は空間ダイバーシチを利用するために、各データストリームにコンステレーションマッピング処理を施した後に各データストリームを並べ替えてプリコーダ５６に入力する。プリコーダ５６はコンステレーションマッピングされた信号に対して所定の方式に従ってプリコーディング処理を施す。以降の手順は第１実施例と同じである。
受信側において、MIMO処理部６６はMMSEデコードアルゴリズム或いはZFデコードアルゴリズムに従ってデータストリームを分離し、以後、分離されたデータストリート毎に送信側の処理と逆順に空間デインタリーブ処理、コンステレーションデマッピング処理、周波数デインタリーブ処理、デパンクチュア処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
第７実施例によれば、回路規模は第４実施例より少し増えるが、各データストリームが並列処理できるため回路動作周波数を低減することが可能であり、また、空間ダイバーシチを利用することができる。
【００４５】
（Ｉ）第８実施例
図１２は混成STBC空間多重OFDMシステムにおけるMIMO−OFDM通信装置の構成図で、送信側はNt 個の送信アンテナを有し、受信側はNr 個の受信アンテナを有している。なお、送信データストリーム数をNs であるとすれば、Ns≦min[Nr, Nt]である。
送信側において、FEC エンコーダ１５１は畳み込み符号化により入力データストリームに誤り検出訂正用の符号化処理を施し、パンクチュア部１５２は符号化されたデータ列に対して符号化率に応じてパンクチュア処理を施し、空間インタリーブ部（Spatial Interl
eaver）１５３はパンクチュアしたデータビット列を複数個、図１２ではNs個のデータストリームに分ける。具体的に、空間インタリーブ部１５３はsビットの連続ブロックを周期的に各データストリームに送る（図２参照)。
周波数インタリーブ部(Frequency Interleaver)１５４1〜１５４Nsは入力された直列データ（サブキャリア信号成分）の位置を入れ替える周波数インタリーブを行なう。コンステレーションマッパー(Constellation Mapper）１５５1〜１５５Nsは所定のデータ変調方式、たとえば、QPSK変調に従ってサブキャリア数の各信号成分をコンステレーションマッピングする。
【００４６】
ついで、STBC部１５６のSTBCエンコーダ１５６1〜１５６Nsは、入力する直列データにSTBCエンコード処理を施し、Ns個のデータ系列をNt個のデータ系列に変換してNt本の送信アンテナにマッピングする。図１３は２データストリームを４アンテナにマッピングする例、同様に図１４は２データストリームを３アンテナにマッピングする例である。図１３の例では、コンステレーションマッパー１５５1から入力された周期Tの連続する２シンボルデータ［s11，s12］は、STBCエンコーダ１５６1において２系列のシンボルデータ列［s11，−s12 *］、［s12，s11*］に変換されて送信アンテナ１６０1、１６０2にマッピングされる。また、コンステレーションマッパー１５５2から入力された周期Tの連続する２シンボルデータ［s21，s22］は、STBCエンコーダ１５６2において２系列のシンボルデータ列［s21，−s22 *］、［s22，s21*］に変換されて送信アンテナ１６０3、１６０4にマッピングされる。図１４の例では、コンステレーションマッパー１５５1から入力された周期Tの連続する２シンボルデータ［s11，s12］は、STBCエンコーダ１５６1において２系列のシンボルデータ列［s11，−s12 *］、［s12，s11*］に変換されて送信アンテナ１６０1、１６０2にマッピングされるが、コンステレーションマッパー１５５2から入力された周期Tの連続する２シンボルデータ［s21，s22］はSTBCエンコード処理されずにそのまま送信アンテナ１６０3にマッピングされる。
【００４７】
IFFT(Inverse Fast Fourier Transform)部１５７1〜１５７NtはSTBC部１５６から直列入力するデータをサブキャリア数の並列データに変換し、しかる後、該並列データにIFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。図１５はIFFT処理説明図である。STBCエンコーダ１５６1は、図１５(A)に示すように、コンステレーションマッパー１５５1から入力する連続する２OFDMシンボルに対して上記STBCエンコード処理を行い、それぞれ２系列発生する。IFFT部１５７1のシリパラ変換部１５７SPは図１５(B)に示すように、STBCエンコーダ１５６1から入力する直列のSTBCエンコードデータを時刻t0におけるｋ個の並列データ、時刻t1におけるｋの並列データに変換し、IFFT処理部１５７IFFTは各並列データに順次IFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。IFFT部１５７2もIFFT部１５７1と同様の処理を行う。また、STBCエンコーダ１５６2、IFFT部１５７3〜１５７4はSTBCエンコーダ１５６1、IFFT部１５７1〜１５７2と同様に動作する。
図１２に戻って、ガードインターバル挿入部１５８1〜１５８Ntは、IFFT 部から入力するOFDM信号にガードインターバルを挿入し、送信部(TX)１５９1〜１５９Ntはガードインターバルが挿入されたOFDM信号をDA変換し、ついで、OFDM信号の周波数をベースバンドから無線帯域に変換し、高周波増幅してアンテナ１６０1〜１６０Ntより送信する。
【００４８】
送信アンテナ１６０1〜１６０Ntから送信された信号は、フェージングチャネル（伝搬路）を経て、受信装置の受信アンテナ１６１1〜１６１Nrにより受信され、受信回路（Rx）１６２1〜１６２Nrはアンテナにより受信されたRF信号をベースバンド信号に変換し、該ベースバンド信号をディジタルにAD変換して出力する。シンボル切出し部１６３1〜１６３NrはガードインターバルGIを削除すると共にFFTタイミングでOFDMシンボルを切出してFFT部１６４1〜１６４Nrに入力する。FFT部１６４1〜１６４Nrは切り出されたOFDMシンボル毎にFFT処理を行ない、周波数領域のサブキャリアサンプルに変換する。
図１３(A)のSTBCエンコード処理する場合、サブキャリアf1で送信される信号s11，s12，s21，s22の受信信号y11，y12，y21，y22は以下の式で表現できる。
【数２９】
0005397428.tif 000038
(29)式において、hi,jは図１３(Ｂ)に示すようにｊ番目の送信アンテナからi番目の受信アンテナの周波数f1におけるチャンネル応答である。si,jはi番目のデータストリームにおけるj番目のシンボルを表す。yi,jはi番目の受信アンテナにおけるj番目の受信シンボルを表す。wi,jはi番目の受信アンテナにおけるj番目の受信シンボルAWGN雑音を表す。
同様にサブキャリアf2で送信される信号s13，s14，s23，s24の受信信号y13，y14，y23，y24は以下の式で表現できる。なお、hi,j′はｊ番目の送信アンテナからi番目の受信アンテナの周波数f2におけるチャンネル応答である。
【数３０】
0005397428.tif 000039
以上より、各サブキャリア成分の受信信号は(29)式と同等の式で表現できる。
(29)式から簡単な変換で次式が成立する。この変換は、(29)式の左辺における２行２列のマトリクスを４行１列のマトリクスにするための変換である。なお、以下の(31)〜(40)式ではサブキャリアのインデックスｋを付加しない一般表現になっているが、サブキャリア毎に成立する式である。
【数３１】
0005397428.tif 000040
ただし、
【数３２】
0005397428.tif 000041
【数３３】
0005397428.tif 000042
【数３４】
0005397428.tif 000043
【数３５】
0005397428.tif 000044
であり、(.)Tはマトリックス転置を表す。
また、Ｈｅのチャンネル相関プロファイルは以下の特性を持つ。
【数３３ａ】
0005397428.tif 000045
その中に、
【数３３ｂ】
0005397428.tif 000046
【数３３ｃ】
0005397428.tif 000047
同様に、図１４から簡単な変換で、以下の式
【数３６】
0005397428.tif 000048
【数３７】
0005397428.tif 000049
が得られる。yeとweはそれぞれ(34)式、(35)式と同じである。以上では、Ns=2，Nt=4およびNs=2，Nt=３の場合であるが、２以上の送信データストリーム、３以上の送信アンテナにも適用することができることは明らかである。
【００４９】
受信機においてyeとHeが既知であるから、MIMO処理部１６６は送信信号を求めるために、Maximum likelihood (ML)を使用することができる。しかし、複雑度が膨大となるため複雑度の小さい線形デコーダを採用する。線形デコーダにおいては次式
【数３８】
0005397428.tif 000050
が成立する。ここで、Geは2Ns×2Nrのマトリックスである。ZF(Zero Forcing)線形デコーダを使う場合、Geは次式
【数３９】
0005397428.tif 000051
により与えられる。ただし、(.)+はマトリックスpseudo-inverse を表す。(.)'は、複素共役転置を表す。
式３９の中の逆マトリックスは以下の式で表現できる。
【数３９ａ】
0005397428.tif 000052
式３９aに示したように、実際に逆マトリックスの計算は必要がなくなり、複雑度も減らすことができる。
一方、MMSE (Minimum Mean Square Error) 線形デコーダを使う場合、Geは次式
【数４０】
0005397428.tif 000053
により与えられる。ただし、(.)-1はマトリックス反転、( )′は複素共役転置を表す。また、Noはアンテナ受信雑音電力、εsはk番目のサブキャリアのトータルの送信電力、I2Nsは2Ns×2Nsの単位行列である。
式39aと同様に式４０の中の逆マトリックスは以下の式で表現できる。
【数４０ａ】
0005397428.tif 000054
式40aに示したように、実際に逆マトリックスの計算は必要がなくなり、複雑度も減らすことができる。
【００５０】
SINR計算部１７２は後述するようにチャネルマトリックスHeを用いてSINRj(k)を計算して復号部１７１に入力する。コンステレーションデマッパー（Constellation Demapper）１６７1〜１６７Ns、周波数デインタリーブ部(Frequency Deinterleaver)１６８1〜１６８Ns、空間デインタリーブ部（Spatial Deinterleaver）１６９、デパンクチュア部１７０、復号部７１は、(38),(39)式により、あるいは(38)、(40)式によりMIMO処理部１６６が分離した
0005397428.tif 000055
の各データストリームに対して、送信側の処理と逆順にDemapping処理、周波数デインタリーブ処理、空間デインタリーブ処理、復号処理(重み付の誤り訂正デコーダ処理)を行う。
復号部１７１はSINR計算部１７２から入力されたSINRを信号の信頼度とみなし、該SINRをViterbiデコーディングのpath metricに重み付け乗算し、重み付の誤り訂正デコーダ処理を実行する。
(1) ZFデコーダを使用する場合のSINRの計算及びビタビ復号
SINR計算部１７２は、送信信号分離にZFデコーダ(ZFデコードアルゴリズム)を使用する場合、各サブキャリアｋのSINRを以下の式で計算して復号部１７１に入力する。
【数４１】
0005397428.tif 000056
【数４２】
0005397428.tif 000057
なお、jは変換されたチャンネルの各レイヤのインデックスである。ここで、レイヤのインデックスとは、(33)式あるいは(37)式の正方形マトリクスにおける列のインデックスを意味する。すなわち、上記正方形マトリクスの列はレイヤと定義する。また、[・]j,jはマトリックスのｊ行ｊ列の要素である。復号部１７１は計算されたSINRを次式で示すようにビタビ（Viterbi）のパスメトリック（Path metric）に重み付け(乗算)する。
【数４３】
0005397428.tif 000058
ただし、nはデマップしたビットのインデックスである。例えば、16QAM変調の場合、b1〜b4の4ビットにデマップされるから、n=1であればb1を意味し、n=2であればb2を意味する。iはデータストリームのインデックスである。また、
0005397428.tif 000059
はソフトDecisionの値であり、
0005397428.tif 000060
は参考値である。以上のようにビタビのパスメトリックをSINRで重み付けすることにより、ZFデコードアルゴリズムによる雑音増強を低減して、性能を改善することができる。
【００５１】
(2)MMSEデコーダを使用する場合のSINRの計算及びビタビ復号
SINR計算部１７２は、送信信号分離にMMSEデコーダ(MMSEデコードアルゴリズム)を使用する場合には、各サブキャリアｋのSINRを以下の式で計算して復号部１７１に入力する。
【数４４】
0005397428.tif 000061
【数４５】
0005397428.tif 000062
復号部１７１は計算されたSINRを次式で示すようにビタビ（Viterbi）のパスメトリック（Path metric）に重み付け(乗算)する。
【数４６】
0005397428.tif 000063
以上のようにビタビのパスメトリックをSINRで重み付けすることにより、MMSEデコードアルゴリズムによる雑音増強を低減して、性能を改善することができる。
なお、MMSEデコーダされた信号に(41)式で計算した
0005397428.tif 000064
を重みとして使用することできる。逆に、ZFデコーダされた信号に(44)式で計算した
0005397428.tif 000065
を重みとして使用することができる。
また、複雑度を減らすために、SINRがある閾値以下の場合、SINR=0にし、SINRがある閾値以上の場合、SINRを1としてもよい。
線形デコーダ ZFと線形デコーダ MMSEを使用すると雑音増強となるが、本発明ではビタビパスメトリックにSINRを重み付けすることにより、これらの雑音増強の低減が可能になる。すなわち、本発明では変換されたチャネルの各レイヤ、各サブキャリアのSINRを計算し、このSINRを復調した信号の信頼度として採用し、式(43),(46)のように、重み付きビタビ復号を行う。これにより、SINRが小さい場合(雑音増強が大きい)、ビタビのパスメトリックの貢献度を小さくし、逆に、SINRが大きい場合(雑音増強が小さい)、ビタビのパスメトリックの貢献度を大きくし、結果的に、性能改善を図ることができる。
【００５２】
（Ｊ）第９実施例
図１６は混成SＦBC空間多重OFDMシステムにおけるMIMO−OFDM通信装置の構成図であり、図１２の第８実施例と同一部分には同一符号を付している。異なる点は、STBC部１５６に変えてSFBC(Space Frequency Block Coding)部１８１を設けた点である。
SFBC部１８１のSFBCエンコーダ１８１1〜１８１Nsは、入力する直列データにSFBCエンコード処理を施し、Nsデータ系列をNtデータ系列に変換してNt本の送信アンテナにマッピングする。図１７は２データストリームを４アンテナにマッピングする例、同様に図１８は２データストリームを３アンテナにマッピングする例である。図１７の例では、コンステレーションマッパー１５５1から入力された連続する２サブキャリアデータ［s11，s12］は、SFBCエンコーダ１８１1においてサブキャリア[f1，f2]のシンボルデータ列［s11，s12 ］、［−s12*，s11*］に変換されて送信アンテナ１６０1、１６０2にマッピングされる。また、コンステレーションマッパー１５５2から入力された連続する２サブキャリアデータ［s21，s22］は、SFBCエンコーダ１８１2においてサブキャリア[f1，f2]のシンボルデータ列［s21，s22 ］、［−s22*，s21*］に変換されて送信アンテナ１６０3、１６０4にマッピングされる。図１８の例では、コンステレーションマッパー１５５1から入力された連続する２サブキャリアデータ［s11，s12］は、SFBCエンコーダ１８１1においてサブキャリア[f1，f2]のシンボルデータ列［s11，s12 ］、［−s12*，s11*］に変換されて送信アンテナ１６０1、１６０2にマッピングされるが、コンステレーションマッパー１５５2から入力された連続する２サブキャリアデータ［s21，s22］はSFBCエンコード処理されずにそのまま送信アンテナ１６０3にマッピングされる。
IFFT部１５７1〜１５７NtはSFBC部１８１から直列入力するデータをサブキャリア数の並列データに変換し、しかる後、該並列データにIFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。図１９はIFFT処理説明図である。SFBCエンコーダ１８１1は、図１９(A)に示すように、コンステレーションマッパー１５５1から入力する連続する２OFDMシンボルに対してそれぞれ上記SFBCエンコード処理を行い、時刻t0におけるサブキャリアfi， fi+1成分をそれぞれ２系列発生する。IFFT部１５７1のシリパラ変換部１５７SPは図１９(B)に示すように、SFBCエンコーダ１８１1から入力する直列のSFBCエンコードデータを時刻t0におけるｋ個の並列データに変換し、IFFT処理部１５７IFFTは該並列データにIFFT（逆フーリエ変換）処理を施して合成し、離散時間信号(OFDM信号)にして出力する。IFFT部１５７2もIFFT部１５７1と同様の処理を行う。また、SFBCエンコーダ１８１2、IFFT部１５７3〜１５７4はSFBCエンコーダ１８１1、IFFT部１５７1〜１５７2と同様に動作する。他の部分の動作は第１実施例と同じである。
【００５３】
（Ｋ）第１０実施例
図２０は再送機能を持たせたMIMO−OFDM通信装置の構成図であり、図１２の第８実施例と同一部分には同一符号を付している。異なる点は、誤り検出部１９１、送信アンテナ決定部１９２、ACK/NACK/送信アンテナ情報を送信する送信部１９３、送信アンテナ切替制御部１９４、送信アンテナ切替部１９５を設けた点である。
受信データにエラーが含まれるとき、受信装置は送信装置にデータの再送を要求し、送信装置は再送要求に対してデータの再送を行う。この場合、再送を前回の送信時と同一のアンテナを用いて行なうと再びエラーが発生する可能性が高い。すなわち、MIMOシステムにおいて再送制御を行う場合、再送データを同一のアンテナから送信すると、伝送路状態の悪いアンテナを使い続けることになるため、再送による誤り率の改善が小さくなり、再送制御の利得が得られにくくなる。特に、フェージングやマルチパス環境の変化が、再送間隔（Round Trip Time）に比べて遅い場合に、この問題は深刻である。そこで、第１０実施例において送信装置は、再送による誤り率が減少するように各データを再送するアンテナの組み合わせを変更してデータの再送を行う。なお、受信装置は再受信したデータに誤りがある場合、前に受信したデータと合成して復号を行う。
図２０において、誤り検出部１９１は復号結果に誤りが存在するかチェックし、誤りが存在する場合には、送信アンテナ決定部１９２は後述するアンテナ並べ替え選択基準に従って各再送データと送信アンテナの最適な組み合わせ決定し、送信部１９３は誤りの有無（ACK/NACK）と誤り有りの場合には送信アンテナの組み合わせ情報を送信装置に送信する。送信アンテナ制御部１９４は受信装置から送られてくる送信アンテナの組み合わせ情報に基づいて送信アンテナ切替部１９５を制御してSTBC部１５６から出力するSTBCデータ系列を所定のアンテナに振り分ける。
図１３の４送信アンテナの場合、アンテナ並べパターンの全ての組み合わせは次の６個がある。
｛(1,2)(3,4)]、 [(1,2)(4,3)]、[(1,3)(2,4)]、[(1,3)(4,2)]、[(1,4)(2,3)]、[(1,4)(3,2)]
このため、組み合わせ情報として、各サブキャリアのフィードバックビット数は３ビットが必要である。図１４の３送信アンテナの場合は、アンテナ並べパターンの全ての組み合わせには次の３個がある。
[1,2,3]、[1, 3,2]、[3,1,2]。
したがって、組み合わせ情報として、各サブキャリアのフィードバックビット数は２ビットが必要である
【００５４】
以下に、送信アンテナの最適な組み合わせを決定するためのアンテナ並べ替え選択基準について説明する。
(１) 容量(Capacity )選択基準
受信装置は、容量(Capacity )が最大となるように、各送信データと該送信データを送信するアンテナの全ての組み合わせの中から最良の組み合わせを選択する。 すなわち、次式を満足するアンテナ組み合わせを選択する。
【数４７】
0005397428.tif 000066
Hs は、送信データと該送信データを送信するアンテナの組み合わせに応じたチャンネルマトリックス、H は全組み合わせに応じた変換されたチャネルマトリックスの集合、detは行列のdeterminantである。 容量とは、例えば、スループット(単位時間当たりのデータ転送量)である。
【００５５】
(２)特異値(Singular value)選択基準
受信装置は、最小特異値(Minimum singular value)を最大化するように（最小SINRの最大化と同等）、各送信データと該送信データを送信するアンテナの全ての組み合わせの中から最良の組み合わせを選択する。すなわち、次式を満足するアンテナ並べパターンを選択する。
【数４８】
0005397428.tif 000067
λmin｛Hs｝は、チャネルマトリックスHｓを特異値分解して得られる対角マトリックスの対角要素の最小値(最小特異値)である。
【００５６】
(３)条件数(Condition Number)選択基準：
受信装置は、チャンネルの条件数を最小化するように、各送信データと該送信データを送信するアンテナの全ての組み合わせの中から最良の組み合わせを選択する。すなわち、次式を満足するアンテナ並べパターンを選択する。
【数４９】
0005397428.tif 000068
上式において、分子はチャネルマトリックスHsを特異値分解して得られる対角マトリックスの対角要素の最大値(最大特異値)であり、分母は対角要素の最小値(最小特異値)である。チャネル条件数が小さいほど各伝搬路のチャネル(伝搬特性)は均一であり、受信状態が良好であるといえる。なお、信頼性の通信を確保するため、場合により、データストリーム数を減らすことも可能である。
【００５７】
(４) MMSE(Minimum mean square error )選択基準
受信装置のMSE (Mean square error )は、以下の式で表現できる。
【数５０】
0005397428.tif 000069
MSEのTrace(tr)或いはDeterminant(det)を最小化するように、以下の式により、各送信データと該送信データを送信するアンテナの全ての組み合わせの中から最良の組み合わせを選択する。
【数５１】
0005397428.tif 000070
【数５２】
0005397428.tif 000071
注目すべきは、原理上、最大容量選択基準が最小Det(MSE)基準と同等であることである。
【符号の説明】
【００５８】
１５１  FEC エンコーダ
１５２  パンクチュア部
１５３  空間インタリーブ部
１５４1〜１５４Ns  周波数インタリーブ部
１５５1〜１５５Ns  コンステレーションマッパー
１５６  STBC部
１５６1〜１５６Ns  STBCエンコーダ
１５７1〜１５７Nt  IFFT部
１５８1〜１５８Nt  ガードインターバル挿入部
１５９1〜１５９Nt  送信部
１６０1〜１６０Nt  送信アンテナ
１６１1〜１６１Nr  受信アンテナ
１６２1〜１６２Nr  受信回路
１６３1〜１６３Nr  シンボル切出し部
１６４1〜１６４Nr  FFT部
１６５  チャネル推定部
１６６  MIMO処理部
１６７1〜１６７Ns  コンステレーションデマッパー
１６８1〜１６８Ns  周波数デインタリーブ部
１６９  空間デインタリーブ部
１７０  デパンクチュア部
１7１  復号部
１７２  SINR計算部
【図１】
0005397428.tif 000072
【図２】
0005397428.tif 000073
【図３】
0005397428.tif 000074
【図４】
0005397428.tif 000075
【図５】
0005397428.tif 000076
【図６】
0005397428.tif 000077
【図７】
0005397428.tif 000078
【図８】
0005397428.tif 000079
【図９】
0005397428.tif 000080
【図１０】
0005397428.tif 000081
【図１１】
0005397428.tif 000082
【図１２】
0005397428.tif 000083
【図１３】
0005397428.tif 000084
【図１４】
0005397428.tif 000085
【図１５】
0005397428.tif 000086
【図１６】
0005397428.tif 000087
【図１７】
0005397428.tif 000088
【図１８】
0005397428.tif 000089
【図１９】
0005397428.tif 000090
【図２０】
0005397428.tif 000091
【図２１】
0005397428.tif 000092
【図２２】
0005397428.tif 000093
【図２３】
0005397428.tif 000094
【図２４】
0005397428.tif 000095