JP4212453B2 - 尤度演算方法及び分散演算回路 - Google Patents
尤度演算方法及び分散演算回路 Download PDFInfo
- Publication number
- JP4212453B2 JP4212453B2 JP2003377852A JP2003377852A JP4212453B2 JP 4212453 B2 JP4212453 B2 JP 4212453B2 JP 2003377852 A JP2003377852 A JP 2003377852A JP 2003377852 A JP2003377852 A JP 2003377852A JP 4212453 B2 JP4212453 B2 JP 4212453B2
- Authority
- JP
- Japan
- Prior art keywords
- variance
- qam
- qam demodulated
- phase fluctuation
- phase
- 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 - Fee Related
Links
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Description
本発明は軟判定復号に用いる尤度(軟判定値)の演算方法及び尤度を演算するために使用する分散演算回路に係わり、特に、軟判定ＱＡＭ復調信号の分散を該軟判定ＱＡＭ復調信号に作用させて受信シンボルの尤度を出力する尤度演算方法及び分散演算回路に関する。
一般に誤り訂正を行なう場合、復号器には受信シンボルをそのまま入力して復号処理を行なう。近年、誤り訂正能力を向上させる復号方法、例えば軟判定ビタビ（Viterbi）復号法やターボ復号法が良く用いられるようになっている。その中でもターボ復号法のように、受信シンボルの尤度情報を復号器の入力として用いるケースが増えてきている。
受信シンボルの尤度とは、例えば＋１／−１ のバイポーラ信号の受信信号を観測し、受信信号がどのくらい＋１ である可能性が高いか、またはどのくらい＋１ではない可能性が高いか、という「確からしさ」の指標を意味する。このような受信シンボル値から尤度を求める方法として、Sklar の文献 "Turbo codeconcepts made easy, or how I learned to concatenate and reiterate" BernardSklar, 1997, IEEE がある。この文献によれば、受信シンボルr、受信シンボル列の分散σ２ が判れば、この受信シンボルの対数尤度L が、次式で得られるとされている。
このとき、受信シンボル列の分散σ２は、通常受信機入力端付近で付加される熱雑音に起因するものであり、加法性白色ガウス雑音(AWGN: Additive White Gausian Noise)として捕らえることができ、下式
σ２＝kT NF (2)
で得られる。ここでk はボルツマン定数、T は周辺の絶対温度、NF は受信機装置の雑音指数である。
受信シンボルの尤度とは、例えば＋１／−１ のバイポーラ信号の受信信号を観測し、受信信号がどのくらい＋１ である可能性が高いか、またはどのくらい＋１ではない可能性が高いか、という「確からしさ」の指標を意味する。このような受信シンボル値から尤度を求める方法として、Sklar の文献 "Turbo codeconcepts made easy, or how I learned to concatenate and reiterate" BernardSklar, 1997, IEEE がある。この文献によれば、受信シンボルr、受信シンボル列の分散σ２ が判れば、この受信シンボルの対数尤度L が、次式で得られるとされている。
σ２＝kT NF (2)
で得られる。ここでk はボルツマン定数、T は周辺の絶対温度、NF は受信機装置の雑音指数である。
また、分散σ２は受信シンボルr を長期間観測することによっても得ることができる。一般にN 個の信号ri (ｉ=1,2,…Ｎ)の分散は、下式
で得られる。
ここでｒavgはriベクトル振幅の平均値である。この式から明らかなように、分散はN を大きくとるほど真値に近づくものであり、経験的に数百〜数千個以上の受信シンボルから求めるのが望ましい。これら2 つの方法で得られる分散値は、他のゆらぎ要因が無い限り等しいものとなる。
ここでｒavgはriベクトル振幅の平均値である。この式から明らかなように、分散はN を大きくとるほど真値に近づくものであり、経験的に数百〜数千個以上の受信シンボルから求めるのが望ましい。これら2 つの方法で得られる分散値は、他のゆらぎ要因が無い限り等しいものとなる。
図１１は従来の受信装置の要部構成図であり、軟判定復号器の前段までの構成が示されている。ダウンコンバータ1において、受信キャリア発生部１ａは受信キャリアを発生し、位相シフト部１ｂは受信キャリアの位相をπ／２シフトし、乗算器１ｃ，１ｄはべースバンド信号に受信キャリアを乗算してI-CH信号とQ-CH 信号からなるべースバンドの多値QAM 信号を出力し、ローパスフィルタ１ｅ，１ｆは多値QAM 信号の帯域を制限して軟判定QAM 復調器２に入力する。軟判定QAM 復調器２は発振器３から出力するクロック信号に基づいて多値QAM 信号を直列信号にデマッピング(並び替え)して軟判定のQAM 復調信号を発生して尤度演算回路４に入力する。
尤度演算回路４は軟判定ＱＡＭ復調信号の分散σ2を推定し、該分散σ2を軟判定ＱＡＭ復調信号に作用させて受信シンボルの尤度を計算して軟判定復号器に入力する。すなわち、分散推定部４ａは(2)式により受信シンボル列の分散σ2を推定し、演算部４ｂは２／σ2を計算し、乗算器４ｃは(1)式により軟判定QAM復調信号の尤度を計算して軟判定復号器に入力する。
図１１の受信部は(2)式に従って動作するものであるが、実際の受信装置では、熱雑音以外のゆらぎも存在することが明らかである。このため、図１１の構成をとることは殆どなく、(3)式に従って受信シンボル列の分散σ2を推定して尤度を計算して軟判定復号器に入力する図１２の構成が採用されている。図１２において、図１１と異なる点は分散推定部４ａ′が受信シンボル列の分散σ2を(3)式に従って推定する点で、他の構成、動作は同じである。
図１１の受信部は(2)式に従って動作するものであるが、実際の受信装置では、熱雑音以外のゆらぎも存在することが明らかである。このため、図１１の構成をとることは殆どなく、(3)式に従って受信シンボル列の分散σ2を推定して尤度を計算して軟判定復号器に入力する図１２の構成が採用されている。図１２において、図１１と異なる点は分散推定部４ａ′が受信シンボル列の分散σ2を(3)式に従って推定する点で、他の構成、動作は同じである。
次に、QAM 信号と付加雑音の分散について説明する。
一般的に多値QAM を用いるシステムでは、受信機入力端付近で付加される熱雑音が、受信QAM 信号に不確定さを与え、通信システム全体のパフォーマンスを左右する。熱雑音は、各QAM 信号点で大きさの等しいゆらぎとして観測され、更にこのゆらぎの同相(I-CH)成分Viと直交(Q-CH)成分Vqは同量となる。図１３ に16QAM を例として、AWGNを付加した時のコンスタレーションを示す。各信号点において同相(I-CH)成分と直交(Q-CH)成分は等しく、又、各信号点における同相(I-CH)成分Vi a〜Vi dは等しく、さらに各信号点における直交(Q-CH)成分Vq a〜Vq dは等しい。以下においてこのような性質を持った熱雑音に起因する分散をσ２thermal と記す。
4値QAM（QPSK）のようにQAM 信号ベクトル振幅が等しい変調方式においては、QAM 復号シンボルを観測して(3)式を適用することで比較的容易に分散σ２thermalを推定することが可能である。また、16QAM や64QAM のようにQAM 信号ベクトルが複数種類存在するQAM でも、同じQAM信号ベクトルの振幅から生成されたQAM 復調シンボルの集合を観測しそれらの平均値を取ることで、分散σ２thermalを推定することが可能である。熱雑音付加環境下では、このように得られた分散σ２thermalが受信シンボルの分散となっており、(1)式を適用することで、QAM 復調信号の尤度を求めることができる。
一般的に多値QAM を用いるシステムでは、受信機入力端付近で付加される熱雑音が、受信QAM 信号に不確定さを与え、通信システム全体のパフォーマンスを左右する。熱雑音は、各QAM 信号点で大きさの等しいゆらぎとして観測され、更にこのゆらぎの同相(I-CH)成分Viと直交(Q-CH)成分Vqは同量となる。図１３ に16QAM を例として、AWGNを付加した時のコンスタレーションを示す。各信号点において同相(I-CH)成分と直交(Q-CH)成分は等しく、又、各信号点における同相(I-CH)成分Vi a〜Vi dは等しく、さらに各信号点における直交(Q-CH)成分Vq a〜Vq dは等しい。以下においてこのような性質を持った熱雑音に起因する分散をσ２thermal と記す。
4値QAM（QPSK）のようにQAM 信号ベクトル振幅が等しい変調方式においては、QAM 復号シンボルを観測して(3)式を適用することで比較的容易に分散σ２thermalを推定することが可能である。また、16QAM や64QAM のようにQAM 信号ベクトルが複数種類存在するQAM でも、同じQAM信号ベクトルの振幅から生成されたQAM 復調シンボルの集合を観測しそれらの平均値を取ることで、分散σ２thermalを推定することが可能である。熱雑音付加環境下では、このように得られた分散σ２thermalが受信シンボルの分散となっており、(1)式を適用することで、QAM 復調信号の尤度を求めることができる。
さて、シンボルの尤度を求める時、受信信号に付加された雑音の分散σ２ の精度が重要となる。分散σ２ の精度が熱雑音だけに起因する場合には前述の方法で分散σ２(＝σ２thermal)を精度よく推定することができる。しかし、一般に受信QAM 信号で観測される分散σ２ のゆらぎは熱雑音だけが原因でなく、位相ゆらぎも原因の一つになっている。位相ゆらぎは、ダウンコンバータで使用される発振器(キャリア発生器１ａ)の位相雑音、QAM 復調器２で使用される発振器３の位相雑音またはQAM 復調をディジタル信号処理で行なう場合におけるクロックのジッタであり、これらが同時に発生する場合には各位相雑音やジッタを複合したものとなる。
この位相ゆらぎの影響が大きくなると、QAM 信号のコンスタレーションが同相(I)軸と直交(Q)軸の交点Ｏを中心として左右に回転するようなゆらぎが生じる。従って、精度良く分散を推定するには位相ゆらぎを考慮しなければならない。この位相ゆらぎを考慮して尤度を算出する従来技術がある(たとえば特許文献１参照)。この従来技術はQPSK変調信号（4値QAM変調信号）の位相ずれに応じて尤度計算することを特徴とするものである。
特開平５−１４２１３号公報
この位相ゆらぎの影響が大きくなると、QAM 信号のコンスタレーションが同相(I)軸と直交(Q)軸の交点Ｏを中心として左右に回転するようなゆらぎが生じる。従って、精度良く分散を推定するには位相ゆらぎを考慮しなければならない。この位相ゆらぎを考慮して尤度を算出する従来技術がある(たとえば特許文献１参照)。この従来技術はQPSK変調信号（4値QAM変調信号）の位相ずれに応じて尤度計算することを特徴とするものである。
・問題点1
図１４は16QAMを例として、位相雑音を付加した時のコンスタレーション説明図であり、受信QAM 信号点は円弧状のゆらぎを示し、その角度範囲θは、全信号点でほぼ等しくなる。図１５は、16QAM を例として、AWGN と位相雑音を合成付加した時のコンスタレーション説明図であり、図１３ や図１４ と比較して明らかなように、信号点A, B, C, D のゆらぎの形状は同一にならず、I-CHとQ-CH の成分もそれぞれ異なる。図１１及び図１２に示す従来技術の分散推定法では、I-CH とQ-CH でゆらぎが異なることを前提としないため、このような位相ゆらぎが加わると、熱雑音と位相ゆらぎを総合した平均的な分散値を取得することとなり、分散の推定精度が低下してしまう。その結果、尤度演算の精度が低下する。
また、前記文献の従来技術では位相ずれに応じて尤度計算をしているけれど、4値QAM (QPSK)に限定され、4値以外の多値ＱＡＭにおける分散の推定に適用できない。なぜならば、４値ＱＡＭでは信号ベクトルの大きさが同一であるが、4値以外の多値ＱＡＭでは信号ベクトルの大きさが信号点位置に応じて変化し、図１４から明らかなように、位相ゆらぎに起因する分散が信号ベクトルの大きさ(信号点位置)に応じて変化するからである。すなわち、従来技術のQPSK では、信号点が同心円上に配置されているため、信号ベクトルの大きさが一定であり、位相ゆらぎが信号ベクトルの大きさ(信号点位置)に応じて変化する現象は起こらない。このため、従来技術は分散が信号ベクトルの大きさ(信号点位置)に応じて変化する4値以外の多値ＱＡＭにおける分散の推定に適用できない問題がある。
図１４は16QAMを例として、位相雑音を付加した時のコンスタレーション説明図であり、受信QAM 信号点は円弧状のゆらぎを示し、その角度範囲θは、全信号点でほぼ等しくなる。図１５は、16QAM を例として、AWGN と位相雑音を合成付加した時のコンスタレーション説明図であり、図１３ や図１４ と比較して明らかなように、信号点A, B, C, D のゆらぎの形状は同一にならず、I-CHとQ-CH の成分もそれぞれ異なる。図１１及び図１２に示す従来技術の分散推定法では、I-CH とQ-CH でゆらぎが異なることを前提としないため、このような位相ゆらぎが加わると、熱雑音と位相ゆらぎを総合した平均的な分散値を取得することとなり、分散の推定精度が低下してしまう。その結果、尤度演算の精度が低下する。
また、前記文献の従来技術では位相ずれに応じて尤度計算をしているけれど、4値QAM (QPSK)に限定され、4値以外の多値ＱＡＭにおける分散の推定に適用できない。なぜならば、４値ＱＡＭでは信号ベクトルの大きさが同一であるが、4値以外の多値ＱＡＭでは信号ベクトルの大きさが信号点位置に応じて変化し、図１４から明らかなように、位相ゆらぎに起因する分散が信号ベクトルの大きさ(信号点位置)に応じて変化するからである。すなわち、従来技術のQPSK では、信号点が同心円上に配置されているため、信号ベクトルの大きさが一定であり、位相ゆらぎが信号ベクトルの大きさ(信号点位置)に応じて変化する現象は起こらない。このため、従来技術は分散が信号ベクトルの大きさ(信号点位置)に応じて変化する4値以外の多値ＱＡＭにおける分散の推定に適用できない問題がある。
・問題点２
QAM 復調で扱う信号に位相ゆらぎが加わることに対応して、QAM 信号点ごとに分散を正確に求めようとすると、信号点ごとに分散を演算する必要がある。例えば図１２の従来方法の回路に、各QAM 信号点を判別して、それぞれで(3)式の処理を行なう改良が必要となり、回路規模が大きくなり、しかも、分散推定に必要とされるシンボル数が多くなって制御のスピードが低下する問題が生じる。
以上から本発明の目的は4値以外の多値ＱＡＭにおける分散を、位相ゆらぎが加わっても正確に推定して精度の高い尤度を算出できるようにすることである。
本発明の別の目的は構成が簡単で、しかも、シンボル毎に位相ゆらぎを考慮した分散の推定ができるようにすることである。
QAM 復調で扱う信号に位相ゆらぎが加わることに対応して、QAM 信号点ごとに分散を正確に求めようとすると、信号点ごとに分散を演算する必要がある。例えば図１２の従来方法の回路に、各QAM 信号点を判別して、それぞれで(3)式の処理を行なう改良が必要となり、回路規模が大きくなり、しかも、分散推定に必要とされるシンボル数が多くなって制御のスピードが低下する問題が生じる。
以上から本発明の目的は4値以外の多値ＱＡＭにおける分散を、位相ゆらぎが加わっても正確に推定して精度の高い尤度を算出できるようにすることである。
本発明の別の目的は構成が簡単で、しかも、シンボル毎に位相ゆらぎを考慮した分散の推定ができるようにすることである。
第1の発明は多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する尤度演算方法において、受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、該多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出し、該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正し、補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する尤度演算方法である。
第２の発明は、第1の発明の尤度演算方法において、熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する尤度演算方法である。
第２の発明は、第1の発明の尤度演算方法において、熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する尤度演算方法である。
第３の発明は、第１の発明の尤度演算方法において、前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する尤度演算方法である。
第４の発明は、多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、熱雑音に起因する第1の分散を推定する手段、前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段を備えた分散演算回路である。
第5の発明は、多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、前記軟判定QAM復調信号を用いて統計処理により第1の分散を推定する手段、前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、を備えた分散演算回路である。
第４の発明は、多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、熱雑音に起因する第1の分散を推定する手段、前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段を備えた分散演算回路である。
第5の発明は、多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、前記軟判定QAM復調信号を用いて統計処理により第1の分散を推定する手段、前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、を備えた分散演算回路である。
本発明によれば、位相ゆらぎを考慮して正確に分散を演算でき、該分散を用いることにより尤度演算を正確に行なえ、誤り訂正の精度向上に寄与することができる。
また、本発明によれば、多値ＱＡＭの全信号点のシンボル共通に位相ゆらぎ分散算出部を設け、該位相ゆらぎ分散算出部で位相ゆらぎ分散を算出するように構成することで回路規模を小さくすることができ、しかも、シンボル単位で位相ゆらぎを考慮した分散の推定ができる。
また、本発明によれば、多値ＱＡＭの全信号点のシンボル共通に位相ゆらぎ分散算出部を設け、該位相ゆらぎ分散算出部で位相ゆらぎ分散を算出するように構成することで回路規模を小さくすることができ、しかも、シンボル単位で位相ゆらぎを考慮した分散の推定ができる。
受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、該多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散σ２ wobble）を算出し、該位相ゆらぎ分散を用いて別途計算されている多値QAM復調シンボルの分散を補正し、補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する。
熱雑音に起因する分散σ２ thermalを前記別途計算されている多値QAM復調シンボルの分散とし、この分散に位相ゆらぎ分散σ２ wobbleを加算して補正し、該補正後の分散を用いて(1)式により尤度を計算する。あるいは、前記軟判定QAM復調信号を用いて統計処理により算出した分散を、前記別途計算されている多値QAM復調シンボルの分散とし、この分散に位相ゆらぎ分散を加算して補正し、該補正後の分散を用いて尤度を計算する。
位相ゆらぎ分散σ２ wobbleは位相ゆらぎ情報に応じた分散値σ2 θと、位相ゆらぎの角度範囲θと、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出し、
該同相成分と直交成分にデマッピング処理を施して出力する
本発明によれば、位相ゆらぎがコンスタレーションに回転方向のゆらぎを起こす場合であっても、分散を位相ゆらぎを考慮して正確に推定することができ、該分散を用いるて尤度演算を正確に行なえ、誤り訂正の精度向上に寄与することができる。
熱雑音に起因する分散σ２ thermalを前記別途計算されている多値QAM復調シンボルの分散とし、この分散に位相ゆらぎ分散σ２ wobbleを加算して補正し、該補正後の分散を用いて(1)式により尤度を計算する。あるいは、前記軟判定QAM復調信号を用いて統計処理により算出した分散を、前記別途計算されている多値QAM復調シンボルの分散とし、この分散に位相ゆらぎ分散を加算して補正し、該補正後の分散を用いて尤度を計算する。
位相ゆらぎ分散σ２ wobbleは位相ゆらぎ情報に応じた分散値σ2 θと、位相ゆらぎの角度範囲θと、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出し、
該同相成分と直交成分にデマッピング処理を施して出力する
本発明によれば、位相ゆらぎがコンスタレーションに回転方向のゆらぎを起こす場合であっても、分散を位相ゆらぎを考慮して正確に推定することができ、該分散を用いるて尤度演算を正確に行なえ、誤り訂正の精度向上に寄与することができる。
（A）全体の構成
図1は第1実施例の尤度演算回路を備えた受信装置の要部構成図である。
ダウンコンバータ1１において、受信キャリア発生部１１ａは受信キャリアを発生し、位相シフト部１１ｂは受信キャリアの位相をπ／２シフトし、乗算器１１ｃ，１１ｄはべースバンド信号に受信キャリアを乗算してI-CH信号とQ-CH 信号からなるべースバンドの多値QAM 信号（多値ＱＡＭ復調シンボル）を出力し、ローパスフィルタ１１ｅ，１１ｆは多値QAM 信号の帯域を制限して軟判定QAM 復調器１２に入力する。
軟判定QAM 復調器１２は発振器１３から出力するクロック信号CLKに基づいて多値QAM 信号を直列信号に並び替えて（デマッピングして）軟判定のQAM 復調信号にして尤度演算回路１４に入力する。なお、軟判定QAM 復号器１２は、アナログ的な動作でもディジタル的な動作でもどちらで構成してもよい。また、QAM 信号のI-CH とQ-CH の情報(QAM復調シンボル)をどのように並び替えるか（デマッピングするか）については、送受信で整合が取れていれば、それ以上に特別な定義を必要としない。
図1は第1実施例の尤度演算回路を備えた受信装置の要部構成図である。
ダウンコンバータ1１において、受信キャリア発生部１１ａは受信キャリアを発生し、位相シフト部１１ｂは受信キャリアの位相をπ／２シフトし、乗算器１１ｃ，１１ｄはべースバンド信号に受信キャリアを乗算してI-CH信号とQ-CH 信号からなるべースバンドの多値QAM 信号（多値ＱＡＭ復調シンボル）を出力し、ローパスフィルタ１１ｅ，１１ｆは多値QAM 信号の帯域を制限して軟判定QAM 復調器１２に入力する。
軟判定QAM 復調器１２は発振器１３から出力するクロック信号CLKに基づいて多値QAM 信号を直列信号に並び替えて（デマッピングして）軟判定のQAM 復調信号にして尤度演算回路１４に入力する。なお、軟判定QAM 復号器１２は、アナログ的な動作でもディジタル的な動作でもどちらで構成してもよい。また、QAM 信号のI-CH とQ-CH の情報(QAM復調シンボル)をどのように並び替えるか（デマッピングするか）については、送受信で整合が取れていれば、それ以上に特別な定義を必要としない。
位相ゆらぎ分散生成回路１５は、I-CH，Q-CHの多値QAM信号(多値QAM復調シンボル)と受信処理に用いる発振器１１ａ，１３の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）σ2 wobbleを算出して出力する。なお、位相ゆらぎ情報とは、ダウンコンバータで使用される発振器(キャリア発生器１１ａ)の位相雑音、軟判定QAM 復調器に使用される発振器１３の位相雑音またはQAM 復調をディジタル信号処理で行なう場合におけるクロックCLKのジッタあるいは、これらが同時に発生する場合には各位相雑音やジッタを複合した情報である。
尤度演算回路１４において、熱雑音分散推定部１４ａは軟判定ＱＡＭ復調信号の熱雑音に起因する分散σ2 thermalを(2)式に基づいて推定し、分散補正回路１４ｂは次式
σ2 total＝σ2 thermal＋σ2 wobble
により、分散σ2 thermalを位相ゆらぎ分散σ2 wobbleで補正して全体のσ2 totalを出力する。演算部１４ｃは次式
２／σ2 total
を演算し、乗算器14dは(1)式の演算を行なってQAM復調信号の尤度(軟判定値)を計算して次段の軟判定復号器へ出力する。
尤度演算回路１４において、熱雑音分散推定部１４ａは軟判定ＱＡＭ復調信号の熱雑音に起因する分散σ2 thermalを(2)式に基づいて推定し、分散補正回路１４ｂは次式
σ2 total＝σ2 thermal＋σ2 wobble
により、分散σ2 thermalを位相ゆらぎ分散σ2 wobbleで補正して全体のσ2 totalを出力する。演算部１４ｃは次式
２／σ2 total
を演算し、乗算器14dは(1)式の演算を行なってQAM復調信号の尤度(軟判定値)を計算して次段の軟判定復号器へ出力する。
（B）位相ゆらぎ分散生成回路
図2は位相ゆらぎ分散生成回路の構成図である。位相ゆらぎ分散生成回路１５は、極座標変換器１５ａ、位相ゆらぎ→θ変換器１５ｂ、位相ゆらぎ→σ２ θ変換器１５ｃ、位相ゆらぎ分散演算器１５ｄ、デマッパ１５ｅを備え、入力信号として多値QAM信号(多値QAM復調シンボル)のI-CH とQ-CH、および位相ゆらぎ情報が入力されている。
極座標変換器１５ａは、多値QAM復調シンボルのI-CH とQ-CH をそれぞれ実数と虚数としたときの複素数ベクトルから、該ベクトルの大きさR と偏角φを計算して出力する。位相ゆらぎ→θ変換器１５ｂは、位相ゆらぎ情報からゆらぎ量θ（ゆらぎの角度範囲）を得るもので、いくつかの実現方法が存在するがどの方法によるかは特に限定しない。例えば、実測結果を元に作成した後述する変換テーブル(位相ゆらぎ情報とゆらぎの角度範囲θの対応テーブル)を用いたり、パイロット信号を利用した計測方法が考えられる。なお、ゆらぎ量θは、位相ゆらぎのとりうる値の例えば99%を占める角度範囲のような値である。ここで示した99%という数値は例であり、動作を拘束する条件ではない。
図2は位相ゆらぎ分散生成回路の構成図である。位相ゆらぎ分散生成回路１５は、極座標変換器１５ａ、位相ゆらぎ→θ変換器１５ｂ、位相ゆらぎ→σ２ θ変換器１５ｃ、位相ゆらぎ分散演算器１５ｄ、デマッパ１５ｅを備え、入力信号として多値QAM信号(多値QAM復調シンボル)のI-CH とQ-CH、および位相ゆらぎ情報が入力されている。
極座標変換器１５ａは、多値QAM復調シンボルのI-CH とQ-CH をそれぞれ実数と虚数としたときの複素数ベクトルから、該ベクトルの大きさR と偏角φを計算して出力する。位相ゆらぎ→θ変換器１５ｂは、位相ゆらぎ情報からゆらぎ量θ（ゆらぎの角度範囲）を得るもので、いくつかの実現方法が存在するがどの方法によるかは特に限定しない。例えば、実測結果を元に作成した後述する変換テーブル(位相ゆらぎ情報とゆらぎの角度範囲θの対応テーブル)を用いたり、パイロット信号を利用した計測方法が考えられる。なお、ゆらぎ量θは、位相ゆらぎのとりうる値の例えば99%を占める角度範囲のような値である。ここで示した99%という数値は例であり、動作を拘束する条件ではない。
位相ゆらぎ→σ２ θ変換器１５ｃは、位相ゆらぎ情報あるいはゆらぎ量θを入力とし、その分散値σ2 θを得るものである。この変換器もいくつかの実施例を考えられるが、実測結果を元に作成した後述する変換テーブル（位相ゆらぎ情報あるいはゆらぎ量θと分散値σ２ θとの対応テーブル）を使用するのが望ましい。
位相ゆらぎ分散変換器１５ｄは、上述のR、φ、θ、σ２ θ を入力情報として用いてI-CH とQ-CH の位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑを演算して出力する。デマッパ１５ｅは、入力されたI-CH とQ-CH の位相ゆらぎ分散を多値QAM復調シンボル のデマッピングと同様の方法で並び替えて出力する。
位相ゆらぎ分散変換器１５ｄは、上述のR、φ、θ、σ２ θ を入力情報として用いてI-CH とQ-CH の位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑを演算して出力する。デマッパ１５ｅは、入力されたI-CH とQ-CH の位相ゆらぎ分散を多値QAM復調シンボル のデマッピングと同様の方法で並び替えて出力する。
以下、位相ゆらぎ分散生成回路の動作を説明する。
多値QAM復調シンボルのI-CH とQ-CH が位相ゆらぎ分散生成回路１５に入力すると、極座標変換器１５ａはそれぞれを実数と虚数とする複素数ベクトルの大きさRと偏角φを発生する。一方、位相ゆらぎ情報が入力されると、位相ゆらぎ→θ変換器１５ｂは変換テーブルより位相ゆらぎのゆらぎ量θを発生し、位相ゆらぎ→σ２ θ変換器１５ｃは変換テーブルよりゆらぎ量θに基づいて位相ゆらぎの分散σ２ θを発生する。位相ゆらぎ分散演算器１５ｄは、これらのR,、φ、θ、σ２ θを用いて、位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑを演算する。位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑは異なる値を示し、下式
Ｉ-CH：σ2 wobble ｉ＝σ２ θＲθsinφ (4)
Q-CH：σ2 wobble ｑ＝σ２ θＲθcosφ (5)
によって得られる。なお、(4)、(5)式については後述する。このようにして得られたσ2 wobble ｉ，σ2 wobble ｑをデマッパ１５ｅは多値QAM復調シンボル のデマッピングと同様に並び替えて出力し、位相ゆらぎに起因する位相ゆらぎ分散σ2 wobbleを得る。以上により第1実施例では、多値QAM復調シンボルの全ての信号点を別々に観測することなく、共通の位相ゆらぎ分散生成回路１５より位相ゆらぎ分散σ2 wobbleを発生することができる。
多値QAM復調シンボルのI-CH とQ-CH が位相ゆらぎ分散生成回路１５に入力すると、極座標変換器１５ａはそれぞれを実数と虚数とする複素数ベクトルの大きさRと偏角φを発生する。一方、位相ゆらぎ情報が入力されると、位相ゆらぎ→θ変換器１５ｂは変換テーブルより位相ゆらぎのゆらぎ量θを発生し、位相ゆらぎ→σ２ θ変換器１５ｃは変換テーブルよりゆらぎ量θに基づいて位相ゆらぎの分散σ２ θを発生する。位相ゆらぎ分散演算器１５ｄは、これらのR,、φ、θ、σ２ θを用いて、位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑを演算する。位相ゆらぎ分散σ2 wobble ｉ，σ2 wobble ｑは異なる値を示し、下式
Ｉ-CH：σ2 wobble ｉ＝σ２ θＲθsinφ (4)
Q-CH：σ2 wobble ｑ＝σ２ θＲθcosφ (5)
によって得られる。なお、(4)、(5)式については後述する。このようにして得られたσ2 wobble ｉ，σ2 wobble ｑをデマッパ１５ｅは多値QAM復調シンボル のデマッピングと同様に並び替えて出力し、位相ゆらぎに起因する位相ゆらぎ分散σ2 wobbleを得る。以上により第1実施例では、多値QAM復調シンボルの全ての信号点を別々に観測することなく、共通の位相ゆらぎ分散生成回路１５より位相ゆらぎ分散σ2 wobbleを発生することができる。
第1実施例によれば、位相ゆらぎを考慮して正確に分散を演算でき、該分散を用いることにより尤度演算を正確に行なえ、誤り訂正の精度向上に寄与することができる。
又、第1実施例によれば、長期間のシンボル観測を必要とせずに分散の推定を行なうことが可能である。
また、第1実施例によれば、多値ＱＡＭの全信号点のシンボル共通に位相ゆらぎ分散算出部を設け、該位相ゆらぎ分散算出部で位相ゆらぎ分散を算出するように構成したからは回路規模を小さくすることができる。
以下において、ゆらぎ量θとσ２ θの求め方、位相ゆらぎ分散σ2 wobbleの算出法、デマッパの動作について補足する。
又、第1実施例によれば、長期間のシンボル観測を必要とせずに分散の推定を行なうことが可能である。
また、第1実施例によれば、多値ＱＡＭの全信号点のシンボル共通に位相ゆらぎ分散算出部を設け、該位相ゆらぎ分散算出部で位相ゆらぎ分散を算出するように構成したからは回路規模を小さくすることができる。
以下において、ゆらぎ量θとσ２ θの求め方、位相ゆらぎ分散σ2 wobbleの算出法、デマッパの動作について補足する。
（Ｃ）ゆらぎ量θとσ２ θ
ゆらぎ量θとσ２ θの算出方法ついて説明する。
ゆらぎ量θは、位相ゆらぎ（位相雑音）の大きさによって変化するものであり、例えば発振器の位相雑音が大きくなるほど、θも大きくなる。また、ゆらぎ量θの確率密度分布は正規分布で近似することが可能である。角度方向のゆらぎの確率密度p(θ)は、下式
のようになり、θ＜1 で正規分布に近似可能である（John G.Proakis 著Digital Communications を参照）。ここで、Esはシンボルエネルギ、σr２はI 軸とQ 軸方向に加える雑音密度で
である。ただし、N0は単位周波数当りの雑音電力である。
ゆらぎ量θとσ２ θの算出方法ついて説明する。
ゆらぎ量θは、位相ゆらぎ（位相雑音）の大きさによって変化するものであり、例えば発振器の位相雑音が大きくなるほど、θも大きくなる。また、ゆらぎ量θの確率密度分布は正規分布で近似することが可能である。角度方向のゆらぎの確率密度p(θ)は、下式
図３に正規分布p(x)とゆらぎ量θの確率密度分布p(θ)との比較、および確率密度分布p(θ)が位相ゆらぎ情報(位相雑音PN)でどのように変化するかを示す。なお、位相ゆらぎ情報として位相雑音PNのみを想定している。位相雑音PNを変化させて確率密度分布p(θ)が正規分布p(x)と同一形状になる位相雑音PNと分散σを求めて対応づけることにより位相ゆらぎ(位相雑音)→θ、位相ゆらぎ→σ２ θの変換テーブルを得ることができる。正規分布p(x)と確率密度分布p(θ)とが同一形状になる分散σと位相雑音PNを対応づけると、図3（A）,(B)における正規分布p(x) Ａ〜Ｄと確率密度分布p(θ) Ａ′〜Ｄ′とが同一形状となり、図4に示す対応テーブルが得られる。この対応テーブルを用いて、位相ゆらぎ→θ変換器１５ｂ、位相ゆらぎ→σ２ θ変換器１５ｃは、位相ゆらぎ情報(位相雑音)に応じたゆらぎ量θおよびσ２ θを求めることができる。
（Ｄ）位相ゆらぎ分散σ2 wobbleの演算
受信多値QAM 信号のゆらぎは、熱雑音と位相ゆらぎの合成で考えることができ、その合成結果を示す分散σ2 totalは、熱雑音に起因するゆらぎの分散σ2 thermalと位相ゆらぎ分散σ2 wobbleの和で得られる。
熱雑音は加法性白色ガウス雑音(AWGN)として捕らえることができる。このAWGN が多値QAM に与えるゆらぎを考察する。図５に16QAM を例にQAM の第１象現で観測される4つの信号点Ａ〜ＤについてAWGN 付加時のゆらぎを模式的に示す。なお、信号点Ａは＋１＝“１”，−１＝“０”とすれば、”１０１１”の信号点、信号点Bは”１1１１”の信号点、信号点Ｃは”１０１０”の信号点、信号点Ｄは”１１１０”の信号点である。
QAM の性質上、AWGN はI-CH 成分とQ-CH 成分でほぼ同量で加わり、それはQAM 信号点によらず一定となる。そのため、ゆらぎの振幅範囲Ｖ(σ2 thermal)は送信信号点を中心にほぼ円形となり、I-CH上でのゆらぎ振幅範囲Ｖ(σ2 thermal ｉ)とQ-CH 上でのゆらぎ振幅範囲Ｖ(σ2 thermalｑ)は等しくなる。ここで、Ｖ(・)は、振幅方向のゆらぎ量が.・の関数であることを示す。この円は、正規分布するI-CH 信号とQ-CH 信号の振幅のうち、例えば99%を包括する範囲を示すものであり、AWGN の分散に比例して円の面積は広がる。このように、QAM 復調信号の各シンボルに見られる分散は一定となり、これらを観測することでQAM 復調信号の分散を正確に求めることができる。
受信多値QAM 信号のゆらぎは、熱雑音と位相ゆらぎの合成で考えることができ、その合成結果を示す分散σ2 totalは、熱雑音に起因するゆらぎの分散σ2 thermalと位相ゆらぎ分散σ2 wobbleの和で得られる。
熱雑音は加法性白色ガウス雑音(AWGN)として捕らえることができる。このAWGN が多値QAM に与えるゆらぎを考察する。図５に16QAM を例にQAM の第１象現で観測される4つの信号点Ａ〜ＤについてAWGN 付加時のゆらぎを模式的に示す。なお、信号点Ａは＋１＝“１”，−１＝“０”とすれば、”１０１１”の信号点、信号点Bは”１1１１”の信号点、信号点Ｃは”１０１０”の信号点、信号点Ｄは”１１１０”の信号点である。
QAM の性質上、AWGN はI-CH 成分とQ-CH 成分でほぼ同量で加わり、それはQAM 信号点によらず一定となる。そのため、ゆらぎの振幅範囲Ｖ(σ2 thermal)は送信信号点を中心にほぼ円形となり、I-CH上でのゆらぎ振幅範囲Ｖ(σ2 thermal ｉ)とQ-CH 上でのゆらぎ振幅範囲Ｖ(σ2 thermalｑ)は等しくなる。ここで、Ｖ(・)は、振幅方向のゆらぎ量が.・の関数であることを示す。この円は、正規分布するI-CH 信号とQ-CH 信号の振幅のうち、例えば99%を包括する範囲を示すものであり、AWGN の分散に比例して円の面積は広がる。このように、QAM 復調信号の各シンボルに見られる分散は一定となり、これらを観測することでQAM 復調信号の分散を正確に求めることができる。
次に、位相ゆらぎがQAM に与えるゆらぎを考察する。図６に16QAMを例にQAM の第１象現で観測される位相ゆらぎを模式的に示す。位相ゆらぎは、図１４ に示したのと同様にI 軸とQ 軸の交点を中心として左右に回転するようなゆらぎであり、ゆらぎの角度範囲θは全信号点で同じ量に観測される。上述のように、ゆらぎの角度範囲θは位相のゆらぎ量に対して正規分布で近似され、分散σ２ θを得ることができる。そのゆらぎの角度範囲を例えば99%を包括する範囲を選ぶと図６のような円弧形状となり、これをI 軸とQ 軸に投影したものが、I-CH にあらわれるゆらぎとQ-CH にあらわれるのゆらぎとなり、それぞれが正規分布に近似できる。円弧長Larc ＝Ｒθであるから、I-CH の分散とQ-CH の分散σ2 wobble ｉ，σ2 wobble ｑは、それぞれ(4)式と(5)式で近似できる。
熱雑音と位相ゆらぎが同時に起こる場合、トータルの分散はそれぞれの分散を複合したものとなる。図７に16QAM を例にQAM の第１象現で観測されるゆらぎを模式的に示す。このときのゆらぎの振幅範囲は図５と図６で確認したゆらぎ範囲の複合として扱うことができる。その結果、I-CH とQ-CH のゆらぎの振幅範囲は下式
I-CH のゆらぎ振幅範囲：
Ｖ(σ2 total ｉ)＝Ｖ(σ2 thermalｉ)＋Ｖ(σ2 wobble ｉ) (7)
Q-CH のゆらぎ振幅範囲：
Ｖ(σ2 total ｑ)＝Ｖ(σ2 thermalｑ)＋Ｖ(σ2 wobble ｑ) (8)
となる。つまり、QAM 復調シンボルの分散は、QAM 信号点およびI-CH とQ-CH のそれぞれで検討しなければならない。このようにして得た分散σ2 totalを用いて尤度演算を行なう。
I-CH のゆらぎ振幅範囲：
Ｖ(σ2 total ｉ)＝Ｖ(σ2 thermalｉ)＋Ｖ(σ2 wobble ｉ) (7)
Q-CH のゆらぎ振幅範囲：
Ｖ(σ2 total ｑ)＝Ｖ(σ2 thermalｑ)＋Ｖ(σ2 wobble ｑ) (8)
となる。つまり、QAM 復調シンボルの分散は、QAM 信号点およびI-CH とQ-CH のそれぞれで検討しなければならない。このようにして得た分散σ2 totalを用いて尤度演算を行なう。
(Ｅ) デマッパの動作
図2のデマッパ１５ｅは、マッピングのルールに従ってI-CH 位相ゆらぎ分散σ2 wobble ｉとQ-CH位相ゆらぎ分散σ2 wobble ｑを並び替えて出力する。16QAM のマッピングが、前半2 ビットをI-CH、後半2 ビットをQ-CHと定義していると仮定する。このときの信号点ＡにおけるQAM 信号、軟判定QAM 復調信号、および硬判定QAM 復調信号は以下の例のようになる。
デマッパ１５ｅは、マッピングのルールに従ってI-CH 位相ゆらぎ分散σ2 wobble ｉとQ-CH位相ゆらぎ分散σ2 wobble ｑを並び替えて出力する(デマッピング)。上記の例では、デマッピングにより以下に示す位相ゆらぎ分散
［σ2 wobble ｉσ2 wobble ｉσ2 wobble ｑσ2 wobble ｑ］
を出力する。
図2のデマッパ１５ｅは、マッピングのルールに従ってI-CH 位相ゆらぎ分散σ2 wobble ｉとQ-CH位相ゆらぎ分散σ2 wobble ｑを並び替えて出力する。16QAM のマッピングが、前半2 ビットをI-CH、後半2 ビットをQ-CHと定義していると仮定する。このときの信号点ＡにおけるQAM 信号、軟判定QAM 復調信号、および硬判定QAM 復調信号は以下の例のようになる。
［σ2 wobble ｉσ2 wobble ｉσ2 wobble ｑσ2 wobble ｑ］
を出力する。
図8は第2実施例の尤度演算回路を備えた受信装置の要部構成図であり、図１の第1実施例と同一部分には同一符号を付している。第1実施例と異なる点は、(1)尤度演算回路１４において熱雑音分散推定部１４ａの替わりに分散推定部５１を設け、軟判定ＱＡＭ復調器１２から出力する軟判定QAM復調信号を用いて統計処理により熱雑音に起因する分散を演算する点、および(2)位相ゆらぎ分散生成回路１５による位相ゆらぎ分散生成方法である。
すなわち、第1実施例と第2実施例の動作の違いは、熱雑音の分散σ2 thermalの求め方にある。第1実施例 では、軟判定QAM 復調信号を用いずに、純粋な熱雑音の分散σ2 thermalを(2)式に基づいて算出して尤度演算を実現している。一方、第2実施例では、位相ゆらぎを含んだ軟判定QAM復調信号を用いて(3)式に基づいて統計的に熱雑音の分散を求めている。このため、第2実施例では位相ゆらぎの影響を含んだ熱雑音の分散
を用いることになる。
すなわち、第1実施例と第2実施例の動作の違いは、熱雑音の分散σ2 thermalの求め方にある。第1実施例 では、軟判定QAM 復調信号を用いずに、純粋な熱雑音の分散σ2 thermalを(2)式に基づいて算出して尤度演算を実現している。一方、第2実施例では、位相ゆらぎを含んだ軟判定QAM復調信号を用いて(3)式に基づいて統計的に熱雑音の分散を求めている。このため、第2実施例では位相ゆらぎの影響を含んだ熱雑音の分散
図５〜図7を参照して第１実施例、第2実施例の熱雑音
の関係を考察する。図５ と図７を比較して分かるように、「熱雑音のみの分散」より、「熱雑音に位相ゆらぎを加えた分散」が大きくなる。この関係は信号点やI-CH、Q-CH によらないものであり、以下
σ2 total ｉ＞σ2 thermal ｉ (10)
σ2 total ｑ＞σ2 thermal ｑ) (11)
の関係で示される。また、σ2 total ｉやσ2 total ｑを詳しくみると、信号点によって異なる値を示すことがわかる。例えば、信号点A と信号点D を比べたとき、以下
σ2 total_ｉ_Ａ＞σ2 total_ｉ_D (12)
σ2 total_ｑ_Ａ＜σ2 total_ｑ_D (13)
の関係が成り立つ。ここで、σ2 total_ｉやσ2 total_ｑの添え字_A，D は、信号点の区別を示す。
σ2 total ｉ＞σ2 thermal ｉ (10)
σ2 total ｑ＞σ2 thermal ｑ) (11)
の関係で示される。また、σ2 total ｉやσ2 total ｑを詳しくみると、信号点によって異なる値を示すことがわかる。例えば、信号点A と信号点D を比べたとき、以下
σ2 total_ｉ_Ａ＞σ2 total_ｉ_D (12)
σ2 total_ｑ_Ａ＜σ2 total_ｑ_D (13)
の関係が成り立つ。ここで、σ2 total_ｉやσ2 total_ｑの添え字_A，D は、信号点の区別を示す。
以上のように位相ゆらぎが加わったときの軟判定QAM 復調信号の分散は、シンボル毎に異なる値を示し、従って、第2実施例のように軟判定QAM 復調信号から求められる(9)式の分散は、QAM の第一象現にある各信号点の分散の平均値となり、図７の16QAM では下式
のようになる。
図9は第2実施例における位相ゆらぎ分散生成回路１５の構成図であり、図2の第１実施例における位相ゆらぎ分散生成回路１５と同一部分には同一符号を付している。異なる点は、分散補正係数KMを生成する分散補正係数生成器１５ｆ、−KMσ２ θθを演算する演算器１５ｇ、デマッパ１５ｅから出力する分散σ2 wobbleに−KMσ２ θθを加算して位相ゆらぎ分散
を出力する加算器１５ｈを設けた点である。
一般に位相ゆらぎの影響が無い場合、純粋に熱雑音のみで観測される分散σ2 thermalと、QAM復調信号を観測して得られる(9)式の分散は、ほぼ等しくなる。しかし、位相ゆらぎの影響を受けたQAM復調信号を観測して得られる(9)式の分散は下式
のようになる。ここで、
は分散補正係数であり、QAMの多値数に依存する係数となり、図１０に示すようになる。この数値の根拠については後述する。また、σ２ θやθは、図４の変換テーブルより得られる。
は分散補正係数であり、QAMの多値数に依存する係数となり、図１０に示すようになる。この数値の根拠については後述する。また、σ２ θやθは、図４の変換テーブルより得られる。
第2実施例では、下式を実現させようとしていることに特徴がある。
この式に、上記の(16)式を代入すると、下式
が得られる。第2実施例の位相ゆらぎ分散生成回路１５は上式の後半部分
を演算して出力する。
以上第2実施例によれば、第1実施例の効果に加えて、軟判定QAM復調信号を用いて統計的に熱雑音に起因する分散を算出する場合であっても、位相ゆらぎを考慮して分散を補正演算することが出来、結果的に正しく尤度を算出することができる。
以上第2実施例によれば、第1実施例の効果に加えて、軟判定QAM復調信号を用いて統計的に熱雑音に起因する分散を算出する場合であっても、位相ゆらぎを考慮して分散を補正演算することが出来、結果的に正しく尤度を算出することができる。
以下において、図１０の分散補正係数の求め方について補足する説明する。
位相ゆらぎの影響を受けている時の分散推定回路５１の動作は、「直交座標の信号点や実数/虚数成分毎に得られる分散」の平均値を求めるようになる。
一般的な多値QAMは、次の条件が成り立つ。
(1)QAMシンボルの生起確率は、直交座標の信号点のすべてで等確率である。
(2)直交座標の第１象現、第２象現、第３象現、第４象現で観測される分散は、対称なものとなる。すなわち、全体から求めた分散の平均値と、第１象現のみから求めた分散の平均値は等しくなる。
(3) I軸に見られる分散とQ軸に見られる分散は、対称なものとして観測される。すなわち、全体から求めた分散の平均値と、I軸(実数)成分のみから求めた分散の平均値は等しくなる。
これらの条件は、QAMでは一般的な概念である。これにそぐわないQAMは非常に特殊な方式となる。以上の条件によって、第１象現の受信シンボルのI軸上のばらつきを評価することにより、QAM全体に対する分散の平均値を求めることが可能となる。
位相ゆらぎの影響を受けている時の分散推定回路５１の動作は、「直交座標の信号点や実数/虚数成分毎に得られる分散」の平均値を求めるようになる。
一般的な多値QAMは、次の条件が成り立つ。
(1)QAMシンボルの生起確率は、直交座標の信号点のすべてで等確率である。
(2)直交座標の第１象現、第２象現、第３象現、第４象現で観測される分散は、対称なものとなる。すなわち、全体から求めた分散の平均値と、第１象現のみから求めた分散の平均値は等しくなる。
(3) I軸に見られる分散とQ軸に見られる分散は、対称なものとして観測される。すなわち、全体から求めた分散の平均値と、I軸(実数)成分のみから求めた分散の平均値は等しくなる。
これらの条件は、QAMでは一般的な概念である。これにそぐわないQAMは非常に特殊な方式となる。以上の条件によって、第１象現の受信シンボルのI軸上のばらつきを評価することにより、QAM全体に対する分散の平均値を求めることが可能となる。
16QAMを例に、上記の条件を考慮して、図８の分散推定部５１の出力
がどのようになるかを示したものが(14)式である。すなわち、I軸上に投影される第１象現の4点の分散を求め、その平均値を取ると、
となる。ここで、
は、信号点Xにおける熱雑音と位相ゆらぎを分離せず、またI軸に投影された分散を意味する。1/4は、16QAMの第１象現に配置される信号点の数で平均を取ることを意味している。
(14)式の右辺各項について、位相ゆらぎに依存しない成分（即ち、熱雑音に依存する成分）と、依存する成分に分離すると次式
のようになる。この結果を用いると、(14)式は以下
になる。さらに
は、(4)式で置き換えることができ、下式
が得られる。
(22)式を(21)に代入することで、下式
が得られる。上式右辺の（ ）内について検討する。
一般な16QAMでは、シンボルエネルギの平均値を1とするように、直交座標の値（電圧または電流）を決める。ただし、最小振幅または最大振幅のシンボルエネルギを1とする考えもあるが、その場合RA〜RDが違う値となるので、KMが変わってくる。
シンボルエネルギの平均値を1とするような16QAMの場合、RA〜RD，sinφA~sinφDは、以下のように求められる。
一般な16QAMでは、シンボルエネルギの平均値を1とするように、直交座標の値（電圧または電流）を決める。ただし、最小振幅または最大振幅のシンボルエネルギを1とする考えもあるが、その場合RA〜RDが違う値となるので、KMが変わってくる。
シンボルエネルギの平均値を1とするような16QAMの場合、RA〜RD，sinφA~sinφDは、以下のように求められる。
同様に64QAMでは下式
となる。更に256QAMでは
となる。
・付記
付記１．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する尤度演算方法において、
受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、
該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、
該多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出し、
該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正し、
補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする尤度演算方法。
付記２．前記位相ゆらぎ情報に応じた分散値と、位相ゆらぎの角度範囲と、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出し、
該同相成分と直交成分にデマッピング処理を施して前記位相ゆらぎ分散を出力する、
ことを特徴とする付記１記載の尤度演算方法。
付記３．熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記１記載の尤度演算方法。
付記４．前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記１記載の尤度演算方法。
付記５．予め前記位相ゆらぎ情報に対応させて前記分散値と位相ゆらぎの角度範囲とを記憶しておく、
ことを特徴とする付記２記載の尤度演算方法。
付記６．受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、前記多値QAM復調シンボルの分散を用いて前記軟判定QAM復調信号の尤度情報を演算して出力する尤度演算装置において、
前記多値QAM復調シンボルと、受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出する位相ゆらぎ分散算出部、
該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正する分散補正部、
補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する尤度演算部、
を備えたことを特徴とする尤度演算装置。
付記７．前記位相ゆらぎ分散算出部は、前記位相ゆらぎ情報に応じた分散値と位相ゆらぎの角度範囲と、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出する手段、
該同相成分と直交成分にデマッピング処理を施して前記位相ゆらぎ分散を出力する手段、
を備えたことを特徴とする付記６記載の尤度演算装置。
付記８．前記分散補正部は、熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、前記尤度演算部は、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記６または７記載の尤度演算装置。
付記９．前記分散補正部は、前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、前記尤度演算部は、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記６または７記載の尤度演算装置。
付記１０．予め前記位相ゆらぎ情報に対応させて前記分散値と位相ゆらぎの角度範囲とを記憶する記憶手段を備え、前記位相ゆらぎ分散算出部は該記憶手段より前記位相ゆらぎ情報に応じた分散値と位相ゆらぎの角度範囲を取得する、
ことを特徴とする付記７記載の尤度演算装置。
付記１１．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
熱雑音に起因する第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。
付記１２．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
前記軟判定QAM復調信号を用いて統計処理により第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。
付記１．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する尤度演算方法において、
受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、
該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、
該多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出し、
該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正し、
補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする尤度演算方法。
付記２．前記位相ゆらぎ情報に応じた分散値と、位相ゆらぎの角度範囲と、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出し、
該同相成分と直交成分にデマッピング処理を施して前記位相ゆらぎ分散を出力する、
ことを特徴とする付記１記載の尤度演算方法。
付記３．熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記１記載の尤度演算方法。
付記４．前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記１記載の尤度演算方法。
付記５．予め前記位相ゆらぎ情報に対応させて前記分散値と位相ゆらぎの角度範囲とを記憶しておく、
ことを特徴とする付記２記載の尤度演算方法。
付記６．受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、前記多値QAM復調シンボルの分散を用いて前記軟判定QAM復調信号の尤度情報を演算して出力する尤度演算装置において、
前記多値QAM復調シンボルと、受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出する位相ゆらぎ分散算出部、
該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正する分散補正部、
補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する尤度演算部、
を備えたことを特徴とする尤度演算装置。
付記７．前記位相ゆらぎ分散算出部は、前記位相ゆらぎ情報に応じた分散値と位相ゆらぎの角度範囲と、前記多値QAM復調シンボルとを用いてシンボル点における位相ゆらぎ分散の同相成分と直交成分を算出する手段、
該同相成分と直交成分にデマッピング処理を施して前記位相ゆらぎ分散を出力する手段、
を備えたことを特徴とする付記６記載の尤度演算装置。
付記８．前記分散補正部は、熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、前記尤度演算部は、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記６または７記載の尤度演算装置。
付記９．前記分散補正部は、前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、前記尤度演算部は、該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする付記６または７記載の尤度演算装置。
付記１０．予め前記位相ゆらぎ情報に対応させて前記分散値と位相ゆらぎの角度範囲とを記憶する記憶手段を備え、前記位相ゆらぎ分散算出部は該記憶手段より前記位相ゆらぎ情報に応じた分散値と位相ゆらぎの角度範囲を取得する、
ことを特徴とする付記７記載の尤度演算装置。
付記１１．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
熱雑音に起因する第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。
付記１２．多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
前記軟判定QAM復調信号を用いて統計処理により第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。
１１ ダウンコンバータ
１１ａ 受信キャリア発生部
１２ 軟判定QAM 復調器
１３ 発振器
１４ 尤度演算回路
１４ａ 熱雑音分散推定部
１４ｂ 分散補正回路
１５ 位相ゆらぎ分散生成回路
１１ａ 受信キャリア発生部
１２ 軟判定QAM 復調器
１３ 発振器
１４ 尤度演算回路
１４ａ 熱雑音分散推定部
１４ｂ 分散補正回路
１５ 位相ゆらぎ分散生成回路
Claims (5)
- 多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する尤度演算方法において、
受信信号に直交復調処理を施して多値QAM復調シンボルを発生し、
該多値QAM復調シンボルにデマッピング処理を施して軟判定QAM復調信号を発生し、
該多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する分散（位相ゆらぎ分散）を算出し、
該位相ゆらぎ分散を用いて前記多値QAM復調シンボルの分散を補正し、
補正後の分散を用いて前記軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする尤度演算方法。 - 熱雑音に起因する分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする請求項１記載の尤度演算方法。 - 前記軟判定QAM復調信号を用いて統計処理により算出した分散を前記多値QAM復調シンボルの分散として求め、該分散に前記位相ゆらぎ分散を加算して前記多値QAM復調シンボルの分散を補正し、
該補正後の分散を前記軟判定QAM復調信号に作用させて軟判定QAM復調信号の尤度情報を演算する、
ことを特徴とする請求項１記載の尤度演算方法。 - 多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
熱雑音に起因する第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。 - 多値QAM復調シンボルの分散を用いて軟判定QAM復調信号の尤度情報を演算する受信装置の分散演算回路において、
前記軟判定QAM復調信号を用いて統計処理により第1の分散を推定する手段、
前記多値QAM復調シンボルと受信処理に用いる発振器の位相雑音に応じた位相ゆらぎ情報とを用いて、位相ゆらぎに起因する第2の分散を算出する手段、
前記第１の分散と第２の分散を加算して前記多値QAM復調シンボルの分散を算出する手段、
を備えたことを特徴とする分散演算回路。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003377852A JP4212453B2 (ja) | 2003-11-07 | 2003-11-07 | 尤度演算方法及び分散演算回路 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003377852A JP4212453B2 (ja) | 2003-11-07 | 2003-11-07 | 尤度演算方法及び分散演算回路 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005142880A JP2005142880A (ja) | 2005-06-02 |
| JP4212453B2 true JP4212453B2 (ja) | 2009-01-21 |
Family
ID=34688425
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003377852A Expired - Fee Related JP4212453B2 (ja) | 2003-11-07 | 2003-11-07 | 尤度演算方法及び分散演算回路 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4212453B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11162574B2 (en) | 2017-10-16 | 2021-11-02 | Mitsuba Corporation | Speed reduction mechanism and motor with speed reducer |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5034588B2 (ja) * | 2007-03-20 | 2012-09-26 | 住友電気工業株式会社 | 軟判定復号装置および軟判定復号プログラム |
-
2003
- 2003-11-07 JP JP2003377852A patent/JP4212453B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11162574B2 (en) | 2017-10-16 | 2021-11-02 | Mitsuba Corporation | Speed reduction mechanism and motor with speed reducer |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005142880A (ja) | 2005-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5955481B2 (ja) | 軟判定値生成装置及び軟判定値生成方法 | |
| US7480342B2 (en) | Soft value calculation for multilevel signals | |
| Kenarsari-Anhari et al. | Performance analysis for BICM transmission over Gaussian mixture noise fading channels | |
| AU2015207321B2 (en) | Communication apparatus, demodulation apparatus, carrier reproduction apparatus, phase error compensation apparatus, phase error compensation method, and storage medium on which phase error compensation program has been stored | |
| Simon et al. | On the optimality of bit detection of certain digital modulations | |
| US6452983B1 (en) | Demodulator for improving reception performance in satellite and Mobile communications | |
| US8705665B2 (en) | Process for performing log-likelihood-ratio clipping in a soft-decision near-ML detector, and detector for doing the same | |
| US7680219B2 (en) | Apparatus and method for decoding a bit sequence from QPSK or QAM symbols | |
| JP2007166592A (ja) | 反復デコーダーのための高次変調方式に適する軟判定デマップ方法及びそれを利用したエラー訂正装置 | |
| JP4212453B2 (ja) | 尤度演算方法及び分散演算回路 | |
| US9531577B2 (en) | Bit-likelihood calculating apparatus and bit-likelihood calculating method | |
| JP2008502246A (ja) | ターボ復号化器の位相シフトの補償システム | |
| Park | Bitwise log-likelihood ratios for quadrature amplitude modulations | |
| EP1343285B1 (en) | Soft value calculation for multilevel signals | |
| US7684522B2 (en) | Method and system for determining a log-likelihood ratio (LLR) corresponding to each bit of a symbol | |
| US20020122507A1 (en) | System for carrier phase tracking of multi-dimensional coded symbols | |
| US6795512B1 (en) | System for carrier phase tracking of coded symbols using reliability metrics for symbol estimates | |
| JP3817614B2 (ja) | 位相雑音耐性の高い受信装置用の最適検波器 | |
| JPH1127156A (ja) | 誤り訂正復号装置 | |
| WO2002098090A1 (en) | Method and apparatus for soft information generation | |
| JP3237866B2 (ja) | 軟判定復号方法 | |
| TW201815139A (zh) | 符元判斷方法、符元判斷電路以及數位接收電路 | |
| KR100903876B1 (ko) | 비트 대칭 그레이 코드를 이용하여 위상 편이 방식으로변조된 수신 심볼 신호를 비트 정보로 분할하는 방법 및 그장치 | |
| JP2002330187A (ja) | 復号装置および復号方法 | |
| JP2003273751A (ja) | ビット尤度算出方法及び復調装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060810 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20081022 |
|
| 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: 20081028 |
|
| 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: 20081028 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111107 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: 20111107 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121107 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121107 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131107 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |