JP2006067070A - Ｍｉｍｏシステム受信方法及びその装置 - Google Patents

Ｍｉｍｏシステム受信方法及びその装置 Download PDF

Info

Publication number
JP2006067070A
JP2006067070A JP2004245098A JP2004245098A JP2006067070A JP 2006067070 A JP2006067070 A JP 2006067070A JP 2004245098 A JP2004245098 A JP 2004245098A JP 2004245098 A JP2004245098 A JP 2004245098A JP 2006067070 A JP2006067070 A JP 2006067070A
Authority
JP
Japan
Prior art keywords
transmission
signal
transmission symbol
likelihood
estimating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2004245098A
Other languages
English (en)
Inventor
Bokee Vladimir
ボケー ウラジミール
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2004245098A priority Critical patent/JP2006067070A/ja
Priority to GB0506429A priority patent/GB2417651B/en
Priority to US11/105,435 priority patent/US7492830B2/en
Publication of JP2006067070A publication Critical patent/JP2006067070A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/02—Arrangements for detecting or preventing errors in the information received by diversity reception
    • H04L1/06—Arrangements for detecting or preventing errors in the information received by diversity reception using space diversity
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
    • H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
    • H04L25/03057—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a recursive structure
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03184—Details concerning the metric
    • H04L25/03197—Details concerning the metric methods of calculation involving metrics
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03248—Arrangements for operating in conjunction with other apparatus
    • H04L25/03292—Arrangements for operating in conjunction with other apparatus with channel estimation circuitry
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03312—Arrangements specific to the provision of output signals
    • H04L25/03324—Provision of tentative decisions

Abstract

【課題】 複数のアンテナを利用して信号を受信する受信装置において、復号処理の複雑さを抑えつつ良好な復号性能を得る。
【解決手段】 受信信号から送信シンボルを推定する（ステップＳ４）。推定した送信シンボルの確からしさを計算する（ステップＳ５、Ｓ６）。確からしさを表す情報から送信シンボルの平均値および分散を計算し、さらにそれらに基づいてチャネルを補償するためのＭＭＳＥ係数を計算する（ステップＳ１〜Ｓ３）。ＭＭＳＥ計算を利用して再び送信シンボルを推定する（ステップＳ４）。推定処理は、２回以上繰り返し実行される。
【選択図】図３

Description

本発明は、複数のアンテナを利用して無線通信を行うシステムに係わり、特に、ＭＩＭＯシステムの受信装置においてチャネルを補償する技術に係わる。
高速無線通信を実現するための技術としてＭＩＭＯ（Multi Input Multi Output）が知られている。ＭＩＭＯは、空間分割多重通信方式の１つであり、送信装置および受信装置にそれぞれ複数のアンテナを設け、複数の電波を空間で隔てながら互いに並列に伝送する技術である。ここで、複数のアンテナから送信される複数の電波は、周波数利用効率を上げるため一般には同じ周波数が用いられる。
図１９は、ＭＩＭＯシステムの基本構成図である。図１９において、送信装置は、Ｎt 本の送信アンテナを利用して信号を送信し、受信装置は、Ｎr 本の受信アンテナを利用して送信装置から出力された信号を受信する。この場合、送信装置と受信装置との間のチャネルの特性は、Ｈ0,0〜ＨNr-1,Nt-1により表される。なお、「Ｈ0,0」は、送信アンテナ０から受信アンテナ０へ信号を伝送するためのチャネルの特性（例えば、channel frequency response）を表示する。また、「ＨNr-1,Nt-1」は、送信アンテナNt-1から受信アンテナNt-1へ信号を伝送するためのチャネルの特性を表示する。
そうすると、受信装置において各受信アンテナを介して得られる信号ｙは、下記（１）式により得られる。なお、送信シンボルｘは、送信アンテナＮ0〜ＮNt-1を利用して送信されるシンボルｘ0〜ｘNt-1から構成される。同様に、受信信号ｙは、受信アンテナＮ0〜ＮNr-1を利用して受信される信号ｙ0〜ｙNr-1から構成される。そして、送信シンボルｘおよび受信信号ｙは、それぞれ複素ベクトルで表される。また、チャネル状態情報Ｈは、Ｎt ×Ｎr の行列で表される。さらに、「ｎ」は、各受信アンテナにおけるノイズである。
Figure 2006067070
 なお、数式（１）および後続のいくつかの数式で登場する「ｒ」は、受信信号を表しており、「ｙ」と表記されるものと同じものを表している。また、数式（１）および後続のいくつかの数式でピリオド印を使用するが、これは、乗算を表している。
上記（１）式で記述されるＭＩＭＯ処理は、シンボル毎に実行される。なお、ＭＩＭＯシステムは、例えば、特許文献１に記載されている。
＜非符号化ＭＩＭＯ＞
空間時間ブロック符号（Space Time Block Coding）または空間周波数ブロック符号（Space Frequency Block Coding）などの符号化ＭＩＭＯとは異なり、非符号化ＭＩＭＯは、階層化処理を利用することでデータレートを高めることができる。なお、これらの技術については、下記の非特許文献１に記載されている。
[1] A. Van Zelst, “Space division multiplexing algorisms” Proc. 10th Med. Electrotechnical Conference 2000, pp. 1218-1221
５つの一般的なＭＩＭＯ信号処理技術について記載する前に、まず、ＭＩＭＯ送信におけるチャネル容量について簡単に説明する。なお、ここでは、「Ｈ」は、複数のチャネルの状態を表現する行列である。
受信装置においてチャネル行列Ｈが既知であるものとすると、空間分割多重システムのシャノン限界は、下記（２）式で表される。ここで、「ＩNr」は、Ｎr 次元の単位行列である。また、「ρ」は、各送信アンテナから送信された信号のＳＮＲ（Signal Noise Ratio）である。さらに、「ＨH」は、Ｈのエルミート行列（Hermitian Matrix）である。
Figure 2006067070
 ＭＩＭＯ処理を実行するために、下記の２つの仮定を導入する。
仮定１：システムの設計について
数学的に解を得るためには、ＭＩＭＯシステムは、「Ｎr >=Ｎt 」という条件で設計される必要がある。即ち、受信アンテナの数は、送信アンテナの数以上であるものとする。
仮定２：チャネルの状態について
各送信アンテナまたは受信アンテナは、全部または一部が非相関でなければならない。すなわち、チャネル行列Ｈの全要素が互いに一致していることは許されない。なお、この条件は、フェージングが存在する環境、或いは散乱が大きな環境においては満たされることになる。
次に、５つの一般的なＭＩＭＯ信号処理技術について記載する。
＜Zero Forcingアルゴリズム＞
Zero Forcingアルゴリズムでは、受信装置において、送信信号の推定値ｘestは、下記（３）式で得られる。ここで、「Ｈ+」は、Ｈの擬似逆行列（pseudo inverse）である。なお、Ｈが正方行列であれば、「Ｈ+＝Ｈ-1」である。
Figure 2006067070
 ＜最小平均二乗誤差（ＭＭＳＥ：Minimum Mean Square Error）アルゴリズム＞
ＭＭＥＳアルゴリズムでは、受信装置において、送信信号の推定値ｘestは、下記の（４）式で得られる。（４）式において、「α」は「１／ＳＮＲ」を表し、「Ｉ」は単位行列を表す。また、「Ｈ*」は、Ｈの共役転置行列である。
Figure 2006067070
 なお、ＭＭＳＥアルゴリズムについては、下記の非特許文献２に記載されている。
[2] J. Proakis, Digital Communications, 3rd ed. McGraw-Hill 1995
＜ＺＦ−ＢＬＡＳＴ（Zero Forcing Bell Laboratories Layered Space Time）アルゴリズム＞
ＢＬＡＳＴは、複数の送信アンテナを利用して高速データ送信を実現することを目的としている。このアルゴリズムによれば、Ｎt 本の送信アンテナを利用して同期送信されたシンボルが検出される。干渉を抑制／キャンセルすることで、各送信アンテナから送信された信号が選択されて推定される。このキャンセルステップは、Ｎt −１回繰り返し実行される。以下に、復号手順について記載する。
１回目の復号ステップ（ｕ＝１）においては、下記の条件が設定される。なお、「〜」印は、推定値を意味する。また、Ｈの上に付されている「＝」印は、Ｈが行列であることを表す。
Figure 2006067070
 各ステップｕ毎に、ヌリング行列（nulling matrix）Ｇが計算される。ここで、ｕ回目の復号ステップのヌリング行列Ｇu は、下記（５）式で表されるものとする。すなわち、ｕ回目の復号ステップのヌリング行列Ｇu は、ｕ回目の復号ステップにおけるＨ行列の推定値の擬似逆行列である。
Figure 2006067070
 ＢＬＡＳＴアルゴリズムでは、各ステップｕ毎に、ヌリング行列Ｇから任意のレイヤを選択することができる。ただし、上記非特許文献２によれば、ＳＮＲに対して最大のポスト検出信号を示すレイヤを選択した場合に最高のパフォーマンスが得られる。この条件は、ヌリング行列Ｇからその対角上において最小ノルムを持った行を選択することと等価である。
ｕ回目の復号ステップにおいて、ヌリング行列Ｇu から最小ノルムを持った列を選択する処理は、下記（６）式または（７）式で表される。なお、「ｋu 」は、ヌリング行列Ｇu から選択すべき列を識別すると共に、対応する送信アンテナに相当する。また、「ｗku」は、ヌリング行列Ｇu のｋu 番目の列の各成分から構成されるベクトルを意味する。
Figure 2006067070
 （７）式で得られる「ｗku」に受信信号ベクトルｙuを掛けると、送信アンテナｋu に対応するレイヤを除くすべてのレイヤが抑圧される。このとき、下記（８）式が得られる。ここで、「ｘku」は送信信号を表し、「ｅｓｔ」は推定値であることを示す。
Figure 2006067070
 １つのレイヤが検出されると、以降のレイヤを検出する処理を進めることができる。受信信号ベクトルから検出された信号の一部を差し引くことにより、次のステップで対象となるレイヤ数が１つ小さくできる。即ち、受信ベクトルの変化は、下記（９）式で表される。なお、「Ｑ」は受信シンボルのスライス値を表す。「ｙu」および「ｙu+1」はそれぞれｕ回目のステップにおける受信信号、ｕ＋１回目のステップにおける受信信号を表している。「（Ｈu）ku」は、ｕ回目のステップにおけるＨ行列のｋu 番目の列（即ち、送信アンテナｋu に係わるチャネルの状態を表示する係数から構成される列）を表す。
Figure 2006067070
 この後、チャネル「（Ｈu）ku」のｋu 番目の列はもはや不要になるので、この列は取り除かれる。この処理は、下記（１０）式で表すことができる。
Figure 2006067070
 ＢＬＡＳＴアルゴリズムでは、「ｕ」を１ずつインクリメントしながら、上述の（６）式〜（１０）式の処理が「ｕ＝Ｎt 」を満たすまで繰り返し実行される。すなわち、上述の（６）式〜（１０）式の処理は、送信アンテナの数だけ繰り返し実行される。なお、ＢＬＡＳＴアルゴリズムは、下記の非特許文献３に記載されている。
[3] P. W. Wolniansky and al., “V-BLAST: An architecture for realizing very high data rates over the rich scattering wireless channel”, in Proc. Int. Symposium on Advanced Radio Technologies, Boulder, CO, Sept. 1998.
＜ＭＭＳＥ−ＢＬＡＳＴアルゴリズム＞
ＭＭＳＥ−ＢＬＡＳＴアルゴリズムの原理は、ＺＦ−ＢＬＡＳＴアルゴリズムと同じである。ただし、ＭＭＳＥ−ＢＬＡＳＴアルゴリズムにおいては、上述の（５）式の代わりに下記（１１）式が使用される。なお、「α」はノイズの分散値を表す。
Figure 2006067070
 ＜最尤復号（Maximum Likelihood Decoding）アルゴリズム＞
最尤復号アルゴリズムでは、受信信号ｙと可能なすべての送信ベクトル（具体的には、チャネルＨにより変化した全送信ベクトル）とを比較し、最尤原理により送信データを推定する。この方法は、下記（１２）式で表すことができる。
Figure 2006067070
 最尤復号アルゴリズムは、送信データを再生するための最適な方法の１つである。しかし、このアルゴリズムでは、受信信号から推定される送信ベクトルと、存在し得るすべての送信ベクトルとを比較する必要があるので、計算量が大きくなるという問題がある。例えば、変調方式がＢＰＳＫであるものとすると、送信アンテナが４本のシステムにおいては１６（＝２4）種類の送信ベクトルが存在し、送信アンテナが５本のシステムにおいては３２（＝２5）種類の送信ベクトルが存在する。すなわち、受信信号から推定される送信ベクトルと比較すべきベクトルの数は、送信アンテナの数の増加に対して指数関数的に増加することとなる。
特開２００２−４４０５１号公報
ＭＩＭＯシステムに上述の５つの従来技術を適用すると、以下の問題が生じる。すなわち、Zero ForcingアルゴリズムおよびＭＭＥＳアルゴリズムでは、送信シンボルの推定精度が低く、その結果、復号性能が低くなる。
ＢＬＡＳＴアルゴリズムおよび最尤復号アルゴリズムは、送信シンボルの推定精度は高くなる。特に、最尤復号アルゴリズムは、復号性能に注目すれば最適なアルゴリズムである。しかし、これらのアルゴリズムは、計算が複雑であり、受信装置の負荷が重くなる。例えば、最尤復号アルゴリズムの計算量は、状態数（送信アンテナの数および変調方式に依存する）に対して指数関数的に大きくなる。
本発明の目的は、複数のアンテナを利用して信号を受信する受信装置において、復号処理の複雑さを抑えつつ良好な復号性能を得ることである。
本発明の受信装置は、送信装置からの信号を伝送するチャネルの状態に基づいて複数のアンテナを利用して受信した信号から送信シンボルを推定する推定手段と、上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、上記確からしさを表す情報を上記推定手段に与えながら上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有する。そして、上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する。
上記受信装置においては、推定した送信シンボルの確からしさを表す情報を利用して新たに送信シンボルを推定する処理が繰り返し実行される。このため、繰返し回数を増やせば送信シンボルの推定精度は高くなる。
上記推定手段は、最小平均二乗誤差アルゴリズムで送信シンボルを推定するようにしてもよい。この場合、１回当たりの推定処理の計算量は比較的少ないので、その推定処理を繰り返し実行しても受信装置の負荷はさほど重くならない。
本発明によれば、複数のアンテナを利用して信号を受信する受信装置において、復号処理の複雑さを抑えたアルゴリズムで良好な復号性能を得ることができる。
図１は、本発明の実施形態に係るＭＩＭＯシステムの構成を示す図である。ＭＩＭＯシステム１は、複数の送信アンテナを備える送信装置２および複数の受信アンテナを備える受信装置３から構成される。
送信装置２は、Ｎt 本の送信アンテナを利用して複数のシンボルを並列に送信する。図１に示す例では、シンボル０〜シンボルＮt −１が並列に送信される。なお、並列に送信される複数のシンボルは、１つのデータ列から生成されてもよいし、互いに異なるデータ列から生成されてもよい。また、各送信シンボルは、例えば、変調方式がＢＰＳＫである場合は１ビットにデータに対応し、ＱＰＳＫである場合は２ビットのデータに対応する。
受信装置３は、Ｎr 本の受信アンテナを利用して送信装置２から送信された信号を受信する。ここで、受信アンテナの数は、送信装置２に設けられている送信アンテナの数以上であるものとする。すなわち、「Ｎr >=Ｎt 」であるものとする。また、受信装置３は、複数の受信アンテナを介して受信した信号から送信シンボルを再生する信号再生部４を備えている。
ＭＩＭＯシステム１において、各送信アンテナから送信された信号は、各受信アンテナにより受信される。すなわち、例えば、送信アンテナ（０）から送信された信号は、受信アンテナ（０）〜（Ｎr −１）により受信され、同様に、送信アンテナ（Ｎt −１）から送信された信号も、受信アンテナ（０）〜（Ｎr −１）により受信される。このとき、第ｉ番目の送信アンテナから第ｊ番目の受信アンテナへ信号を伝送するチャネルの状態（または、特性）を、「Ｈij」と表すことにする。そうすると、送信装置２と受信装置３との間の複数のチャネルは、Ｎt ×Ｎr 行列で表すことができる。以下では、この行列のことを「チャネル行列Ｈ」と呼ぶことがある。
受信装置３は、通信に先だってチャネル行列Ｈの各成分を検出する。受信装置３は、例えば、以下の手順でチャネル行列Ｈの各成分を検出する。送信装置２は、まず、送信アンテナ（０）を利用してパイロット信号を送信する。受信装置３は、送信アンテナ（０）から送信されたパイロット信号を受信アンテナ（０）〜（Ｎr −１）を介して受信する。このとき、パイロット信号の送信電力は予め決められているので、受信装置３は、受信アンテナ（０）〜（Ｎr −１）を介して受信した信号の電力やＳＮＲ等をモニタすることにより、チャネル行列の第１列目の各成分（Ｈ0,0〜Ｈ0,Nr-1）を得る。以下、同様に、各送信アンテナから送信されるパイロット信号を利用して、チャネル行列の２〜Ｎt 列の各成分を得ることができる。
上記ＭＩＭＯシステム１において、送信装置２からシンボルｘ（ｘ0〜ｘNt-1）が送信されると、受信装置３において検出される信号ｙ（ｙ0〜ｙNr-1）は、上述した（１）式により表される。したがって、ノイズｎがゼロであるものとすると、受信装置３は、チャネル行列Ｈを検出することにより、正しい送信シンボルを得ることができる。すなわち、受信装置３は、チャネル行列Ｈの各成分による信号への影響を逆算することにより、受信信号ｙから送信シンボルｘを得ることができる。以下では、この逆算処理を含む概念のことを「補償」又は「チャネル補償」と呼ぶことがある。しかし、一般に、ノイズｎはゼロではない。また、チャネル行列Ｈを完全に正確に求めることは出来ない。よって、受信装置３は、受信信号ｙおよびチャネル行列Ｈから送信シンボルを推定し、その推定値の誤差を最小にするためのアルゴリズムを導入する。
図２は、受信装置３が備える信号再生部４について説明する図である。信号再生部４には、受信信号（ｙ：ｙ0〜ｙNr-1）、チャネル状態情報（ＣＳＩ：Channel State Information）、ノイズの分散値が与えられる。ここで、受信信号は、図１に示す受信アンテナ（０）〜（Ｎr −１）を介してそれぞれ受信した信号である。また、チャネル状態情報は、図１を参照しながら説明したチャネル行列Ｈの各成分である。さらに、ノイズの分散値は、例えば、「１／ＳＮＲ」であり、復調に先だって予め求めておくものとする。
信号再生部４は、受信信号、チャネル状態情報、ノイズの分散値、前回の推定値の確からしさ(尤度)に基づいて、各送信シンボルを推定すると共に、その推定値の確からしさ（尤度）を更新する。推定値の確からしさは、この例では、対数尤度比（ＬＬＲ：Log Likelihood Ratios）により表される。
図３および図４は、信号再生部４の処理を示すフローチャートである。ただし、図３および図４は、同じ動作を記述している。
実施形態の信号再生アルゴリズムは、図３または図４に示すステップＳ１〜Ｓ６を２回以上繰り返し実行することにより、送信シンボルを推定する。なお、これらのフローチャートにおいては、チャネル行列Ｈおよびノイズの分散値αは、予め検出されているものとする。また、以降の説明において、「（ｉ）」は第ｉ回目の繰返しサイクルを意味する。「ｋ」は送信アンテナを識別する。「＾」印は推定値を意味する。
ステップＳ１では、後述するステップＳ６で得られる対数尤度比ＬＬＲを利用して、送信シンボルの平均値を計算する。ただし、第１回目の繰返しサイクルでは、送信シンボルの平均値は、予め決められた初期値が使用される。例えば、変調方式がＢＰＳＫである場合は、送信シンボルは「＋１」又は「−１」なので、その平均値の初期値として「０」が設定される。
ステップＳ２では、ステップＳ１で得られた送信シンボルの平均値を利用して、送信シンボルの分散が計算される。ただし、第１回目の繰返しサイクルでは、送信シンボルの分散も、予め決められた初期値が使用される。送信シンボルの分散の初期値は、変調方式に応じて設定される。例えば、変調方式がＢＰＳＫである場合は、その分散の初期値として「１」が設定される。
ステップＳ３では、ステップＳ２で得られた送信シンボルの分散、チャネル行列Ｈ、ノイズの分散値αを利用して、ＭＭＳＥ係数が計算される。ＭＭＳＥ係数は、最小平均二乗誤差アルゴリズムで受信信号から送信シンボルを推定するために使用する行列（または、その成分）であり、チャネル補償のための係数である。
ステップＳ４では、ステップＳ３で得られたＭＭＳＥ係数を利用して、受信信号から送信シンボルを再生する。
ステップＳ５では、ステップＳ４で得た推定値を利用して対数尤度比ＬＬＲの変量ΔＬを計算する。このとき、誤差はガウス分布であるものとする。
ステップＳ６では、前回の繰返しサイクルで計算された対数尤度比ＬＬＲを、変量ΔＬを利用して更新する。ここで、対数尤度比ＬＬＲは、送信シンボルの推定値の確からしさを表すパラメータである。なお、第ｉ回目の繰返しサイクルにおいて新たに得られる対数尤度比ＬＬＲは、第ｉ＋１回目のサイクルで使用される。
実施形態の信号再生アルゴリズムは、上述のステップＳ１〜Ｓ６を２回以上繰り返し実行する。そして、繰返しサイクル毎に、対数尤度比ＬＬＲ、送信シンボルの平均値、送信シンボルの分散が更新され、それらの更新結果に従って送信シンボルが推定される。したがって、送信シンボルの推定精度は、ステップＳ１〜Ｓ６の処理の繰返し回数が多いほど高くなる。
繰返し回数は、特に限定されるものではないが、シミュレーション等により所望の復号特性が得られるように決定されるようにしてもよい。あるいは、復号エラーの発生率をモニタし、その値が所定の閾値よりも小さくなった時点で繰返し処理を停止するようにしてもよい。
なお、図３または図４に示すフローチャートにおいて、ステップＳ１〜Ｓ４の処理は、既存の最小平均二乗誤差アルゴリズムと同じである。ただし、既存の最小平均二乗誤差アルゴリズムにおいては、送信シンボルの推定値から対数尤度比ＬＬＲを計算する処理、および対数尤度比ＬＬＲを利用して送信シンボルを繰り返し推定する手順は存在しない。すなわち、実施形態のアルゴリズムは、対数尤度比ＬＬＲを利用して送信シンボルを繰り返し推定することにより、その推定精度を高めたことを１つの特徴とする。
次に、実施形態のアルゴリズムの具体的な計算方法を詳しく説明する。なお、以下の説明では、「Ｅ｛ａ｝」は「ａ」の平均値を意味する。また、「＾」印が付された「ｘk」は、ｋ番目の送信アンテナから送信されたシンボルの推定値を表す。
最小平均二乗誤差（ＭＭＳＥ：Minimum Mean Square Error）アルゴリズムでは、一般に、コスト関数を最小にする解が得られる。実施形態の信号再生処理に最小平均二乗誤差アルゴリズムを適用すると、コスト関数は下記の通りである。
Figure 2006067070
 このコスト関数を最小化して最小平均二乗誤差アルゴリズムの解を得ることは、下記の（１３）式および（１４）式と等価である。ここで、（１３）式は、ｋ番目の送信アンテナから送信されるシンボルについて記述し、（１４）式は、ＭＩＭＯシステム全体について記述している。
Figure 2006067070
 なお、「Ｃｏｖ」は、共分散（Covariance）を計算する演算子である。例えば、確率変数Ｘ１、Ｘ２の平均値をμ１、μ２とするとき、この演算子は下記を意味する。
Ｃｏｖ（Ｘ１，Ｘ２）＝Ｅ｛（Ｘ１−μ１）（Ｘ２−μ２）｝
また、最小平均二乗誤差アルゴリズムにおいて、（１３）式または（１４）式を得るための手順は、例えば、下記の非特許文献４に記載されている。
[4] M. Tuchler and al., “Minimum mean square error equalization using a priori information”, IEEE Trans. Comm., March 2002.
（１４）式に含まれている「Ｃｏｖ（ｘ，ｙ）」および「Ｃｏｖ（ｙ，ｙ）」は、上述の（１）式を利用して単純化すると、それぞれ下記（１５）式および（１６）式で記述される。
Figure 2006067070
 さらに、下記（１７）式に示す行列Ｖを定義する。行列Ｖは、「Ｃｏｖ（ｘ，ｘ）」により得られる対角行列である。よって、この行列Ｖの対角線上の各成分は、それぞれ、対応する送信アンテナから送信されるシンボルの分散に相当する。
Figure 2006067070
 そして、（１５）式および（１６）式を（１４）式に代入すると、下記の（１８）式が得られる。
Figure 2006067070
 このように、実施形態のアルゴリズムによれば、送信シンボルは、送信シンボルの平均値、行列Ｇ、受信信号ｙ、チャネル行列Ｈを利用して推定される。ここで、行列Ｇは、送信シンボルの分散を含んでいる。すなわち、送信シンボルの推定は、各送信シンボルの平均および分散に依存する。よって、実施形態のアルゴリズムでは、この推定のためにガウス分布を導入する。そうすると、条件付き確率は、下記（２０）式で記述される。このとき、「μk」は、ｋ番目の送信アンテナから送信されたシンボルの推定値の平均値である。また、「σk 2」は、ｋ番目の送信アンテナから送信されたシンボルの推定値の分散である。
Figure 2006067070
 また、送信シンボルの対数尤度比ＬＬＲは、（２０）式の条件付き確率を利用して下記（２２）で定義される。ここで、「ΔＬ（ｘk）」は、ｋ番目の送信アンテナから送信されるシンボルが「＋１」であるときの確率と「−１」であるときの確率の比の対数により定義され、あるサイクルにおける対数尤度比ＬＬＲとその次のサイクルにおける対数尤度比ＬＬＲとの変量に相当する。
Figure 2006067070
 そうすると、ｉ＋１回目のサイクルにおける対数尤度比ＬＬＲは、下記（２３）式に示すように、ｉ回目のサイクルにおける対数尤度比ＬＬＲを「ΔＬ」で更新することにより得られる。
Figure 2006067070
 ここで、実施形態のＭＩＭＯシステムの変調方式がＢＰＳＫであるものとする。そうすると、送信シンボルの取り得る値は「＋１」および「−１」なので、下記（２４）式の関係が得られる。
Figure 2006067070
 そして、（２０）式および（２４）式を（２２）式に代入すると、下記（２５）式が得られる。
Figure 2006067070
 さらに、（２５）式を計算すると、対数尤度比ＬＬＲの変量として下記（２６）式が得られる。なお、「Ｒｅ（）」は、かっこ内に記述される虚数の実部を取り出す演算子を意味している。
Figure 2006067070
 ここで、ＭＩＭＯシステムの変調方式がＢＰＳＫであるものとする。そして、ｋ番目の送信アンテナから送信されるシンボルの平均値および分散が、それぞれ「０」「１」であるものとする。そうすると、（２１）式は、下記（２７）で記述される。さらに、この場合、対数尤度比ＬＬＲの変量ΔＬは、下記（２８）式で記述される。なお「ｈjk」及び「ｈkj」は、チャネル行列Ｈの成分である。また、「ｇjk」及び「ｇkj」は、（１９）式で表される行列Ｇの成分である。
Figure 2006067070
 ところで、ｉ回目のサイクルにおける送信シンボルｘの平均値は、ＭＩＭＯシステムの変調方式がＢＰＳＫであるものとすると、下式で定義することができる。
Figure 2006067070
 そうすると、この定義式に上述の（２２）式を代入することにより、下記（２９）式が得られる。さらに、この（２９）式から（３０）式を得ることができる。
Figure 2006067070
 図５は、実施形態のアルゴリズムをまとめて記載したものである。なお、各計算式は、（１３）式から（３０）式へ至る過程で登場したものと実質的に同じである。
ここで、図３または図４に示したフローチャートと、上述の計算との対応関係を説明する。
＜第１回目の繰返しサイクル＞
ステップＳ１の平均値およびステップＳ２の分散は、変調方式に応じた初期値として与えられる。ステップＳ３のＭＭＳＥ係数（行列Ｇ）は、（１９）式により計算される。ここで、ＭＭＳＥ係数を計算するためのチャネル行列Ｈ、ノイズの分散値αは、予め検出されている。また、送信シンボルの分散は、ステップＳ２で得られている。
ステップＳ４では、（１８）式により送信シンボルの推定値が計算される。受信信号ｙは、複数の受信アンテナを介して受信する。また、ＭＭＳＥ係数（行列Ｇ）は、ステップＳ３で得られた値を使用する。
ステップＳ５では、対数尤度比ＬＬＲの変量ΔＬが計算される。この変量ΔＬは、一般的には、（２２）式により得られ、ＭＩＭＯシステムの変調方式がＢＰＳＫである場合には、（２６）式により得られ、さらに送信シンボルの平均値および分散としてそれぞれ「０」「１」が与えられた場合には、（２８）式により得られる。
ステップＳ６では、（２３）式によって対数尤度比ＬＬＲが計算される。ただし、第１回目の繰返しサイクルにおいて得られた対数尤度比ＬＬＲは、第２回目の繰返しサイクルにおいて使用される。
＜第２回目の繰返しサイクル＞
ステップＳ１において、送信シンボルの平均値は、例えばＭＩＭＯシステムの変調方式がＢＰＳＫである場合は、（２９）式により得られる。このとき、対数尤度比ＬＬＲは、第１回目の繰返しサイクルのステップＳ６で得られた値を使用する。続いて、ステップＳ２において、送信シンボルの分散は、例えばＭＩＭＯシステムの変調方式がＢＰＳＫである場合は、（３０）式により得られる。
ステップＳ３〜Ｓ６の処理は、第１回目の繰返しサイクルと同じである。ただし、第２回目の繰返しサイクルにおいては、第３回目の繰返しサイクルにおいて使用すべき対数尤度比ＬＬＲが計算される。
＜第３回目以降の繰返しサイクル＞
基本的に、第２回目の繰返しサイクルと同じである。なお、繰返し回数は、例えば、シミュレーション等により予め決められていている。或いは、所定の条件が満たされた時点で処理を終了するようにしてもよい。そして、最終的に得られる推定値が、再生された送信シンボルとして出力される。
次に、本発明の実施形態のアルゴリズムと従来技術とを比較する。以下では、復号性能およびアルゴリズムの複雑さについて比較する。
＜復号性能＞
図６および図７は、実施形態のアルゴリズムおよび従来技術のエラー率についてのシミュレーション結果を示す図である。なお、図６は、送信アンテナおよび受信アンテナがそれぞれ４本である場合の結果であり、図７は、送信アンテナおよび受信アンテナがそれぞれ８本である場合の結果である。また、シミュレーションは、以下の条件に基づいて行った。
キャリア周波数：２．４ＧＨｚ
帯域幅：２０ＭＨｚ
変調方式：ＢＰＳＫ
符号：畳込み符号
符号化率：１／２
データサブキャリア数：５２（＝４８＋４）
ガードインターバル長：１６サンプル
チャネルモデル：5-path exponential decay
最大遅延拡散：０．２０μ秒
端末速度：０．５ｍ／秒
サンプル周期：０．５μ秒
本発明の実施形態のアルゴリズムは、Zero Forcingアルゴリズム、ＭＭＥＳアルゴリズム、ＺＦ−ＢＬＡＳＴアルゴリズムと比較して良好な受信性能が得られる。ただし、ＭＭＳＥ−ＢＬＡＳＴアルゴリズムと比較すると、やや劣っている。
また、本発明の実施形態のアルゴリズムでは、繰返し回数が２回であるときよりも３回のときの方がエラー率は改善する。なお、繰返し回数が１回であるときの特性は、実質的に、ＭＭＳＥアルゴリズムと同じになる。
＜アルゴリズムの複雑さ＞
ここでは、本発明の実施形態のアルゴリズムとＢＬＡＳＴとを比較する。また、送信シンボルを推定するために実行される演算回数について検討する。さらに、以下で検討するＭＩＭＯシステムは、送信アンテナ数＝受信アンテナ数＝Ｎであるものとする。
ＢＬＡＳＴの演算回数は、下記のように表すことができる。なお、「Ｆ（Ｎ3）」は、「Ｎ」についての３次の関数である。
掛け算回数＝４Ｎ4＋Ｆ（Ｎ3）
足し算回数＝４Ｎ4＋Ｆ（Ｎ3）
一方、本発明の実施形態のアルゴリズムにおける演算回数は、下記のように表すことができる。ここで、「Ｋ」は、繰返し回数である。
掛け算回数＝６Ｎ3＋６（Ｋ−１）Ｎ3＝６ＫＮ3
足し算回数＝６Ｎ3＋６（Ｋ−１）Ｎ3＝６ＫＮ3
図８、図９（ａ）、図９（ｂ）は、掛け算回数についての比較結果である。また、図１０、図１１（ａ）、図１１（ｂ）は、足し算回数についての比較結果である。なお、図９（ｂ）は、図９（ａ）の一部を拡大した図であり、図１１（ｂ）は、図１１（ａ）の一部を拡大した図である。なお、ＢＬＡＳＴの掛け算／足し算回数を計算する式において、Ｎが大きくなると、Ｎ4項と比較して「Ｆ（Ｎ3）」の寄与度が小さくなるので、図８〜図１１においては、「Ｆ（Ｎ3）」を無視して計算している。
このように、送信／受信アンテナの数が増加すると、本発明の実施形態のアルゴリズムよりもＢＬＡＳＴの方が、送信シンボルを推定するために必要な計算量が多くなる。上述の例では、送信アンテナの数が４本以上になると、本発明の実施形態のアルゴリズムの方が有利になる。なお、一般に、ＭＩＭＯシステムでは、送信アンテナの数が４本以上である。
上述のように、本発明の実施形態のアルゴリズムによれば、Zero Forcingアルゴリズム、ＭＭＥＳアルゴリズムと比較して良好な受信性能が得られる。また、本発明の実施形態のアルゴリズムは、ＢＬＡＳＴアルゴリズムと比較して送信シンボルを推定するために必要な計算量が少なくなる。なお、本発明の実施形態のアルゴリズムの計算量が最尤復号アルゴリズムと比較して少なくなることは言うまでもない。
実施例１では、本発明のアルゴリズムがＭＩＭＯ−ＯＦＤＭシステムに適用される。ＯＦＤＭについては、例えば、下記の非特許文献５、６に記載されており、ＭＩＭＯ−ＯＦＤＭについては下記の非特許文献７に記載されているが、ここで概要を説明する。
[5] R. Van Nee and R. Prasad, “OFDM For Wireless Multimedia Communications”, Artech House Publishers, 2000
[6] L. J. Cimini, “Analysis and simulation of digital mobile channel using orthogonal frequency division multiple access”, IEEE trans. Comm., pp.665-675, 1995
[7] A. Van Zeist and al., “Implementation of a MIMO OFDM-based wireless LAN system”, IEEE Trans. Signal. Process. 52, No 2, 2004, pp. 483-494
ＯＦＤＭでは、各サブキャリアのレートを低くし、それらの低速サブキャリアの数を多くすることで高速送信を実現する技術である。周波数帯域幅は、複数の小さな領域に分割され、それらがそれぞれサブキャリアとして使用される。複数のサブキャリアは互いに直交している。この特性を得るためには、各サブキャリアの周波数は、シンボル期間の逆数の整数倍離れている必要がある。また、チャネルは、各サブキャリアについて推定されて訂正される必要がある。
＜送信機＞
図１２は、ＭＩＮＯ−ＯＦＤＭ送信機の構成図である。以下、この図を参照しながら送信機の構成および動作を説明する。
高速バイナリデータは、符号化（畳込み符号）、インターリーブ、変調（ＢＰＳＫ、ＱＰＳＫ、１６ＱＡＭなど）される。つづいて、データは、並列低速変調データ列に変換され、Ｎc 個のサブキャリアに載せられる。変調信号は、Ｎ点逆離散フーリエ変換（ＩＤＦＴ）により周波数分割多重される。結果として得られるＯＦＤＭ信号は、Ｄ／Ａコンバータによりアナログ信号に変換され、ＲＦ帯にアップコンバートされて出力される。
ここで、Ｎ個のサブキャリアを利用してＮ個のシンボルを送信するＯＦＤＭシステムを想定する。そうすると、第ｉ番目のＯＦＤＭシンボル期間に送信されるＯＦＤＭ信号は、下記（３１）で表される。ここで、「ｓｉ」は、第ｉ番目の送信シンボル（データ、ゼロパディングを含む）を表す。「Ｎ」は、ブロックサイズを表す。「ｊ」は、「−１」の平方根である。
Figure 2006067070
 ＯＦＤＭでは、ＩＳＩ（Inter Symbol Interference）およびＩＣＩ（Inter channel Interference）の影響を排除するために、各ＯＦＤＭシンボル間にガードインターバルが挿入される。ガードインターバルはcyclic prefixなので、ＯＦＤＭシンボルの遅延波の遅延時間がカードインターバルよりも小さければ、ＯＦＤＭシンボルはその遅延により劣化することはない。ここで、もし、ガードインターバルが出力シーケンスの循環拡張であれば、ガードインターバルが挿入された送信信号は下記（３２）式で表される。尚、「Ｇ」はガードインターバルの長さを意味する。
Figure 2006067070
 伝搬チャネルは、異なる振幅および時間遅延を持ったＰ個の離散パスから構成されるものとすると、ベースバンドインパルス応答は、下記（３３）式で表される。尚、「αp」および「τp」は、それぞれ、ｐ番目の離散パスの複素チャネル利得および時間遅延を意味する。
Figure 2006067070
 さらに、総チャネルパワーを「１」とし、τmax＝maxτj＜Ｎであるものとする。そして、コヒーレント検出が完全（すなわち、タイミングが完全）であれば、１つのＯＦＤＭシンボルのみが干渉する。そうすると、周波数領域において、チャネル転送関数Ｈは、下記（３４）式になる。
Figure 2006067070
 ＜受信機＞
図１３は、ＭＩＮＯ−ＯＦＤＭ受信機の構成図である。以下、この図を参照しながら受信機の構成および動作を説明する。
受信信号は、バンドパスフィルタを通過した後、ダウンコンバートされる。ダウンコンバードされた信号は、Ａ／Ｄコンバータでサンプリング（即ち、デジタル信号に変換）され、ガードインターバルが除去される。サンプリングされたデータは、Ｎ点離散フーリエ変換器（ＩＤＦＴ）で処理された後、サブキャリアチャネルデータに分離される。チャネル歪の影響を減らすために、最小平均二乗誤差またはZero Forcingの係数を利用して、サブキャリア毎にチャネルが周波数領域で補償される。最後に、各サブキャリアのデータが、復調、逆インターリーブ、復号される。
ここで、ガードインターバルが最大遅延拡散よりも大きければ、受信ベースバンド信号は、下記（３５）式で表すことができる。ここで、「ｒ」および「ｎ」は、ガードインターバルが除去された後であって且つ処理前における、ＯＦＤＭシンボルの受信信号および付加された白色ガウス雑音（ＡＷＧＮ：Additive White Gaussian Noise）である。
Figure 2006067070
 この後、復調処理が実行され、情報ビットについて判定が行われる。
上述したように、チャネル変動は、サブキャリア毎に訂正される必要がある。例えば、ＩＥＥＥ８０２．１１ａ、１１ｇにおいては、ＩＦＦＴ／ＦＦＴのディメンジョンは６４であるが、４８個のサブキャリアのみがデータシンボルを送信するために使用される。他の４個のサブキャリアが周波数領域におけるデータフレームの同期のために使用される。よって、１２個のサブキャリアは、受信処理過程でキャンセルされる。この場合、データを送信しているサブキャリアのみについてチャネルの推定および訂正が行われる。
＜本発明のＭＩＭＯ−ＯＦＤＭへの適用＞
送信機は、Ｎt 本の送信アンテナを備えている。また、各送信アンテナに対して互いに独立したチャネルエンコーダが設けられる。そして、各エンコーダの出力は、それぞれ複数のサブキャリアに載せられて、対応する送信アンテナを介して送信される。
受信機は、Ｎr 本の受信アンテナを備える。これらの受信アンテナを介して受信した信号は、それぞれフーリエ変換された後、ＭＩＭＯ信号処理部１１に与えられる。ここで、ＭＩＭＯ信号処理部１１は、本発明の実施形態のアルゴリズムを実行する信号再生部４を含み、予め記述されたプログラムを実行するプロセッサにより実現される。そして、ＭＩＭＯ信号処理部１１の出力は、Ｎt 個の独立したデータ列であり、各データ列について独立して復調および復号が行われる。
実施例２では、本発明のアルゴリズムがマルチキャリアＣＤＭＡシステムにおいて使用される。マルチキャリアＣＤＭＡ技術は、マルチアクセスのためのＯＦＤＭ変調に取り込まれている。なお、マルチキャリアＣＤＭＡについては、例えば、下記の非特許文献８、９に記載されている。
[8] Yee N., J. P. Linnartz and G. Fettweis, “Multi-carrier CDMA in indoor wireless radio networks”, IEICE Trans. Comm., E77-B pp: 900-904, July 1994
[9] S. Abeta, H. Atarashi and M. Sawahasgi, “Performance of Coherent Multi-Carrier/DS-CDMA and MC-CDMA for Broadband Packet Wireless Access”, IEICE Trans. Commun., vol. E84-B, No.3 March 2001
＜送信機＞
図１４は、マルチキャリアＣＤＭＡ送信機の構成図である。以下、この図を参照しながら送信機の構成および動作を説明する。
マルチキャリアＣＤＭＡ信号は、以下のようにして生成される。すなわち、バイナリデータが符号化（畳込み符号またはターボ符号）、変調（ＢＰＳＫ、ＱＰＳＫ、１６ＱＡＭなど）される。Ｎpl個のデータシンボルが同時に処理される。データシンボルは、ＳＦ個のコピーが作成される。この場合、データは、Ｎpl×ＳＦ個のサブキャリアを利用して送信される。各サブキャリア上で、データに拡散コードが乗算される。拡散コードは、長さがＳＦチップの直交符号である。
送信機は、拡散コードを利用して信号を周波数領域で拡散させる。その後、逆高速フーリエ変換（ＩＦＦＴ）が実行される。ＩＦＦＴのディメンジョンは、「Ｎpl×ＳＦ」である。パラレルデータをシリアルデータに変換した後、シンボル間干渉を回避するためのガードインターバルが付加される。
＜受信機＞
図１５は、マルチキャリアＣＤＭＡ受信機の構成図である。以下、この図を参照しながら受信機の構成および動作を説明する。
受信信号は、ガードインターバルが除去され、さらにパラレルデータに変換されて高速フーリエ変換器（ＦＦＴ）に与えられる。パイロット信号は、サブキャリア毎のチャネル推定のために使用できる。チャネルの補償は、高速フーリエ変換器（ＦＦＴ）の出力で行われる。他の結合方法として、ＯＲＣ（Orthogonal restoring combining）、ＥＧＣ（Equal gain combining）、ＭＲＣ（Maximum ratio combining）、ＭＭＳＥＣ（Minimum mean square error combining）などが適用される。この後、信号から対応するユーザのデータが抽出される。データの抽出は、拡散コード間の直交性が利用される。この処理は、サブキャリア毎に行われる。そして、周波数領域において、ＳＦチップ間隔でデータがコヒーレント的に加算される。
なお、ＯＦＤＭ送信形態は、マルチキャリアＣＤＭＡの特殊ケース（すなわち、ＳＦ＝１）として考えることができる。
＜本発明のマルチキャリアＣＤＭＡへの適用＞
送信機は、Ｎt 本の送信アンテナを備えている。各送信アンテナを介して送信される信号は、それぞれ複数のサブキャリアを含んでおり、各サブキャリアの信号について拡散処理が行われる。
受信機は、Ｎr 本の受信アンテナを備える。本発明のアルゴリズムは、サブキャリア毎に、ＭＩＭＯ検出部１２により逆拡散処理の前に実行される。ＭＩＭＯ検出部１２は、信号再生部４を含み、予め記述されたプログラムを実行するプロセッサにより実現される。ここで、ＭＩＭＯ検出部１２を設ける位置は、公知の構成と同じであり、従来の受信機においてＭＭＳＥアルゴリズムまたはＢＬＡＳＴアルゴリズム等を実行する処理部が設けられている位置と同じである。
実施例３では、本発明のアルゴリズムは、単一キャリア通信システムにおいて時間領域イコライザとして使用される。この場合、送信機は、図１６（ａ）に示すように、入力データからＮt 個のデータ列を生成するシリアル／パラレル変換器、各データ列をそれぞれ符号化するＮt 個のエンコーダ、各データ列をそれぞれ変調するＮt 個の変調器、そして各変調信号を送信するＮt 本の送信アンテナを備える。
受信機は、図１６（ｂ）に示すように、Ｎr 本の受信アンテナ、ＭＩＭＯ処理部１３、ＭＩＭＯ処理部１３の出力を復調するＮt 個の復調器、復調信号を復号するＮt 個のデコーダを備える。そして、ＭＩＭＯ処理部１３は、本発明のアルゴリズムを実行することにより時間領域イコライザとして動作し、Ｎr 個の受信信号からＮt 個のシンボルの推定値を出力する。
実施例４では、本発明のアルゴリズムは、単一キャリア通信システムにおいて周波数領域イコライザとして使用される。この場合、送信機は、図１７（ａ）に示すように、第３の実施例の送信機にガードインターバルＧＩを挿入する機能を追加することにより得られる。なお、ガードインターバルは、Ｎデータ毎に挿入される。
受信機では、図１７（ｂ）に示すように、各受信アンテナを介して受信した信号は、ガードインターバルＧＩが除去された後、それぞれＮ個のパラレルデータに変換される。すなわち、Ｎr 組のＮパラレルデータが生成される。そして、各Ｎパラレルデータは、それぞれ高速フーリエ変換器（ＦＦＴ）により周波数領域の信号に変換されてＭＩＭＯ検出部１４に与えられる。
ＭＩＭＯ検出部１４は、本発明のアルゴリズムを実行することにより周波数領域イコライザとして動作し、Ｎr 組のＮパラレルデータからＮt 組のＮパラレルデータの推定値を出力する。そして、ＭＩＭＯ検出部１４の出力は、逆高速フーリエ変換器（ＩＦＦＴ）により時間領域の信号に変換された後、復調および復号される。
実施例５では、本発明のアルゴリズムは、単一キャリアＣＤＭＡに適用される。この場合、送信機は、図１８（ａ）に示すように、第３の実施例の送信機に拡散コードを乗算する拡散器を設けることにより得られる。なお、送信信号は、拡散コードを乗算することにより、時間領域で拡散される。また、受信機は、図１８（ｂ）に示すように、第３の実施例の受信機に拡散コードを乗算する逆拡散器を設けることにより得られる。
本発明は、上述の実施例１〜５に限定されるものではなく、下記の実施形態も可能である。
（１）本発明の受信装置は、空間キャンセラを備えていてもよい。空間キャンセラは、特定の方向から到来する無線信号をキャンセルする回路である。なお、本発明のアルゴリズムを実行する信号再生部４は、例えば、空間キャンセラの出力から送信シンボルを推定するようにしてもよい。或いは、信号再生部４の出力が空間キャンセラに与えられるようにしてもよい。この構成によれば、不要な無線信号の影響を抑制できるので、受信特性がさらに改善する。
（２）本発明の受信装置は、最尤復号器（ＭＬＤ：Maximum Likelihood Decoder）を備えていてもよい。この場合、例えば、本発明のアルゴリズムを実行する信号再生部４の出力が最尤復号器に与えられる。この構成によれば、受信特性がさらに改善する。
（３）上述の実施例は、ＭＩＭＯシステムを前提として記載したが、本発明はこれに限定されるものではなく、１本の送信アンテナから送信された信号を複数の受信アンテナで受信するシステムにも適用可能である。また、上述の実施例は、１本の送信アンテナから送信された信号を１本の受信アンテナで受信するシステムに適用する形態を排除するものではない。
（４）上述の実施例では、送信シンボルを推定するアルゴリズムとしてＭＭＳＥを導入しているが、本発明はこの構成に限定されるものではない。すなわち、本発明の１つの特徴は、送信シンボルの推定結果の確からしさを示す情報を次回の推定処理に利用しながら推定処理を繰り返すことであり、推定アルゴリズム自体はＭＭＳＥに限定されるものではない。よって、ＭＭＳＥアルゴリズムの代わりに他のアルゴリズム（例えば、Zero Forcing アルゴリズム、ＢＬＡＳＴアルゴリズム、最尤復号アルゴリズム等）で送信シンボルを推定してもよい。
（５）上述の実施例では、ＭＩＭＯシステムの変調方式がＢＰＳＫであるものとして説明したが、本発明のアルゴリズムを導入するうえで、変調方式は特に限定されるものではない。
（付記１）送信装置からの信号を伝送するチャネルの状態に基づいて、複数のアンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
ことを特徴とする受信装置。
（付記２）付記１に記載の受信装置であって、
上記推定手段は、
上記確からしさを表す情報に基づいて、上記チャネルを介した伝送による信号への影響を補償するための補償係数を導出する手段と、
上記補償係数を利用して上記受信信号から送信シンボルを推定する手段、
を有する。
（付記３）付記１に記載の受信装置であって、
上記推定手段は、
上記確からしさを表す情報に基づいて送信シンボルの平均値を計算する手段と、
その平均値に基づいて送信シンボルの分散を計算する手段と、
その分散に基づいて上記チャネルを介した伝送による信号への影響を補償するための補償係数を導出する手段と、
上記補償係数を利用して上記受信信号から送信シンボルを推定する手段、
を有する。
（付記４）付記１に記載の受信装置であって、
上記推定手段は、最小平均二乗誤差アルゴリズムで複数のアンテナを利用して受信した信号から送信シンボルを推定する。
（付記５）付記１に記載の受信装置であって、
上記尤度計算手段は、繰返しサイクル毎に、前回の繰返しサイクルで得られた確からしさを表す情報を更新する。
（付記６）付記１に記載の受信装置であって、
上記尤度計算手段は、推定された送信シンボルについての対数尤度比を計算する。
（付記７）付記６に記載の受信装置であって、
上記尤度計算手段は、繰返しサイクル毎に、前回の繰返しサイクルで得られた対数尤度比を更新する。
（付記８）付記１に記載の受信装置であって、
単一キャリア通信システムにおいて時間領域イコライザとして使用される。
（付記９）付記１に記載の受信装置であって、
単一キャリア通信システムにおいて周波数領域イコライザとして使用される。
（付記１０）付記１に記載の受信装置であって、
単一キャリアＣＤＭＡにおいて使用される。
（付記１１）付記１に記載の受信装置であって、
ＯＦＤＭシステムにおいて使用される。
（付記１２）付記１に記載の受信装置であって、
マルチキャリアＣＤＭＡシステムにおいて使用される。
（付記１３）付記１に記載の受信装置であって、
さらに空間キャンセラを備え、
上記推定手段は、上記空間キャンセラにより処理された受信信号から送信シンボルを推定する。
（付記１４）付記１に記載の受信装置であって、
上記推定手段による推定結果について最尤復号を行う最尤復号回路をさらに備える。
（付記１５）複数のアンテナを利用して信号を受信する受信装置であって、
送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数のアンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
ことを特徴とする受信装置。
（付記１６）複数の送信アンテナを備える送信装置および複数の受信アンテナを備える受信装置から構成されるＭＩＭＯシステムであって、
上記受信装置は、
上記送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数の受信アンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
ことを特徴とするＭＩＭＯシステム。
（付記１７）ＭＩＭＯシステムの受信装置において送信シンボルを推定する方法であって、
送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数のアンテナを利用して受信した信号から送信シンボルを推定し、
上記推定手段により推定された送信シンボルの確からしさを計算し、
上記確からしさを表す情報を上記推定手段に与えながら上記推定処理および計算処理を繰り返し実行させ、
２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
ことを特徴とする送信シンボル推定方法。
本発明の実施形態に係るＭＩＭＯシステムの構成を示す図である。 受信装置が備える信号再生部について説明する図である。 信号再生部の処理を示すフローチャート（その１）である。 信号再生部の処理を示すフローチャート（その２）である。 実施形態のアルゴリズムをまとめて記載したものである。 エラー率についてのシミュレーション結果を示す図（その１）である。 エラー率についてのシミュレーション結果を示す図（その２）である。 掛け算回数を比較した表である。 掛け算回数を比較した図である。 足し算回数を比較した表である。 足し算回数を比較した図である。 ＭＩＮＯ−ＯＦＤＭ送信機の構成図である。 ＭＩＮＯ−ＯＦＤＭ受信機の構成図である。 マルチキャリアＣＤＭＡ送信機の構成図である。 マルチキャリアＣＤＭＡ受信機の構成図である。 （ａ）は第３の実施例の送信機、（ｂ）は第３の実施例の受信機である。 （ａ）は第４の実施例の送信機、（ｂ）は第４の実施例の受信機である。 （ａ）は第５の実施例の送信機、（ｂ）は第５の実施例の受信機である。 ＭＩＭＯシステムの基本構成図である。
符号の説明
１ ＭＩＭＯシステム
２ 送信装置
３ 受信装置
４ 信号再生部
１１ ＭＩＭＯ信号処理部
１２、１４ ＭＩＭＯ検出部
１３、１５ ＭＩＭＯ処理部

Claims (10)

  1. 送信装置からの信号を伝送するチャネルの状態に基づいて、複数のアンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
    上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
    上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
    上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
    ことを特徴とする受信装置。
  2. 請求項１に記載の受信装置であって、
    上記推定手段は、
    上記確からしさを表す情報に基づいて、上記チャネルを介した伝送による信号への影響を補償するための補償係数を導出する手段と、
    上記補償係数を利用して上記受信信号から送信シンボルを推定する手段、
    を有する。
  3. 請求項１に記載の受信装置であって、
    上記推定手段は、
    上記確からしさを表す情報に基づいて送信シンボルの平均値を計算する手段と、
    その平均値に基づいて送信シンボルの分散を計算する手段と、
    その分散に基づいて上記チャネルを介した伝送による信号への影響を補償するための補償係数を導出する手段と、
    上記補償係数を利用して上記受信信号から送信シンボルを推定する手段、
    を有する。
  4. 請求項１に記載の受信装置であって、
    上記推定手段は、最小平均二乗誤差アルゴリズムで複数のアンテナを利用して受信した信号から送信シンボルを推定する。
  5. 請求項１に記載の受信装置であって、
    上記尤度計算手段は、推定された送信シンボルについての対数尤度比を計算する。
  6. 請求項１に記載の受信装置であって、
    ＯＦＤＭシステムにおいて使用される。
  7. 請求項１に記載の受信装置であって、
    マルチキャリアＣＤＭＡシステムにおいて使用される。
  8. 複数のアンテナを利用して信号を受信する受信装置であって、
    送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数のアンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
    上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
    上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
    上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
    ことを特徴とする受信装置。
  9. 複数の送信アンテナを備える送信装置および複数の受信アンテナを備える受信装置から構成されるＭＩＭＯシステムであって、
    上記受信装置は、
    上記送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数の受信アンテナを利用して受信した信号から送信シンボルを推定する推定手段と、
    上記推定手段により推定された送信シンボルの確からしさを計算する尤度計算手段と、
    上記確からしさを表す情報を上記推定手段に与え、上記推定手段および尤度計算手段に上記処理を繰り返し実行させる繰返し手段、を有し、
    上記推定手段は、２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
    ことを特徴とするＭＩＭＯシステム。
  10. ＭＩＭＯシステムの受信装置において送信シンボルを推定する方法であって、
    送信装置からの信号を伝送するチャネルの状態に基づいて、上記複数のアンテナを利用して受信した信号から送信シンボルを推定し、
    上記推定手段により推定された送信シンボルの確からしさを計算し、
    上記確からしさを表す情報を上記推定手段に与えながら上記推定処理および計算処理を繰り返し実行させ、
    ２回目以降の推定処理においては、上記確からしさを表す情報を利用して送信シンボルを推定する
    ことを特徴とする送信シンボル推定方法。
JP2004245098A 2004-08-25 2004-08-25 Ｍｉｍｏシステム受信方法及びその装置 Pending JP2006067070A (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2004245098A JP2006067070A (ja) 2004-08-25 2004-08-25 Ｍｉｍｏシステム受信方法及びその装置
GB0506429A GB2417651B (en) 2004-08-25 2005-03-30 Signal receiving method in MIMO system and device therefor
US11/105,435 US7492830B2 (en) 2004-08-25 2005-04-14 Signal receiving method in MIMO system and device therefor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004245098A JP2006067070A (ja) 2004-08-25 2004-08-25 Ｍｉｍｏシステム受信方法及びその装置

Publications (1)

Publication Number Publication Date
JP2006067070A true JP2006067070A (ja) 2006-03-09

Family

ID=34567606

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004245098A Pending JP2006067070A (ja) 2004-08-25 2004-08-25 Ｍｉｍｏシステム受信方法及びその装置

Country Status (3)

Country Link
US (1) US7492830B2 (ja)
JP (1) JP2006067070A (ja)
GB (1) GB2417651B (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8199844B2 (en) 2008-09-10 2012-06-12 Samsung Electronics Co., Ltd. Apparatus and method for canceling inter-antenna interference of orthogonally coded signals
JP2013229763A (ja) * 2012-04-25 2013-11-07 Nippon Hoso Kyokai <Nhk> 受信装置及びプログラム

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0222555D0 (en) 2002-09-28 2002-11-06 Koninkl Philips Electronics Nv Packet data transmission system
KR100943610B1 (ko) * 2004-07-20 2010-02-24 삼성전자주식회사 다중 시공간 블록 부호화 방식을 사용하는 다중 입력 다중 출력 시스템에서 안테나 셔플링 정보 피드백 장치 및 방법
JP4387282B2 (ja) * 2004-10-20 2009-12-16 株式会社エヌ・ティ・ティ・ドコモ 信号分離装置及び信号分離方法
US7953039B2 (en) * 2005-04-21 2011-05-31 Samsung Elecronics Co., Ltd. System and method for channel estimation in a delay diversity wireless communication system
US7961800B2 (en) * 2005-09-07 2011-06-14 Nec Corporation Adaptive radio/modulation apparatus, receiver apparatus, wireless communication system, and wireless communication method
KR100922936B1 (ko) * 2006-02-07 2009-10-22 삼성전자주식회사 통신 시스템에서 데이터 송수신 방법 및 시스템
JP2007295549A (ja) * 2006-03-31 2007-11-08 Matsushita Electric Ind Co Ltd Ｍｉｍｏ受信装置およびｍｉｍｏ通信システム
WO2008021008A2 (en) * 2006-08-07 2008-02-21 Interdigital Technology Corporation Method, apparatus and system for implementing multi-user virtual multiple-input multiple-output
US8265212B2 (en) 2006-12-01 2012-09-11 Apple Inc. Antenna selection for MIMO decoding
JP4871801B2 (ja) * 2007-07-09 2012-02-08 キヤノン株式会社 通信装置、通信方法、当該通信方法をコンピュータに実行させるためのコンピュータプログラム
KR101454027B1 (ko) * 2007-08-10 2014-10-24 한국전자통신연구원 병렬 구조를 가지는 시분할 다중화 통신 시스템 및 방법
US7778340B2 (en) * 2007-09-06 2010-08-17 Hong Kong Applied Science And Technology Research Institute Co., Ltd. Accurate channel quality indicator for link adaptation of MIMO communication systems
GB0720559D0 (en) 2007-10-19 2007-11-28 Fujitsu Ltd MIMO wireless communication system
US8391421B2 (en) * 2008-06-19 2013-03-05 Simon Fraser University EDA-based detection of communication signals
US8457240B2 (en) * 2008-08-25 2013-06-04 Daniel Lee Methods of selecting signal transmitting, receiving, and/or sensing devices with probabilistic evolutionary algorithms in information conveyance systems
US8605831B1 (en) * 2008-10-15 2013-12-10 Marvell International Ltd. Demodulation that accounts for channel estimation error
EP2461532A2 (en) 2009-07-31 2012-06-06 Hong Seok Seo Method for processing a received signal to which channel state information is applied in an orthogonal frequency-division multiplexing system
US8761317B2 (en) * 2012-04-04 2014-06-24 Telefonaktiebolaget L M Ericsson (Publ) Soft-bit generation and channel estimation for interference cancellation using decoder output
KR20150117155A (ko) * 2014-04-09 2015-10-19 한국전자통신연구원 다중입력 다중출력 통신 시스템의 연판정 검출 방법 및 장치
US20230117788A1 (en) * 2021-10-15 2023-04-20 Cypress Semiconductor Corporation Systems, methods, and devices for channel interference detection in wireless devices

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7106813B1 (en) * 2000-03-16 2006-09-12 Qualcomm, Incorporated Method and apparatus for combined soft-decision based interference cancellation and decoding
US7068628B2 (en) 2000-05-22 2006-06-27 At&T Corp. MIMO OFDM system
US6731700B1 (en) * 2001-01-04 2004-05-04 Comsys Communication & Signal Processing Ltd. Soft decision output generator
US7154936B2 (en) * 2001-12-03 2006-12-26 Qualcomm, Incorporated Iterative detection and decoding for a MIMO-OFDM system
US7386057B2 (en) * 2003-02-20 2008-06-10 Nec Corporation Iterative soft interference cancellation and filtering for spectrally efficient high-speed transmission in MIMO systems
AU2003293212A1 (en) * 2003-04-14 2004-11-19 Bae Systems Information And Electronic Systems Integration Inc. Joint symbol, amplitude, and rate estimator
GB2406759B (en) * 2003-10-02 2006-06-07 Toshiba Res Europ Ltd Signal decoding methods and apparatus
GB2407007B (en) * 2003-10-09 2006-06-28 Toshiba Res Europ Ltd Estimator for mimo receiver

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8199844B2 (en) 2008-09-10 2012-06-12 Samsung Electronics Co., Ltd. Apparatus and method for canceling inter-antenna interference of orthogonally coded signals
JP2013229763A (ja) * 2012-04-25 2013-11-07 Nippon Hoso Kyokai <Nhk> 受信装置及びプログラム

Also Published As

Publication number Publication date
GB0506429D0 (en) 2005-05-04
GB2417651A (en) 2006-03-01
GB2417651B (en) 2007-10-10
US20060045200A1 (en) 2006-03-02
US7492830B2 (en) 2009-02-17

Similar Documents

Publication Publication Date Title
US7492830B2 (en) Signal receiving method in MIMO system and device therefor
KR100624504B1 (ko) Ｍｉｍｏ시스템에서의 스펙트럼 효율 고속 송신을 위한 반복적 소프트 간섭 소거 및 필터링
US7974350B2 (en) Propagation path estimation method and apparatus
US8462899B2 (en) Multi-user detecting method and device of OFDM transmission signal
JP5344121B2 (ja) シングルキャリア伝送方式における無線通信方法および装置
KR100843253B1 (ko) 반복 수신 방법 및 반복 수신기
Haene et al. A real-time 4-stream MIMO-OFDM transceiver: system design, FPGA implementation, and characterization
Zemen et al. Iterative joint time-variant channel estimation and multi-user detection for MC-CDMA
JP4147193B2 (ja) マルチキャリヤ拡散スペクトル信号の受信
JP2006262039A (ja) 伝搬路推定方法及び伝搬路推定装置
KR100849338B1 (ko) 직교주파수분할다중 방식의 이동통신시스템에서시공간-주파수 부호화/복호화 장치 및 방법
US20090141834A1 (en) Receiver and receiving method for rf signals
KR20070114803A (ko) 추정된 복소수 심벌들의 유연한 복조를 위한 방법
KR100866195B1 (ko) 직교주파수분할다중 방식의 이동통신시스템에서 시공간-주파수 부호화/복호화 장치 및 방법
JP2008258899A (ja) 受信装置および受信方法
JP5030312B2 (ja) 受信機、受信方法および集積回路
Linsalata et al. Design of DGT-based linear and non-linear equalizers for GFDM transmission
JP5143533B2 (ja) 受信装置、及び信号処理方法
Yuan et al. Low-complexity iterative detection in multi-user MIMO ISI channels
JP5241437B2 (ja) 受信装置、及び信号処理方法
US8718162B2 (en) Reception device, reception method, and reception program
Chang et al. Cancellation of ICI by Doppler effect in OFDM systems
JP2011139294A (ja) 送信装置および受信装置
Venkataraman et al. An improved detection technique for receiver oriented MIMO-OFDM systems
JP2008211325A (ja) 受信装置、送信装置、無線送受信システム及び無線受信方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061127

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090310

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090623

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20091020