JP4490265B2 - 同一チャネル干渉受信機 - Google Patents

同一チャネル干渉受信機 Download PDF

Info

Publication number
JP4490265B2
JP4490265B2 JP2004522953A JP2004522953A JP4490265B2 JP 4490265 B2 JP4490265 B2 JP 4490265B2 JP 2004522953 A JP2004522953 A JP 2004522953A JP 2004522953 A JP2004522953 A JP 2004522953A JP 4490265 B2 JP4490265 B2 JP 4490265B2
Authority
JP
Japan
Prior art keywords
module
signal
user
receiver
received
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP2004522953A
Other languages
English (en)
Other versions
JP2005534226A (ja
Inventor
タイラー、マシュー・エー
マクレオド、ロバート・ビー
ラーンド、レイチェル・イー
ニードズウィッキー、ジョシュア・ディー.
ブロマー、カール・ディー
マッセルワイン、トーマス・ピー
Original Assignee
ビーエイイー・システムズ・インフォメーション・アンド・エレクトロニック・システムズ・インテグレイション・インコーポレーテッド
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US10/228,787 external-priority patent/US6947502B2/en
Application filed by ビーエイイー・システムズ・インフォメーション・アンド・エレクトロニック・システムズ・インテグレイション・インコーポレーテッド filed Critical ビーエイイー・システムズ・インフォメーション・アンド・エレクトロニック・システムズ・インテグレイション・インコーポレーテッド
Publication of JP2005534226A publication Critical patent/JP2005534226A/ja
Application granted granted Critical
Publication of JP4490265B2 publication Critical patent/JP4490265B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/7103Interference-related aspects the interference being multiple access interference
    • H04B1/7105Joint detection techniques, e.g. linear detectors
    • H04B1/71055Joint detection techniques, e.g. linear detectors using minimum mean squared error [MMSE] detector
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/10Polarisation diversity; Directional diversity
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0242Channel estimation channel estimation algorithms using matrix methods
    • H04L25/0246Channel estimation channel estimation algorithms using matrix methods with factorisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0242Channel estimation channel estimation algorithms using matrix methods
    • H04L25/0248Eigen-space methods
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03203Trellis search techniques
    • H04L25/03216Trellis search techniques using the M-algorithm
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03203Trellis search techniques
    • H04L25/03222Trellis search techniques using the T-algorithm
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03305Joint sequence estimation and interference removal
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03312Arrangements specific to the provision of output signals
    • H04L25/03318Provision of soft decisions
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03331Arrangements for the joint estimation of multiple sequences
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L2025/0335Arrangements for removing intersymbol interference characterised by the type of transmission
    • H04L2025/03426Arrangements for removing intersymbol interference characterised by the type of transmission transmission using multiple-input and multiple-output channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • H04L25/0228Channel estimation using sounding signals with direct estimation from sounding signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03184Details concerning the metric
    • H04L25/03197Details concerning the metric methods of calculation involving metrics
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03248Arrangements for operating in conjunction with other apparatus
    • H04L25/03286Arrangements for operating in conjunction with other apparatus with channel-decoding circuitry

Description

(省益の説明書)
本発明の部分は、政府の財政的支援とともに行われた可能性があり、政府には特定の権利がある可能性がある。
(関連出願)
本願は、2002年7月24日に出願された米国仮出願第60/398,451号の利点を主張する。本願は、2002年4月16日に出願された米国仮出願第60/372,956号に対する優先権を主張する2002年8月26日に提出された米国出願第10/228,787号の一部継続出願である。本願は2002年3月25日に提出された米国出願第10/105,918号の一部継続出願でもある。本願は「マルチユーザ検出のための周波数不一致補償(Frequency Mismatch Compensation for Multiuser Detection)」(弁護士事件整理番号20020028−US)と題される2003年4月25日に出願された米国出願(番号未定)に関連する。本願は、「過飽和状態の通信用の延期されたDecorrelating決定フィードバック検出器(Deferred Decorrelating Decision−Feedback Detector For Supersaturated Communications)」(弁護士事件整理番号D4625−US)と題される2003年4月25日に提出された米国出願(番号未定)に関連している。これらの出願のそれぞれは全体として参照によりここに組み込まれる。
本発明は無線通信に関し、さらに詳しくは、同様に変調された干渉信号からの同一チャネル干渉により壊されたデジタル信号の回復が望ましい応用例用の無線デジタル信号受信機に関する。
無線ネットワークは、携帯電話とコンピュータのような多様な装置間の通信のために利用されている。バイナリ位相偏移変調された信号及び直角位相偏移変調された信号などのデジタル変調された信号はネットワークのノード間で送信される。例は、端末が衛星トランスポンダを通って送信する衛星通信ネットワーク、端末が中継塔を通して伝送する地上システム、及び端末が中央中継要素を通って送信する屋内ローカルエリアネットワークを含む。
これらのネットワークに接続されるコンピュータ要素は、種々のユーザサービスを提供する。例はデジタル音声符号化を使用する電話トラフィック、ビデオ会議、広域コンピュータネットワーク接続性、及びインターネットサービスを含む。このような応用例では、干渉と雑音の存在下で指定された帯域幅でネットワークトラフィック容量を最大にすることが望ましい。
その目的のため、ネットワーク要素を通信チャネルに効率的に分割するための種々の方式が存在している。例えば、周波数領域多元接続(FDMA)方式は各ネットワーク端末を別個の重複しない周波数帯に割り当てる。時間領域多元接続(TDMA)方式は、各端末を別個の重複しないタイムスロットに割り当てる。符号分割多元接続(CDMA)方式は、各端末間の相互相関がごくわずかであるように各端末を別個の変調波形に割り当てる。
このような方式は、十分な信号対雑音比またはコード化冗長性を考慮すると、干渉の存在下で信号を復調するための手段だけではなく、過剰な信号マージンを検出するための手段も備える場合に、さらに多くの通信機構が割り当てられた帯域幅を使用できるという点で効率的ではない。要するに、無線送信及び受信における進歩にも関わらず、従来のシステムは、干渉とマルチパスの問題などの信号の劣化を被る現実世界の無線通信信号を適切に捕らえていない。
さらに具体的には、現実世界のマルチユーザシステムは同時に信号を送信する多数の独立したユーザを含んでいる。これらの伝送のそれぞれは、受信された信号の中で現れるマルチパスと同一チャネル(co-channel)干渉というリアルタイムの問題に関連付けられている。マルチパスは、受信機がさまざまな無作為に変化する遅延と振幅を有するエコーに遭遇するように、一つではなく多くの経路に沿って信号が受信機に進むときに発生する。同一チャネル干渉は、他のユーザから受信される信号を指す。
マルチユーザ検出(MUD)受信機は、同一チャネル干渉デジタル信号をいっしょに復調するために使用できる。一般的には、MUDは非直交マルチプレックスでのデータの検出を指す。MUD処理はチップごとに使用可能な情報ビット数または干渉が制限されたシステムのためのシグナリング次元を増加する。最大尤度原則に基づいた最適MUDは、受信機で受信される波形を生じさせるために、送信機群で発生したかもしれない可能性の全体数に、受信された信号を比較することにより動作する。
しかしながら、さらに大きな容量の信号を調べるマルチユーザ検出器にとって、計算は
複雑であり、時間を要し、従ってリアルタイム動作を実行不可能にする。従来のツリープ
ルーニング(tree-pruning)アルゴリズムに基づいた複雑さが低減された手法はある程度ま
で役に立つ。しかしながら、このようなマルチユーザ検出アルゴリズムの性能は、パラメ
ータM(プルーニング(pruning)係数)が減少するにつれて劣化するが、Mは必要とされ
る計算数を支配する。したがって、不適切なプルーニングに対抗するには、基本的なツリ
ープルーニングが、Mが十分に大きいことを確実にしなければならない。その結果、従来
のプルーニング方法は、特に干渉信号の数が中程度から大であるときに、依然として計算量の多い高い複雑性と関連している。
したがって、必要とされているのは、複数のユーザが同じ通信チャネルで動作できるよ
うにする技法である。このような技法は、同一チャネル信号を正確に分離し、複雑な処理
を削減しなければならない。
〔発明が解決しようとする課題〕
本発明は、ツリープルーニング・アルゴリズムに基づくマルチユーザ検出と比較して、より計算量の少ないマルチユーザ検出法を提供することを目的とする。
本発明の一実施形態は、同一チャネル干渉受信機を提供する。該受信機はK個の同一チャネル干渉信号からの情報を含む合成(complex)信号を受信するように適応されるマルチユーザ検出器モジュールを含む。該受信機はさらに、該合成信号を受信し、該K個の同一チャネル干渉信号のそれぞれの推定シグナチャ(signature)波形を生成するように適応されるパラメータ推定モジュールを含む。該推定されたシグナチャ波形はマルチユーザ検出器に与えられ、それによりK個の同一チャネル干渉信号の復調を可能にする。マルチユーザ検出器モジュール及びパラメータ推定モジュールがそれぞれ複数の合成信号を受信するように適応できることに注意せよ。
マルチユーザ検出器モジュールは、例えば、部分的に量子化された過去の情報を有するアルゴリズムと動作するように構成される。代替として、マルチユーザ検出器モジュールは事前白色化された(prewhitened)データに基づいたアルゴリズムと動作するように構成される。代りに、受信機は両方のアルゴリズムと動作するように構成されてよい。1つのこのような実施形態では、部分的に量子化された過去の情報を有するアルゴリズムは低複雑性線形(low complexity linear)MMSEアルゴリズムである。別のこのような実施形態では、事前白色化されたデータに基づくアルゴリズムは、M−アルゴリズムとT−アルゴリズムの1つである。
受信機は、さらに、マルチユーザ検出器モジュールとパラメータ推定モジュールに動作可能に結合されるアナログフロントエンド(analog front end)を含んでもよい。該アナログフロントエンドは、それぞれ複数の送信機からの1つまたは複数の合成波形信号を受信し、それぞれの受信された合成波形信号を対応する合成信号に変換するように適応される。1つのこのような実施形態では、アナログフロントエンドは、それぞれが受信された合成波形をデジタル波形に変換するように適応される1または複数のアナログ/デジタル変換器を含み、それぞれがそれぞれのアナログ/デジタル変換器に動作可能に結合され、受信された合成波形に関連する周波数をさらに低い周波数に変換するように適応される1または複数のダウンコンバータも含む。
アナログフロントエンドは、さらに、それぞれが複数の送信機から対応する合成波形信号を受信するように構成される1または複数のアンテナを含んでもよい。1つのこのような実施形態では、該1または複数のアンテナは単一偏向アンテナである。代替として、該1または複数のアンテナは2つの偏向ポートと適応された二重偏向アンテナであり、それにより偏向(polarization)として、該1または複数のアンテナは、それぞれが2つの偏向ポートで適応される2または3以上の二重偏向アンテナを含み、それにより空間ダイバシティと偏向ダイバシティを提供する。
マルチユーザ検出器モジュールが部分的に量子化された過去の情報を有する低複雑性線形MMSEアルゴリズムと動作するように構成される実施形態では、マルチユーザ検出器モジュールは、ターボMUDモジュールと、コンバイナモジュールと、エラー訂正モジュールと、スレッショルド(thresholding)モジュールとを含む。該ターボMUDモジュールはK個の同一チャネル干渉信号のそれぞれに個々のビットの推定値を提供するように適応される。ビットに関連する誤り率が所定の数字以下に低下するまで、該推定値はエラー訂正モジュールを含むフィードバックループで反復適用される。該コンバイナモジュールはターボMUDに動作可能なように結合され、ターボMUDにより出力される再計算されたビット推定値を次の反復で量子化されたビット値と結合するように適応される。該スレッショルドモジュールは、エラー訂正モジュールの出力に動作可能なように結合され、所定の閾値を超える各ビット推定値に量子化された値を割り当て、コンバイナモジュールまでそれらの量子化されたビット値を通過し、それにより部分的に量子化された過去の情報を有効にする。1つのこのような実施形態では、以後のそれぞれの反復プロセスでの該エラー訂正モジュールは、ターボMUDにより出力された再計算されたビット推定値とスレッショルドモジュールにより出力される量子化されたビット値の組み合わせを処理し、スレッショルドモジュールを通してターボMUDにその出力を戻し、それにより反復ごとの不確かなビット推定値の数を削減する。
マルチユーザ検出器モジュールが事前白色化されたデータに基づいた低複雑性Mアルゴリズムと動作するように構成される実施形態では、マルチユーザ検出器モジュールは整合(matched)フィルタと、ホワイトナデザイナ(whitener designer)モジュールと、非同期ホワイトナモジュールと、シンボル仮説試験(symbol hypothesis testing)モジュールとを含む。該整合フィルタは受信機により受信される合成信号を事前白色化するように適応され、それにより多元接続干渉からユーザを部分的に減結合する(decoupling)。該ホワイトナデザイナモジュールはパラメータ推定器に動作可能に結合され、パラメータ推定器からのパラメータ推定値に基づいた各受信合成信号のモデルを作成し、整合フィルタにより出力されるフィルタリングされたデータを白色化する非同期ホワイトナモジュールを計算するように適応される。該シンボル仮説試験モジュールはホワイトナデザイナモジュールに動作可能に結合され、非同期ホワイトナモジュールにより出力される白色化されたデータを受信するように構成される。シンボル仮説試験モジュールは、仮説の尤度を特徴付ける測定基準(metric)の連続した評価に基づきシンボル仮説試験を実施するように適応される。1つのこのような実施形態では、ホワイトナデザイナモジュールは平方根因数分解を利用する。例えば、ホワイトナデザイナモジュールは白色化モジュールでの白色化のために使用され、シンボル仮説試験モジュールの仮説試験でも使用される対角線上で(diagonally)でロードされる(loaded)コレスキー因数分解(Cholesky Factorization)を計算するためにパラメータ推定モジュールにより提供される相関行列(correlation matrix)を活用する。代替として、ホワイトナデザイナモジュールは、ハウスホルダ変換(Householder transformations)を使用するQR因数分解を利用する。さらに、ホワイトナデザイナモジュールは受信されたエネルギー及び/または位相がシンボル期間の間でしか変化しないときに非同期ホワイトナモジュールを効率的に更新するために双曲線ハウスホルダ(Hyperbolic Householder)変換を利用できる。
パラメータ推定モジュールのパラメータ推定値は、それぞれの受信された複合信号と関連付けられるチャネルをモデル化するために使用でき、それにより整合フィルタの適用及び非同期デコレレイト(decorrelating)フィルタバンクの作成を可能にする。一実施形態では、パラメータ推定モジュールはトレーニングシーケンスロケータモジュールと、雑音推定器モジュールと、シグナチャ波形推定器モジュールと、アクティブユーザテスタモジュールと、変換行列リビルダ(rebuilder)モジュールとを含む。該トレーニングシーケンスロケータモジュールは、受信された合成信号の各フレーム内のトレーニングシーケンスロケーション指数を推定するように適応される。該雑音推定器モジュールは、該トレーニングシーケンスロケーション指数に従って該受信された合成信号内の平均雑音電力の推定値を計算するように適応される。該シグナチャ波形推定器モジュールは、該トレーニングシーケンスロケーション指数と変換行列に従って該受信された合成信号内の各ユーザに一意のシグナチャ波形を推定するように適応される。アクティブユーザテスタモジュールは、雑音推定器モジュールの出力と、シグナチャ波形推定器モジュールの出力とに動作可能に結合され、受信された合成信号に関連付けられるアクティブユーザの数を判断するように適応される。該変換行列リビルダモジュールは、ユーザテスタモジュールと、ユーザごとに事前に記憶された既知のトレーニングシーケンスとに動作可能に結合され、シグナチャ波形推定モジュールにより使用される変換行列を生成するように適応される。
本発明の別の実施形態は、K個の同一チャネル干渉信号からの情報を含む合成信号を受信するための方法を提供する。該方法は、該K個の同一チャネル干渉信号のそれぞれのシグナチャ波形を推定することと、部分的に量子化された過去の情報を有する低複雑性線形MMSEアルゴリズムと、事前白色化されたデータに基づいた低複雑性Mアルゴリズムの少なくとも1つでシグナチャ波形に基づいた合成信号を処理することを含む。1つのこのような実施形態では、部分的に量子化された過去の情報を有する低複雑性線形MMSEアルゴリズムは、それぞれの処理反復から所定の閾値を上回る推定値を有するそれらのビットの検討を排除することを含み、閾値を上回るビット推定値は確実と考えられる。別のこのような実施形態では、事前白色化されたデータに基づいた低複雑性Mアルゴリズムは合成信号をフィルタリングすることを含み、それによりユーザを多元干渉から部分的に切り離し(decoupling)、事前白色化されたデータを提供する。方法は、パラメータ推定値に基づいて受信された合成信号のモデルを作成することと、事前白色化されたデータを白色化するためのモデルに基づき非同期ホワイトナを計算することと、仮説の尤度を特徴付ける測定基準の連続した評価に基づいて記号仮説試験を実施することに進む。
ここに説明されている特徴及び優位点は包括的ではなく、特に多くの追加の機能及び優位点は、図面、明細書及び請求項を考慮して当業者に明らかとなるであろう。さらに、明細書で使用される言語がおもに読みやすさと教育上の目的のために選択され、本発明の主題の範囲を制限するためではないことが注意されなければならない。
今日の通信システムでは、セントラルコントローラは、通常、通信チャネルセットアップ期間の間に1つの通信機構をそれぞれのチャネルに割り当てる。チャネルはタイムスロット、周波数及び拡散コードのなんらかの組み合わせであってよい。大部分のシステムにおいて、チャネルは離れた地域で再利用され、それにより同一チャネル(co-channel)干渉を引き起こす。再利用が発生する地域の間に大きな距離がある場合には、信号減衰は同一チャネル干渉を許容できるレベルまで削減する。従来の受信機は、かなりの同一チャネル干渉の存在下では信号を復調することができないため、これは必要な結果である。
本発明の原理に従って構成される同一チャネル干渉受信機は、同じチャネル上で送信される2つまたは3つ以上の信号をいっしょに復調できる。同一チャネル干渉受信機の実施形態を活用するシステムは類似するチャネルセットアップフォーマットを使用できるであろうが、チャネル割り当ては任意の指定領域内のチャネルごとに単一通信機構(communicator)に制限されないであろう。さらに具体的には、すべての使用可能なチャネルがチャネルごとに一人のユーザで満たされると、セントラルコントローラは新しい通信機構をすでに占有されているチャネルに追加することによって新しいチャネル要求を満たし始めることができる。
これは第1の通信機構のビット誤り率をわずかに劣化させるであろう。セントラルコントローラはオプションでこのチャネル上の送信機にビット誤り率を下げるために送信電力を増加するように命令してもよい。代わりに、セントラルコントローラは、オプションでこのチャネル上の送信機と受信機に帯域幅を減少し、ビット誤り率を下げるように命令してもよい。
さらに、本発明の受信機は空間的にダイバースアンテナ及びまたは偏向ダイバースアンテナを使用してデータを受信するように構成される。単一アンテナの使用だけではなく、両方の概念とも本発明の範囲内にある。複数の偏向ポートを包含することにより、送信機は両方の偏向で、及びかなりのマルチパスまたは電磁散乱の場合にも送信できるようになり、たとえ1つだけが送信されたとしても両方の偏向とも受信、処理されてもよい。偏向及び/または空間ダイバシティが信号空間の次元数を増加し、配列点間の距離を効果的に拡大する。結果として、ビット誤り率は改善される。
本発明の実施形態は複数の応用例で使用できる。1つのこのような応用例は、システムが、ほぼ同一の搬送周波数で同様に変調されたデータを送信する複数の独立した通信網から信号を同時に受信する空地(airborne−to−ground)通信用同一チャネル通信システムである。主要なビームエリアカバレージは低密度(例えば、10,000平方キロメートル)または高密度(例えば100平方キロメートル)である場合がある。例えば、28の基地局をカバーする主要ビームと、基地局あたり8人のユーザがいて、基地局の内の4つが同じ周波数で動作している、7という周波数再利用係数を考えてみる。このようなシステムは224人までのユーザのダウンリンク信号を受信し、多様なソースから生じる種々の信号タイプは単一受信機により処理される必要があるであろう。
同一チャネルシステムの別の応用例は、周波数再利用を利用する同じ通信網の要素から信号を同時に受信する地上周波数通信受信機である。通信網は例えばパケット無線ネットワーク、携帯電話網、または無線ローカルエリアネットワークである場合がある。ネットワーク要素の不注意な位置決めのために、ネットワークは干渉により質が落ちる。
ここでの説明は無線通信応用例を説明しているが、本発明の原理は通信回線に接続されたケーブルシステム及びローカルエリアネットワーク、ディスクドライブまたは他の記憶装置の読み書き動作、衛星通信、及び多くの複数のソースの中からデジタルデータを操作することから恩恵を受けるどんなアプリケーションにも等しく適用可能である。
同一チャネル干渉受信機
図1を参照すると、マルチユーザ検出(MUD)受信機12と、同じ名目(nominal)データレートを有する、デジタル位相偏移変調(phase shift key)(PSK)または直角位相振幅変調(QAM)信号などのすべて同種の変調方式を使用して、同一チャネル干渉デジタル信号をすべて同じ周波数で同時に送信している複数のユーザ送信機11から11を備える通信システム10のシステムブロック図が示されている。
ユーザ送信機11から11のそれぞれは一意の既知のトレーニングシーケンスを有する。該トレーニングシーケンスは、受信機アンテナ13で受信されるように大まかに位置合わせされ、その結果トレーニングシーケンスは大部分重複する。この種の同期は、通常、MUD受信機12と同じ場所に配置される装置から送信される同期信号を使用することにより通信システム内で実現される。シンボル遷移の位置合わせは必要とされない。
この特定の実施形態では、MUD受信機12、アンテナ13、信号サンプラ14、ダウンコンバータ16を含むアナログフロントエンド、及びダウンコンバータ16の出力(例えばベースバンド信号または受信信号の他のより低い周波数バージョン)がマルチユーザ検出器18と、ユーザごとにシグナチャ波形を推定するパラメータ推定器20にフィードされる。
ユーザ送信機11〜11からのK個の信号はアンテナ13により送信機11〜11からの信号の合計として受信される。アンテナ13は、信号サンプラ14への信号接続を備える単一偏向アンテナである。この接続は、例えば伝送線路または1のアンテナ13から1の信号サンプラ14に接続する導波管22によって行われる。
信号サンプラ14は、アナログ/デジタル変換器(A/D)により実現されてもよい。信号サンプラ14の出力は、アンテナ13から受信されたサンプリング波形(R)のスナップショット15であり、このスナップショット15は少なくともデータの2個のフレーム内のサンプル数から構成される。代わりに、スナップショット15はデータの複数のフレーム内のサンプル数から構成されてもよい。スナップショット15は、通常、受信信号Rの周波数をベースバンドに変換するためにデジタル無線で使用されるダウンコンバータ16にフィードされる。ダウンコンバータ16の出力17は、同じ周波数と時空間のK個すべての同一チャネル干渉信号からの情報を含む複合ベースバンド信号、r(n,1)である。
ベースバンド信号r(n,1)は、パラメータ推定器20に送信される。マルチユーザ検出器18は、パラメータ推定器20により提供される情報を使用して、同一チャネル干渉デジタル信号をいっしょに復調する。パラメータ推定器20はユーザ送信機11から11により送信され、該K個の信号のシグナチャ波形について解決するためにコンポジット受信信号r(n,1)内に含まれる一意のトレーニングシンボルの(メモリ19に記憶される)知識を使用する。用語「シグナチャ波形」はここで、信号が通過するチャネルのインパルス応答を示すために使用される。用語「チャネル」はここで、伝搬チャネルとアンテナの影響だけではなく、送信機11から11と受信機12フロントエンドで使用されるどんなフィルタリングも含むために使用される。さらに直接シーケンススペクトル拡散システムでは、それは拡散符号も含むであろう。
最適マルチユーザ検出器18は、受信信号と、それらの一意のシグナチャ応答により変換される各ユーザの送信データシンボルの考えられるすべての組み合わせの間の平均平方誤差を最小限に抑える検出器である。この最適マルチユーザ検出器18は、以下のとおりに数学的に表現することが可能であり、
Figure 0004490265
ここでΩ=送信されたデータシンボルの考えられるすべての組み合わせの制約セットである。多くの低複雑性MUDアルゴリズムがここに説明される。
パラメータ推定器20の目的はマルチユーザ検出器18に式1を解くために必要とされる情報を与えることである。各ユーザ及び各ダイバシティポートに一意であるシグナチャ波形30は、各ユーザの送信済みのシンボルの変換を、それらが送信機11〜11から受信機12に伝搬するに従って説明する。これは、送信機11から11でのパルス形状フィルタリングと、受信機12での受信機フィルタリングを含む。いくつかのマルチユーザ検出器は同期のためにデータの各フレーム内のトレーニングシーケンスのロケーションについての情報も必要とし、それらはユーザごとに送信されたシンボルのさらに優れた推定を行うために受信信号内での雑音電力についての情報もまた必要としてもよい。パラメータ推定器20はこのようなパラメータを計算するように構成されてもよく、したがってこれらの入力を必要とするどんなマルチユーザ検出器18とも動作する。
一実施形態では、パラメータ推定器20はスナップショットあたり1度発生する出力を生成し、そのスナップショットの中のデータのフレームごとにパラメータ推定値を含む。これらのパラメータ推定値は、各ダイバシティポート(p)、フレーム(m)、及びアクティブユーザ(k)について、推定されたシグナチャ波形30、
Figure 0004490265
を含む。出力はまた、雑音の平均電力を表すスカラーである推定雑音電力26
Figure 0004490265
及びスナップショット15の各フレームのトレーニングシーケンスのロケーションに対するポインタであるトレーニングシーケンス指数28τTSを含む。出力はまた各ユーザkの状態を含むアクティブユーザベクタ29(u(k))も含む。状態は、ユーザが「アクティブに送信している」または「送信していない」ことを指す。
パラメータ推定器20の出力は、r(n,1)ベースバンド信号17も受信し、送信機11〜11により送信されるK個の同一チャネル干渉信号のそれぞれに対応する信号1、信号2、信号Kまでについて、送信機1シンボル39から送信機Kシンボル38の別々のストリームを生成するMUD18に送信される。システムは、さらに、周波数不一致補償モジュール、ブロックエラー復号モジュール、逆多重化モジュールまたはでデパケッタイジング(depacketizing)モジュール、及びルーティングモジュールなどのMUD18の出力を受け取るように適応された追加のMUD後(post−MUD)処理構成部品(不図示)も含んでもよい。
偏向ダイバシティMUD
図2は、本発明の別の実施形態に従って二重偏向アンテナを有するマルチユーザ検出受信機を描いている。二重偏向アンテナ40を含むと、さらに優れたシンボル決定を行い、それによりシンボル誤り率を削減するためにマルチユーザ検出器18にさらに多くの情報が与えられる。この余分な情報は、直交偏向されたアンテナポートにより受信される信号が実質的には異なるチャネルを通って移動するという事実から引き出される。パラメータ推定器20は、マルチユーザ検出器18に、両方のアンテナポート(p=1,2)が受信するユーザ「k」ごとのシグナチャ波形30を提供する。したがって、マルチユーザ検出器18に渡されるKx2のシグナチャ波形30があるであろう。したがって、MUD18はただ1つの偏向と比較して2倍多くの式を有するが、解くための同数のシンボルを有するであろう。
二重偏向アンテナの使用は、例えば以下の2つのケースで有利である。つまり、第1に、信号が二重直交偏向で送信されるケース、第2に、ただ1つの偏向だけが送信されたとしても、電磁散乱によりかなりの相互偏向されたエネルギーが受信アンテナ40で受信されるケースである。
空間と偏向ダイバシティMUD
図3は、本発明の別の実施形態に従って各アンテナが1つまたは2つの偏向を有する複数のアンテナを有するマルチユーザ検出受信機を描いている。余分なアンテナポートを含むと、MUD18になおさらに多くの情報が与えられ、さらによいシンボル決定を下し、それによりシンボル誤り率を削減することを可能にする。追加の利点を提供するために、余分なアンテナは空間ダイバースでなければならない。言い換えると、アンテナは大幅に異なる伝搬チャネルを提供するほど十分に離間されなければならない。このケースでは、パラメータ推定器20が、p=1,2...Pである信号r(n,p)を処理し、MUD18に渡すためのKxPのシグナチャ波形30を生成する。
フレーム構造及びトレーニングシーケンススライド式検索ウィンドウ
図4は、各フレームセグメントの同じロケーションで受信された複合トレーニングシーケンスβ(n,p)を有するフレーム(framed)セグメントf(n,p)のシーケンスを示す、指定されたダイバシティポート信号「p」について複数(K個)の同一チャネル干渉信号のケースについて受信されたベースバンド信号r(n,p)の下にあるフレーム構造の図であり、トレーニングシーケンスロケータで使用するためのトレーニングシーケンススライド式検索ウィンドウ(l(τ,p))もまた示す。
一実施形態では、受信された複合トレーニングシーケンス、β(n,p)は、付加的なガウス白色雑音、w(n、p)を加えたそのそれぞれのシグナチャ波形S(n,p)で(アスタリスクによって示される)畳み込まれる(convolved)各ユーザのトレーニングシーケンスの合計b(n)の複雑なベースバンドバージョンとして定義される。この関係は数学的には以下のとおりに定義される。
Figure 0004490265
図4は、トレーニングシーケンスロケータ内で使用するためのトレーニングシーケンススライド式検索ウィンドウ1(τ、p)もまた示す。特に、トレーニングシーケンススライド式検索ウィンドウは、トレーニングシーケンスロケータ56(図6)の一部である検出統計計算機90により使用される。これらのスライド式検索ウィンドウはLサンプル長であり、ここでLは受信された複合トレーニングシーケンスβ(n,p)でのサンプル数である。各スライド式検索ウィンドウの第1の指数はF個のサンプルによって分離され、ここでFはデータのフレーム内のサンプル数である。各スライド式検索ウィンドウ1m(τ,p)は受信データ、f(n,p)の対応するフレームを横切って移動し、F個のサンプルシフトの合計について一度に1サンプルづつシフトされる。サンプルシフトτごとに、各スライド式検索ウィンドウ1(τ、p)内のデータは、検出統計の対応する値を計算するために検出統計計算機90により使用される。
パラメータ推定
図5は、本発明の一実施形態に従って構成されるパラメータ推定器のブロック図である。パラメータ推定器20は、トレーニングシーケンスロケータ56と、雑音推定器52と、信号推定器ループ57と、初期変換行列ビルダ63とを含むソフトウェアモジュールを備えて図示されている。該信号推定器ループ57は、さらにシグナチャ波形推定器58と、アクティブユーザテスタ60と、変換行列セレクタ61と、変換行列リビルダ62とを含む。代替実施形態及び構成は類似する機能性を実現してよい。
トレーニングシーケンスロケータ56は、複合受信トレーニングシーケンスβ(n,p)の受信データf(n,p)の各フレーム内でロケーション指数τTSを推定するために使用され、雑音推定器52はダイバシティポートpごとに受信信号r(n,p)内の平均雑音電力
Figure 0004490265
の推定値を計算するために使用される。シグナチャ波形推定器58は、各ユーザKに一意の特徴的なシグナチャ波形
Figure 0004490265
及び受信されたスナップショット内のフレームmごとに各ダイバシティポートpを推定するために使用される。
シグナチャ波形推定器58の出力は、どのユーザ信号が指定されたスナップショットに存在するのかを検出するアクティブユーザテスタ60にフィードされ、シグナチャ波形推定器58で使用される変換行列
Figure 0004490265
を再構築する変換行列リビルダ62に出力を提供する。この行列はアクティブユーザテスタ60により計算されるようにアクティブユーザのトレーニングシーケンス、b(n)だけを使用することにより再構築される。
変換行列リビルダ62の出力は,シグナチャ波形推定器58に送信するために初期変換行列ビルダ63から出力
Figure 0004490265
または変換行列リビルダ62から出力
Figure 0004490265
を選択する変換行列セレクタ61にフィードされる。この特定の実施形態では、変換行列セレクタ61は、指定されたスナップショット内でのシグナチャ波形の初期推定値につねに
Figure 0004490265
を選択し、データの同じスナップショットのシグナチャ推定値のすべての以後の再計算に
Figure 0004490265
を選択する。これにより、シグナチャ波形推定器58は、アクティブユーザテスタ60により決定されるようにアクティブなユーザのためだけに特徴的なシグナチャ波形
Figure 0004490265
のさらによい推定値を計算できるようになる。
シグナチャ波形推定器58を動作する、アクティブユーザテスタ60を実行する、及び変換行列リビルダ62を動作するこのプロセスはシグナチャ推定ループ57と呼ばれる。シグナチャ推定ループ57は、過去の反復で計算されたアクティブユーザテスタ60の出力が現在の反復でのアクティブユーザテスタ60の出力に等しくなるまで繰り返すことができる。パラメータ推定20構成要素内のシグナチャ推定器ループ57の最大数を設定することもまた可能である。
このループ57を通して反復のたびに、シグナチャ波形推定器58の出力でのシグナチャ波形の数が過去の反復で計算されたアクティブユーザの数に等しいことに注意せよ。さらに、第1の反復で、シグナチャ推定の数が考えられるユーザの総数Kに等しいことに注意せよ。いったんアクティブユーザの最終的なシグナチャ推定値が計算されると、結果として生じる波形は、どのユーザが現在のスナップショットでアクティブであるのかを報告するユーザ状態ベクタ、u(k)とともにパラメータ推定器20の出力として渡される。
初期変換行列ビルダ63は、ユーザごとに、例えばマルチユーザ検出受信機12のメモリ19に事前に記憶されてよい、既知のトレーニングシーケンスデータ、b(n)を受け取る。各ユーザのトレーニングシーケンスデータは、変換行列セレクタ61にフィードされる初期変換行列
Figure 0004490265
を構築するために使用される。
雑音推定器52は、すべてのp=1,2,...P、ダイバシティポートについて入信信号、r(n,p)の雑音電力を推定し、アクティブユーザテスタ60とMUD18に情報をフィードする。この推定は、通常、スナップショットあたり1回行われ、該スナップショットは少なくとも、2つのフレームの中のサンプル数であるが、雑音電力がゆっくりと変化する、あるいはまったく変化しない場合には頻繁に行われる必要はない。雑音推定器52の精度は、複合トレーニングシーケンス推定値の数、
Figure 0004490265
が増加するにつれ改善する。複合トレーニングシーケンス推定値の数を増加するためには、受信された合成ベースバンド信号r(n,p)の中のフレームの数、f(n,p)が増加しなければならない。この結果スナップショットサイズが拡大する。代わりに、トレーンニングシーケンスセレクタ56は、受信されたデータの複数のスナップショットについて複合トレーニングシーケンス推定値、
Figure 0004490265
を記憶し、記憶された複合トレーニングシーケンス推定値
Figure 0004490265
の合計数を使用して推定される雑音電力を計算しなければならない。
トレーニングシーケンスロケータ56は、受信されたスナップショットベクタ、r(n,p)の各フレーム、f(n,p)内のトレーニングシーケンスの位置を決定し、サンプル指数、トレーニングシーケンスロケーション指数28と呼ばれるτTSの形式でこの情報をMUD18にフィードする。加えて、受信されたスナップショットの中のトレーニングシーケンスの位置は雑音推定器52に、及びそれぞれ平均雑音電力推定値
Figure 0004490265
及びシグナチャ推定値
Figure 0004490265
を決定するためにr(n,p)内の各フレーム、f(n,p)のどのセクションを処理するのかを決定することが使用されるシグナチャ波形推定器58にフィードされる。シグナチャ波形推定器58は、ダイバシティポートpごとに、コンポジット受信入力信号、r(n,p)内のそれぞれK個の個々の同一チャネル干渉信号の各フレームmのシグナチャ波形
Figure 0004490265
を推定し、アクティブユーザテスタ60とMUD18にこの情報を出力する。
図6は、本発明の一実施形態に従って構成されるトレーニングシーケンスロケータ構成要素のフローチャートである。トレーニングシーケンスロケータ56構成要素は、受信データの2つのフレームの最小値を処理する、各ダイバシティポート(p=1,2...P)から受信される、受信合成ベースバンド信号、r(n,p)内でトレーニングシーケンスの位置を見つけるように構成される。データ(F)の送信フレーム内のサンプル数は受信機に既知であるため、各受信フレーム、r(n、p)内のf(n,p)はデータの送信フレームと同じタイムレートで反復する。したがって、複合トレーニングシーケンス、β(n,p)のロケーションは、受信されたデータの各フレーム、f(n,p)について同じとなる。しかしながら、データの送信フレームと受信フレームは当初同期していないため、データの送信フレームの始まりは受信されたデータのフレーム、f(n,p)の始まりで開始しない可能性がある。これは、受信されたデータの各フレームの中のトレーニングシーケンスのロケーションもまた各フレームの始まりで開始しない可能性があることを意味する。
トレーニングシーケンスロケータ56は、受信されたデータの各フレームの中のトレーニングシーケンスのロケーションを見つける。これを行うために、Lサンプル長(受信された複合トレーニングシーケンスと同じ長さ)であるスライド式検索ウィンドウベクタ、l(τ,p)が受信データの各フレームを通して同時に適用され、ウィンドウ表示のフレームの各組み合わせの間の相互関係が計算され、次に検出統計計算機90で平均化される。結果は、まさに受信データのフレームの長さ(Fサンプル長)である検出統計、d(τ)である。ペイロードデータはフレームごとに無相関であるため、スライド式検索ウィンドウが各フレームのペイロードデータの上にあるときには検出統計は非常に低い値を有する。他方、複合トレーニングシーケンス、β(n,p)はフレームごとにきわめて相関性がある。したがって、スライド式検索ウィンドウが各フレーム内の複合トレーニングシーケンス上にあるとき、検出統計は非常に高くなる。したがって、検出統計d(τ)のピークのロケーションτは、各フレームシーケンスf(n,p)の中のトレーニングシーケンスのロケーションとなろう。
トレーニングシーケンスロケータ56構成要素に対する入力は、各ダイバシティポート、(p=1,2,...P)からの合成ベースバンド受信信号r(n,p)である。トレーニングシーケンスロケーション指数τの推定値は、検出統計計算機90によりダイバシティポート信号r(n,p)ごとに別個に計算される。トレーニングシーケンスロケーション指数τを推定する上での第1のステップは、その受信信号を使用して検出統計d(τ)を計算するために検出統計計算機90に受信信号r(n,p)を提供することである。前述されたように、この検出統計の各要素は、指定されたトレーニングシーケンスサンプル指数τのためのスライド式検索ウィンドウのそれぞれの組み合わせについて、相互関係係数
Figure 0004490265
を計算することにより生成される。いったん相互関係係数の各組み合わせが計算されると、それらは平均され、τという特定された値について、検出統計d(τ)の値として出力される。本発明の一実施形態に従ってこのプロセスを実行するために必要とされる段階的な計算は以下のとおりである。
ステップ1.指定された検索ウィンドウサンプル指数τについて、指定された信号r(n,p)の中の受信データの各フレームのためのスライド式検索ウィンドウ、l(τ,p)を定める。
Figure 0004490265
ステップ2.各スライド式検索ウィンドウ、l(τ,p)のエネルギーe(τ,p)を計算する。
Figure 0004490265
ステップ3.スライド式検索ウィンドウのそれぞれの組み合わせについて、相互関係係数
Figure 0004490265
を計算する。
Figure 0004490265
ステップ4.対応する相互関係係数を平均化することにより、ダイバシティポートpのための指定された検索ウィンドウサンプル指数、τについて検出統計d(τ)を計算する。
Figure 0004490265
このプロセス(ステップ1から4)は、検索ウィンドウサンプル指数{τ=1,2...F}ごとに、及びダイバシティポート{p=1,2...P}ごとに繰り返される。
図6を参照すると、いったんダイバシティポート受信信号r(n,p)ごとの検出統計d(τ)が計算されるたびに、それは、ダイバシティポート信号ごとにトレーニングシーケンスの推定されたロケーションτが、検出統計を最大にするスライド式検索ウィンドウ指数τを検出することにより計算されるトレーニング指数ファインダ92にフィードされる。トレーニング指数ファインダ92は、検出統計d(τ)ごとに次の式を計算する。
Figure 0004490265
次に、信頼測定基準計算機94が、各検出統計から信頼測定基準を計算する。これは、各検出統計のピーク対rms値を計算することにより行われる。このプロセスは、例えば検出統計d(τ)ごとに次の計算を実行することにより実現できる。
Figure 0004490265
前述されたように、この検出プロセス全体が各ダイバシティポートpの受信信号r(n,p)に別個に適用される。いったん信号ごとにトレーニングシーケンスロケーション、τが推定され、信頼測定基準cが計算されると、どの推定値を使用するのかを決定するために決定試験が適用される。比較器(comparator)96は、各信頼測定基準の値を比較し、最高の信頼測定基準cを有する推定されたトレーニングシーケンスτに等しい出力トレーニングシーケンスロケーションτTSを設定することによりこの決定試験を実行する。このプロセスは、以下の式の中にて数学的に説明される。
Figure 0004490265
図7は、本発明の一実施形態に従って構成される雑音推定器構成要素のブロック図である。パラメータ推定器20の雑音推定器52は、各ダイバシティポートpから受信信号r(n,p)内の平均雑音電力の正確な推定値を計算する。通常、おもな雑音ソースは、信号サンプラ14の中の低雑音増幅器(LNAs)の第1段により生成される熱雑音である。この雑音は合成ゼロ平均付加ガウス白色雑音(AWGN)として正確にモデル化できる。
各ダイバシティポートpはLNAの独自の段を含むため、受信信号r(n,p)の雑音電力はポートごとにわずかに異なる場合がある。したがって、ダイバシティポートpごとに雑音電力の別々の推定値を取得することが重要である。雑音は付加的であるため、複合トレーニングシーケンス、β(n,p)を含む受信データ、r(n,p)の各フレームのセクションは付加的な白色雑音、w(n,p)を加えられた複合トレーニングシーケンスとしてモデル化できる。また、この雑音は、受信されたデータとは統計的に無関係であると考えられる。前記に説明されたように、トレーニングシーケンスロケータ56構成要素は、複合トレーニングシーケンスを含む受信データの各フレームf(n,p)の中のロケーション指数を決定する。これにより、雑音推定器52は、受信された複合トレーニングシーケンス
Figure 0004490265
の推定値を含む、各フレームのセクションを抽出できる。
トレーニングシーケンスセレクタ70は、受信データのフレームごとの受信トレーニングシーケンスに対応するその信号のセクションを出力するルーチンの中にダイバシティポートpに対する合成ベースバンド信号r(n,p)を受け入れる。トレーニングシーケンスセレクタ70はまた、トレーニングシーケンス指数τTS28と、このような出力を選択する上で使用するための受信信号r(n,p)内で捕捉されたデータのフレーム数Mを受信する。Mは、受信信号r(n,p)のサイズとフレームあたりの既知のサンプル数、Fに基づいて計算される。
トレーニングシーケンスセレクタ70は、(m=1,2,...M)について、複合トレーニングシーケンス
Figure 0004490265
の優れた表現を獲得するために、ベクタ
Figure 0004490265
のそれぞれの平均値を計算するアベレージャ72にフィードされる、
Figure 0004490265
を出力する。この特定の実施形態では、アベレージャ(averager)72は、加算ルーチン74と1/M乗算器(multiplier)76を備える。
推定された信号
Figure 0004490265
は、各ベクタ
Figure 0004490265
の中に含まれる、推定される雑音信号、
Figure 0004490265
を取得するために加算器(summer)78、78、...78内の各ベクタ
Figure 0004490265
から差し引かれる。次に、各雑音信号、
Figure 0004490265
が、各フレームMの中の平均雑音電力
Figure 0004490265
の推定値を取得するために各雑音信号の分散が計算される分散計算機80,80,...80にフィードされる。各フレーム内の平均雑音電力
Figure 0004490265
は、ダイバシティポートpについて受信された信号r(n,p)内の平均雑音電力
Figure 0004490265
の推定値を決定するアベレージャ82にフィードされる。複数のダイバシティポートがある場合には、各ダイバシティポート(p=1,2...P)によって受信される受信信号r(n,p)の雑音電力推定値を取得するために前記プロセスは繰り返される。
受信された複合トレーニングシーケンス
Figure 0004490265
の正確な推定値は、各ダイバシティポートpから受信信号r(n,p)の雑音電力推定値を別々に取得するために使用される。ダイバシティポートpの受信された信号、r(n,p)の中の雑音電力を推定する際の第1のステップは、トレーニングシーケンスセレクタ70が、トレーニングシーケンスロケーション指数(τTS)、フレーム(F)あたりのサンプル数、処理するための各受信トレーニングシーケンス内のサンプル数(N)、及び使用するための各受信トレーニングシーケンス(δ)へのオフセットに基づいて、受信データの各フレームの中の推定された複合トレーニングシーケンス
Figure 0004490265
を抽出することである。Nとδの両方が,トレーニングシーケンスセレクタ70のメモリに記憶されるパラメータであり、したがって受信された複合トレーニングシーケンスの任意のセクションを選択するために修正できることに注意せよ。これらの値は、通常、各複合受信トレーニングシーケンスの中心部分が受信信号r(n,p)から抽出されるように設定されるであろう。該トレーニングシーケンスセレクタ70は、数学的に以下のとおりに説明される。
Figure 0004490265
いったん受信データのフレームmごとの推定された受信トレーニングシーケンス
Figure 0004490265
が受信信号r(n,p)から抽出されると、それらはすべて、受信されたトレーニングシーケンスの各推定値に追加される雑音ベクタ、w(n,p)の影響を最小限にするためにアベレージャ72によって互いと平均化される。これにより、実際の受信されたトレーニングシーケンス
Figure 0004490265
のより正確な推定値が生じる。このプロセスは、数学的に以下のとおりに説明される。
Figure 0004490265
いったんこれが行われると、推定された受信トレーニングシーケンス
Figure 0004490265
はそれぞれに含まれる雑音信号
Figure 0004490265
の推定値を取得するために各ベクタ
Figure 0004490265
から差し引かれる。
次に、各フレームの中の平均雑音電力の推定値
Figure 0004490265
を取得するために、各雑音信号の分散が分散計算機80から80によって計算される。この計算は以下のとおりに表現される。
Figure 0004490265
次に、これらの雑音電力推定値
Figure 0004490265
のそれぞれは受信信号r(n,p)の中の平均雑音電力の推定値
Figure 0004490265
を取得するためにアベレージャ82の中で平均化される。この平均化は、数学的には以下のとおりに実行される。
Figure 0004490265
各受信信号r(n,p)の雑音電力推定値を取得するために、このプロセス全体がダイバシティポートpごとに繰り返される。
図8は、本発明の一実施形態に従って構成されるシグナチャ波形推定器構成要素のフローチャートである。シグナチャ波形推定器58構成要素は、受信データの各フレームmの受信トレーニングシーケンス
Figure 0004490265
を行列乗算器(Matrix Multiplier)66に提供するトレーニングシーケンスセレクタ64を含む。シグナチャ波形推定器58は、各ダイバシティポートp及び受信データの各フレームmについて各ユーザの伝送信号を、それが送信機から受信機に伝搬するにつれて変換する特徴的なシグナチャ波形
Figure 0004490265
を推定する。この変換は、数学的には以下のとおりに表現される。
Figure 0004490265
この式は、ダイバシティポート(p)から受信された合成ベースバンド信号が、ユーザk、及び付加的な白色雑音、w(n,p)が追加されるダイバシティポートpにとって一意であるサンプルあたりT秒でサンプリングされた対応する特徴的なシグナチャ波形、s(nT,p)で畳み込まれる(convolved)(アスタリスクにより示される)各ユーザの伝送信号、d(n)の合計であることを示している。この式は以下のように行列形式で書き直すことができ、
Figure 0004490265
ここで
Figure 0004490265
Figure 0004490265
Figure 0004490265
である。
指定されたダイバシティポートpの場合、これらのシグナチャ応答を推定するために使用される手法は、複合トレーニングシーケンスβ(n,p)を含む受信信号の部分を、各ユーザk=1,2,...Kにより送信される実際の既知のトレーニングシーケンス、b(n)と比較することである。これは、各ユーザから送信されたトレーニングシーケンスb(n)が受信機12により知られているため、及び受信された複合トレーニングシーケンス、β(n,p)の推定値がトレーニングシーケンスロケータ56を使用して受信された信号、r(n,p)から抽出できるために達成できる。このケースでは、
Figure 0004490265
である。
各ユーザk、及び各ダイバシティポートpの特徴的なシグナチャ波形
Figure 0004490265
の最大尤度推定値は、受信された複合トレーニングシーケンス、β(n,p)とその対応するシグナチャ波形
Figure 0004490265
で畳み込まれる各ユーザのトレーニングシーケンスの合計の間の平方誤差を集合的に最小限に抑えるものである。最大尤度推定値は、数学的に以下のとおりに行列形式で表され、
Figure 0004490265
ここでΩはの考えられるすべての組み合わせの集合に等しい。
最大尤度手法を使用すると、シグナチャ推定値はゼロ強制(zero forcing)基準を使用することにより解決できる。これは、以下のとおりにゼロに等しい最小の内側で式を設定することにより行われる。
Figure 0004490265
いったんこれが行われると、特徴的なシグナチャ波形ベクタは、以下のとおりにsML(p)の一次(linear)方程式の前記集合を解くことにより計算できる。
Figure 0004490265
前記最大尤度式の解に基づき、第1のステップは、シグナチャ波形推定器58内のトレーニングシーケンスセレクタ64が、フレームmごとに、そのフレームの中の受信された複合トレーニングシーケンス
Figure 0004490265
を含む受信信号r(n,p)の部分を抽出することである。これは、トレーニングシーケンスセレクタ(ステップ64)で行われ、トレーニングシーケンスτTSのロケーション、フレーム、Fあたりのサンプル数、選択するための受信トレーニングシーケンスのサンプル数、Nβ、及び使用するための各受信トレーニングシーケンスδβへのオフセットに基づいている。
βが(N+N−1)に等しく、Nはシグナチャ推定値ごとに使用するためのサンプル数であり、Nは既知のトレーニングシーケンスの中のサンプル数であることに注意せよ。さらに、Nとδβがこのトレーニングシーケンスセレクタ64のメモリに記憶されるパラメータであり、したがって各フレーム内の受信複合トレーニングシーケンスの任意の部分を選択するために修正できることに注意せよ。これらの値は、通常、複合受信トレーニングシーケンス全体が受信信号r(n,p)から抽出されるように設定されるであろう。トレーニングシーケンスセレクタステップ64は、以下のとおりに数学的に説明される。
Figure 0004490265
シグナチャ波形を推定する際の次のステップ66から66は、以下のとおりに行列乗算器ステップ66から66を使用して、変換行列セレクタ61から受信される変換行列
Figure 0004490265
を、m=1,2,...Mのフレームmごとに複合受信トレーニングシーケンス推定値、
Figure 0004490265
を含む受信された合成ベースバンド信号の部分で乗算することである。
Figure 0004490265
Figure 0004490265
ここで
Figure 0004490265
Aがアクティブなユーザの総数に等しく、kがa番目のアクティブなユーザの指数に等しいことに注意せよ。したがって、最初の反復の場合、K人すべてのユーザがアクティブであると仮定されるため、シグナチャ推定ループを通る最初の反復の場合、k=k、k−1,2...Kの場合はA=Kである。
シグナチャ波形推定値の初期計算時、変換行列、
Figure 0004490265
は、変換行列セレクタ61ルーチンを経由して初期変換行列ビルダ63からシグナチャ波形推定器58に渡される。指定された合成ベースバンド受信信号r(n,p)のためのシグナチャ波形のすべての以後の推定値では、変換行列
Figure 0004490265
が、変換行列セレクタルーチン61を経由して変換行列リビルダ62からシグナチャ波形推定器58に渡される。これは、アクティブなユーザ(k)のシグナチャ推定値だけが計算されるように行われる。変換行列
Figure 0004490265
の次元が、それぞれの特徴的なシグナチャ推定値、
Figure 0004490265
の中のサンプル数(Ns)、現在アクティブなユーザの推定数(A)、及びそれぞれの受信複合トレーニングシーケンス推定値
Figure 0004490265
のサンプル数(Nβ)の関数であることに注意せよ。
図9は、本発明の一実施形態に従って構成されるアクティブユーザテスタ構成要素のフローチャートである。アクティブユーザテストの目的は、受信されたデータ、r(n,p)の各スナップショットの情報をどのユーザが送信しているのかを判断することである。この試験は、シグナチャ推定ループの反復のたびに一度実行される。これは、各ダイバシティポートpにわたって、受信されたデータの各フレームmの中のユーザkごとのシグナチャ推定値
Figure 0004490265
に基づいてユーザ信号の中の平均受信電力を計算することにより行われる。これらの電力推定値
Figure 0004490265
は、各フレームmにわたり平均され、結果、
Figure 0004490265
が推定された雑音フロア
Figure 0004490265
に関して設定される所定の閾値rthに比較される。
アクティブユーザテストを実行するために、シグナチャ推定値
Figure 0004490265
は、ダイバシティポート「p」用のアクティブテストと呼ばれるステップ100pによって別々にダイバシティポートpごとに処理される。このステップの出力は、長さ(K)要素である1とゼロのシーケンス、u(k)である。u(k)=0の場合には、ユーザkは「非アクティブ」であると推定され、u(k)=1の場合には、ユーザkはダイバシティポートpのためのシグナチャ推定値に基づいて「アクティブに送信している」と推定される。この出力シーケンスu(k)は、ダイバシティポートpのアクティブユーザテスト結果シーケンスと呼ばれる。
いったん各ダイバシティポートについてこのテスト結果シーケンスが計算されると、該シーケンスは論理「OR」演算子102に渡される。この論理「OR」関数はPu(k)シーケンスのどれかが、k={1,2,...K}である場合、kの各値について1に等しい場合に、結合されたアクティブユーザテスト結果シーケンス、u(k)を1に等しく設定する。したがって、結合されたアクティブユーザテスト結果シーケンスu(k)は、u(k)シーケンスのどれかが1に等しい場合に1に等しくなり、それ以外の場合、各ユーザk=1,2...Kの場合u(k)は0に等しくなる。
図10は、本発明の一実施形態に従ってダイバシティポート「p」用のアクティブユーザテストのフローチャートである。指定されたダイバシティポートpについて、アクティブユーザテスト結果、u(k)が、以下のプロセスを実行することにより、シグナチャ推定ループの過去の反復でアクティブであるために計算された指定されたユーザkについて計算される。
第1のステップ104から104は、図10に示されるようにk=k,k...kである場合にフレームmごと及びユーザkごとに推定された特徴的なシグナチャ応答、
Figure 0004490265
を使用してユーザ(k)ごとに平均受信信号電力を推定することであり、それは以下のように計算され、
Figure 0004490265
ここでFsymはシンボルあたりのサンプル数に等しい。
いったん信号電力
Figure 0004490265
がユーザkごとに推定されると、各フレームからの結果は以下のとおりに平均される。
Figure 0004490265
次のステップ106から106では、ユーザ、(k=k,k,...k)ごとのこれらの推定されたシグナチャ電力が、推定された雑音フロア、
Figure 0004490265
に関して検出閾値、rthに比較される。ユーザkのための推定されたシグナチャ電力、
Figure 0004490265
が、推定雑音フロア
Figure 0004490265
で相対閾値rthの積以上である場合には、そのユーザkについてアクティブユーザ試験結果、u(k)は1に設定される。それ以外の場合、それは0に設定される。このテストは数学的に以下のとおりに表される。
Figure 0004490265
結合結果ステップ108では、u(k)の結果のすべてが以下のとおりにオリジナルユーザ状態ベクタu(k)と結合される。
Figure 0004490265
図5を参照すると、変換行列リビルダ62は、各ユーザk、ダイバシティポートp及びフレームmについて、シグナチャ応答を推定するためにシグナチャ波形推定器58により使用される変換行列
Figure 0004490265
を再構築するためにアクティブユーザテスタ60の結果を使用する。この変換行列は、トレーニングシーケンス行列、
Figure 0004490265
内の、非アクティブユーザに対応する部分行列、
Figure 0004490265
を削除することにより再構築される。これは、アクティブユーザkaのシグナチャ応答を推定するためだけに使用される更新された変換行列
Figure 0004490265
を作成するために実行される。このプロセスは、数学的に以下のとおりに表現され、
Figure 0004490265
Figure 0004490265
ここでkは、以下のアルゴリズムを使用して定義できる。
Figure 0004490265
更新された変換行列
Figure 0004490265
は、変換行列セレクタ61を経由してシグナチャ波形推定器58に渡される。シグナチャ波形推定器58の内側では、更新された変換行列
Figure 0004490265
が、アクティブユーザのためだけにさらに正確なシグナチャ波形推定値を計算するために、各ダイバシティポートp及び各フレームmについて推定された受信トレーニングシーケンス
Figure 0004490265
に適用し直される。
初期変換行列ビルダ63は、各考えられるユーザ(k=1,2,...K)について長さNbサンプルである既知のトレーニングシーケンス(b(n))に基づき初期の変換行列
Figure 0004490265
を計算する。この変換行列は以下のとおりに計算される。
最初に、既知のトレーニングシーケンス合成積行列(convolution matrix)
Figure 0004490265
が求められる。
Figure 0004490265
Figure 0004490265
第2に、変換行列
Figure 0004490265
は以下のように計算される。
Figure 0004490265
初期変換行列
Figure 0004490265
は、変換行列セレクタ61を経由してシグナチャ波形推定器58に渡され、各ダイバシティポートpにわたり考えられる各ユーザk及び受信データの各フレームmについて初期シグナチャ波形推定値
Figure 0004490265
を計算するために使用される。また、既知のトレーニングシーケンス合成積行列は、部分行列
Figure 0004490265
が、シグナチャ推定ループの反復のたびに再生される必要がないように、変換行列リビルダ62に渡される。
変換行列セレクタ61構成要素は、どの変換行列がシグナチャ波形推定器58に渡されるのかを選択するために使用される。合成ベースバンド受信信号r(n,p)の指定されたスナップショットについて、初期シグナチャ波形推定値
Figure 0004490265
は、初期変換行列
Figure 0004490265
を使用して計算される。したがって、このケースでは、変換行列セレクタは、その出力
Figure 0004490265

Figure 0004490265
に等しく設定することにより、シグナチャ波形推定器58に
Figure 0004490265
を渡す。いったんシグナチャ波形推定器58がユーザごとにシグナチャ波形を推定すると、結果はどのユーザが現在アクティブであるのかを判断するためにアクティブユーザテスタ60に渡される。これらの結果は、次に、アクティブなユーザについて既知のトレーニングシーケンス合成積行列
Figure 0004490265
だけを使用して変換行列を再構築するために変換行列リビルダ62に渡される。したがって、初期シグナチャ波形推定値が計算された後、変換行列セレクタ61は、その出力
Figure 0004490265

Figure 0004490265
に設定することにより、再構築された変換行列
Figure 0004490265
をシグナチャ波形推定器58に渡す。
低複雑性MUDアルゴリズム
過去の情報を備える従来の最小平均二乗誤差(Minimum Mean Squared Error)(MMSE)マルチユーザ検出器が一般にターボ設定で現れる。ユーザの集合全体(通信チャネル)について同時に線形演算子を計算するために、代数的な手段が使用される。これは、時間の特定の瞬間ごとに、過去の情報、または情報の各ユーザのビットのありそうな値に関する知識を活用することにより行われる。しかしながら、このような従来のマルチユーザ検出器は、オーバロードされた通信システム向けに設計されていなかったため、不合理な計算の複雑性だけではなく数値の不安定性にもさらされる。
本発明の原理を使用して、従来の検出器の不安定性の問題は解決される。この解決策は現在利用可能な、あるいは出現しつつある多くの種々のマルチユーザ検出器で役立ち、ターボ検出器のマルチユーザ検出器ブロックで必要とされる計算の数を劇的に削減し、その結果信頼できる演算がリアルタイムインプリメンテーションで達成できる。一実施形態では、MUDモジュールはMMSE検出器内の確率(probability)データの部分的な量子化を利用する。MMSE検出及びハイブリッド取り消しは、ターボ検出器での可変サイズ行列演算だけではなく、同じマルチユーザ検出ブロックでも可能にされる。
部分的に量子化された過去の情報を備える低複雑性MUDの例の実施形態は、図11から図15Cを参照して説明される。事前白色化された(prewhitened)データに基づいた低複雑性MUDの別の例の実施形態は図16から図17に関して説明される。
部分的に量子化された過去の情報を備える低複雑性MUD
図11は、本発明のある実施形態に従ったターボMUD内の入信データストリームの反復処理の図表示である。各ストリームは異なるユーザに対応する。例えば、任意の1つのストリームは同じ周波数上で動作し、同時に送信しており、干渉信号を生じさせる複数のCDMA携帯電話からである場合がある。ストリームは、再び同じ周波数上で、同時にヘッドエンドシステムと通信している多くのケーブルボックスからである場合もある。代りに、ストリームが802.11無線ローカルエリアネットワークの複数のユーザから、あるいは超ワイドバンドネットワークからである場合がある。代りに、ストリームは、読み取られるヘッドが極めて濃密なトラックの1つまたは複数の重複しているデジタル記憶装置ディスク上の隣接するトラックからの干渉信号である場合がある。とにかく、発生する事実は、110でここで描かれている多くのパケットまたはデータの塊が同時に送信され、各パケットが例えば数百ビットを含んでいるということである。
干渉信号を析出するためには、入信信号は時間間隔単位で入信信号を処理し、116で描かれているように、ユーザごとに個々のビットの最も妥当だとされる数字、つまり推定値を提供するターボMUD114に、矢印112により描かれるように適用される。これらの最も妥当だとされる数字または推定値は、次に、特定のビットがそうなると考えられているものとなる尤度を高めるために、フィードバックループ内でMUDに適用される。これは、誤り率のような時間が所定値を下回ることが統計的に既知となるまで反復プロセスで行われる。それが下回ると、デジタル入信ストリームに対応するビットは、入信データがターボMUD114により分離された干渉信号で回復されるように出力される。
図12は、本発明の一実施形態に従って構成されたターボマルチユーザ検出器を描くブロック図である。この実施形態は、例えば、図1から3のマルチユーザ検出器18を実現するために使用できるであろう。ここでは、MUDモジュールからのソフト出力がエラー訂正器を通して、特定の確実度を有するそれらの推定値が確かめられ、フィードバックループ内を通ってMUDに戻されないようにスレッショルドされる(thresholded)最良の推定値平均と結合し直される。むしろ、それらの推定値は、相対的な確実度のそれらのビットを無視することによって実質的に削減されるMDUによる処理を使用してコンバイナ内で活用される。
さらに詳細には、MUD114は、L時間間隔の間に例えば数百ビットがあることを示すためにLが使用されるビットストリーム内の独立したビットを処理するように適応される。該L時間間隔は、波形122により描かれている干渉信号の入信データストリームの一部を描写するために120で描かれている。MUD114は、ここで124によって示されるように一度にこれらの入信信号のビットサイズ部分だけを考え、その瞬間にユーザのそれぞれからビットの最善の推定値の数を生成するためにこれらの信号を処理する。ここでは、第1のユーザの最善の推定値は列124に、第2のユーザは列126で示されている。
「ソフト」MUD出力と呼ばれる実数に関してのこれらの推定値は、各ユーザの先行するビットと続いて起こるビットのウィンドウに基づいて、一実施形態においてビットが1または0のどちらかであるという尤度を判断する129で示されるように、エラー訂正機能を実行するエラー訂正装置128に結合される。したがって、ビット131は、先行するウィンドウ133と続いて起こるウィンドウ135内のビットにより決定される。エラー訂正装置の出力は、ある特定のビットがそれがそうなるといわれているものとなる尤度を高めるためのビットの反復処理を可能にするように、MUD114に適用し直される。
この反復処理は、ユーザのそれぞれからのビットのそれぞれについての膨大な数の可能性のために状況の中の最良な状況の元でさえ法外な時間量を要することである。本発明の一実施形態は、MUDの次の反復の間の考慮から、ほとんど確実であると判断されるそれらのビットを排除する。
例えば、10個のビットの中から、3個が確実であると仮定すると、MUDによる次の反復での処理は、前回の反復の間の処理の12%に削減できる。これを行うためには、スレッショルドモジュール130がMUD114に戻されることから、それらが一実施形態では−1または+1のどちらかを割り当てられるほど十分に確実であると判断されるそれらのビットを排除する。このスレッショルドモジュール130は、線路132上の相対的に不確実なビット推定値を通過し、線路134と136上でそれらの量子化された値、+1または−1として、次の反復でMUD114からの再計算された最も妥当とされる数字とこれらの確実なビットを結合する結合装置138まで相対的に確実なビット推定値を通過する。
結果は、次の反復でエラー訂正モジュール128に提供される情報が該確実なビットと、ビットの残りについて再計算された最も妥当とされる数字も考慮に入れる。エラー訂正装置128はこれらの結合された推定値を処理する。エラー訂正装置128の出力は、スレッショルド装置130を通して再びMUD114にフィードバックされ、その結果、この次の通過で、なおさらに多くのビットが確実と宣言でき、次の反復は残された不確実なビット推定値だけを検討し、再びこれ以後の反復で劇的な処理の優位性をもたらす。エラー訂正装置128の出力が所定の数未満の誤り率を有すると予測されるときには、反復は停止して、出力は分離され復調された信号として線路140で使用できる。
図13は、図12のシステムにより実行される例のスレッショルド(thresholding)の図表示である。この特定の例では、−1と−.95の間にあるビットには−1の値が割り当てられるのに対し、.95と+1の間にあるビットには+1の値が与えられる。他のビット推定値が、矢印132により示されるように通過する一方、該スレッショルド装置30は、例えば、それぞれ−1と+1を戻す、正しいという−1と−.95の確率の間、及び正しいという+.95と+1の間の確率などのそれらのビット推定値を求める。分かるように、スレッショルドは、144と146で示されている閾値のどちらかの側に存在するそれらのビット推定値だけが、コンバイナ138に直接的に結合されるそれらの「確実な」ビットとなるほどであるのに対し、残りのビット推定値はMUD114まで通過する推定値である。
図14は、図12のシステムにより実行される例のスレッショルドと結合を描くブロック図である。分かるように、K掛けるLの行列を考慮すると、本発明のある実施形態に従ったスレッショルドはMUD114にとっての処理優位点を生じさせる。K掛けるLの行列の所定の部分集合150は限られた数のユーザと限られた数のビットを反映する。
図15Aは、確実な値を有する図14の例のK掛けるLの行列内のビットの内の2つの概略図である。分かるように、この部分集合は3人のユーザk、k及びkと、3つの時間間隔l、l及びlを含んでいる。スレッショルド装置130のスレッショルド動作のために、kビット(152と示される)とkビット(154と示される)が確実であると見なされる。MUD114による次の反復では、不確実なままであるビット推定値がMUD114により処理されるのに対し、量子化され、MUD114によってさらに処理されない(つまり未処理)ビット推定値は通過する。
行列からの値のこの処理は、MUD14がフィードバックループの(160として示される)最良のビット推定値の7つを処理することを示す図15Bにさらに明確に示されている。2つの確実な値(162として示される)は通過する。図15Cは、エラー訂正装置130に適切な推定値を与えるために、未処理信号と結合される次の反復で処理されるMUD信号の概略図である。
処理されたビットと未処理のビットは、エラー訂正器128に結合される結果が164で図示されるようであるようにそれ以後138で結合される。理解されるように、スレッショルドの後、エラー訂正器128へ次の反復で通過されたものが、確実なビットとともに、ビットがどうなるべきかの新しい最善の推定値である。これにより、エラー訂正器128は、コンバイナ138からの出力をさらによく処理し、スレッショルド回路130を通してさらに優れた推定値をMUD114に提供できるようになる。結果的に、最終的なビット推定値がより迅速に確かめられ、処理時間を劇的に短縮することになる。したがって、オリジナルの入信干渉信号は、より少ない反復とより少ない計算でさらに迅速に回復される。
量子化ターボMUD
部分的に量子化された過去の情報アルゴリズムを有する低複雑性Mudに根拠を与える数学的な背景がここに説明される。ターボMUDプロセッサのMUD処理ブロックが、それが以下の尤度関数を直接的に最大にするときにその究極の性能に達することは、文献の中で周知である。
Figure 0004490265
この尤度の2つの部分は、候補ビット(第1の項)に対するデータの適合、及びビットが何であったのか(第2の項)に関する過去の知識のための尤度貢献のためである。ターボMUDでは、これらの過去の尤度は未知数(等しくありそうな0または1)として開始し、以後ループが反復するにつれ、最後の反復からの尤度にすぎなくなる。第2の項の尤度貢献が相互に無関係であると仮定することは典型的であり、それは以下のように書き直すことができる。
Figure 0004490265
部分的に量子化された過去の情報を備えた低複雑性Mudでは、式36はその極値に対する第2の尤度貢献を量子化することによりほぼ解かれる。実際に、ビットの過去の尤度が完全に未知である場合には、それが1または0である確率は等しく(式36に影響を与えない)、過去の尤度が完全に既知である場合には、正しい選択の場合確率は1であり、間違った選択の場合は0である。これが、そのビットのありそうもないあるいは間違った推測について−∞の、式36に対する貢献を強制する。
尤度関数の式36のこの単純化または量子化は、ターボ反復でのその点におけるビット決定が2つの部分S=[S]、b=[b ]、つまり「既知」と「未知」に分けることができることを数学的に暗示している。これは尤度を以下のように簡略化し、
Figure 0004490265
ここで最適化が未知のビット(2の下付き数字)に対してだけ実施され、既知のビット(1の下付き数字)は入力として計算に渡されることが理解される。さらに、尤度式、したがって未知のビットのための式38によって暗示される処理は、通常のMUD処理に同等であり、データを修正されたデータr−Sで置換する。
この実施形態では、任意の数のMUDアルゴリズムが置換できるが、MMSE MUDアルゴリズムが式38を計算するために使用される。MMSE演算子は
Figure 0004490265
となり、以下のアルゴリズム的なステップを使用して計算できる。
1)
Figure 0004490265
を入力し、ビット推定値が正しい確率を決定し、確実なビットと不確実なビットとを生成する
2)r、SAを入力し、残留
Figure 0004490265
を計算する。
3)
Figure 0004490265
を計算する。
4)
Figure 0004490265
を計算する。
5)
Figure 0004490265
を計算する。
6)
Figure 0004490265
を通過する。
部分的に量子化された過去の情報アルゴリズムを有する低複雑性Mudは、例えば図12に描かれている実施形態によって実施できる。
事前白色化されたデータに基づく低複雑性M−アルゴリズム
MUDにより処理されるデータは、多元接続干渉からユーザを部分的に切り離すように事前白色化(prewhitened)することが可能で、これが性能を劇的に改善する。図16は、本発明のある実施形態に従って事前白色化されたデータに基づくM−アルゴリズムを実現する。この実施形態のモジュール205、208、210及び213は、例えば、図1から図3のマルチユーザ検出器を実現するために使用できるであろう。
ここで図16の装置がさらに詳細に説明される。潜在的に複雑なデータストリーム200はなんらかのソースから受信される。符号分割多元接続(CDMA)通信方式のケースでは、データストリームはチップレートのなんらかの倍数でサンプリングされる。TDMA通信方式の場合、データストリームはシンボルレートのなんらかの倍数でサンプリングされる。
データ200は、なんらかのレート(例えば、シンボルレート)で転送されるデータのベクタを表す。このデータ200は整合フィルタ205に送信される。加えて、同じベクタ200は、パラメータ推定モジュール292に渡される。パラメータ推定モジュールの目的は、伝送チャネルのタイミング、信号振幅、位相、偏向及び識別を推定することである。パラメータの推定値は整合フィルタバンク205を設計するために渡され、パラメータの推定値は対応するホワイトナ208を設計するためにも渡される。
図16の213で実施されるシンボル仮説試験は、以下として数学的に表される最大尤度検出器を含むことができるであろう。
Figure 0004490265
式39は、受信されたサンプルを送信されたビット、b及びSで定義されたチャネルのモデルに関して定義する、
Figure 0004490265
により描かれている受信されたサンプルの線形モデルに基づいている。最大尤度検出器は、網羅的な検索を必要とする強引な手法である。最大尤度の解は、ビットシーケンスの考えられるあらゆる仮説を使用して、受信されたサンプルとサンプルの線形モデルの間のユークリッド距離を網羅的に評価することから構成される。多数のユーザまたはマルチパスからの厳しいシンボル間干渉の有する問題のためにはそれは計算集約的すぎる。ここで検討される該手法は、過飽和の環境で同じ性能をほぼ達成し、計算数で大きな節約を行う最大尤度検出器の簡略化されたバージョンから構成される。この手法は、以下のとおりに説明される。
最大尤度解は以下のように書き直され
Figure 0004490265
ここでΣは雑音、Nの共分散を表す。雑音が白色であるとき、式41の加重最小二乗法の解は式39の最大尤度検出器に同一である。任意の行列Wについて、式40の加重最小二乗法の解は以下の解に同一であり、
Figure 0004490265
ここでw=Wr及びΣ=WΣWである。受信データの線形組み合わせを探求する動機は、特定の変換がより可能性が高いビット仮説のさらに効率的な検索を可能にするためである。式41から42の表記法はすべてのユーザの存在に基づいた受信波形のサンプルのための線形行列表現に基づいている(式40を参照すること)。
Figure 0004490265
とし、次に行列Wで定められるフィルタバンクはデコレレイタ(decorrelator)受信機で使用されるフィルタバンクである。該デコレレイタ受信機は、それが最適に多元接続干渉を削減するが、有色雑音を説明しないため魅力的である。具体的には、既知の相互関係行列を仮定して多元接続干渉は排除されるが、白色雑音成分は(シグナチャ波形が直交でないときに)着色されている。さらに適切なフィルタバンクは相互関係行列の平方根の逆数を含む。この組み合わせの結果、多元接続干渉を部分的に減結合するが、無相関の雑音成分を維持するフィルタバンクが生じる。平方根フィルタと整合フィルタのカスケードは、シグナチャ波形に対して最小二乗の意味で最も近いフィルタの直交セットを表す。相関行列のコレスキー(Cholesky)因数分解を使用して定められる平方根フィルタバンクはより魅力的な平方根因数分解の1つである。相関行列は、以下
Figure 0004490265
により式43で表現される。
相関行列Hのコレスキー因数分解は、
Figure 0004490265
により定められ、ここでFは上三角行列である。したがって、白色化フィルタは、下三角行列であるF−Hとして定められる。したがって、式43でデコレレイトフィルタバンクを活用するより、
Figure 0004490265
として定められる以下の部分的なデコレレイトフィルタバンクは、加重最小二乗法の解の効率的な検索により適している。この特定の平方根因数分解の魅力は、以下
Figure 0004490265
で示されている。
同一チャネル干渉の部分的な減結合は、Fbであると定められる白色化された出力、wによって描かれており、ここでFは上三角(upper triangular)行列である。列ベクタbを、一番上の行が第1のユーザを表し、一番下の行がK番目のユーザを表すようにユーザにより並べ替えさせる。Fbの中の項を計算すると、K番目のユーザが他のすべてのユーザのビット仮説から完全に切り離されていることを示す。また、(K−1)番目のユーザビット仮説はユーザKのビット仮説とだけ結合される。用語部分的な減結合は、(K−m)人のユーザのための決定が、第1の1人から(K−m−1)人のユーザの知識が後の(K−m)人のユーザに対する決定を下すために必要とされないように、他のユーザのどれかから切り離されるため、使用される。
このようにして続行すると、任意のユーザのための測定が任意の「将来の」ユーザの実際のビットからどのようにして切り離されたのかが示される。ユーザkの項「将来」が、すべてのユーザ1からk−1までを指すことに注意せよ。
雑音は、
Figure 0004490265
によって定められる部分的なデコレレイタを使用して白色化される。白色雑音は、
Figure 0004490265
によって示され、ここでIは識別行列を表し、Eはランダム変数の期待値を表す。斜めの共分散(diagonal covariance)行列は、雑音が式48に定められる部分的にデコレレイトするフィルタバンクを使用して白色化されたことを証明している。式48の中のデコレレイトフィルタバンクを式42に代入すると、同じ最大尤度解が生じる。
式42の最大尤度表現は
Figure 0004490265
によって示される測定基準Ω(b)に関して書き直され、
ここで、
Figure 0004490265
式51を使用すると、ビットの最適集合のための検索が、ユーザkのためのビット仮説の尤度を特徴付ける測定基準、bが構成要素、
Figure 0004490265
で表現されるデシジョンツリー(decision tree)に関して再公式化することができる。
項bは、ユーザkのビット仮説を表し、項wは、ユーザKによって使用されるシグナチャ波形に一致されたフィルタkのフィルタバンク出力を表す。項Fjkはユーザjとユーザkについて定められたコレスキー係数を表す。ML測定基準の連続的な性質は以下の式によりさらに明確に示される。総和(式52を参照)の各構成要素は、デシジョンツリーのK個の段階の1つと見なすことができる。以下の式は、デシジョンツリーの最初の3つの段階の構成要素に相当する距離測定基準の最初の3つの項を描いている。
Figure 0004490265
これらの最初の3つの段階で測定基準を評価することは、QPSKアルファベットについて図17で示されている。第1の段階では、第1の構成要素
Figure 0004490265
がユーザKのための4つの考えられる状態の尤度を評価する。図17により示されるように、この尤度測定基準は他のユーザのビットに関する他の決定に依存していない。第2の段階では、(K−1)番目のユーザのビットb(K−1)の推定値は、
Figure 0004490265
に基づいている。第2段階での決定は、ユーザb(K−1)に対応するフィルタリングされたデータだけに依存し、決定はK番目のユーザのために下される。同様に、(K−2)番目のユーザの推定値は、過去の2つの決定に依存する距離測定基準の第3の構成要素を使用して第3の段階で評価される。これは、これらの構成要素のすべてKについて続行する。傾向は、ユーザのビットに関する初期の決定が、残りのユーザのビットのために下される「将来の」決定から切り離されることを示している。
ML測定基準全体がこれらの構成要素の内のKの総和であるため、デシジョンツリーに対するK個の段階がある。結合して最適な決定は、デシジョンツリーの枝のそれぞれの枝を探索することを必要とする。デシジョンツリー手法は、マルチパスがないQPSKケースについて4の仮説が評価されることをまた必要とする。デシジョンツリーに関して問題を表現することにより、M−アルゴリズムまたはT−アルゴリズムなどのプルーニング技法の探求が可能になる。
過飽和状態の通信の問題に対しては、統計学的に独立した次元より多くのユーザが存在するため、相関行列が半正値(positive semidefinite)(つまり、反転可能ではない)となるように、過少決定(underdetermined)された問題がある。相関行列は反転可能ではないため、式48の部分的なデコレレイタ(decorrelator)は、チャネル間の相互関係のためにもはや存在していない。
マルチ共直線性(multicolinearity)に対抗するための回帰分析で一般的に使用される技法がリッジ回帰(ridge regression)と呼ばれる。マルチ共直線性は、独立変数間の高い相互関係から生じ、この特定のケースではユーザの送信済み波形に相当する。単に置かれた(ASA)は大きな非対角(off-diagonal)項を有し、高い条件数の不安定な相関行列を生じさせる(つまり、最大固有値対最低固有値の比)。回帰分析では、これは非常に高い共分散の推定値を生じさせる。これは、推定値の中の分散を最小限にするために小さなバイアスを受け入れることで解決される。これは、相関行列を(ASA)から(ASA+γI)に変換することにより行われ、ここでγは対角ロードを表す。これが最小の固有値を上昇させ、回帰分析のバイアスがかけられた推定値を犠牲にして、さらに安定した相関行列を生じさせる。
ここで検討される手法は、デシジョンツリーで固有の非線形決定プロセスにより重大となる推定値の分散を削減することを目的として、相関行列を対角でロードし、このようにして意図的にバイアスを生じさせることにより部分的にデコレレイト(decorrelating)するフィルタバンクを改変することである。具体的には、新しい部分的なデコレレイトフィルタバンクは、
Figure 0004490265
で定められ、ここで対角にロードされた相関行列のコレスキー因数分解は、
Figure 0004490265
ほどである。
前述されたように、リッジ回帰の概念は、推定値の分散を削減する意図をもってバイアスを意図的に導入することから構成されている。式54で定められる部分的にデコレレイトするフィルタバンクを受信されたサンプルに適用すると、K個の要素の列ベクタwにより定められるフィルタバンクの中からベクタK個のサンプルが生じる。式54を受信されたサンプルに適用することは、
Figure 0004490265
で表される。
式40に定められる受信されたサンプルの行列モデル表現を代入すると、結果として生じる簡略化が実現される。
Figure 0004490265
ここで相関行列Hは式44で定められる。式57に基づき、雑音分散を使用して対角ロードを通して意図的に導入されたバイアスが、
Figure 0004490265
であることは明らかである。
対角にロードする相関行列のコレスキー因数分解は、以下
Figure 0004490265
により定められる式57の雑音の項の共分散により示されるように雑音を完全に白色化しない。
妥当な信号対雑音比、及び穏当な対角ロードのケースでは、対角にロードされた基底付きの(based)部分的なデコレレイタを適用した後の雑音共分散は
Figure 0004490265
により近似される。加えて、小さな対角ロードにとって、そのバイアスは小さいと考えられるこれらの近似に基づいて、加重最小2乗法の解は
Figure 0004490265
により近似される。
前記のように、最大尤度の解は、
Figure 0004490265
で表され、ここで
Figure 0004490265
である。
式62を評価するための手順は、
Figure 0004490265
により記述される個々の測定基準がデシジョンツリーの各ノードで評価される図17に説明される手法と一貫している。図17を参照すると、デシジョンツリーの過去の段階から現在の段階への遷移は、デシジョンツリーの経路に沿った個別測定基準の発生から構成されている。M−アルゴリズムとT−アルゴリズムのような多様な次善のプルーニング手法は、許容できる誤り率のデシジョンツリーを効率的に横切るために使用できる。
図16は、式55、56、61及び62に記述されるアルゴリズムを実現する。式56は、整合フィルタバンクの適用に続いて対角ロードされた相関行列のコレスキー因数分解を適用することから構成される。したがって、式56に記述される整合フィルタバンクは205で適用される。部分的にデコレレイトするフィルタバンクを完成するために、オーバロードした非同期ホワイトナ(whitener)応用210がフィルタリングされたデータ207に適用される。210の非同期ホワイトナは、式55に記述される対角にロードされた相関行列のコレスキー因数分解を適用することから構成される。
パラメータ推定モジュール202の目的は、アクティブチャネルのタイミング、信号振幅、位相、偏向及び識別を推定することである。該パラメータの推定値は、整合フィルタの適用及び非同期デコレレイトフィルタバンクの作成に必要とされるチャネルをモデル化するために使用される。パラメータ推定モジュールはチャネルモデル206及び相関行列203を提供する。
該208の目的は、ボックス210で使用されるホワイトナを設計することである。ボックス208は、式55に記述される対角でロードされるコレスキー因数分解を計算するためにパラメータ推定モジュールによって提供される相関行列203を活用する。因数分解はボックス210で白色化に使用され、ボックス213で仮説試験に使用される。この手法はデコレレイト決定フィードバック検出の拡張であるため、「受信エネルギーを減少する」ことによりユーザを順序付けるという概念は、最強のユーザが最初に評価されるように維持される。この順序付けは、平方根行列の型を定め、仮説試験モジュール213で維持されなければならない。変形物は、SNRベース及び尤度ベースなどの他の順序付け技法を含む。ボックス208で使用される手法は、相関行列のブロックバンド(block−banded)構造を利用する非同期因数分解を含む。記号同期の問題は退化したケースであり、したがって本発明の実施形態により対処され得る。ウィンドウ化ベースの技法もまたここで利用されてよい。
ボックス208は、語長を削減し、処理速度を改善するためにも適用されてよい。例えば、対角行列で増補されたシグナチャ波形行列から構成される、行列で実現されるハウスホルダ変換を使用するQR因数分解は、対角にロードされる相関行列のコレスキー因数分解より少ない語長を必要とする。少ない語長を使用するアルゴリズムインプリメンテーションは、固定点処理ハードウェア構成により適している。加えて、本発明の一実施形態は、受信エネルギー及び/または位相だけがシンボル期間の間に変化するときにホワイトナを効率的に更新するためにボックス208での双曲線ハウスホルダ(Hyperbolic Householder)変換の使用を含む。
ボックス208内で生成される相関行列の平方根因数分解はボックス210でデータを白色化するために使用される。本発明の実施形態は、整合フィルタリングされたデータを白色化する2つの手法をカバーするように適用できる。第1の手法は、平方根行列の共役転置の逆数により定められるフィルタのバンクを適用することに基づくことができる。平方根は三角形構造で定められているため、後方置換を使用する白色化手順を実現できる。この代替手法が必要とする演算の数はより少ない。
白色化されたデータストリームはボックス210を出て、線路212上でボックス213に渡される。ボックス208に定められる平方根因数分解は、線路211上でボックス213に渡される。この因数分解は、M−アルゴリズムまたはT−アルゴリズムなどの横幅第1(breadth-first)技法を使用して実現できるデシジョンツリー内でビット仮説を連続して評価するための測定基準で使用される。相関行列因数分解及び白色化フィルタを定義するために使用されるユーザ順序付けがシンボル仮説試験モジュール、ボックス213で維持されることを思い出せ。
ボックス213の目的は、K人すべてのユーザについてより可能性が高いビット仮説を効率的に調査することである。一実施形態では、ボックス213で実施されるシンボル仮説試験が、図17に記述される仮説の尤度を特徴付ける測定基準の順次評価に基づいている。この評価は等式61と62に記述される測定基準に基づいている。
デシジョンツリー内の段階である特定のユーザのビット仮説に対応する測定基準がここで言葉で説明される。該測定基準は、線路212に沿った白色化フィルタの1つの出力と問題のユーザのためのビット仮説と過去に試験されたユーザのために選択された仮説に対応する平均信号エネルギーに基づいて仮定される平均信号エネルギーの間のユークリッド距離から構成される。この平均信号エネルギーは、ボックス208で計算され、線路211の上で提供される対角にロードされる相関行列のコレスキー因数分解に基づいている。図17に描かれているデシジョンツリーの各ノードでの測定基準は、過去の決定に対応する測定基準の蓄積を含む。
決定フィードバック手法とは異なり、決定は直ちに下されない。しかしながら、ここで検討される手法は、決定フィードバック技法の汎用化であり、したがって決定フィードバック技法を含む。この連続した概念は、すべてのフィルタバンクサンプルと任意のビット仮説に基づいた仮定された平均信号の間のユークリッド距離の項を拡大することにより、図17と式53で観察された。ボックス213の出力は、線路214で報告される重要なシンボル期間でのすべてのユーザのためのシンボルの抑制された推定値である。
多様な効率的なデシジョンツリー検索戦略は、ボックス213で利用できる。例えば、M−アルゴリズムは各段階での仮説の数を固定数に制限する1つのこのような手法である。T−アルゴリズムは事実上M−アルゴリズムに類似しているが、それは発生した測定基準を閾値に比較することで仮説の数を制限する。デシジョンツリー検索に対するこの手法及び他の効率的な手法の変形物はこの開示を鑑みて明らかとなるであろう。
本開示を鑑みて明らかになるであろう、本発明の原理に従って構成される同一チャネル干渉受信機のための多数の代替実施形態がある。例えば、MUDは多様なターボMUD設計に準拠し、MUDコアとエラー訂正デコーダのためのさまざまな組み合わせを可能にする。MUDコアは、多様なM−アルゴリズムインプリメンテーションだけでなくMMSE MUD、ウィンドウ化した最適MUDを備えることがある。さらに、ターボループ内の畳み込みエラー訂正デコーダは、とりわけソフト出力ビタビ(Viterbi)またはソフト出力BCJRを備える場合がある。
本発明の実施形態は、ソフトウェア、ハードウェア、ファームウェアまたはその任意の組み合わせで実現されてよい。例えば、パラメータ推定モジュールとMUDモジュールは、デジタル信号プロセッサで実行中の命令の集合で、あるいは他の適当な処理環境で実現されてよい。代わりに、モジュールは専用のシリコン(例えば、ASIC、FPGA、または他のチップ上システム設計)で実現できる。
本発明の実施形態の前記説明は、図解及び説明のために提示されてきた。網羅的となる、あるいは開示されている正確な形式に本発明を制限することは意図されていない。本開示を鑑みて多くの修正及び変形が考えられる。本発明の範囲がこの詳細な説明によってではなくむしろここに添付される請求項によって制限されることが意図される。
本発明の一実施形態に従って構成されるマルチユーザ検出受信機を有する通信システムのシステムブロック図である。 本発明の別の実施形態に従って二重偏向アンテナを有するマルチユーザ検出受信機を示す図である。 本発明の別の態様に従ってそれぞれのアンテナが1つまたは2つの偏向を有する、複数のアンテナを有するマルチユーザ検出受信機を示す図である。 指定されたダイバシティポート信号のための2つまたは3つ以上の同一チャネル干渉信号のケースについて受信されたベースバンド信号の下にあるフレーム構造の図であり、本発明の一実施形態に従ってトレーニングシーケンスロケータで使用するためのトレーニングシーケンススライド式検索ウィンドウも示す図である。 本発明の一実施形態に従って構成されるパラメータ推定器のブロック図である。 本発明の一実施形態に従って構成されるトレーニングシーケンスロケータ構成要素のフローチャートである。 本発明の一実施形態に従って構成される雑音推定器の構成要素のブロック図である。 本発明の一実施形態に従って構成されるシグナチャ波形推定器構成要素のフローチャートである。 本発明の一実施形態に従って構成されるアクティブユーザテスタ構成要素のためのフローチャートである。 本発明の一実施形態に従ってダイバシティポート「p」用のアクティブユーザテストのフローチャートである。 本発明の一実施態様に従ったターボMUDの中の入信データストリームの反復処理の図表示である。 本発明の一実施形態に従って構成されるターボMUDシステムを描くブロック図である。 図12のシステムにより実行される例のスレッショルド(thresholding)の図表示である。 図12のシステムにより実行される例のスレッショルド及び結合を描くブロック図である。 特定の値を有する図14の例のK掛けるLの行列内のビットの2つの概略図である。 図14の行列からの値の例の処理の概略図である。 図14のエラー訂正装置に適切な推定値を提供するために、未処理信号と結合される次の反復で処理されるMUD信号の概略図である。 本発明の実施形態に従って事前白色化されたデータに基づいたMアルゴリズムを実現する図である。 本発明のある実施形態に従って仮説の尤度を特徴付ける測定基準の順次評価に基づいたシンボル仮説試験を立証する図である。

Claims (19)

  1. K個の同一チャネル干渉合成信号を受信する同一チャネル干渉受信機であって、
    前記受信機により受信される前記合成信号を事前に白色化するように構成され、多元接続干渉からユーザを部分的に切り離すように動作可能な整合フィルタと、
    K個の同一チャネル干渉信号からの情報を含む前記合成信号を受信するように構成され、前記K個の同一チャネル干渉信号を復調するために、前記事前に白色化されたデータに基づいた計算量の少ないM−アルゴリズム動作するように構成されたマルチユーザ検出器モジュールと、
    前記合成信号を受信し、前記K個の同一チャネル干渉信号のそれぞれの推定シグナチャ波形を生成して、前記K個の同一チャネル干渉信号を復調する前記マルチユーザ検出器モジュールに、前記推定シグナチャ波形を提供するように構成されるパラメータ推定モジュールと、
    前記パラメータ推定モジュールに動作可能に結合され、前記パラメータ推定モジュールからのパラメータ推定値に基づいてそれぞれの受信合成信号のモデルを作成し、前記整合フィルタにより出力されるフィルタリングされたデータを白色化する非同期ホワイトナモジュールを計算するように構成されるホワイトナデザイナモジュールと、
    前記ホワイトナデザイナモジュールに動作可能に結合され、前記非同期ホワイトナモジュールにより出力される白色化されたデータを受信するように構成され、仮説の尤度を特徴付ける測定基準の連続した評価に基づいてシンボル仮説試験を実施するように構成されるシンボル仮説試験モジュールと、
    を含、同一チャネル干渉受信機。
  2. 前記ホワイトナデザイナモジュールは、前記パラメータ推定モジュールにより提供される相関行列を活用し、前記非同期ホワイトナモジュールでの白色化のために使用され、前記シンボル仮説試験モジュールでの仮説試験にも使用される対角にロードされるコレスキー因数分解を計算する、請求項に記載の受信機。
  3. 前記ホワイトナデザイナモジュールは、ハウスホルダ変換を使用するQR因数分解を利用する、請求項に記載の受信機。
  4. 前記ホワイトナデザイナモジュールは、受信されたエネルギー及び/または位相だけがシンボル期間の間に変化するときに前記非同期ホワイトナモジュールを効率的に更新する
    ために双曲線ハウスホルダ変換を利用する、請求項に記載の受信機。
  5. 前記ホワイトナデザイナモジュールは平方根因数分解を利用する、請求項に記載の受信機。
  6. 前記非同期ホワイトナモジュールはデータを白色化するために前記ホワイトナデザイナモジュールにより生成される相関行列平方根因数分解を利用する、請求項1乃至5の何れか1項に記載の受信機。
  7. 前記非同期ホワイトナモジュールは、前記相関行列平方根因数分解の共役転置の逆数により定められるフィルタのバンクを利用する、請求項に記載の受信機。
  8. 前記相関行列平方根因数分解は三角形構造を有し、前記非同期ホワイトナモジュールは後方置換を利用する、請求項に記載の受信機。
  9. 前記シンボル仮説試験モジュールは、横幅第1技法を使用して実現されるデシジョンツリー内で前記ビット仮説を連続して評価するために測定基準前記ホワイトナデザイナモジュールにより生成される相関行列平方根因数分解を利用する、請求項に記載の受信機。
  10. 前記パラメータ推定モジュールパラメータ推定値を提供し、前記パラメータ推定値は、各受信合成信号と関連付けられるチャネルをモデル化するために使用され、前記整合フィルタの適用及び非同期デコレレイトフィルタバンクの作成を可能にする、請求項1乃至9の何れか1項に記載の受信機。
  11. 前記パラメータ推定モジュールは、
    前記受信合成信号の各フレーム内でトレーニングシーケンスロケーション指数を推定するように構成されるトレーニングシーケンスロケータモジュールと、
    前記トレーニングシーケンスロケーション指数に従って前記受信合成信号の平均雑音電力の推定値を計算するように構成される雑音推定器モジュールと、
    前記トレーニングシーケンスロケーション指数及び変換行列に従って前記受信合成信号内で各ユーザに一意のシグナチャ波形を推定するように構成されるシグナチャ波形推定器モジュールと、
    前記雑音推定器モジュールの出力と、前記シグナチャ波形推定器モジュールの出力とに動作可能に結合され、前記受信合成信号に関連付けられるアクティブユーザの数を決定するように構成されるアクティブユーザテスタモジュールと、
    前記アクティブユーザテスタモジュールと、各ユーザの事前に記憶された既知のトレーニングシーケンスとに動作可能に結合され、前記シグナチャ波形推定モジュールにより使用される前記変換行列を生成するように構成される変換行列リビルダモジュールと、
    を含む、請求項1乃至10の何れか1項に記載の受信機。
  12. 前記マルチユーザ検出器モジュール及び前記パラメータ推定モジュールに動作可能に結合され、複数の送信機からそれぞれ1つまたは複数の複合波形信号を受信し、それぞれの受信された複合波形信号を対応する合成信号に変換するように構成されるアナログフロントエンドをさらに含む、請求項1乃至11の何れか1項に記載の受信機。
  13. 前記アナログフロントエンドは、
    それぞれが受信された複合波形をデジタル波形に変換するように構成される1または複数のアナログ/デジタル変換器と、
    それぞれがそれぞれのアナログ/デジタル変換器に動作可能に結合され、受信された複合波形に関連付けられる周波数をさらに低い周波数に変換するように構成される1または複数のダウンコンバータと、
    を含む、請求項12に記載の受信機。
  14. 前記アナログフロントエンドは、
    それぞれが前記複数の送信機から対応する複合波形信号を受信するように構成された1または複数のアンテナを含む、請求項12又は13に記載の受信機。
  15. 前記1または複数のアンテナは単一偏向アンテナである、請求項14に記載の受信機。
  16. 前記1または複数のアンテナは、2つの偏向ポートと構成される二重偏向アンテナであり、それにより偏向ダイバシティを提供する、請求項15に記載の受信機。
  17. 前記1または複数のアンテナは、それぞれが2つの偏向ポートと構成される以上の二重偏向アンテナを含み、それにより空間ダイバシティ及び偏向ダイバシティを提供する、請求項14に記載の受信機。
  18. 前記パラメータ推定モジュールは、アクティブチャネルのタイミング、信号振幅、位相、偏向及び識別のうちの少なくとも1つを含む信号パラメータを推定するように構成される、請求項1乃至17の何れか1項に記載の受信機。
  19. K個の同一チャネル干渉信号からの情報を含む少なくとも1つの合成信号を受信する方法であって、
    a)前記少なくとも1つの合成信号をフィルタリングして、多元接続干渉からユーザを部分的に切り離し、事前に白色化されたデータを提供するステップと、
    b)前記K個の同一チャネル干渉信号のそれぞれのシグナチャ波形を推定するステップと、
    c)事前に白色化されたデータに基づく計算量の少ないM−アルゴリズムで前記推定されたシグナチャ波形に基づいて前記少なくとも1つの合成信号を処理するステップと、
    d)前記少なくとも1つの合成信号のモデルをパラメータ推定値に基づいて作成するステップと、
    e)前記事前に白色化されたデータを白色化するための前記モデルに基づいて非同期ホワイトナを計算するステップと、
    f)仮説の尤度を特徴付ける測定基準の連続した評価に基づいてシンボル仮説試験を実施するステップと、
    g)K個の同一チャネル干渉信号についての情報を提供するために、前記少なくとも1つの合成信号を復調するステップと
    を含む、方法
JP2004522953A 2002-07-24 2003-04-25 同一チャネル干渉受信機 Expired - Lifetime JP4490265B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US39845102P 2002-07-24 2002-07-24
US10/228,787 US6947502B2 (en) 2002-04-16 2002-08-26 Parameter estimator for a multiuser detection receiver
PCT/US2003/012917 WO2004010572A1 (en) 2002-07-24 2003-04-25 Co-channel interference receiver

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2009290622A Division JP4958021B2 (ja) 2002-07-24 2009-12-22 同一チャネル干渉受信機

Publications (2)

Publication Number Publication Date
JP2005534226A JP2005534226A (ja) 2005-11-10
JP4490265B2 true JP4490265B2 (ja) 2010-06-23

Family

ID=30772568

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2004522953A Expired - Lifetime JP4490265B2 (ja) 2002-07-24 2003-04-25 同一チャネル干渉受信機
JP2009290622A Expired - Lifetime JP4958021B2 (ja) 2002-07-24 2009-12-22 同一チャネル干渉受信機

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2009290622A Expired - Lifetime JP4958021B2 (ja) 2002-07-24 2009-12-22 同一チャネル干渉受信機

Country Status (4)

Country Link
EP (1) EP1554798B1 (ja)
JP (2) JP4490265B2 (ja)
AU (1) AU2003223740A1 (ja)
WO (1) WO2004010572A1 (ja)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003272709A1 (en) * 2003-09-25 2005-05-11 Bae Systems Information And Electronic Systems Integration Inc. System for decreasing processing time in an iterative multi-user detector system
JP4403010B2 (ja) * 2004-02-03 2010-01-20 株式会社エヌ・ティ・ティ・ドコモ 信号分離装置
CN1965501B (zh) * 2004-02-13 2012-07-04 日本电气株式会社 无线通信系统、接收设备、解调方法
US7233777B2 (en) * 2004-08-31 2007-06-19 L-3 Integrated Systems Company Separation of AM cochannel signals in an overloaded signal environment
US7711035B2 (en) 2004-09-17 2010-05-04 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for suppressing communication signal interference
JP4666150B2 (ja) 2005-05-31 2011-04-06 日本電気株式会社 Mimo受信装置、受信方法、および無線通信システム
JP4545663B2 (ja) * 2005-09-06 2010-09-15 日本電信電話株式会社 マルチユーザ受信装置およびマルチユーザ受信方法
US7817101B2 (en) 2006-10-24 2010-10-19 Com Dev International Ltd. Dual polarized multifilar antenna
US7876865B2 (en) 2007-06-08 2011-01-25 COM DEV International Ltd System and method for decoding automatic identification system signals
PL2329629T3 (pl) * 2008-08-20 2013-08-30 Qualcomm Inc Przesunięcie wieloramek pomiędzy sąsiednimi komórkami dla poprawy wydajności sygnalizacji sieci GERAN przy stosowanym sterowaniu mocą sygnalizacji
CA2743442A1 (en) * 2008-11-13 2010-05-20 Nortel Networks Limited Reduced complexity channel estimation for uplink receiver
CN101420708B (zh) * 2008-11-28 2010-12-08 华为技术有限公司 一种多用户异步信息检测方法和装置
JP5475276B2 (ja) * 2008-12-24 2014-04-16 三星電子株式会社 受信装置、及び信号検出方法
PL2430869T3 (pl) 2009-05-11 2014-05-30 Qualcomm Inc Przesuwanie czasowe transmisji danych wspólnego kanału w celu zmniejszenia zakłóceń wspólnego kanału
US8331495B2 (en) * 2009-08-17 2012-12-11 Telefonaktiebolaget L M Ericsson (Publ) Method and apparatus for an improved wireless communication receiver
US8780788B2 (en) 2009-09-25 2014-07-15 Com Dev International Ltd. Systems and methods for decoding automatic identification system signals
JP2011211335A (ja) * 2010-03-29 2011-10-20 Japan Radio Co Ltd 無線パケット受信装置
US9331774B2 (en) 2010-06-09 2016-05-03 Exactearth Ltd. Systems and methods for segmenting a satellite field of view for detecting radio frequency signals
KR101317593B1 (ko) 2011-10-07 2013-12-19 경상대학교산학협력단 비동기 합동 검출 방법 및 장치
US9015567B2 (en) 2012-04-12 2015-04-21 Com Dev International Ltd. Methods and systems for consistency checking and anomaly detection in automatic identification system signal data

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5905742A (en) * 1995-12-27 1999-05-18 Ericsson Inc. Method and apparauts for symbol decoding
US6172970B1 (en) * 1997-05-05 2001-01-09 The Hong Kong University Of Science And Technology Low-complexity antenna diversity receiver
US5867478A (en) 1997-06-20 1999-02-02 Motorola, Inc. Synchronous coherent orthogonal frequency division multiplexing system, method, software and device
US6144710A (en) * 1998-04-23 2000-11-07 Lucent Technologies, Inc. Joint maximum likelihood sequence estimator with dynamic channel description
US6377611B1 (en) * 1999-02-01 2002-04-23 Industrial Technology Research Institute Apparatus and method for a DSSS/CDMA receiver
GB2354676B (en) * 1999-09-20 2001-10-10 Lucent Technologies Inc Turbo equalization with a decision feedback equalizer

Also Published As

Publication number Publication date
AU2003223740A1 (en) 2004-02-09
JP2010093852A (ja) 2010-04-22
EP1554798A4 (en) 2009-03-18
EP1554798A1 (en) 2005-07-20
EP1554798B1 (en) 2018-05-30
WO2004010572A1 (en) 2004-01-29
JP4958021B2 (ja) 2012-06-20
JP2005534226A (ja) 2005-11-10

Similar Documents

Publication Publication Date Title
JP4958021B2 (ja) 同一チャネル干渉受信機
US7092452B2 (en) Co-channel interference receiver
EP1471655B1 (en) Deferred decorrelating decision-feedback detector for supersaturated communications
US5889827A (en) Method and apparatus for digital symbol detection using medium response estimates
US7583757B2 (en) Joint symbol, amplitude, and rate estimator
US7010069B2 (en) Method for co-channel interference identification and mitigation
RU2484582C2 (ru) Способ и устройство для последовательного вычитания помех с помощью обработки матрицы корня ковариации
US7031671B2 (en) Multi-user detection method
US8421654B2 (en) Method of decoding a signal implementing a progressive construction of a decoding tree, corresponding computer program and decoding device
US20080095259A1 (en) Pre-coding for multiple-input-multiple-output communications
US7817754B2 (en) M-algorithm with prioritized user ordering
US20070076669A1 (en) Multi-channel communication method and apparatus using plural Markov Chain Monte Carlo Simulations
Li et al. Timing estimation and resynchronization for amplify-and-forward communication systems
KR101216343B1 (ko) 간섭에 대한 최소 제곱 추정을 이용한 수신
JP2002111563A (ja) マルチユーザ検波方法およびマルチユーザ検波装置
KR20070118837A (ko) 다중 안테나 시스템에서 개선된 파노 알고리즘을 이용한신호 검출 장치 및 방법
US8018985B2 (en) CDMA receivers for the golden code
Ghauri et al. Blind and semi-blind single user receiver techniques for asynchronous CDMA in multipath channels
Correal Efficient interference cancellation techniques for advanced DS-CDMA receivers
Wan et al. Joint channel estimation and Markov Chain Monte Carlo detection for frequency-selective channels
KR100816057B1 (ko) 다중 송수신 시스템에서의 신호검출 방법 및 수신 장치
Slimane et al. An improved successive interference cancellation multiuser detector for narrowband signals
Zhou et al. A soft tree pruning based fixed-complexity sphere decoder for interference-limited MIMO systems
Dai et al. Multiuser receiver design
Gu et al. Study of Interference Cancellation and Relay Selection Algorithms Using Greedy Techniques for Cooperative DS-CDMA Systems

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080115

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20080402

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20080409

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080603

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090630

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090925

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20091002

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091222

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20100302

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100401

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130409

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130409

Year of fee payment: 3

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130409

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130409

Year of fee payment: 3

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130409

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140409

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250