(57)【特許請求の範囲】
【請求項１】
複数のサブキャリアを介して少なくとも2ユーザの信号を送受する通信システムにおいて、
異なるコードで拡散された2ユーザの信号をコード多重し、複数のサブキャリアを用いて送信する送信装置、
信号を受信し、受信信号を送信側コードで逆拡散し、非線形処理して得られる各ユーザの受信データの軟判定を行う受信部を含む受信装置、
各受信部における軟判定の対象値を他方の受信部に入力する手段、
を備え、一方の受信部は他方の受信部から入力された軟判定対象値を用いて、自身の軟判定対象値を調整し、該軟判定対象値に基づいて受信データを判定する、
ことを特徴とする通信システム。
【請求項２】
前記各受信部は、
受信したデータが2値のうち一方である確率と他方である確率との差を前記軟判定対象値として演算する手段、
他方の受信部から入力された軟判定対象値を用いて、自身の軟判定対象値を調整する手段、
該軟判定対象値に基づいて受信データを判定する判定部
を備えたことを特徴とする請求項1記載の通信システム。
【請求項３】
前記受信装置は、
信号を各サブキャリアにおける信号に変換するFFT部、
FFT部からの出力信号を用いて各ユーザデータを判定する受信部を備え、
各受信部は、
各サブキャリアの受信シンボルに、送信側拡散コードと同じコードパターンと各サブキャリアに関連するチャネル推定値とを結合した信号を乗算して逆拡散する4つの乗算器、
振幅リミッター特性に応じた伝達関数を備えた非線形ユニット、
閾値に基づいて送信データについて判定を行う判定部、
前記逆拡散信号の2つの第1加算部、
一方の第1加算部からの出力信号と他ユーザデータの軟判定対象値との加算を行い、加算結果を前記非線形ユニットに入力する第2加算部、
該非線形ユニットの出力信号と前記他方の第1加算部からの出力信号を加算して軟判定対象値を出力する手段、
を備えたことを特徴とする請求項1記載の通信システム。
【請求項４】
複数のサブキャリアを介して少なくとも2ユーザの信号を送受する通信システムの受信装置において、
信号をサブキャリア毎の信号に変換するFFT部、
サブキャリア毎の信号に送信側と同じコードを用いて逆拡散処理、非線形処理を施してユーザデータの判定対象値を発生し、該判定対象値の軟判定を行う受信部、
各受信部における前記判定対象値を他方の受信部に入力する手段、
を備え、一方の受信部は他方の受信部から入力された軟判定対象値を用いて、自身の軟判定対象値を調整し、該軟判定対象値に基づいて受信データを判定する、
ことを特徴とする通信システムの受信装置。
【請求項５】
複数のサブキャリアを介して少なくとも2ユーザの信号を送受する通信システムの受信装置において、
信号を各サブキャリアにおける信号に変換するFFT部、
FFT部からの出力信号を用いて各ユーザデータを判定する受信部を備え、
各受信部は、
各サブキャリアの受信シンボルに、送信側拡散コードと同じコードパターンと各サブキャリアに関連するチャネル推定値とを結合した信号を乗算して逆拡散する4つの乗算器、
振幅リミッターの特性に応じた伝達関数を備えた非線形ユニット、
閾値に基づいて送信データについて判定を行う判定部、
前記逆拡散信号の2つの第1加算部、
一方の第1加算部からの出力信号と他ユーザデータの軟判定対象値との加算を行い、加算結果を前記非線形ユニットに入力する第2加算部、
該非線形ユニットの出力信号と前記他方の第1加算部からの出力信号を加算して軟判定対象値を出力する手段、
を備えたことを特徴とする通信システムの受信装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、2以上のサブチャネルを採用する通信システム、特にDMT(discrete multi-tone)変調やフィルタバンク(filter-bank)変調、などの幾つかのサブチャネルを備えた通信システム、あるいは、マルチキャリア通信システム、あるいは、マルチキャリア変調によりバンドが独立した狭バンドに分割され、各サブバンドが、多重ユーザのアクセスをサポートするために、独立してCDMA方式により変調される通信システムに関する。
【０００２】
この10年の間に、マルチキャリア(MC)と符号分割多重アクセス（CDMA）の組合わせに基づいた多くの新しい変調方式が考えられてきた。これら方式の共通の特徴は、拡散データシーケンス(チップ)の並列送信であり、チップ周期は、マルチパスフェージングに対するロバスト性を得るために実際に増大させられてきた。マルチキャリア(MC)と符号分割多重アクセス（CDMA）の組合わせは、異なったアプローチで実現されている。対応する方式は、MC-CDMA、直接拡散CDMA (DS-CDMA)、マルチトーンCDMA (MT-CDMA)が知られている。全方式において、ユーザは同時に利用可能なバンド幅を共用する。ユーザ信号の分離は、DS-CDMAに類似した方法でコードドメインにおいて実行される。
周波数選択性フェージング環境は、サブチャネルにおいてマルチアクセス干渉(MAI)またはチャネル間干渉(ICI)を引き起こし、結果的にBER性能を劣化する(たとえば、非特許文献1、２)。
【０００３】
本発明では、MC-CDMAべースの通信システムのために最適なシンボルワイズレシーバ(symbol-wisereceiver)の構成を導出する。Markov チェーンアプローチ（Markov chain approach）を使って、干渉が隣接サブチャネルに影響を与える状況において最適レシーバをシンセサイズする。新しい非線形信号処理技術が、すなわち、最大事後確率推定に基づいた非線形信号処理技術が提案される。最大事後確率推定は、システムのBER(bit error rate)性能に対する多重アクセス干渉（MAI）あるいはチャネル間干渉(ICI)の効果をかなり緩和する。ほとんどの研究において、MAIまたはICIを零平均ガウス分布確率変数（Gaussian-distributed randomvariable with a zero mean）として扱うことが、普通に行われている。これに対して、我々は、MAIまたはICIを記述するために有限状態離散マルコフプロセスモデル(a finite-statediscrete Malkov process model)を採用する。2ユーザデータの送受信について説明するが、いくらかの修正により、以後に説明される本発明のアイデアは、幾人かのユーザがMAIあるいはICIにより互いに影響されるより複雑なケースに発展可能である。
【０００４】
【従来の技術】
DMTあるいはOFDMによるMC-CDMA方式を考察する。MC-CDMA方式は、DS拡散とMC変調との直列連鎖に基づいている。議論を単純化するために、IFFT(図6)とFFT(図７)は理想的に動作し、付加的なICIを引き起こさないものとする(たとえば、非特許文献3)。IFFTとFFTは、MC-CDMAシンボルを変調信号に変換する際、あるいは逆の変換をする際、システムにおいて実施される。
図6はMC-CDMA送信機の構成図であり、ユーザ数が2で、サブキャリア数が２の場合である。ユーザ０のコピー部CPY0はユーザ0の送信データD0をコピーして２つのサブキャリア用のデータを出力する。乗算器MPL00は第１複製データD0に第1拡散コード列｛C00，C01｝の第1コードC00を乗算し、乗算器MPL01は第2複製データD0に第２コードC01を乗算し、それぞれサブチャネル0，1側に出力する。
【０００５】
ユーザ１のコピー部CPY1はユーザ1の送信データD1をコピーして２つのサブキャリア用のデータを出力する。乗算器MPL10は第1複製データD１に第２拡散コード列｛C10，C11｝の第1コードC10を乗算し、乗算器MPL11は第２複製データD1に第２コードC11を乗算し、それぞれサブチャネル0，1側に出力する。
サブチャネル0(第1サブキャリア)側の加算器ADD0は
D0・C00＋D1・C10
を演算して出力し、サブチャネル１(第２サブキャリア)側の加算器ADD1は
D0・C01＋D1・C11
を演算して出力する。
逆高速フーリエ変換部IFFTは2つのサブキャリア信号に逆高速フーリエ変換（あるいは逆離散フーリエ変換）を施して実部及び虚部の時間軸信号に変換する。各時間軸信号は図示しないGI挿入部、QPSK変調部を介してアンテナより送信される。
【０００６】
図７はMC-CDMA受信機の構成図であり、ユーザ数が2で、サブキャリア数が２の場合である。フーリエ変換部FFTは図示しない受信部、復調部、GI除去部を介して入力されたべースバンドの受信信号にフーリエ変換処理を施してサブチャネル0,1の2つのサブキャリア信号S0，S1を出力する。チャネル推定部CEUはサブキャリア信号S0，S1に含まれるパイロット信号を用いて各サブキャリアのチャネルを推定し、チャネル補償部CC0,CC1はチャネル推定値k0，ｋ1をサブキャリア信号S0，S1に乗算してチャネル補償(フェージング補償)を行う。
【０００７】
乗算器MP00は第１キャリアデータＳ0に第1拡散コード列｛C00，C01｝の第1コードC00を乗算し、乗算器MP01は第2キャリアデータＳ1に第２コードC01を乗算し、それぞれサブチャネル0の加算器AD0に入力する。サブチャネル0(第1サブキャリア)側の加算器AD0は次式
Ｄ0＝Ｓ0・C00＋Ｓ1・C01
を演算して送信データＤ0を復調する。
同様に、乗算器MP10は第１キャリアデータＳ0に第2拡散コード列｛C10，C11｝の第1コードC10を乗算し、乗算器MP11は第2キャリアデータＳ1に第２コードC11を乗算し、それぞれサブチャネル1の加算器AD1に入力する。サブチャネル１(第2サブキャリア)側の加算器AD1は次式
Ｄ1＝Ｓ0・C10＋Ｓ1・C11
を演算して送信データＤ1を復調する。
【０００８】
IFFTとFFTが理想的に動作し、付加的なICIを引き起こさないものとすれば、以降の説明においてそれらを省略することができる。IFFTを省略した場合の2ユーザのMC-CDMA用送信機は図８に示され、FFTを省略した場合の2ユーザのMC-CDMA用受信機は図９に示され、それぞれ図6、図７と同一部分には同一符号を付している。図８において、増幅部AF0、AF1の係数k0，k1は振幅に関しての第1、第2サブキャリア０，１のチャネル伝達関数であり、伝送環境(チャネル)におけるサブキャリアの減衰を意味している。
これからの説明において、記号を簡単化にするために、信号Di(t)=±１やSi(t)の時間依存性を省略し、新しい記号としてDi(t)をDｉで、Si(t)をSｉで表記する。
【０００９】
サブキャリア(サブチャネル)数Nは２である。各サブキャリアシーケンスDｉ(ｉはサブキャリア番号)は、拡散ファクタSFの直交ショートコードCiにより拡散され、ついで、同一のサブキャリアを共用するために他ユーザからの対応するサブキャリアと結合される。図８において、C0＝｛C00，C01｝は第1ユーザ(ユーザ番号0)の拡散コードであり、C1＝｛C10，C11｝は第２ユーザ(ユーザ番号1)の拡散コードである。直交する一組の拡散コード列C0、C1を発生するために使用可能な多くの異なるシーケンスが存在するが、ウォルシュ(Walsh)やハダマード(Hadamard)シーケンスはCDMAにおける有益な組を与える。ここで、ハダマードマトリクスから拡散コードを選択する。表記において、ｃの第1インデックスはユーザ番号(ユーザ数2の場合には0又は１)を表わし、第2インデックスはサブキャリア(サブチャネル)番号を示す。図８に示されるMC-CDMAシステムは、サブキャリア数Ｎ，拡散ファクタSFがすべて等しいと仮定しており、N=SF=2である。図８のモデルは、MAIやICIを引き起こす物理プロセスを理解するために有益である。
【００１０】
MC伝送時、各サブキャリア上に周波数非選択フェージングが本質的に存在し、また、一般的にN=SFとする必要はない。実際に、初期のシンボルレートが周波数選択フェージングの支配を受けるのに十分に高ければ、信号を周波数ドメインで拡散する前に、直列／並列変換する必要がある。
MC-CDMAレシーバは、逆拡散を成功させるためにコヒーレントな検出を必要とする。本明細書では、理想ガウシャン白色ノイズチャネル(AWGN)と理想的なサブチャネルのフレーム／時間同期を仮定する。レシーバにおいて、各サブキャリアSiの各受信シンボルは、拡散シーケンスCiの対応するチップパターン(コードパターン)の信号によって乗算され、推定されたチャネル推定値ｋが各サブキャリアに関連付けられる。それらは図９に示されるように周波数ドメインにおいてSFチップの期間にわたってコヒーレントに積算される。図９は2ユーザの場合の従来技術（たとえば、非特許文献２，３）で報告されたMC-CDMAを示している。このレシーバは”クラシカル”MC-CDMAとして本明細書で後に参照される。理想的なチャネル推定ができたと仮定し、各サブキャリアのチャネル推定値がk0，k1であるとしている。
【００１１】
図式的例示として、ユーザへの拡散コード割り当てのためにウォルシュコードＣ＝｛w0，w1｝が使用されたものとする。最初のユーザ(第1ユーザ)の拡散コードはＣ0＝｛＋１，＋１｝となり、第2ユーザの拡散コードはＣ1＝｛＋１，−１｝となる。送信機における第1、第2サブキャリア信号は、
S0＝D0＋D１
S1＝D0−D１
となり、チャネル出力は
S0＝ｋ0・（D0＋D１）
S1＝ｋ1・（D0−D１）              (1)
となる。レシーバにおいて、各サブキャリアSiの信号は逆拡散されてコヒーレントに結合される。フラットフェージングの場合、逆拡散及び結合後、データ推定は(2)式の様に得られる。
【数１】
0004187555.tif 000002
(2)式において、k0，k1は、第1、第2サブチャネルの振幅に関するチャネル伝達係数である。フラットフェージングの場合には、k0＝k1と仮定することができる。この簡単な例から明らかなように、フラット周波数フェージングの場合、適当な拡散コードを選択することにより、受信シンボルDiに付加する歪を発生するICI又はMAIは存在しない。これは、MC-CDMAべース通信システムにおいてユーザデータ分離の基本的な原理である。
【００１２】
以下では、もはやフラットフェージングを仮定しない。このため、k0≠k1である。(2)式から明らかなように、
【数２】
0004187555.tif 000003
が成立する。(3)式において、被乗数(k02−k12)をもったメンバーは、周波数選択性フェージングによるMAIを表わしている。(3)式におけるマルチアクセス干渉MAIは、ICIとしても扱われる。実際、(3)式からわかるように、被乗数(k02−k12)はクロスフィード係数又は2つのサブチャネル間の漏れを表わしている。
C0＝｛ｋ0・ｗ00，ｋ1・ｗ01｝，C1＝｛ｋ0・ｗ10，ｋ1・ｗ11｝としての拡散コードパターンに対して、あるいは図6に示す記号を用いた(4)式の拡散コードパターンに対して、サブチャネル伝達係数k0，k1が再計算可能であることを示すことは容易である。
【数３】
0004187555.tif 000004
(4)式により拡散コードを修正して、周波数選択性フェージングをフラットフェージングとして扱うことができる。一般的な場合におけるC＝｛ｗ0，ｗ1｝と対比すると、(4)式のコードは非直交である。図10及び図11はこのアプローチによるMC-CDMA送信機及びMC-CDMA受信機の構成図である。
【非特許文献１】
S.Hara, P Prasad “Overview of Multicarrier CDMA”, IEEE CommunicationMagazine, Dec.1997, pp.126-133.
【非特許文献２】
S.Hara “Design and Performance of Multicarrier CDMA System inFrequency-Selective Raleigh Fading”, IEEE Transaction on Vehicular Technology,Vol.48, Sep. 1999.
【非特許文献３】
S.Abeta, H.Atarashi, N.Sawahashi and F.Adachi “Performance of CoherentMulticarrier/DS-CDMA and MC-CDMA for Broadband Packet Wireless Access”, IEICETrans. Commun., Vol. EB-84, No.3, March 2001, pp.406-414.
【００１３】
【発明が解決しようとする課題】
(3)式からk0=k1のとき、MAI又はICIが存在しなくなる。実際、この場合、零に等しいクロスフィード係数及びサブチャネル間の漏れは除去される。他のどのようなkの選択であってもMAIあるいはICIの原因になる。この簡単な例は、注意しないとk0≠k1の時、周波数選択フェージングがサブチャネルにMAIあるいはICIを引き起こし、BERを劣化する。
従って、本発明の目的は、MC-OFDMの通信システムのBER性能を改善することである。そのような改善は非線形信号処理の実施により得られる。すなわち、非線形処理後に一つのサブチャネルから導出される情報が他サブチャネルの推定事後確率を精練するために使用されるとき、BERの改善が得られる。
【００１４】
【課題を解決するための手段】
本発明は、マルチアクセス方式と結合された2-サブチャネルマルチキャリア通信システムであり、サブチャネルを介して独立にデータを発生する2つのデータソースを備え、これらデータのシーケンスは周波数ドメインにおいて拡散ファクタSF=2の直交コードシーケンスにより拡散され、ついで、IFFT変換され、周波数選択性環境を通して送信される。IFFTは各サブキャリアにおけるシンボルを信号に変換するために使用される。
すなわち、本発明は、複数のサブキャリアを介して少なくとも2ユーザの信号を送受する通信システムであり、異なるコードで拡散された2ユーザの信号をコード多重し、複数のサブキャリアを用いて送信する送信装置、信号を受信し、受信信号を送信側コードで逆拡散し、非線形処理して得られる各ユーザの受信データの軟判定を行う受信部を含む受信装置、各受信部における軟判定の対象値を他方の受信部に入力する手段を備え、一方の受信部は他方の受信部から入力された軟判定対象値を用いて、自身の軟判定対象値を調整し、該軟判定対象値に基づいて受信データを判定する。
【００１５】
受信装置は、好ましくは、▲１▼信号を各サブキャリアにおける信号に変換するFFTプロセッサ、▲２▼各サブキャリアのシンボルに、コードパターンと各サブキャリアに関連するチャネル推定値とを結合した信号を乗算する4つの乗算器、▲３▼振幅リミッター特性に応じた伝達関数を備えた非線形ユニット、▲４▼閾値に基づいて送信データについて判定を行う判定部、▲５▼前記逆拡散信号の2つの第1加算部、▲６▼一方の第1加算部からの出力信号と他ユーザデータの軟判定対象値との加算を行い、加算結果を前記非線形ユニットに入力する第2加算部、▲７▼該非線形ユニットの出力信号と前記他方の第1加算部からの出力信号を加算して軟判定対象値を出力する手段、を備えている。本発明によれば、BER性能を改善することができる。
【００１６】
【発明の実施の形態】
（Ａ）サブキャリアの受信シンボル復調アルゴリズム
図１はMC-CDMA送信部の構成図であり、ユーザ数が2で、サブキャリア数(サブチャネル数)が２の場合である。なお、IFFTは省略し図示していない。
ユーザ０のコピー部11はユーザ0の送信データD0をコピーして２つのサブキャリア用のデータを出力する。乗算器12は第１複製データD0に第1拡散コード列｛C00，C01｝の第1コードC00を乗算し、乗算器13は第2複製データD0に第２コードC01を乗算し、それぞれサブチャネル0，1側に出力する。
ユーザ１のコピー部21はユーザ1の送信データD1をコピーして２つのサブキャリア用のデータを出力する。乗算器22は第1複製データD１に第1拡散コード列｛C10，C11｝の第1コードC10を乗算し、乗算器23は第２複製データD1に第２コードC11を乗算し、それぞれサブチャネル0，1側に出力する。
サブチャネル0(第1サブキャリア)側の加算器14は
D0・C00＋D1・C10
を演算して出力し、サブチャネル１(第２サブキャリア)側の加算器24は
D0・C01＋D1・C11
を演算して出力する。
【００１７】
伝送中に第1サブキャリア信号にノイズn0(t)が付加されてｙ0(t)となって受信機31に受信される。また、第２サブキャリア信号にノイズn1(t)が付加されてｙ1(t)となって受信機に受信される。ノイズn0(t)，n1(t)は統計的に独立である。
バイナリ情報(2値情報) Diがサブチャネル０，１を介して信号S*ij(t)として送信されるものとする。なお、S*ij(t)におけるインデックスｉはユーザ番号(ｉ＝０，１)を示し、インデックスjは情報シンボルDiの符号により決定される。すなわち、
Di=＋1ならばj=0
Di=−1ならばj=１
である。以後、表記を簡単にするために、式においてS*ij (t)の時間依存性を省略する。すなわち、S*ij(t)をS*ij と表記する。
【００１８】
送信情報シンボルD1、D2は統計的に独立で(相関がなく)、且つ、等分布された確率変数であるとする。図１から、ICIの影響を受けた信号は、サブチャネル０，１で送信された信号S*0j，S*1j (j=0，１)の線形結合として表現される。
【数４】
0004187555.tif 000005
ICIの導入後は、２つの送信信号Diに替わって(5)，(6)式に従って、サブチャネル０の送信信号としてSｉ(ｉ=0,1,2,3)を使用し、サブチャネル１の送信信号としてSｉ(ｉ=4,5,6,7)を使用する。送信信号Siにおける番号ｉは、第1、第2サブチャネルにおけるシンボルD1,D2をペア(対)にすることにより決定される信号番号を示す。
【００１９】
サブチャネルｉの受信情報シンボルDｉがjとなる事後確率P(Di=j/y(t))は、
Dｉがjである信号を受信する事後確率として得ることができる。この事後確率
P(Di =j/y(t))は、2つのサブチャネルにおいてDi情報シンボルをそれぞれ送信する事後確率の乗積の和として得られる。すなわち、第1サブチャネル(サブチャネル０)の受信情報シンボルD0が”0”(=+1)となる事後確率P(D0=+1/y(t))、および、第１サブチャネルの受信情報シンボルD0が”1”(=−1)となる事後確率P(D0=−1/y(t))はそれぞれ(7)式で求めることができる。
【数５】
0004187555.tif 000006
(7)式より、D0=+1となる事後確率は、第1サブチャネルの信号がS0で、第2サブチャネルの信号がS4である確率の積と、第1サブチャネルの信号がS1で、第2サブチャネルの信号がS5である確率の積の和である。同様に、D0=−1となる事後確率は、第1サブチャネルの信号がS2で、第2サブチャネルの信号がS6である確率の積と、第1サブチャネルの信号がS3で、第2サブチャネルの信号がS7である確率の積の和である。
【００２０】
受信信号がSiである事後確率P(Si/y(t)) は、次式
【数６】
0004187555.tif 000007
により与えられる。ただし、
k0は正規化因子、
ｉはサブチャネル番号(ｉ=1,2)、jは信号番号(j=0,1,2,3)、
y(t)は第ｉサブチャネルにおいて、ICIを伴う信号系列とスペクトルパワー強度N0を有する白色ガウスノイズn(t)との合成信号（y(t)= Sｉ +n(t)）、
Papr(Si)は受信信号Siの事前確率、
P(y(t)/ Si)は条件付き確率であり、受信語がy(t)であった時、送られた符号語がSiであったという確率である。
【００２１】
(7)式に(8)式を用いて変形すると、(9)式が得られる。
【数７】
0004187555.tif 000008
事前確率Papr(Si)(ｉ=0〜7)は(5)(6)式より(10)式で示すように2つの隣接サブチャネルで送信されるデータの交差積として表現される。
【数８】
0004187555.tif 000009
第１サブチャネル(サブチャネル0)の受信情報シンボルD0が”0”(=+1)、”1”(=−1)となる事後確率P(D0=＋1/y(t))、P(D0=−1/y(t))が求まれば、それらの大小比較により、あるいはそれらの対数(logarithm)の差と閾値との比較により受信情報シンボルの符号(+1又は−1)を決定できる。
【００２２】
・大小比較による判定
第1サブチャネルの情報シンボルD0が+1であるか、−１であるかは、まず、
【数９】
0004187555.tif 000010
【００２３】
・対数の差による判定
第1サブチャネルの情報シンボルD0が+1であるか、−１であるかは、まず、
【数１０】
0004187555.tif 000011
を演算し(lnはeを底とする対数)、しかる後、その正負により判定する。すなわち、
【数１１】
0004187555.tif 000012
であればD１＝＋１と判定し、
【数１２】
0004187555.tif 000013
であればD１＝−１であると判定する。
(11c),(11d)式における各対数項は(12a),(12b)式に示すようになる。但し、表記を簡単にするために、P(Si/y(t))→P(Si)とする。
【数１３】
0004187555.tif 000014
(12a),(12b)式のln P(D0=+1/y(t))，ln P(D0=−1/y(t))を求めると(13)〜(14)式になる。ただし、ln P(D0=+1)＝ln P(D0=+1/y(t))，ln P(D0=−1)＝ln P(D0=−1/y(t))である。
【００２４】
【数１４】
0004187555.tif 000015
lnは対数、chはcosh（ハイパボリク コサイン）を意味する。
【００２５】
さて、(8)式より、
【数１５】
0004187555.tif 000016
が成立しており、両辺の対数を演算すると
【数１６】
0004187555.tif 000017
また、(5)、(6)式より、
S0=−S3，  S1=−S2，  S4=−S7，  S5=−S6
であるから、エネルギーに関して
Ｅ0=Ｅ3，  Ｅ1=Ｅ2，  Ｅ4=Ｅ7，  Ｅ5=Ｅ6        (16)
が成立する。従って、(17)，(18)式が成立する。
【数１７】
0004187555.tif 000018
【数１８】
0004187555.tif 000019
次式
【数１９】
0004187555.tif 000020
に、(13),(14)式を代入して(17)、(18)式を考慮すると(19)式が導出される。
【００２６】
【数２０】
0004187555.tif 000021
上式において、ΔlnP0，ΔlnP１は、第1、第2ユーザのソフト判定値で
【数２１】
0004187555.tif 000022
である。ｙ0(t)，ｙ1(t)は、第1サブチャネル、第2サブチャネルの受信信号で
ｙ0(t)=Si(t)+n0(t)（ｉ＝０〜３）
ｙ1(t)=Si(t)+n1(t)（ｉ＝４〜７）
である。又、
ΔE0=E0−E1，ΔE4=E4−E5
である。
【００２７】
（Ｂ）受信機の構成
図2は(19)式の判定アルゴリズムを実現する受信機の構成図であり、ソフト判定値ΔlnP0演算部50、FFT演算部51、ソフト判定値ΔlnP1演算部70を有している。第1サブチャネルのソフト判定値ΔlnP0演算部50のみの詳細を示すが第2サブチャネルの構成も同様である。
(19)式における(S0+S1)，(S4+S5)，(S0−S1)，(S4−S5)は(5)，(6)式より、
(S0＋S1)=2・C00，
(S4＋S5)=2・C01，
(S0−S1)=2・C10，            (20)
(S4−S5)=2・C11
である。尚、C00〜C11は(4)式で与えられる。
【００２８】
(20)式の関係およびFFT演算部51から
【数２２】
0004187555.tif 000023
が出力されるものとすると、乗算器52は(19)式の右辺第1項
【数２３】
0004187555.tif 000024
の演算を行い、乗算器53は(19)式の右辺第2項
【数２４】
0004187555.tif 000025
の演算を行う。乗算器54は(19)式の右辺第３項、第4項における
【数２５】
0004187555.tif 000026
の演算を行い、乗算器55は(19)式の右辺第３項、第4項における
【数２６】
0004187555.tif 000027
の演算を行う。加減算器56〜58は(19)式の右辺第３項のln coshの［］内の演算を行い、加減算器56〜57,59は(19)式の右辺第4項のln coshの［］内の演算を行う。なお、
【数２７】
0004187555.tif 000028
である。ln cosh演算部60は、(19)式の右辺第３項の演算を行い、ln cosh演算部61は、(19)式の右辺第4項の演算を行う。加算器62〜64は(19)式右辺の演算を行ってソフト判定値ΔlnP0を出力する。判定部65はこのソフト判定値ΔlnP0の正負に基づいて受信信号D0の”0”，”1”を判定して出力する。
【００２９】
ソフト判定値ΔlnP1演算部70も同様にソフト判定値ΔlnP1を演算し、該ソフト判定値ΔlnP1の正負に基づいて受信信号D1の”0”，”1”を判定して出力する。以上では、ln cosh演算部60,61で(19)式の右辺第３項、第4項の演算を行う例を示したであるが、図３に示す振幅リミッターの特性（伝達関数）を備えた非線形ユニットを設けることにより、ln cosh演算部60,61を省略することができる。なお、図3において、縦軸は、
【数２８】
0004187555.tif 000029
である。
【００３０】
以上要約すれば、図2のMA-CDMA用のレシーバは、▲１▼受信MC-CDMA信号を各サブキャリアにおけるCDMA信号に変換するFFTプロセッサ、▲２▼サブキャリアの各受信MC-CDMAシンボルに、チップパターン(コードパターン)と各サブキャリアに関連するチャネル推定値とを結合した信号((4)式参照)を乗算する4つの乗算器、▲３▼振幅リミッター特性に応じた伝達関数を備えた非線形ユニット、▲４▼閾値に基づいて送信データについて判定を行う判定部、▲５▼前記逆拡散信号の2つの第1加算部、▲６▼一方の第1加算部からの出力信号と他ユーザデータの軟判定対象値との加算を行い、加算結果を前記非線形ユニットに入力する第2加算部、▲７▼該非線形ユニットの出力信号と前記他方の第1加算部からの出力信号を加算して軟判定対象値を出力する手段を備えている。
【００３１】
本発明の受信データの復調アルゴリズムはターボ符号のターボデコーダに類似している。ターボデコーダにおいて、各デコーダは情報を他のデコーダに渡し、そして、他のデコーダにより導き出された情報を用いて順番に推定された事後確率を精練する。同様に、本発明のアルゴリズムにおいても、非線形処理後に導出される第1ユーザデータの軟判定対象値が、第2ユーザデータの推定事後確率を精練するために使用され、非線形処理後に導出される第2ユーザデータの軟判定対象値が、第1ユーザデータの推定事後確率を精練するために使用され、BERの改善効果が得られる。
【００３２】
（C）ノイズイミュニティとシミューレーション結果
一般的な場合、本発明レシーバは事後確率の推定を改善するために非線形ユニットとクロスフィードを採用した。この構成は、BER性能を解析することが非常に難しいことを意味している。それゆえ、非線形信号処理がクラシカルなレシーバより優れていることを証明するためにコンピュータシミューレーションに頼った。図４において特性Ａは、本発明レシーバの平均BERパフォーマンスを2Eb/N0の関数として示している。ただし、k0＝0.5，k1＝1.0である。また、参考として特性Bは、クラシカルなレシーバでMAIが存在する時のBERパフォーマンスを示す。特性ＣはMAI(ICI)のない本発明のBER特性である。特性Ｄは、(21)式の公式を使って計算したマッチトフィルタレシーバのICIが存在しない時のBERシミューレーション結果を表示する。
【数２９】
0004187555.tif 000030
である。
コンピュータシミューレーションにより得られた本発明のBERパフォーマンスと(21)式で計算されたBERパフォーマンスとはかなり良く一致している。又、図４のプロットから明らかなようにMAIが存在しない本発明レシーバのBERは、従来のマッチトフィルタべースレシーバの(21)式で得られたBERと差異がない。尚、後者のBERは図４で”Reference”として示している。
【００３３】
図5はk1=1.0に固定し、k0を0〜1,0の範囲で可変した時の本発明レシーバのBER特性ＡとクラシカルなレシーバのBER特性Ｂである。ただし、2Eb/N0=10dBである。図５より、本発明レシーバは広い範囲にわたって良好なBERパフォーマンスが得られた。しかしながら、BERにおける最大の改善は大きな2Eb/N0で発生する(図4参照)。これは、2Eb/N0が小さいと、入力ノイズがMAI(ICI)を支配し、事後確率の推定を改善するために利用できる信頼性のあるデータが存在しないからである。これに対して、2Eb/N0が大きくなると、MAI(ICI)がノイズを支配するからである。以上のことはｋ0の範囲に対しても同様である。相対的に伝達係数ｋ0(ｋ0<0.5)を小さくして第1サブチャネル信号のEb/N0比を低くすると、推定後に得られたデータは全く信頼性が低い。しかし、第1サブチャネル信号のEb/N0比を高くすると、得られたデータの信頼性を高くできる。
【００３４】
【発明の効果】
以上の本発明によれば、一方の受信部は他方の受信部から入力された軟判定対象値を用いて、自身の軟判定対象値を調整し、該軟判定対象値に基づいて受信データを判定するように構成したから、BER性能を改善することができる。
【図面の簡単な説明】
【図１】  MC-CDMA送信部の構成図である。
【図２】本発明の判定アルゴリズムを実現するMC-CDMA受信機の構成図である。
【図３】非線形ユニットの特性図である。
【図４】本発明のBERパフォーマンス説明図である。
【図５】  k1=1.0に固定し、k0を0〜1,0の範囲で可変した時の本発明レシーバのBER特性である。
【図６】従来のMC-CDMA送信機の構成図である。
【図７】従来のMC-CDMA受信機の構成図である。
【図８】従来のMC-CDMA送信機の第1変形図である。
【図９】従来のMC-CDMA受信機の第1変形図である。
【図１０】従来のMC-CDMA送信機の第２変形図である。
【図１１】従来のMC-CDMA受信機の第２変形図である。
【符号の説明】
５０  ソフト判定値ΔlnP0演算部
５１  FFT演算部
６０，６１  ln cosh演算部
６５  判定部
７０  ソフト判定値ΔlnP1演算部
【図１】
0004187555.tif 000031
【図２】
0004187555.tif 000032
【図３】
0004187555.tif 000033
【図４】
0004187555.tif 000034
【図５】
0004187555.tif 000035
【図６】
0004187555.tif 000036
【図７】
0004187555.tif 000037
【図８】
0004187555.tif 000038
【図９】
0004187555.tif 000039
【図１０】
0004187555.tif 000040
【図１１】
0004187555.tif 000041