(57)【要約】送信シンボルの離散的な信号点配置を一様連続的な分布で近似することによって、Bayesian型ＭＭＳＥ等化による送信シンボルの推定値Ｘ＾の算出に要する総和演算の回数を低減して演算量を低減すると共に特性劣化を抑制する。送信ベクトルを構成する各送信シンボルについて行うBayesian型ＭＭＳＥ等化における各送信シンボルについての総和演算において、送信シンボルの内の一つの送信シンボルについて、離散的な信号点配置を一様連続な分布で近似することによって当該送信シンボルの総和演算に代えて積分を含まない閉形式とし、送信シンボルの内の近似を行わない残りの送信シンボルについて総和演算を行う。
2012066981.tif 000002
【特許請求の範囲】
【請求項１】
ＭＩＭＯ空間多重通信の受信において、Bayesian型ＭＭＳＥ等化（最小平均二乗誤差）によって信号分離を行う信号分離方法であって、
送信ベクトルを構成する各送信シンボルについて行うBayesian型ＭＭＳＥ等化の総和演算において、
当該送信シンボルの内の一つの送信シンボルについては、離散的な信号点配置を一様連続な分布で近似することによって当該送信シンボルの総和演算に代えて積分を含まない閉形式とし、
前記送信シンボルの内の近似を行わない残りの送信シンボルについては、総和演算を行うことを特徴とする、信号分離方法。
【請求項２】
前記一様連続な分布は、信号点配置を表す複素平面上で、実軸および虚軸においてGauss分布であることを特徴とする、請求項１に記載の信号分離方法。
【請求項３】
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする少なくとも一つの円上において一様連続な分布であることを特徴とする、請求項１に記載の信号分離方法。
【請求項４】
前記送信シンボルは、Ｍ−ＰＳＫの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする単位円上において一様連続な分布であることを特徴とする、請求項３に記載の信号分離方法。
【請求項５】
前記送信シンボルは、Ｍ−ＡＰＳＫの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とし、径を異にする複数の各円上において一様連続な分布であることを特徴とする、請求項３に記載の信号分離方法。
【請求項６】
前記送信シンボルは、ＱＡＭの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする正方形領域内において一様連続な分布であることを特徴とする、請求項１に記載の信号分離方法。
【請求項７】
Bayesian型ＭＭＳＥ等化の期待値演算を事後確率密度関数で表したときの演算式において、前記事後確率密度関数を一様連続な分布関数によって近似することによって、前記送信シンボルの離散的な信号点配置を一様連続な分布で近似することを特徴とする、請求項１に記載の信号分離方法。
【請求項８】
前記分布関数はGauss関数であることを特徴とする、請求項７に記載の信号分離方法。
【請求項９】
前記分布関数は、円上において一様連続な分布を表す関数であることを特徴とする、請求項７に記載の信号分離方法。
【請求項１０】
前記分布関数は、正方形領域内において一様連続な分布を表す関数であることを特徴とする、請求項７に記載の信号分離方法。
【請求項１１】
前記一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行う送信シンボルであることを特徴とする、請求項１から１０の何れか一つに記載の信号分離方法。
【請求項１２】
前記一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行わない送信シンボルであることを特徴とする、請求項１から１０の何れか一つに記載の信号分離方法。
【請求項１３】
ＭＩＭＯ空間多重通信の受信において、Bayesian型ＭＭＳＥ等化（最小平均二乗誤差）によって信号分離を行う信号分離装置であって、
Bayesian型ＭＭＳＥ等化を行う演算部を備え、
当該演算部は、
送信ベクトルを構成する各送信シンボルについて行うBayesian型ＭＭＳＥ等化における各送信シンボルについての総和演算において、
当該送信シンボルの内の一つの送信シンボルについて、離散的な信号点配置を一様連続な分布で近似することによって当該送信シンボルの総和演算に代えて積分を含まない閉形式とし、
前記送信シンボルの内の近似を行わない残りの送信シンボルについて総和演算を行うことを特徴とする、信号分離装置。
【請求項１４】
前記一様連続な分布は、信号点配置を表す複素平面上で、実軸および虚軸においてGauss分布であることを特徴とする、請求項１３に記載の信号分離装置。
【請求項１５】
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする少なくとも一つの円上において一様連続な分布であることを特徴とする、請求項１３に記載の信号分離装置。
【請求項１６】
前記送信シンボルは、Ｍ−ＰＳＫの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする単位円上において一様連続な分布であることを特徴とする、請求項１５に記載の信号分離装置。
【請求項１７】
前記送信シンボルは、Ｍ−ＡＰＳＫの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とし、径を異にする複数の各円上において一様連続な分布であることを特徴とする、請求項１５に記載の信号分離装置。
【請求項１８】
前記送信シンボルは、ＱＡＭの信号点配置であり、
前記一様連続な分布は、信号点配置を表す複素平面上で、原点を中心とする正方形領域内において一様連続な分布であることを特徴とする、請求項１３に記載の信号分離装置。
【請求項１９】
Bayesian型ＭＭＳＥ等化の期待値演算を事後確率密度関数で表したときの演算式において、前記送信シンボルの離散的な信号点配置を一様連続な分布で近似することを特徴とする、請求項１３に記載の信号分離装置。
【請求項２０】
前記分布関数はGauss関数であることを特徴とする、請求項１９に記載の信号分離装置。
【請求項２１】
前記分布関数は、円上において一様連続な分布を表す関数であることを特徴とする、請求項１９に記載の信号分離装置。
【請求項２２】
前記分布関数は、正方形領域内において一様連続な分布を表す関数であることを特徴とする、請求項１９に記載の信号分離装置。
【請求項２３】
前記一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行う送信シンボルであることを特徴とする、請求項１３から２２の何れか一つに記載の信号分離装置。
【請求項２４】
前記一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行わない送信シンボルであることを特徴とする、請求項１３から２２の何れか一つに記載の信号分離装置。
【請求項２５】
前記演算部は、前記事後確率密度関数を一様連続の分布で近似する複数種の分布関数を備え、送信シンボルの分布に応じて選択して適用することを特徴とする、請求項１９に記載の信号分離装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号分離方法および信号分離装置に関し、特にＭＩＭＯ（Multi-Input Multi-Output）による空間多重通信(ＳＭ: spatial multiplexing)の受信側で使用される信号分離方法および信号分離装置に関する。
【背景技術】
【０００２】
ＭＩＭＯ(multiple-input multiple-output)空間多重 (ＳＭ: spatial multiplexing)は複数の送信アンテナから異なる情報を送ることで、帯域の拡張を必要とせずにスループットを倍増できる技術である。ＭＩＭＯ−ＳＭによる通信は異なる信号が伝搬路で線形に重ね合わされて受信されるため、受信した信号を適切に分離することが必要である。信号分離法の精度はＭＩＭＯ−ＳＭのスループット増大効果を決定する主要因となるため、マルチアンテナ技術の広がりと共に様々な手法が提案されている（例えば、非特許文献１）。【０００３】
最も単純な信号分離法として、ＭＭＳＥ (minimum mean square error)やＺＦ(zero forcing)規範に基づく線形等化が挙げられる。線形等化による信号分離法は、線形処理のみで信号分離を行う簡易な方法であり、単純に実装できるという利点がある反面、達成できるダイバーシティオーダに制限があり、受信アンテナ数が少ない場合は非常に特性が悪い（例えば、非特許文献２）。【０００４】
逐次干渉キャンセラ(ＳＩＣ: successive interference canceller)はＭＩＭＯ−ＳＭの先駆的検討となったV-BLASTで用いられている構成であり、線形等化と干渉信号レプリカの除去を繰り返すことで信号分離の精度を高める方法である。ＳＩＣは線形等化に若干の処理量の追加で実装できるが、判定帰還を必要とし、複数のシンボルを順に処理していくため同時処理ができず遅延が問題となる。また、良い特性のためには、はじめに検出するシンボルの精度を高める工夫や誤り訂正を含めた判定帰還が必要とされる（例えば、非特許文献２）。 【０００５】
近年では、ＬＲ (lattice reduction)法と線形等化や ＳＩＣの組み合わせ（例えば、非特許文献３，４，５）が知られている。 【０００６】
ＬＲ法は、信号点の格子を通信路行列が好ましい形になるように変換する手法であり、変換された格子上でシンボルの判定を実行することで、通信路行列が悪条件 (ill conditioned)の場合でも線形等化やＳＩＣで安定した信号分離が可能となる。しかし、変換された格子上でシンボルの判定をしなければならないために軟判定値の出力が困難であり、シンボル点を硬判定しなければならないため、誤り訂正符号との相性が悪いという問題がある。さらに、ＰＳＫのように格子構造に基づかない信号点配置には適用できないといった欠点がある。【０００７】
また、ＬＲ法はＱＲ分解よりもさらに複雑であり、さらに、適用できる信号点配置はＱＡＭに限られる。【０００８】
他の信号分離法として、低演算量化された最尤 (ＭＬ: maximum likelihood)検出や最大事後確率 (ＭＡＰ: maximum a posteriori)検出が挙げられる。最尤検出（ＭＬＤ）は、誤り率を最小化できるが、送信アンテナ数をＮｔ、変調多値数をＭとするとＭＮｔのオーダーという膨大な計算量を必要とするという問題がある。【０００９】
また、通信路行列をＱＲ分解した後に、Ｍアルゴリズムやsphere decodingを用いてＭＬ検出を近似する方法が知られている（非特許文献２）。【００１０】
ＱＲ分解とＭアルゴリズムとを組み合わせた方法は、ＱＲ分解を施した通信路行列に基づいて実装可能な計算量でＭＬＤを近似するものであるが、ＳＩＣと同様に逐次処理であるため処理遅延が問題となる。各段階で残す候補数をＳとすると、Ｓ×Ｍの計算をＮt回繰り返さなければならず、計算量を低減するためにＳを小さくすると特性が劣化する。さらに、フェージング係数が変化する度にＱＲ分解の実行が必要である。 【００１１】
ＱＲ分解とSphere decodingとを組み合わせた方法は、ＱＲ分解とＭアルゴリズムとによる方法と同様に実装可能なＭＬＤの近似であり、上記した問題の他に、残す候補の判定が難しいという問題がある。【００１２】
また、ＱＲ分解後にＳＭＣ (sequential Monte-Carlo)法でシンボルの事後確率を計算し、実現可能な計算量でＭＡＰ推定を行う方法も提案されている（例えば、非特許文献６）。計算量はＮｔに対し線形のオーダーであるが、各ステージにおける処理は比較的複雑である。【００１３】
上記したように、ＱＲ分解に基づく方法は逐次処理であるために、各アンテナからの信号の検出を同時処理できないという問題があり、また、演算量は指数関数的増加を示さないものの、１シンボルの検出当たりに要する処理は依然として大きいという問題がある。【００１４】
また、Bayesian型ＭＭＳＥ等化による信号分離方法が知られている（例えば、非特許文献７，８）。【００１５】
ＭＭＳＥ規範は元の送信シンボルと推定値との平均二乗誤差(ＭＳＥ)を最小化するものであり、信号分離処理の処理を線形演算に制限しない場合は、次式（１）で表される推定値によって最小ＭＳＥが達成される。【００１６】
2012066981.tif 000003
【００１７】
なお、式では推定値はＸの上部に“＾”の記号を付しているが、以下の文書中では、Ｘ＾によって表記する。【００１８】
この推定値は、受信ベクトルがＹ＝ｙとなるような送信ベクトルＸの平均値を推定値とすることを意味しており、ＭＳＥが最小化される。式 (１)で定義される信号分離法をBayesian型 ＭＭＳＥ等化と呼ぶ（非特許文献７，８）。【００１９】
式(１)において、期待値演算を事後pdfpX|Y (x|y)を用いて書き下せば、次式（２）となる。【００２０】
2012066981.tif 000004
【００２１】
式（２）は、ＭＡＰ検出のように事後pdfを最大化するものを検索するのではなく、全通りのシンボルパターンに対応する事後pdfで重み付けすることを意味している。その点で、Bayesian型ＭＭＳＥ等化は軟判定のＭＡＰ検出と見なすことができる。Bayesian型ＭＭＳＥ等化の推定結果Ｘを硬判定すれば、ＭＡＰ検出と等価となる。また、送信するシンボル点が等確率で選ばれるならば、ＭＬ検出とも等価である。 
【先行技術文献】
【非特許文献】
【００２２】
【非特許文献１】E.G.Larsson,“MIMO detection methods: How they work,” IEEE Sig. Proc. Mag., vol. 26, pp. 91. 95, May 2009.
【非特許文献２】樋口，田岡 ,“マルチアンテナ無線伝送技術その3 MIMO多重法における信号分離技術,” NTT DoCoMoテクニカル・ジャーナル ,  vol. 14, pp. 66.75, 4月 2006.
【非特許文献３】H. Yao and G. W. Wornell, “Lattice-reductionaided detectors for MIMO  communications systems”, in Proc. IEEE Globecom'02, pp. 424.428, Nov. 2002.
【非特許文献４】C. Windpassinger and R. F. H. Fischer, “Lowcomplexity near-maximum-likelihood detection and precoding for MIMIO systems using lattice reduction”,; in Proc. ITW'03, Apr. 2003.
【非特許文献５】D. W¨ubben, R. B¨ohnke, V. K¨uhn, and K. D. Kammeyer,“Near-maximum-likelihood detection of MIMO systems using MMSE-based lattice reduction,” in Proc. IEEE ICC'04, pp. 798.802, June 2004.
【非特許文献６】B. Dong, X. Wang, and A. Doucet, “A new class of soft MIMO demodulation algorithms,” IEEE Trans. Sig. Proc., vol. 51, pp. 2752.2763, Nov. 2003.
【非特許文献７】S. Chen, S. McLaughlin, B. Mulgrew, and P. M. Grant, “Adaptive Bayesian decision feedback equalizer for dispersive mobile radio channels,”; IEEE Trans. Commun., vol. 43, pp. 1937.1946, May 1995.
【非特許文献８】S. K. Patra and B. Mulgrew, “Efficient architecture for Bayesian equalization using fuzzy filters,” IEEE Trans. Circuits Syst. II, Analog Digit. Sig. Proc., vol. 45, pp. 812.820, July 1998.
【発明の概要】
【発明が解決しようとする課題】
【００２３】
ＭＩＭＯ空間多重通信は、複数のアンテナから送信された信号を受信側において分離することが必要であり、信号分離の精度を高めるにはその計算量が膨大となるという問題がある。【００２４】
Bayesian型ＭＭＳＥ等化による信号分離方法において、式（２）で表される定義式をそのまま実行すると、ＭＬやＭＡＰ検出と同様にＭＮｔ個のシンボルパターンを検査しなければならない。【００２５】
例えば送信ベクトルが２つの送信シンボルｘ１、ｘ２である場合には、式（２）の送信シンボルＸ＾は以下の式（３）で表される。【００２６】
2012066981.tif 000005
【００２７】
上記式（３）は、雑音の分散をＮ０とすると、ベイズの定理によって次式（４）に変形される。【００２８】
2012066981.tif 000006
【００２９】
ただし、λ（ｘ１、ｘ２）は次式（５）で表される。式（５）において、“：＝”は定義を表す等号である。【００３０】
2012066981.tif 000007
【００３１】
上記式から、送信シンボルの個数が２の場合にはＭ２パターンの［ｘ１ ｘ２］について２乗距離を計算する必要があり、送信シンボルの個数がＮｔの場合には、ＭＮｔパターンの［ｘ１ ｘ２ …ｘＮｔ］について２乗距離を計算する必要がある。【００３２】
式（２）において、推定値Ｘ＾を算出するに要する総和演算の回数は、送信シンボルのパターン数と同数が必要である。【００３３】
図１４は、上記したBayesian型ＭＭＳＥ等化による信号分離の演算手順を説明する図であり、２つの送信シンボルｘ１、ｘ２の例について示している。【００３４】
式（４）において、はじめに一方の送信シンボルｘ１について総和演算を行う。ここでは、式（４）の分子において送信シンボルｘ１について行った総和演算をＡとし、式（４）の分母において送信シンボルｘ１について行った総和演算をＢとしている（S100）。【００３５】
次に、他方の送信シンボルｘ２について総和演算を行う。ここでは、式（４）の分子において送信シンボルｘ２について行った総和演算をＣとし、式（４）の分母において送信シンボルｘ２について行った総和演算をＤとしている（S101）。送信シンボルの推定値Ｘ＾は、Ｃ／Ｄの演算によって求める（S102）。【００３６】
したがって、本発明は上記課題を解決して、Bayesian型ＭＭＳＥ等化による信号分離の演算量を低減することを目的とする。
【課題を解決するための手段】
【００３７】
本願発明は、送信シンボルの離散的な信号点配置を一様連続的な分布で近似することによって、Bayesian型ＭＭＳＥ等化による送信シンボルＸ＾の算出に要する総和演算の回数を低減して演算量を低減すると共に、特性劣化を抑制するものである。【００３８】
本発明は、ＭＩＭＯ空間多重通信の受信において、Bayesian型ＭＭＳＥ等化（最小平均二乗誤差）によって信号分離を行う信号分離であり、信号分離方法の形態と信号分離装置の形態とすることができる。【００３９】
（本発明の信号分離方法の形態）
本発明の信号分離方法の形態は、送信ベクトルを構成する各送信シンボルについて行うBayesian型ＭＭＳＥ等化の総和演算において、送信シンボルの内の一つの送信シンボルについては、離散的な信号点配置を一様連続な分布で近似することによって当該送信シンボルの総和演算に代えて積分を含まない閉形式とし、送信シンボルの内の近似を行わない残りの送信シンボルについては、総和演算を行う。【００４０】
一つの送信シンボルの離散的な信号点配置を一様連続な分布で近似することによって、各送信シンボルについて行う総和演算の内の一つの総和演算を積分を含まない閉形式に置き換え、これによって総和演算の回数を低減し、Bayesian型ＭＭＳＥ等化による信号分離の演算量を低減する。【００４１】
一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行う送信シンボル、あるいは、送信シンボルの内で推定を行わない送信シンボルとすることができる。【００４２】
本発明の信号分離方法において、送信シンボルの離散的な信号点配置を近似する一様連続な分布は、種々の分布を適用することができる。【００４３】
一様連続な分布の一例は、信号点配置を表す複素平面上で、実軸および虚軸においてGauss分布であり、一様連続な分布の他の例は、信号点配置を表す複素平面上において原点を中心とする少なくとも一つの円上において一様連続な分布である。【００４４】
円上の一様連続な分布の例は、信号点配置を表す複素平面上で、原点を中心とする単位円上において一様連続な分布であり、Ｍ−ＰＳＫの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００４５】
また、円上の一様連続な分布の別の例は、信号点配置を表す複素平面上で、原点を中心とし径を異にする複数の各円上において一様連続な分布であり、Ｍ−ＡＰＳＫの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００４６】
さらに、一様連続な分布の別の例は、信号点配置を表す複素平面上において面状の一様連続な分布であり、原点を中心とする正方形領域内において一様連続な分布とすることができ、ＱＡＭの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００４７】
本発明は、Bayesian型ＭＭＳＥ等化の期待値演算を事後確率密度関数で表したときの演算式において、送信シンボルの離散的な信号点配置を一様連続な分布で近似する。【００４８】
分布関数は、Gauss関数、円上において一様連続な分布を表す関数、正方形領域内において一様連続な分布を表す関数とすることができる。円上において一様連続な分布を表す関数は、Ｍ−ＰＳＫの離散的な信号点配置や、Ｍ−ＡＰＳＫの離散的な信号点配置を一様連続な分布で近似する際に適用することができ、正方形領域内において一様連続な分布を表す関数は、ＱＡＭの離散的な信号点配置の一様連続な分布で近似する際に適用することができる。【００４９】
（本発明の信号分離装置の形態）
本発明の信号分離装置の形態は、送信ベクトルを構成する各送信シンボルについて行うBayesian型ＭＭＳＥ等化の総和演算を行う演算部を備え、当該演算部は、送信シンボルの内の一つの送信シンボルについては、離散的な信号点配置を一様連続な分布で近似することによって当該送信シンボルの総和演算に代えて積分を含まない閉形式とし、送信シンボルの内の近似を行わない残りの送信シンボルについては、総和演算を行う。【００５０】
演算部は、一つの送信シンボルの離散的な信号点配置を一様連続な分布で近似することによって、各送信シンボルについて行う総和演算の内の一つの総和演算を積分を含まない閉形式に置き換え、これによって総和演算の回数を低減し、Bayesian型ＭＭＳＥ等化による信号分離の演算量を低減する。【００５１】
一様連続な分布で近似する送信シンボルは、送信シンボルの内で推定を行う送信シンボル、あるいは、送信シンボルの内で推定を行わない送信シンボルとすることができる。【００５２】
本発明の信号分離装置において、送信シンボルの離散的な信号点配置を近似する一様連続な分布は、種々の分布を適用することができる。【００５３】
一様連続な分布の一例は、信号点配置を表す複素平面上で、実軸および虚軸においてGauss分布であり、一様連続な分布の他の例は、信号点配置を表す複素平面上において原点を中心とする少なくとも一つの円上において一様連続な分布である。【００５４】
円上の一様連続な分布の例は、信号点配置を表す複素平面上で、原点を中心とする単位円上において一様連続な分布であり、Ｍ−ＰＳＫの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００５５】
また、円上の一様連続な分布の別の例は、信号点配置を表す複素平面上で、原点を中心とし径を異にする複数の各円上において一様連続な分布であり、Ｍ−ＡＰＳＫの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００５６】
さらに、一様連続な分布の別の例は、信号点配置を表す複素平面上において面状の一様連続な分布であり、原点を中心とする正方形領域内において一様連続な分布とすることができ、ＱＡＭの離散的な信号点配置の送信シンボルを一様連続な分布に近似する場合に適用することができる。【００５７】
本発明の信号分離装置の演算部は、Bayesian型ＭＭＳＥ等化の期待値演算を事後確率密度関数で表したときの演算式において、送信シンボルの離散的な信号点配置を一様連続な分布で近似する。【００５８】
分布関数は、Gauss関数、円上において一様連続な分布を表す関数、正方形領域内において一様連続な分布を表す関数とすることができる。円上において一様連続な分布を表す関数は、Ｍ−ＰＳＫの離散的な信号点配置や、Ｍ−ＡＰＳＫの離散的な信号点配置を一様連続な分布で近似する際に適用することができ、正方形領域内において一様連続な分布を表す関数は、ＱＡＭの離散的な信号点配置の一様連続な分布で近似する際に適用することができる。【００５９】
演算部は、近似を行う分布関数を複数種備えておき、送信シンボルの形態に応じて選択して適用することができる。
【図面の簡単な説明】
【００６０】
【図１】ＭＩＭＯ通信モデルおよび本発明の信号分離装置を説明するための図である。
【図２】本発明による信号分離の概略を説明するための図である。
【図３】Ｍ−ＰＳＫの変調および近似分布を説明するための図である。
【図４】Ｍ−ＱＡＭの変調および近似分布を説明するための図である。
【図５】Ｍ−ＡＰＳＫの変調および近似分布を説明するための図である。
【図６】Uniform ring近似によるＱＰＳＫのビット誤り率特性を示すグラフである。
【図７】Uniform ring近似による８−ＰＳＫのビット誤り率特性を示すグラフである。
【図８】Uniform ring近似による１６−ＰＳＫのビット誤り率特性を示すグラフである。
【図９】Uniform ring近似による８−ＰＳＫのビット誤り率特性を示すグラフである。
【図１０】Uniform ring近似による１６−ＡＰＳＫのビット誤り率特性を示すグラフである。
【図１１】Uniform square近似による１６−ＱＡＭのビット誤り率特性を示すグラフである。
【図１２】Gauss近似による８−ＰＳＫの特性およびUniform ring近似との比較を示すグラフである。
【図１３】Gauss近似による１６−ＱＡＭの特性およびUniform square近似との比較を示すグラフである。
【図１４】Bayesian型ＭＭＳＥ等化による信号分離の演算手順を説明する図である。
【発明を実施するための形態】
【００６１】
以下、本発明の実施の形態について、図を参照しながら詳細に説明する。
【００６２】
Bayesian型 ＭＭＳＥ等化によれば、定義式のまま実行するとＭＬ検出とほぼ同じ特性を達成できるものの同等の計算量が必要となる。本発明は、このBayesian型 ＭＭＳＥ等化において、送信シンボルの離散的な信号点配置を一様連続な分布で近似することによって、分離処理の計算量を低減する。
【００６３】
本発明は、変調形式としてＭ−ＰＳＫ、Ｍ−ＡＰＳＫ、ＱＡＭに適用することができ、また、送信アンテナ数も任意の個数とすることができる。例えば、信号点配置がＰＳＫである場合には、信号分離の計算量をＭＮｔオーダーからＭＮｔ−１のオーダーに低減することができる。計算量は依然として指数関数的のままであるが、冪数を１低減することができる。
【００６４】
例えば、Ｎｔ＝２に限定した場合には、計算量は線形のオーダーとなる。小型端末はサイズの制限から搭載可能なアンテナ数が少なく、バッテリの節約から電力増幅器の非線形性に耐性のあるＰＳＫ変調が用いられる事情を考慮すると、送信アンテナの個数Ｎｔを２としＰＳＫ変調とする通信形態の設定は、小型端末から基地局への伝送において十分に想定できるものである。 
【００６５】
図１は、ＭＩＭＯ通信モデルおよび本発明の信号分離装置を説明するための図である。ＭＩＭＯ通信モデルは、送信アンテナ１，伝送路２，受信アンテナ３、チャンネル推定部４，および検出器５を備える。
【００６６】
ＭＩＭＯ通信モデルでは、送信アンテナ１から別々の送信ベクトルＸが同時に同一の周波数で送信される。これらの送信ベクトルＸは、伝送路２を経て複数の受信アンテナで受信される。検出器５では、各受信アンテナ３で受信された受信ベクトルＹに基づいて、送信された複数の送信ベクトルＸを個々の信号に分離する。
【００６７】
チャンネル推定部４は、パイロット信号を検出することによって伝送路２の特性を表わす通信路行列Ｈを推定する。検出器５は、チャンネル推定部４で推定した通信路行列Ｈを用いて、信号分離を行う。
【００６８】
本発明の検出器５は、受信した送信ベクトルを構成する各送信シンボルについて、Bayesian型ＭＭＳＥ等化の演算を行う演算部５Ａ、および演算部５Ａの演算結果に基づいて送信シンボルの推定値を算出する算出部５Ｂを備える。演算部５Ａは、Bayesian型ＭＭＳＥ等化の演算で行う各送信シンボルの総和演算において、送信シンボルの内の一つの送信シンボルについては、離散的な信号点配置を一様連続な分布で近似することによってその送信シンボルの総和演算に代えて積分を含まない閉形式とし、送信シンボルの内の近似を行わない残りの送信シンボルについては、総和演算を行うことによって演算量を低減した信号分離を行う。
【００６９】
演算部５Ａは、近似に用いる近似分布関数を記憶手段に格納し、送信ベクトルの変調方式に応じて記憶手段から適応する近似分布関数を読み出して、信号分離の演算に用いることができる。記憶手段は、演算部５Ａ内に設ける構成とする他、別に設ける構成としてもよい。また、近似分布関数は、記憶手段に予め格納する他、図示しない入力部によって格納する構成とすることもできる。
【００７０】
以下、大文字（Ｘ１，Ｘ２，Ｙ１，…）を確率変数、その取り得る値（実現値）を小文字（ｘ１，ｘ２，ｙ１，…）で表現する。
【００７１】
送信アンテナ１のアンテナ数をＮｔ、受信アンテナ３のアンテナ数をＮｒとし、図１は送信アンテナ１のアンテナ数をＮｔが２である例を示している。各送信アンテナ１は送信シンボル［Ｘ1、…、ＸＮt］を有する送信ベクトルを送信する。受信アンテナ３は、ＭＩＭＯチャネルＨを伝搬した受信ベクトルＹを受信する。受信ベクトルＹは受信シンボル［Ｙ1、…、ＹＮｒ］で表される。
【００７２】
通信路の入出力関係は次式（６）で表される。
【００７３】
2012066981.tif 000008
【００７４】
ここで、送信ベクトルＸ、受信ベクトルＹ、雑音ベクトルＷ、通信路行列Ｈは、以下の式（７）〜（１０）で表される。“：＝”は定義を表す等号である。
【００７５】
2012066981.tif 000009
2012066981.tif 000010
2012066981.tif 000011
2012066981.tif 000012
【００７６】
各送信シンボル Ｘｉ,ｉ＝１，…，Ｎtは互いに独立で同じ分布を持つ確率変数の列（i.i.d.）であり、それぞれＭ要素の信号点の集合χＭから任意の値をとる。雑音ベクトルの各要素も互いに独立で同じ分布を持つ確率変数の列（i.i.d.）であり、それぞれ平均０、分散Ｎ０の複素ガウス分布に従う。
【００７７】
ＭＩＭＯ−ＳＭにおいて、信号分離によって受信ベクトルＹから対応する送信ベクトルの推定値Ｘ＾：＝［Ｘ＾１  …  Ｘ＾Ｎｔ］を求める。
【００７８】
上記したＭＩＭＯ通信モデルにおいて、信号点配置、送信アンテナ数、および受信アンテナ数は任意とすることができる。例えば、図１において、Ｍ要素の信号点の集合χＭをＭ−ＰＳＫ、Ｍ−ＡＰＳＫ，Ｍ−ＱＡＭ等の出力点集合とし、送信アンテナ数をＮｔ＝２とすることができる。
【００７９】
複素信号Ｘ１，Ｘ２を送信してＹ１，Ｙ２，…，ＹＮｒを受信したとすれば、上記式（６）は次式（１１）で表される。
【００８０】
2012066981.tif 000013
【００８１】
ここで、ｈｋ，ｎは送信アンテナｎから受信アンテｋへの伝搬路応答であり、Ｗ１，Ｗ２，…，ＷＮｒは各受信アンテナから入力される信号の雑音成分である。各送信シンボルＸ１，Ｘ２は、それぞれＭ要素の信号点の集合χＭから任意の値をとる。
【００８２】
集合χＭは、例えば、Ｍ−ＱＡＭ(ただしＭは２の偶数乗)、Ｍ−ＰＳＫの場合には、それぞれ次式（１２）、（１３）で表される。
【００８３】
2012066981.tif 000014
2012066981.tif 000015
【００８４】
集合χＭがＭ−ＡＰＳＫの場合には、リング数をＫ、各リングの半径をρ１，…，ρＫ、各リング上の信号点数をＭ（ｋ），ｋ＝１，…,Ｋ（ただし、ΣＫｋ＝１Ｍ（ｋ）＝Ｍ）とすると、次式（１４）で表される。
【００８５】
2012066981.tif 000016
【００８６】
なお、ＭーＰＳＫは、この定義においてＫ＝１，ρ１＝１，Ｍ（１）＝Ｍの場合と見なすことができる
【００８７】
Bayesian型ＭＭＳＥ等化において、式（１）から式（２），式（３）を導く際には、各送信シンボルＸｉの事前確率密度関数pdfは離散的であるとしている。本発明は、送信シンボルＸｉの離散的な値を一様連続的な分布で近似する。
【００８８】
本発明は、式（４）において送信シンボルＸを一様連続な分布f(x)に従う確率変数で近似することによって計算量を低減する。
【００８９】
送信シンボルＸを一様連続な分布で近似する際、推定を行う送信シンボルＸに対して一様連続な分布の近似を適用する態様の他に、推定を行わない送信シンボルＸに対して一様連続な分布の近似を適用する態様とすることができる。以下、推定を行う送信シンボルＸに対して一様連続な分布の近似を適用する態様を第１のタイプとし、推定を行わない送信シンボルＸに対して一様連続な分布の近似を適用する態様を第２のタイプとして説明する。
【００９０】
（第１のタイプ）
第１のタイプの態様では、送信シンボルＸ１を推定する場合において、送信シンボルＸ１を一様連続な分布f(x)で近似することは、推定を行う所望の送信シンボルの分布を近似することを意味する。所望のシンボルを近似して導かれる第１のタイプの推定式は、式（４）を一様連続な分布f(x)を用いて書き換えることによって以下の式（１５）〜(１７)で表される。
【００９１】
2012066981.tif 000017
2012066981.tif 000018
2012066981.tif 000019
ここで、∫ｃは複素平面全体に渡る積分を表している。
【００９２】
関数f(x) には、
１）元の信号点配置をできるだけ忠実に再現していること、および
２)α(x),β(x)で与えられる積分を、積分を含まない形式である閉形式とすること
が求められる。閉形式とすることによって、式（１５）で示されるように、Ｎｔ＝２であるときに通常Ｍ２パターンの［ｘ１ ｘ２］の検査を必要とするBayesian型ＭＭＳＥ等化においてx１についての総和演算が無くなり、Ｍパターンのx２のみの総和演算によって、信号分離の算出を行うことができる。
【００９３】
（第２のタイプ）
送信シンボルＸ１ではなく送信シンボルＸ２に対して一様連続な分布f(x)の近似を適用することも可能である。
【００９４】
推定を行う所望の送信シンボルＸ１からみて送信シンボルＸ２は干渉シンボルに相当する。したがって、第２のタイプの態様は、干渉となる送信シンボルに対して一様連続な分布f(x)の近似を適用することを意味している。
【００９５】
干渉となる送信シンボルを一様連続な分布を表す関数f(x)で近似して導かれる第２のタイプの推定式は、式（４）を一様連続な分布を表す関数f(x)を用いて書き換えることによって以下の式（１８）で表される。
【００９６】
2012066981.tif 000020
【００９７】
ただし、β´(x1)はβ(x2)において 
2012066981.tif 000021
と変数を交換したものに相当する。
【００９８】
第２のタイプは、第１のタイプと同様に計算量を削減することができるが、その特性は第１のタイプ と比較して劣ることが確認されている。
【００９９】
以下、様々な形の一様連続な分布関数f(x)について、第１のタイプおよび第２のタイプの推定式の導出を説明する。
【０１００】
なお、各推定式は、Ｘ＾１ではなくlogＸ＾１を与える式として表記されている。これは計算機で実装する際に、対数領域で計算をすることでオーバーフローを避けるためである。したがって、最終的な推定結果はＸ＾１＝exp（logＸ＾１）として得られる。
【０１０１】
図２は、本発明による信号分離の概略を説明するための図である。図２において、信号分離の演算は、Bayesian型ＭＭＳＥ等化の式において、複数の送信ベクトルの内の一つの送信ベクトルの送信シンボルに係わる総和演算を、一様連続な分布関数f(x)を用いることによって閉形式の積分（α（ｘ），β（ｘ））とし、残りの送信ベクトルの送信シンボルについて総和演算を行い（S1）、得られた演算結果から送信ベクトルの推定値Ｘ＾を算出する（S2）。
【０１０２】
S1の工程では、Bayesian型ＭＭＳＥ等化で行う総和演算の内、送信ベクトルの一つを一様連続な分布関数f(x)で近似することで、この送信ベクトルの各送信シンボルについての総和演算を不要とし、これによって計算量を低減することができる。
【０１０３】
ここで、以下の推定式に共通して現れる変数を定義しておく。 
第１のタイプの推定式で用いられる変数は以下の通りである。
【０１０４】
2012066981.tif 000022
2012066981.tif 000023
2012066981.tif 000024
2012066981.tif 000025
2012066981.tif 000026
2012066981.tif 000027
【０１０５】
第２のタイプの推定式で用いられる変数は以下の通りである。第２のタイプの推定式の元となるβ′(x１)はβ(x２)において式(１９)の変数交換が施されたものに相当する。したがって、第２のタイプにおいては上記した式（２０）〜（２３）の変数を以下のものと読み替えればよい。
【０１０６】
2012066981.tif 000028
2012066981.tif 000029
2012066981.tif 000030
2012066981.tif 000031
【０１０７】
また、推定式において上記した変数の他に、導出過程中のみに現れる最終的な結果に影響を及ぼさない定数をＣで表している。
【０１０８】
対数領域で表記された推定式において、ΣR,log,ΣC,logで表記される総和演算が現れる。これらはそれぞれ、実数領域および複素数領域のヤコビアン対数に基づく総和演算であり、通常の加算を以下の式（３０）、（３０’）、（３１）、（３１’）で示すように、別の演算に置き換えて実行することを意味する。
【０１０９】
2012066981.tif 000032
2012066981.tif 000033
2012066981.tif 000034
2012066981.tif 000035
ここで、maxＲｅ(a,b),minＲｅ(a,b)は、aとbのうち実数部が大きい(小さい)方を表している。 
【０１１０】
式（３０’）および式（３１’）は、式（３０）および式（３１）中の対数演算の項を省略した式である。式（３０’）および式（３１’）によって得られる特性は、後に示す図６〜図１１に示すように、式（３０）および式（３１）によって得られる特性と比較して劣化の程度は低いことが確認される。また、対数演算の項を省略することによって計算量を削減することができる。
【０１１１】
以下、送信ベクトルを一様連続な分布で近似する例として、Gauss分布による近似、ＱＡＭに適用する正方形領域の一様連続分布による近似、ＰＳＫやＡＰＳＫに適用する円（リング）上の一様連続分布による近似の各例について説明する。
【０１１２】
［Gauss分布による近似］
送信ベクトルを近似する一様連続な分布としてガウス分布を用いる。ガウス分布による近似において、一様連続な分布の分布関数として、以下のガウス関数ｆ（ｘ）を用いる。
ｆ(x)=exp(−|x|２)                  … （３２）
【０１１３】
図３（ｄ）は、ガウス分布の一例を示し、ガウス関数ｆ（ｘ）は複数の平面の実軸方向および虚軸方向についてガウス分布するものとする。ガウス分布は、実際のシンボルの分布と大きく異なるが、積分を簡易な閉形式にできる点で好ましい。 
【０１１４】
式（３２）に式（１６），（１７）を代入して計算すると、α、βは以下の式（３３）、（３４）で表される。
【０１１５】
2012066981.tif 000036
2012066981.tif 000037
【０１１６】
（第１のタイプ）
推定する送信ベクトルに対してガウス分布による近似を適用すると、式（１５）に式（３３），（３４）を代入することによって、推定送信シンボルＸ＾１は以下の式（３５）で表される。
【０１１７】
2012066981.tif 000038
【０１１８】
対数領域では、以下の式（３６）で表される。
2012066981.tif 000039
【０１１９】
（第２のタイプ）
推定を行わない送信ベクトルに対してガウス分布による近似を適用すると、式（１８）に式（３３），（３４）を代入することによって、推定送信シンボルＸ＾１は以下の式（３７）で表される。
【０１２０】
2012066981.tif 000040
【０１２１】
対数領域では、以下の式（３８）で表される。
2012066981.tif 000041
【０１２２】
［正方形領域の一様連続分布による近似］
信号点配置がＭ−ＱＡＭ(ただしＭは２の偶数乗)である場合を近似する好適な例として、複素平面内の正方形領域における一様連続分布による近似（uniform square近似）を設定することができる。
【０１２３】
図４（ａ）は１６ＱＡＭの例を示し、図４（ｂ）は正方形領域における一様連続分布の近似例を示している。
【０１２４】
このuniform square近似は、格子状に離散的に分布するＱＡＭの信号点配置を、信号点数が無限大のＱＡＭ(∞−ＱＡＭ)で近似するものであり、原点を中心とし、頂点がκ+jκ,−κ+jκ,−κ−jκ,κ−jκの正方形領域Ｄκに、無数のシンボル点が一様連続的に分布している。
【０１２５】
このuniform square近似を表す分布関数f(x)は、ユニットステップ関数U(・)を用いて以下の式（３９）で表される。
【０１２６】
2012066981.tif 000042
【０１２７】
式（３９）において、κは、ＭーＱＡＭの信号点の集合においてＭを無限大としたときの極限を求めることで求めることができ、κ＝√3/2となる。
【０１２８】
上記式（３９）で表される分布関数ｆ（ｘ）を式（１７）、（１６）に代入することによって、式（４０）、（４１）で表されるβおよびαが得られる。
【０１２９】
2012066981.tif 000043
2012066981.tif 000044
【０１３０】
ここで、
2012066981.tif 000045
である。 
【０１３１】
（第１のタイプ）
推定する送信ベクトルに対してuniform square近似による近似を適用すると、式（１５）に式（４０），（４１）を代入することによって、推定送信シンボルＸ＾１は以下の式（４６）で表される。なお、式（４１）において、
2012066981.tif 000046
および
2012066981.tif 000047
は非常に小さいので無視している。
【０１３２】
2012066981.tif 000048
【０１３３】
式（４６）において、ｅＩ，ｅＱは非負であるので、式（４６）は以下の式（４７）の対数領域で計算することができる。
【０１３４】
2012066981.tif 000049
【０１３５】
（第２のタイプ）
推定を行わない送信ベクトルに対してuniform square近似を適用すると、式（１８）に式（３３）を代入することによって、推定送信シンボルＸ＾１は以下の式（４８）で表される。
【０１３６】
2012066981.tif 000050
【０１３７】
対数領域では、以下の式（４９）で表される。
2012066981.tif 000051
【０１３８】
［円（リング）上の一様連続分布による近似］
信号点配置がＭ−ＡＰＳＫ、Ｍ−ＰＳＫである場合を近似する好適な例として、複素平面内において各円（リング）上で一様連続分布による近似（uniform ring近似）を設定することができる。なお、Ｍ−ＰＳＫは、Ｍ−ＡＰＳＫにおいて円（リング）の個数が１つの場合に相当している。
【０１３９】
図３（ａ）〜（ｃ）は、Ｍ−ＰＳＫにおいて、それぞれＢＰＳＫ、ＱＰＳＫ、および８ＰＳＫの例を示し、図３（ｅ）は単一の円（リング）上の一様連続分布の近似を示している。また、図５（ａ）はＭ−ＡＰＳＫにおいてＤＡＰＳＫの例を示し、図５（ｂ）は２つの円（リング）上の一様連続分布の近似を示している。
【０１４０】
このuniform ring近似は、円（リング）上に離散的に分布するＭ−ＡＰＳＫ、Ｍ−ＰＳＫの信号点配置を、信号点数が無限大の円（リング）の分布で近似するものである。
【０１４１】
このuniform ring近似を表す分布関数f(x)は、以下の式（５０）で表される。
2012066981.tif 000052
【０１４２】
式（５０）を用いると共に、式（１６）、（１７）の被積分変数ｘをｘ＝ρκｅｊθの極座標表示し、φ０を任意の実数定数とすることによって、以下の式（５１）、（５２）が得られる。
【０１４３】
2012066981.tif 000053
2012066981.tif 000054
【０１４４】
ここで、Ｉ０（・），Ｉ１（・）はそれぞれ０次と１次の第１種変形ベッセル関数である。これらのベッセル関数は、以下に示すように指数関数で高精度に近似できることが知られている。
【０１４５】
2012066981.tif 000055
【０１４６】
（第１のタイプ）
推定する送信ベクトルに対してuniform ring近似による近似を適用すると、式（５１），（５２）においてベッセル関数を指数関数で近似したα（ｘ２）、β（ｘ２）を式（１５）に代入することによって、推定送信シンボルＸ＾１は以下の式（５３）で表される。
【０１４７】
2012066981.tif 000056
【０１４８】
対数領域では、以下の式（５４）で表される。
2012066981.tif 000057
【０１４９】
ＰＳＫは、ＡＰＳＫにおいてＫ＝１，ρ１＝１とした単一リングに相当する。このＰＳＫでは、上記式（５４）は以下の式（５５）で表される。
【０１５０】
2012066981.tif 000058
【０１５１】
（第２のタイプ）
推定を行わない送信ベクトルに対してuniform ring近似を適用すると、式（１８）に式（５１）を代入することによって、推定送信シンボルＸ＾１は以下の式（５６）で表される。
【０１５２】
2012066981.tif 000059
【０１５３】
対数領域では、以下の式（５７）で表される。
2012066981.tif 000060
【０１５４】
ＰＳＫは、ＡＰＳＫにおいてＫ＝１，ρ１＝１とした単一リングに相当する。このＰＳＫでは、上記式（５７）は以下の式（５８）で表される。
【０１５５】
2012066981.tif 000061
【０１５６】
表１に導出した推定式を示す。
【表１】
2012066981.tif 000062
【０１５７】
［実施例］
以下に、本発明の信号分離による計算機シミュレーションの例を示す。
ここで、通信路行列Ｈのフェージング係数hｋ、ｎ,k=１,…,Ｎｒ、ｎ＝１，２は全て独立で、各々平均０、分散１の複素ガウス分布に従うと仮定する。図６〜図１１において、式（３０’）、式（３１’）を用いて得られる第１のタイプおよび第２のタイプのシミュレーション結果については、それぞれtype-Iおよびtype-IIで示し、式（３０’）、式（３１’）を用いて対数演算の項を省略したmax演算のみで得られる第１のタイプおよび第２のタイプのシミュレーション結果については、それぞれtype-I(max）およびtype-II(max）で示している。
【０１５８】
図６〜図８は、受信アンテナ数がＮｒ＝２の場合において、uniform ring近似に基づいて本発明の信号分離法で受信したＱＰＳＫ，８−ＰＳＫ，１６−ＰＳＫのビット誤り率(ＢＥＲ)特性を示している。
【０１５９】
第１のタイプ（黒の三角印で示す）と第２のタイプ（白の三角印で示す）で示すビット誤り率を比較すると、第１のタイプは第２のタイプよりも優れた特性を有することが確認される。つまり、干渉となるシンボルを近似するよりも、推定を行う所望のシンボルを近似する方が好ましいことを示している。 
【０１６０】
ＢＥＲ＝１０−４において、第１のタイプと線形ＭＭＳＥ等化の特性差は１１−１２dBである。また、ＢＥＲを最小化できるＭＬ検出と比較すると、同じＢＥＲにおける差は３dB以内に収まっている。
【０１６１】
図９は、８−ＰＳＫで受信アンテナ数Ｎｒを増加させたときのＢＥＲ特性を示している。第２のタイプは第１のタイプよりも特性が劣ることを既に確認済みであるため、ここでは第１のタイプの結果のみを示している。図９から、受信アンテナ数Ｎｒを増加させるとＭＬＤとの差が縮まることが確認できる。 
【０１６２】
また、受信アンテナ数Ｎｒ＝３においてＭＬＤとの差はわずか０．８dBであり、受信アンテナ数Ｎｒ＝４の場合ではグラフ上で視認できる差は無い。
【０１６３】
図１０は、受信アンテナ数がＮｒ＝２の場合において、１６−ＡＰＳＫを用いたときのＢＥＲを示している。ここで用いた１６−ＡＰＳＫは振幅が異なる二つの８−ＰＳＫの重ね合わせであり、Ｋ＝２，ρ1＝√２/５，ρ２＝２√２／５，Ｍ（１）＝Ｍ（２）＝８である。これらのリング半径は、よく用いられるリング比１：２に対応している。図１０より 、ＰＳＫの場合と同様の傾向のＢＥＲ特性が確認できる。 
【０１６４】
ＡＰＳＫの推定式はＰＳＫの推定式と比べ複雑度が高いが、第１のタイプとＭＬＤの差はＰＳＫの場合より大きい。
【０１６５】
図１１は、受信アンテナ数がＮｒ＝２の場合において、１６−ＱＡＭをuniform square近似に基づいて受信したときのＢＥＲ特性を示している。 
【０１６６】
ＰＳＫおよびＡＰＳＫの場合と同様に、第１のタイプの方が第２のタイプよりも特性が優れているが、ＢＥＲ特性の傾きから、達成できるダイバーシティオーダがＰＳＫおよびＡＰＳＫの場合より低いことが確認できる。その結果、ＢＥＲ＝１０−４において線形ＭＭＳＥ等化と比較したときのＥｂ／Ｎ０の改善量は、第１のタイプでも３．５dB程度である。 
【０１６７】
図１２，１３はGauss近似の特性を示している。なお、Gauss近似に基づく推定式は、ＰＳＫおよびＱＡＭの両方に適用し、uniform ring近似およびuniform square近似の特性を比較している。 
【０１６８】
１６−ＱＡＭの結果より、ＱＡＭにおけるGauss近似の特性はuniform square近似の特性とほぼ同じであることが分かる。一方、ＰＳＫにおいては、uniform ring近似がGauss近似よりも大幅に優れていることが確認できる。この結果より、シンボル点分布の近似法が特性に大きな影響を与えると言える。 
【０１６９】
処理の複雑度の評価として、表２に計算機上に実装したプログラムの実行時間を示す。
【表２】
2012066981.tif 000063
【０１７０】
第１のタイプおよび第２のタイプの処理時間は、簡易な手法である線形ＭＭＳＥ法には及ばないものの、ＭＬＤと比べて大幅に短縮されていることが確認できる。 
【０１７１】
本発明の信号分離によれば、Bayesian型 ＭＭＳＥ等化において、何れか一つの送信アンテナから到来する送信シンボルの分布を一様連続分布で近似することによって、Bayesian型 ＭＭＳＥ等化の推定式に含まれる積分について閉形式とすることができるため、Ｍ２のオーダーをＭのオーダーの計算量に低減して信号分離を実行することができる。また、シミュレーション結果から、ＭＬ検出に近い特性が達成できることも確認できる。ＡＰＳＫ、ＰＳＫに対するuniform ring、ＱＡＭに対する uniform square近似、および Gauss近似の三種類の近似を、シミュレーション結果により比較した場合には、uniform ring近似が優れた特性を示すことが確認される。
【０１７２】
なお、本発明は前記各実施の形態に限定されるものではない。本発明の趣旨に基づいて種々変形することが可能であり、これらを本発明の範囲から排除するものではない。
【産業上の利用可能性】
【０１７３】
本発明は、Ｍ−ＡＰＳＫ、Ｍ−ＰＳＫ、Ｍ−ＱＡＭの各変調方式に適用することができる。
【符号の説明】
【０１７４】
１    送信アンテナ
２    伝送路
３    受信アンテナ
４    チャンネル推定部
５    検出器
５Ａ    演算部
５Ｂ    算出部
Ｘ    送信ベクトル
Ｙ    受信ベクトル
Ｈ    通信路行列
Ｗ    雑音ベクトル
【図１】
2012066981.tif 000064
【図２】
2012066981.tif 000065
【図５】
2012066981.tif 000066
【図６】
2012066981.tif 000067
【図７】
2012066981.tif 000068
【図８】
2012066981.tif 000069
【図９】
2012066981.tif 000070
【図１０】
2012066981.tif 000071
【図１１】
2012066981.tif 000072
【図１２】
2012066981.tif 000073
【図１３】
2012066981.tif 000074
【図１４】
2012066981.tif 000075
【図３】
2012066981.tif 000076
【図４】
2012066981.tif 000077
【国際調査報告】
2012066981.tif 000078
2012066981.tif 000079
2012066981.tif 000080