JP2008514075A - Ｑａｍ変調における伝送チャネルで送信される観測データの移送とゲインの推定方法 - Google Patents

Abstract

ＱＡＭ変調で伝送される観測データのゲイン又は位相を推定するための方法に関する。当該方法は、（Ａ）特定の位相とゲインの法則に基づいて位相及び/又はゲインを反復推定すること及び、（Ｂ）位相及び/又はゲインパラメータを推定するための適応型プロセスを実行する、なお、適応型プロセスは、ＱＡＭ変調のシンボルを構成するビット組に関して各観測データ（ｙ_ｋ）について、少なくとも、log尤度確率で表現される尤度確率値に基づいて推定する少なくとも１つの関数を有する。本発明は、単一キャリア又はマルチキャリアのデジタル伝送に適用される。
【選択図】図２ａ

Description

現在のデジタル伝送システムにおいて、伝送デジタル信号は時間的連続のアナログ信号に変換される。該アナログ信号は伝送チャネルと呼ばれる物理的伝送路を通して伝送される。例えば、無線波に対しては空中、光波に対しては光ファイバを通して伝送される。受信時には、受信信号は伝送チャネルとの物理的相互作用を受け、処理されデジタル形式に変換される。

通常、伝送の処理ステップは次を含む。
―伝送される２値、ビットの組合せを複素数シンボル(complex symbols)の組合せへ変換すること。前記複素数シンボルの組合せは有限のアルファベットで、複素数面におけるコンステレーションの形式で表現される；
―前記シンボルの組を、そのスペクトルが周波数０の周辺に集中するベースバンドの時間連続の波形へ変換すること；
―キャリア周波数の周辺で周波数を変換することによりシフトすること；
受信時には、対の動作（paired operations）が実行される：
―コンプレックス デモジュレータ（complex demodulator）を使って逆シフト操作（reverse shifting）してベースバンドに復帰すること；
―ベースバンドの時間連続波形を、複素数値の組に変換すること；
―伝送される２値を再構成すること。

送信時における周波数シフト操作と受信時におけるベースバンド復帰の操作は、２つの独立の発振器により制御される。１つは発振に対し、他のものは受信に対するものである。周波数は、まして位相は完全に一致することはない。位相エラーは一般的に存在し、時間的に変化するものである。

前記エラーの大部分はアナログ装置を使って通常のやり方で補償されるが、周波数変動の回路レベルに適用するＰＬＬのような回路については、常に、ベースバンドキャリア周波数の残留(residue)及びコンプレックス デモジュレータの出力において最終的に位相誤差がある。

送受信時に発振器に影響を与える周波数及び/又は位相偏移（difference）は擾乱要素を構成し、コンプレックス デモジュレータから出力される信号の観測における寄生的位相ずれ誤差を招く。

別の複数の要因が前記寄生的位相ずれ誤差を強める。特に伝送チャネルを通過する信号が必要とする伝播時間、発振器・受信器の間の相対的動き(relative movement)はドップラー効果を生じ、それは擾乱の位相ずれ誤差を招く。

受信信号を適正に処理し、送信シンボルを抽出し、十分な精度を識別するためには位相ドリフトの補償が不可欠である。

公知の位相ドリフト補償技術をベースバンドの信号について説明する。ベースバンドの信号に対して、寄生的位相ずれ誤差θ_ｋの効果は観測データｙ_ｋに付加される。観測データｙ_ｋは複素数データ（compleｘ data）で、受信器のコンプレックス デモジュレータから出力されるもので、次の関係式で表現される。
ｙ_ｋ＝ａ_ｋｅ^ｉθｋ+ｂ_ｋ

この関係式中、ａ_ｋは送信された複素数シンボルを意味し、Ｍ個の要素を持つＱＡＭ変調におけるＭ要素の有限のアルファベット｛Ｑ１，Ｑ２…ＱＭ｝に属している。ＱＡＭとはQuadrature Amplitude Modulationである。但し、
ａ_ｋ∈｛Ｑ１，Ｑ_２…Ｑ_Ｍ｝、Ｍ＝２^Ｎ、
Ｎは複素数シンボルａ_ｋを構成するのに使うバイナリーのパケットの長さを意味する。
-ｂ_ｋは付加的雑音で、ガウス雑音、白色雑音、サーキュラー雑音及び中心雑音であると仮定される。

寄生的位相ずれ誤差θｋを推定し、訂正することを可能にする公知技術の中で、最も洗練された推定方法は最も負荷の大きいデジタル処理、即ち、マルコフ連鎖等によるモンテカルロ法等である。この方法は受信した観測データ全体を同時に処理するものである。

しかし、これら方法を実際に使うことは困難であり、不可能である。何故なら前記方法はリアルタイムで膨大な計算を必要とするからである。

使用するときの簡便性から、シーケンシャルに、連続的に受信された観測データを処理するためには、位相ロックループ、Phase Locked Loop:ＰＬＬ技術が望ましい。

通常、位相ロックループはデジタルの反復アルゴリズムで、これにより位相の値、従って、寄生的位相ずれ誤差を推定することができる。前述のアルゴリズムとデジタル処理は変調のタイプに密接に依存している。

例示としてあげるが、２状態位相偏移変調ＭＤＰ２、又、Binary Phase Shift Keying:BPSKの場合、送信時の送信シンボルは、-１、又は+１である。伝送チャネルがもたらす寄生的位相ずれ誤差θｋのために、コンプレックスデモジュレータの出力として受信時に得られる複数の観測データは、もはや-１、又は+１に対応せず、図１で表わされる、位相のずれた値である。

従来の位相ロックループを使って、真の位相、従ってBPSK変調の場合の寄生的位相ずれ誤差θｋを推定することができる。従来の位相ロックループは、COSTASループであり、これによれば、現行観測データ(current observation datum)ｙ_ｋの位相φ_ｋを、現行観測データｙ_ｋと位相の以前の推定(previous estimate)φ_ｋ−1を使って次の反復公式に基づいて推定できる。
φ_ｋ＝φ_ｋ−１+γIm（ｙ^２ _ｋｅ^{−ｉ２φｋ―１}） （２）

別の変調タイプとして、特に４つ（quadrature）の内２のキャリアの変調、ＱＡＭと呼ばれるQuadrature Amplitude Modulationが用いられる。

一般的に、どのような変調タイプが使用されようとも、位相ロックループは次の関係式を満たす。
φ_ｋ＝φ_ｋ-１+γF（ｙ_ｋ、φ_ｋ-１） （３）

このタイプの全ての位相ロックループは、Ｍ−ＱＡＭ変調タイプに密接に関係する関数Ｆを使って、以前の位相φ_k−1に基づいて現在の位相φ_ｋを算出するように設計されている。

更に、パラメータγは、第２次フィルタ機能、比例と積分補正、又は更に高度次数のフィルタ関数により形成される。

前記位相ロックループは、従来のアナログモデルを採用しているが、この位相ロックループには大きな限界がある。即ち、位相推定φ_ｋは本質的に以前の推定値φ_ｋ-１に基づくという事実、及び１つの又は複数の過去の観測データ、及び/又は現在の観測データ及び１又は複数の過去の推定の関数に基づくという事実により限界がある。

このために、現在の位相の推定及びその現在の位相の補正は、大半、準最適な（sub-optimal）ものである。

例えば、エラー訂正の必要性に対して観測データブロックが記憶される時には、過去の観察に依存する原因推定（causal estimation）に対する制限は、もはや必要ではない。

特許出願PCT WO 2004/036753（2004年4月29日公開）は、ＢＰＳＫ又はＱＡＭ変調に基づいて送信チャネルを使って送信された観測データの位相推定プロセスを説明している。この操作は観測データブロックに対し実行されるが、少なくとも１つの位相ロックループを使って、前記ブロックから抽出された観測データの所定のシーケンスで実行される。

本書に記載のプロセスによれば、実質的に原因推定に対する上記制限を実質的に除去することができる。なぜならば、位相φ_ｋについての現在の評価をするとき、ブロックによる処理は同一ブロック内の以前の観測データのみならず、以後の観測データを考慮に入れることが可能であるからである。

よって、上記のプロセスの実施モードにおいて、この制限は、一方向とその逆方向における観測データを読み取ることにより、従来の位相ロックループを含む第１と、第２反復プロセスを使うことにより克服される。

しかしながら、上記プロセスは、特にターボデコーダを備える受信機だけでしか利用できない情報を使用するという不利がある。なお、その場合、シンボルが+１か、−１かに等しいBPSKタイプの変調の場合においても、観測データの信頼度の情報は使用できる。

この情報は、ターボデコーダにおいて、ソフト情報形式で(in the form of soft information)使用できる。即ち、各シンボルについてのアプリオリの情報である。

しかしながら、マトリックスの状態数が２に限定される場合、上記の複数状態と対応する複数シンボルに関する各観測データの各位相の誤差範囲は±Π/２である。ＢＰＳＫ変調で伝送される観測データの検出のために使われる先端技術を使った実際の位相ロックは正確に機能するが、ＢＰＳＫ変調の上記ソフト情報に基づく追加的訂正を導入することは、結局有用性を半減する。

具体的には、本発明は、特に、フレキシブルな誤差訂正システムを持つ受信機にＱＡＭ変調で伝送されるデジタル信号処理に対応するデジタル受信機のための位相推定プロセスを提案することを目的とする。又は、更に一般的に、所謂ターボ方法として反復プロセスを用いる全受信機に対応するものである。この方法は従来エラー訂正符号化（ターボコード）、等価（ターボ等価）または同期（ターボ同期）に使われていたものである。

本発明の別の目的は、自動ゲインコントロール（AGC）を有するデジタル受信機に対するゲインの推定プロセスを提案することである。これは特にソフト誤差訂正システムを備えた受信機にＱＡＭ変調で伝送されるデジタル信号処理に適応するものである。更に一般的に、本発明の位相推定プロセスに関係し前述した反復方法を使う全受信機に対応する。

本発明の別の目的は、デジタル受信機に対してデジタル信号処理に特に適応する、位相とゲインとを合わせて推定するプロセスを提案することである。特にＱＡＭ変調で伝送されるデジタル信号処理に適合するものである。更に一般的に、本発明の位相又はゲインの推定プロセスに関係し前述した反復方法を使う全受信機に対応する。

本発明の別の目的は、１キャリア及び/又はマルチキャリアの受信機において、本発明の位相とゲインとの共同推定プロセス（joint process for estimate of phase and gain）における位相推定プロセスとゲインの推定プロセスを実施することである。

最後に、本発明の別の目的は、特定の位相ロックループ構成を実施することである。この構成は、変化する誤差の訂正システムを備えるデジタル受信機、又は更に一般的に、等価な反復方法を使う受信機において、位相とゲインとの共同推定において位相とゲインのそれぞれの推定精度を高めるものである。

メモリに記憶された観測データパラメータの位相及び/又はゲインの推定プロセスは、ＱＡＭ変調で伝送チャネルを介して伝送されるビットの組で形成される一連のデジタルシンボルに対応するものであり、本発明の対象である。同推定プロセスは、一連の観測データから、少なくとも、位相及び/又はゲインのパラメータの反復推定を実行するステップを有している点で注目すべきである。この反復推定は特定の位相及び/又はゲインの関係に基づいて実行される。この関係は前記シーケンスにおける観測データの推定位相及び/又はゲインをリンクするものである。これら観測データから推定される、連続した位相及び/又はゲイン値の少なくとも１つに基づいて前記位相及び/又はゲインパラメータを推定する少なくとも１つの適応プロセスをイニシャライズし、尤度確率値（likelihood probability）に基づいて、これら位相及び/又はゲインパラメータを推定する少なくとも１つの機能を含む該適応型推定プロセスを実行する。なお、前記尤度確率値（likelihood probability）は、これらシンボルを構成するビットのフルセットに関して各観測データのlog尤度確率で表現されるものである。

本発明のプロセスは、「ターボ」タイプのデコード構成を持つデジタル信号受信機を使う応用分野に使われる。特に、多数状態を持つＱＡＭ変調を使う大量のデジタル信号のための受信機を使う応用分野に使われる。

従来例に関する図１を除いた以下の図面の説明を読むことにより、本発明はより深く理解される。

本発明のプロセスの詳細な説明は、図２ａと図２ｂとに関連して与えられる。

指摘すべきことは、一般的に、本発明のメモリに記憶されている観測データに対する位相及び/又はゲインパラメータを推定するプロセスが、伝送チャネルを使って伝送されるＱＡＭ変調における１組のビットにより形成される一連のデジタルシンボルに対応するデータに適用されることである。

“メモリに記憶されている観測データ”とは、メディア上のメモリに記憶されている観測データｙ_ｋの組を意味する。

特に、本発明のプロセスの有利な実施例において、本発明のプロセスはブロックとしてメモリに記憶されている観測データに対して実施される。特に、ターボデコード装置を備えた受信機により処理される観測データに対して実施される。これについては後述する。

一般的に、図２ａに関して指摘されていることは、本発明のプロセスは、ステップＡでは、メモリに記憶されている観測データから位相及び/又はゲインパラメータの反復推定を行うことを含む。

前記の反復推定は、選択された一連の観測データの推定位相に関連する特定の位相及び/又はゲイン関係をベースにして行われる。

反復推定ステップＡの後で、勿論、ステップＡから生じる複数の推定される位相及び/又はゲイン値を利用することができる。

ステップＡの後のステップＢは、ステップＡで得られた観測データから推定される連続した位相及び/又はゲイン値の少なくとも１つに基づいて前記位相及び/又はゲインパラメータを推定する適応プロセスをイニシャライズする(initialize)。

前記イニシャライズした後で、適応推定プロセスが実行される。そして、本発明が関係するプロセスの特に注目すべきアスペクトに対応して、これは、その尤度確率値に基づいて位相及び/又はゲインパラメータを推定する少なくとも１つの機能を含む。尤度確率値は対象とする変調におけるシンボルを構成するビットの組に関する各観測データに対してlog尤度で表現されている。

一般的に、図２ａで示される本発明のプロセスの実施例、及び、図２ａ以降の図面で示され、以下で説明する実施例において、以下の点が指摘される。

―プロセスが１つ又は複数の以前の事例で推定される同一変数のためのパラメータ値で、現在の変数のためのパラメータ値を評価し、特に推定することができる場合、プロセスは反復的であるといわれる。

―逆に、プロセスが、例えば外的物理法則に関係するパラメータ値における変化の推定又は評価に関し現在の変数のパラメータを評価し、特に推定するプロセスである時、プロセスは適応型であるといわれる。

―指摘すべきことは、本発明によるプロセスを実行する例において、ステップＡで実行される反復推定は少なくとも以前の観測データに対する位相とゲインパラメータの記憶された値を使って、ランクｋの現行観測データｙ_ｋに対する位相又はゲインパラメータ値を得る点である。

逆に、ステップＢで実施される適応型プロセスは位相とゲインパラメータの値の反復性の概念のみならず、外部変数を使う。この外部変数は外的に得られた尤度確率値に依存する位相及び/又はゲインパラメータの推定に対応する。

図２ａのステップＡの実施に関して、次の点が指摘される。即ち、特定の位相及び/又はゲイン関係は次の式を満足する。
（４） φ_ｋ＝φ_ｋ−１+γＦ（ｙ_ｋ,φ_ｋ−１）；Ｇ_ｋ＝Ｇ_ｋ−１＋γＧ（ｙ_ｋ,φ_ｋ−１）
上記関係式において、
φ_ｋ,φ_ｋ−１はランクｋ及びｋ−１の観測データy_ｋ、y_ｋ−１の推定位相値を表わす。
Ｇ_ｋ，G_ｋ−１はランクｋ及びｋ−１の観測データy_ｋ、y_ｋ−１の推定ゲイン値を表わす。
ＦとＧはそれぞれＱＡＭ変調のタイプに依存する特定関数を表わす。
γは予め決めたフィルター関数である。

逆に、図２ａのステップＢにおける、位相及び/又はゲインパラメータを推定するための適応型プロセスは、
ＡＥΦ（φ_ｋ，φ_ｋ−１,Ｌ^ｋ）、
ＡＥＧ（Ｇ_ｋ,Ｇ_ｋ−１，Ｌ^ｋ）
により表わされる。
前記位相及び/又はゲインパラメータは、対象のＱＡＭ変調に対するシンボルを構成するビットの組に関する各観測データに対するlog尤度Ｌ^ｋで表現される尤度確率値に依存する、少なくとも１つの位相及び/又はゲインパラメータを推定する関数を含んでいる。

図２ａに関連して、留意すべきは、ステップＢが現行データブロックに対して実行されることである。ステップＢが実行された後に、本発明が関係する観測データに対する位相及び/又はゲインパラメータの推定プロセスは、当然次のデータブロックの実行を行う。この操作は図２ａにおいて、ドットとしてリターン矢印で示されており、関係式Ｂ＝Ｂ+１を通って次のブロックに進む。

更に具体的には、シンボルを構成するビットの組に対する各観測データｙ_ｋの尤度で表現される尤度確率が考慮され、ステップＢで実行される適応型プロセスが、ＱＡＭ変調に対する全シンボルに関連してランクｋの観測データｙ_ｋに対する位相又はゲインを推定する反復関数を含まれることは好ましいことであると指摘されている。

図２ｂを参照すると、イニシャライズステップ(initialization stage)は、次式を検証することにより、適応型プロセス及び位相又はゲインを推定する反復関数のイニシャライズを含んでいる。
ＡＥΦ：φ_ｋ＝φ_ｋ−１+ＣＡｒｇ_ｋ（Im_ｋ,Ｗ_ｊ）
ＡＥＧ：Ｇ_ｋ＝Ｇ_ｋ−１+ＣＭ_ｋ（Re_ｋ,Ｗ_ｊ）

図２ｂのステップＢ２は当然、適応型プロセス及び特に、ＱＡＭ変調に対する全シンボルに関連してランクｋの観測データｙ_ｋに対する位相又はゲイン位相又はゲインを推定する反復関数のイニシャライズを含む。

一般に、図２ｂを参照すると、適応型プロセスを定義できる位相推定のための反復関数は上記式ＡＥΦを満たすことが指摘されている。

該反復関数は次を含む。
―ランクｋ−１の前の観測データｙ_ｋ−１に対する推定される位相φ_ｋ−１を含む、推定される位相角(phase argument),
―ＣArg_ｋ（Im_ｋ,Ｗ_ｊ）で表現される補正位相角（correcting phase argument）。この補正位相角は本発明によるプロセスの注目すべきアスペクトに応じて、（ＱＡＭ変調に対する全シンボルに関連して現行観測データに対する確率で表現される）重み値又は信頼値により重み付けられた現行観測データｙ_ｋに対して測定された位相角に比例するものである。

上記検討及び上記補正位相角の特定値を考慮すると、各観測データの位相を推定するための反復関数は次の関係式を満たす：

Ｗｊ（ｙ_ｋ，Ｌ^ｋ，φ_ｋ−１）は、現行観測データｙ_ｋに関するシンボルＱ_ｊに（attributed to）尤度確率で表現される重み値又は信頼値を表わす。

ゲイン推定の適応型に関しては、反復機能ＡＥＧは以下のものを含む。
―以前のする観測データｙ_ｋ−１に対する推定されるゲインＧ_ｋ−１、
―ランクｋの現行観測データｙ_ｋに関連するゲインと以前の観測データに対する推定ゲインＧ_ｋ−１との差値に比例するゲイン要素補正項ＣＭ_ｋ（Re_ｋ,Ｗ_ｊ）。
ランクｊの所与のシンボルＱ_ｊに対し、この差値はシンボルの組に関してシンボルＱ_ｊに寄与する尤度で重みまたは信頼度により重み付けされる。

前記補正ゲイン項の関数定義に対するコメントに関して、反復関数ＡＥＧは次の式を満たす。

Ｗｊ（ｙ_ｋ，Ｌ^ｋ，Ｇ_ｋ−１）は、現行観測データｙ_ｋに関するシンボルＱ_ｊに帰する（attributed to）尤度確率で表現される重み値又は信頼値を表わす。

当然、反復関数が位相推定のために実行される時には、観測データｙ_ｋに関してシンボルＱ_ｊに帰する尤度で表現される重み値又は信頼性値は推定位相に依存する。又、反対に、適応的プロセスの反復関数がゲイン推定のために実行される時には推定ゲインに依存する。

位相推定のために、現行観測データｙ_ｋに関するシンボルＱ_ｊに帰する確率で表現される重み又は信頼性値は次の関係式を満たす。

上記式において、expは指数関数を表わす。
ｑ^ｊ _ｍはＱＡＭ変調におけるシンボルＱの第ｍ番目ビットを表わす。
Ｌⁿ _mは、Ｎ+１シンボルの変調に対する第ｎ番ＱＡＭシンボルの第ｍビットの現行観測データのｌｏｇ尤度値を表わす。
Ｌ^ｎは全ビットに対するｌｏｇ尤度値のリストを表わす。Ｌｎ＝（Ｌｎ１，…ＬｎＮ）である。
σ_ｂ ^２は対象の伝送チャネルに対するノイズ電力を意味する。
θは対象の観測データに対する推定位相角を表わす。

同様に、ゲイン推定の場合、観測データに関するシンボルＱ_ｊに帰する尤度で表現された重み又は信頼性値は次の関係式を満たす。

上記関係式において、関係式（７）と同一のパラメータは、パラメータＧを除き、同一のパラメータを表わす。パラメータＧは観測データｙ_ｎに対する推定ゲイン値を意味する。

加えて、上記関係式（７）、（８）により導入される上記重み又は信頼性値は制限するものでない。実際、ｙ_ｋに対するシンボルＱ_ｊに帰する確率を使って表現される重み又は信頼性値はソフトデマッパー又はターボデコーダーによりシンボル毎に提供される各シンボルの全体の値に有利に対応する。

この変形例において、ｙ_ｋに対するシンボルＱ_ｊに帰する尤度を使って表現された重み又は信頼性値は位相についての次の関係式を満たす。

また、ゲインについて次の関係式を満たす。

上記の関係式において、Ｌn,jは、ソフトウェアデマッパー又はターボデコーダによるシンボル毎に提供される尤度比ＩｎP(a_j)/Ｐ（ref）を表わす。但し、Inはネピア対数で、Ｐ（a_j）は複素シンボルａ_ｊの確率、Ｐ（ref）はレファレンス値の確率を表わす。

図３ａとそれに続く図を使って、本発明の位相及び/又はゲイン推定プロセスを実行する種々のモードについて以下に説明する。

図３ａは、第１の制限されない例で、本発明のプロセスが使われ、位相パラメータの推定を行う、特にＱＡＭ変調の全てのタイプに適合する例である。

上記の図３ａを参照すると、プロセスはメモリに観測データｙ_ｋを複数ブロックとして記憶し、図２ａで説明したステップＡを実行する。留意すべき事は、このステップが、一連の連続する観測データに対する推定位相を結ぶ特定位相関係を基礎に位相パラメータの反復推定を実行することを含んでいることである。又、留意すべき事は、メモリに記憶するこのステップが、コンプレックス デモジュレータにより生成される出力信号、又はターボデコーディングモジュールに適合した受信機の別の構成要素からの出力に基づいて、Ｋ+１観測値ｙ₀乃至ｙ_ｋを記憶することを含むことである。前記応用において、本発明のプロセスは２００から２０００の観測データのブロックに有利に適用される。各ブロックを構成する観測データ数は、観測データの伝送で使用されるＱＡＭ変調の応用とタイプに応じて選択される。

一旦、これらが前記メモリに記憶されると、図２ａのステップＡは、メモリに記憶された対象の観測データブロックの観測データの所定のシーケンスで実行される。特に、最初であればどのようなシーケンスも構成することができる。例えば、観測データｙ_ｋのデータが受信される順序であるクロノロジカル シーケンスを構成することができる。前記の反復推定を実行するのに使われる関数は図２ａにおけるステップＡを実行するために明細書に記載された関数である。

図３ａを参照すると、ステップｏはステップｂ１を伴う。ステップｂ１は第１の適応型プロセス（ＡＥΦ１と呼ぶ）をイニシャライズを含む。最後の推定位相値のような値から第１の適応型反復プロセスの第１の値を決定する。

イニシャライズステップ（initialization stage）ｂ１は第１適応型プロセスＡＥΦ１を実行することを含むステップｂ２を伴う。この適応型プロセスは、尤度値に依存する位相パラメータを推定するための、少なくとも１つの関数を有する。前記尤度値はＱＡＭ変調コンステレーションにおけるシンボルの全構成ビットに関連して、各観測データに対するlog尤度で表現される。

尤度値、又はソフト情報(soft information)は外部変数(external variable)を構成する。外部変数により前記プロセスが、本明細書で前記した定義と共に適応的になる。この第１の前記の適応型プロセスは、ランクｋの観測データｙ_ｋを例えば前方向に読取ることにより、第１組の連続した中間の推定位相ずれ誤差値φ_０乃至φ_Ｎを生成することができる。

ステップｂ１において実行されるイニシャライズは、第１適応型プロセスに対する第１値を固定することを可能にする。特に位相の場合において、推定されるパラメータが１つの観測データブロックから他の観測データブロックへの連続性を持つとき、第１適応型プロセスが、例えば、前のブロックに対する最後の推定値を考慮して、有利にイニシャライズされることが好ましい。特に、通常の伝送チャネルについては、観測データのブロックを伝送する期間、位相パラメータは伝送の安定性のために緩慢に変化するパラメータであることが分かる。

ステップｂ２において前記第１適応型プロセスを実行することにより、図３ｃの上の左から右に向かう矢印で示されるように、一シーケンスの推定位相値φ₀…φ_ｋ…φ_ｋを構成できる。

ステップｂ２の後で、ステップｂ３の本発明のプロセスが続く。ステップｂ３では第２適応型プロセスＡＥΦ_２をイニシャライズする。このプロセスでは第１適応型プロセスＡＥΦ_１の実行の結果得られる最終中間的位相ずれ誤差値から、第２適応型プロセスＡＥΦ_２の第１値を決定する。

第２適応型プロセスＡＥΦ_２の第１値、即ち、φ′_ｋが第１適応型プロセスＡＥΦ_１により算出される最後の数値φ_ｋでイニシャライズされることが好ましい。この操作は図３ｃの下の右から左への矢印で示される。

ステップｂ３の後、ステップｂ４が続く。ステップｂ４では第２の適応型プロセスを実行する。これは当然、シンボルを構成するビットの組に関する各観測データに対するlog尤度で表現される尤度値に依存する位相パラメータを推定するための関数を含む。第２の適応型プロセスＡＥΦ_２により、推定された中間的位相ずれ誤差値φ′_ｋ-１乃至φ′₀の第２の組を生成することができる。

次に本発明のプロセスでは、ステップｂ５において、下記の関係式に従って、中間的位相ずれ誤差値のランクｋの第１と第２の組合せとして、ランクｋの各観測データｙ_ｋに対する推定位相ずれ誤差φ″_ｋに対する最終値を算出する。
（９） φ″_ｋ＝ｇ（φ_ｋ、φ′_ｋ）

一般的に次のことに留意すべきである。即ち、中間的位相ずれ誤差値の第１と第２の組を組合せる関係式はＱＡＭ変調のタイプに関係して選択される関数である。

ある実施例において、例えば、ｇは、第１と第２の連続的推定中間的位相ずれ誤差値のリニア結合として、最終推定位相ずれ誤差値を表現するように選択される。

特定の選択では、リニア結合係数Ａ_ｋ＝Ｂ_ｋ＝１／２が選択され、リニア結合は次の式で表わされる。
（１０） φ″_ｋ＝Ａ_ｋφ_ｋ+Ｂｋφ′_ｋ

更に、観測データのランクｋをベースに２つの適応型プロセスの内の１つを優先するために、リニア結合係数Ａ_ｋ、Ｂ_ｋは可変係数であってよい。従って、図３ｃで説明するブロックの右部分において、第１の適応型プロセスＡＥΦ１を優先してリニア結合のための重み付けを選択することができる。又、逆に、ブロックの左側部分において第２の適応型プロセスにより重みを付けることもできる。この操作モードは反復を行った適応型プロセスを優先することを可能にし、位相算出時に、より高い精度を要求することを可能にする。

第１と第２の適応型プロセスＡＥΦ１とＡＥΦ２の実施に関して留意すべきことは次の点である。関係式１１を満足する第１と第２の位相ループを使って有利なやり方で、これら（第１と第２の適応型プロセスＡＥΦ１とＡＥΦ２）が使用されてもよい。

特に、前記関係式については、第１の適応型プロセスは、前記関係式においてε＝-１で使用され、第２の適応型プロセスは、ε＝+１で使用される。

第１と第２の位相ループを実装し、第１と第２の適応型プロセスを連続して実行する時、観測データｙ_ｋに関してシンボルＱ_ｊに帰するlog尤度で表現される重み付け又は信頼性値の表現は、前記の関係式７を満たす。

図３aを使って説明された実施例モードにおける本発明の位相推定プロセスは、全タイプのＱＡＭ変調に適しており、又、例えば、ターボデコーダーにより提供されるシンボルのソフト情報（soft information）のような情報を利用する。

実際問題として次のことに留意すべきである。即ち、位相ループを表わす適応型関数を表わす式におけるパラメータγは、補正する位相ずれ誤差モデルに基づいて位相又はゲイン角（argument）の補正項を形成するために、それぞれ位相角又はゲイン要素に適用するデジタルフィルタ関数を含んでもよい。単純な位相ずれ誤差に対しては単純な比例補正を使うことが出来るが、更に複雑なケースでは積分補正、又は、より高次のフィルタを使うと有利である。フィルタ関数γは、関係式１２を満たす２次のデジタルフィルタを使って実行すると有利である。

上記関係式において、ｚはｚ変換を意味する。

観測データｙ_ｋを受信する受信機のゲインを推定する本発明のプロセスを実施する第２の例は図３ｂに関連して説明する。この受信機は、例えば、オートゲインコントロールループを備えるものである。

一般則として、ＱＡＭタイプ変調に対して、受信機においてＱＡＭシンボルの正確なデマッピング処理が行えるように、チャネルゲインを推定することは必要である。

このような状況で、受信された観測データｙ_ｋは下記関係式１３によるシンボルａ_ｋと関係していると仮定する。
（１３） ｙ_ｋ＝Ｇａ_ｋ+ｂ_ｋ、ｋ∈[0、Ｋ]

上記関係式において、ａ_ｋはＱＡＭコンステレーションのシンボルを意味し、Ｇは伝送チャネルにより提供されるゲイン、又は減衰を表わす。

図３ａに関係する実施例において位相推定だけの場合と同様に、ゲイン推定だけの実施例は、メモリにおける複数観測データブロックｙ_ｋを記憶するステップｏ及び図２ａのステップＡを使うことを含む。

図２ａのステップＡを実行する場合において、当然、一連の前記の観測データにおける連続した観測データに対する推定ゲインに関する、特定ゲイン関係式の反復関数を使う。

ステップｏの後、ゲインパラメータ推定のための本発明のプロセスはステップｃ１を呼び出すこと（calling）を含む。ステップｃ１は、第１適応型プロセスをイニシャライズし、第１適応型プロセスの第１値、例えば、最後の推定ゲイン値を決定する。第１の適応型プロセスは図３ｂではＡＥＧ１と表記される。前記イニシャライズは、位相のみを推定するためのプロセスを実行するときに行われるものと同様の条件の下で、行われてもよい。特に、補正がしたい観測データｙ_ｋの大きさの一時的変化のためのノーマルモデル（noramal model）に関係して要素γを算出するためのデジタルフィルタの使用に関しては、前記イニシャライズは行われてもよい。位相だけを推定する場合と同様にして、単純な比例補正器を単純なモデルに対して使っても良い。これに対し、更に複雑な場合は、積分補正、又は、より高次のフィルタを使うことができる。ゲインを推定する適応型プロセスのイニシャライズは第１ゲイン値Ｇ₀を提供することを含み、位相だけを推定するための本発明のプロセスを実行する例において前記したように、フィルタγは、２次のデジタルフィルタを使って提供される。

前記したステップｃ１の後、ステップｃ２が続く。ステップｃ２では第１適応型プロセスＡＥＧ１を実行する。これは当然、ＱＡＭ変調コンステレーションのシンボルの構成ビットの組に関する各観測データに対し、log尤度で表現される尤度値に依存するゲインパラメータ推定関数を少なくとも１つ含む。第１適応型プロセスＡＥＧ１は、図３ｃの例についてランクｋの観測データｙ_ｋを順方向に読むことによって、Ｇ₀乃至Ｇ_ｋで表される連続的中間的ゲイン値（successive intermediate gain value）の第１の組を生成できる。

位相だけを推定する第１の例と同様なやり方で、ゲインだけを推定するための本発明のステップｃ３におけるプロセスは、ＡＥＧ２と呼ばれる第２適応型プロセスをイニシャライズすることを含む。このプロセスでは第１適応型プロセスの結果得られた最終推定ゲイン値Ｇｋに基づいて第２適応型プロセスの第１値を決定する。

前記のイニシャライズステップの後、ステップｃ４が続く。ステップｃ４では第２適応型プロセスを実行する。これは当然、ＱＡＭ変調のシンボルを構成するビットの組に関する各観測データに対し、log尤度で表現される尤度値に依存するゲインパラメータを推定するための関数を少なくとも１つ含む。

第２適応型プロセスＡＥＧ２を実行すると、逆方向にランクｋの観測データｙ_ｋを読むことにより、連続的中間的ゲイン値Ｇ′_ｋ-1乃至Ｇ′₀の第２の組を生成できる。

位相だけを推定する本発明のプロセスの実施例と同様なやり方で、ステップｃ５が呼び出されて（is called）、ランクｋの全ての観察データｙ_ｋに対して最終推定ゲイン値Ｇ"ｋが算出される。この最終ゲイン値は、関係式１４によりランクｋの第１と第２の中間ゲイン値の組合せにより表現される。
（１４） Ｇ″_ｋ＝ｇ（Ｇ_ｋ，Ｇ′_ｋ）

勿論、ステップｃ６によりブロック プロセシングの最終を処理することができ、処理は図３ａの場合と同様なやり方で、各ブロックに対し戻り時にＢ＝Ｂ+１でリピートしてもよい。

位相だけを推定するプロセスと同様なやり方で、ゲインだけを推定する実施例における本発明のプロセスは、有利に行われ、次の関係式１５を満足する第１と第２のゲインループを経て、第１と第２の適応型プロセスＡＥＧ１、ＡＥＧ２が実行される。

上記関係式において、εは、第１適応プロセスＡＥＧ１を実行する場合-１を取り、第２適応プロセスＡＥＧ２を実行する場合１を取る。

本発明の、位相とゲインの推定プロセスの非制限的な実施例（第１と第２の）について留意すべき点は以下の通りである。即ち、本発明のプロセスの実行を可能にする関係(relationship)における位相とゲインのパラメータの完全な対称性は、位相及び/又はゲインの変数（前記適応型プロセスを表わす関数における位相とゲイン値の表現を決定している）の独立性に関係している。

上記の図３ａ及び３ｂに対応する位相だけ又はゲインだけを推定する実施例を実施するために、他のパラメータの正確な知識を持つことが望ましい、つまり、位相に対するゲインの知識、又はその逆である。

前記パラメータの１又は他の知識に関連する問題は、実際に同時に相互に存在するのであるが、各推定は通常、別の結果に基づかなければならないことである。複数パラメータの内の１つが正確に推定されなかった場合、伝送誤差のリスクは避けられない。

本発明のプロセスを実行する第３の実施例によると、以下の条件で上記伝送誤差のリスクを解決できる。この第３の実施例はこれら２つのパラメータを組合せて推定することを含んでいる。

以下、本発明の要旨の位相とゲインの推定に使われる組合せ推定によるプロセスを図４ａに関連して説明する。

上記の図４ａを参照すると、本発明によるプロセスは、観測データｙ_ｋに基づいてゲイン及び/又は位相パラメータの反復的推定を行うステップｏ′を実行することを含んでいる。この反復的推定はゲインー位相ループを使って実行され、対象のＱＡＭ変調のためのシンボルに関する各観測データのゲインと位相の推定を可能にする。

特に、前記のステップｏ′は観測データｙ_ｋのブロックを記憶することを含む。この記憶には、図２ａにおけるステップＡに実質的に対応するステップａ′の実行が伴う。特にステップＡ′はステップＡにおける反復関数及び同ステップにおけるゲインに関連する反復関数の実行に対応する。

ステップＡ′を実施するために、又、好ましい非制限的な実施例に従って、位相とゲインとの２つの反復関数の使用は、適応型位相又はゲインプロセスの呼び出し(calling)と、有利に置換される。それら関数において、全ビットに対するlog尤度値のリストＬ^ｋの値は、任意に０と等しくされる。これら条件の下に、ゲインと位相を推定するための適応型プロセスが各シンボルに関して実行される。ＱＡＭ変調に対する各シンボルは等しく起こりえる（equally likely）と考えられる。この仮定はイニシャライズステップのみに対し容認できる精度を十分保証するものである。

前記ステップｏ′の後、ゲインと位相のパラメータの共同推定（joint estimation）のための本発明のプロセスはステップｄ１を呼び出すことを含む。このプロセスについては図４ａで説明する。ステップｄ１では、第１適応型プロセスがイニシャライズされ、推定ゲイン値Ｇ₀と位相値φ₀が決定される。

ゲインと位相の共同推定を実行するために留意すべきことは、本発明のプロセスは次の関係式１６を満たす観測データｙ_ｋを処理するステップを含む点である。
（１６）ｙ_ｋ＝ａ_ｋＧｅ^ｉθｋ+ｂ_ｋ ｋ∈[0、Ｋ]

この関係式において、ａ_ｋはＱＡＭシンボルを表わし。Ｇは、例えば、伝送チャネルにより提供される寄生的ゲインを表す。θｋは処理されるべき寄生的位相ずれ誤差を意味する。特に第１適応型処理ＡＥＧΦ１を実行するために、ステップｄ１は第１推定ゲイン値Ｇ₀及び位相値φ₀をイニシャライズし、決定することができる。

上記イニシャライズ ステップの後、ステップｄ２が続く。ステップｄ２では第１適応型処理ＡＥＧΦ１を実行する。イニシャライズ ステップは、ＱＡＭ変調のシンボルに対する構成ビットの組に対し各観測データの尤度確率（log尤度で表現される）に依存するゲインと位相のパラメータを推定する少なくとも１つの関数を含む。

第１適応型処理ＡＥＧΦ１を実行することにより、ゲインＧ₀乃至Ｇ_ｋ及びφ₀乃至φ_ｋに対し、それぞれ連続的中間値の第１の組を生成することができる。

ステップｄ２の後にステップｄ３が続く。ステップｄ３では第２適応型ゲイン・位相処理ＡＥＧΦ２のイニシャライズを実行し、第１適応型処理ＡＥＧΦ１の実行の結果得られる最終推定ゲインと位相値に基づいて第２適応型プロセスに対する第１値を決定する。

その後、第２適応型プロセスＡＥＧΦ２が実行される。この適応型プロセスは、ＱＡＭ変調のシンボルに対する構成ビットの組に対し各観測データの尤度確率（log尤度で表現される）に依存する、少なくとも１つの関数を含む。

第２適応型プロセスＡＥＧΦ２を実行して、逆方向にランクｋの観測データｙ_kを読むことにより連続する中間的ゲイン・位相値Ｇ′_ｋ+１乃至Ｇ′₀及びφ′_ｋ+１乃至φ′₀の第２の組を生成することができる。

上記ステップｄ４の後にステップｄ５が続く。ステップｄ５では、関係式１７によるランクｋの第１と第２の中間的ゲインと位相値の組み合わせとして、ランクｋの観測データｙ_ｋに対し、それぞれ最終ゲインと位相を算出する。
（１７） Ｇ″_ｋ＝ｇ（Ｇ_ｋ、Ｇ′_ｋ）
φ″_ｋ＝ｇ（φ_ｋ、φ′_ｋ）
上記式において、Ｇは特定の関数（specific function）を表わす。

図４に関して説明した位相とゲインの組合せの推定を実行するために、第１と第２の適応型プロセスＡＥＧΦ１とＡＥＧΦ２が、次の関係式１８、１９を満たすゲイン−位相ループを使って実行される。

上記の関係式において、留意すべきことは、位相又はゲインだけをそれぞれ推定するときには、パラメータεは第１適応型プロセスＡＥＧΦ１に対しては-１に等しくし、第２適応型プロセスＡＥＧΦ２に対しては+１に等しくすることである。

更にパラメータγ１とγ２は、ＱＡＭ変調タイプに従って選択された特定のフィルタ関数を表わす。位相だけ、ゲインだけの推定を実行する例においてパラメータγを選択するのと同様なやり方で選択してよい。

位相とゲインの共同推定を実行するために、観測データに関するシンボルＱｊに帰するlog尤度で表現される重み又は信頼性値は、位相及びゲインパラメータについて同時に算出されねばならない。

このような仮定の上に、上記重み又は信頼性値は次の関係式（２０）を満たす。

上記関係式において、Ｇは対象の観測データｙ_ｎに対する推定ゲインを表わし、θは対象の観測データｙ_ｎに対する推定位相を表わす。

図３ａ、３ｂと４ａで説明した本発明のプロセスが実施される全ての場合において、適応型プロセスの実行ステップｂ２、ｂ４；ｃ２、ｃ４及びｄ２、ｄ４は反復プロセスｂ′２、ｂ′４、ｃ′２、ｃ′４及びｄ′２、ｄ′４の数だけ繰り返してもよい。これらは計算の精度を高めるためにループＲ＝Ｒ+１、Ｒ′＝Ｒ′+１によりシンボル化されるものであり、反復回数は３、４に達する。

本発明によるゲイン位相ループを、図４ｂに関連して以下に説明する。

上記ゲイン位相ループは加算器４０を有し、その入力端で、以前の観測データと位相補正項CArgk（Im、Wｊ）、又、ゲインCMｋ（Rek、Wj）のために推定される、位相パラメータφ_k＋εとゲインパラメータＧ_k＋εを受け取り、ランクｋの現行観測データγ_ｋのために推定される位相パラメータγ_ｋ及び/又はゲインパラメータＧ_kを生成する。位相ＦとゲインＧの角（argument）の機能モジューラス（functional modulus）４１が備えられ、以前の観測データに対し推定されるランクｋの現行観測データｙ_ｋ、位相とゲインのパラメータφ_k＋ε、Ｇ_k＋ε、ＱＡＭシンボルのアルファベットからの複数シンボルＱ１乃至ＱＭ、及び全ビットＬ^ｋに対するlog尤度値のリストを受信すると、機能モジューラス４１は位相及び/又はゲイン角を提供する。位相及び/又はゲイン角に適用されるフィルタ４２（γ１、γ２）は位相及び/又はゲイン補正項を加算器４０に渡す。図４ｂにおける位相ループはソフトウェアモジュール又は専用計算器を有しても良い。

本発明のプロセスは、特に位相とゲインとを共に推定する場合、単一キャリア、マルチキャリア伝送のときに実施されると有利である。

特に理解しておくことは、精度レベルのために、位相とゲインとを共に推定する場合、このプロセスは広範囲の周波数レンジに適用することができることである。特に、マルチキャリアの場合そうである。実際に、伝送周波数に依存する伝送チャネルの変化性（variability）は搬送周波数に関して同一ではない。本発明のプロセスの実施のフレキシビリティは、伝送チャネルとマルチキャリア伝送条件の変化性にも拘らず、結果の精度が高いために特に魅力的である。

特に、マルチキャリア伝送の場合、本発明のプロセスは、マルチキャリアシステムを形成する少数のサブキャリア上で独立に実施してもよい。最終的に、独立性の概念は、サブキャリア周波数値に基づいて、例えばフィルタパラメータγのようなパラメータを限定する（qualify）ことになる。

更に、このような場合、本発明のプロセスは、マルチキャリアシステムの他のサブキャリアに対し、対象サブキャリア周波数値に応じてゲイン及び/又は位相値を補間してよい。勿論、このような補間により、最適精度を達成するようにゲイン及び/又は位相値を調整することができる。

図４ａに示す本発明のゲインと位相の共同推定のプロセスは、図５ａで説明する受信機システムにおいて実施される。即ち、観測データｙ_ｋを提供する後者のコンプレックスデモジュレータからの出力において実施される。

前記受信機に関して、コンプレックス デモジュレータからの出力において（図５ａでは図示していないが）、後者はゲイン・位相処理モジュール１を有している。同モジュールを通して、ランクｋの観測データｙ_ｋに対する推定位相ずれ誤差値φ_ｋが推定ゲイン値Ｇ_ｋと共に与えられる。推定ゲイン値Ｇ_ｋは受信機のＡＧＣループ（自動ゲインコントロールループ）で与えられる。ゲイン・位相処理モジュール１にはフレキシブル デマッパー２が接続される。デマッパー２にはターボデコーダ３が接続される。ターボデコーダ３はフレキシブル情報、即ち、Ｌ^ｋを提供する。又、Ｌ^kは、全ビットに対するlog尤度値のリストを意味するものであって、観測データに関して、シンボルＱjに帰するlog尤度で表現される重み又は信頼性値を算出する。この重み又は信頼性値は、関係式７及び８に従って、位相及びゲインに対して提供される。これら式については前述したとおりである。次に、前記Ｌ^ｋ値は、ＡＥＧΦ関数を実行するゲインー位相ループに渡される。反復プロセスはゲイン位相ループで、符号４で表わされており、対象の観測データｙ_ｋに関する推定される位相とゲインのパラメータφ_ｋとＧkを位相処理モジュール１に与える。

勿論、位相とゲインの組合せ検出に基づいて本発明のプロセスを実行し、ｏ′ステップをイニシャライズするためには、反復イニシャライズ プロセス(iterative initialization procedure)がモジュール４を介して実行される。そのモジュール４において、Ｌ^ｋ＝０がｋに関係なくセットされる。次にモジュール４はイニシャライズ値（initialization values）を単なる反復プロセスＡＥＧΦ（Ｌ^ｋ＝０）から提供する。次に、２つの適応型プロセスＡＥＧΦ１とＡＥＧΦ２がそれぞれイニシャライズされ、実行される。

振幅・位相処理モジュール１、フレキシブル デマッパー２、ターボデコーダー３は、この技術分野では公知であるから、これらのモジュールについては詳しく説明しない。

ゲインー位相ループモジュール４はデジタル処理モジュールで、関係式１８、１９及び２０を満たすゲインー位相ループを実行する関数を有している。これら関係式は図４ｂに関係して既に説明した。

最後に、図５ｂについて説明する。ＱＡＭシンボルがインターリーブされるように、同シンボルが伝送ビットシーケンスのインターリービングで伝送される場合には、本発明の受信機は、ターボデコータ３の上流に配置されるデインターリーバモジュール３ａ及び、ターボデコータ３からの出力（即ち、ゲインー位相ループ４の上流）に配置されるインターリーバモジュール３ｂを有している。

上記複数インターリーバの操作プロセスとインターリーバモジュールについては、この技術分野では公知であるから、詳しく説明しない。観測データとシンボルについて、ｑπ/２（ｑ：整数）の位相の曖昧さを除去することができる。

従来技術のＢＰＳＫの説明図である。 例示的に示す、本発明のプロセスを実行する際の本質的ステップのフロー図である。 例示的に示す、〔図２ａ〕で示した本発明のプロセスによる適応型プロセスの実行の後のイニシャライズステップの実装を示すフロー図である。 例示的に示す、位相推定に関する第１の例における本発明のプロセスを実装する際の重要ステップのフローチャートである。 例示的に示す、ゲイン推定に関する第２の例における本発明のプロセスを実装する際の重要ステップのフローチャートである。 例示的に示す、図３ａと３ｂに示される本発明のプロセスを実装するための順方向及び逆方向に観測データを読取るクロノグラム（chronogram）を示す。 例示的に示す、ゲインと位相の組合せ推定に関する第３の好ましい例における本発明のプロセスを実装する際の重要ステップのフローチャートである。 本発明の対象による位相ゲインループである。 〔図４ａ〕で説明した発明によるプロセスを実施する本発明の対象に応じて、ゲインー位相ループを備えたターボデコード受信機を機能ブロックを使って表わした概略図である。 位相の明確化を行うことのできるターボデコード受信機を機能ブロックを使って表わした図である。

Claims (20)

1. 伝送チャネルにより伝送されるＱＡＭ変調におけるビットの組により形成される連続するデジタルシンボルに対応する、メモリに記憶される観測データの位相及び/又はゲインパラメータを推定するための方法であって、以下のステップを有する方法。
   ａ）一シーケンスの観測データに基づいて、前記位相及び/又はゲインパラメータの反復推定を行うステップであって、前記反復推定は前記シーケンスにおける連続した観測データの推定位相をリンクする特定の位相及び/又はゲイン関係を使って得られることを特徴とするステップ。
   ｂ）前記観測データから推定される、連続した位相及び/又はゲイン値の少なくとも１つに基づいて、前記位相及び/又はゲインパラメータを推定する少なくとも１つの適応型プロセスのイニシエーションを行い、前記シンボルを構成するビットのセットに関し、各観測データのlog尤度確率で表現される尤度確率値（likelihood probability）に依存する前記位相及び/又はゲインパラメータを推定する機能を含む少なくとも１つの適応型推定プロセスを実行するステップ。
2. 前記特定の位相及び/又はゲイン関係が次の式を満足する請求項１に記載の方法。
   φ_ｋ＝φ_ｋ−１+γＦ（ｙ_ｋ,φ_ｋ−１）；
   Ｇ_ｋ＝Ｇ_ｋ−１＋γＧ（ｙ_ｋ,Ｇ_ｋ−１）
   但し、
   φ_ｋ,φ_ｋ−１は、ランクｋ及びｋ−１の観測データy_ｋ、y_ｋ−１の推定位相値を表わ す。
   Ｇ_ｋ，Ｇ_ｋ−１は、ランクｋ及びｋ−１の観測データy_ｋ、y_ｋ−１の推定ゲイン値を表 わす。
   ＦとＧは、それぞれＱＡＭ変調のタイプに依存する特定関数を表わす。
   γは、予め決めたフィルタ関数を表わす。
3. 前記適応型方法は、シンボルを構成するビットの組に対する各観測データのlog尤度で表現される尤度確率が考慮され、対象となるＱＡＭ変調の全シンボルに関するランクｋの各観測データｙ_ｋの位相又はゲインを推定する反復関数を含むことを特徴とする請求項１又は２に記載の方法。
4. 各観測データの位相を推定するための反復関数は次の関係式を満たすことを特徴とする請求項３に記載の方法。
   

   

   
5. 各観測データのゲインを推定するための前記反復関数は次の式を満たすことを特徴とする請求項３又は４に記載の方法。
   

   

   

   Ｗｊ（ｙ_ｋ，Ｌ^ｋ，Ｇ_ｋ−１）は、現行観測データｙ_ｋに関するシンボルＱ_ｊに帰する尤度確率で表現される重み値又は信頼値を表わす。
6. 位相推定時に、現行観測データに関するシンボルＱｊに帰する尤度で表現される重み又は信頼性値は次の関係式を満たす請求項１乃至５の何れか一項記載の方法。
   

   

   上式において、
   expは指数関数を表わす。
   ｑ^ｊ _ｍはＱＡＭ変調におけるシンボルＱ_ｊの第ｍ番目ビットを表わす。
   Ｌⁿ _mは、Ｎ+１シンボルの変調に対する第ｎ番ＱＡＭシンボルの第ｍビットの現行観測データのｌｏｇ尤度値を表わす。
   Ｌ^ｎは、全ビットに対するｌｏｇ尤度値のリストを表わす。Ｌｎ＝（Ｌ^ｎ _１，…Ｌ^ｎ _Ｎ）である。
   σ_ｂ ^２は、対象の伝送チャネルのノイズ電力を意味する。
   θは、推定位相角を表わす。
7. ゲイン推定時に、観測データに関するシンボルＱｊに帰する尤度で表現された重み又は信頼性値が次の関係式を満たすことを特徴とする請求項１乃至６の何れか一項に記載の方法。
   

   

   但し、
   expは指数関数を表わす。
   ｑ^ｊ _ｍはＱＡＭ変調におけるシンボルＱの第ｍ番目ビットを表わす。
   Ｌⁿ _mは、Ｎ+１シンボルの変調に対する第ｎ番ＱＡＭシンボルの第ｍビットの現行観測データのｌｏｇ尤度値を表わす。
   Ｌ^ｎは、全ビットに対するｌｏｇ尤度値のリストを表わす。Ｌｎ＝（Ｌ^ｎ _１，…Ｌ^ｎ _Ｎ）である。
   σ_ｂ ^２は、対象の伝送チャネルに対するノイズ電力を意味する。
   Ｇは、推定ゲイン値を表わす。
8. 位相パラメータを推定するために、前記シーケンスにおける連続する観測データの推定位相にリンクする特定の位相関係に基づいて、以下のステップを含む位相パラメータの反復推定を行うステップｏの後に、次のステップを含むことを特徴とする請求項１乃至４の何れか一項に記載の方法。
   -ｂ１）第１適応型プロセスをイニシャライズし、第１適応型プロセスの複数の第１値を、例えば、最後の推定位相値を決定するステップ、
   -ｂ２）第１適応型プロセスを実行するステップであって、前記第１適応型プロセスは、ランクｋの観測データｙ_ｋを順方向に読むことにより、連続する中間的推定位相ずれ誤差値φ₀乃至φ_ｋの第１の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記位相パラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ、
   -ｂ３）第２適応型プロセスをイニシャライズし、前記第１の適応型プロセスの実行の結果得られた最終中間的位相ズレ誤差値に基づいて第２適応型プロセスの第１値（複数）を決定するステップ、
   -ｂ４）第２適応型プロセスを実行するステップであって、第２適応型プロセスはランクｋの観測データｙ_ｋを逆方向に読むことにより、連続する中間的推定位相ずれ誤差値φ′_ｋ-1乃至φ′₀の第２の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記位相パラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ。
   -ｂ５）下記の関係式に従って、ランクｋに対する第１と第２の位相ずれ誤差値組合せとして、ランクｋの各観測データｙ_ｋに対する最終推定位相ずれ誤差φ″_ｋを算出するステップ。
   φ″_ｋ＝ｇ（φ_ｋ、φ′_ｋ）
9. 第１と第２の適応型プロセスが、次式を満足する第１と第２の位相ループを介して実行されることを特徴とする請求項４と８に記載の方法。
   

   

   但し、第１適応型プロセスに対してはε＝-１、第２適応型プロセスに対しては、ε＝+１で使用される。
10. ゲインパラメータを推定するために、前記シーケンスにおける連続する観測データの推定ゲインにリンクする特定のゲイン関係に基づいて、以下のステップを含む位相パラメータの反復推定を行うステップｏの後に、次のステップを含むを特徴とする請求項１乃至５の何れか一項に記載の方法。
    -ｃ１）第１適応型プロセスをイニシャライズし、第１適応型プロセスの複数の第１値を、例えば、最後の推定ゲイン値を決定するステップ、
    -ｃ２）前記第１適応型プロセスを実行するステップであって、前記第１適応型プロセスは、ランクｋの観測データｙ_ｋを順方向に読むことにより、連続する中間的推定ゲイン誤差値Ｇ₀乃至Ｇ_ｋの第１の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記ゲインパラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ、
    -ｃ３）第２適応型プロセスをイニシャライズし、前記第１の適応型プロセスの実行の結果得られた最終推定ゲイン値に基づいて第２適応型プロセスの第１値（複数）を決定するステップ、
    -ｃ４）第２適応型プロセスを実行するステップであって、第２適応型プロセスはランクｋの観測データｙ_ｋを逆方向に読むことにより、連続する中間的ゲイン値Ｇ′_ｋ-1乃至Ｇ′₀の第２の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記ゲインパラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ、
    -ｃ５）下記の関係式に従って、ランクｋに対する第１と第２のゲイン値組合せとして、ランクｋの各観測データｙ_ｋに対する最終推定ゲイン値Ｇ″_ｋを算出するステップ。
    Ｇ″_ｋ＝ｇ（Ｇ_ｋ，Ｇ′_ｋ）
11. 第１と第２の適応型プロセスが、次式を満足する第１と第２の位相ループを介して実行されることを特徴とする請求項５と１０に記載の方法。
    

    

    但し、第１適応型プロセスに対してはε＝-１、第２適応型プロセスに対しては、ε＝+１で使用される。
12. ゲインと位相パラメータの組合せ推定のために次のステップ含むことを特徴とする請求項１乃至７の何れか一項に記載の方法であって方法。
    -ｏ′）観測データの１シーケンスから前記ゲイン及び/又は位相のパラメータの反復推定を行うステップであって、前記反復推定はゲインー位相ループを使って実行され、前記反復推定ステップはＱＡＭ変調シンボルの各観測データのゲインと位相を推定することを可能にすることを特徴とするステップ、
    ｄ１）第１適応型プロセスをイニシャライズし、第１推定ゲインと位相値Ｇ₀とφ₀を決定するステップ、
    ｄ２）前記第１適応型プロセスを実行するステップであって、前記第１適応型プロセスは、ランクｋの観測データｙ_ｋを順方向に読むことにより、連続する中間的推定ゲイン誤差値Ｇ₀乃至Ｇ_ｋの第１の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記ゲイン位相パラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ、
    ｄ３）第２適応型プロセスをイニシャライズし、前記第１の適応型プロセスの実行の結果得られた最終推定ゲインと位相の値に基づいて第２適応型プロセスの第１値（複数）を決定するステップ、
    ｄ４）第２適応型プロセスを実行するステップであって、第２適応型プロセスはランクｋの観測データｙｋを逆方向に読むことにより、連続する中間的ゲイン値Ｇ′_ｋ-1乃至Ｇ′₀の第２の組を生成するために前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現される尤度確率に依存する前記ゲインと位相のパラメータを推定する少なくとも１つの関数を含むことを特徴とするステップ、
    ｄ５）下記の関係式に従って、ランクｋに対する第１と第２のゲインと位相値の組合せとして、ランクｋの各観測データｙ_ｋに対する最終ゲインと位相の値を算出するステップ。
    Ｇ″ｋ＝ｇ（Ｇｋ，Ｇ′ｋ）
    φ″ｋ＝ｇ（φｋ、φ′ｋ）
    なお、ｇは特定関数を意味する。
13. 第１と第２の適応型プロセスが、次の関係式を満たすゲイン−位相ループを使って実行される請求項１２に記載の方法。
    

    

    

    但し、第１適応型プロセスに対してはε＝-１、第２適応型プロセスに対してはε＝+１であり、
    γ_１とγ_２は、ＱＡＭ変調のタイプに基づいて、予め決めたフィルタ関数である。
14. ゲインと位相の組合せ推定のために、観測データに関するシンボルＱｊに帰する尤度で表現された重み又は信頼性値は次の関係式を満たす請求項１２又は１３に記載の方法。
    

    

    但し、
    expは指数関数を表わす。
    ｑ^ｊ _ｍはＱＡＭ変調におけるシンボルＱの第ｍ番目ビットを表わす。
    Ｌⁿ _mは、Ｎ+１シンボルの変調に対する第ｎ番ＱＡＭシンボルの第ｍビットの現行観測データのｌｏｇ尤度値を表わす。
    Ｌ^ｎは全ビットに対するｌｏｇ尤度値のリストを表わす。Ｌｎ＝（Ｌ^ｎ _１，…Ｌ^ｎ _Ｎ）である。
    σ_ｂ ^２は対象の伝送チャネルに対するノイズ電力を意味する。
    Ｇは推定ゲイン値を表わす。
15. 適応型プロセスの各ステップが特定数の反復を行うことを特徴とする請求項８乃至１４の何れか一項に記載の方法。
16. 伝送チャネルを通ってＱＡＭ変調で伝送されるシンボルの組に関する以前の観測データのゲイン及び/又は位相に関する、現行観測データのゲイン及び/又は位相の適応型推定のためのゲイン位相ループであって、
    該ゲイン位相ループは、前記シンボルを構成するビットの組に関する各観測データのlog尤度で表現された尤度確率値に基づく推定のためのゲイン位相ループであって、少なくとも以下を有することを特徴とするゲイン位相ループ。
    ―入力端で、以前の観測データに対し位相とゲインパラメータ、及びゲインとゲイン偏角を補正する項を受け取り、現行観測データに対して推定位相及びゲインのパラメータを出力する加算手段、
    ―フィルタリング モジュールとカスケードに接続する位相及び/又はゲイン角のための機能モジュールであって、
    前記位相及び/又はゲイン機能モジュールは以前の観測データの位相、ゲインパラメータ、前記現行観測データ及びＱＡＭ変調の各シンボルの観測データの全ビットに対するlog尤度のリストを受信し、
    ＱＡＭ変調シンボルの組の現行観測データの重み又は信頼性値により重み付けされた現行観測データの測定位相角の値に比例する値、及び、
    現行観測データと重み又はlog尤度の項で表現される信頼値により重み付けされた、ＱＡＭ変調シンボルの以前の観測データの推定ゲインとの間のゲイン差値に比例する値を、前記フィルタリングモジュールに出力する機能モジュール。
    なお、前記log尤度は前記シンボルの組に対するシンボルに帰し、フィルタは位相角補正項及び/又はゲイン要素を加算手段に出力する。
17. 伝送チャネルを通ってＱＡＭ変調で伝送されるデジタル観測データのための受信機で、
    コンプレックスデモジュレータに対する入力部において、少なくとも組合せとして、現行観測期間の推定の位相ずれ誤差とゲインの値を与えられる位相処理モジュールと、ソフトデマッパーと、観測データのシンボルＱjに帰するlog尤度のリストとしてソフト情報を出力するターボデコーダー、及び、請求項１６のゲインー位相ループを有する受信機。
18. 伝送されるＱＡＭシンボルのインターリーブされた観測データ伝送において、以下を有する請求項１７に記載の受信機。
    ―ターボデコーダの上流側に設けられるデインターリーバモジュール、
    ―受信される観測データの位相の曖昧さを除去することのできる前記ゲインー位相ループの上流側に設けられるインターリーバモジュール。
19. 単一又はマルチキャリア伝送における請求項１乃至１３のいずれか一項に基づく、位相及び/又はゲインを推定するための方法の使用。
20. マルチキャリア伝送において以下の方法による請求項１９に記載の使用。
    α）マルチキャリアシステムを構成する低減された数のサブキャリアに基づいて独立的に使用される方法。
    β）ゲイン及び/又は位相値が、対象サブキャリアの周波数値に応じて、マルチキャリアシステムのサブキャリアに対し補間される。

Images