JPWO2011158932A1 - 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法 - Google Patents

周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法 Download PDF

Info

Publication number
JPWO2011158932A1
JPWO2011158932A1 JP2012520505A JP2012520505A JPWO2011158932A1 JP WO2011158932 A1 JPWO2011158932 A1 JP WO2011158932A1 JP 2012520505 A JP2012520505 A JP 2012520505A JP 2012520505 A JP2012520505 A JP 2012520505A JP WO2011158932 A1 JPWO2011158932 A1 JP WO2011158932A1
Authority
JP
Japan
Prior art keywords
frequency
frequency offset
received signal
offset
estimation
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.)
Granted
Application number
JP2012520505A
Other languages
English (en)
Other versions
JP5404926B2 (ja
Inventor
中川　匡夫
匡夫 中川
松井　宗大
宗大 松井
浩一 石原
浩一 石原
理一 工藤
理一 工藤
小林　孝行
孝行 小林
山崎　悦史
悦史 山崎
佐野　明秀
明秀 佐野
英二 吉田
英二 吉田
匡人 溝口
匡人 溝口
宮本　裕
宮本　　裕
泰司 鷹取
泰司 鷹取
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2012520505A priority Critical patent/JP5404926B2/ja
Publication of JPWO2011158932A1 publication Critical patent/JPWO2011158932A1/ja
Application granted granted Critical
Publication of JP5404926B2 publication Critical patent/JP5404926B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/10—Polarisation diversity; Directional diversity
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H04B10/616—Details of the electronic signal processing in coherent optical receivers
    • H04B10/6164—Estimation or correction of the frequency offset between the received optical signal and the optical local oscillator
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H04B10/65—Intradyne, i.e. coherent receivers with a free running local oscillator having a frequency close but not phase-locked to the carrier signal
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H04B17/318—Received signal strength
    • H04B17/327—Received signal code power [RSCP]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
    • H04L27/14—Demodulator circuits; Receiver circuits
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • H04L27/2659—Coarse or integer frequency offset determination and synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset
    • H04L2027/0032—Correction of carrier offset at baseband and passband
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0044—Control loops for carrier regulation
    • H04L2027/0046—Open loops

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Circuits Of Receivers In General (AREA)
  • Optical Communication System (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Superheterodyne Receivers (AREA)

Abstract

受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する場合に、受信信号の周波数オフセットを適切に推定し補償する、周波数オフセット推定装置を提供する。周波数オフセット推定装置は、所定のサンプリング周波数で予めサンプリングされた受信信号を、Ｎ個の周波数成分を持つ周波数スペクトルに変換し、周波数スペクトルの、１からＮ/２までの負の周波数成分とＮ/２＋１からＮまでの正の周波数成分とを、それぞれ周波数帯域制限し、周波数帯域制限された周波数スペクトルの正の周波数成分と負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、それぞれの電力の電力差の絶対値が予め設定されたしきい値以下になるまで周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、しきい値以下になるまで移動させた移動量に基づいて周波数オフセットを推定する。

Description

本発明は、ディジタルコヒーレント光受信機および無線通信受信機における周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法に関する。
本願は、２０１０年６月１７日に日本へ出願された特願２０１０−１３８４０２号、および、２０１０年１１月１０日に日本へ出願された特願２０１０−２５１８６８号に基づき優先権を主張し、それらの内容をここに援用する。
光通信の分野において、周波数利用効率を飛躍的に向上させる同期検波方式とディジタル信号処理を組み合わせたディジタルコヒーレント通信システムが注目されている。ディジタルコヒーレント通信システムは、直接検波により構築されていたシステムと比較すると、受信感度を向上させることができるだけでなく、送信信号をディジタル信号として受信することで、光ファイバ伝送によって受ける波長分散や偏波モード分散による送信信号の波形歪みを補償することができることが知られており、次世代の光通信技術として導入が検討されている。
コヒーレント受信機において受信された信号光は、局部発振光と掛け合わされてベースバンド信号に変換される。信号光の搬送波や局部発振光を生成するレーザ発振器は、無線通信用の発振器で一般に用いられている位相同期ループによる周波数安定化が困難であり、送信機のレーザ発振器の出力周波数と受信機のレーザ発振器の出力周波数との間に大きな周波数オフセットが生じる。実際の光通信システムにおいて周波数オフセットは、±5GHzに達する。コヒーレント通信システムでは搬送波の位相に情報を乗せているため、受信機において周波数オフセットを推定し、補償する必要がある。
また無線通信においては、送信機と受信機に用いている基準発振器の発振周波数の誤差や送信機と受信機の移動に伴うドップラーシフトにより周波数オフセットが生じる。この場合も受信機において周波数オフセットを推定し、補償する必要がある。
従来の周波数オフセットの推定には、既知のパイロットシンボルを用いる方法がある（非特許文献１参照）。しかしながらこの方法には、情報伝送に寄与しない既知のパイロットシンボルを送信信号に付加することによる伝送速度の低下や、既知のパイロットシンボルを検出するための回路や手順が必要となるという欠点がある。
一方、既知のパイロットシンボルを必要としない周波数オフセット推定方法には、1シンボル周期におけるシンボルの位相変化情報を利用する位相増加アルゴリズムや（非特許文献２参照）、周波数スペクトルを利用する方法（非特許文献３参照）が知られている。
図１７はM-PSK（M-Phase Shift Keying）変調信号に対して位相増加アルゴリズムを用いる従来の周波数オフセット推定装置の構成例を示すブロック図である。図１７に示される周波数オフセット推定装置は、１シンボル遅延部１０１、複素共役部１０２、乗算部１０３、M乗部１０４、加算部１０５、位相検出部１０６を備えている。
入力信号I+jQは、受信信号が所定のサンプリング周波数で予めサンプリングされた複素信号である。この入力信号は二分岐され、分岐された一方の信号は１シンボル遅延部１０１および複素共役部１０２を通過し、分岐された他方の信号と乗算部１０３で乗算されて、１シンボル間の位相変化情報を持った複素信号となる。この複素信号をM乗部１０４でM（正の整数）乗することにより、データの変調に起因する位相変化を取り除く。この位相変化を取り除いた信号を加算部１０５でN（正の整数）シンボルに渡って加算することにより、位相に関する平均化を行い、瞬時的な変化を取り除く。この加算後の信号から位相検出部１０６で位相を取り出し、さらに、M乗部１０４のM乗の演算によって１シンボル間の位相変化のM倍になっている位相を1/M倍にしている。この結果、周波数オフセットによって生じる１シンボル間の位相変化Δθを得る。周波数オフセット推定値Δfは次式によって算出される。ただし、この式においてRSはシンボルレートである。
Figure 2011158932
図１８は周波数スペクトルを利用する従来の周波数オフセット推定装置の構成例を示すブロック図である。図１８の周波数オフセット推定装置は、乗算部１０７、ＦＦＴ（fast Fourier transform）部１０８、周波数誤差検出部１０９、ＮＣＯ（numerically-controlled oscillator）１１０を備えている。
入力信号I+jQは、受信信号が所定のサンプリング周波数で予めサンプリングされた複素信号である。この入力信号はＮＣＯ１１０の出力信号と乗算部１０７で乗算されて、その周波数が変更される。この周波数が変更された信号はＦＦＴ部１０８に入力されて周波数領域の周波数スペクトルに変換される。周波数誤差検出部１０９は周波数スペクトルを測定し周波数誤差信号を出力する。この周波数誤差信号を元に、ＮＣＯ１１０はその出力信号の周波数を所定のステップで変化させる。以上のループの演算を周波数誤差信号がほぼ０となるまで繰り返し、周波数誤差信号がほぼ０となって収束した時点で周波数オフセット推定が完了する。
図１７に示される位相増加アルゴリズムは、コンスタレーション上で信号点が等位相間隔に配置されているM-PSK変調信号に対してのみ正確に動作し、直交振幅変調（QAM: quadrature amplitude modulation）で変調された信号に対しては動作しない。QAMで変調された信号に対して、既知のパイロットシンボルを必要としない周波数オフセット推定方法が非特許文献４や非特許文献５に開示されている。これら文献に記載されている推定方法を表す数式は以下の通りである。
Figure 2011158932
 ここで、y(p,t)は受信信号であり、偏波pと時間tの関数である。また、Nは推定に用いるシンボル数、RSはシンボルレートである。
Figure 2011158932
 ここで、数式２の周波数オフセットの推定範囲は、数式３により制限されている。
次に、この推定方法の動作について、図１９から図２２を参照して説明する。
図１９から図２２は、上述の非特許文献４や非特許文献５に係る周波数オフセット推定方法の動作を示す説明図である。
図１９に示すものは、64QAMで変調された信号に対して、０またはRS/4の整数倍の周波数オフセットを持っており、また位相オフセットが残留している場合のコンスタレーションである。信号点の周期はシンボルレートの逆数1/RSであり、周波数オフセットがRS/4であればコンスタレーション上の信号点は周波数オフセットが０の時の位置から丁度π/2回転したところに位置にする。すなわち周波数オフセットが０の時とRS/4の整数倍の周波数オフセットの時とでは同じコンスタレーション配置を示す。QAMで変調された信号はπ/2の位相対称性を持っているため、ある１つの信号点の位相を位相θとすると、原点からの距離がこの信号点と原点との距離に等しい点であって、且つこの信号点に対してπ/2の整数倍の位相差を持つ点が他に３点ある。図１９においては、これらの点を黒丸ｋ１からｋ４で示している。これら４点の黒丸ｋ１，ｋ２，ｋ３，ｋ４の位相βは、次の数式４で表される。
Figure 2011158932
 これらの信号を４乗した場合の位相は４βであり、次の数式５で表される。
Figure 2011158932
すなわち、黒丸の４点ｋ１，ｋ２，ｋ３，ｋ４は４乗すると複素平面上の同一点に収束する。同様に他の信号点も４乗すると原点からの距離が等しくπ/2の整数倍の位相差を持つ４点ごとに複素平面上の同一点に収束する。
図２０は、図１９のコンスタレーションの信号点を４乗した場合の信号点配置図である。図１９の６４点は図２０では１６点に収束している。これらの信号点は実軸（横軸）および虚軸（縦軸）に対して非対称であり、信号点の加算値または平均値は０ではない値を取る。
図２１に示すものは、64QAMで変調された信号に対して、RS/4の整数倍以外の周波数オフセットが有る場合のコンスタレーションである。原点からの距離が等しい信号点が円周上に配置されている。図２２は、図２１のコンスタレーションの信号点を４乗した場合の信号点配置図である。これらの信号点は実軸（横軸）および虚軸（縦軸）に対して対称であるため、信号点の加算値または平均値は０になる。
上述した数式２には、次に示す項が存在する。
Figure 2011158932
 この項を評価関数φc(f)とすると、これは、受信信号y(p,t)を４乗した後に、周波数fの逆回転演算子exp(-j2πft)を作用させて時間平均を取る操作である。逆回転演算子exp(-j2πft)の作用は、元の信号の周波数を−ｆだけ周波数変換する。従って、受信信号y(p,t)の周波数オフセットがfoの時、評価関数φc(f)の中の周波数fがf≠4fo+kRS(kは整数)の時に評価関数φc(f)は０となる。また、周波数fがf＝4fo+kRS (kは整数)の時に評価関数φc(f)は０ではない値を取る。
さらに、数式３の周波数オフセットの推定範囲の制限によってkRS(kは整数)の不確定性はなくなるため、評価関数φc(f)が最大値となるfを求めて1/4倍することにより、周波数オフセットfoを推定することができる。
M. K. Nezami et al., "DFT-based frequency acquisition algorithm for large carrier offsets in mobile satellite receivers," Electronics Letters, volume 37, pp. 386-387, March 2001. A. Leven et al., "Frequency estimation in intradyne reception," IEEE Photonics Technology Letters, volume 19, pp. 366-368, March 2007. K. Piyawanno et al., "Fast and accurate automatic frequency control for coherent receivers," ECOC2009, paper7.3.1, Sept., 2009. P. Ciblat et al., "Blind NLLS carrier frequency-offset estimation for QAM, PSK, and PAM modulations: performance at low SNR," IEEE Transactions on Communications, volume 54, pp. 1725-1730, Oct. 2006. M. Selmi et al., "Accurate digital frequency offset estimator for coherent PolMux QAM transmission systems," ECOC2009, paperP3.08, Sept.,2009.
しかしながら図１７に示す位相増加アルゴリズムを用いる従来の周波数オフセット推定装置では、位相の不確定性によって推定できる周波数オフセットの範囲が限定される。位相検出部１０６において、位相を取り出す演算「arg(・)」によって取り出すことのできる位相の範囲は[-π〜π]である。この取り出した位相をさらに1/M倍にするために、位相検出部１０６が検出することのできる位相の範囲は[-π/M〜π/M]である。従って図１７に示す周波数オフセット推定装置が推定できる周波数範囲は数式１によって、[-RS/2M〜RS/2M]に限定されるという問題がある。
また、図１８に示す周波数スペクトルを利用する従来の周波数オフセット推定装置では、ＮＣＯ１１０の出力信号の周波数を所定のステップで変えながら、それぞれの周波数で数百から数千のＦＦＴする際のデータ数に等しいサンプル数の入力信号を周波数オフセット推定装置に取り込む必要があるため、推定処理が収束するまでに時間がかかるという問題がある。さらに周波数スペクトルを利用する周波数オフセット推定装置では、伝送路中のバンドパスフィルタや受信機中のローパスフィルタのカットオフ周波数が小さく、かつ、周波数オフセットが大きい場合に周波数スペクトルの片側が削られて非対称となることにより、周波数オフセットの推定精度が劣化するという問題点がある。
また、数式２および数式３で表される、従来のQAMで変調された信号に対する周波数オフセットの推定方法（非特許文献４、非特許文献５）では、推定可能な周波数範囲が４乗した周波数帯で[-RS/2〜RS/2]、４乗する前の元の周波数帯で[-RS/8〜RS/8]に限定されるという問題がある。すなわち、この範囲を超えた周波数オフセットについては、誤って検出されることとなり問題である。さらに、位相雑音がある場合、コンスタレーション上の信号点が円周に沿って移動し、周波数オフセットが有る場合の図２１と同様の状態となるため正しく推定できないという問題がある。これらの問題について、図２３から図２５を参照して具体的に説明する。図２３から図２５は、従来の問題点を説明するためのシミュレーション結果の一例を示す説明図である。
なお、これらのシミュレーション結果においては、変調はシンボルレート28GBaudのDP-64QAM（偏波多重64直交振幅変調）、OSNR（光信号対雑音比）は10dB、推定に用いるシンボル数Nは1028であり、また周波数オフセットとして5GHzを与えて、周波数fと、次の数式６で表される評価関数φc(f)をプロットしている。
Figure 2011158932
図２３に示すものは、周波数fを-64GHzから64GHzと広帯域に掃引した場合のシミュレーション結果である。なお、本シミュレーションでは、位相雑音を付加していない。この図２３では、周波数オフセットの４倍である20GHzから、RS=28GHzの周波数間隔ごとに評価関数が大きくなる不確定性を示している。従って、非特許文献４や非特許文献５に記載された周波数オフセット推定方法では、評価関数のピークが常に２つ以上含まれない周波数範囲は[-RS/2〜RS/2]に制限されるという問題があることが分かる。
また、図２４および図２５に示すものは、レーザの線幅として10MHzを仮定してシミュレーションした結果である。なお、本シミュレーションでは、位相雑音を付加している。これらの図２４と図２５では、推定に用いている1028個のシンボル系列が互いに異なっている。図２４では、周波数オフセットの４倍である20GHzで最も大きな値を示しているが、他の周波数でも局所的なピークが多数見られる。図２５では、20GHzにもピークはあるが、20GHzとは異なる局所的なピークの方が全体の中で最も大きな値を示している。すなわち、図２５から評価関数φc(f)を最大とするfを求めると、誤った値を周波数オフセットとして推定することになる。このように、非特許文献４や非特許文献５に記載された周波数オフセット推定方法では、位相雑音や熱雑音の状態によって、しばしば誤った推定結果を出力するという問題がある。
本発明は、このような事情に鑑みてなされたもので、その目的は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する場合に、受信信号の周波数オフセットを適切に推定することができる、周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法を提供することにある。
この発明は上述した課題を解決するためになされたもので、本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換部と、前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限部と、前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御部と、を備える周波数オフセット推定装置である。
本発明の周波数オフセット推定装置において、前記周波数オフセット推定制御部は、前記周波数スペクトルの前記全ての周波数成分を前記周波数軸上で循環して移動させる場合に、前記正の周波数成分の電力が前記負の周波数成分の電力より大きい場合は、前記周波数スペクトルの全周波数成分を負の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号が１未満の場合は周波数成分番号にＮを加算し、前記正の周波数成分の電力が前記負の周波数成分の電力以下の場合は、前記周波数スペクトルの全周波数成分を正の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号がＮを越える場合は周波数成分番号からＮを減算するようにしてもよい。
本発明は、本発明の周波数オフセット推定装置と、前記周波数オフセット推定装置によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する第１の周波数オフセット補償部と、前記第１の周波数オフセット補償部によって補償された前記受信信号に対して、前記周波数オフセットを位相増加アルゴリズムに基づいて推定する、位相増加周波数オフセット推定部と、前記位相増加周波数オフセット推定部によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記周波数オフセットを補償する第２の周波数オフセット補償部と、を備える受信装置である。
本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する周波数オフセット推定装置において用いられる、周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換手順と、前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限手順と、前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御手順と、を備える周波数オフセット推定方法である。
本発明の周波数オフセット推定方法では、前記周波数オフセット推定制御手順において、前記周波数スペクトルの前記全ての周波数成分を前記周波数軸上で循環して移動させる場合に、前記正の周波数成分の電力が前記負の周波数成分の電力より大きい場合は、前記周波数スペクトルの全周波数成分を負の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号が１未満の場合は周波数成分番号にＮを加算し、前記正の周波数成分の電力が前記負の周波数成分の電力以下の場合は、前記周波数スペクトルの全周波数成分を正の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号がＮを越える場合は周波数成分番号からＮを減算することにより、前記周波数スペクトルを循環して移動するようにしてもよい。
本発明は、本発明の周波数オフセット推定方法による手順と、前記周波数オフセット推定方法による手順によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する第１の周波数オフセット補償手順と、前記第１の周波数オフセット補償手順によって補償された前記受信信号に対して、前記周波数オフセットを位相増加アルゴリズムに基づいて推定する、位相増加周波数オフセット推定手順と、前記位相増加周波数オフセット推定手順によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記周波数オフセットを補償する第２の周波数オフセット補償手順と、を備える受信方法である。
本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定部と、前記周波数オフセット粗推定部の粗推定値に基づいて掃引周波数の範囲を決定する掃引周波数範囲制御部と、前記掃引周波数範囲制御部により決定された前記掃引周波数の範囲において前記受信信号の周波数オフセットを推定する周波数オフセット精推定部と、を備え、前記周波数オフセット精推定部は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後、該受信信号に対して該受信信号の周波数から前記掃引周波数を減算する周波数変換を行う第１の演算部と、前記第１の演算部における各偏波の演算結果に対して、Ｎ（Ｎ：正の整数）シンボルの平均化または加算を行った後、絶対値または絶対値のべき乗の演算を行う第２の演算部と、前記第２の演算部の後段側において、前記２つの偏波の演算結果を加算、またはピーク値の大きい方の偏波の演算結果を選択するとともに、前記ＮシンボルからなるフレームについてＵ（Ｕ：正の整数）フレームの加算または平均化を行う第３の演算部と、前記第３の演算部の演算結果が最大値となる掃引周波数を検出し、当該掃引周波数を１／Ｗ倍して前記周波数オフセットを推定する第４の演算部と、を有する周波数オフセット推定装置である。
本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定手順と、前記周波数オフセット粗推定手順により推定された粗推定値に基づいて掃引周波数範囲を決定する掃引周波数範囲制御手順と、前記掃引周波数範囲制御手順により決定された前記掃引周波数の範囲において前記受信信号の周波数オフセットを推定する周波数オフセット精推定手順と、を備え、前記周波数オフセット精推定手順は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後、該受信信号に対して該受信信号の周波数から前記掃引周波数を減算する周波数変換を行う第１の演算手順と、前記第１の演算手順における各偏波の演算結果に対して、Ｎ（Ｎ：正の整数）シンボルの加算または平均化を行った後、絶対値または絶対値のべき乗の演算を行う第２の演算手順と、前記第２の演算手順の後において、前記２つの偏波の演算結果を加算、またはピーク値の大きい方の偏波の演算結果を選択するとともに、前記ＮシンボルからなるフレームについてＵ（Ｕ：正の整数）フレームの加算または平均化を行う第３の演算手順と、前記第３の演算手順の演算結果が最大値となる掃引周波数を検出し、当該掃引周波数を１／Ｗ倍して前記周波数オフセットを推定する第４の演算手順と、を含む周波数オフセット推定方法である。
本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定部と、前記受信信号または前記受信信号の分散が補償された信号に対して周期的な周波数オフセット推定特性を有する周波数オフセット精推定部と、前記周波数オフセット粗推定部によって推定された前記周波数オフセットの値に基づいて、前記周波数オフセット精推定部によって推定された周波数オフセットの周波数不確定性を除去し、前記周波数オフセットを推定する周波数不確定性除去制御部とを備える周波数オフセット推定装置である。
本発明の周波数オフセット推定装置において、前記周波数オフセット精推定部は、前記受信信号または前記受信信号の分散が補償された信号に対して位相増加アルゴリズムに基づいて前記周波数オフセットを推定するようにしてもよい。
本発明の周波数オフセット推定装置において、前記周波数オフセット精推定部は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後に周波数スペクトルに変換し、その変換結果に対して絶対値または絶対値のべき乗の演算を行い、これら２つの偏波の周波数スペクトルを加算、またはピーク値の大きい方の偏波の周波数スペクトルを選択するとともに、Ｎ（Ｎ：正の整数）シンボルからなるフレームの周波数スペクトルについてＵ（Ｕ：正の整数）フレームの加算または平均化を行い、この演算結果が最大値となる周波数を検出するようにしてもよい。
本発明の周波数オフセット推定装置において、前記周波数不確定性除去制御部は、前記周波数不確定性を除去し前記周波数オフセットを推定する際に、前記周波数オフセット精推定部によって推定された前記周波数オフセットに基づいて、前記周波数不確定性を含む前記周波数オフセットの候補となる周波数を算出し、周波数軸上で隣り合う前記周波数オフセットの候補となる周波数の中点を、隣り合う前記周波数オフセットの候補の境界とし、前記周波数軸上で前記境界に基づいた領域の中から、前記周波数オフセット粗推定部によって推定された前記値が含まれる領域を選択し、前記選択した領域に含まれる前記周波数オフセットの候補となる周波数を、前記周波数オフセットの推定値として選択するようにしてもよい。
本発明の周波数オフセット推定装置において、前記周波数オフセット粗推定部は、前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換部と、前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限部と、前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御部と、を備えるようにしてもよい。
本発明は、本発明の周波数オフセット推定装置と、前記周波数オフセット推定装置によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する周波数オフセット補償部と、前記周波数オフセット補償部によって補償された前記受信信号に対して、位相を補償する位相補償部と、前記位相を補償された前記受信信号のシンボルの判定を行う判定部と、を備える受信装置である。
本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定手順と、前記受信信号または前記受信信号の分散が補償された信号に対して周期的な周波数オフセット推定特性を有する周波数オフセット精推定手順と、前記周波数オフセット粗推定手順によって推定された前記周波数オフセットの値に基づいて、前記周波数オフセット精推定手順によって推定された周波数オフセットの周波数不確定性を除去し、前記周波数オフセットを推定する周波数不確定性除去制御手順と、を備える周波数オフセット推定方法である。
本発明の周波数オフセット推定方法において、前記周波数オフセット精推定手順は、前記受信信号または前記受信信号の分散が補償された信号に対して位相増加アルゴリズムに基づいて前記周波数オフセットを推定するようにしてもよい。
本発明の周波数オフセット推定方法において、前記周波数オフセット精推定手順は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後に周波数スペクトルに変換し、その変換結果に対して絶対値または絶対値のべき乗の演算を行い、これら２つの偏波の周波数スペクトルを加算、またはピーク値の大きい方の偏波の周波数スペクトルを選択するとともに、Ｎ（Ｎ：正の整数）シンボルからなるフレームの周波数スペクトルについてＵ（Ｕ：正の整数）フレームの加算または平均化を行い、この演算結果が最大値となる周波数を検出するようにしてもよい。
本発明の周波数オフセット推定方法において、前記周波数不確定性除去制御手順は、前記周波数不確定性を除去し前記周波数オフセットを推定する際に、前記周波数オフセット精推定手順によって推定された前記周波数オフセットに基づいて、前記周波数不確定性を含む前記周波数オフセットの候補となる周波数を算出し、周波数軸上で隣り合う前記周波数オフセットの候補となる周波数の中点を、隣り合う前記周波数オフセットの候補の境界とし、前記周波数軸上で前記境界に基づいた領域の中から、前記周波数オフセット粗推定手順によって推定された前記値が含まれる領域を選択し、前記選択した領域に含まれる前記周波数オフセットの候補となる周波数を、前記周波数オフセットの推定値として選択するようにしてもよい。
本発明の周波数オフセット推定方法において、前記周波数オフセット粗推定手順は、前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換手順と、前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限手順と、前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御手順と、を備えるようにしてもよい。
本発明は、本発明の周波数オフセット推定方法による手順と、前記周波数オフセット推定方法によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する周波数オフセット補償手順と、前記周波数オフセット補償手順によって補償された前記受信信号に対して、位相を補償する位相補償手順と、前記位相を補償された前記受信信号のシンボルの判定を行う判定手順と、を備える受信方法である。
この発明によれば、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する場合に、受信信号の周波数スペクトルの中の周波数帯域制限された正の周波数成分の電力値と負の周波数成分の電力値との電力差の絶対値が予め設定されたしきい値以下になるまで全周波数成分を移動し、その移動量に基づいて周波数オフセットを推定することにより、受信信号の周波数オフセットを適切に推定することができる。
すなわち、上記のように受信信号の周波数スペクトルの中の周波数帯域制限された正の周波数成分の電力値と負の周波数成分の電力値との電力差の絶対値が予め設定されたしきい値以下になるまで全周波数成分を移動し、その移動量に基づいて周波数オフセットを推定するため、位相増加アルゴリズムを用いる周波数オフセット推定装置と比較して、広帯域に周波数オフセットを推定することができる。また、既知のパイロットシンボルを必要としない方法であって、さらにＮＣＯも不要であるため、広帯域な他の周波数スペクトルを利用する周波数オフセット推定装置に比較して、高速かつ高精度に周波数オフセットを推定することができる。
また、この発明によれば、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する場合に、粗推定した周波数オフセットの値に基づいて、広帯域な周波数範囲において推定するため、受信信号の周波数オフセットを適切に推定することができる。
本発明の第１実施形態による周波数オフセット推定装置の構成例を示すブロック図である。 本第１実施形態の動作の概要を示す説明図であり、周波数オフセットがないときの周波数スペクトルを示す説明図である。 本第１実施形態の動作の概要を示す説明図であり、周波数オフセットがあるときの周波数スペクトルを示す説明図である。 本第１実施形態の動作をシミュレーションした結果の一例を示した説明図である。 従来の位相増加アルゴリズムを用いる周波数オフセット推定装置の動作をシミュレーションした結果の一例を示した説明図である。 本発明の第２実施形態による受信装置の構成例を示すブロック図である。 本発明の第３実施形態による受信装置の構成例を示すブロック図である。 本第３実施形態の動作の概要を示す説明図である。 第４の実施形態による周波数オフセット推定装置の構成例を示すブロック図である。 第４の実施形態による周波数オフセット推定装置の動作をシミュレーションした結果の一例を示す説明図である。 第５実施形態による周波数オフセット推定装置の構成例を示すブロック図である。 第６実施形態による受信装置の構成例を示すブロック図である。 第７実施形態による周波数オフセット推定装置の構成例を示すブロック図である。 第７実施形態による周波数オフセット推定装置が、周波数オフセットを推定する動作の概要を示す説明図である。 第７実施形態の実験結果の一例を示す説明図である。 第８実施形態による受信装置の構成例を示すブロック図である。 従来の位相増加アルゴリズムを用いる周波数オフセット推定装置の構成例を示すブロック図である。 従来の周波数スペクトルを利用する周波数オフセット推定装置の構成例を示すブロック図である。 従来の周波数オフセット推定方法の動作を示す第１の説明図である。 従来の周波数オフセット推定方法の動作を示す第２の説明図である。 従来の周波数オフセット推定方法の動作を示す第３の説明図である。 従来の周波数オフセット推定方法の動作を示す第４の説明図である。 従来の問題点を説明するためのシミュレーション結果の一例を示す第１の説明図である。 従来の問題点を説明するためのシミュレーション結果の一例を示す第２の説明図である。 従来の問題点を説明するためのシミュレーション結果の一例を示す第３の説明図である。
以下、本発明の実施形態を、図面を参照して説明する。なお、各図において、同一または対応する構成には同一の符号を用いている。
＜第１実施形態＞
以下、図面を参照して、本発明の第１実施形態について説明する。図１は、本発明の第１実施形態による、受信信号の搬送波周波数と受信側の局部発振器の出力信号の周波数との差である周波数オフセットを推定する周波数オフセット推定装置１５の構成例を示す概略ブロック図である。
周波数オフセット推定装置１５は、ＦＦＴ部１、ＳＰＤＴ（Single Pole, Double Throw）スイッチ２、バンドパスフィルタ３、および周波数オフセット推定制御部１２を備えている。さらに、周波数オフセット推定制御部１２は、第１の２乗加算部４、第２の２乗加算部５、減算部６、第１の判定部７、第２の判定部８、周波数スペクトル循環移動部９、カウンタ１０、および周波数スペクトル移動距離検出部１１を備えている。
図１において、入力信号I+jQは、受信信号が所定のサンプリング周波数で予めサンプリングされた複素信号である。
ＦＦＴ部1は、この入力信号を周波数領域の周波数スペクトルに変換する。例えば、ＦＦＴ部1がＦＦＴする際のデータ数をＮ（Ｎは任意の自然数）とすると、周波数の大きさの順に１からＮまで順序付けしたＮ個の番号の周波数成分を持つ周波数スペクトルを出力する。ここで、この周波数スペクトルの１からＮ/２までの周波数成分番号の周波数成分を負の周波数成分、周波数スペクトルのＮ/２＋１からＮまでの周波数成分番号の周波数成分を正の周波数成分としてこの後の信号処理を説明する。
ＳＰＤＴスイッチ２は、後述する第１の判定部７からの制御によりバンドパスフィルタ３に入力する信号を切り替える。ＳＰＤＴスイッチ２は、第１の判定部７の判定がＮＯの場合はＦＦＴ部1の出力信号を、ＹＥＳの場合は周波数スペクトル循環移動部９の出力信号を、それぞれバンドパスフィルタ３に入力するように切り替える。
初期設定は第１の判定部７はＮＯの判定で、ＳＰＤＴスイッチ２の共通端子は上側の接点に接続されており、ＳＰＤＴスイッチ２はＦＦＴ部1の出力信号を通過させてバンドパスフィルタ３に入力する。
バンドパスフィルタ３は、ＳＰＤＴスイッチ２を通過して入力された信号に対して周波数帯域制限をする。ここで、バンドパスフィルタ３は、ＦＦＴ部1で周波数スペクトルに変換された受信信号の正の周波数成分の信号と負の周波数成分の信号とに対して、１からＮの周波数成分における中心周波数に対して正と負が対称となるように周波数帯域制限をし、通過帯域が正の周波数成分の信号を第1の２乗加算部４へ、また通過帯域が負の周波数成分の信号を第２の２乗加算部５へ出力する。
第1の２乗加算部４および第２の２乗加算部５は、それぞれバンドパスフィルタ３の通過帯域の正の周波数成分の電力値PPおよび負の周波数成分の電力値PMを算出する。減算部６は、これら２つの電力値PPと電力値PMとの間の電力差を算出する。
ここで、図２および図３は、この周波数スペクトルに変換された受信信号に対して、バンドパスフィルタ３の通過帯域の正の周波数成分の電力値PPおよび負の周波数成分の電力値PMを算出する動作の概要を示す説明図である。図２と図３の横軸は周波数、縦軸は電力である。
バンドパスフィルタ３はＦＦＴ部1で周波数領域の周波数スペクトルに変換された信号に対して、その各周波数成分の信号にバンドパスフィルタ３の伝達関数を乗算する。例えばバンドパスフィルタ３は矩形波フィルタを用いた場合に、通過帯域の周波数成分には１、阻止帯域の周波数成分には０を乗算する。第1の２乗加算部４および第２の２乗加算部５は、通過帯域の周波数成分を２乗してその総和を取ることで、バンドパスフィルタ３の通過帯域の正負それぞれの周波数成分の電力の総和である電力値PPおよび電力値PMを算出する。
図２に示すように、周波数オフセットがない場合には、周波数スペクトルの対称性によって正の周波数成分の電力値PPと負の周波数成分の電力値PMは等しくなる。一方、図３に示すように、周波数オフセットがある場合には、周波数スペクトルは正の方向あるいは負の方向のいずれかに偏り、正の周波数成分の電力値PPと負の周波数成分の電力値PMとの間には電力差が生じる。
次に、図１の第１の判定部７は、電力値PPと電力値PMとの電力差の絶対値と、予め定められたしきい値THとを比較して大小の判定をする。第１の判定部７は、電力値PPと電力値PMの電力差の絶対値がしきい値THより大きく、ＹＥＳと判定した場合は、第２の判定部８の処理を開始させ、また、ＳＰＤＴスイッチ２の接点を切り替え、バンドパスフィルタ３への入力信号を、ＦＦＴ部1からの出力信号から、周波数スペクトル循環移動部９からの出力信号に変更する。
第２の判定部８は、電力値PPと電力値PMとを比較して大小の判定をする。第２の判定部８は、その判定した結果すなわち周波数スペクトルが正の方向あるいは負の方向のいずれに偏っているかを判定した結果に基づいて、周波数スペクトル循環移動部９およびカウンタ１０の中の移動制御を選択して通知する。周波数スペクトル循環移動部９は、第２の判定部８で判定した周波数成分の正の方向あるいは負の方向のいずれかの偏りが小さくなる方向に、全周波数成分を周波数軸上で移動させる。カウンタ１０は、正の方向に移動した場合と負の方向に移動した場合とのそれぞれの回数を計数し保持する。
ここで、全周波数成分を周波数軸上で移動するための制御について説明する。
この移動は、ＦＦＴ部１で変換された周波数スペクトルの周波数成分番号が１からＮとすると、移動後も周波数成分番号が１からＮになるような循環した移動とする。すなわち移動後の周波数成分番号が１未満の場合は周波数成分番号にＮを加算し、移動後の周波数成分番号がＮを越える場合は周波数成分番号からＮを減算する。
第２の判定部８において電力値PPが電力値PMより大きいと判定された場合、周波数スペクトル循環移動部９は、負の方向に周波数成分を１ずつ循環して移動させ、カウンタ１０の値dを１増やす。反対に、第２の判定部８において電力値PPが電力値PMより小さいと判定された場合は、周波数スペクトル循環移動部９は、正の方向に周波数成分を１ずつ循環して移動させ、カウンタ１０の値dを１減らす。
また、このときＳＰＤＴスイッチ２の接点は、周波数スペクトル循環移動部９からの出力信号がバンドパスフィルタ３に入力されるように変更されているため、周波数スペクトル循環移動部９によって循環移動された周波数スペクトルが、再びバンドパスフィルタ３に入力され、第1の２乗加算部４、第２の２乗加算部５、および減算部６の順に信号処理が実行され、電力値PPと電力値PMとの電力差が算出され第１の判定部７に入力される。このループの信号処理は、第１の判定部７がＮＯと判定するまで、繰り返し実行される。すなわち、このループの信号処理は、電力値PPと電力値PMとの電力差の絶対値が予め定められたしきい値TH以下と判定されるまで繰り返し実行される。
第１の判定部７がＮＯと判定した場合にこのループの信号処理は収束し、周波数スペクトル移動距離検出部１１は、収束時点でのカウンタ１０の値dの値を検出する。この値は周波数スペクトルの周波数軸上での移動距離の総和であり、この移動距離の総和である移動量に対応する周波数に基づいて周波数オフセットを推定する。
以上で周波数オフセット推定は終了し、第１の判定部７は、ＳＰＤＴスイッチ２の接点をＦＦＴ部１の出力信号がバンドパスフィルタ３に入力されるように切り替え、ＦＦＴ部１から新たに周波数スペクトルデータをこのループ信号処理に入力し、周波数オフセットの推定は開始される。
図４は本第１実施形態の動作をシミュレーションした結果の一例を示した説明図であり、図５は比較のために、従来の位相増加アルゴリズムを用いる周波数オフセット推定装置の動作をシミュレーションした結果の一例を示した説明図である。
図４、図５共に、横軸は周波数オフセットの設定値、縦軸は周波数オフセット推定結果である。変調方式は112Gbit/s偏波多重ＱＰＳＫ（Quadrature Phase Shift Keying）であり、シンボルレートRSは28Gbaudである。
図５に示す従来の位相増加アルゴリズムを用いる周波数オフセット推定装置では、±RS/2M=±3.5GHz（位相数M＝４）の範囲内で精度良く推定できているが、その範囲を越えると位相の不確定性により誤った推定結果となっている。一方、図４に示す本第１実施形態では、±5GHzの広帯域に渡って精度良く推定できていることが分かる。
図１から図４を用いて説明したように、本第１実施形態の周波数オフセット推定装置１５は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する場合に、受信信号の周波数スペクトルの中の周波数帯域制限された正の周波数成分の電力値PPと負の周波数成分の電力値PMとの電力差の絶対値が予め設定されたしきい値以下になるまで全周波数成分を循環して移動し、その移動量に基づいて周波数オフセットを推定する。
これにより、本第１実施形態の周波数オフセット推定装置１５は、図１７に示した従来の位相増加アルゴリズムを用いる周波数オフセット推定装置が推定できる周波数オフセットの範囲に比較して、広帯域に周波数オフセットを推定することができる。また、本第１実施形態による方法は、既知のパイロットシンボルを必要としない方法であって、さらに、図１８に示した従来の周波数オフセット推定装置で用いていたＮＣＯは不要な方法である。そのため本実施形態によれば、ステップ的な周波数変化の各周波数で数百から数千のＦＦＴする際のデータ数に等しいサンプル数の入力信号を取り込む必要がない。よって本実施形態によれば、周波数オフセットを高速かつ高精度に推定することができる。
また、本第１実施形態の周波数オフセット推定装置１５は、周波数成分を循環して移動する大きさは１に限る必要はない。そのため、電力値PPと電力値PMとの電力差の絶対値が大きい時には循環して移動する大きさを大きくし、電力値PPと電力値PMとの電力差の絶対値がしきい値THに近づいた時に循環して移動する大きさを小さくすることができる。これにより、さらに周波数オフセットを高速に推定することができる。
なお、第１の判定部７の予め定めるしきい値THの値は、入力信号の電力によって決めれば良い。例えばＦＦＴ部１の出力信号の周波数スペクトルが一様な矩形であると仮定し、また周波数スペクトルの全電力の総和を電力値PTOTALとして周波数成分を負の方向に１だけ循環移動すると、電力値PPは移動前に比べてPTOTAL/N 小さくなり、電力値PMはPTOTAL/N 大きくなる。従って電力値PPと電力値PMとの差は2PTOTAL/N変化する。このとき、しきい値THを2PTOTAL/N程度とすれば、ほぼ1周波数成分の精度で周波数オフセットを推定することができる。推定精度を犠牲にしても高速な推定処理時間が要求される場合には、しきい値THを大きくすれば良い。
なお、ＦＦＴ部１がＦＦＴする際のデータ数Ｎは、要求される推定精度と推定処理時間とによって決めることができる。一般にＮが小さいほど推定精度は劣化するが、推定処理時間は短縮される。
また、周波数スペクトル移動距離検出部１１より出力される周波数オフセット推定値は、複数個を蓄積した後に相加平均しても良いし、忘却係数を用いた平均化を行って時間的な変動に追随させても良い。
なお、入力信号を周波数スペクトルに変換するＦＦＴ部１は、周波数解析方法として高速フーリエ変換に限らず、短時間フーリエ変換、離散フーリエ変換、ウエーブレット変換、あるいはその他の周波数解析方法を用いても良い。
また、バンドパスフィルタ３は、周波数帯域を制限するフィルタとして、バンドパスフィルタに限らず、受信する信号の特性に基づいて、ローパスフィルタ、ハイパスフィルタ、あるいはその他の周波数帯域を制限するフィルタを用いても良い。
＜第２実施形態＞
以下、図面を参照して、本発明の第２実施形態について説明する。本第２実施形態は、第１実施形態に記載した周波数オフセット推定装置１５により周波数オフセットを粗く推定し、その推定値で周波数オフセットを補償した後に、従来の位相増加アルゴリズムを用いる周波数オフセット推定装置で周波数オフセットを推定し、さらに補償するという２種類の周波数推定装置を組み合わせた構成を持つ受信装置である。
図６は、本発明の第２実施形態による受信装置の構成例を示す概略ブロック図である。
受信装置は、周波数オフセット推定装置１５、第１の周波数オフセット補償部２８、偏波分離部１９、位相増加周波数オフセット推定装置２０、第２の周波数オフセット補償部２９、第１の位相補償部２４、第２の位相補償部２５、第１の判定部２６、および第２の判定部２７を備えている。ここで、位相増加周波数オフセット推定装置２０は、図１７に構成例を示した位相増加アルゴリズムを用いる従来の周波数オフセット推定装置と同様の構成を有している。図６に示す周波数オフセット推定装置１５は、図１に示す第１実施形態に記載の構成に対応しており、その説明を省略する。
第１の周波数オフセット補償部２８は、第1のＮＣＯ１６、第1の乗算部１７、および第２の乗算部１８を備えている。
第２の周波数オフセット補償部２９は、第２のＮＣＯ２１、第３の乗算部２２、および第４の乗算部２３を備えている。
入力信号I+jQおよび入力信号I’+jQ’はそれぞれX偏波とY偏波の受信信号が所定のサンプリング周波数で予めサンプリングされた複素信号である。
本第２実施形態では、まず初めに、周波数オフセット推定装置１５は、入力信号I+jQおよび入力信号I’+jQ’に対して、周波数オフセットを広帯域に粗く推定する。ただし、周波数オフセット推定装置１５は、次のような条件を満たすように設計されている。その条件とは、周波数オフセット推定装置１５がこの周波数オフセットを粗く推定する際は、この粗く推定された周波数オフセットを補償した入力信号に対して、位相増加周波数オフセット推定装置２０が周波数オフセットを推定する場合に、周波数の不確定性なく推定できる周波数の範囲の中に入るような条件である。
次に、周波数オフセット推定装置１５で粗く推定された周波数オフセット推定値は、第１の周波数オフセット補償部２８に入力される。そして第１の周波数オフセット補償部２８は、入力信号の周波数オフセットを補償する。
つまり、粗く推定された周波数オフセットの推定値に基づいて、第１の周波数オフセット補償部２８が有する第1のＮＣＯ１６の発振周波数は調整される。そして、入力信号I+jQおよび入力信号I’+jQ’と、調整された第1のＮＣＯ１６の発振周波数を持った信号とが、第１の周波数オフセット補償部２８が有する第1の乗算部１７および第２の乗算部１８において乗算され、入力信号I+jQおよび入力信号I’+jQ’の周波数オフセットは補償される。
偏波分離部１９は、第１の周波数オフセット補償部２８で入力信号I+jQおよび入力信号I’+jQ’の周波数オフセットを補償された信号に対して、偏波分離および残留分散の除去を実行する。位相増加周波数オフセット推定装置２０は、偏波分離部１９から出力された信号を入力し、周波数オフセットの推定を実行する。
この、位相増加周波数オフセット推定装置２０で推定した周波数オフセット推定値は、第２の周波数オフセット補償部２９に入力される。そして第２の周波数オフセット補償部２９は、入力信号の周波数オフセットを補償する。
つまり、位相増加周波数オフセット推定装置２０で推定した周波数オフセット推定値に基づいて、第２の周波数オフセット補償部２９が有する第２のＮＣＯ２１の発振周波数は調整される。そして、第２の周波数オフセット補償部２９への入力信号と、調整された第２のＮＣＯ２１の発振周波数を持った信号とが、第２の周波数オフセット補償部２９が有する第３の乗算部２２および第４の乗算部２３において乗算され、第２の周波数オフセット補償部２９への入力信号の周波数オフセットは補償される。
次に、第２の周波数オフセット補償部２９において周波数オフセットが補償された信号に対して、第１の位相補償部２４および第２の位相補償部２５は位相を補償し、それぞれ第１の判定部２６および第２の判定部２７へ出力する。第１の判定部２６および第２の判定部２７は、シンボルの判定をして復調信号を生成する。
このように、本第２実施形態の周波数オフセットの推定装置を備える受信装置は、まず周波数オフセット推定装置１５で周波数オフセットを粗く推定した値に基づいて入力信号を補償してから、従来の位相増加周波数オフセット推定装置２０に入力する。また、偏波分離および残留分散の除去を合わせて実行することで、周波数オフセットおよび分散値が共に小さい入力信号に対して、位相増加アルゴリズムを用いて周波数オフセットを推定することができる。これにより、本第２実施形態の周波数オフセットの推定装置を備える受信装置は、受信信号に対して、広帯域かつ高精度に周波数オフセットを推定して補償し、復調することができる。
ところで、周波数スペクトルを利用する周波数オフセット推定装置において、伝送路中のバンドパスフィルタや受信機中のローパスフィルタのカットオフ周波数が小さく、かつ周波数オフセットが大きい場合に周波数スペクトルの片側が削られて非対称となることがある。これに対して、本第２実施形態の周波数オフセットの推定装置においては、周波数オフセット推定装置１５で周波数オフセットを粗く推定し補償した後に、さらに従来の位相増加周波数オフセット推定装置２０で周波数オフセットを推定する。そのため、本第２実施形態の周波数オフセットの推定装置においては、伝送路中のバンドパスフィルタや受信機中のローパスフィルタのカットオフ周波数が小さく、かつ周波数オフセットが大きい場合であっても、周波数オフセットの推定精度が劣化してしまうことを防ぐことができる。
なお、位相増加周波数オフセット推定装置２０に入力される信号は、偏波分離部１９からの出力信号に代えて、偏波分離部１９に入力される前の、第１の周波数オフセット補償部２８で周波数オフセットを補償された信号であっても良い。ただし、偏波分離部１９において偏波分離および残留分散の除去をされた信号を、位相増加周波数オフセット推定装置２０へ入力した方が、より高精度に周波数オフセットを推定できる。
＜第３実施形態＞
以下、図面を参照して、本発明の第３実施形態について説明する。
上記に説明したように、第２実施形態の受信装置においては、周波数オフセット推定装置１５が粗く推定した値に基づいて入力信号を補償している。これに対して、本第３実施形態の受信装置は、周波数オフセット推定装置１５が粗く推定した値に基づいて、入力信号は補償せずに、従来の位相増加アルゴリズムを用いて検出する際の位相の不確定性を除去して周波数オフセットを推定する構成を持つ受信装置である。なお、本実施形態における受信装置では、第１実施形態の周波数オフセット推定装置に代えて、非特許文献３に記載の発明を利用することも可能である。
図７は、本発明の第３実施形態による受信装置の構成例を示す概略ブロック図である。
受信装置は、周波数オフセット推定装置１５、周波数不確定性除去制御部３０、偏波分離部１９、位相増加周波数オフセット推定装置２０、周波数オフセット補償部２９、第１の位相補償部２４、第２の位相補償部２５、第１の判定部２６、および第２の判定部２７を備えている。周波数オフセット補償部２９は、ＮＣＯ２１、第１の乗算部２２、および第２の乗算部２３を備えている。
周波数オフセット推定装置１５は、第１実施形態に記載の構成であり説明を省略する。
図７において図６の各部に対応する構成には同一の符号をつけており、その説明を省略する。
図７において、入力信号I+jQおよび入力信号I’+jQ’はそれぞれX偏波とY偏波の受信信号が所定のサンプリング周波数で予めサンプリングされた複素信号である。
本第３実施形態では、まず、周波数オフセット推定装置１５は、入力信号I+jQおよび入力信号I’+jQ’に対して、周波数オフセットを広帯域に粗く推定し、その推定値を周波数不確定性除去制御部３０に出力する。
また、位相増加周波数オフセット推定装置２０は、偏波分離部１９から出力された信号に対して周波数オフセットを推定し、その推定値を周波数不確定性除去制御部３０に出力する。
次に、周波数不確定性除去制御部３０は、周波数オフセット推定装置１５、および位相増加周波数オフセット推定装置２０から入力されたそれぞれの周波数オフセットの推定値に基づいて、周波数オフセットを推定し、周波数オフセット補償部２９のＮＣＯ２１に出力し、ＮＣＯ２１の発振周波数を調整する。
周波数オフセット補償部２９で、入力信号I+jQおよび入力信号I’+jQ’の周波数オフセットが補償されたあと、第１の位相補償部２４、第２の位相補償部２５、第１の判定部２６、および第２の判定部２７で復調される。
ここで、周波数不確定性除去制御部３０が、位相増加周波数オフセット推定装置２０の周波数オフセット推定結果の不確定性を除去し、周波数オフセットを推定する際の動作について説明する。
まず、周波数不確定性除去制御部３０は、位相増加周波数オフセット推定装置２０から入力された周波数オフセットの推定値に基づいて、複数の周波数オフセットの候補の周波数を算出する。次に、周波数不確定性除去制御部３０は、算出した複数の周波数オフセットの候補において、周波数軸上で隣り合う候補と候補の中点を算出し境界とする。さらに、周波数不確定性除去制御部３０は、周波数オフセット推定装置１５から入力された周波数オフセットの推定値が、上記周波数オフセットの候補の隣り合う境界と境界の間の領域にあることを検出する。そして、周波数不確定性除去制御部３０は、上記検出された領域の中にある周波数オフセットの候補を、周波数オフセットの推定値として選択する。
次に、図８を用いて、本第３実施形態における周波数オフセットを推定する動作の概要について説明する。
周波数オフセット推定装置１５の推定値に基づいて、周波数不確定性除去制御部３０は、位相増加周波数オフセット推定装置２０の周波数オフセット推定結果の不確定性を除去する。
この周波数オフセット推定結果の不確定性は、位相増加周波数オフセット推定装置２０において、位相増加アルゴリズムで検出する位相の不確定性に起因するものであり、周期性を持っている。位相増加周波数オフセット推定装置２０の周波数オフセット推定値を周波数fMthとすると、推定値が周波数fMthとなる可能性のある周波数オフセットの候補の周波数fCND(k)は、次の数式７で与えられる。
Figure 2011158932
実際の周波数オフセットを周波数fCND(Ｌ)とする。すなわち不確定性を表す係数kをk=Ｌとする（Ｌは整数）。この時、周波数オフセット推定装置１５の周波数オフセット推定値を周波数fCOとすると、周波数fCOが次の数式８の関係式を満足するように周波数オフセット推定装置１５を設計することで、不確定性を除去することができる。
Figure 2011158932
図８は、周波数軸上に周波数fMthからRS/M周期に周波数fCND(k)が並んでおり、隣り合う周波数fCND(k)の周波数の中点を、隣り合う周波数オフセットの候補の境界としていることを示している。図８では、周波数fCOが周波数オフセットの候補の周波数fCND(1)の領域にあるため、周波数不確定性除去制御部３０は、実際の周波数オフセットをfCND(1)と推定する。
このように、本第３実施形態の周波数オフセットの推定装置を備える受信装置は、まず周波数オフセット推定装置１５で周波数オフセットを粗く推定する。次に、周波数不確定性除去制御部３０は、その推定値に基づいて位相増加アルゴリズムで検出する位相の不確定性を除去するとともに、周波数オフセットを推定する。これにより、本第３実施形態の周波数オフセットの推定装置を備える受信装置は、受信信号に対して、広帯域、かつ高精度に周波数オフセットを推定して補償し、復調することができる。また第２実施形態と比較すると、周波数オフセット推定装置１５の推定値で周波数オフセットを補償するためのＮＣＯや乗算器が不要となり、回路規模を低減することができる。
なお、位相増加周波数オフセット推定装置２０に入力される入力信号は、本第３実施形態の構成例の図７では、偏波分離部１９からの出力信号であるが、偏波分離部１９に入力される前の信号であっても良い。
＜第４実施形態＞
以下、図面を参照して、本発明の第４実施形態について説明する。
図９は、本発明の第４の実施形態による周波数オフセット推定装置１５０の構成例を示すブロック図である。
この周波数オフセット推定装置は、周波数オフセット精推定部４１、周波数オフセット粗推定部５７、および掃引周波数範囲制御部５８を備えている。さらに、周波数オフセット精推定部４１は、第１の４乗回路４２、第２の４乗回路４３、第１の乗算回路４４、第２の乗算回路４５、ＮＣＯ（numerically-controlled oscillator）４６、第1のＮシンボル加算回路４７、第２のＮシンボル加算回路４８、第１の１／Ｎ除算回路４９、第２の１／Ｎ除算回路５０、第１の絶対値２乗回路５１、第２の絶対値２乗回路５２、偏波加算回路５３、Ｕフレーム加算回路５４、最大値掃引周波数検出回路５５、および１／４除算回路５６を備えている。
周波数オフセット粗推定部５７は、所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する。周波数オフセット粗推定部５７には、第２実施形態と同様、第１実施形態に記載の周波数スペクトルを利用する周波数オフセット推定装置１５を用いることができる。また、第１実施形態に記載の周波数オフセット推定装置１５に代えて、非特許文献３に記載の周波数オフセット推定装置も同様に用いることができる。掃引周波数範囲制御部５８は、周波数オフセット粗推定部５７の粗推定値に基づいて掃引周波数範囲を決定する。周波数オフセット精推定部４１は、掃引周波数範囲制御部５８により決定された掃引周波数範囲において受信信号の周波数オフセットを推定する。ここで、受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗとする。なお、Ｗは変調方式によって定まる数（正の整数）であり、QAM変調であればＷ＝４、PAM（pulse amplitude modulation）変調であればＷ＝２、M-PSK変調であればＷ＝Ｍである。本第４実施形態においては、Ｗ＝４の場合として、説明する。
周波数オフセット精推定部４１は、第１の演算部４１ａ、第２の演算部４１ｂ、第３の演算部４１ｃ、および第４の演算部４１ｄを備えている。第１の演算部４１ａは、受信信号における２つの偏波をそれぞれ４乗した後、該受信信号に対して該受信信号の周波数から前記掃引周波数を減算する周波数変換を行う。そして、第２の演算部４１ｂは、第１の演算部４１ａにおける各偏波の演算結果に対して、Ｎ（Ｎ：正の整数）シンボルの平均化（平均値の算出）を行った後、絶対値のべき乗の演算（平均値に対して絶対値のべき乗を算出）を行う。さらに、第３の演算部４１ｃは、第２の演算部４１ｂによる２つの偏波の演算結果を加算した後、ＮシンボルからなるフレームについてＵ（Ｕ：正の整数）フレームの加算処理を行う。そして、第４の演算部４１ｄは、第３の演算部４１ｃの演算結果が最大値となる掃引周波数を検出し、この掃引周波数を１／４倍して出力して周波数オフセットを推定する。
なお、第１の演算部４１ａは、第１の４乗回路４２、第２の４乗回路４３、第１の乗算回路４４、第２の乗算回路４５、およびＮＣＯ４６を備えている。第２の演算部４１ｂは、第1のＮシンボル加算回路４７、第２のＮシンボル加算回路４８、第１の１／Ｎ除算回路４９、第２の１／Ｎ除算回路５０、第１の絶対値２乗回路５１、および第２の絶対値２乗回路５２を備えている。また、第３の演算部４１ｃは、偏波加算回路５３およびＵフレーム加算回路５４を備えている。第４の演算部４１ｄは、最大値掃引周波数検出回路５５および１／４除算回路５６を備えている。
以下、本実施形態における周波数オフセット精推定部４１の具体的な動作例を説明する。
この第４実施形態の動作を表す数式を数式９および数式１０に示す。
Figure 2011158932
Figure 2011158932
 ここで、Ｎは1フレームの推定に用いるシンボル数、uはＮ個のシンボル系列からなるフレームのフレーム番号、Ｕは推定に用いる総フレーム数、y(u,p,t)は受信信号であり、フレームu、偏波p、および時間tの関数である。また、RSはシンボルレートであり、fcoarseは周波数オフセット粗推定部５７の粗推定周波数である。
図９において、入力信号I+jQおよびI'+jQ'は、それぞれ、数式９の受信信号y(u,p,t)が所定のサンプリング周波数で予めサンプリングされたX偏波およびY偏波の複素信号である。周波数オフセット精推定部４１は、この入力信号を第１の４乗回路４２および第２の４乗回路４３で４乗し、ＮＣＯ４６の出力信号exp(-j2πft)と第１の乗算回路４４および第２の乗算回路４５において乗算する。この乗算は、４乗した受信信号の周波数から周波数ｆを減算する周波数変換を行う演算に相当する。
ＮＣＯ４６の掃引周波数は、掃引周波数範囲制御部５８により制御される。具体的には、掃引周波数範囲制御部５８は、周波数オフセット粗推定部５７により推定された粗推定周波数fcoarseとシンボルレートRSから数式１０を満足するようにＮＣＯ４６の掃引周波数の周波数範囲を決め、数式１０の下限周波数から上限周波数まで所定のステップで周波数ｆを変化させる。周波数オフセット精推定部４１は、各々の周波数ｆにおいて、第１のＮシンボル加算回路４７および第２のＮシンボル加算回路４８でのＮ個のシンボルに対する演算結果を加算し、第１の１／Ｎ除算回路４９および第２の１／Ｎ除算回路５０においてＮで除算してＮ個の平均値とし、第１の絶対値２乗回路５１および第２の絶対値２乗回路５２において絶対値の２乗を求め、偏波加算回路５３において２つの偏波に対する演算結果を加算する。ここで、周波数オフセット精推定部４１は、Ｎ個のシンボル系列を１フレーム目、次のＮ個のシンボル系列を２フレーム目とみなし、Ｕフレームの演算結果をＵフレーム加算回路５４において加算する。当該加算処理により、周波数オフセット精推定部４１は雑音成分を除去する。ここまでの演算により、次の数式１１で示す評価関数φn(f)が求まる。
Figure 2011158932
周波数オフセット推定装置１５０は、この評価関数φn(f)が最大値となる周波数fを最大値掃引周波数検出回路５５で算出し、1/4除算回路５６で1/4倍することにより、周波数オフセットを推定することができる。
本第４実施形態では、周波数オフセット粗推定部５７と掃引周波数範囲制御部５８とがＮＣＯ４６の掃引周波数範囲を制御することにより、周波数オフセットの推定範囲を広帯域にすることができる。この点について、図２３を用いて説明する。
例えば、周波数オフセットが5GHzである図２３において、粗推定周波数fcoarseが4GHzであったとすると、シンボルレートRSは28GHzであることから、数式１０よりＮＣＯ４６の掃引周波数範囲は[2GHz〜30GHz]となる。このため、周波数オフセット推定装置１５０は、20GHzのピークのみを検出してその両隣のピークは検出しないため、不確定性が無く20GHzを最大値掃引周波数として検出し、その1/4である5GHzを周波数オフセット推定値として出力することができる。
本実施形態の周波数オフセット推定装置１５０は、推定可能な周波数範囲が [fcoarse-RS/8〜fcoarse+RS/8]であり、また粗推定周波数fcoarseが可変であることから、広帯域に推定が可能である。また本実施形態では、真の周波数オフセット値と粗推定周波数fcoarseとの誤差には±RS/8と大きな値が許容できるため、周波数オフセット粗推定部５７には高精度な特性が要求されない。なお、従来の技術では、推定範囲が[-14GHz〜14GHz]となるため、-8GHzを最大値掃引周波数として検出し、その1/4である-2GHzを周波数オフセット推定値として出力することとなる。
図１０は、図２４および図２５と同じ周波数オフセット5GHz、レーザの線幅10MHz、OSNR10dBを仮定し、1フレームの推定に用いるシンボル数Ｎを１０２８、推定に用いる総フレーム数Ｕを８０としてシミュレーションを行って、数式１１の評価関数φn(f)を求めた結果を示している。
この図１０を図２４および図２５と比べると、周波数オフセットの４倍である20GHzのピークが強調され、フレームごとに一時的に現れていた局所的なピークは相対的に低く抑えられている。従って、周波数オフセット推定装置１５０は、位相雑音や熱雑音によりフレームごとに一時的に現れていた局所的なピークに影響されることなく、正しい推定結果を算出することが可能となる。
以上のように、本第４実施形態の周波数オフセット推定装置１５０は、周波数オフセット粗推定部５７が粗く推定した結果に基づいて、周波数オフセット精推定部４１が周波数オフセットを推定する際の周波数不確定性を除去するように、掃引周波数範囲制御部５８がＮＣＯ４６の掃引周波数の周波数範囲を決め、周波数オフセットを推定する。また、周波数オフセット推定装置１５０は、周波数オフセット精推定部４１が周波数オフセットを推定する演算処理において、受信信号の雑音成分を除去する。これにより、本第４実施形態の周波数オフセット推定装置１５０は、受信信号に対して広帯域かつ高精度に周波数オフセットを推定することができる。
＜第５実施形態＞
以下、図面を参照して、本発明の第５実施形態について説明する。
図１１は、本第５実施形態による周波数オフセット推定装置１５０の構成例を示すブロック図である。なお、図１１において図９の各部に対応する構成には同一の符号をつけており、その説明を省略する。
なお、本第５実施形態による周波数オフセット精推定部４１において、第１の演算部４１ａ、および第４の演算部４１ｄは、第４実施形態による第１の演算部４１ａ、および第４の演算部４１ｄと同様の構成である。一方、本第５実施形態による周波数オフセット精推定部４１において、第２の演算部４１ｂ、および第３の演算部４１ｃは、第４実施形態による第２の演算部４１ｂ、および第３の演算部４１ｃと異なる構成である。本第５実施形態による周波数オフセット精推定部４１において、第２の演算部４１ｂは、第1のＮシンボル加算回路４７、第２のＮシンボル加算回路４８、第１の絶対値回路５９、および第２の絶対値回路６０を備えている。また、第３の演算部４１ｃは、第１のＵフレーム加算回路６１,第２のＵフレーム加算回路６２、第１の１／Ｕ除算回路６３、第２の１／Ｕ除算回路６４、および偏波加算回路５３を備えている。
本実施形態において周波数オフセット精推定部４１は、第２の演算部４１ｂにおいて、第１の演算部４１aにおける各偏波の演算結果に対して、Ｎシンボルの加算を行った後、絶対値の演算（加算した値に対して絶対値を算出）を行う。さらに、周波数オフセット精推定部４１は、第３の演算部４１ｃにおいて、第２の演算部４１ｂで得られた演算結果に対して、ＮシンボルからなるフレームについてＵフレームの平均化処理を行った後、２つの偏波の演算結果を加算する。そして、周波数オフセット精推定部４１は、第４の演算部４１ｄにおいて、第３の演算部４１ｃの演算結果が最大値となる掃引周波数を検出し、この掃引周波数を１／４倍して出力して周波数オフセットを推定する。
本実施形態における周波数オフセット精推定部４１は、図９に示す第４実施形態の構成と比較して、第１の１／Ｎ除算回路４９および第２の１／Ｎ除算回路５０が削除されている。このように第１の１／Ｎ除算回路４９および第２の１／Ｎ除算回路５０が削除されることにより、評価関数φn(f)の各周波数成分はＮ倍になるものの、周波数間の相対的な大小関係は変わらない。よって、本実施形態の周波数オフセット精推定部４１は、第４実施形態の周波数オフセット精推定部４１と同様に周波数オフセットを推定する。ここで、Ｎシンボル加算回路と１／Ｎ除算回路との両方により処理される場合は平均化処理、１／Ｎ除算回路を省いてＮシンボル加算回路により処理される場合は加算処理となる。
また、第４実施形態の周波数オフセット精推定部４１における第１の絶対値２乗回路５１および第２の絶対値２乗回路５２は、本実施形態の周波数オフセット精推定部４１においては、第１の絶対値回路５９および第２の絶対値回路６０に置き換えられている。このように２乗を１乗や他のべき乗に変えても、評価関数φn(f)の各周波数成分の相対的な大小関係は変わらない。そのため、本実施形態の周波数オフセット精推定部４１は、第４実施形態の周波数オフセット精推定部４１と同様に周波数オフセットを推定する。
また、第４実施形態の周波数オフセット精推定部４１の構成と比較して、本実施形態の周波数オフセット精推定部４１においては、第１の１／Ｕ除算回路６３および第２の１／Ｕ除算回路６４が追加されている。これによって評価関数φn(f)の各周波数成分は１／Ｕ倍になるものの、周波数間の相対的な大小関係は変わらない。そのため、本実施形態の周波数オフセット精推定部４１は、第４実施形態の周波数オフセット精推定部４１と同様に周波数オフセットを推定する。ここで、Ｕフレーム加算回路と１／Ｕ除算回路との両方により処理される場合は平均化処理、１／Ｕ除算回路を省いて、Ｕフレーム加算回路により処理される場合は加算処理となる。
さらに、第４実施形態の周波数オフセット精推定部４１における偏波加算回路５３からＵフレーム加算回路５４への順番は、本実施形態の周波数オフセット精推定部４１においては、第１のＵフレーム加算回路６１（または、第２のＵフレーム加算回路６２）から、１／Ｕ除算回路を介して偏波加算回路５３への順番に入れ替えられている。このようにフレーム加算と偏波加算との順番が入れ替えられても評価関数φn(f)の各周波数成分の相対的な大小関係は変わらない。そのため、本実施形態の周波数オフセット精推定部４１は、第４実施形態の周波数オフセット精推定部４１と同様に周波数オフセットを推定する。なお、推定に用いる総フレーム数Ｕが大きければ、偏波加算を行わず、ピーク値の大きい方の偏波の評価関数のみを選択することも可能である。
これにより、本第５実施形態の周波数オフセット推定装置１５０は、上述の周波数オフセット精推定部４１の構成においても、第４実施形態の周波数オフセット推定装置１５０の周波数オフセット精推定部４１の構成と同様に、受信信号に対して広帯域かつ高精度に周波数オフセットを推定することができる。
また、第４実施形態や第５実施形態はQAMで変調された信号に限るものではなく、PSKやPSKの振幅を多段階とした変調方式、あるいはPAM変調方式についても適用できる。この場合、周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗとすると、第１の４乗回路４２および第２の４乗回路４３をＷ乗回路、１／４除算回路５６を１／Ｗ除算回路に変更すれば良い。例えばPAM変調であればＷ＝２、8PSK変調であればＷ＝８となる。
＜第６実施形態＞
以下、図面を参照して、本発明の第６実施形態について説明する。
図１２は、本発明の第６実施形態による受信装置２００の構成例を示すブロック図である。
この受信装置２００は、周波数オフセット精推定部４１、周波数オフセット粗推定部５７、掃引周波数範囲制御部５８、偏波分離部６６、周波数オフセット補償部６７、第1の位相補償部６８、第２の位相補償部６９、第1の判定部７０、および第２の判定部７１を備えている。ここで、図１２に示す、周波数オフセット精推定部４１、周波数オフセット粗推定部５７、および掃引周波数範囲制御部５８は、図９および図１１に示す同一符号の各部とそれぞれ対応している。
周波数オフセット補償部６７は、周波数オフセット推定装置１５０によって推定された受信信号の周波数オフセットの値に基づいて、受信信号の周波数オフセットを補償する。第１の位相補償部６８および第２の位相補償部６９は、周波数オフセット補償部６７によって補償された受信信号に対して、位相を補償する。第1の判定部７０および第２の判定部７１は、前記位相を補償された前記受信信号のシンボルの判定を行う。
本第６実施形態の受信装置２００は、まず周波数オフセット粗推定部５７で周波数オフセットを粗く推定した値に基づいて、数式１０を満足するように、周波数オフセット精推定部４１の掃引周波数を掃引周波数範囲制御部５８において決める。
偏波分離部６６は、入力信号I+jQおよび入力信号I’+jQ’に対して偏波分離および残留分散の除去を実行する。周波数オフセット精推定部４１は、偏波分離部６６から出力された信号を入力し、周波数オフセットの推定を実行する。
この周波数オフセット精推定部４１で推定した周波数オフセット推定値は、周波数オフセット補償部６７に入力され、入力信号I+jQおよび入力信号I’+jQ’の周波数オフセットを補償する。次に、周波数オフセット補償部６７において周波数オフセットが補償された信号に対して、第１の位相補償部６８および第２の位相補償部６９は位相を補償し、それぞれ第１の判定部７０および第２の判定部７１へ出力する。第１の判定部７０および第２の判定部７１は、シンボルの判定をして復調信号を生成する。
以上のように、周波数オフセットの推定と共に、偏波分離および残留分散の除去を合わせて実行することで、周波数オフセットおよび分散値が共に小さい入力信号に対して、周波数オフセット精推定部４１を動作させることができる。これにより、本第６実施形態の周波数オフセット推定装置１５０を備える受信装置２００は、受信信号に対して、広帯域かつ高精度に周波数オフセットを推定して補償し、復調することができる。
なお、周波数オフセット精推定部４１に入力される信号は、偏波分離部６６からの出力信号に代えて、偏波分離部６６に入力される前の信号を入力しても良い。ただし、偏波分離部６６において偏波分離および残留分散の除去をされた信号を周波数オフセット精推定部４１へ入力した方が、より高精度に周波数オフセットを推定することができる。
＜第７実施形態＞
以下、図面を参照して、本発明の第７実施形態について説明する。
図１３は、本発明の第７実施形態による周波数オフセット推定装置１５０の構成例を示すブロック図である。
この周波数オフセット推定装置１５０は、周波数オフセット粗推定部５７、周波数オフセット精推定部４１、および周波数不確定性除去制御部９７を備えている。周波数オフセット精推定部４１は、第１のＷ乗回路８８、第２のＷ乗回路８９、第１のＦＦＴ部９０、第２のＦＦＴ部９１、第１の絶対値２乗回路９２、第２の絶対値２乗回路９３、偏波スペクトル加算回路９４、Ｕフレームスペクトル加算回路９５、および最大値周波数検出回路９６を備えている。なお、Ｗは変調方式によって定まる数（正の整数）であり、周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗとする。QAM変調であればＷ＝４、PAM変調であればＷ＝２、M-PSK変調であればＷ＝Ｍである。
周波数オフセット粗推定部５７には、第１実施形態に記載の周波数スペクトルを利用する周波数オフセット推定装置１５を使うことができる。また、第１実施形態に記載の周波数オフセット推定装置１５に代えて、非特許文献３に記載の発明も同様に使うことができる。本実施形態は、時間平均である評価関数が離散フーリエ変換にもなっていることに着目した構成である。
本実施形態において周波数オフセット精推定部４１は、受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗと定義した場合に、受信信号における２つの偏波をそれぞれＷ乗した後に周波数スペクトルに変換する。そして、周波数オフセット精推定部４１は、その変換結果に対して絶対値または絶対値のべき乗の演算を行い、これら２つの偏波の周波数スペクトルを加算するとともに、Ｎ（Ｎ：正の整数）シンボルからなる各フレームの周波数スペクトルに対してＵ（Ｕ：正の整数）フレーム加算または平均化を行い、この演算結果が最大値となる周波数を検出する。
まず、周波数オフセット精推定部４１は、入力信号I+jQ およびI'+jQ'を第１のＷ乗回路８８および第２のＷ乗回路８９でＷ乗した後、第１のＦＦＴ部９０および第２のＦＦＴ部９１でＮ個のシンボル系列を一括して周波数スペクトルに変換する。次に、周波数オフセット精推定部４１は、この周波数スペクトルのデータに対して、第１の絶対値２乗回路９２および第２の絶対値２乗回路９３で絶対値の２乗を求め、偏波スペクトル加算回路９４において２つの偏波の周波数スペクトルを加算する。ここで、周波数オフセット精推定部４１は、Ｎ個のシンボル系列を1フレーム目、次のＮ個のシンボル系列を２フレーム目とみなし、Ｕフレームの周波数スペクトルをＵフレームスペクトル加算回路９５において加算する。当該加算処理により、周波数オフセット精推定部４１は、雑音成分を除去する。
この加算された周波数スペクトルから、最大値周波数検出回路９６は、最大値を取る周波数ｆを算出する。第４実施形態や第５実施形態では、ＮＣＯ４６の掃引周波数を所定のステップで変化させて、各々の周波数において評価関数φn(f)を求めた後にφn(f)が最大となる周波数ｆを検出した。これに対し本第７実施形態では、周波数スペクトルを求め、周波数領域で周波数スペクトルが最大となる周波数ｆを検出する。
また、第４実施形態や第５実施形態では、周波数不確定性のない周波数オフセットを推定するために、ＮＣＯ４６の掃引周波数の周波数範囲を予め掃引周波数範囲制御部５８で制限していた。一方、本第７実施形態で用いる周波数スペクトルには周波数不確定性が含まれている。この周波数不確定性を周波数不確定性除去制御部９７で除去する。周波数不確定性除去制御部９７は、周波数オフセット粗推定部５７によって推定された周波数オフセットの値に基づいて、周波数オフセット精推定部４１によって推定された周波数オフセットの周波数不確定性を除去し、周波数オフセットを推定する。
例えば、周波数不確定性除去制御部９７は、周波数不確定性を除去し周波数オフセットを推定する際に、周波数オフセット精推定部４１によって推定された周波数オフセットに基づいて、周波数不確定性を含む前記周波数オフセットの候補となる周波数を算出する。さらに、周波数不確定性除去制御部９７は、周波数軸上で隣り合う周波数オフセットの候補となる周波数の中点を、隣り合う前記周波数オフセットの候補の境界とし、周波数軸上で境界に基づいた領域の中から、周波数オフセット粗推定部５７によって推定された値が含まれる領域を選択する。そして、周波数不確定性除去制御部９７は、前記選択した領域に含まれる周波数オフセットの候補となる周波数を、周波数オフセットの推定値として選択する。
以下、周波数不確定性を除去して周波数オフセットを推定する本実施形態の動作の具体例について、図１４および数式１２〜数式１６を用いて説明する。図１４は、第７実施形態による周波数オフセット推定装置１５０が周波数オフセットを推定する動作の概要を示す説明図である。
Figure 2011158932
Figure 2011158932
Figure 2011158932
Figure 2011158932
Figure 2011158932
数式１２の絶対値記号の中は、yW(u,p,t)に関する離散フーリエ変換となっている。すなわちyW(u,p,t)に関する周波数スペクトルとなっている。この周波数スペクトルの最大値を求める周波数範囲は、数式１３によって、折り返しの影響を受けない第１ナイキストゾーンである[-RS/2からRS/2]に制限している。受信信号のＷ乗の周波数スペクトルが最大となる周波数をfMAX とすると、周波数オフセットの推定値がfMAXとなる可能性のある周波数オフセットの候補の周波数fCND(k)は数式１４で表される。図１４はＷ＝４の場合について示している。この図１４にはfMAX を縦軸に、またfCND(k)を横軸にプロットしている。
周波数オフセット粗推定部９７の粗推定値をfCO とし、隣り合うfCND(k)の中点を判定境界とする。この時、２つの判定境界に挟まれる領域のうち、fCO を含む領域が実際の周波数オフセットを与える。数式１５ではk=i の領域がfCO を含んでいることを表し、数式１６により周波数オフセットの推定値を確定することができる。この場合、図１４に示している例では、次の数式１７が成り立つ。
Figure 2011158932
以上のように、本第７実施形態の周波数オフセット推定装置１５０は、周波数オフセット粗推定部５７によって推定された周波数オフセットの値に基づいて、周波数オフセット精推定部４１によって推定された周波数オフセットの周波数不確定性を除去し、周波数オフセットを推定する。これにより、本第７実施形態の周波数オフセット推定装置１５０は、受信信号に対して広帯域かつ高精度に周波数オフセットを推定することができる。また、第４実施形態や第５実施形態と比較すると、掃引周波数範囲制御部５８およびＮＣＯ４６をＦＦＴ部９０および９１に置き換えて一括して周波数スペクトルに変換して推定するため、推定に要する時間を短縮することができる。
なお、周波数スペクトルに変換するＦＦＴ部９０およびＦＦＴ部９１は、周波数解析方法として高速フーリエ変換に限らず、短時間フーリエ変換、離散フーリエ変換、ウエーブレット変換、あるいはその他の周波数解析方法を用いても良い。また、第１の絶対値２乗回路９２および第２の絶対値２乗回路９３の２乗を１乗や他のべき乗に変えても周波数スペクトルの各周波数成分の相対的な大小関係は変わらず、動作に影響しない。また偏波スペクトル加算回路９４において偏波加算を行わず、ピーク値の大きい方の偏波の周波数スペクトルを選択することも可能である。
図１５は、本第７実施形態の実験結果の一例を示す説明図である。
この実験結果は、シンボルレートRSが10GHzの64QAMで変調された信号を受信した時のＵフレームスペクトル加算回路９５の出力信号の周波数スペクトルを示している。また周波数オフセット粗推定部５７の推定値fCO は3.7GHzであった。図１５より、数式１２のfMAX は-4GHzである。数式１４およびＷ＝４ よりfCND(k)は、-1GHz(k=0)、1.5GHz(k=1)、4GHz(k=2)、6.5GHz(k=3)、…となり、このうち数式１５を満たすｉは２である。従って、数式１６より周波数オフセットは4GHzと定めることができる。以上のように本実施形態により、RS/8=1.25GHzを越える周波数オフセットを推定することが可能となる。
＜第８実施形態＞
以下、図面を参照して、本発明の第８実施形態について説明する。
図１６は、本発明の第８実施形態による受信装置２００の構成例を示すブロック図である。なお、図１６において図１２の各部に対応する構成には同一の符号をつけており、その説明を省略する。
偏波分離部６６は、入力信号I+jQおよび入力信号I’+jQ’に対して偏波分離および残留分散の除去を実行する。周波数オフセット精推定部４１は、偏波分離部６６から出力された信号を入力し、周波数オフセットの推定を実行する。
周波数不確定性除去制御部９７は、周波数オフセット粗推定部５７の推定値および周波数オフセット精推定部４１の推定値がそれぞれ入力され、数式１２から数式１６に従って周波数オフセット推定値を確定する。そして、周波数オフセット補償部６７は、周波数不確定性除去制御部９７から出力される確定した周波数オフセットの推定値が入力され、入力信号の周波数オフセットを補償する。
次に、周波数オフセット補償部６７において周波数オフセットが補償された信号に対して、第１の位相補償部６８および第２の位相補償部６９は位相を補償し、それぞれ第１の判定部７０および第２の判定部７１へ出力する。第１の判定部７０および第２の判定部７１は、シンボルの判定をして復調信号を生成する。
以上のように、本第８実施形態の受信装置２００は、周波数オフセット推定装置１５０による周波数オフセットの推定と共に、偏波分離および残留分散の除去を合わせて実行することで、周波数オフセットおよび分散値が共に小さい入力信号に対して、周波数オフセット精推定部４１を動作させることができる。これにより、本第８実施形態の受信装置２００は、受信信号に対して、広帯域かつ高精度に周波数オフセットを推定して補償し、復調することができる。
なお、周波数オフセット精推定部４１に入力される信号は、偏波分離部６６からの出力信号に代えて、偏波分離部６６に入力される前の信号を入力しても良い。ただし、偏波分離部６６において偏波分離および残留分散の除去をされた信号を周波数オフセット精推定部４１へ入力した方が、より高精度に周波数オフセットを推定することができる。
また、図１、図６、図７、図９，図１１，図１２、図１３、および図１６における周波数オフセットの推定装置を構成する各部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各構成による処理が実行されてもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳ(operating system)や周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ＷＷＷ（world wide web）システムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（read only memory）、ＣＤ（compact disc）−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
例えば、周波数不確定性除去制御部を備えた受信装置または周波数オフセット推定装置において、周波数オフセットの精推定を実現するための構成は、位相増加周波数オフセット推定装置（部）２０（図７）や周波数オフセット精推定部４１（図１３，図１６）の構成に限定されるものではなく、これらに代えて、図５や図１４にのこぎり波状に示されているような周期的な周波数オフセット推定特性を持つ従来の一般的な構成を用いてもよい。
本発明は、例えば、ディジタルコヒーレント光受信機や無線通信受信機に利用することができる。本発明によれば、受信信号の周波数オフセットを広帯域，高速かつ高精度に推定することができる。
１・・・ＦＦＴ部、３・・・バンドパスフィルタ、１２・・・周波数オフセット推定制御部、１５、１５０・・・周波数オフセット推定装置、２０、６７・・・位相増加周波数オフセット推定装置、２４、２５、６８、６９・・・位相補償部、２６，２７、７０，７１・・・判定部、２８、２９・・・周波数オフセット補償部、３０、９７・・・周波数不確定性除去制御部、４１・・・周波数オフセット精推定部、４１ａ・・・第１の演算部、４１ｂ・・・第２の演算部、４１ｃ・・・第３の演算部、４１ｄ・・・第４の演算部、５７・・・周波数オフセット粗推定部、５８・・・掃引周波数範囲制御部、２００・・・受信装置

Claims (20)

  1. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する周波数オフセット推定装置であって、
    所定のサンプリング周波数で予めサンプリングされた前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換部と、
    前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限部と、
    前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御部と、
    を備える周波数オフセット推定装置。
  2. 前記周波数オフセット推定制御部は、前記周波数スペクトルの前記全ての周波数成分を前記周波数軸上で循環して移動させる場合に、
    前記正の周波数成分の電力が前記負の周波数成分の電力より大きい場合は、前記周波数スペクトルの全周波数成分を負の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号が１未満の場合は周波数成分番号にＮを加算し、
    前記正の周波数成分の電力が前記負の周波数成分の電力以下の場合は、前記周波数スペクトルの全周波数成分を正の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号がＮを越える場合は周波数成分番号からＮを減算する請求項１に記載の周波数オフセット推定装置。
  3. 請求項１あるいは請求項２に記載の周波数オフセット推定装置と、
    前記周波数オフセット推定装置によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する第１の周波数オフセット補償部と、
    前記第１の周波数オフセット補償部によって補償された前記受信信号に対して、前記周波数オフセットを位相増加アルゴリズムに基づいて推定する、位相増加周波数オフセット推定部と、
    前記位相増加周波数オフセット推定部によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記周波数オフセットを補償する第２の周波数オフセット補償部と、
    を備える受信装置。
  4. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差である周波数オフセットを推定する周波数オフセット推定装置において用いられる、周波数オフセット推定方法であって、
    所定のサンプリング周波数で予めサンプリングされた前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換手順と、
    前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限手順と、
    前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御手順と、
    を備える周波数オフセット推定方法。
  5. 前記周波数オフセット推定制御手順において、前記周波数スペクトルの前記全ての周波数成分を前記周波数軸上で循環して移動させる場合に、
    前記正の周波数成分の電力が前記負の周波数成分の電力より大きい場合は、前記周波数スペクトルの全周波数成分を負の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号が１未満の場合は周波数成分番号にＮを加算し、
    前記正の周波数成分の電力が前記負の周波数成分の電力以下の場合は、前記周波数スペクトルの全周波数成分を正の方向に予め定められた大きさだけ移動し、移動後の周波数成分番号がＮを越える場合は周波数成分番号からＮを減算することにより、前記周波数スペクトルを循環して移動する請求項４に記載の周波数オフセット推定方法。
  6. 請求項４あるいは請求項５に記載の周波数オフセット推定方法による手順と、
    前記周波数オフセット推定方法による手順によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する第１の周波数オフセット補償手順と、
    前記第１の周波数オフセット補償手順によって補償された前記受信信号に対して、前記周波数オフセットを位相増加アルゴリズムに基づいて推定する、位相増加周波数オフセット推定手順と、
    前記位相増加周波数オフセット推定手順によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記周波数オフセットを補償する第２の周波数オフセット補償手順と、
    を備える受信方法。
  7. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、
    所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定部と、
    前記周波数オフセット粗推定部の粗推定値に基づいて掃引周波数の範囲を決定する掃引周波数範囲制御部と、
    前記掃引周波数範囲制御部により決定された前記掃引周波数の範囲において前記受信信号の周波数オフセットを推定する周波数オフセット精推定部と、
    を備え、
    前記周波数オフセット精推定部は、
    前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後、該受信信号に対して該受信信号の周波数から前記掃引周波数を減算する周波数変換を行う第１の演算部と、
    前記第１の演算部における各偏波の演算結果に対して、Ｎ（Ｎ：正の整数）シンボルの平均化または加算を行った後、絶対値または絶対値のべき乗の演算を行う第２の演算部と、
    前記第２の演算部の後段側において、前記２つの偏波の演算結果を加算、またはピーク値の大きい方の偏波の演算結果を選択するとともに、前記ＮシンボルからなるフレームについてＵ（Ｕ：正の整数）フレームの加算または平均化を行う第３の演算部と、
    前記第３の演算部の演算結果が最大値となる掃引周波数を検出し、当該掃引周波数を１／Ｗ倍して前記周波数オフセットを推定する第４の演算部と、
    を有する周波数オフセット推定装置。
  8. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、
    所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定手順と、
    前記周波数オフセット粗推定手順により推定された粗推定値に基づいて掃引周波数の範囲を決定する掃引周波数範囲制御手順と、
    前記掃引周波数範囲制御手順により決定された前記掃引周波数の範囲において前記受信信号の周波数オフセットを推定する周波数オフセット精推定手順と、
    を備え、
    前記周波数オフセット精推定手順は、
    前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後、該受信信号に対して該受信信号の周波数から前記掃引周波数を減算する周波数変換を行う第１の演算手順と、
    前記第１の演算手順における各偏波の演算結果に対して、Ｎ（Ｎ：正の整数）シンボルの加算または平均化を行った後、絶対値または絶対値のべき乗の演算を行う第２の演算手順と、
    前記第２の演算手順の後において、前記２つの偏波の演算結果を加算、またはピーク値の大きい方の偏波の演算結果を選択するとともに、前記ＮシンボルからなるフレームについてＵ（Ｕ：正の整数）フレームの加算または平均化を行う第３の演算手順と、
    前記第３の演算手順の演算結果が最大値となる掃引周波数を検出し、当該掃引周波数を１／Ｗ倍して前記周波数オフセットを推定する第４の演算手順と、
    を含む周波数オフセット推定方法。
  9. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、
    所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定部と、
    前記受信信号または前記受信信号の分散が補償された信号に対して周期的な周波数オフセット推定特性を有する周波数オフセット精推定部と、
    前記周波数オフセット粗推定部によって推定された前記周波数オフセットの値に基づいて、前記周波数オフセット精推定部によって推定された周波数オフセットの周波数不確定性を除去し、前記周波数オフセットを推定する周波数不確定性除去制御部と
    を備える周波数オフセット推定装置。
  10. 前記周波数オフセット精推定部は、前記受信信号または前記受信信号の分散が補償された信号に対して位相増加アルゴリズムに基づいて前記周波数オフセットを推定する請求項９に記載の周波数オフセット推定装置。
  11. 前記周波数オフセット精推定部は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後に周波数スペクトルに変換し、その変換結果に対して絶対値または絶対値のべき乗の演算を行い、これら２つの偏波の周波数スペクトルを加算、またはピーク値の大きい方の偏波の周波数スペクトルを選択するとともに、Ｎ（Ｎ：正の整数）シンボルからなるフレームの周波数スペクトルについてＵ（Ｕ：正の整数）フレームの加算または平均化を行い、この演算結果が最大値となる周波数を検出する請求項９に記載の周波数オフセット推定装置。
  12. 前記周波数不確定性除去制御部は、前記周波数不確定性を除去し前記周波数オフセットを推定する際に、
    前記周波数オフセット精推定部によって推定された前記周波数オフセットに基づいて、前記周波数不確定性を含む前記周波数オフセットの候補となる周波数を算出し、
    周波数軸上で隣り合う前記周波数オフセットの候補となる周波数の中点を、隣り合う前記周波数オフセットの候補の境界とし、
    前記周波数軸上で前記境界に基づいた領域の中から、前記周波数オフセット粗推定部によって推定された前記値が含まれる領域を選択し、
    前記選択した領域に含まれる前記周波数オフセットの候補となる周波数を、前記周波数オフセットの推定値として選択する請求項９〜１１のいずれか１項に記載の周波数オフセット推定装置。
  13. 前記周波数オフセット粗推定部は、
    前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換部と、
    前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限部と、
    前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御部と、
    を備える請求項９〜１２のいずれか１項に記載の周波数オフセット推定装置。
  14. 請求項７、請求項９、請求項１０、請求項１１、請求項１２、または請求項１３に記載の周波数オフセット推定装置と、
    前記周波数オフセット推定装置によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する周波数オフセット補償部と、
    前記周波数オフセット補償部によって補償された前記受信信号に対して、位相を補償する位相補償部と、
    前記位相を補償された前記受信信号のシンボルの判定を行う判定部と、
    を備える受信装置。
  15. 受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、
    所定のサンプリング周波数で予めサンプリングされた２つの偏波から成る受信信号の周波数スペクトルから周波数オフセットを推定する周波数オフセット粗推定手順と、
    前記受信信号または前記受信信号の分散が補償された信号に対して周期的な周波数オフセット推定特性を有する周波数オフセット精推定手順と、
    前記周波数オフセット粗推定手順によって推定された前記周波数オフセットの値に基づいて、前記周波数オフセット精推定手順によって推定された周波数オフセットの周波数不確定性を除去し、前記周波数オフセットを推定する周波数不確定性除去制御手順と、
    を備える周波数オフセット推定方法。
  16. 前記周波数オフセット精推定手順は、前記受信信号または前記受信信号の分散が補償された信号に対して位相増加アルゴリズムに基づいて前記周波数オフセットを推定する請求項１５に記載の周波数オフセット推定方法。
  17. 前記周波数オフセット精推定手順は、前記受信信号の周波数オフセットの無いコンスタレーション上の信号点の持つ回転対称性を２π／Ｗ（Ｗ：正の整数）と定義した場合に、前記受信信号における前記２つの偏波をそれぞれＷ乗した後に周波数スペクトルに変換し、その変換結果に対して絶対値または絶対値のべき乗の演算を行い、これら２つの偏波の周波数スペクトルを加算、またはピーク値の大きい方の偏波の周波数スペクトルを選択するとともに、Ｎ（Ｎ：正の整数）シンボルからなるフレームの周波数スペクトルについてＵ（Ｕ：正の整数）フレームの加算または平均化を行い、この演算結果が最大値となる周波数を検出する請求項１５に記載の周波数オフセット推定方法。
  18. 前記周波数不確定性除去制御手順は、前記周波数不確定性を除去し前記周波数オフセットを推定する際に、
    前記周波数オフセット精推定手順によって推定された前記周波数オフセットに基づいて、前記周波数不確定性を含む前記周波数オフセットの候補となる周波数を算出し、
    周波数軸上で隣り合う前記周波数オフセットの候補となる周波数の中点を、隣り合う前記周波数オフセットの候補の境界とし、
    前記周波数軸上で前記境界に基づいた領域の中から、前記周波数オフセット粗推定手順によって推定された前記値が含まれる領域を選択し、
    前記選択した領域に含まれる前記周波数オフセットの候補となる周波数を、前記周波数オフセットの推定値として選択する請求項１５〜１７のいずれか１項に記載の周波数オフセット推定方法。
  19. 前記周波数オフセット粗推定手順は、
    前記受信信号を周波数変換し、周波数の大きさの順に１からＮ（Nは任意の自然数）まで順序付けしたＮ個の周波数成分を持つ周波数スペクトルを出力する周波数変換手順と、
    前記周波数スペクトルの、１からＮ/２までの周波数成分番号を持った周波数成分である負の周波数成分とＮ/２＋１からＮまでの周波数成分番号を持った周波数成分である正の周波数成分とを、それぞれ周波数帯域制限する周波数帯域制限手順と、
    前記周波数帯域制限された前記周波数スペクトルの前記正の周波数成分と前記負の周波数成分とをそれぞれ２乗加算してそれぞれの電力を算出し、前記正の周波数成分の電力と前記負の周波数成分の電力とから算出した電力差の絶対値が、予め設定されたしきい値以下になるまで前記周波数スペクトルの全ての周波数成分を周波数軸上で循環して移動させ、前記しきい値以下になるまで移動させた移動量に基づいて前記周波数オフセットを推定する周波数オフセット推定制御手順と、
    を備える請求項１５〜１８のいずれか１項に記載の周波数オフセット推定方法。
  20. 請求項８、請求項１５、請求項１６、請求項１７、請求項１８、または請求項１９に記載の周波数オフセット推定方法による手順と、
    前記周波数オフセット推定方法による手順によって推定された前記受信信号の前記周波数オフセットの値に基づいて、前記受信信号の前記周波数オフセットを補償する周波数オフセット補償手順と、
    前記周波数オフセット補償手順によって補償された前記受信信号に対して、位相を補償する位相補償手順と、
    前記位相を補償された前記受信信号のシンボルの判定を行う判定手順と、
    を備える受信方法。
JP2012520505A 2010-06-17 2011-06-17 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法 Active JP5404926B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012520505A JP5404926B2 (ja) 2010-06-17 2011-06-17 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2010138402 2010-06-17
JP2010138402 2010-06-17
JP2010251868 2010-11-10
JP2010251868 2010-11-10
JP2012520505A JP5404926B2 (ja) 2010-06-17 2011-06-17 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法
PCT/JP2011/063903 WO2011158932A1 (ja) 2010-06-17 2011-06-17 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法

Publications (2)

Publication Number Publication Date
JPWO2011158932A1 true JPWO2011158932A1 (ja) 2013-08-19
JP5404926B2 JP5404926B2 (ja) 2014-02-05

Family

ID=45348324

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012520505A Active JP5404926B2 (ja) 2010-06-17 2011-06-17 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法

Country Status (5)

Country Link
US (1) US8781029B2 (ja)
EP (2) EP2903173B1 (ja)
JP (1) JP5404926B2 (ja)
CN (2) CN102934378B (ja)
WO (1) WO2011158932A1 (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5545892B2 (ja) 2009-05-18 2014-07-09 日本電信電話株式会社 信号生成回路および信号受信回路、信号生成回路、光信号送信装置、信号受信回路、光信号同期確立方法、および光信号同期システム
US9264278B2 (en) * 2012-10-19 2016-02-16 Apple Inc. Robust scalable and adaptive frequency estimation and frequency tracking for wireless systems
JP6082125B2 (ja) * 2012-12-14 2017-02-15 テレフオンアクチーボラゲット エルエム エリクソン（パブル） マルチキャリア変調信号のための受信機
WO2014097583A1 (ja) * 2012-12-20 2014-06-26 日本電気株式会社 周波数オフセット補償装置および周波数オフセット補償方法
CN104104624B (zh) * 2013-04-10 2018-02-16 富士通株式会社 频差估计方法、装置和系统
ES2644734T3 (es) * 2013-07-31 2017-11-30 Huawei Technologies Co., Ltd. Receptor OFDM multi-longitud de onda, método de compensación del desplazamiento de frecuencia y sistema
JP6458730B2 (ja) * 2013-08-21 2019-01-30 日本電気株式会社 周波数偏差補償方式、周波数偏差補償方法及びプログラム
EP2849362A1 (en) * 2013-09-11 2015-03-18 Alcatel Lucent Blind frequency offset estimation for pulse-shaped signals
WO2015052874A1 (ja) * 2013-10-09 2015-04-16 日本電信電話株式会社 光伝送システム
US9866364B2 (en) * 2014-08-29 2018-01-09 Huawei Technologies Co., Ltd. System and method for semi-orthogonal multiple access
HUE048813T2 (hu) 2014-09-24 2020-08-28 Guangdong Oppo Mobile Telecommunications Corp Ltd Eljárás és vezeték nélküli kommunikációs eszköz vett jel frekvenciaeltolásának becslésére
US9497048B2 (en) * 2015-01-23 2016-11-15 Huawei Technologies Co., Ltd. Method and apparatus with improved estimation accuracy and robustness for fine frequency offset estimation in coherent receivers
WO2016161638A1 (zh) * 2015-04-10 2016-10-13 华为技术有限公司 一种相干光源频偏估计和补偿的相干接收机、方法和系统
US9819405B2 (en) * 2015-06-12 2017-11-14 Huawei Technologies Co., Ltd. Systems and methods for transmission pairing mixed transmission modes
CN105141564B (zh) * 2015-07-28 2018-10-12 广东顺德中山大学卡内基梅隆大学国际联合研究院 一种高子载波数高阶调制水平ofdm采样频率同步方法
EP3360476B1 (en) * 2017-02-13 2020-09-30 Stichting IMEC Nederland A method and a device for detecting of a vital sign
CN111406385B (zh) 2017-11-24 2021-12-28 华为技术有限公司 一种用于网络接入节点生成相位补偿后的调制符号的处理设备
JP7064141B2 (ja) * 2018-09-05 2022-05-10 日本電信電話株式会社 光受信装置、及び周波数オフセット推定方法
CN110995632B (zh) * 2019-11-29 2023-03-21 深圳市统先科技股份有限公司 卫星通信带宽复用电路
US11251766B2 (en) * 2020-01-13 2022-02-15 Maxim Integrated Products, Inc. Ultra-wide band frequency offset estimation systems and methods for analog coherent receivers
CN114422314B (zh) * 2021-12-22 2024-05-24 江苏科大亨芯半导体技术有限公司 多载波鉴频方法及系统
CN115333603B (zh) * 2022-07-14 2023-06-27 航天恒星科技有限公司 载波同步方法、装置、电子设备及存储介质
CN115333654B (zh) * 2022-10-13 2023-02-03 成都爱旗科技有限公司 一种频偏检测方法、系统及电子设备

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5307515A (en) 1991-08-05 1994-04-26 Ford Motor Company Adjacent channel controller for radio receiver
EP1482699A3 (en) * 1998-06-08 2005-02-09 Telefonaktiebolaget LM Ericsson (publ) Burst structure for multicarrier transmission, and synchronisation of bursts, symbols and frequency
US7221917B2 (en) 2002-05-01 2007-05-22 Ibiquity Digital Corporation Adjacent channel interference mitigation for FM digital audio broadcasting receivers
US8428198B2 (en) * 2006-03-15 2013-04-23 Qualcomm Incorporated Frequency tracking which adapts to timing synchronization
JP4983365B2 (ja) 2006-05-16 2012-07-25 ソニー株式会社 無線通信装置
CN102611673B (zh) * 2007-01-05 2015-01-21 Lg电子株式会社 在考虑了频率偏移的情况下设定循环移位的方法
CN101442364B (zh) * 2007-11-19 2011-10-19 富士通株式会社 光相干接收机、光相干接收机用频差估计装置及方法
EP2187531A4 (en) 2008-01-21 2016-09-07 Nec Corp COMMUNICATION APPARATUS, COMMUNICATION SYSTEM, CONTROL METHOD, AND CONTROL PROGRAM
US8649750B2 (en) * 2008-03-11 2014-02-11 Intel Corporation Techniques for efficient carrier recovery for passband communciation systems
JP5545892B2 (ja) * 2009-05-18 2014-07-09 日本電信電話株式会社 信号生成回路および信号受信回路、信号生成回路、光信号送信装置、信号受信回路、光信号同期確立方法、および光信号同期システム

Also Published As

Publication number Publication date
WO2011158932A1 (ja) 2011-12-22
US20130070874A1 (en) 2013-03-21
EP2584720A4 (en) 2013-11-27
EP2903173A1 (en) 2015-08-05
EP2584720B1 (en) 2015-03-11
CN104539342A (zh) 2015-04-22
EP2584720A1 (en) 2013-04-24
CN102934378A (zh) 2013-02-13
CN102934378B (zh) 2015-01-14
US8781029B2 (en) 2014-07-15
CN104539342B (zh) 2018-01-09
JP5404926B2 (ja) 2014-02-05
EP2903173B1 (en) 2019-09-04

Similar Documents

Publication Publication Date Title
JP5404926B2 (ja) 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法
JP5304083B2 (ja) 周波数オフセットモニタリング装置及びコヒーレント光受信機
JP5391660B2 (ja) 光コヒーレント受信機、周波数オフセット推定装置、及び光コヒーレント受信機の方法
JP5651990B2 (ja) デジタルコヒーレント受信器および受信方法
US8908809B2 (en) Complexity reduced feed forward carrier recovery methods for M-QAM modulation formats
US7203254B2 (en) Method and system for synchronizing in a frequency shift keying receiver
EP2583424B1 (en) Method for phase and oscillator frequency estimation
US8463121B2 (en) Ultra wide-range frequency offset estimation for digital coherent optical receivers
US8514922B2 (en) Filter coefficient control apparatus and method
KR0157500B1 (ko) 자동주파수조절방법 및 그 장치
US9042741B2 (en) System and method for blind frequency recovery
US7430247B2 (en) Carrier frequency detection for N-ary phase modulated signal
US8989316B1 (en) Method for estimating and correcting a carrier frequency offset over dispersive but spectrally flat channels
JP2020010195A (ja) 周波数推定装置および追尾受信機
Zhou et al. DFT-based carrier recovery for Satellite DVB Receivers
JP5656750B2 (ja) シングルキャリア受信装置
JP6809814B2 (ja) 信号検出装置及び信号検出方法
JP2014017609A (ja) 受信装置および受信方法
KR20040038247A (ko) 디지털 ｔｖ의 락 검출 장치
JP2012147094A (ja) シンボル推定回路及び復調回路

Legal Events

Date Code Title Description
RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20130606

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20130809

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: 20131022

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131029

R150 Certificate of patent or registration of utility model

Ref document number: 5404926

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150