(57)【特許請求の範囲】
【請求項１】
時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、分解前の条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象である送信信号ベクトルを含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いることを特徴とする空間多重信号検出方法。
【請求項２】
前記複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）及び前記分解可能な条件付確率を複数のグループに分割する処理において、各グループの条件付確率の条件となる事象を含むグループの条件付確率を簡易推定する処理（オーダリング）を含み、該簡易推定する処理によって、前記グループ間の順序付けを前記グループの条件付確率が高く簡易推定される順序で処理可能なように分割することを特徴とする請求項１記載の空間多重信号検出方法。
【請求項３】
前記複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）としてＱＲ分解を用いたことを特徴とする請求項１または請求項２記載の空間多重信号検出方法。
【請求項４】
前記複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）としてブロック三角化分解を用いたことを特徴とする請求項１または請求項２記載の空間多重信号検出方法。
【請求項５】
前記複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）として三重対角化手法を用いたことを特徴とする請求項１または請求項２項記載の空間多重信号検出方法。
【請求項６】
前記二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法は、前記半環（セミリング）としてＭＡＸ−ｌｏｇ領域における半環を用い、和を最大値演算（ＭＡＸ）、積を通常の加算として、
【数１】
0004803384.tif 000214
の如く行われることを特徴とする請求項１から請求項５のいずれか記載の空間多重信号検出方法。
【請求項７】
前記二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法は、前記半環（セミリング）としてＬｏｇ領域における半環を用い、和をヤコビアン対数に基づく演算、積を通常の加算として、
【数２】
0004803384.tif 000215
ここで、
【数３】
0004803384.tif 000216
の如く行われることを特徴とする請求項１から請求項５のいずれか記載の空間多重信号検出方法。
【請求項８】
時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象を含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際しては、自グループにおける条件付確率の条件となる事象を含む前記グループ内の最大条件付確率を示す送信系列を条件付確率の条件として計算するか、あるいは予め簡易推定された複数個の事象（送信系列）を条件として複数個計算し、前記グループ間の順序付けにしたがって各グループ内の条件付確率の計算を前段の最大条件付確率を示す送信系列のもとに計算していく処理を有することを特徴とし、
軟判定出力である該送信系列のビット尤度は、請求項６記載の半環（セミリング）を用いて
（１）対象となるビットが推定対象として存在する該グループにおいて対象ビットと前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理並びに該対象ビットに対する排反事象と該前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とを条件とする複数個の条件付確率をメトリックベースで検出する処理、
（２）次段のグループにおいて対象ビットを含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理並びに対象ビットに対する排反事象を含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理からなる処理、
（３）同様に対象となるビットが存在するグループ以降の段におけるグループにおいても該前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件として複数個の条件付確率をメトリックベースで検出する処理並びに対象ビットに対する排反事象を含む（１）の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とそれに伴って検出される前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とを条件とする複数個の条件付確率をメトリックベースで検出する処理、
（４）以降最終段のグループまで（３）を繰り返す処理、
（５）（４）を完了後、対象ビットの対数尤度比を、前記（１）における対象ビットと前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とその対象ビットに対する排反事象と前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とをメトリックベースで引き算し、前記（２）から（４）の対象となるビットが存在するグループ以降の段のグループにおいても前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率と該対象ビットに対する排反事象を含む前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところの対数尤度比として検出する処理、
を用いて計算する処理を含むことを特徴とする空間多重信号検出方法。
【請求項９】
時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象を含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際しては、自グループにおける条件付確率の条件となる事象を含む前記グループ内の最大条件付確率を示す送信系列を条件付確率の条件として計算し、前記グループ間の順序付けにしたがって各グループ内の条件付確率の計算を前段の最大条件付確率を示す送信系列のもとに計算していく処理を有し、前記尤度を最大化する送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いることが特徴であって、最終段における処理が完了した後の復元抽出（リサンプリング）として、対象となるビットが推定対象として計算された前記条件付確率の組合せの集合からメトリックベースの最大尤度を選択する処理、並びに該対象ビットに対する排反事象が推定対象として計算された前記条件確率の組合せの集合からメトリックベースの最大尤度を選択する処理を有し、両メトリックの差分をとって対象ビットの軟判定出力とする処理を含むことを特徴とする空間多重信号検出方法。
【請求項１０】
時空間多重信号分離において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、分解前の条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象である送信信号ベクトルを含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いる軟入力軟出力検出器と軟入力軟出力復号器とを有し、該軟入力軟出力復号器は符号化前の情報ビット系列に対する対数尤度比（以下、ＬＬＲとする）を出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を元に前記軟入力軟出力検出器のアプリオリ（a priori）入力を作ることを特徴とする時空間反復復号器。
【請求項１１】
時空間検出部と復号部よりなる時空間多重信号分離装置において、該時空間検出部は前記軟入力軟出力検出器と前記軟入力軟出力符号器とを含み、該復号部は前記軟入力軟出力復号器を含む形の構成であって、前記軟入力軟出力復号器へのアプリオリ（a  priori）入力は前記軟入軟出符号器出力に差し引く形で作用して前記軟入力軟出力検出器へのアプリオリ（a  priori）入力を形成し、前記軟入力軟出力検出部へのアプリオリ（a  priori）入力は前記軟入力軟出力検出器出力に差し引く形で作用して前記軟入力軟出力復号器へのアプリオリ（a  priori）入力を形成することを特徴とする請求項１０記載の時空間反復復号器。
【請求項１２】
時空間多重信号分離において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、分解前の条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象である送信信号ベクトルを含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いる軟入力軟出力検出器と軟入力軟出力復号器とを有し、該軟入力軟出力復号器は符号化前の情報ビット系列に対する対数尤度比（以下、ＬＬＲとする）を出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を元に前記軟入力軟出力検出器のソフトレプリカ入力を作ることを特徴とする時空間反復復号器。
【請求項１３】
前記軟入力軟出力符号器は送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素であるｑを法とする加算に代わって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1  と対数尤度比ＬＬＲ2とに対して前記ｑを法とする加算結果ａ1＋ａ2（mod  q）を事後値に持つ対数尤度比ＬＬＲを出力することを特徴とする請求項１０から請求項１２のいずれか記載の時空間反復復号器。
【請求項１４】
前記軟入力軟出力符号器は送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素である２を法とする加算に代わって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1と対数尤度比ＬＬＲ2とに対して、
【数４】
0004803384.tif 000217
を出力する手段、あるいはその近似値を出力する手段を含むことを特徴とする請求項１３記載の時空間反復復号器。
【請求項１５】
前記軟入力軟出力符号器は送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素である２を法とする加算に代わって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1と対数尤度比ＬＬＲ2とに対してその絶対値の比較を行い、値の小さい方を選択する手段を含み、その選択結果に対しての極性付けを前記対数尤度比ＬＬＲ1  及び対数尤度比ＬＬＲ2  のＭＳＢ（Most  Significant  Bit）に対する前記２を法とする加算結果に基づいて行うことを特徴とする請求項１４記載の時空間反復復号器。
【請求項１６】
時空間多重信号分離において、軟入力軟出力検出器と軟入力軟出力復号器とを有し、前記軟入力軟出力検出器は請求項６記載のＭａｘ−ｌｏｇ領域における半環（セミリング）を用いた空間多重検出器であって、その軟出力に当たる対数尤度比に対して重み付けを行い、次段へのアプリオリ（a  priori）入力として用いることを特徴とする時空間反復復号器。
【請求項１７】
前記重み付けを０．７５として用いることを特徴とする請求項１６記載の時空間反復復号器。
【請求項１８】
前記重み付け０．７５をシフト加算によって実現したことを特徴とする請求項１７記載の時空間反復復号器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は空間多重信号検出方法及びそれを用いる時空間反復復号器に関し、特に移動通信において空間多重信号検出方法と反復逐次処理を用いた時空間反復復号器に関する。
【背景技術】
【０００２】
移動通信における電波伝搬路では送信アンテナから到来した電波が周囲の地形等に応じて反射や散乱を受け、一群の素波の集まりとなって受信機に到着する。そのため、品質の高い移動通信を実現する上で常に障害となっていたのが、これらの結果から生じるところのフェージング現象である。このフェージングによる劣悪な電波伝搬環境の克服が長年の移動通信技術者における課題であり、色々な対応策が実用化されてきている。【０００３】
近年、このフェージング現象を悪者扱いするのではなく、逆にフェージングを移動通信における電波伝搬に内在する可能性を秘めた環境資源として見直す動きが活発化している（例えば、非特許文献１，２参照）。【０００４】
また、近年、Multi-USER  Diversityと呼ばれるフェージング変動における空間的位置独立性を利用して電波伝搬路に内在する環境資源を活用する動きもあり、これも同様の動向の一つとも言える。【０００５】
上記の非特許文献１，２には、内在する伝搬路資源を活用する手だてとして空間多重化処理された信号を効率的に活用するＢＬＡＳＴ（Bell  Labs  Layered  Space-Time）と呼ばれる空間伝送処理が開示されている。また、このＢＬＡＳＴの空間多重分離を低複雑度で実現するアーキテクチャとして線形フィルタリングと干渉キャンセラとを組み合わせたＶ−ＢＬＡＳＴと呼ばれる手法も開示されている。【０００６】
線形フィルタリングとしては、干渉成分を抑圧（ヌリング：nulling）するＺＦ（Zero-Forcing）規範のもの、あるいは最小平均自乗誤差（ＭＭＳＥ：Minimum  Mean  Square  Error）規範のものが一般的である。【０００７】
ＺＦ規範にしたがって抑圧を行う線形変換としてはＭＰ（Moore-Penrose）の一般逆行列が知られており、干渉キャンセラの特性向上を目的に検出後のＳＮＲ（Signal  to  Noise  Ratio）が最も高いと簡易推定される順に検出する順序付け処理（オーダリング）がなされる。このシンボルの順序付けを行う操作としてＭＰ一般逆行列の重みベクトルに相当する最小ノルムを有する列ベクトルを優先して使うことが知られている。【０００８】
あるいは、さらに低複雑度化の手法としてＱＲ分解による方法がある。すなわち、通信路行列ＨをＱＲ分解によって、Ｈ＝Ｑ・Ｒとした後、ｎT次元の送信アンテナ信号ベクトル
【数１】
0004803384.tif 000002
と、ｎR次元の受信アンテナ信号ベクトル
【数２】
0004803384.tif 000003
との間に、
ＱH・Ｙ＝Ｒ・Ｘ＋ＱH・ｖ
という関係が成り立つ。【０００９】
ここで、
【数３】
0004803384.tif 000004
はユニタリー行列、
【数４】
0004803384.tif 000005
は上三角行列であり、雑音成分ベクトル
【数５】
0004803384.tif 000006
はユニタリー変換されるので、雑音強調は生ぜず、信号点間距離を維持したまま変換されることになる。【００１０】
このＱＲによる分解過程でＳＮＲが高い順序で処理できるように行列内ベクトルの並び替えが可能で、ＳＮＲが最大化されるような順序（オーダリング）で検出するステップ処理が実現できる。このような方法はＺＦ規範によるヌリング処理にあたるため、本質的に受信アンテナの本数ｎRが送信アンテナの本数ｎTと同数か、それ以上であることが前提である。【００１１】
しかしながら、これらの方法の欠点は、初回のステップでのヌリングによる線形処理でｎT−１次のヌル生成を行うため、ダイバシティ利得が、
ｎR−ｎT＋１
のオーダしか得られない。したがって、初回のステップにおける検出誤りが起こりやすく、その影響が後段の検出誤りを引き起こす誤り伝搬が生じる。【００１２】
一方、最適検出を行うには、
【数６】
0004803384.tif 000007
という式におけるＭＬＤ（Maximum  Likelihood  Decoding：最尤検出）になる。【００１３】
このため、アンテナの本数と変調信号点のサイズＡ、
【数７】
0004803384.tif 000008
とに対して指数関数的に複雑度が増大し、符号化を考慮に入れると、ＭＬＤは事実上不可能である。【００１４】
そこで、低複雑度化の手法としてターボ原理に基づく手法等が検討されている。上記の式は検出器のみについてのＭＬＤであるが、この複雑度を回避するため、及び上述のＶ−ＢＬＡＳＴにおける初段から後段への誤り伝搬による特性劣化、言い換えるとフェージング環境におけるダイバシティ利得を得る目的で、ＳＤ（Sphere  Decoding：球内復号）と呼ばれる復号法の適用が提案されている。【００１５】
ＳＤの基本的な考え方は受信信号点を中心とした適当な半径ｒの球に含まれる信号点について尤度計算を行い、限定された範囲内でＭＬＤを行うといったもので、半径ｒの選び方によって効率が決まってくる。あるいは、信号点の数を尤度の大きさによって限定することによって複雑度を回避する方法もある。【００１６】
【非特許文献１】“Layered  space-time  architecture  for  wireless  communications  in  a  fading  environment  when  using  multiple  antennas”（１９９６年、ベル研究所テクニカル・ジャーナル、ボリューム６、ナンバ２、４１〜５９頁）
【非特許文献２】“Capacity  of  multi-antenna  Gaussian  channels”（１９９９年１１月／１２月、ヨーロピアン・トランザクション・オン・テレコミュニケーション、５８５〜５９５頁）
【発明の開示】
【発明が解決しようとする課題】
【００１７】
新世代移動通信システムにおけるシステムスループット実現のためには、通信路容量の拡大策である時空間信号多重における信号分離を高性能で、しかも低複雑度で実現する必要がある。【００１８】
しかしながら、上述したＶ―ＢＬＡＳＴでは、方式自体に内在する誤り伝搬のため、低複雑度で実現できる代わりに、特性劣化を引き起こす。一方、最適検出であるＭＬＤでは、高性能化は実現できるものの、高複雑度で実現することになり、採用することができない。【００１９】
また、球内復号であるＳＤは半径ｒ内に入る信号点数によってその複雑度が変動し、装置化に適さない。あるいは、対象とする信号点の数を尤度によって限定するとしても、高性能を実現するためには複雑度の増加が必要で、後段に続く復号器との尤度のやりとりも考慮されていない。すなわち、最終的にはターボ符号やＬＤＰＣ（Low  Density  Parity  Check）といった符号化を含む最適受信を考えることが必要である。【００２０】
その復号器出力で性能が決まるのであるから、ターボ復号器やＬＤＰＣへ渡す尤度情報は対数尤度代数（Log  Likelihood  Algebra）にしたがった正確なものが必要で、これと軟入力軟出力復号器のExtrinsic情報による時空間反復復号を行うことによって、符号化を考慮した最尤復号を低複雑度で実現することになる。しかしながら、そのための正確な尤度情報が得られる構成とはなっていなかったという問題がある。【００２１】
また、所要システムスループットによってはアンテナ本数を増加させる必要が生じる。上述した説明よりアンテナ本数が増えるにしたがってその複雑度が増加するので、できるだけ増加分を押さえたいが、従来方法では難しいという問題がある。【００２２】
さらに、空間多重信号分離を行う軟入力軟出力空間多重信号検出器とターボ復号器やＬＤＰＣといった軟入力軟出力復号器との連接をターボ原理に基づく反復復号によって実現した場合、検出器と復号器との間でExtrinsic情報の受け渡しが行われる。【００２３】
しかしながら、従来行われているやり方は外符号用復号器に相当する軟入力軟出力復号器においてExtrinsic情報抽出のしにくさからIntrinsic情報を含んだままExtrinsic情報として使われたり、あるいは一部のExtrinsic情報で賄われたりしている。このため、ターボ原理に基づく反復回数による特性向上が数回で飽和してしまうという問題がある。【００２４】
さらにまた、既製の軟入力軟出力復号器を用いる場合、情報ビット系列に対する対数尤度比を出力としており、符号語（シンボル系列）に対する対数尤度比は出力されていない場合が多い。【００２５】
一方、軟入力軟出力検出器は、符号語に対してＭＬＤ検出となるべく作用するので、軟入力軟出力復号器から軟入力軟出力検出器へのExtrinsic情報は符号語に対するものが必要となる。その結果、軟入力軟出力復号器内部のトレリス上遷移確率伝搬を演算している箇所を変更することになり、既製コアブロックの変更を余儀なくされるという問題がある。【００２６】
そこで、本発明の目的は上記の問題点を解消し、ターボ原理に基づく時空間反復復号の特性を向上させることができる空間多重信号検出方法及びそれを用いる時空間反復復号器を提供することにある。
【課題を解決するための手段】
【００２７】
本発明による空間多重信号検出方法は、時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象を含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いている。【００２８】
本発明による他の空間多重信号検出方法は、時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象を含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際しては、自グループにおける条件付確率の条件となる事象を含む前記グループ内の最大条件付確率を示す送信系列を条件付確率の条件として計算するか、あるいは予め簡易推定された複数個の事象（送信系列）を条件として複数個計算し、前記グループ間の順序付けにしたがって各グループ内の条件付確率の計算を前段の最大条件付確率を示す送信系列のもとに計算していく処理を有することを特徴とし、
軟判定出力である該送信系列のビット尤度は、請求項６記載の半環（セミリング）を用いて
（１）対象となるビットが推定対象として存在する該グループにおいて対象ビットと前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理並びに該対象ビットに対する排反事象と該前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とを条件とする複数個の条件付確率をメトリックベースで検出する処理、
（２）次段のグループにおいて対象ビットを含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理並びに対象ビットに対する排反事象を含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する処理からなる処理、
（３）同様に対象となるビットが存在するグループ以降の段におけるグループにおいても該前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件として複数個の条件付確率をメトリックベースで検出する処理並びに対象ビットに対する排反事象を含む（１）の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とそれに伴って検出される前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）とを条件とする複数個の最大条件付確率をメトリックベースで検出する処理、
（４）以降最終段のグループまで（３）を繰り返す処理、
（５）（４）を完了後、対象ビットの対数尤度比を、前記（１）における対象ビットと前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とその対象ビットに対する排反事象と前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とをメトリックベースで引き算し、前記（２）から（４）の対象となるビットが存在するグループ以降の段のグループにおいても前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率と該対象ビットに対する排反事象を含む前段までの最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところの対数尤度比として検出する処理、
を用いて計算する処理を含んでいる。【００２９】
本発明による別の空間多重信号検出方法は、時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象を含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際しては、自グループにおける条件付確率の条件となる事象を含む前記グループ内の最大条件付確率を示す送信系列を条件付確率の条件として計算し、前記グループ間の順序付けにしたがって各グループ内の条件付確率の計算を前段の最大条件付確率を示す送信系列のもとに計算していく処理を有し、前記尤度を最大化する送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いることが特徴であって、最終段における処理が完了した後の復元抽出（リサンプリング）として、対象となるビットが推定対象として計算された前記条件付確率の組合せの集合からメトリックベースの最大尤度を選択する処理、並びに該対象ビットに対する排反事象が推定対象として計算された前記条件確率の組合せの集合からメトリックベースの最大尤度を選択する処理を有し、両メトリックの差分をとって対象ビットの軟判定出力とする処理を含んでいる。【００３０】
本発明による時空間反復復号器は、時空間多重信号分離において、軟入力軟出力検出器と軟入力軟出力復号器とを有し、該軟入力軟出力復号器は符号化前の情報ビット系列に対する対数尤度比（以下、ＬＬＲとする）を出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を元に前記軟入力軟出力検出器のアプリオリ（a  priori）入力を作っている。【００３１】
本発明による他の時空間反復復号器は、時空間多重信号分離において、軟入力軟出力検出器と軟入力軟出力復号器とを有し、該軟入力軟出力復号器は符号化前の情報ビット系列に対する対数尤度比（以下、ＬＬＲとする）を出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を元に前記軟入力軟出力検出器のソフトレプリカ入力を作っている。【００３２】
本発明による別の時空間反復復号器は、時空間多重信号分離において、軟入力軟出力検出器と軟入力軟出力復号器とを有し、前記軟入力軟出力検出器は請求項６記載のＭａｘ−ｌｏｇ領域における半環（セミリング）を用いた空間多重検出器であって、その軟出力に当たる対数尤度比に対して重み付けを行い、次段へのアプリオリ（a  priori）入力として用いている。【００３３】
すなわち、本発明の空間多重信号検出方法は、上記の問題に鑑みなされたものであり、空間多重信号分離における軟入力軟出力検出器で尤度を最大化する送信系列推定に半環（セミリング）を用いたメトリック演算手法を用いることによって、最適検出であるＭＬＤ（Maximum  Likelihood  Decoding：最尤検出）に近い高性能化を低複雑度で実現するものである。【００３４】
また、本発明の空間多重信号検出方法では、検出器に続く後段のターボ復号器やＬＤＰＣ（Low  Density  Parity  Check）といった軟入力軟出力復号器へ渡す尤度情報も対数尤度代数（Log  Likelihood  Algebra）にしたがった正確なものであり、最終的な性能である復号器出力の誤り率特性に対して劣化を生じることなく、本来の復号器性能を引き出すことが可能となる。【００３５】
さらに、本発明の空間多重信号検出方法では、尤度を最大化する送信系列を推定するために用いられる半環（セミリング）を用いたメトリック演算手法として、Ｍａｘ−ｌｏｇ領域における半環を用いることで、対象ビットの軟判定出力として複数のグループに分割した条件付確率の各グループ間における前段最大条件付確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）に対する対象ビットを含む最大条件付確率とその排反事象を含む最大条件付確率とのメトリックベースの差分とそれ以降の同最大条件付確率の差分との総和として表すことができるので、所要システムスループットに対応するため、アンテナ本数を増加させる場合でも複雑度を抑えた形で対応することが可能なようにしたものである。【００３６】
本発明の空間多重信号検出方法では、軟入力軟出力検出器と軟入力軟出力復号器とのターボ原理に基づく時空間反復復号構成をとったとしても、軟入力軟出力復号器におけるExtrinsic情報抽出を正確に行うことによって、ターボ原理に基づく反復回数による特性向上が数回で飽和してしまうといったことがないように構成するものである。【００３７】
また、本発明の空間多重信号検出方法では、既製の符号語（シンボル系列）に対する対数尤度比を軟出力として持たない軟入力軟出力復号器を用いる場合でも、前段の軟入力軟出力検出器へのアプリオリ（a  priori）を供給すべく、符号語に対応したExtrinsic情報を生成することができるように構成するものである。【００３８】
つまり、本発明の空間多重信号検出方法は、時空間多重信号分離における軟入力軟出力検出方法であり、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する手段を有し（ファクタライゼーション）、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際してはグループにおける条件付確率の条件となる事象を含むグループを先に処理する如く確率計算するグループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いることを特徴としている。【００３９】
また、本発明の空間多重信号検出方法は、上述した複数の条件付確率の積で表せるように分解する手段（ファクタライゼーション）及び分解可能な条件付確率を複数のグループに分割する手段において、各グループの条件付確率の条件となる事象を含むグループの条件付確率を簡易推定する簡易推定手段を有し（オーダリング）、該簡易推定手段によって、グループ間の順序付けをグループの条件付確率が高く簡易推定される順序で処理できるように分割することを特徴としている。【００４０】
本発明の空間多重信号検出方法では、上述した複数の条件付確率の積で表せるように分解する手段（ファクタライゼーション）として、ＱＲ分解を用いたことを特徴としている。【００４１】
または、本発明の空間多重信号検出方法では、上述した複数の条件付確率の積で表せるように分解する手段（ファクタライゼーション）として、ブロック三角化分解を用いたことを特徴としている。【００４２】
あるいは、本発明の空間多重信号検出方法では、上述した複数の条件付確率の積で表せるように分解する手段（ファクタライゼーション）として、三重対角化手法を用いたことを特徴としている。【００４３】
一方、本発明の空間多重信号検出方法では、上記二つの排反する事象の尤度の比によって送信系列を推定するために用いられる半環（セミリング）を用いたメトリック演算手法を、和を最大値演算（ＭＡＸ）、積を通常の加算として、
【数８】
0004803384.tif 000009
の如く行われることを特徴としている。【００４４】
また、本発明の空間多重信号検出方法では、上記二つの排反する事象の尤度の比によって送信系列を推定するために用いられる半環（セミリング）を用いたメトリック演算手法を、和をヤコビアン対数に基づく演算、積を通常の加算として、
【数９】
0004803384.tif 000010
ここで、
【数１０】
0004803384.tif 000011
の如く行われることを特徴としている。【００４５】
本発明の空間多重信号検出方法では、時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する手段を有し（ファクタライゼーション）、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際してはグループにおける条件付確率の条件となる事象を含むグループを先に処理する如く確率計算するグループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際しては、自グループにおける条件付確率の条件となる事象を含むグループ内の最大条件付確率を示す事象（送信系列）を条件付確率の条件として計算するか、あるいは予め簡易推定された複数個の事象（送信系列）を条件として複数個計算し、グループ間の順序付けにしたがって各グループ内の条件付確率の計算を行う手段を有することを特徴とし、軟判定出力である該送信系列のビット尤度は、和を最大値演算（ＭＡＸ）、積を通常の加算とする上記の半環（セミリング）を用いて、
（１）対象となるビットが推定対象として存在する該グループにおいて対象ビットと前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する手段、並びに該対象ビットに対する排反事象と該前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個のを大条件付確率をメトリックベースで検出する手段（対象となるビットが存在するグループ）
（２）次段のグループにおいて対象ビットを含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する手段、並びに対象ビットに対する排反事象を含む前段最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する手段（次段のグループ）
（３）同様に、対象となるビットが存在するグループ以降の段におけるグループにおいても該前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する手段、並びに対象ビットに対する排反事象を含む（１）の最大確率事象とそれに伴って検出される前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで検出する手段
（４）以降最終段のグループ迄、（３）を繰り返す手段（それ以降の段のグループ）
（５）（４）を完了後、対象ビットの対数尤度比を、（１）における対象ビットと前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率と、その対象ビットに対する排反事象と前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率とをメトリックベースで引き算し、（２）から（４）の対象となるビットが存在するグループ以降の段のグループにおいても前段迄の最大確率事象を条件とする最大条件付確率と該対象ビットに対する排反事象を含む前段迄の最大確率事象か、あるいは予め簡易推定された複数個の事象（送信系列）を条件とする複数個の条件付確率をメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところの対数尤度比として検出する手段
を用いて計算する手段を有することを特徴としている。【００４６】
本発明の空間多重信号検出方法では、時空間多重信号分離における軟入力軟出力検出方法において、空間多重された送信系列を送信したとするときに受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する手段を有し（ファクタライゼーション）、該分解可能な条件付確率を複数のグループに分割し、該尤度計算に際してグループにおける条件付確率の条件となる事象を含むグループを先に処理する如く確率計算するグループ間の順序付けが可能であって、該各グループ内の条件付確率計算に際して、自グループにおける条件付確率の条件となる事象を含むグループ内の最大条件付確率を示す事象（送信系列）を条件付確率の条件として計算し、グループ間の順序付けにしたがって各グループ内の条件付確率の計算を前段の最大条件付確率を示す事象（送信系列）を基に計算していく手段を有し、尤度を最大化する送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いることが特徴であって、最終段における処理が完了した後の復元抽出（リサンプリング）として、対象となるビットが推定対象として計算された条件付確率の組合せの集合からメトリックベースの最大尤度を選択する手段、並びに該対象ビットに対する排反事象が推定対象として計算された条件確率の組合せの集合からメトリックベースの最大尤度を選択する手段を有し、両メトリックの差分をとって対象ビットの軟判定出力とする手段を特徴として構成される。【００４７】
本発明の時空間反復復号器は、時空間多重信号分離において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、分解前の条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象である送信信号ベクトルを含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いる軟入力軟出力検出器と軟入力軟出力復号器を有し、該軟入力軟出力復号器は符号化前の情報ビット系列に対する対数尤度比（ＬＬＲ：Log Likelihood Ratio）を出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を元に軟入力軟出力検出器のアプリオリ（a priori）入力を作ることを特徴として構成される。
【００４８】
また、本発明の時空間反復復号器は、時空間検出部と復号部とからなる時空間多重信号分離装置において、該時空間検出部が軟入力軟出力検出器と軟入力軟出力符号器とを含み、該復号部が軟入力軟出力復号器を含む構成であって、軟入力軟出力復号器へのアプリオリ（a  priori）入力が軟入軟出符号器の出力に差し引く形で作用して軟入力軟出力検出器へのアプリオリ入力を形成し、軟入力軟出力検出部へのアプリオリ（a  priori）入力が軟入力軟出力検出器出力に差し引く形で作用して軟入力軟出力復号器へのアプリオリ（a  priori）入力を形成することを特徴として構成される。【００４９】
さらに、本発明の時空間反復復号器は、時空間多重信号分離において、空間多重された送信系列を送信したとする時に受信される信号系列に対して得られる尤度と呼ばれる条件付確率を、複数の条件付確率の積で表せるように分解する処理（ファクタライゼーション）を有し、分解前の条件付確率を複数のグループに分割し、該尤度計算に際しては前記グループにおける条件付確率の条件となる事象である送信信号ベクトルを含む前記グループを先に処理する如く確率計算する前記グループ間の順序付けが可能であって、該グループ内の確率計算に際し、二つの排反する事象の尤度の比によって送信系列を推定するために半環（セミリング）を用いたメトリック演算手法を用いる軟入力軟出力検出器と軟入力軟出力復号器とを有し、該軟入力軟出力復号器が符号化前の情報ビット系列に対するＬＬＲを出力とし、該対数尤度比を入力として符号化後の符号語系列に対する対数尤度比を出力とする軟入力軟出力符号器を有し、該軟入力軟出力符号器の出力を基に軟入力軟出力検出器のソフトレプリカ入力を作ることを特徴として構成される。
【００５０】
本発明の時空間反復復号器は、上記の軟入力軟出力符号器が送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素であるｑを法とする加算に代わって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1と対数尤度比ＬＬＲ2とに対してｑを法とする加算結果［＝ａ1＋ａ2（mod  q）］を事後値に持つＬＬＲを出力することを特徴として構成される。【００５１】
本発明の時空間反復復号器は、上記の軟入力軟出力符号器が送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素である２を法とする加算にかわって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1と対数尤度比ＬＬＲ2とに対して、
【数１１】
0004803384.tif 000012
を出力する手段、あるいはその近似値を出力する手段を特徴として構成される。【００５２】
本発明の時空間反復復号器は、上記の軟入力軟出力符号器が送信側の符号器と同じ構成を軟判定データに対応させたものであって、その構成要素である２を法とする加算に代わって二つの事後値ａ1と事後値ａ2とに対する対数尤度比ＬＬＲ1と対数尤度比ＬＬＲ2とに対してその絶対値の比較を行い、値の小さい方を選択する手段を有し、その選択結果に対して極性付けを対数尤度比ＬＬＲ1及び対数尤度比ＬＬＲ2のＭＳＢ（Most  Significant  Bit）に対する２を法とする加算結果に基づいて行うことを特徴として構成される。【００５３】
本発明の時空間反復復号器は、時空間多重信号分離において、軟入力軟出力検出器と軟入力軟出力復号器とを有し、軟入力軟出力検出器が上述したＭａｘ−ｌｏｇ領域における半環を用いた空間多重検出器であって、その軟出力に当たる対数尤度比に対して重み付けを行い、次段へのアプリオリ（a  priori）入力として用いることを特徴として構成される。
本発明の時空間反復復号器では、上記の重み付けを０．７５として用いることを特徴として構成される。
本発明の時空間反復復号器では、上記の重み付け０．７５をシフト加算によって実現したことを特徴として構成される。【００５４】
上記のように、本発明は、時空間反復復号方式並びに空間多重信号検出方法で、空間多重信号分離のための軟入力軟出力検出の尤度最大化による送信系列推定に半環（セミリング）を用いたメトリック演算を用いることによって、最適検出であるＭＬＤに近い性能でしかも低複雑度で実現可能な方法を提供するものであり、後段に続くターボ復号器やＬＤＰＣといった軟入力軟出力復号器に対しても特性劣化を引き起こすことのないように正確な尤度情報を渡すことが可能な方法を提供するものである。【００５５】
また、本発明は、半環（セミリング）を用いたメトリック演算で和を最大値演算、積を通常加算としたセミリングを適用すると、検出対象ビットの軟判定出力が分割した複数のグループにおける対象ビットを含む最大条件付確率とその排反事象を含む最大条件付確率のメトリックベースの差分、かつそれ以降の同最大条件付確率のメトリックベースの差分の総和として表すことが可能となるので、所要システムスループットに対応するためにアンテナ本数を増やし、処理段数を増加させた場合でも、絞り込まれた最低限の最大条件付確率事象のみによる追加グループ内処理で対処することが可能となるので、増加複雑度を抑えた形でアンテナ数増による処理に対処することが可能な方法を提供するものである。【００５６】
本発明の時空間反復復号器は、軟入力軟出力検出器と軟入力軟出力復号器との間で独立にアプリオリ（a  priori）情報のやりとりが可能となるので、軟入力軟出力復号器におけるExtrinsic情報の抽出を容易とし、Extrinsic情報抽出のしにくさから来るIntrinsic情報を含んだままのExtrinsic情報となったり、あるいは一部のExtrinsic情報で賄われたりすることがない。したがって、ターボ原理に基づく反復回数による特性向上が数回で飽和してしまうといった劣化を生じることのない時空間反復復号器を提供することが可能となる。【００５７】
また、本発明の時空間反復復号器は、既製の軟入力軟出力復号器出力である情報ビット系列に対する対数尤度比を符号語（シンボル系列）に対するＬＬＲに外部で変換可能となるので、既製コアブロックの内部変更をすることなく、実現可能な時空間反復復号器を提供することが可能となる。
【発明の効果】
【００５８】
本発明は、以下に述べるような構成及び動作とすることで、ターボ原理に基づく時空間反復復号の特性を向上させることができるという効果が得られる。
【図面の簡単な説明】
【００５９】
【図１】本発明の第１の実施例による時空間多重信号分離装置の構成を示すブロック図である。
【図２】本発明の第１の実施例による送信側装置の構成を示すブロック図である。
【図３】本発明の第５の実施例によるＱＲ分解を用いたグループ分割手段（factorization）の一例を示す図である。
【図４】本発明の第６の実施例による三重対角マトリクス分解を用いたグループ分割手段（factorization）の一例を示す図である。
【図５】本発明の第７の実施例によるブロック三角化分解を用いたグループ分割手段（factorization）の一例を示す図である。
【図６】ベイズ統計学における周辺化（マージナライゼーション：ｍａｒｇｉｎａｌｉｚａｔｉｏｎ）処理のイメージ図である。
【図７】グループ分割を条件付確率の積の分割と同じにした場合の一例を示す図である。
【図８】１６ＱＡＭにおける外円上信号点や内円上信号点等を説明するための図である。
【図９】貪欲法（ｇｒｅｅｄｙ）によるステップ処理で最大条件付確率事象の選択過程を示す図である。
【図１０】貪欲法（ｇｒｅｅｄｙ）によるステップ処理で複数の条件付確率の積でグループを構成した時の最大条件付確率事象選択過程を示す図である。
【図１１】本発明の第１８の実施例による時空間反復復号器の構成を示すブロック図である。
【図１２】本発明の第１９の実施例による時空間反復復号器の構成を示すブロック図である。
【図１３】本発明の第２０の実施例による符号器の構成例を示す図である。
【図１４】本発明の第２１の実施例による符号器の構成例を示す図である。
【図１５】図１４に使用されている軟入力軟出力要素符号器の内容を示す図である。
【図１６】軟入力軟出力要素符号器で使われる２を法とする加算をＬＬＲで実行する近似回路を示す図である。
【図１７】通常のターボ復号における符号化前の情報系列（Ｉ）に対するＬＬＲ算出過程を示すトレリス線図である。
【図１８】ターボ復号における符号化後のパリティ系列系列（Ｐ）に対するＬＬＲ算出過程を示すトレリス線図である。
【符号の説明】
【００６０】
１，５，７  軟入力軟出力検出器
２，６，８  軟入力軟出力復号器
３，５３，５５  デインタリーバ
４，３２，５９  インタリーバ
１１−１〜１１−ｎ，
３４−１〜３４−ｎ，
５１−１〜５１−ｎ，
７１−１〜７１−ｎ  アンテナ
１２  デテクタ
１３，１４  変換器
２１  デコーダ
３１  エンコーダ
３３  時空間マッパ
５２  時空間検出器
５４，５８  減算器
５６  加算器
５７，７３  軟入力軟出力符号器
６１，８１  軟入力軟出力復号器
６２，８２  判定器
７２  空間多重信号検出器
７４  ソフトレプリカ生成器
７２１  干渉キャンセラ
７２２  線形フィルタリング
【発明を実施するための最良の形態】
【００６１】
次に、本発明の実施例について図面を参照して説明する。図１は本発明の第１の実施例による時空間多重信号分離装置の構成を示すブロック図である。図１において、本発明の第１の実施例による時空間多重信号分離装置は軟入力軟出力検出器１と、軟入力軟出力復号器２と、デインタリーバ３と、インタリーバ４とから構成されている。【００６２】
軟入力軟出力検出器１はアンテナ１１−１〜１１−ｎと、デテクタ（Ｄｅｔｅｃｔｏｒ）１２と、変換器１３，１４とを備え、空間多重信号分離を行う。軟入力軟出力検出器１の処理方法としては、以下に述べるような空間多重信号検出方法が使われる。軟入力軟出力復号器２はデコーダ（Decoder）２１を備えている。【００６３】
図２は本発明の第１の実施例による送信側装置の構成を示すブロック図である。図２において、本発明の第１の実施例による送信側装置はエンコーダ（Encoder：符号器）３１と、インタリーバ３２と、時空間マッパ（Constellation  mapper）３３と、アンテナ３４−１〜３４−ｎとを備えている。【００６４】
この送信側装置において、送信対象となる情報系列はまずエンコーダ３１に入力し、ここで符号語（シンボル系列）に変換された後、インタリーバ３２によって攪拌され、時空間マッパ３３によって各信号点及び各アンテナ３４−１〜３４−ｎにマッピングされた後、送信系列として無線伝送路（図示せず）で空間多重される。【００６５】
図１に示す受信側では、上記のように空間多重された送信系列であるところの信号を軟入力軟出力検出器１で分離・抽出し、符号語（シンボル系列）の尤度情報としてインタリーバ３２の逆過程であるデインタリーバ３によって元の並びに入れ替えられた後、軟入力軟出力復号器２に入力される。【００６６】
軟入力軟出力復号器２は符号化後の符号語系列に対するExtrinsic情報を軟出力として発生する機能を有しており、送信系列の順序に合わせるため、インタリーバ４を介して軟入力軟出力検出器１へアプリオリ（a  priori）として入力される。【００６７】
この軟入力軟出力検出器１→デインタリーバ３→軟入力軟出力復号器２→インタリーバ４→軟入力軟出力検出器１で形成されたループ内を、上記の尤度情報を複数回伝達することによって、飛躍的に検出能力や復号能力を向上させる手法がターボ原理であり、そのためにターボ原理にしたがったExtrinsic情報を正確に抽出する必要がある。【００６８】
本実施例による空間多重信号検出方法は、上記のような構成において用いられるもので、以下、この空間多重信号検出方法について説明する。【００６９】
今、送信系列をｎT本の送信アンテナで送った場合の送信信号ベクトルをＸとし、それをｎR本の受信アンテナで受けた場合の受信信号ベクトルをＹとすると、ＭＬＤ（Maximum  Likelihood  Decoding：最尤検出）の場合、
【数１２】
0004803384.tif 000013
となる。【００７０】
したがって、条件付確率ｐ（Ｙ｜Ｘ）を計算することになる。ちなみに、上記の式の右辺のｐ（Ｘ）はアプリオリ（a  priori）に相当する部分である。また、送信信号ベクトルＸ及び受信信号ベクトルＹは、
【数１３】
0004803384.tif 000014
【数１４】
0004803384.tif 000015
である。また、上記の条件付確率ｐ（Ｙ｜Ｘ）を、
【数１５】
0004803384.tif 000016
として表し、送信信号系列ｘi  のとりうる信号点の数Ａを、
【数１６】
0004803384.tif 000017
とすると、条件付確率ｐ（Ｙ｜Ｘ）をメトリックベースで表した種類の合計は、
【数１７】
0004803384.tif 000018
個となる。【００７１】
次に、本発明の第２の実施例について説明する。本発明の第２の実施例において、送信側及び受信側の各装置の構成は、上述した本発明の第１の実施例と同様である。本発明の第２の実施例は、条件付確率の積で表せるように分解する手段（factorization）としてＱＲ分解を用いている点が上述した本発明の第１の実施例と異なる。このファクタライゼーションによる低複雑化を図った場合の空間多重信号検出方法について説明する。【００７２】
通信路行列を、
【数１８】
0004803384.tif 000019
とすると、ＱＲ分解によって、
Ｈ＝Ｑ・Ｒ
となるので、
【数１９】
0004803384.tif 000020
となる。【００７３】
ここで、
【数２０】
0004803384.tif 000021
はユニタリー行列、
【数２１】
0004803384.tif 000022
は上三角行列であり、雑音成分ベクトル
【数２２】
0004803384.tif 000023
はユニタリー変換された雑音ベクトルであるから、雑音強調は生ぜず、信号点間距離を維持したままの変換である。尚、このＱＲによる分解過程でＳＮＲ（Signal  to  Noise  Ratio）が高い順序で処理できるように行列内ベクトルの並び替えが可能で、ＳＮＲが最大化されるような順序（オーダリング）で検出するステップ処理を実現することもできる。【００７４】
上記の式を行列の要素レベルまで展開して再記すると、
【数２３】
0004803384.tif 000024
となる。ここで、
【数２４】
0004803384.tif 000025
はユニタリー変換後の雑音ベクトル要素で、統計的に独立であるから、ユニタリー変換後の受信信号ベクトルＺに対する条件付確率は、
【数２５】
0004803384.tif 000026
となって、送信系列の要素に対応した条件付確率の積で表せるように分解することができる。【００７５】
その送信系列の部分集合による送信信号ベクトルを、
【数２６】
0004803384.tif 000027
というように定義すると、上記の条件付確率は、
【数２７】
0004803384.tif 000028
となって送信信号系列ｘi  のとりうる信号点の数をＡとすると、条件付確率ｐ（Ｙ｜Ｘ）をメトリックベースで表した種類の合計が、
【数２８】
0004803384.tif 000029
となる。ここで、信号点の数Ａは、
【数２９】
0004803384.tif 000030
で表される。よって、本実施形態では、上述したＭＬＤの場合に比べて、約１／ｎRに複雑度が減少し、しかも誤り率特性がＭＬＤの場合と変わらないといった効果がある。【００７６】
続いて、本発明の第３の実施例について説明する。本発明の第３の実施例において、送信側及び受信側の各装置の構成は、上述した本発明の第１の実施例と同様である。本発明の第３の実施例は、条件付確率の積で表せるように分解する手段（factorization）として三重対角マトリクス分解（Tridiagonal  matrices  factorization）を用いている点が上述した本発明の第１の実施例と異なる。この三重対角マトリクス分解によるファクタライゼーションにて低複雑化を図った場合の空間多重信号検出方法について説明する。【００７７】
三重対角マトリクス分解の場合、行列の要素レベルまで展開して記すと、
【数３０】
0004803384.tif 000031
となる。したがって、Ｚに対する条件付確率は、
【数３１】
0004803384.tif 000032
となる。【００７８】
送信信号系列ｘi  のとりうる信号点の数をＡとすると、条件付確率ｐ（Ｚ｜Ｘ）をメトリックベースで表した種類の合計が、
【数３２】
0004803384.tif 000033
となる。ここで、信号点の数Ａは、
【数３３】
0004803384.tif 000034
で表される。【００７９】
よって、本実施形態では、上述したＭＬＤの場合に比べて、約Ａ0、
【数３４】
0004803384.tif 000035
に複雑度が減少し、しかも誤り率特性がＭＬＤの場合と変わらないといった効果がある。【００８０】
さらに、本発明の第４の実施例について説明する。本発明の第４の実施例において、送信側及び受信側の各装置の構成は、上述した本発明の第１の実施例と同様である。本発明の第４の実施例は、条件付確率の積で表せるように分解する手段（factorization）としてブロック三角化分解を用いている点が上述した本発明の第１の実施例と異なる。このブロック三角化分解によるファクタライゼーションにて低複雑化を図った場合の空間多重信号検出方法について説明する。【００８１】
ブロック三角化分解の場合、行列の要素レベルまで展開して記すと、
【数３５】
0004803384.tif 000036
となる。したがって、Ｚに対する条件付確率は、３次のブロックを使ったとすると、
【数３６】
0004803384.tif 000037
となって、上記の条件付確率は、
【数３７】
0004803384.tif 000038
となる。【００８２】
送信信号系列ｘi  のとりうる信号点の数をＡとすると、条件付確率ｐ（Ｚ｜Ｘ）をメトリックベースで表した種類の合計が、
【数３８】
0004803384.tif 000039
となる。ここで、信号点の数Ａは、
【数３９】
0004803384.tif 000040
で表される。よって、本実施形態では、上述したＭＬＤの場合に比べて、約３／ｎRに複雑度が減少し、しかも誤り率特性がＭＬＤの場合と変わらないといった効果がある。【００８３】
以上説明した方法は、最適検出であるＭＬＤの場合と同じ誤り率特性を示し、複雑度を軽減して実現した実施形態である。次に、貪欲法（greedy）の適用によってマルコフ連鎖とみなして分割し、さらに低複雑度化を図った実施例について説明する。【００８４】
本発明の第５の実施例では、まず、上述した受信信号ベクトルＺに対する条件付確率ｐ（Ｚ｜Ｘ）を送信系列の要素に対応した条件付確率の積で表せるように分解し、複数のグループに分けた後（ファクタライゼーション）、ステップ毎に条件に相当する候補をグループ内で絞り込む。この絞り込み具合によって、大小の差はあるものの、原理的にグループ間の誤り伝搬による特性劣化を伴うが、許容範囲内に収めて複雑度を軽減する。【００８５】
図３は本発明の第５の実施例によるＱＲ分解を用いたグループ分割手段（factorization）の一例を示す図である。図３においては、分割の手段（ファクタライゼーション）として上述のＱＲ分解を用いた場合で、三つのグループに分けて３ステップで処理を行った場合の処理を示している。【００８６】
図３において、実線で囲まれた部分が第一のステップで処理される部分、荒い破線で囲まれた部分が第二のステップで処理される部分、細かい破線で囲まれた部分が第三のステップで処理される部分である。【００８７】
図３に示す式において、上記と同様に、
【数４０】
0004803384.tif 000041
はユニタリー変換後の雑音ベクトル要素で統計的に独立であるから、ユニタリー変換後の受信信号ベクトルＺに対する条件付確率は、
【数４１】
0004803384.tif 000042
となる。【００８８】
これを図３に示すグループ毎に独立に各ステップの処理を行う。各ステップの処理で行う条件付確率の計算は以下の通りである。
＜ステップ１＞
【数４２】
0004803384.tif 000043
＜ステップ２＞
ステップ１で得られた全ての候補を条件として、
【数４３】
0004803384.tif 000044
という条件付確率の計算を行う。
＜ステップ３＞
ステップ１，２で得られた絞り込み候補を条件として、
【数４４】
0004803384.tif 000045
という条件付確率を計算する。【００８９】
図４は本発明の第６の実施例による三重対角マトリクス分解を用いたグループ分割手段（factorization）の一例を示す図である。図４においては、分解手段（factorization）として三重対角マトリクス分解を用いた場合で、三つのグループに分けて３ステップで処理を行った場合の処理を示している。【００９０】
図４において、実線で囲まれた部分が第一のステップで処理される部分、荒い破線で囲まれた部分が第二のステップで処理される部分、細かい破線で囲まれた部分が第三のステップで処理される部分である。受信信号ベクトルＺに対する条件付確率は、
【数４５】
0004803384.tif 000046
となる。【００９１】
これを図４に示すグループ毎に独立に各ステップの処理を行う。各ステップの処理で行う条件付確率の計算は以下の通りである。
＜ステップ１＞
【数４６】
0004803384.tif 000047
＜ステップ２＞
ステップ１で得られた全ての候補を条件として、
【数４７】
0004803384.tif 000048
という条件付確率の計算を行う。
＜ステップ３＞
ステップ１，２で得られた絞り込み候補を条件として、
【数４８】
0004803384.tif 000049
という条件付確率を計算する。【００９２】
図５は本発明の第７の実施例によるブロック三角化分解を用いたグループ分割手段（factorization）の一例を示す図である。図５においては、分解手段（factorization）としてブロック三角化分解を用いた場合で、三つのグループに分けて３ステップで処理を行った場合の処理を示している。【００９３】
図５において、実線で囲まれた部分が第一のステップで処理される部分、荒い破線で囲まれた部分が第二のステップで処理される部分、細かい破線で囲まれた部分が第三のステップで処理される部分である。受信信号ベクトルＺに対する条件付確率は、３次のブロックを使ったとすると、
【数４９】
0004803384.tif 000050
となって、上記の条件付確率は、
【数５０】
0004803384.tif 000051
となる。【００９４】
これを図５に示すグループ毎に独立に各ステップの処理を行う。各ステップの処理で行う条件付確率の計算は以下の通りである。
＜ステップ１＞
【数５１】
0004803384.tif 000052
＜ステップ２＞
ステップ１で得られた全ての候補を条件として、
【数５２】
0004803384.tif 000053
という条件付確率の計算を行う。
＜ステップ処理３＞
ステップ１，２で得られた絞り込み候補を条件として、
【数５３】
0004803384.tif 000054
という条件付確率を計算する。【００９５】
尚、上述した条件付確率の積で表せるように分解する手段として使ったマトリクス分解過程で各グループの条件付確率が高い順序でステップ処理ができるように分割することが可能で、例えば、各送信アンテナから送られてくる互いに直交したパイロット信号を用いて受信端で一番ＳＮＲの高い順に下から送信信号ベクトル要素の並べ替えを行い、それに対応した形で通信路行列Ｈの並び替えを行った後、マトリクス分解を行えば、上述したステップ１の条件付確率が他のグループより高くなることが期待できる。【００９６】
あるいは、通信路行列Ｈの各列ベクトルのノルムを計算して、ノルムの高い列ベクトル順に対応する送信信号ベクトル要素を下から並び替え、それに応じて通信路列Ｈの並び替えを行えば、上述したステップ１の条件付確率が他のグループより高くなることが期待できるし、ステップ順に条件付確率が並ぶように期待できる。【００９７】
以上のようにして得られた絞り込み候補からは、符号語（シンボル系列）のシンボルｘi  に対する確率ｐ（ｘi｜Ｚ）を計算する必要がある。これはベイズ統計学における周辺化（マージナライゼーション：marginalization）処理であり、
【数５４】
0004803384.tif 000055
というようにして求めることができる。これらの実際の演算はメトリック領域で行われる。したがって、乗算は全て加算に置き換わり、指数演算も必要ない。具体例な実施形態については後述する。【００９８】
メトリック領域で行われる演算としては、有力な低複雑化手法に半環（セミリング：Semi-ring）がある。上記の式に半環（セミリング）を適用すると、
【数５５】
0004803384.tif 000056
というようになる。ここで、ｆ（・）は確率ｐ（・）に対応したメトリックである。【００９９】
半環（セミリング）を用いた低複雑化手法を説明するに当たり、まず半環（セミリング：Semi-ring）の概念について説明する。尚、上記のような記載が、「１９９６年７月、アイ・イー・イー・イー・トランザクション・オン・インフォメーション・セオリ、第４２巻、第４号、１０７２〜１０９２頁（IEEE  Transactions  on  information  theory,Vol.42,No.4,pp.1072-1092）」（参考文献１）に開示されている。【０１００】
半環（semiring）は２種類の演算、
【数５６】
0004803384.tif 000057
が定義された集合semiringで構成され、任意の３つの元、
【数５７】
0004803384.tif 000058
の間で、
（閉包則）
【数５８】
0004803384.tif 000059
（結合則）
【数５９】
0004803384.tif 000060
（単位元）
【数６０】
0004803384.tif 000061
（ゼロ元）
【数６１】
0004803384.tif 000062
という関係を満たす。尚、半環（semiring）は環（ring）と異なり、逆元を有しない。【０１０１】
以上の性質を持つ半環（semiring）の代表的な二つの例（Max-Log領域及びＬｏｇ領域）、つまり（Ａ）Max-Log領域（Max-Log  Domain）での半環（セミリング）と、（Ｂ）Ｌｏｇ領域（Max-Log  Domain）での半環（セミリング）とを以下に示す。【０１０２】
（Ａ）において、Max-Log領域における半環（セミリング）は、
【数６２】
0004803384.tif 000063
という式で対応付けることができる。【０１０３】
すなわち、半環上の加算は最大値演算（ＭＡＸ）、半環上の積は通常の加算として定義される。尚、このような演算方法については、上記の参考文献１に記載されている。【０１０４】
これによって、上述したｍは−∞に相当し、「１の下線付き」は通常のゼロに相当することになる。この場合、上記の式の、
【数６３】
0004803384.tif 000064
はメトリックの最大値検出によって簡単に実行することができる。上記の方法は簡易なものであるが、その計算結果は近似値である。【０１０５】
（Ｂ）において、Ｌｏｇ領域における半環（セミリング）は、
【数６４】
0004803384.tif 000065
という式で対応付けることができる。【０１０６】
ここで、最初の演算はヤコビアンロガリズム（Jacobian  Logarithm）によって、
【数６５】
0004803384.tif 000066
というようにして行うことができる。この場合、補正項ｆ（｜＊｜）は簡単なルックアップテーブル（look-up  table）で実現することができる。【０１０７】
この補正が理想的に行われると、
【数６６】
0004803384.tif 000067
のメトリック演算によって、
【数６７】
0004803384.tif 000068
という周辺化処理を、メトリック領域上で正確に行うことができる。【０１０８】
本発明の第８の実施例は、以上の半環（セミリング）による低複雑化手法を、上述した貪欲法（greedy）を適用したステップ毎に処理する方法に適用したものである。以下の説明で用いる半環（セミリング）は、Max-Log領域の実施形態で、条件付確率の積で表せるように分解する手段（factorization）としてＱＲ分解を用いた場合である。【０１０９】
ユニタリー変換後の受信信号ベクトルＺに対する条件付確率は、上述したＱＲ分解を用いた場合の式から、
【数６８】
0004803384.tif 000069
となる。これにベイズ統計学における周辺化（マージナライゼーション：marginalization）処理を行うと、
【数６９】
0004803384.tif 000070
という式を得る。この式に半環（セミリング）を適用すると、
【数７０】
0004803384.tif 000071
というようになる。【０１１０】
これに対しては、図３に示すグループ毎に独立な各ステップの処理を行う。各ステップの処理で行う条件付確率の計算は以下の通りである。
＜ステップ１＞
【数７１】
0004803384.tif 000072
Max-Log領域の実施形態の場合、ステップ１のグループ内最大値になる。ここで、予め簡易推定された複数個の事象（送信系列）を条件として用いた場合、それぞれに対して複数個計算することになるが、周辺化処理に対する計算はそれぞれの最大値として簡略化される。この処理はステップ１のグループ内周辺化処理に相当する。図６にそのイメージを示す。
＜ステップ２＞
ステップ１で得られた最大値を持つ候補を条件として用いるか、あるいは予め簡易推定された複数個の事象（送信系列）を条件として用い、
【数７２】
0004803384.tif 000073
というメトリック演算を行う。予め簡易推定された複数個の事象（送信系列）を条件として用いると、複数個計算することになるが、Max-Log領域の実施形態の場合、ステップ２のグループ内最大値になって簡略化される。この処理はステップ２のグループ内周辺化処理に相当する。図６にそのイメージを示す。
＜ステップ３＞
ステップ１，２で得られた最大値を持つ候補を条件として用いるか、あるいは予め簡易推定された複数個の事象（送信系列）を条件として用い、
【数７３】
0004803384.tif 000074
というメトリック演算を行う。予め簡易推定された複数個の事象（送信系列）を条件として用いると、複数個計算することになるが、Max-Log領域の実施形態の場合、ステップ３のグループ内最大値になって簡略化される。この処理はステップ３のグループ内周辺化処理に相当する。図６にそのイメージを示す。【０１１１】
尚、上記の各ステップを通じて自グループ内に対象となる要素、すなわちｘi  ＝ａがない場合、そのステップでは条件のない最大値となる。したがって、半環（セミリング）を適用した周辺化処理の式は、
【数７４】
0004803384.tif 000075
というようになる。【０１１２】
ここで、各ステップの演算結果を以下のようにα，β，γと置く。
ステップ１：
【数７５】
0004803384.tif 000076
ステップ２：
【数７６】
0004803384.tif 000077
ステップ３：
【数７７】
0004803384.tif 000078
したがって、
【数７８】
0004803384.tif 000079
となる。【０１１３】
あるいは予め簡易推定された複数個の事象（送信系列）を条件として用いた場合、“’”をその複数個中の任意の送信系列に対応した演算結果や事象として用いると、各ステップの演算結果は複数個となり、α’，β’，γ’を任意として、
ステップ１：
【数７９】
0004803384.tif 000080
ステップ２：
【数８０】
0004803384.tif 000081
ステップ３：
【数８１】
0004803384.tif 000082
したがって、
【数８２】
0004803384.tif 000083
となる。この場合、複数個の事象に対して得られた結果の中で最大値をとる事象が最も確からしい送信系列となる。【０１１４】
Max-Log領域の実施形態の場合、前段の拘束条件の下で計算した各ステップのグループ内最大メトリックα，β，γの和が求めるｆ（ｘi＝ａ｜Ｚ）になる。尚、最大値のみの単一候補を用いた場合、ステップ内に対象となる要素、すなわちｘi  ＝ａがない場合、そのステップでは条件のない最大値となるので、他との共有化が可能で、さらに低複雑度で実現することができる。【０１１５】
以上の条件付確率の積で表せるように分解し（factorization）、分解可能な条件付確率を複数のグループに分割して各ステップの処理を行い、各ステップの処理における前段の最大条件付確率を示す事象、すなわちこの実施形態の場合、メトリックの最大値となる候補Ｘを次段の条件とすべく処理する過程を具体的なメトリック計算による実施形態で以下説明する。尚、本実施例では説明のし易さから最大値のみの単一候補を用いた場合を例にとって説明しているが、簡易推定によって複数個の事象（送信系列）に絞り込んだ候補を使った場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１１６】
以下の実施形態では、送信信号系列の要素ｘi  のとりうる信号点の数Ａを「１６」とした１６ＱＡＭ（Quadrature  Amplitude  Modulation）の場合において、説明のし易さからグループ分割を条件確率の積の分割と同じ場合について説明する。ここで、信号点の数Ａは、
【数８３】
0004803384.tif 000084
と表される。【０１１７】
この分割の様子を図７に示す。図７において、ステップ１に相当するメトリック処理は、
【数８４】
0004803384.tif 000085
という式によって行われる。【０１１８】
ここで、最終的な対象ビットの軟判定出力は対数尤度比（ＬＬＲ：Log  Likelihood  Ratio）であり、メトリックベースの演算では対象ビットに対するメトリックとその排反事象に対するメトリックとの差となるので、共通項であるｚ、
【数８５】
0004803384.tif 000086
は予め削除されている。【０１１９】
本発明の第９の実施例では、Max-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値をもって最大確率事象となる。【０１２０】
１６ＱＡＭでは一つの信号点に対して四ビット割り付けられる。そこで、
【数８６】
0004803384.tif 000087
と表すことにする。また、各信号点は、図８に示すように、外円上に信号点のあるグループ（outer  signal  point）と、内円上に信号点のあるグループ（inner  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（１／３）にあるグループ（±ｔａｎ-1（１／３）（ｍｏｄ  π）  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（３）にあるグループ（±ｔａｎ-1（３）  signal  point)とに分けることができる。【０１２１】
それぞれのグループにおけるメトリック計算は、outer  signal  point（３2＋３2＝１８）  ａ＝±１，ｂ＝±１，
【数８７】
0004803384.tif 000088
inner  signal  poinｔ（１2＋１2＝２）
【数８８】
0004803384.tif 000089
±ｔａｎ-1（１／３）（ｍｏｄ  π）  signal  point（３2＋１2＝１０）
【数８９】
0004803384.tif 000090
±ｔａｎ-1（３）  signal  point（１2＋３2＝１０）
【数９０】
0004803384.tif 000091
となる。【０１２２】
このメトリック計算によって得られた全ての組合せからMax-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値と最大確率事象、つまりその時の送信系列、
【数９１】
0004803384.tif 000092
を選び出す。尚、上記の例では各ビットに対して単一候補に絞り込んだ場合を例にとって説明しているが、簡易推定によって複数個に絞り込んだ場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１２３】
【数９２】
0004803384.tif 000093
以上の結果、得られた上記の結果がステップ１における
【数９３】
0004803384.tif 000094
となる。但し、この実施例の場合、グループ分割と条件確率の積の分割とが同じなので、
【数９４】
0004803384.tif 000095
となっている。この処理で特徴的なのは、全ての組合せを得るために加減算のみしか使っていない点であり、低複雑化がなされている。【０１２４】
ステップ２の処理は、前段の最大確率事象、すなわち各ビットに対するビット条件付最大事象、
【数９５】
0004803384.tif 000096
とその最大排反事象とに対して、それぞれ、
【数９６】
0004803384.tif 000097
というメトリック処理によって行われる。もちろん、簡易推定によって複数個の事象に絞りこんだ候補を使った場合でも、上記と同様の手段を講じることが可能である。【０１２５】
ここで、最終的な対象ビットの軟判定出力はＬＬＲであり、メトリックベースの演算では対象ビットに対するメトリックとその排反事象に対するメトリックとの差となるので、共通項であるｚ0、
【数９７】
0004803384.tif 000098
は予め削除されている。【０１２６】
本発明の第１０の実施例では、上述したのと同様に、Max-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値をもって最大確率事象となる。【０１２７】
同様に、１６ＱＡＭで一つの信号点に対して四ビット、
【数９８】
0004803384.tif 000099
と表すことにするが、この例の場合、最大値のみの単一候補に絞り込んだ場合を例にとっているので、前段の最大確率事象に基づく信号点であることを明確にするため、
【数９９】
0004803384.tif 000100
と表記することにする。【０１２８】
前段では、各ビットに対するビット条件付最大事象とその最大排反事象とが決定されているので、それぞれの事象に対して上記のメトリック処理を行っている。すなわち、
【数１００】
0004803384.tif 000101
の全てに対して計算を行う。尚、この条件事象の中には、無条件最大事象が必ず存在する。【０１２９】
各信号点は、図８に示すように、外円上に信号点のあるグループ（outer  signal  point）と、内円上に信号点のあるグループ（inner  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（１／３）にあるグループ（±ｔａｎ-1（１／３）（ｍｏｄ  π）  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（３）にあるグループ（±ｔａｎ-1（３）  signal  point）とに分ける。【０１３０】
それぞれのグループにおけるメトリック計算は
【数１０１】
0004803384.tif 000102
となる。【０１３１】
このメトリック計算によって得られた全ての組合せからMax-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値と最大確率事象、つまりその時の送信系列
【数１０２】
0004803384.tif 000103
を選び出す。尚、上記の例では各ビットに対して単一候補に絞り込んだ場合を例にとって説明しているが、簡易推定によって複数個に絞り込んだ場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１３２】
ここで、前段の無条件最大事象が、
【数１０３】
0004803384.tif 000104
であったとすると（複数個に絞り込まれた候補を使った場合、その候補それぞれが前段の条件事象となる）、
【数１０４】
0004803384.tif 000105
また、前段に対象ビットが存在している事象に対しては、その事象を条件とする最大メトリックを選び出し、
【数１０５】
0004803384.tif 000106
というようにする。【０１３３】
選択の過程を図９に示す。左からステップ１、ステップ２、・・・の処理となる。図９では図示の都合上、条件事象を省略して記載している。この例では、各ステップにおける無条件最大事象が、
【数１０６】
0004803384.tif 000107
の順となっている。尚、上記の例では説明のし易さから最大値のみの単一候補に絞り込んだ場合を例にとって説明しているが、簡易推定によって複数個に絞り込んだ場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１３４】
以上、得られた結果がステップ２における、
【数１０７】
0004803384.tif 000108
となる。但し、この実施例の場合、グループ分割と条件確率の積の分割とが同じなので、
【数１０８】
0004803384.tif 000109
となっている。この処理で特徴的なのは、全ての組合せを得るために加減算のみしか使っていない点であり、低複雑化がなされている。【０１３５】
次に、任意のステップｍの処理は、前段の最大確率事象、すなわち各ビットに対するビット条件付最大事象ｘm＝（αm，βm，γm，δm）と、その最大排反事象とに対して、それぞれ、
【数１０９】
0004803384.tif 000110
というメトリック処理によって行われる。もちろん、簡易推定によって複数個の事象に絞りこんだ候補を使った場合でも、上記と同様の手段を講じることが可能である。【０１３６】
ここで、最終的な対象ビットの軟判定出力はＬＬＲであり、メトリックベースの演算では対象ビットに対するメトリックとその排反事象に対するメトリックとの差となるので、共通項である、
【数１１０】
0004803384.tif 000111
は予め削除されている。【０１３７】
本発明の第１１の実施例では、上述したのと同様に、Max-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値をもって最大確率事象となる。【０１３８】
また、上記と同様に、１６ＱＡＭで一つの信号点に対して四ビットｘm＝（αm，βm，γm，δm）と表すことにするが、この例の場合、最大値のみの単一候補に絞り込んだ場合を例にとっているので、前段の最大確率事象に基づく信号点であることを明確にするため、
【数１１１】
0004803384.tif 000112
と表記することにする。【０１３９】
前段までの間には、各ビットに対するビット条件付最大事象とその最大排反事象とが決定されているので、それぞれの事象に対して上記のメトリック処理を行っている。すなわち、例えばステップ３の状態を例にすると、
【数１１２】
0004803384.tif 000113
の全てに対して計算を行う。尚、この条件事象の中には、無条件最大事象が必ず存在する。【０１４０】
各信号点は、図８に示すように、外円上に信号点の有るグループ（outer  signal  point）と、内円上に信号点のあるグループ（inner  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（１／３）にあるグループ（±ｔａｎ-1（１／３）（ｍｏｄ  π）  signal  point）と、信号点が水平軸に対して±ｔａｎ-1（３）にあるグループ（±ｔａｎ-1（３）  signal  point）とに分けることができる。【０１４１】
それぞれのグループにおけるメトリック計算は、
【数１１３】
0004803384.tif 000114
となる。【０１４２】
このメトリック計算によって得られた全ての組合せからMax-Log領域における半環（セミリング）の実施形態をとっているので、各候補の最大条件付確率、すなわちメトリックの最大値と最大確率事象、つまりその時の送信系列
【数１１４】
0004803384.tif 000115
を選び出す。尚、上記の例では各ビットに対して単一候補に絞り込んだ場合を例にとって説明しているが、簡易推定によって複数個に絞り込んだ場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１４３】
ここで、説明の都合上、ステップ３の状態として説明すると、前段までの無条件最大事象が、
【数１１５】
0004803384.tif 000116
であったとして（複数個に絞り込まれた候補を使った場合、その候補それぞれが前段の条件事象となる）、
【数１１６】
0004803384.tif 000117
また、前段に対象ビットが存在している事象に対しては、その事象を条件とする最大メトリックを選び出し、
【数１１７】
0004803384.tif 000118
というようにする。【０１４４】
選択の過程を図９に示す。左からステップ１、ステップ２、ステップ３、・・・の処理で、上記の例は三番目の場合である。図９では、図示の都合上、条件事象を省略して記載している。この例の各ステップにおける無条件最大事象は、
【数１１８】
0004803384.tif 000119
である。尚、上記の例では説明のし易さから最大値のみの単一候補を用いた場合を例にとって説明しているが、簡易推定によって複数個の事象（送信系列）に絞り込んだ候補を使った場合でも、上記と同様の手段を講じることが可能なことはいうまでもない。【０１４５】
以上、得られた結果が任意のステップｍ、この例ではステップ３における
【数１１９】
0004803384.tif 000120
となる。但し、この例の場合、グループ分割と条件確率の積の分割とが同じなので、
【数１２０】
0004803384.tif 000121
となっている。この処理で特徴的なのは、全ての組合せを得るために加減算のみしか使っていない点であり、低複雑化がなされている。【０１４６】
本発明の第１２の実施例では、このようにして得られた各ステップのグループ内ビット条件付最大メトリックを基に、空間多重信号検出における対象ビットのＬＬＲを求める。例えば、ステップ３までで完了する構成の場合において、
【数１２１】
0004803384.tif 000122
に対する軟判定出力であるところのＬＬＲを求める場合、
【数１２２】
0004803384.tif 000123
の関係から、
【数１２３】
0004803384.tif 000124
となる。【０１４７】
ここで、Ｌ（ｘ，ｘ，ｘ，ｘ｜１，ｘ，ｘ，ｘ）はステップ２における前段最大確率事象、
【数１２４】
0004803384.tif 000125
を条件とする最大メトリック、Ｌ（ｘ，ｘ，ｘ，ｘ｜１，ｘ，ｘ，ｘ，ｘ，ｘ，ｘ，ｘ）はステップ３における前前段最大確率事象、
【数１２５】
0004803384.tif 000126
と、それを条件とするステップ３での最大確率事象、
【数１２６】
0004803384.tif 000127
とを条件とする最大メトリックである。【０１４８】
Max-Log領域の実施形態の場合、前段の拘束条件の下で計算した各ステップのグループ内最大メトリックα，β，γの和が求めるｆ（ｘi  ＝ａ｜Ｚ）になるから、
【数１２７】
0004803384.tif 000128
として計算することができる。次に、この排反事象である
【数１２８】
0004803384.tif 000129
に対しては、
【数１２９】
0004803384.tif 000130
となる。したがって、
【数１３０】
0004803384.tif 000131
となる。【０１４９】
これによって、
【数１３１】
0004803384.tif 000132
に対する軟判定出力であるところのＬＬＲは、
【数１３２】
0004803384.tif 000133
となる。【０１５０】
すなわち、対象ビットが存在するグループにおけるステップ処理で、対象ビットと、前段までの最大確率事象を条件とする最大条件確率と、その排反事象と、前段までの最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、対象ビットが存在するグループ以降のグループにおいても前段までの対象ビットを含む最大確率事象を条件とする最大条件付確率と、その排反事象を含む最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとしているのである。尚、この例の場合、対象ビットがステップ１の処理に含まれているので、対象ビットが存在するグループにおいて前段までの最大確率事象を条件とする必要はない。【０１５１】
あるいは、予め簡易推定によって絞り込まれた複数個の事象（送信系列）を条件とした場合には、上記と同様に、ステップ３までで完了する構成の場合において、
【数１３３】
0004803384.tif 000134
に対する軟判定出力であるところのＬＬＲを求める場合、
【数１３４】
0004803384.tif 000135
の関係から、
【数１３５】
0004803384.tif 000136
となる。【０１５２】
ここで、“’”は予め簡易推定された複数個の事象（送信系列）を基に得られた結果を示しており、Ｌ（ｘ’，ｘ’，ｘ’，ｘ’｜１，ｘ’，ｘ’，ｘ’）はステップ２における前段確率事象、
【数１３６】
0004803384.tif 000137
を条件とするメトリック、Ｌ（ｘ’，ｘ’，ｘ’，ｘ’｜１，ｘ’，ｘ’，ｘ’，ｘ’，ｘ’，ｘ’，ｘ’）はステップ３における前段確率事象、
【数１３７】
0004803384.tif 000138
と、それを条件とするステップ３における前段確率事象、
【数１３８】
0004803384.tif 000139
とを条件とするメトリックである。【０１５３】
Max-Log領域の実施形態の場合、前段の拘束条件の下で計算した各ステップのグループ内最大メトリックα，β，γの和が求めるｆ（ｘi＝ａ｜Ｚ）になるから、
【数１３９】
0004803384.tif 000140
として計算することができる。【０１５４】
次に、この排反事象である
【数１４０】
0004803384.tif 000141
に対しては、
【数１４１】
0004803384.tif 000142
となる。したがって、
【数１４２】
0004803384.tif 000143
となる。【０１５５】
これによって、
【数１４３】
0004803384.tif 000144
に対する軟判定出力であるところのＬＬＲは、
【数１４４】
0004803384.tif 000145
となる。【０１５６】
すなわち、対象ビットが存在するグループにおけるステップ処理で、対象ビットと、前段までの予め簡易推定された複数個の事象（送信系列）とを基に得られた確率事象を条件とする条件確率と、その排反事象と、前段までの予め簡易推定された複数個の事象（送信系列）とを基に得られた確率事象を条件とする条件付確率とをメトリックベースで引き算し、対象ビットが存在するグループ以降のグループにおいても、前段までの予め簡易推定された複数個の事象（送信系列）を基に得られた確率事象を条件とする条件付確率と、その排反事象を含む確率事象を条件とする条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとしているのである。【０１５７】
本発明の第１３の実施例は、対象ビットがステップ２のグループに含まれている場合を示している。例えば、ステップ３までで完了する構成の場合において、
【数１４５】
0004803384.tif 000146
に対する軟判定出力であるところのＬＬＲを求める場合、
【数１４６】
0004803384.tif 000147
の関係から、
【数１４７】
0004803384.tif 000148
となる。【０１５８】
ここで、
【数１４８】
0004803384.tif 000149
におけるαは排反事象と共通なので、予め削除されている。また、
【数１４９】
0004803384.tif 000150
はステップ２におけるビット条件付最大事象、
【数１５０】
0004803384.tif 000151
における最大メトリックで、前段の条件
【数１５１】
0004803384.tif 000152
はステップ１における無条件最大確率事象である。また、Ｌ（ｘ，ｘ，ｘ，ｘ｜ｘ，１，ｘ，ｘ，  ｘ，ｘ，１，ｘ）はステップ３における前前段無条件最大確率事象、
【数１５２】
0004803384.tif 000153
と、それを条件とするステップ２でのビット条件付最大確率事象、
【数１５３】
0004803384.tif 000154
とを条件とする無条件最大メトリックである。【０１５９】
Max-Log領域の実施形態の場合、前段の拘束条件の下で計算した各ステップのグループ内最大メトリックα，β，γの和（この場合、αは共通項なので削除され、βとγとの和となる）が求めるｆ（ｘi＝ａ｜Ｚ）になるので、
【数１５４】
0004803384.tif 000155
として計算することができる。【０１６０】
次に、この排反事象である
【数１５５】
0004803384.tif 000156
に対しては、
【数１５６】
0004803384.tif 000157
となる。したがって、
【数１５７】
0004803384.tif 000158
となる。【０１６１】
これによって、
【数１５８】
0004803384.tif 000159
に対する軟判定出力であるところのＬＬＲは、
【数１５９】
0004803384.tif 000160
となる。【０１６２】
すなわち、対象ビットが存在するグループにおけるステップ処理で、対象ビットと前段までの最大確率事象を条件とする最大条件確率と、その排反事象と、前段までの最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、対象ビットが存在するグループ以降のグループにおいても、前段までの対象ビットを含む最大確率事象を条件とする最大条件付確率と、その排反事象を含む最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとしているのである。【０１６３】
あるいは、予め簡易推定によって絞り込まれた複数個の事象（送信系列）を条件とした場合で、“’”を予め簡易推定された複数個の事象を基に得られた結果として、上記と同様に、ステップ３までで完了する構成の場合において、
【数１６０】
0004803384.tif 000161
に対する軟判定出力であるところのＬＬＲを求める場合、
【数１６１】
0004803384.tif 000162
の関係から、
【数１６２】
0004803384.tif 000163
となる。【０１６４】
ここで、
【数１６３】
0004803384.tif 000164
はステップ２におけるビット条件付き事象、
【数１６４】
0004803384.tif 000165
におけるメトリックで、前段の条件
【数１６５】
0004803384.tif 000166
はステップ１における予め簡易推定された複数個の事象（送信系列）を基に得られた確率事象である。また、Ｌ（ｘ’，ｘ’，ｘ’，ｘ’｜ｘ’，ｘ’，ｘ’，ｘ’，ｘ’，ｘ’，１，ｘ’）はステップ３における前段確率事象、
【数１６６】
0004803384.tif 000167
と、それを条件とするステップ２でのビット条件付確率事象、
【数１６７】
0004803384.tif 000168
とを条件とするメトリックである。【０１６５】
Max-Log領域の実施形態の場合、前段の拘束条件の下で計算した各ステップのグループ内最大メトリックα，β，γの和が求めるｆ（ｘi  ＝ａ｜Ｚ）になるので、
【数１６８】
0004803384.tif 000169
として計算することができる。【０１６６】
次に、この排反事象である
【数１６９】
0004803384.tif 000170
に対しては、
【数１７０】
0004803384.tif 000171
となる。したがって、
【数１７１】
0004803384.tif 000172
となる。【０１６７】
これによって、
【数１７２】
0004803384.tif 000173
に対する軟判定出力であるところのＬＬＲは、
【数１７３】
0004803384.tif 000174
となる。【０１６８】
ここで、上式第一項の
【数１７４】
0004803384.tif 000175
は、
【数１７５】
0004803384.tif 000176
を求めるに当たっての
【数１７６】
0004803384.tif 000177
と、
【数１７７】
0004803384.tif 000178
を求めるに当たっての
【数１７８】
0004803384.tif 000179
でその差分である。【０１６９】
すなわち、対象ビットが存在するグループにおけるステップ処理で、対象ビットと、前段までの予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に得られた条件確率と、その排反事象と、同じくその前段までの予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に得られた条件付確率とをメトリックベースで引き算し、対象ビットが存在するグループ以降のグループにおいても、前段までの対象ビットを含む確率事象を条件とする条件付確率と、その排反事象を含む確率事象を条件とする条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとしているのである。【０１７０】
尚、以上の実施形態は、説明のし易さから、グループ分割を条件確率の積の分割と同じにしているが、複数の条件確率の積でグループ化しても、上記と同様に実施することがで
きる。【０１７１】
本発明の第１４の実施例は、図１０に示す複数の条件確率の積を初段に用いた場合の例である。この例の場合も説明のし易さから最大値に候補を絞った場合を例にとって説明するが、簡易推定によって複数個に絞り込んだ場合も、上記と同様の手段が成り立つことはいうまでもない。初段の事象を上述した表記法に合わせて記すと、以下のようになる。【０１７２】
１６ＱＡＭで一つの信号点に対して四ビット、図１０に示す例ではステップ１の処理として二つの信号点を用いているので、
【数１７９】
0004803384.tif 000180
と表すことにする。したがって、次段の中に対象ビットが含まれている場合の選択対象は、
【数１８０】
0004803384.tif 000181
となり、この全てに対して計算を行う。この条件事象の中には、無条件最大事象が必ず存在する。この処理は上述したグループ分割を条件確率の積の分割と同じとした場合のステップ３の処理と同じものである。【０１７３】
図１０の選択の過程は、左からステップ１、ステップ２、ステップ３の処理で上記の例は三列目のステップ２の処理である。図１０では、図示の都合上、条件事象を省略して記載している。この例の各ステップにおける無条件最大事象は、
【数１８１】
0004803384.tif 000182
である。尚、通常、縮退があるので、サンプル数は図中以下の数となる。【０１７４】
本発明の第１５の実施例の場合、初段に複数アンテナ分の処理を同時に行うので、その分ダイバーシティ利得が得られ、次段以降への誤り伝搬の影響を改善することができる。【０１７５】
以上の処理を纏めると、以下のようになる。
（１）対象となるビットが存在するグループ：対象となるビットが推定対象として存在する該グループにおいて、対象ビットと前段までの最大確率事象を条件とする対象ビット条件付最大条件付確率をメトリックベースで検出するとともに、その対象ビットに対する排反事象と該前段までの最大確率事象を条件とする排反ビット条件付最大条件付確率とをメトリックベースで検出する。
（２）次段のグループ：次段のグループにおいて、対象ビットを含む前段最大確率事象を条件とする無条件ビット最大条件付確率をメトリックベースで検出するとともに、対象ビットに対する排反事象を含む前段最大確率事象を条件とする無条件ビット最大条件付確率をメトリックベースで検出する。
（３）同様に、対象となるビットが存在するグループ以降の段におけるグループにおいても、該前段までの最大確率事象を条件とし、無条件ビット最大条件付確率をメトリックベースで検出するとともに、対象ビットに対する排反事象を含む（１）の最大確率事象と、それに伴って検出される前段までの最大確率事象を条件とする無条件ビット最大条件付確率とをメトリックベースで検出する。
（４）それ以降の段のグループ：以降、最終段のグループまで（３）を繰り返す手段。（５）（４）を完了後、対象ビットのＬＬＲを、（１）における対象ビットと、前段までの最大確率事象を条件とする最大条件付確率と、その対象ビットに対する排反事象と、前段までの最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、（２）から（４）の対象となるビットが存在するグループ以降の段のグループにおいても、前段までの最大確率事象を条件とする最大条件付確率と、該対象ビットに対する排反事象を含む前段までの最大確率事象を条件とする最大条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとして検出する。
また、対象ビットが存在するグループの前後のグループで、上記の例ではすべてメトリックの計算対象としているが、適当な範囲で打ち切ることによって、低複雑化が可能であり、そのような使用形態も可能である。【０１７６】
尚、本発明の空間多重検出方法は、反復処理を用いない場合でも、低複雑度化に有効となり、そのような使用形態も適宜なされるのはいうまでもない。あるいは、予め簡易推定された複数個の事象（送信系列）を条件とした場合について上記の処理を纏めると、（１）対象となるビットが存在するグループ：予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に、対象となるビットが推定対象として存在する該グループにおいて、対象ビットと前段までの確率事象を条件とする対象ビット条件付確率をメトリックベースで検出するとともに、その対象ビットに対する排反事象と該前段までの確率事象を条件とする排反ビット条件付確率とをメトリックベースで検出する。
（２）次段のグループ：予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に、次段のグループにおいて、対象ビットを含む前段確率事象を条件とする無条件ビット条件付確率をメトリックベースで検出するとともに、対象ビットに対する排反事象を含む前段確率事象を条件とする条件付確率をメトリックベースで検出する。
（３）同様に、予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に、対象となるビットが存在するグループ以降の段におけるグループにおいても、該前段までの確率事象を条件とし、条件付確率をメトリックベースで検出するとともに、対象ビットに対する排反事象を含む（１）の確率事象と、それに伴って検出される前段までの確率事象を条件とする条件付確率とをメトリックベースで検出する。
（４）それ以降の段のグループ：予め簡易推定によって絞り込まれた複数個の事象（送信系列）を基に、以降、最終段のグループまで（３）を繰り返す手段。（５）（４）を完了後、対象ビットのＬＬＲを、（１）における対象ビットと、前段までの確率事象を条件とする条件付確率と、その対象ビットに対する排反事象と、前段までの確率事象を条件とする条件付確率とをメトリックベースで引き算し、（２）から（４）の対象となるビットが存在するグループ以降の段のグループにおいても、前段までの確率事象を条件とする条件付確率と、該対象ビットに対する排反事象を含む前段までの確率事象を条件とする条件付確率とをメトリックベースで引き算し、それぞれの総和をもって対象ビットの軟判定出力であるところのＬＬＲとして検出する。また、対象ビットが存在するグループの前後のグループで、上記の例ではすべてメトリックの計算対象としているが、適当な範囲で打ち切ることによって、低複雑化が可能であり、そのような使用形態も可能である。【０１７７】
以上、メトリック領域で行われる半環（セミリング：Semi-ring）を用いた低複雑化方法について説明したが、最大値単一候補に絞り込んだ場合で、さらに特性向上のために復元抽出（リサンプリング）がある。【０１７８】
この復元抽出（リサンプリング）は貪欲法（greedy）によるステップ処理のために生じる誤り伝搬の軽減を狙ったもので、最終的に推定対象となるビットとその排反事象を含む上述したメトリックベースの条件付確率の集合とが得られるが、その中からビット条件付の最大メトリックを再選択することによって実現するものである。【０１７９】
本発明の第１６の実施例は、排反事象に対しても、上記と同様に、選択されたものを採用し、両メトリックの差分を用いて、上述した方法と同様な方法によって対象ビットの軟判定出力であるところのＬＬＲを検出するものである。例えば、図９に示す場合、最終的に縮退を無視すると、
【数１８２】
0004803384.tif 000183
個のサンプル数を得ることができる。同例の場合、
【数１８３】
0004803384.tif 000184
を代入すると、
【数１８４】
0004803384.tif 000185
個のサンプルとなり、この集合から上述の条件にある最大メトリックを復元抽出（リサンプリング）する。【０１８０】
本発明の第１７の実施例は、上記と同様に、初段に複数アンテナ分の処理を同時に行い、ダイバーシティ利得を得た段以降への誤り伝搬の影響を改善している。図１０において、最終的に縮退を無視すると、
【数１８５】
0004803384.tif 000186
個のサンプル数を得ることができる。上記と同様に、
【数１８６】
0004803384.tif 000187
を代入すると、
【数１８７】
0004803384.tif 000188
個のサンプルとなり、この集合から上述した条件にある最大メトリックを復元抽出（リサンプリング）する。【０１８１】
以上、本発明の空間多重信号検出方法を用いた軟入力軟出力検出方法について説明したが、次に、時空間多重信号分離における軟入力軟出力検出器と軟入力軟出力復号器とを用いた時空間反復復号器について説明する。【０１８２】
図１１は本発明の第１８の実施例による時空間反復復号器の構成を示すブロック図である。図１１において、本発明の第１８の実施例による時空間反復復号器は軟入力軟出力検出部５と、軟入力軟出力復号部６とから構成されている。【０１８３】
軟入力軟出力検出部５はアンテナ５１−１〜５１−ｎと、時空間検出器５２と、デインタリーバ５３，５５と、減算器５４，５８と、加算器５６と、軟入力軟出力符号器５７と、インタリーバ５９とから構成され、軟入力軟出力復号部６は軟入力軟出力復号器６１と、判定器６２とを備えている。【０１８４】
軟入力軟出力復号部６は符号化前の情報ビット系列に対するＬＬＲを出力としている。軟入力軟出力検出部５内にある軟入力軟出力符号器５７はそのＬＬＲを入力とし、符号化後の符号語系列に対するＬＬＲを出力する。軟入力軟出力検出部５はこの符号化後の符号語に対するＬＬＲを基に、空間多重信号検出用のアプリオリ（a  priori）情報を生成し、このアプリオリ情報を基にターボ原理に基づいて軟入力軟出力検出部５と軟入力軟出力復号部６との間で反復復号を実行している。【０１８５】
軟入力軟出力復号部６は軟入力軟出力復号器６１を含む構成であって、軟入力軟出力復号器６１へのアプリオリ（a  priori）入力は軟入力軟出力符号器５７の軟判定出力に減算器５８によって差し引く形で作用し、軟入力軟出力検出部５内にある時空間検出器５２へのアプリオリ（a  priori）入力を形成している。この時、Intrinsic情報も加算器５６を介して減算器５８によって差し引くように作用している。【０１８６】
また、時空間検出器５２へのアプリオリ（a  priori）入力は時空間検出器５２の軟判定出力に減算器５４によって差し引く形で作用し、軟入力軟出力復号器６１へのアプリオリ（a  priori）入力を形成している。【０１８７】
尚、図２に示すように、送信側である符号器エンコーダ３１と時空間マッパ３３との間にはインタリーバ３２があり、これがターボ原理に基づく反復復号に重要な役割を果たしている。すなわち、図１に示すように、受信側でこのインタリーバ３２の逆過程であるデインタリーバ３及びインタリーバ４を入れることになり、軟入力軟出力検出処理と軟入力軟出力復号処理とを統計的に独立させ、Extrinsic情報、転じてアプリオリ（a  priori）とすることによって検出能力や復号能力を飛躍的に向上させる、所謂、ターボ原理を用いた手法であり、図１１に示すデインタリーバ５３，５５及びインタリーバ５９がその役割を担っている。【０１８８】
図１１に示す例の場合、軟入力軟出力検出部５から軟入力軟出力復号部６へのIntrinsic情報とアプリオリ（a  priori）情報とを別々に受け渡しているので、デインタリーバ５３，５５が二つとなっている。反復処理によって最終的に得られたＬＬＲ結果は、軟入力軟出力復号部６の中にある判定器６２によって硬判定され、データとして出力される。【０１８９】
図１２は本発明の第１９の実施例による時空間反復復号器の構成を示すブロック図である。本発明の第１９の実施例による時空間反復復号器はソフトキャンセラによる時空間反復復号器であり、本来のターボ原理に基づく反復復号ではなく、ソフトレプリカによる干渉キャンセラと線形フィルタリングとによる干渉抑圧（ヌリング：nulling）を基に構成されている。【０１９０】
図１２において、本発明の第１９の実施例による時空間反復復号器は軟入力軟出力検出部７と、軟入力軟出力復号部８とから構成されている。軟入力軟出力検出部７はアンテナ７１−１〜７１−ｎと、空間多重信号検出器７２と、軟入力軟出力符号器７３と、ソフトレプリカ生成器７４とから構成され、空間多重信号検出器７２内には干渉キャンセラ７２１と線形フィルタリング７２２とを備えている。また、軟入力軟出力復号部８は軟入力軟出力復号器８１と、判定器８２とを備えている。【０１９１】
本実施例は、本来のExtrinsic情報をやり取りするターボ原理に基づいたものではないので、反復処理による特性向上が数回で飽和してしまう欠点がある。しかしながら、本実施例は、簡易に実現することができるという利点がある。【０１９２】
図１２における時空間多重信号分離は空間多重信号検出器（軟入力軟出力検出器）７２を含む軟入力軟出力検出部７と、軟入力軟出力復号器８１を含む軟入力軟出力復号部８とから構成されている。【０１９３】
軟入力軟出力復号器８１は符号化前の情報ビット系列に対するＬＬＲを出力としている。軟入力軟出力検出部７の中には軟入力軟出力符号器７３があって、そのＬＬＲを入力とし、符号化後の符号語系列に対するＬＬＲを出力する。【０１９４】
ソフトレプリカ生成器７４はその符号語系列に対するＬＬＲを基にソフトレプリカを生成し、空間多重信号検出器７２へ送る。空間多重信号検出器７２は送られてきたソフトレプリカを基に干渉キャンセラ７２１にてキャンセリングを行い、線形フィルタリング７２２によるヌリングを行って空間多重信号を分離する。【０１９５】
ソフトレプリカ生成器７４は軟入力軟出力符号器７３からの符号語系列に対するＬＬＲをソフトレプリカに変換するもので、ＢＰＳＫ（Binary  Phase  Shift  Keying）の場合、
E[xi]=(+1)・p(xi  =+1|y)+(-1)p(xi  =-1|y)
という関係から実現することができる。【０１９６】
ここで、符号語系列に対するＬＬＲをＬ（ｘi）と置くと、
【数１８８】
0004803384.tif 000189
となる。同様に、
【数１８９】
0004803384.tif 000190
となる。したがって、
【数１９０】
0004803384.tif 000191
となる。【０１９７】
さらに、
【数１９１】
0004803384.tif 000192
として符号語系列に対するＬＬＲ［Ｌ（ｘi）］からソフトレプリカＥ［ｘi  ］を求めることができる。【０１９８】
以上の処理はＢＰＳＫの場合についてであるが、ＱＰＳＫ（Quadrature  Phase  Shift  Keying）の場合には２ビット構成で、Ｅ［ａ2i］，Ｅ［ａ2i+1］とすると、シンボルレプリカは、
E[xi  ]=E[a2i]+jE[a2i+1]
となる。【０１９９】
あるいは、１６ＱＡＭの場合には４ビット構成で、Ｅ［ａ4i］，Ｅ［ａ4i+1］，Ｅ［ａ4i+2］，Ｅ［ａ4i+3］とすると、例えば、
E[xi  ]=E[a4i]・(2-E[a4i+2])+jE[a4i+1]・(2-E[a4i+3])
といったマッピングで、シンボルレプリカＥ［ｘi］を構成することができる。【０２００】
以上の反復キャンセリング処理によって、最終的に得られたＬＬＲの結果は、軟入力軟出力復号部８の中にある判定器８２によって硬判定され、データとして出力される。【０２０１】
次に、上述した処理で使われてきた軟入力軟出力符号器の構成例について説明する。【０２０２】
軟入力軟出力符号器は送信側の符号器、例えば図２に示すエンコーダ３１と同じ構成を軟判定データに対応させたもので、その構成要素であるｑを法とする加算、すなわち、ａ1＋ａ2（ｍｏｄ  ｑ）といった計算をＬＬＲで実行したものである。【０２０３】
図１３は本発明の第２０の実施例による符号器の構成例を示す図である。図１３においては、送信側で使用される符号器の一例を示しており、ターボ符号で用いられる再帰的組織的畳み込みによる並列連接という手法を示している。【０２０４】
図１３において、使用されている信号は１ビットであり、使用されている加算器も２を法とした加算で、ａ1  ＋ａ2  （ｍｏｄ  ２）といった計算になる。すなわち、本実施例による符号器は、排他的論理和で簡単に実現することができる。【０２０５】
図１４は本発明の第２１の実施例による符号器の構成例を示す図である。図１４においては、上記の図１３に示す符号器と同じ構成を、軟判定データに対応させたものである。【０２０６】
図１５は図１４に使用されている軟入力軟出力要素符号器の内容を示す図である。この軟入力軟出力要素符号器の構成は図１３に示す送信側で使用されている符号器と同じものであるが、使用されている信号は複数ビットの軟判定データであり、２を法とする加算器もＬＬＲで実行したものである。【０２０７】
以下、２を法とする加算をＬＬＲで実行する場合の計算方法について説明する。この場合、ＬＬＲは、
【数１９２】
0004803384.tif 000193
という式で定義する。したがって、
【数１９３】
0004803384.tif 000194
となる。よって、
【数１９４】
0004803384.tif 000195
となる。【０２０８】
この関係から、
【数１９５】
0004803384.tif 000196
となるので、分母分子を、
【数１９６】
0004803384.tif 000197
で割ると、
【数１９７】
0004803384.tif 000198
となる。【０２０９】
ここで、
【数１９８】
0004803384.tif 000199
であるから、
【数１９９】
0004803384.tif 000200
となり、さらに、
【数２００】
0004803384.tif 000201
であるから、整理すると、
【数２０１】
0004803384.tif 000202
となる。【０２１０】
すなわち、
【数２０２】
0004803384.tif 000203
であるから、書き直して、
【数２０３】
0004803384.tif 000204
となって、２を法とする加算をＬＬＲで実行することができる。【０２１１】
この演算をそのままテーブルを用いて実装しても良いが、
【数２０４】
0004803384.tif 000205
という近似を使って計算することができる。あるいは、書き直して、
【数２０５】
0004803384.tif 000206
となる。【０２１２】
つまり、二つの対数尤度比ＬＬＲ1，ＬＬＲ2に対して、その絶対値の比較を行い、値の小さい方を選択し、その選択結果に対して極性付けを二つの対数尤度比ＬＬＲ1，ＬＬＲ2のＭＳＢ（Most  Significant  Bit）に対する２を法とする加算結果に基づいて行えばよいことになる。実際に回路を構成した例を図１６に示す。【０２１３】
以上、具体的な例として図１３に示す送信側で使用される符号器の一例としてターボ符号で用いられる再帰的組織的畳み込みによる並列連接という手法で構成した場合の受信側での本発明の第２１の実施例である軟入力軟出力要素符号器について説明したが、このターボ符号に限らず、ＬＤＰＣ（Low  Density  Parity  Check）等でも同様に本発明を用いることができることはいうまでもない。【０２１４】
本発明の軟入力軟出力符号器の意義を明確にするため、通常行われる符号語系列に対するＬＬＲ算出方法について説明する。図１７は通常のターボ復号における符号化前の情報系列（Ｉ）にたいするＬＬＲの算出過程を示したトレリス線図である。図１７において、情報系列に対するＬＬＲ（Ｉ）は、
【数２０６】
0004803384.tif 000207
というようにして計算される。【０２１５】
上記の式における条件式、
【数２０７】
0004803384.tif 000208
に相当する箇所を図１７の実線で、条件式、
【数２０８】
0004803384.tif 000209
に相当する箇所を図１７の点線で示す。【０２１６】
図１８は符号化後のパリティ系列（Ｐ）に対するＬＬＲの計算過程を示したトレリス線図である。図１８において、パリティ系列に対するＬＬＲ（Ｐ）は、
【数２０９】
0004803384.tif 000210
というようにして計算される。【０２１７】
上記の式における条件式、
【数２１０】
0004803384.tif 000211
に相当する箇所を図１８の実線で、条件式、
【数２１１】
0004803384.tif 000212
に相当する箇所を図１８の点線で示す。【０２１８】
半環（セミリング）による上記の演算をＬｏｇ領域の半環で行うと、ＬＬＲは理論値通りの値となる。また、Max-Log領域の半環で行うと、近似ではあるが、計算量が少なくてすむ。このようにして得られたＬＬＲ（Ｉ）とＬＬＲ（Ｐ）とが符号化後の符号語系列のＬＬＲとなる。【０２１９】
しかしながら、これらのトレリス線図に基づく計算は軟入力軟出力復号器内部のトレリス上遷移確率伝搬を演算しているコアブロック内部で行われている。従って、機能追加を行うには既製のコアブロックの変更を余儀なくされる。【０２２０】
また、Max-Log領域の半環を使った場合、Ｌｏｇ領域よりＬＬＲの値が高くなる性質がある。そこで、特に図中に記載していないが、Extrinsic情報に重み付けを行い、元に戻して次段へのアプリオリ（a  priori）入力としている。この時、実装のしやすさからシフト加算で重み付けを行っている。【０２２１】
本発明の場合、１／２＋１／４＝０．７５の重み付けを行っている。１／２は１ビットシフトで、１／４は２ビットシフトで実現することができるので、乗算器を使わずに実現することができる。【０２２２】
また、上記の説明では、式の煩雑さを防ぐため、アプリオリ（a  priori）を省略して説明している。半環上の演算では、積の形で対象ビットのメトリックに加わる。特に、但し書きがない場合には、各メトリックにアプリオリが含まれているとして記載している。【０２２３】
このように、本発明では、上述した空間多重信号検出方法と時空間反復復号器とによって、空間多重信号検出方法では必要最小限のサンプル数で最大の効果を得ようというものである。サンプル数が少ないということは条件付確率によってグループ分けされた次段の複雑度をそれだけ軽減することができるということである。【０２２４】
この場合には、上記のようにサンプル数を絞っても、対象ビットに対する最大条件付確率とその排反事象に対する最大条件確率とをメトリックベースで演算しているので、対数尤度代数（Log  Likelihood  Algebra）にしたがったＬＬＲを軟判定出力として提供することができる。すなわち、空間多重信号分離を行う軟入力軟出力検出器とターボ復号器やＬＤＰＣといった軟入力軟出力復号器の連接に際し、Extrinsic情報の受け渡しが行われるが、その元となる正しいＬＬＲを提供することによって、ターボ原理に基づく時空間反復復号の特性を向上させるものである。【０２２５】
つまり、本発明は、空間多重信号分離における軟入力軟出力検出器で尤度を最大化する送信系列推定に半環（セミリング）を用いたメトリック演算手法を用いることによって、最適検出であるＭＬＤに近い高性能を低複雑度で実現するものである。【０２２６】
また、本発明では、尤度を最大化する送信系列を推定するために用いられる半環（セミリング）を用いたメトリック演算手法として、和を最大値演算（ＭＡＸ）、積を通常の加算の
【数２１２】
0004803384.tif 000213
のごとく行うと（Max-log領域における半環）、対象ビットの軟判定出力として複数のグループに分割した条件付確率の各グループ間における前段最大条件付確率事象に対する対象ビットを含む最大条件付確率とその排反事象を含む最大条件付確率とのメトリックベースの差分と、それ以降の同最大条件付確率の差分との総和として表すことができるので、所要システムスループットに対応するため、アンテナ本数を増加させる場合でも、複雑度を抑えた形で対応することができる。あるいは、予め簡易推定された複数個の事象（送信系列）を条件とした場合でも、その条件を基に対象ビットを含む条件付確率とその排反事象を含む条件確率のメトリックベースの差分とそれ以降の条件付き確率の差分の総和として表すことができる。【０２２７】
また、本発明では、空間多重信号検出の軟入力軟出力検出器とターボ復号器やＬＤＰＣといった軟入力軟出力復号器の連接に際し、Extrinsic情報の受け渡しが行われるが、従来行われている外符号用復号器に相当する軟入力軟出力復号器においてExtrinsic情報抽出のしにくさからIntrinsic情報を含んだままExtrinsic情報として使われていたり、或いは一部のExtrinsic情報で賄われたりしていた問題を解消しターボ原理に基づく反復回数による特性向上効果的に実現したものである。【０２２８】
さらに、本発明では、既製の軟入力軟出力復号器を用いる場合、情報ビット系列に対する対数尤度比を出力としており、符号語（シンボル系列）に対するＬＬＲが出力されていない場合が多い。【０２２９】
一方、軟入力軟出力検出器は、符号語に対してＭＬＤ検出となるべく作用するので、軟入力軟出力復号器から軟入力軟出力検出器へのExtrinsic情報として符号語に対するものが必要となる。その結果、軟入力軟出力復号器内部のトレリス上の遷移確率伝搬を演算している箇所を変更することになり、既製コアブロックの変更を余儀なくされるといった問題がある。これに対し、本発明では、軟入力軟出力符号器を用いた時空間反復復号器を用いることで、既製コアブロックの変更をすることなく、ターボ原理に基づく時空間反復復号を実現することができる。【０２３０】
以上、纏めると、本発明では、低複雑化と特性向上との要因として空間多重信号分離のための軟入力軟出力検出器における条件付確率を、複数の条件付確率の積の形でグループ分けし、各グループ間の処理順序付が可能であるので、複雑度の軽減並びにグループ間の順序付けを条件付確率の高い順にすることによって、低複雑化と特性向上とを図ることができる（ファクタライゼーション、オーダリング）。【０２３１】
さらに、本発明では、送信系列を推定するための半環（セミリング）を用いたメトリック演算手法による低複雑度化と、Max-Log領域における半環（セミリング）として、前段までの最大事象を考慮した対象ビットが存在するグループ内最大条件付確率の検出と、後段への最大条件付確率、並びにその排反事象を含む最大条件確率をメトリックベースの差分として各段の差分の総和をとることによって、ＬＬＲ演算をメトリックベースで計算することができるようになるので、対数尤度代数（Log  Llikelihood  Algebra）にしたがったＬＬＲを軟判定出力として提供することができる。
【図１】
0004803384.tif 000218
【図２】
0004803384.tif 000219
【図３】
0004803384.tif 000220
【図４】
0004803384.tif 000221
【図５】
0004803384.tif 000222
【図６】
0004803384.tif 000223
【図７】
0004803384.tif 000224
【図８】
0004803384.tif 000225
【図９】
0004803384.tif 000226
【図１０】
0004803384.tif 000227
【図１１】
0004803384.tif 000228
【図１２】
0004803384.tif 000229
【図１３】
0004803384.tif 000230
【図１４】
0004803384.tif 000231
【図１５】
0004803384.tif 000232
【図１６】
0004803384.tif 000233
【図１７】
0004803384.tif 000234
【図１８】
0004803384.tif 000235