(57)【特許請求の範囲】
【請求項１】
複数の送信アンテナから、ＣＤＭＡ方式により送信された送信信号を、複数の受信アンテナにより受信する受信装置であって、
各受信アンテナで受信した受信信号を１次復調して各送信アンテナからの送信信号を推定し、該推定結果に基づいてマルチパス環境における受信アンテナ毎の各パスの受信信号を再生するマルチパス受信信号再生手段と、
前記各受信アンテナで受信した受信信号から、着目するパス以外のパスに対応する前記マルチパス受信信号再生手段により再生された再生受信信号を減じるマルチパス干渉キャンセル手段と、
該マルチパス干渉キャンセル手段により、前記再生受信信号が減じられた受信信号を用いて２次復調を行う復調手段と
を備えることを特徴とする受信装置。
【請求項２】
請求項１に記載の受信装置であって、
前記マルチパス受信信号再生手段は、最小平均自乗誤差法（MMSE: Minimum Mean Square Error）を用いて、各受信アンテナで受信した受信信号を１次復調することを特徴とする受信装置。
【請求項３】
請求項１に記載の受信装置であって、
前記マルチパス受信信号再生手段は、最尤検出法（MLD: Maximum Likelihood Detention）を用いて、各受信アンテナで受信した受信信号を１次復調することを特徴とする受信装置。
【請求項４】
請求項１に記載の受信装置であって、
前記マルチパス受信信号再生手段は、ＱＲ分解を利用した最尤検出法を用いて複数パスを一括して前記１次復調を実行することを特徴とする受信装置。
【請求項５】
請求項１に記載の受信装置であって、
前記マルチパス受信信号再生手段は、ＱＲ分解を利用した最尤検出法を用いてパス毎に前記１次復調を実行することを特徴とする受信装置。
【請求項６】
請求項２に記載の受信装置であって、
前記マルチパス受信信号再生手段は、各受信アンテナで受信した受信信号を１次復調することにより推定された送信シンボル系列の確からしさに基づいて、受信信号の振幅を制御することを特徴とする受信装置。
【請求項７】
請求項２に記載の受信装置であって、
前記マルチパス受信信号再生手段は、前記複数の送信アンテナから送信される既知のパイロット信号を用いてチャネル係数を推定することを特徴とする受信装置。
【請求項８】
請求項２に記載の受信装置であって、
所定数の前記マルチパス受信信号再生手段及び前記マルチパス干渉キャンセル手段を多段接続することを特徴とする受信装置。
【請求項９】
請求項８に記載の受信装置であって、
前記マルチパス受信信号再生手段が多段接続される場合に、２段目以降の各段において、前段のマルチパス干渉キャンセル手段により前記再生受信信号が減じられた受信信号を用いて、前記複数の送信アンテナから送信される既知のパイロット信号に基づいて推定されるチャネル係数推定値の更新を行うことを特徴とする受信装置。
【請求項１０】
請求項１に記載の受信装置であって、
前記復調手段は、最尤検出法を用いて２次復調を行うことを特徴とする受信装置。
【請求項１１】
請求項１に記載の受信装置であって、
前記復調手段は、ＱＲ分解を利用した最尤検出法を用いて複数パスを一括して２次復調を行うことを特徴とする受信装置。
【請求項１２】
請求項１に記載の受信装置であって、
前記復調手段は、ＱＲ分解を利用した最尤検出法を用いてパス毎に２次復調を行うことを特徴とする受信装置。
【請求項１３】
請求項１に記載の受信装置であって、
前記複数の送信アンテナから符号多重された送信信号が送信されたときに、
前記マルチパス受信信号再生手段は、各受信アンテナで受信した受信信号を１次復調して、拡散符号毎に受信アンテナ毎の各パスの受信信号を再生し、
前記マルチパス干渉キャンセル手段は、前記各受信アンテナで受信した受信信号から、着目するパス以外のパスの全拡散符号に対応する前記マルチパス受信信号再生手段により再生された再生受信信号を減じた信号を生成し、
前記復調手段は、前記マルチパス干渉キャンセル手段により前記再生受信信号が減じられた受信信号を用いて、拡散符号毎に２次復調を行うことを特徴とする受信装置。
【請求項１４】
複数の送信アンテナから、ＣＤＭＡ方式により送信された送信信号を、複数の受信アンテナにより受信する受信装置の受信方法であって、
各受信アンテナで受信した受信信号から、所定のアルゴリズムを用いて、各送信アンテナからの送信信号を推定する送信信号推定ステップと、
該送信信号推定ステップにより推定された送信信号と、既知のパイロット信号に基づいて推定されるチャネル推定値とを乗算することにより、マルチパス環境における受信アンテナ毎の各パスの受信信号を再生する受信信号再生ステップと、
前記各受信アンテナで受信した受信信号から、着目するパス以外のパスに対応する前記受信信号再生ステップにより再生された再生受信信号を減算する再生受信信号減算ステップと、
該再生受信信号減算ステップにより前記再生受信信号が減算じられた受信信号を用いて復調を行う復調ステップと
を有することを特徴とする受信方法。
【請求項１５】
複数の送信アンテナを備え、各送信アンテナからＣＤＭＡ信号を送信する送信装置と、請求項１に記載の受信装置と
を備えることを特徴とする無線通信システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数アンテナを用いて信号の復調（信号分離）を行う受信装置、受信方法、および無線通信システムに関する。
【背景技術】
【０００２】
第４世代移動体通信の無線通信方式においては、高速の伝送速度を実現することが重要である。かかる観点から、複数の送受信アンテナを用いて信号伝送を行うＭＩＭＯ（Multi-Input Multi-Output、複数入出力）チャネルを用いた信号伝送法において、各送信アンテナより異なる信号を同一時刻、同一周波数を用いて送信するＭＩＭＯ多重法（MIMO Multiplexing）の技術が注目されている。【０００３】
図１４は、ＭＩＭＯ多重法を説明するための図であり、複数のアンテナを用いたＭＩＭＯ通信システムの構成を示す図である。このシステムによれば、送信側の複数のアンテナ１０１１１〜１０１１Ｎから各々異なった信号を同じ周波数を用いて送り、受信側においても複数のアンテナ１０２１１〜１０２１Ｎを用いてそれら全てを同時に受信することにより、伝送帯域を増やさずに、送信アンテナ数に比例した伝送速度の高速化が可能である（例えば、非特許文献１参照）。【０００４】
また、ＭＩＭＯと総称される技術の１つにＢＬＡＳＴ（Bell Labs Layered Space-Time）がある。ＢＬＡＳＴは、異なる情報を複数の送信アンテナから同一周波数において同時に並列伝送し、受信側において干渉抑圧で制御されたダイバーシチ受信とレプリカ減算によって信号分離を行う方式である（例えば、非特許文献１参照）。【０００５】
上記の通り、ＭＩＭＯ多重法では高速の伝送速度を実現することが可能であるが、複数送信アンテナから異なるデータ系列が同一の周波数帯域、時間スロットで送信されるため、受信装置では受信信号から各送信アンテナから送信された送信信号系列を抽出する信号分離が復調のために必要となる。【０００６】
ＭＩＭＯ多重における信号分離法については種々の方法が提案されている。例えば、最小平均自乗誤差法（ＭＭＳＥ）やＺＦ（Zero Forcing）等の線形フィルタを用いた信号分離アルゴリズムは、送信アンテナ数と同じかそれよりも多い数の受信アンテナで受信した複数の受信信号を、着目する送信アンテナ以外の送信アンテナからの受信信号電力を抑圧するように合成する方法であり、受信側の演算量は比較的少なくてすむといった特徴がある。【０００７】
また、最尤検出法（ＭＬＤ）を用いた信号分離アルゴリズムは、各送信アンテナからの受信信号のレプリカ候補を生成し、受信信号と全送信アンテナからの受信信号のレプリカ候補の和とのユークリッド距離が最小になる受信信号レプリカを求めることで、最も確からしい各送信アンテナの送信信号系列を推定する方法であり、上記したＭＭＳＥに比較して、信号分離精度が高く、復調性能は優れるが、信号分離に要する演算量はアンテナ数に比例して指数的に増加するといった欠点がある。そこで、ＭＬＤにおける演算量を低減するために、ＱＲ分解を利用してＭＬＤにおける２乗ユークリッド距離を計算する信号点候補を大幅に削減する方法が提案されている（例えば、非特許文献２参照）。【０００８】
ところで、直接拡散（ＤＳ）ＣＤＭＡは従来の情報データ変調信号を高速レートの拡散符号で拡散する２次変調を行って伝送することで複数の通信者が同一の周波数帯を用いて通信を行う方式であり、ＤＳ−ＣＤＭＡを用いた無線通信では、信号伝送帯域幅の広帯域化により、マルチパスフェージング（周波数選択性フェージング）が生じ、送信信号は互いに伝搬遅延時間の異なる複数のマルチパスに分離されて受信される。【０００９】
ＤＳ−ＣＤＭＡでの受信では、この複数のパスを合成するＲａｋｅ受信により受信品質を改善することができるが、異なるパス間には干渉（以下、マルチパス干渉という）が生じるため、Ｒａｋｅ受信による受信品質改善がオフセットされる。【００１０】
上記マルチパス干渉は、拡散符号を乗算するスピードであるチップレートと情報シンボルのシンボルレートの比で定義される拡散率の逆数に大きさが比例する。このため、情報ビットレートを増大するために拡散率を１に近づけると、Ｒａｋｅ受信効果よりもマルチパス干渉による受信品質の劣化のほうが支配的になり、高速データ伝送時に受信特性が劣化するという問題があった。そこで、このような問題を解決するために、マルチパス干渉キャンセラが提案されている（例えば、非特許文献３参照）。【００１１】
このマルチパス干渉キャンセラは、非特許文献３に記載されているように、仮のＲａｋｅ受信結果で推定された送信信号系列と各受信パスのチャネル係数（伝搬路の複素包絡線）を基にパス毎の受信信号系列を推定し、受信信号からあるパス以外の全ての推定した受信信号系列を差し引くことを、パス数分繰り返すことで得られるマルチパス干渉を低減したパス毎の受信信号を用いて最終的なＲａｋｅ受信を行うことにより、マルチパス環境下での高品質受信を実現している。【００１２】
また、ＤＳ−ＣＤＭＡを用いた無線通信におけるＭＩＭＯ多重でのマルチパス干渉の影響を低減する信号分離法として、他の送信アンテナからの受信信号により生じる干渉の抑圧とマルチパス干渉の抑圧を同時に行う２次元ＭＭＳＥも提案されている。
【非特許文献１】G. J. Foschini, Layered Space-Time Architecture for Wireless Communication in a Fading Environment When Using Multiple Antennas, Bell Labs Technical Journal, Vol. 1, No. 2, autumn 1996, pp 41-59.
【非特許文献２】Bin Dong, Xaodong Wang, and Amaud Doucet, "Sampling-based Near-optimal MIMO demodulation Algorithms," in Proc. 42nd IEEE Conference on Decision and Control, Hawaii, Dec.2003.
【非特許文献３】K. Higuchi, A. Fujiwara, and M. Sawahashi, "Multipath Interference Canceller for High-Speed Packet Transmission With Adaptive Modulation and Coding Scheme in W-CDMA Forward Link, " IEEE J.Select. Areas Commun., Vol.20, No.2, pp.419-432, February 2002.
【非特許文献４】Frederik Petre et. Al, " Combined Space-Time Chip Equalization And Parallel Interference Cancellation For DS-CDMA Downlink With Spatial Multiplexing," in Proc. IEEE PIMRC2002.
【発明の開示】
【発明が解決しようとする課題】
【００１３】
ここで、ＤＳ−ＣＤＭＡを用いた無線通信に情報ビットレートの高速化のために、上記したＭＩＭＯ多重法を適用することを考えると、まず１送信アンテナあたりの情報ビットレートを増大するために拡散率を小さくしなければならない。この場合、図１５に示すように、同一の受信タイミングのパス同士の送信アンテナ間の相互干渉（同図（ａ）参照）に加えて、全送信アンテナの異なる受信タイミングのマルチパスからのマルチパス干渉（同図（ｂ）、（ｃ）参照）が生じる。このため、受信側での信号分離の精度は、先に述べた１アンテナ送信時のＲａｋｅ受信と同様に、マルチパス干渉の影響で大きく劣化すると考えられる。【００１４】
先に提案されたマルチパス干渉キャンセラでは、１アンテナからのみ信号が送信される場合の構成であり、ＭＩＭＯ多重時に適用しても、Ｒａｋｅ受信では高精度な送信信号系列の推定ができない。【００１５】
また、図１５から明らかなように、ＭＩＭＯ多重時におけるマルチパス干渉は同一送信アンテナからの干渉だけでなく、異なる送信アンテナからの受信信号からも生じるため、高精度な受信のためには、これらのマルチパス干渉も抑圧する必要がある。【００１６】
しかしながら、上記したＭＭＳＥは原理的にＭＬＤに比較して信号分離精度が悪く、受信側で最適な性能が得られない。また、２次元ＭＭＳＥにおいては、他の送信アンテナからの受信信号により生じる干渉の抑圧に加えて、マルチパス干渉も抑圧するようにフィルタ係数が制御されるため、信号分離精度の劣化がさらに大きくなるといった問題があった。【００１７】
本発明は、上記のような問題点に鑑みてなされたもので、その課題とするところは、マルチパス干渉による受信特性の劣化を抑圧して高精度な信号分離を実現することのできる受信装置、受信方法、および無線通信システムを提供することである。
【課題を解決するための手段】
【００１８】
本発明は、ＣＤＭＡ方式を用いて信号の受信を行う受信装置において、Ｍ本（Ｍは、正の整数）の送信アンテナから送信された送信信号を、Ｎ本（Ｎは、正の整数）の受信アンテナにより受信しているであって、各受信アンテナで受信した受信信号を１次復調して各送信アンテナからの送信信号を推定し、推定結果に基づいてマルチパス環境における受信アンテナ毎の各パスの受信信号を再生するマルチパス受信信号再生手段と、前記各受信アンテナで受信した受信信号から着目するパス以外のパスの前記再生受信信号を減じるマルチパス干渉キャンセル手段と、前記減じた信号を用いて２次復調を行う復調手段と、を備えることを特徴の１つとしている。【００１９】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、前記１次復調を、最小平均自乗誤差法（MMSE: Minimum Mean Square Error）を用いて実行することを特徴としている。【００２０】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、前記１次復調を、最尤検出法（MLD: Maximum Likelihood Detention）を用いて実行することを特徴としている。【００２１】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、ＱＲ分解を利用した最尤検出法を用いて複数パスを一括して前記１次復調を実行することを特徴としている。【００２２】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、ＱＲ分解を利用した最尤検出法を用いてパス毎に前記１次復調を実行することを特徴としている。【００２３】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、上記の方法を用いて推定される送信シンボル系列の確からしさに基づいて、受信信号の振幅を制御することを特徴としている。【００２４】
また、前記受信装置であって、前記マルチパス受信信号再生手段は、前記Ｍ本の送信アンテナから送信される既知のパイロット信号を用いてチャネル係数を推定することを特徴としている。【００２５】
また、前記受信装置であって、所定数の前記マルチパス受信信号再生手段及び前記マルチパス干渉キャンセル手段を多段接続することを特徴としている。【００２６】
また、前記受信装置であって、前記マルチパス受信信号再生手段が多段接続される場合に、各段において、前記マルチパス干渉キャンセル手段により減じられた信号を用いて、前記Ｍ本の送信アンテナから送信される既知のパイロット信号に基づいて推定されるチャネル係数推定値の更新を行うことを特徴としている。【００２７】
また、前記受信装置であって、前記復調手段は、最尤検出法を用いて２次復調を行うことを特徴としている。【００２８】
また、前記受信装置であって、前記復調手段は、ＱＲ分解を利用した最尤検出法を用いて複数パスを一括して２次復調を行うことを特徴としている。【００２９】
また、前記受信装置であって、前記復調手段は、ＱＲ分解を利用した最尤検出法を用いてパス毎に２次復調を行うことを特徴としている。【００３０】
また、前記受信装置であって、前記Ｍ本の送信アンテナから符号多重された送信信号が送信されたときに、前記マルチパス受信信号再生手段は、各受信アンテナで受信した受信信号を１次復調して拡散符号毎に受信アンテナ毎の各パスの受信信号を再生し、前記マルチパス干渉キャンセル手段は、前記各受信アンテナで受信した受信信号から着目するパス以外のパスの全ての拡散符号に対応する前記再生受信信号を減じた信号を生成し、前記復調手段は、前記減じた信号を用いて拡散符号毎に２次復調を行うことを特徴としている。
【発明の効果】
【００３１】
本発明の実施例によれば、ＣＤＭＡ方式を用いて複数の送信アンテナから異なるデータを同時送信したときに、マルチパス干渉を低減し、異なる送信アンテナから送信された信号の高精度な分離を実現することができる。その結果、マルチパスフェージング環境における、受信品質を大幅に向上させることができる。
【図面の簡単な説明】
【００３２】
【図１】本発明の実施の形態に係る受信装置を含んで構成される無線通信システム１を示す図である。
【図２】マルチパス受信信号再生部の第１の実施形態を示す構成図である。
【図３】マルチパス受信信号再生部の第２の実施形態を示す構成図である。
【図４】マルチパス受信信号再生部の第３の実施形態を示す構成図である。
【図５】マルチパス受信信号再生部の第４の実施形態を示す構成図である。
【図６】チャネル係数推定部の構成を示す図である。
【図７Ａ】送信装置より送信される送信信号のフレーム構成を示す図である。
【図７Ｂ】送信装置より送信される送信信号のフレーム構成を示す図である。
【図８】本発明の実施形態に係る受信装置の第２の実施形態を示す構成図である。
【図９】本発明の実施形態に係る復調部の第１の実施形態を示す構成図である。
【図１０】本発明の実施形態に係る復調部の第２の実施形態を示す構成図である。
【図１１】本発明の実施形態に係る復調部の第３の実施形態を示す構成図である。
【図１２】本発明の実施形態に係る受信装置の第３の実施形態を示す構成図である。
【図１３】本発明についての計算機シミュレーションの結果を示す図である。
【図１４】複数のアンテナを用いたＭＩＭＯ通信システムの構成を示す図である。
【図１５】ＤＳ−ＣＤＭＡにおけるＭＩＭＯ多重では、マルチパス干渉の影響で復調（信号分離）の精度が劣化することを説明するための図である。
【符号の説明】
【００３３】
１    無線通信システム
１０，１０１０    送信装置
１１１〜１１ｍ，１０１１１〜１０１１ｎ    送信アンテナ
２０，２００，６００，１０２０    受信装置
２１１〜２１ｎ，１０２１１〜１０２１ｎ    受信アンテナ
２２，３０〜５０，７０，２１１〜２１３，６１１〜６１４    マルチパス受信信号再生部
２３，２２１１，２２１２，２２２１，２２２２，２２３１，２２３２，６１５〜６１８    マルチパス干渉キャンセル部
２４，３００，４００，５００，６２１，６２２    復調部
３１，４１，５１，７１，１００，３１１，４１１，５１１    チャネル係数推定部
３２    線形フィルタ係数計算部
３３    線形フィルタリング部
３４１〜３４ｍ，４２〜４５，５２〜５５，７２〜７５，３１２〜３１５，４１２〜４１５，５１２〜５１５    逆拡散部
３５１〜３５ｍ    送信シンボル系列推定部
３６１〜３６ｍ，５０１〜５０ｍ，６２１〜６２ｍ，８４１〜８４ｍ    マルチパス受信信号再生処理部
４６，５８，８０，３１６，４１８，５２０    送信シンボル候補生成部
４７，３１７    受信信号レプリカ生成部
４８，６０，８２，３１８，４２０，５２２    尤度計算部
４９１〜４９ｍ，６１１〜６１ｍ，８３１〜８３ｍ    送信シンボル系列推定部
５６，７６，７７，４１６，５１６，５１７    ＱＲ分解部
５７，７８，７９，４１７，５１８，５１９    ＱＨ演算部
５９，８１，４１９，５２１    変換信号レプリカ生成部
１０１〜１０４    相関検出部
１１１〜１１４    レプリカ信号生成部
３１９，４２１，５２３    送信系列推定部
６３１〜６３４，６４１〜６４４    加算器
【発明を実施するための最良の形態】
【００３４】
以下、本発明の実施の形態を図面に基づいて説明する。【００３５】
図１は、本発明の実施の形態に係る受信装置を含んで構成される無線通信システム１を示す図である。本発明に係る無線通信システム１は、ＤＳ−ＣＤＭＡにおけるＭＩＭＯチャネルまたはＭＩＳＯ（Multi-Input Single-Output、複数入力単出力、つまり、受信装置の受信アンテナが１つ）チャネル用いたシステムであり、本実施形態では、ＭＩＭＯチャネルを用いて信号伝送がなされる場合を例にとり、以下説明する。【００３６】
同図において、この無線通信システム１は、送信装置１０と受信装置２０が無線通信において接続可能となっている。送信装置１０は、入力される送信データビット系列をＭ本の送信系列にシリアル・パラレル変換し、同一の周波数帯かつ同一の拡散符号を用いてデータ変調を行い、送信信号としてＭ本の送信アンテナ１１１〜１１ｍから同時に送信する。このようにして送信された送信信号はマルチパス伝搬路を経てＬ個のマルチパス受信信号となって受信装置２０に具備されるＮ本の受信アンテナ２１１〜２１ｎで受信される。受信装置２０は、Ｎ本の受信アンテナ２１１〜２１ｎと、マルチパス受信信号再生部２２と、Ｎ個のマルチパス干渉キャンセル部２３１〜２３ｎと、復調部２４とを具備して構成される。【００３７】
ここで、送信装置１０の送信アンテナｍから送信された送信信号をＳｍ（ｔ）とすると、受信装置２０の受信アンテナｎで受信される受信信号ｒｎ（ｔ）は、次式で表すことができる。【００３８】
【数１】
0004322918.tif 000002
上記式において、ｈｍ，ｎ，ｌ（ｔ）は、送信アンテナｍからの送信信号を受信アンテナｎで受信したときの受信パスｌに対するチャネル係数を表し、τｌは受信パスｌの伝搬遅延時間を表す。【００３９】
次に、本発明に係る受信装置２０の動作を説明する。【００４０】
受信装置２０では、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）が、マルチパス受信信号再生部２２に入力される。マルチパス受信信号再生部２２では、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号を基に、仮の信号分離（１次復調という）を行うことで、受信アンテナ２１１〜２１ｎごとに各送信アンテナ１１１〜１１ｍからの受信パス毎の受信信号系列【００４１】
【数２】
0004322918.tif 000003
を推定して出力する。上記１次復調の方法には、後述する所定のアルゴリズムが用いられる。【００４２】
次に、受信アンテナ２１１〜２１ｎ数分の各マルチパス干渉キャンセル部２３１〜２３ｎでは、当該受信アンテナ２１１〜２１ｎの受信信号と、各送信アンテナ１１１〜１１ｍからの受信パス毎の受信信号系列を入力とし、受信信号から他のパスの受信信号全てを差し引いたマルチパス干渉キャンセル後の受信信号ｒｎ，ｌ（ｔ）を次式にしたがって演算し、出力する。【００４３】
【数３】
0004322918.tif 000004
復調部２４は、上記のようにしてマルチパス干渉キャンセル部２３１〜２３ｋから出力されるＮ×Ｌ個のマルチパス干渉キャンセル後の受信信号ｒｎ，ｌ（ｔ）を入力とし、送信シンボル系列の尤度もしくは、送信シンボル系列によって送信されたビットに対する尤度λｉを２次復調出力として後続のビタビ復号器、ターボ復号器等の誤り訂正（チャネル）復号器に出力して復号処理を実行させる。【００４４】
このように本実施形態によれば、マルチパス受信信号再生部で、１次復調が行われ、送信信号が推定される。そして、その推定された送信信号とチャネル変動値（チャネル係数）とを掛け合わせることでパス毎の受信信号を推定し、マルチパス干渉キャンセル部で、受信信号から着目するパス以外の推定した受信信号が減算される。これにより、復調部では、マルチパス干渉を除去した後の受信信号を用いて復調することが可能となり、高精度に信号分離を行うことができる。【００４５】
すなわち、上りリンクにＤＳ−ＣＤＭＡ方式を用い、ＭＩＭＯ多重法を適用した場合であっても、マルチパス干渉に起因した信号分離精度の劣化を回避することができる。【００４６】
図２は、図１に示すマルチパス受信信号再生部の第１の実施形態を示す構成図である。本実施形態では、マルチパス受信信号再生部は、１次復調方法としてＭＭＳＥアルゴリズムを用いる。【００４７】
同図において、このマルチパス受信信号再生部３０は、チャネル係数推定部３１と、線形フィルタ係数計算部３２と、線形フィルタリング部３３と、Ｍ個の逆拡散部３４１〜３４ｍと、Ｍ個の送信シンボル系列推定部３５１〜３５ｍと、Ｍ個のマルチパス受信信号再生処理部３６１〜３６ｍとから構成される。【００４８】
本実施形態におけるマルチパス受信信号再生部３０では、最初に、チャネル係数推定部３１において、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を入力し、受信アンテナ２１１〜２１ｎと送信アンテナ１１１〜１１ｎ（図１参照）との間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される。【００４９】
次に、線形フィルタ係数計算部３２において線形フィルタの係数計算、及びその求めた線形フィルタ係数を用いた等化のためのフィルタリング処理が行われる。これらの処理は、時間領域の信号処理で実現する方法と、周波数領域の信号処理で実現する方法とが考えられるが、以下では、周波数領域で信号処理を行う方法について例を挙げて説明する。【００５０】
線形フィルタ係数計算部３２では、得られたチャネル係数の推定値【００５１】
【数４】
0004322918.tif 000005
を用いて、送信アンテナ毎に該送信アンテナの信号の遅延パス成分と他の送信アンテナからの受信信号成分を同時に抑圧する線形フィルタの係数を計算する。具体的には、時間領域でパイロットチャネルを用いて推定したチャネル推定値【００５２】
【数５】
0004322918.tif 000006
と、各受信パスの遅延時間【００５３】
【数６】
0004322918.tif 000007
を用いて送信アンテナ１１１〜１１mと受信アンテナ２１１〜２１ｎ間のチャネルのインパルス応答を求める。【００５４】
次に、上記のようにして求められたチャネルのインパルス応答に対して、線形フィルタリングするブロックサイズに相当するチップ数×オーバサンプリング数分の大きさとなるＮｆポイントのＦＦＴを行うことで、送信アンテナ１１１〜１１mからの送信信号の受信アンテナ２１１〜２１ｎでのチャネル変動値の周波数成分である【００５５】
【数７】
0004322918.tif 000008
を推定し、これより、ｎ行ｍ列の行列表記したチャネル行列【００５６】
【数８】
0004322918.tif 000009
の推定値【００５７】
【数９】
0004322918.tif 000010
を算出する。そして、このようにして算出された【００５８】
【数１０】
0004322918.tif 000011
を用いて、ＦＦＴ後の周波数成分毎に、線形フィルタリング係数を計算する。【００５９】
線形フィルタの係数は、ＺＦ基準あるいはＭＭＳＥ基準にしたがって求めることができる。【００６０】
例えば、ＺＦ基準の線形フィルタの係数は、次式により求めることができる.【００６１】
【数１１】
0004322918.tif 000012
また、ＭＭＳＥ基準の線形フィルタの係数は、次式により求めることができる。【００６２】
【数１２】
0004322918.tif 000013
ここで、【００６３】
【数１３】
0004322918.tif 000014
は雑音成分を表す。【００６４】
次に、線形フィルタリング部３３では、ＮｆポイントのＦＦＴにより、Ｎ系列の受信信号それぞれを周波数領域の受信信号Ｙ（ｆ）に変換する。その後、周波数領域の受信信号Ｙ（ｆ）にＷ（ｆ）を乗算することにより、周波数領域のチャネル変動によるコードチャネル間の直交性の崩れ（時間領域で見たＭＰＩI(Multi-Path-Interference)）と送信アンテナ間の干渉を同時に等化（抑圧）したＭ個の送信信号の推定値【００６５】
【数１４】
0004322918.tif 000015
を次式にしたがって生成する。【００６６】
【数１５】
0004322918.tif 000016
なお、【００６７】
【数１６】
0004322918.tif 000017
はＮｆポイントのＩＦＦＴ後に並直列変換することで時間領域の仮の復調信号（１次復調信号）【００６８】
【数１７】
0004322918.tif 000018
毎に再変換される。このようにして得られた送信アンテナ数分の仮の復調系列（１次復調系列）は、逆拡散部３４１〜３４ｍに入力される。逆拡散部３４１〜３４ｍでは、送信時に用いた拡散符号と同一の拡散符号で上記送信アンテナ数分の仮の復調系列を逆拡散し、逆拡散信号ｚｍを送信シンボル系列推定部３５１〜３５ｍに出力する。送信シンボル系列推定部３５１〜３５ｍでは、逆拡散信号ｚｍを硬判定もしくは軟判定することにより送信シンボル系列【００６９】
【数１８】
0004322918.tif 000019
を推定し出力する。【００７０】
例えば、送信シンボル系列推定部３５１〜３５ｍにおいて硬判定を行う場合の実施例は以下のようになる。【００７１】
ここで、送信シンボルｄｍ，ｉ（ｉはシンボル候補番号【００７２】
【数１９】
0004322918.tif 000020
を表す。なお、Ｃは送信シンボル点数を表し、ＱＰＳＫならＣ＝４、１６ＱＡＭならＣ＝１６である）に対して、【００７３】
【数２０】
0004322918.tif 000021
となるｉより、【００７４】
【数２１】
0004322918.tif 000022
とする。【００７５】
また、送信シンボル系列推定部３５１〜３５ｍにおいて軟判定を行う場合の実施例は以下のようになる。【００７６】
まず、次式にしたがって、軟判定のビット系列を求める。【００７７】
【数２２】
0004322918.tif 000023
ここで、【００７８】
【数２３】
0004322918.tif 000024
は、送信アンテナｍのビットｉの対数尤度比である。【００７９】
【数２４】
0004322918.tif 000025
上記式において、Ｓｍｉｎ，νは、第ｉビットが"ν" であるシンボルの集合のなかで、送信信号点ｚｍからのユークリッド距離が最小のシンボル候補を表し、σ２
は、雑音電力を表す。【００８０】
軟判定シンボル【００８１】
【数２５】
0004322918.tif 000026
は、【００８２】
【数２６】
0004322918.tif 000027
を用いて以下のように推定される。【００８３】
【数２７】
0004322918.tif 000028
最後に、マルチパス受信信号再生処理部３６１〜３６ｍは、上記のようにして推定された送信シンボル系列【００８４】
【数２８】
0004322918.tif 000029
を入力し、次式に基づき、その入力された送信シンボル系列【００８５】
【数２９】
0004322918.tif 000030
に拡散符号とチャネル係数を乗算することにより、受信アンテナごとに各送信アンテナからの受信パス毎の受信信号系列【００８６】
【数３０】
0004322918.tif 000031
を次式にしたがって推定して出力する。【００８７】
【数３１】
0004322918.tif 000032
次に、上記マルチパス受信信号再生部の他の実施形態について説明する。【００８８】
図３は、マルチパス受信信号再生部の第２の実施形態を示す構成図である。本実施形態では、マルチパス受信信号再生部は、１次復調方法としてＭＬＤアルゴリズムを用いる。【００８９】
同図において、このマルチパス受信信号再生部４０は、チャネル係数推定部４１と、Ｎ×Ｌ個の逆拡散部４２〜４５と、送信シンボル候補生成部４６と、受信信号レプリカ生成部４７と、尤度計算部４８と、Ｍ個の送信シンボル系列推定部４９１〜４９ｍと、Ｍ個のマルチパス受信信号再生処理部５０１〜５０ｍとから構成される。【００９０】
本実施形態におけるマルチパス受信信号再生部４０では、最初に、チャネル係数推定部４１において、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を入力し、受信アンテナ２１１〜２１ｎと送信アンテナ１１１〜１１ｎとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される。【００９１】
次に、Ｎ×L個の逆拡散部４２〜４５により、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を送信時に用いた拡散符号と同一の拡散符号で各パスの受信タイミングで逆拡散し、Ｎ×Ｌ個の各受信アンテナでの各パスの逆拡散信号ｚｎ，ｌを求める。【００９２】
送信シンボル候補生成部４６は、各送信アンテナに対する送信シンボルｄｍ，ｉ（ｉはシンボル候補番号【００９３】
【数３２】
0004322918.tif 000033
を表し、Ｃは、送信シンボル点数を表し、例えば、ＱＰＳＫならＣ＝４、１６ＱＡＭならＣ＝１６である）を生成して出力する。【００９４】
受信信号レプリカ生成部４７は、送信シンボル候補生成部４６で生成された送信シンボルと、チャネル係数推定部４１で推定されたチャネル係数とを入力して受信信号レプリカ【００９５】
【数３３】
0004322918.tif 000034
を次式にしたがって生成して出力する。【００９６】
【数３４】
0004322918.tif 000035
尤度計算部４８は、受信逆拡散信号ｚｎ，ｌと受信信号レプリカ【００９７】
【数３５】
0004322918.tif 000036
を入力し、次式にしたがって誤差演算を行う。【００９８】
【数３６】
0004322918.tif 000037
送信シンボル系列推定部４９１〜４９ｍでは、各送信アンテナに対応する生成した送信シンボルとそれに対応する誤差信号とを入力し、最小誤差を選択し、当該誤差を与える送信シンボル系列【００９９】
【数３７】
0004322918.tif 000038
を推定する。この送信シンボル系列推定部４９１〜４９ｍにおいて硬判定を行う場合の実施例は以下のようになる。【０１００】
【数３８】
0004322918.tif 000039
また、送信シンボル系列推定部４９１〜４９ｍにおいて軟判定を行う場合の実施例は以下のようになる。【０１０１】
まず、次式にしたがって、軟判定のビット系列を求める。【０１０２】
【数３９】
0004322918.tif 000040
ここで、【０１０３】
【数４０】
0004322918.tif 000041
は、送信アンテナｍのビットｉの対数尤度比である。【０１０４】
【数４１】
0004322918.tif 000042
上記式において、ｅｍｉｎ，νは、第ｉビットが"ν" である【０１０５】
【数４２】
0004322918.tif 000043
の最小値、σ２（ｎ）は雑音電力をあらわす。【０１０６】
軟判定シンボル【０１０７】
【数４３】
0004322918.tif 000044
は、【０１０８】
【数４４】
0004322918.tif 000045
を用いて以下のように推定される。【０１０９】
【数４５】
0004322918.tif 000046
最後に、マルチパス受信信号再生部３６１〜３６ｍは、上記のようにして推定された送信シンボル系列【０１１０】
【数４６】
0004322918.tif 000047
を入力し、次式に基づき、その入力された送信シンボル系列【０１１１】
【数４７】
0004322918.tif 000048
に拡散符号とチャネル係数を乗算することにより、
受信アンテナごとに各送信アンテナからの受信パス毎の受信信号系列【０１１２】
【数４８】
0004322918.tif 000049
を次式にしたがって推定して出力する。【０１１３】
【数４９】
0004322918.tif 000050
図４は、マルチパス受信信号再生部の第３の実施形態を示す構成図である。本実施形態では、マルチパス受信信号再生部は、１次復調方法としてパス一括処理を行うＭＬＤアルゴリズムを用いる。【０１１４】
同図において、このマルチパス受信信号再生部５０は、チャネル係数推定部５１と、Ｎ×Ｌ個の逆拡散部５２〜５５と、ＱＲ分解部５６と、ＱＨ演算部５７と、送信シンボル候補生成部５８と、変換信号レプリカ生成部５９と、尤度計算部６０と、Ｍ個の送信シンボル系列推定部６１１〜６１ｍと、Ｍ個のマルチパス受信信号再生処理部６２１〜６２ｍとから構成される。【０１１５】
本実施形態におけるマルチパス受信信号再生部５０では、最初に、チャネル係数推定部５１において、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を入力し、受信アンテナ２１１〜２１ｎと送信アンテナ１１１〜１１ｎとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される。【０１１６】
次に、Ｎ×L個の逆拡散部５２〜５５により、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を送信時に用いた拡散符号と同一の拡散符号で各パスの受信タイミングで逆拡散し、Ｎ×Ｌ個の各受信アンテナでの各パスの逆拡散信号ｚｎ，ｌを求める。【０１１７】
次に、ＱＲ分解部５６では、チャネル係数からなる下記のＭ行×（Ｎ×Ｌ）列のチャネル行列を生成し、チャネル行列のＱＲ分解を行って、Ｑ行列およびＲ行列を出力する。【０１１８】
【数５０】
0004322918.tif 000051
上記ＱＲ分解部５６から出力されるＱ行列は、（Ｎ×Ｌ）行×Ｍ列のユニタリー行列であり、ＱＨＱ＝Ｉを満たす。【０１１９】
ここで、Ｈは共役複素転置を表し、Ｉは単位行列を表す。また、Ｒ行列はＭ行×Ｍ列の上三角行列となる。【０１２０】
ＱＨ演算部部５７では、次式にしたがった演算が行われる。【０１２１】
【数５１】
0004322918.tif 000052
送信シンボル候補生成部５８は、まず、送信アンテナＭに対する送信シンボルｄＭ，ｉを生成して変換信号レプリカ生成部５９に出力する。変換信号レプリカ生成部５９は、送信シンボル候補生成部５８で生成された送信シンボルと、ＱＲ分解部５６から出力されるＲ行列を入力して、変換信号レプリカ【０１２２】
【数５２】
0004322918.tif 000053
を次式にしたがって生成して出力する。【０１２３】
【数５３】
0004322918.tif 000054
尤度計算部６０は、まず、ｘＭと【０１２４】
【数５４】
0004322918.tif 000055
を用い、次式にしたがって誤差演算を行う。【０１２５】
【数５５】
0004322918.tif 000056
尤度計算部６０は、誤差演算を行った後、送信アンテナＭに対応する生成した送信シンボルとそれに対応する誤差信号とを入力し、誤差の小さいＳＭ個の送信アンテナＭに対応する送信シンボル系列【０１２６】
【数５６】
0004322918.tif 000057
とそのときの誤差【０１２７】
【数５７】
0004322918.tif 000058
を保持する。【０１２８】
次に、送信シンボル候補生成部５８は、送信アンテナＭ−１に対する送信シンボル【０１２９】
【数５８】
0004322918.tif 000059
を生成して出力する。【０１３０】
変換信号レプリカ生成部５９は、送信アンテナＭに対応するＳＭ個の送信シンボル系列と送信アンテナＭ−１に対する送信シンボルとからＲ行列を入力として変換信号レプリカ【０１３１】
【数５９】
0004322918.tif 000060
を次式にしたがって生成して出カする。【０１３２】
【数６０】
0004322918.tif 000061
尤度計算部６０は、次いで、誤差演算を次式により行い、【０１３３】
【数６１】
0004322918.tif 000062
誤差の小さいＳＭ−１個の送信アンテナＭと送信アンテナＭ−１に対応する送信シンボル系列の組み合わせ【０１３４】
【数６２】
0004322918.tif 000063
とそのときの誤差【０１３５】
【数６３】
0004322918.tif 000064
を保持する。【０１３６】
同様にして、送信シンボル候補生成部５８は、送信アンテナｍに対する送信シンボルｄｍ，ｉを生成して出力する。変換信号レプリカ生成部５９は、送信アンテナｍ＋１から、送信アンテナＭに対応するＳｍ+１個の送信シンボル系列と送信アンテナｍに対する送信シンボルとＲ行列を入力として変換信号レプリカ【０１３７】
【数６４】
0004322918.tif 000065
を次式にしたがって生成して出カする。【０１３８】
【数６５】
0004322918.tif 000066
尤度計算部６０は、誤差演算を次式にしたがって行い、【０１３９】
【数６６】
0004322918.tif 000067
誤差の小さいＳｍ個の送信アンテナＭから送信アンテナｍに対応する送信シンボル系列の組み合わせ【０１４０】
【数６７】
0004322918.tif 000068
とそのときの誤差【０１４１】
【数６８】
0004322918.tif 000069
を保持する。【０１４２】
以上の操作を繰り返すことにより、得られたＣ・Ｓ２個の全送信アンテナに対応する送信シンボル系列の組み合わせ【０１４３】
【数６９】
0004322918.tif 000070
とそのときの誤差【０１４４】
【数７０】
0004322918.tif 000071
を得る（ｊ＝１〜Ｃ・Ｓ２の整数）。【０１４５】
送信シンボル系列推定部６１１〜６１ｍでは、各送信アンテナに対応する、生き残り送信シンボルとそれに対応する誤差信号とを入力とし、最小誤差を選択し、当該誤差を与える送信シンボル系列を推定する。【０１４６】
送信シンボル系列推定部６１１〜６１ｍにおいて硬判定を行う場合の実施例は以下のようになる。【０１４７】
【数７１】
0004322918.tif 000072
また、送信シンボル系列推定部６１１〜６１ｍにおいて軟判定を行う場合の実施例は以下のようになる。【０１４８】
まず、次式にしたがって、軟判定のビット系列を求める。【０１４９】
【数７２】
0004322918.tif 000073
ここで、【０１５０】
【数７３】
0004322918.tif 000074
は、送信アンテナｍのビットｉの対数尤度比である。【０１５１】
【数７４】
0004322918.tif 000075
上記式において、ｅｍｉｎ，νは、第ｉビットが"ν" である【０１５２】
【数７５】
0004322918.tif 000076
の最小値、σ２（ｎ）は雑音電力を表す。【０１５３】
軟判定シンボル【０１５４】
【数７６】
0004322918.tif 000077
は、【０１５５】
【数７７】
0004322918.tif 000078
を用いて以下のように推定される。【０１５６】
【数７８】
0004322918.tif 000079
最後に、マルチパス受信信号再生処理部６２１〜６２ｍは、上記のようにして推定された送信シンボル系列【０１５７】
【数７９】
0004322918.tif 000080
を入力し、次式に基づき、その入力された送信シンボル系列【０１５８】
【数８０】
0004322918.tif 000081
に拡散符号とチャネル係数を乗算することにより、受信アンテナごとに各送信アンテナからの受信パス毎の受信信号系列【０１５９】
【数８１】
0004322918.tif 000082
を次式にしたがって推定して出力する。【０１６０】
【数８２】
0004322918.tif 000083
以上説明したように、本実施形態におけるマルチパス受信信号再生部の構成は、図３に示すマルチパス受信信号再生部の構成に比較して、若干の受信信号系列の再生精度の劣化を許容すれば、誤差演算の回数をＣＭ回から【０１６１】
【数８３】
0004322918.tif 000084
回に大幅に低減することが可能であり、受信装置の複雑さを大幅に低減することが可能となる。【０１６２】
図５は、マルチパス受信信号再生部の第４の実施形態を示す構成図である。本実施形態では、マルチパス受信信号再生部は、１次復調方法としてパス毎に処理を行うＭＬＤアルゴリズムを用いる。【０１６３】
同図において、このマルチパス受信信号再生部７０は、チャネル係数推定部７１と、Ｎ×Ｌ個の逆拡散部７２〜７５と、Ｌ個のＱＲ分解部７６、７７と、Ｌ個のＱＨ演算部７８、７９と、送信シンボル候補生成部８０と、変換信号レプリカ生成部８１と、尤度計算部８２と、Ｍ個の送信シンボル系列推定部８３１〜８３ｍと、Ｍ個のマルチパス受信信号再生処理部８４１〜８４ｍとから構成される。【０１６４】
本実施形態におけるマルチパス受信信号再生部７０では、最初に、チャネル係数推定部７１において、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を入力し、受信アンテナ２１１〜２１ｎと送信アンテナ１１１〜１１ｎとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される。【０１６５】
次に、Ｎ×Ｌ個の逆拡散部７２〜７５により、Ｎ本の受信アンテナ２１１〜２１ｎで受信された受信信号ｒｎ（ｔ）を送信時に用いた拡散符号と同一の拡散符号で各パスの受信タイミングで逆拡散し、Ｎ×Ｌ個の各受信アンテナでの各パスの逆拡散信号ｚｎ，ｌを求める。【０１６６】
次に、第ｌ番目のＱＲ分解部では、第ｌ番目のパスのチャネル係数からなる下記のチャネル行列をパス数分生成し、それぞれ、チャネル行列のＱＲ分解を行って、Ｑ行列およびＲ行列を出力する。【０１６７】
【数８４】
0004322918.tif 000085
上記ＱＲ分解部から出力されるＱｌ行列は、Ｎ行×Ｍ列のユニタリー行列であり、ＱｌＨＱｌ＝Ｉを満たす。また、Ｒｌ行列は、Ｍ行×Ｍ列の上三角行列となる。【０１６８】
したがって、第ｌ番目のＱＨ演算部における演算は、【０１６９】
【数８５】
0004322918.tif 000086
と記述することができる。【０１７０】
送信シンボル候補生成部８０は、まず、送信アンテナＭに対する送信シンボルｄＭ，ｉを生成して変換信号レプリカ生成部８１に出力する。変換信号レプリカ生成部８１は、送信シンボルｄＭ，ｉと、Ｒｌ行列を入力して、変換信号レプリカ【０１７１】
【数８６】
0004322918.tif 000087
を次式にしたがって生成して出力する。【０１７２】
【数８７】
0004322918.tif 000088
尤度計算部８２は、まず、ｘＭ，ｌと【０１７３】
【数８８】
0004322918.tif 000089
を用い、次式にしたがって誤差演算を行う。【０１７４】
【数８９】
0004322918.tif 000090
尤度計算部８２は、誤差演算を行った後、送信アンテナＭに対応する生成した
送信シンボルとそれに対応する誤差信号とを入力とし、誤差の小さいＳＭ個の送信アンテナＭに対応する送信シンボル系列【０１７５】
【数９０】
0004322918.tif 000091
とそのときの誤差【０１７６】
【数９１】
0004322918.tif 000092
を保持する。【０１７７】
次に、送信シンボル候補生成部８０は、送信アンテナＭ−１に対する送信シンボル【０１７８】
【数９２】
0004322918.tif 000093
を生成して出力する。【０１７９】
変換信号レプリカ生成部８１は、送信アンテナＭに対応するＳＭ個の送信シンボル系列と送信アンテナＭ−１に対する送信シンボルとからＲ行列を入力として変換信号レプリカ【０１８０】
【数９３】
0004322918.tif 000094
を次式にしたがって生成して出カする。【０１８１】
【数９４】
0004322918.tif 000095
尤度計算部８２は、次いで、誤差演算を次式により行い、【０１８２】
【数９５】
0004322918.tif 000096
誤差の小さいＳＭ−１個の送信アンテナＭと送信アンテナＭ−１に対応する送信シンボル系列の組み合わせ【０１８３】
【数９６】
0004322918.tif 000097
とそのときの誤差【０１８４】
【数９７】
0004322918.tif 000098
を保持する。【０１８５】
同様にして、送信シンボル候補生成部８０は、送信アンテナｍに対する送信シンボルｄｍ，ｉを生成して出力する。変換信号レプリカ生成部８１は、送信アンテナｍ＋１から、送信アンテナＭに対応するＳｍ+１個の送信シンボル系列と送信アンテナｍに対する送信シンボルとＲ行列を入力として変換信号レプリカ【０１８６】
【数９８】
0004322918.tif 000099
を次式にしたがって生成して出カする。【０１８７】
【数９９】
0004322918.tif 000100
尤度計算部８２は、誤差演算を次式にしたがって行い、【０１８８】
【数１００】
0004322918.tif 000101
誤差の小さいＳｍ個の送信アンテナＭから送信アンテナｍに対応する送信シンボル系列の組み合わせ【０１８９】
【数１０１】
0004322918.tif 000102
とそのときの誤差【０１９０】
【数１０２】
0004322918.tif 000103
を保持する。【０１９１】
以上の操作を繰り返すことにより、得られたＣ・Ｓ２個の全送信アンテナに対応する送信シンボル系列の組み合わせ【０１９２】
【数１０３】
0004322918.tif 000104
とそのときの誤差【０１９３】
【数１０４】
0004322918.tif 000105
を得る（ｊ＝１〜Ｃ・Ｓ２の整数）。【０１９４】
送信シンボル系列推定部８３１〜８３ｍでは、各送信アンテナに対応する、生き残り送信シンボルとそれに対応する誤差信号とを入力とし、最小誤差を選択し、当該誤差を与える送信シンボル系列を推定する。【０１９５】
送信シンボル系列推定部８３１〜８３ｍにおいて硬判定を行う場合の実施例は以下のようになる。【０１９６】
【数１０５】
0004322918.tif 000106
また、送信シンボル系列推定部８３１〜８３ｍにおいて軟判定を行う場合の実施例は以下のようになる。【０１９７】
まず、次式にしたがって、軟判定のビット系列を求める。【０１９８】
【数１０６】
0004322918.tif 000107
ここで、【０１９９】
【数１０７】
0004322918.tif 000108
は、送信アンテナｍのビットｉの対数尤度比である。【０２００】
【数１０８】
0004322918.tif 000109
上記式において、ｅｍｉｎ，νは、【０２０１】
【数１０９】
0004322918.tif 000110
の最小値を表し、σ２（ｎ）は雑音電力を表す。【０２０２】
軟判定シンボル【０２０３】
【数１１０】
0004322918.tif 000111
は、【０２０４】
【数１１１】
0004322918.tif 000112
を用いて以下のように推定される。【０２０５】
【数１１２】
0004322918.tif 000113
最後に、マルチパス受信信号再生処理部８４１〜８４ｍは、上記のようにして推定された送信シンボル系列【０２０６】
【数１１３】
0004322918.tif 000114
を入力し、次式に基づき、その入力された送信シンボル系列【０２０７】
【数１１４】
0004322918.tif 000115
に拡散符号とチャネル係数を乗算することにより、受信アンテナごとに各送信アンテナからの受信パス毎の受信信号系列【０２０８】
【数１１５】
0004322918.tif 000116
を推定して出力する。【０２０９】
【数１１６】
0004322918.tif 000117
以上説明したように、本実施形態におけるマルチパス受信信号再生部の構成は、図４に示すマルチパス受信信号再生部の構成と同様に、図３に示すマルチパス受信信号再生部の構成に比較して、若干の受信信号系列の再生精度の劣化を許容すれば、誤差演算の回数をＣＭ回から【０２１０】
【数１１７】
0004322918.tif 000118
回に大幅に低減することが可能であり、受信装置の複雑さを大幅に低減することが可能となる。【０２１１】
図６は、本発明に係るチャネル係数推定部の構成図であり、図７Ａおよび図７Ｂは、当該チャネル係数推定部を用いる場合の送信装置より送信される送信信号の構成例を示す図である。【０２１２】
まず、図７Ａおよび図７Ｂを参照しながら、送信装置から送信される送信信号について説明する。同図に示されるように、本実施形態では、各送信アンテナ（ここでは、送信アンテナ１、２）からの送信信号に、送信アンテナ毎に異なる４シンボル長のパイロットシンボル（斜線部）がデータシンボルに対して周期的に挿入されている。図７Ａに示される例では、各送信アンテナのパイロットシンボルパターンは互いに直交している。【０２１３】
また、図７Ｂに示される例では、各送信アンテナのパイロットシンボルの拡散に用いられる拡散符号（Ｃ１、Ｃ２）が直交している。【０２１４】
パイロット信号の送信は、図７Ａ、図７Ｂのいずれかの方法によって行えばよく、送信アンテナ間のパイロット信号を直交させることで高精度なチャネル推定が可能となる。上記のような直交シンボルパターンもしくは直交拡散符号は、例えば、パイロットシンボル数と同じ長さもしくは、パイロットシンボルの拡散率と同じ長さのＷａｌｓｈ系列を用いることで生成可能である。これ以降、送信アンテナmのパイロット信号系列をpｍ(n)として説明を進める。なお、nはチップ番号を表すものとする。【０２１５】
図６に戻り、本発明に係るチャネル係数推定部の構成を説明する。このチャネル係数推定部１００は、送信アンテナｍと受信アンテナｎ間の各パスのチャネル係数を推定する機能を備える。本例では、３送信アンテナ、４受信アンテナの場合のチャネル係数推定部の構成例を示している。すなわち、このチャネル係数推定部は、３×４の相関検出部１０１〜１０４と、パイロット信号レプリカ生成部１１１〜１１４とを備える。なお、本例では、構成要素が複数有る相関検出部と、パイロット信号レプリカ生成部の符号の末尾に連番を付すものとして図示した。【０２１６】
まず、同図を参照しながら送信アンテナ１と受信アンテナ１間のチャネル係数ｈ１，１，ｌを推定する場合の動作を説明する。【０２１７】
同図において、受信アンテナ１で受信された受信信号ｒ１は相関検出部１０１に入力される。また、パイロット信号レプリカ生成部１１１では送信アンテナ１のパイロットシンボル系列ｐ１を生成し、相関検出部１０１に入力する。【０２１８】
相関検出部１０１では、受信信号ｒ１に送信アンテナ１のパイロットシンボル系列p1の複素共役値をパスｌの受信タイミングを考慮して乗算した値を４パイロットシンボル区間で平均化することにより送信アンテナ１と受信アンテナ１間のチャネル係数ｈ１，１，ｌを次式にしたがって推定し出力する。【０２１９】
【数１１８】
0004322918.tif 000119
ここで、ｒ１（ｎ）は、パイロットシンボルｎが受信されるときの受信信号ｒ１を示す。実際にはｈ１，１，ｌの推定は複数のパイロットシンボル送信区間で得られたチャネル係数推定値を重み付け平均することで求めることも可能である。【０２２０】
同様にして、受信信号ｒ１を入力とする２段目の相関検出部（図示省略）では、受信信号ｒ１と２段目のパイロットシンボルレプリカ生成部（図示省略）で生成された送信アンテナｍのパイロットシンボル系列ｐｍを入力として、チャネル係数ｈ１，ｍ，ｌを推定し出力する。【０２２１】
さらに同様にして、受信信号ｒ４とパイロットシンボルレプリカ生成部１１３で生成される送信アンテナ１のパイロットシンボル系列ｐ１を相関検出部１０３に入力し、相関を求めることでチャネル係数ｈ４，１，ｌを推定し出力する。【０２２２】
以上の動作を繰り返すことで、３送信アンテナと４受信アンテナ間の各パスのチャネル係数を推定することができる。なお、上記では、パイロットシンボルがデータシンボルに時間的に多重される構成を例にあげて説明したが、符号多重を用いた場合も同様の方法でチャネル係数推定値を得ることができる。【０２２３】
図８は、本発明の実施形態に係る受信装置の第２の実施形態を示す構成図である。同図に示すように本実施形態では、受信装置２００は、複数のマルチパス受信信号再生部２１１〜２１３がシリアルにマルチパス干渉キャンセル部２２１１、２２１２、２２２１、２２２２、２２３１、２２３２を介して接続されており（本例では、３段構成）、最終段に復調部２３１が配置されている。初段のマルチパス受信信号再生部２１１には、上述した図２から図５に記載のいずれか任意のマルチパス受信信号再生部の構成を適用することができる。【０２２４】
また、２段目以降のマルチパス受信信号再生部２１２、２１３にも、以上で説明した図２から図５に記載のいずれか任意のマルチパス受信信号再生部の構成を適用することができる。ここで、２段目以降の第ｐ段目のマルチパス受信信号再生部への入力信号【０２２５】
【数１１９】
0004322918.tif 000120
は、受信信号と第p−１段目のマルチパス受信信号再生部の出力信号【０２２６】
【数１２０】
0004322918.tif 000121
を用いて、マルチパス干渉キャンセル部により次式の演算により生成される。【０２２７】
【数１２１】
0004322918.tif 000122
上記式にしたがい演算されて得られたＮ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２２８】
【数１２２】
0004322918.tif 000123
は受信アンテナｎ、パスｌに対応する逆拡散部（図示省略）に入力される。【０２２９】
また、この構成を用いた場合は、２段目以降の第ｐ段目のマルチパス受信信号再生部におけるチャネル係数推定部（図示省略）では、チャネル係数ｈ１１ｌの推定において、受信信号ｒｎ（ｔ）の代わりにマルチパス干渉キャンセル後の受信信号【０２３０】
【数１２３】
0004322918.tif 000124
を用いることでより高精度なチャネル係数推定を行うことができる（構成の簡単化のために、前段のマルチパス受信信号再生部におけるチャネル係数を用いることも同様に可能である）。【０２３１】
このように本実施形態によれば、マルチパス受信信号再生部を複数段設けることで、後段のマルチパス受信信号再生部では、マルチパス干渉キャンセル後の受信信号を用いて高精度にチャネル推定、送信シンボル系列推定を行うことができ、結果としてより高精度な受信アンテナごとの各送信アンテナからの受信パス毎の受信信号系列【０２３２】
【数１２４】
0004322918.tif 000125
を推定することができる。【０２３３】
図９は、本発明の実施形態に係る受信装置で用いられる復調部の第１の実施形態を示す構成図である。本実施形態では、復調部は、復調アルゴリズムとしてＭＬＤを用いる。【０２３４】
同図において、この復調部３００は、チャネル係数推定部３１１と、Ｎ×Ｌ個の逆拡散部３１２〜３１５と、送信シンボル候補生成部３１６と、受信信号レプリカ生成部３１７と、尤度計算部３１８と、送信系列推定部３１９とから構成される。なお、構成要素が複数有る場合には、末尾に連番を付すものとして図示した。【０２３５】
上記のように構成された復調部３００の動作について説明する。【０２３６】
復調部３００に入力される入力信号は、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２３７】
【数１２５】
0004322918.tif 000126
であり、チャネル係数推定部３１１において、受信アンテナと送信アンテナとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される（なお、構成の簡単化のために、前段のマルチパス受信信号再生部におけるチャネル係数を用いることも同様に可能である）。【０２３８】
さらに、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２３９】
【数１２６】
0004322918.tif 000127
は受信アンテナｎ、パスｌに対応する逆拡散部に入力され、逆拡散信号ｚ´ｎ、ｌを得る。送信シンボル候補生成部３１６は、各アンテナに対する送信シンボルｄｍ，ｉを生成して出力する。受信信号レプリカ生成部３１７は、送信シンボルとチャネル係数とを入力として受信信号レプリカ【０２４０】
【数１２７】
0004322918.tif 000128
を次式にしたがって生成して尤度計算部３１８に出力する。【０２４１】
【数１２８】
0004322918.tif 000129
尤度計算部３１８は、逆拡散部３１２〜３１５から出力される逆拡散信号ｚ´ｎ、ｌと受信信号レプリカ【０２４２】
【数１２９】
0004322918.tif 000130
とを入力として、誤差演算を次式にしたがって行う。【０２４３】
【数１３０】
0004322918.tif 000131
送信系列推定部３１９では、各送信アンテナに対応する生成した送信シンボルとそれに対応する誤差信号とを入力し、送信シンボル系列によって送信されたビッに対する尤度λｉを出力する。ここで、誤差信号に基づいたビット尤度の計算法は、既存のいかなる方法も適用可能である。【０２４４】
上述したビット尤度は、チャネル復号器（例えば、ターボ復号器）等に入力されて最終的に情報ビット系列が復元される。【０２４５】
図１０は、本発明の実施形態に係る復調部の第２の実施形態を示す構成図である。本実施形態では、復調部は、復調アルゴリズムとしてパス一括で処理を行うＭＬＤを用いる。【０２４６】
同図において、この復調部４００は、チャネル係数推定部４１１と、Ｎ×Ｌ個の逆拡散部４１２〜４１５と、ＱＲ分解部４１６と、ＱＨ演算部４１７と、送信シンボル候補生成部４１８と、変換信号レプリカ生成部４１９と、尤度計算部４２０と、送信系列推定部４２１とから構成される。なお、構成要素が複数有る場合には、末尾に連番を付すものとして図示した。【０２４７】
上記のように構成された復調部４００の動作について以下説明する。【０２４８】
復調部４００に入力される入力信号は、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２４９】
【数１３１】
0004322918.tif 000132
であり、チャネル係数推定部４１１において、受信アンテナと送信アンテナとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定される（なお、構成の簡単化のために、前段のマルチパス受信信号再生部におけるチャネル係数を用いることも同様に可能である）。【０２５０】
さらに、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２５１】
【数１３２】
0004322918.tif 000133
は受信アンテナｎ、パスｌに対応する逆拡散部に入力され、逆拡散信号ｚ´ｎ、ｌを得る。【０２５２】
次に、ＱＲ分解部４１６では、チャネル係数からなる下記のチャネル行列を生成し、チャネル行列のＱＲ分解を行って、Ｑ行列およびＲ行列をＱＨ演算部４１７に出力する。【０２５３】
【数１３３】
0004322918.tif 000134
上記ＱＲ分解部４１７から出力されるＱ行列は（Ｎ×Ｌ）行×Ｍ列のユニタリー行列であり、ＱＨＱ＝Ｉを満たす。また、Ｒ行列はＭ行×Ｍ列の上三角行列となる。したがって、ＱＨ演算部４１７における演算は、【０２５４】
【数１３４】
0004322918.tif 000135
と記述することができる。【０２５５】
送信シンボル候補生成部４１８は、まず、送信アンテナＭに対する送信シンボルｄＭ，ｊを生成して出力する。変換信号レプリカ生成部４１９は、送信シンボルｄＭ，ｊとＲ行列を入力して変換信号レプリカ【０２５６】
【数１３５】
0004322918.tif 000136
を次式にしたがって生成して尤度計算部４２０に出力する。【０２５７】
【数１３６】
0004322918.tif 000137
尤度計算部４２０はまず、ｘ´Ｍと、【０２５８】
【数１３７】
0004322918.tif 000138
の誤差演算を次式にしたがって行う。【０２５９】
【数１３８】
0004322918.tif 000139
尤度計算部４２０では、送信アンテナＭに対応する生成した送信シンボルとそれに対応する誤差信号とを入力とし、誤差の小さいＳＭ個の送信アンテナＭに対応する送信シンボル系列【０２６０】
【数１３９】
0004322918.tif 000140
とそのときの誤差【０２６１】
【数１４０】
0004322918.tif 000141
次に、送信シンボル候補生成部４１８、送信アンテナＭ−１に対する送信シンボル【０２６２】
【数１４１】
0004322918.tif 000142
を生成して出力する。【０２６３】
変換信号レプリカ生成部４１９は、送信アンテナＭに対応するＳＭ個の送信シンボル系列と送信アンテナＭ−１に対する送信シンボルとからＲ行列を入力として変換信号レプリカ【０２６４】
【数１４２】
0004322918.tif 000143
を次式にしたがって生成して出カする。【０２６５】
【数１４３】
0004322918.tif 000144
尤度計算部４２０は、次いで、誤差演算を次式により行い、【０２６６】
【数１４４】
0004322918.tif 000145
誤差の小さいＳＭ−１個の送信アンテナＭと送信アンテナＭ−１に対応する送信シンボル系列の組み合わせ【０２６７】
【数１４５】
0004322918.tif 000146
とそのときの誤差【０２６８】
【数１４６】
0004322918.tif 000147
を保持する。【０２６９】
同様にして、送信シンボル候補生成部４１８は、送信アンテナｍに対する送信シンボルｄｍ，ｉを生成して出力する。変換信号レプリカ生成部４１９は、送信アンテナｍ＋１から、送信アンテナＭに対応するＳｍ+１個の送信シンボル系列と送信アンテナｍに対する送信シンボルとＲ行列を入力として変換信号レプリカ【０２７０】
【数１４７】
0004322918.tif 000148
を次式にしたがって生成して出カする。【０２７１】
【数１４８】
0004322918.tif 000149
尤度計算部４２０は、誤差演算を次式にしたがって行い、【０２７２】
【数１４９】
0004322918.tif 000150
誤差の小さいＳｍ個の送信アンテナＭから送信アンテナｍに対応する送信シンボル系列の組み合わせ【０２７３】
【数１５０】
0004322918.tif 000151
とそのときの誤差【０２７４】
【数１５１】
0004322918.tif 000152
を保持する。【０２７５】
以上の操作を繰り返すことにより、得られたＣ・Ｓ２個の全送信アンテナに対応する送信シンボル系列の組み合わせ【０２７６】
【数１５２】
0004322918.tif 000153
とそのときの誤差【０２７７】
【数１５３】
0004322918.tif 000154
を得る。【０２７８】
送信系列推定部４２１では、各送信アンテナに対応する、生き残り送信シンボルとそれに対応する誤差信号とを入力とし、送信シンボル系列によって送信されたビットに対する尤度λｉを出力する。なお、誤差信号に基づいたビット尤度の計算法は、既存のいかなる方法も適用可能である。【０２７９】
上述したビット尤度はチャネル復号器（例えば、ターボ復号器）等に入力されて最終的に情報ビット系列が復元される。【０２８０】
以上のように、第２の実施形態に係る復調部の構成（図１０）は、図９において説明した第１の実施形態の復調部の構成に比較にして、若干のビット尤度の推定精度の劣化を許容すれば、誤差演算の回数をＣＭ回から【０２８１】
【数１５４】
0004322918.tif 000155
回に大幅に低減することが可能であり、受信装置の複雑さを大幅に低減することが可能となる。【０２８２】
図１１は、本発明の実施形態に係る復調部の第３の実施形態を示す構成例図である。本実施形態では、復調部は、復調アルゴリズムとしてパス毎に処理を行うＭＬＤを用いる。【０２８３】
同図において、この復調部５００は、チャネル係数推定部５１１と、Ｎ×Ｌ個の逆拡散部５１２〜５１５と、Ｌ個のＱＲ分解部５１６、５１７と、Ｌ個のＱＨ演算部５１８、５１９と、送信シンボル候補生成部５２０と、変換信号レプリカ生成部５２１と、尤度計算部５２２と、送信系列推定部５２３とから構成される。なお、構成要素が複数有る場合には、末尾に連番を付すものとして図示した。【０２８４】
上記のように構成された復調部５００の動作について以下説明する。【０２８５】
復調部５００に入力される入力信号は、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２８６】
【数１５５】
0004322918.tif 000156
であり、チャネル係数推定部５１１において、受信アンテナと送信アンテナとの間各々の各パスのチャネル係数ｈｍ，ｎ，ｌが推定さ
される（なお、構成の簡単化のために、前段のマルチパス受信信号再生部におけるチャネル係数を用いることも同様に可能である）。【０２８７】
さらに、Ｎ×Ｌ個のマルチパス干渉キャンセル後の受信信号【０２８８】
【数１５６】
0004322918.tif 000157
は受信アンテナｎ、パスｌに対応する逆拡散部に入力され、逆拡散信号ｚ´ｎ、ｌを得る。【０２８９】
次に、第ｌ番目のＱＲ分解部では、第ｌ番目のパスのチャネル係数からなる下記のチャネル行列をパス数分生成し、それぞれ、チャネル行列のＱＲ分解を行って、Ｑ行列およびＲ行列を出力する。【０２９０】
【数１５７】
0004322918.tif 000158
上記ＱＲ分解部から出力されるＱｌ行列は、Ｎ行×Ｍ列のユニタリー行列であり、ＱｌＨＱｌ＝Ｉを満たす。ここで、Ｈは共役複素転置を表し、Ｉは単位行列を表す。また、Ｒｌ行列は、Ｍ行×Ｍ列の上三角行列となる。【０２９１】
したがって、第ｌ番目のＱＨ演算部における演算は、【０２９２】
【数１５８】
0004322918.tif 000159
と記述することができる。【０２９３】
送信シンボル候補生成部５２０は、まず、送信アンテナＭに対する送信シンボルｄＭ，ｉを生成して変換信号レプリカ生成部５２１に出力する。変換信号レプリカ生成部５２１は、送信シンボルｄＭ，ｉと、Ｒｌ行列を入力して、変換信号レプリカ【０２９４】
【数１５９】
0004322918.tif 000160
を次式にしたがって生成して出力する。【０２９５】
【数１６０】
0004322918.tif 000161
尤度計算部５２２は、まず、ｘ´Ｍ，ｌと【０２９６】
【数１６１】
0004322918.tif 000162
を用い、次式にしたがって誤差演算を行う。【０２９７】
【数１６２】
0004322918.tif 000163
尤度計算部５２２は、誤差演算を行った後、送信アンテナＭに対応する生成した送信シンボルとそれに対応する誤差信号とを入力とし、誤差の小さいＳＭ個の送信アンテナＭに対応する送信シンボル系列【０２９８】
【数１６３】
0004322918.tif 000164
とそのときの誤差【０２９９】
【数１６４】
0004322918.tif 000165
を保持する。【０３００】
次に、送信シンボル候補生成部５２０は、送信アンテナＭ−１に対する送信シンボルｄＭ−１，ｉを生成して出力する。【０３０１】
変換信号レプリカ生成部５２１は、送信アンテナＭに対応するＳＭ個の送信シンボル系列と送信アンテナＭ−１に対する送信シンボルとからＲ行列を入力として変換信号レプリカ【０３０２】
【数１６５】
0004322918.tif 000166
を次式にしたがって生成して出カする。【０３０３】
【数１６６】
0004322918.tif 000167
尤度計算部５２２は、次いで、誤差演算を次式により行い、【０３０４】
【数１６７】
0004322918.tif 000168
誤差の小さいＳＭ−１個の送信アンテナＭと送信アンテナＭ−１に対応する送信シンボル系列の組み合わせ【０３０５】
【数１６８】
0004322918.tif 000169
とそのときの誤差【０３０６】
【数１６９】
0004322918.tif 000170
を保持する。【０３０７】
同様にして、送信シンボル候補生成部５２０は、送信アンテナｍに対する送信シンボルｄｍ，ｉを生成して出力する。変換信号レプリカ生成部５２１は、送信アンテナｍ＋１から、送信アンテナＭに対応するＳｍ+１個の送信シンボル系列と送信アンテナｍに対する送信シンボルとＲ行列を入力として変換信号レプリカ【０３０８】
【数１７０】
0004322918.tif 000171
を次式にしたがって生成して出カする。【０３０９】
【数１７１】
0004322918.tif 000172
尤度計算部５２２は、誤差演算を次式にしたがって行い、【０３１０】
【数１７２】
0004322918.tif 000173
誤差の小さいＳｍ個の送信アンテナＭから送信アンテナｍに対応する送信シンボル系列の組み合わせ【０３１１】
【数１７３】
0004322918.tif 000174
とそのときの誤差【０３１２】
【数１７４】
0004322918.tif 000175
を保持する。【０３１３】
以上の操作を繰り返すことにより、得られたＣ・Ｓ２個の全送信アンテナに対応する送信シンボル系列の組み合わせ【０３１４】
【数１７５】
0004322918.tif 000176
とそのときの誤差【０３１５】
【数１７６】
0004322918.tif 000177
を得る。【０３１６】
送信系列推定部５２３では、Ｌ個の尤度計算部から得られた各送信アンテナに対応する、生き残り送信シンボルとそれに対応する誤差信号とを入カとし、Ｌ個の誤差信号の和から、送信シンボル系列によって送信されたビットに対する尤度
λｉを出カする。なお、誤差信号に基づいたビット尤度の計算法は、既存のいかなる方法も適用可能である。【０３１７】
上述したビット尤度は、チャネル復号器（例えば、ターボ復号器）等に入力されて最終的に情報ビット系列が復元される。【０３１８】
以上のように、第３の実施形態に係る復調部の構成（図１１）は、図９において説明した第１の実施形態の復調部の構成に比較にして、若干のビット尤度の推定精度の劣化を許容すれば、誤差演算の回数をＣＭ回から【０３１９】
【数１７７】
0004322918.tif 000178
回に大幅に低減することが可能であり、受信装置の複雑さを大幅に低減すること
が可能となる。【０３２０】
図１２は、本発明の実施形態に係る受信装置の第３の実施形態を示す構成図である。本実施の形態による無線通信システムにおいては、送信装置（図示省略）において、送信データビット系列をＭ本の送信系列にシリアル・パラレル変換し、同一の周波数帯かつ同一の拡散符号群を用いてデータ変調を行い、Ｍ本の送信アンテナからＮｃｏｄｅ本のマルチコードチャネルで同時に送信される場合を示している。【０３２１】
本実施形態に係る受信装置６００では、ステージ数分の各コードチャネル（本例では、コードチャネル１、２）に対応するマルチパス受信信号再生部６１１〜６１４と、各コードチャネルに対応する復調部６２１、６２２を用意し、マルチパス干渉キャンセル部６１５〜６１８を介して接続される。  
本実施形態では、初段の各コードチャネルに対応するマルチパス受信信号再生部６１１、６１２への入力信号は、Ｎ本の受信アンテナ（ここでは、Ｎ＝２の例を示している）で受信された受信信号ｒｎ（ｔ）であり、マルチパス受信信号再生部６１１、６１２は、それぞれ当該コードチャネルの各送信アンテナからの受信パス毎の再生受信信号系列【０３２２】
【数１７８】
0004322918.tif 000179
（ｋはコード番号で、【０３２３】
【数１７９】
0004322918.tif 000180
はコードチャネルｋの拡散符号を表す）を出力する。【０３２４】
Ｎ個のマルチパス干渉キャンセル部（本例では、Ｎ＝２）６１５、６１６には、受信信号ｒｎ（ｔ）と全コードチャネル再生受信信号系列【０３２５】
【数１８０】
0004322918.tif 000181
とが加算器６３１〜６３４により加算（合成）されて入力され、マルチパス干渉キャンセル後の受信信号【０３２６】
【数１８１】
0004322918.tif 000182
を次式にしたがって生成する。【０３２７】
【数１８２】
0004322918.tif 000183
第２段の各コードチャネルに対応するマルチパス受信信号再生部６１３、６１４への入力信号は、マルチパス干渉キャンセル後の受信信号【０３２８】
【数１８３】
0004322918.tif 000184
となり、以上の操作をステージ数だけ繰り返して、最終的なマルチパス干渉キャンセル後の受信信号【０３２９】
【数１８４】
0004322918.tif 000185
が生成される。【０３３０】
このようにして生成されたマルチパス干渉キャンセル後の受信信号【０３３１】
【数１８５】
0004322918.tif 000186
は各コードチャネルに対応する復調部６２１、６２２に入力され、各コードチャネルにおける送信シンボル系列によって送信されたビットに対する尤度が出力される。【０３３２】
次に、本発明について具体的に計算機シミュレーションを行った例を、図１３を用いて説明する。同図は、従来のＭＬＤと２次元ＭＭＳＥおよび本発明による受信装置の復調方法を用いた場合の受信アンテナ当たりの平均受信Ｅｂ／Ｎｏ（情報１ビット当りの信号エネルギーに対する雑音電力密度）に対するスループット特性を計算機シミュレーションにより評価した場合の結果を示す図である。【０３３３】
本シミュレーションでは、受信装置のマルチパス受信信号再生部を、２ステージで構成し、第１ステージでは、図２の構成を、第２ステージでは図４の構成を用いるものとした。また、マルチパスとして平均受信電力が等しい２パスモデルを仮定した。拡散率は１６であり、１５コードチャネルを符号多重している（実効的な拡散率は１５／１６である）。図中Ｒはターボ符号化を用いたチャネル符号化の符号化率を表し、実線は、本発明法による特性（ＱＲ−ＭＬＤ  ｗｉｔｈ  ２−ｓｔａｇｅ  ＭＰＩＣ）を、点線は従来におけるＭＭＳＥの特性を、×印は従来におけるＭＬＤの特性を示している。また、●、○はＱＰＳＫ変調を用いたＭＩＭＯ多重（４送信アンテナ、４受信アンテナ）の受信アンテナ当たりの平均受信Ｅｂ／Ｎｏを、■、□は１６ＱＡＭ変調を用いたＭＩＭＯ多重の受信アンテナ当たりの平均受信Ｅｂ／Ｎｏ特性をそれぞれ表している。【０３３４】
同図により、本発明の受信装置の構成を用いることで、従来の信号分離法を用いる場合に比較して、大幅にあるスループットを得るために必要な平均受信電力Ｅｂ／Ｎｏを低減できていることがわかるが、これは、本発明の受信装置の構成を用いることでより大幅に少ない送信電力で従来と同じスループットを実現できることを示している。換言すれば、同じ送信電力であれば従来よりも大幅にスループットを増大できることを示している。【０３３５】
以上説明してきたように、本発明によれば、ＣＤＭＡ方式において情報ビットレートを増大するために複数の送信アンテナから異なるデータを同時送信するＭＩＭＯ多重を適用したときに、マルチパス干渉を低減し、異なる送信アンテナから送信された信号の高精度な分離を実現することができる。その結果、マルチパスフェージング環境における、受信ビット誤り率・受信パケット誤り率の低減、並びにスループット（誤り無く伝送できる情報ビットの伝送レート）を大幅に向上させることができる。【０３３６】
上記各実施形態では、マルチパス受信信号再生部に、ＭＭＳＥやＭＬＤなどのアルゴリズム（例えば、アルゴリズムを記述したプログラム）が採用される場合を例示したが、より好ましくは、実用性の範囲の演算量で高精度な信号分離を実現可能とするＱＲ−ＭＬＤアルゴリズムを採用するとよい。【０３３７】
また、上記各実施形態では、無線通信システムで用いられる受信装置を適用したが、当該受信装置を、移動通信システムで用いている基地局において適用することも勿論可能である。また、無線回線や有線回線を介して上記アルゴリズムを選択的にダウンロードし、無線設備の特性を変化させる、いわゆるソフトウェア無線基地局を本発明に適用することも可能である。
【産業上の利用可能性】
【０３３８】
本発明は、複数アンテナを用いて信号の復調を行う無線通信システムに適用可能である。
【図１】
0004322918.tif 000187
【図２】
0004322918.tif 000188
【図３】
0004322918.tif 000189
【図４】
0004322918.tif 000190
【図５】
0004322918.tif 000191
【図６】
0004322918.tif 000192
【図７Ａ】
0004322918.tif 000193
【図７Ｂ】
0004322918.tif 000194
【図８】
0004322918.tif 000195
【図９】
0004322918.tif 000196
【図１０】
0004322918.tif 000197
【図１１】
0004322918.tif 000198
【図１２】
0004322918.tif 000199
【図１３】
0004322918.tif 000200
【図１４】
0004322918.tif 000201
【図１５】
0004322918.tif 000202