JP6876099B2 - 適応等化器、適応等化方法及び光通信システム - Google Patents

適応等化器、適応等化方法及び光通信システム Download PDF

Info

Publication number
JP6876099B2
JP6876099B2 JP2019115330A JP2019115330A JP6876099B2 JP 6876099 B2 JP6876099 B2 JP 6876099B2 JP 2019115330 A JP2019115330 A JP 2019115330A JP 2019115330 A JP2019115330 A JP 2019115330A JP 6876099 B2 JP6876099 B2 JP 6876099B2
Authority
JP
Japan
Prior art keywords
tap coefficient
polarization state
step size
polarization
received signal
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.)
Active
Application number
JP2019115330A
Other languages
English (en)
Other versions
JP2021002743A (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.)
NTT Electronics Corp
Original Assignee
NTT Electronics 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 NTT Electronics Corp filed Critical NTT Electronics Corp
Priority to JP2019115330A priority Critical patent/JP6876099B2/ja
Priority to PCT/JP2020/022780 priority patent/WO2020255809A1/ja
Priority to US17/596,617 priority patent/US12040834B2/en
Priority to CN202080044772.9A priority patent/CN114026800B/zh
Publication of JP2021002743A publication Critical patent/JP2021002743A/ja
Application granted granted Critical
Publication of JP6876099B2 publication Critical patent/JP6876099B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25Arrangements specific to fibre transmission
    • H04B10/2507Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • H04B10/2569Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to polarisation mode dispersion [PMD]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/04Control of transmission; Equalising
    • H04B3/06Control of transmission; Equalising by the transmitted signal
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25Arrangements specific to fibre transmission
    • H04B10/2507Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/614Coherent receivers comprising one or more polarization beam splitters, e.g. polarization multiplexed [PolMux] X-PSK coherent receivers, polarization diversity heterodyne coherent receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • H04B10/6162Compensation of polarization related effects, e.g., PMD, PDL
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • H04B10/6166Polarisation demultiplexing, tracking or alignment of orthogonal polarisation components
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/04Control of transmission; Equalising

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Optical Communication System (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Description

本発明は、データ通信において伝送路の特性を補償する適応等化器、適応等化方法及び通信装置に関する。
コヒーレント光通信では、受信側において処理タイミングとの同期誤差及び伝送信号の歪をデジタル信号処理により補償することで、数十Gbit/s以上の大容量伝送を実現している。デジタル信号処理では、主に、波長分散補償、周波数制御・位相調整、偏波多重分離及び偏波分散補償等の処理を行っている。
特に、偏波状態SOP(State of Polarization)は、高速で変動を生じる場合がある。そのSOP変動に対するデジタル信号処理の追従性も重要な課題となっている。rad/μsオーダの速度でSOP変動が生じる場合、デジタル信号処理の追従性が低下し、符号誤り率BER(Bit Error Rate)が急激に悪化することが知られている。
偏波多重分離及び偏波分散補償の処理は、主に適応等化によって行われる。デジタル信号処理における適応等化器は、一般的にデジタルフィルタで構成される。そのデジタルフィルタに伝送信号の歪が相殺されるタップ係数を設定することで、伝送信号を補償できる。タップ係数は時間的に変化する状況に適応して逐次更新され、SOP変動に追従した補償が行われる。
適応等化器を構成するデジタルフィルタのタップ係数更新には、一般的には定包絡線基準アルゴリズム(CMA:Constant Modulus Algorithm)のような逐次更新アルゴリズムが用いられる。このタップ係数の適応制御時において、デジタル信号処理の追従性を決定するのがステップサイズμである。ステップサイズμが大きくなると、デジタル信号処理の追従性が向上し、高速なSOP変動に対する受信耐力が向上する。その一方、低速なSOP変動時には、雑音に対する影響を受けてBERが悪化する。
そこで、特許文献1では、ステップサイズμをSOP変動速度に応じて適応制御することが提案されている。
特開2013−223128号公報 特開2018−174413号公報 特開2013−168983号公報
しかしながら、特許文献1で提案されている方法では、低速なSOP変動により小さなステップサイズμが一旦設定されると、高速なSOP変動への追従性が低下するという問題があった。
本発明は、このような課題を解決するためになされたものであり、その目的は、様々なSOP変動に対して常に安定な追従性を実現する適応等化器、適応等化方法及び光通信システムを提供することにある。
このような目的を達成するために、本発明に係る適応等化器(70)は、タップ係数を変更することで受信信号(61)の偏波変動を適応的に補償する適応等化フィルタ(71)と、可変のステップサイズを用いて受信信号(61)の偏波変動に応じたタップ係数を計算し、計算されたタップ係数で適応等化フィルタ(71)のタップ係数を自身で更新する第1のタップ係数更新器(72)と、前記受信信号を入力とし前記適応等化フィルタの出力を受けることなく自身で、固定のステップサイズを用いて受信信号(61)の偏波変動に応じたタップ係数を計算する第2のタップ係数更新器(73)と、第2のタップ係数更新器(73)で計算されたタップ係数を用いて受信信号(61)の偏波状態を推定する偏波状態推定器(74)と、偏波状態推定器(74)で推定された偏波状態に対応したステップサイズを求め、求められたステップサイズで可変のステップサイズを更新するステップサイズ更新器(75)とを備えることを特徴とする。
第1のタップ係数更新器(72)及び第2のタップ係数更新器(73)の各々は、伝送されてきた既知信号と既知信号の真値との比較誤差を最小化するようにタップ係数をステップサイズ毎に更新して求める逐次更新アルゴリズムを用いてタップ係数を計算するように構成されてもよい。
ステップサイズ更新器(75)は、偏波状態変動量及び偏波状態変動速度の少なくとも一方の所定範囲毎に最適な誤り率特性が得られるステップサイズが予め記憶されているメモリ(97)を含んでいてもよい。
偏波状態推定器(74)は、受信信号(61)の偏波状態の推定値の移動平均を求め、求められた平均値をステップサイズ更新器(75)に供給する移動平均回路(98)を含んでいてもよい。偏波状態推定器(74)はまた、推定された偏波状態変動量からヒステリシス特性を利用してピークを検出するピーク検出回路(99)を含んでいてもよい。
また、本発明に係る適応等化方法は、タップ係数を変更することで受信信号(61)の偏波変動を適応的に補償する適応等化フィルタリングを実行するステップと、可変のステップサイズを用いて受信信号(61)の偏波変動に応じたタップ係数を計算し、計算されたタップ係数で適応等化フィルタリングのタップ係数を自身で更新する第1のタップ係数更新ステップと、前記受信信号を入力とし前記適応等化フィルタの出力を受けることなく自身で、固定のステップサイズを用いて受信信号(61)の偏波変動に応じたタップ係数を計算する第2のタップ係数更新ステップと、第2のタップ係数更新ステップで計算されたタップ係数を用いて受信信号(61)の偏波状態を推定する偏波状態推定ステップと、偏波状態推定ステップで推定された偏波状態に対応したステップサイズを求め、求められたステップサイズで可変のステップサイズを更新するステップサイズ更新ステップとを備えることを特徴とする。
第1のタップ係数更新ステップ及び第2のタップ係数更新ステップの各々は、伝送されてきた既知信号と既知信号の真値との比較誤差を最小化するようにタップ係数をステップサイズ毎に更新して求める逐次更新アルゴリズムを用いてタップ係数を計算するステップを含んでいてもよい。
ステップサイズ更新ステップは、偏波状態変動量及び偏波状態変動速度の少なくとも一方の所定範囲毎に最適な誤り率特性が得られるステップサイズが予め記憶されているメモリ(97)を使用して、ステップサイズを選択するステップを含んでいてもよい。
偏波状態推定ステップは、受信信号(61)の偏波状態の推定値の移動平均を求め、求められた平均値をステップサイズ更新ステップに供給するステップを含んでいてもよい。偏波状態推定ステップはまた、推定された偏波状態変動量からヒステリシス特性を利用してピークを検出するステップを含んでいてもよい。
また、本発明に係る光通信システムは、受信信号の波長分散による歪を補償する波長分散補償器(60)と、波長分散補償器(60)の出力信号の偏波変動を補償する上述した適応等化器(70)と、適応等化器(70)の出力信号を復号する復号器とを備えることを特徴とする。
光通信システムは、光信号を受信し、受信された光信号を上記受信信号に変換する光受信機(40)を更に備えていてもよい。
本発明では、適応等化フィルタ(71)のタップ係数を計算するための可変のステップサイズとは別に、固定のステップサイズを利用する。この固定のステップサイズを用いてタップ係数を計算し、このタップ係数を用いて偏波状態を推定する。このため、偏波状態が変動して可変のステップサイズがどのように変化しても、固定のステップサイズに基づいて偏波状態の高速に計算することができる。従って、様々な偏波状態の変動に対して常に安定な追従性を実現することができる。
更に、本発明では、上述したように推定された偏波状態に対応した可変のステップサイズを用いてタップ係数を計算し、そのタップ係数で適応等化フィルタ(71)のタップ係数を更新する。このため、偏波状態の変動が小さいときにも、BERの悪化を防止することができる。
図1は、本発明の実施の形態に係る光通信システムにおける送信側の構成を示すブロック図である。 図2は、本発明の実施の形態に係る光通信システムにおける受信側の構成を示すブロック図である。 図3は、本発明の実施の形態に係る適応等化器の構成を示すブロック図である。 図4は、図3における適応等化フィルタの構成を示すブロック図である。 図5は、図4における偏波状態推定器の出力のイメージを示すグラフである。 図6は、図3におけるステップサイズ更新器の構成を示すブロック図である。 図7は、図3における偏波状態推定器の構成を示すブロック図である。 図8は、適応等化器の比較例のブロック図である。
以下、図面を参照しながら、本発明の実施の形態について詳細に説明する。
[光通信システム]
図1に示すように、本発明の実施の形態に係る光通信システムは、送信側に、送信信号処理器10、及び光送信機20を備えている。
送信信号処理器10は、入力データに対して所定の処理を施す回路である。具体的には、送信信号処理器10は、入力データを水平偏波用データと垂直偏波用データに分け、それぞれのデータに対して誤り訂正用符号化、帯域制限フィルタリング、及び変調用マッピング等の処理を行う。このような処理を施された水平偏波用データ11a及び垂直偏波用データ11bが、光送信機20に出力される。
光送信機20は、水平偏波用データ11a及び垂直偏波用データ11bを光信号に変換し、変換された光信号を送信する回路である。この光送信機20は、信号光源(信号LD)21、2つの90゜合成器22a,22b、及び偏波合成器23を含んでいる。90゜合成器22a,22bは、信号光源21の出力光を水平偏波用データ11a及び垂直偏波用データ11bで変調することによって、これらのデータ11a,11bを光信号に変換する。偏波合成器23は、光信号に変換された水平偏波用データ11aと垂直偏波用データ11bを合成する。合成された光信号は、光ファイバ伝送路30を通して受信側に伝送される。
図2に示すように、本発明の実施の形態に係る光通信システムは、受信側に、光受信機40、AD変換器50、波長分散補償器60、適応等化器70、及び復号器80を備えている。
光受信機40は、光信号を受信し、受信された光信号を電気信号に変換し、変換された電気信号を出力する回路である。この光受信機40は、偏波分離器41、局部発振光源(局発LD)42、2つの90゜ハイブリッド回路43a,43b、及び光電変換器44を含んでいる。偏波分離器41は、光信号を2つの直交偏波成分、すなわちX偏波(水平偏波)とY偏波(垂直偏波)に分離する。90゜ハイブリッド回路43a,43bは、偏波分離器41から出力された光信号の各偏波に局部発振光源42の出力光を合成させ、光信号の各偏波を更に同相(I)成分と直交(Q)成分に分離する。光電変換器44は、90゜ハイブリッド回路43a,43bから出力された光信号の各成分を電気信号に変換し、この電気信号をX偏波信号45a及びY偏波信号45bとして出力する。以下、X偏波信号45a及びY偏波信号45bを受信信号と呼ぶ。
AD変換器50は、光受信機40から出力された受信信号をデジタル信号に変換する。
光ファイバ伝送路30を光信号が伝搬する際に、波長分散によって信号波形が歪む。波長分散補償器60は、その歪の大きさをAD変換器50から出力された受信信号から推定し、波長分散による受信信号の歪みを補償する。
また、送信側においてX偏波信号とY偏波信号を合成して送信し、受信側においてX偏波信号とY偏波信号を分離する際に、偏波モード分散によって偏波変動が生じ、信号波形が歪む。適応等化器70は、波長分散補償器60から出力された受信信号の偏波変動による歪みを補償する等化処理を行う。なお、偏波分離は最初に光受信機40で行われるが、適応等化器70でより完全に偏波分離が処理される。
復号器80は、適応等化器70から出力された受信信号を復号し、元のデータ(すなわち、送信信号処理器10の入力データ)を再生する。
[適応等化器の構成]
次に、本発明の実施の形態に係る適応等化器70の構成について説明する。図3に示すように、適応等化器70は、適応等化フィルタ71、第1のタップ係数更新器72、第2のタップ係数更新器73、偏波状態推定器74、及びステップサイズ更新器75を備えている。
適応等化フィルタ71は、タップ係数を更新することで受信信号61(Xin(n)、Yin(n))の偏波変動による波形歪みを適応的に補償する回路である。この適応等化フィルタ71は、一般的なデジタルフィルタ、特にはFIRフィルタ(Finite Impulse Response Filter)で構成できる。
図4に示す適応等化フィルタ71は、バタフライ型の回路構成を有する。すなわち、入力Xin(n)がFIRフィルタ91及びFIRフィルタ93に与えられる。また、入力Yin(n)がFIRフィルタ92及びFIRフィルタ94に与えられる。FIRフィルタ91の出力とFIRフィルタ92の出力が加算回路95で加算され、その加算結果が出力Xout(n)となる。また、FIRフィルタ93の出力とFIRフィルタ94の出力が加算回路96で加算され、その加算結果が出力Yout(n)となる。入力Xin(n)及び入力Yin(n)は、光受信機40で偏波分離され、波長分散補償器50で波長分散補償されたX偏波信号及びY偏波信号である。FIRフィルタ91、92、93、94に設定されるタップ係数をそれぞれhxx(n)、hxy(n)、hyx(n)、hyy(n)とすると、出力Xout(n)及び出力Yout(n)は下記の式で表される。但し、「*」は畳み込み演算を示す。nはデータのサンプル順を示す。なお、各タップ係数は、複数のタップ係数のグループである。
out(n)=hxx(n)*Xin(n)+hxy(n)*Yin(n)
out(n)=hyx(n)*Xin(n)+hyy(n)*Yin(n)
この適応等化フィルタ71によって、偏波分離や波形歪の補償が実現される。偏波の状態が変動する場合、補償精度を向上するため、適応等化フィルタ71のタップ係数も偏波変動に応じて更新される必要がある。
第1のタップ係数更新器72は、可変のステップサイズμを用いて受信信号61(Xin(n)、Yin(n))の偏波変動に応じたタップ係数を計算し、計算されたタップ係数で適応等化フィルタ71のタップ係数を更新する回路である。適応等化フィルタ71のタップ係数を計算するアルゴリズムとして例えばCMAが用いられ、下式によりタップ係数が更新される。但し、εx(n)及びεy(n)は誤差評価関数である。Xin *(n)及びYin *(n)は、FIRフィルタ91〜94の入力Xin(n)及びYin(n)の複素共役である。
xx(n+1)=hxx(n)+μεx(n)Xout(n)・Xin *(n)
xy(n+1)=hxy(n)+μεx(n)Xout(n)・Yin *(n)
yx(n+1)=hyx(n)+μεy(n)Yout(n)・Xin *(n)
yy(n+1)=hyy(n)+μεy(n)Yout(n)・Yin *(n)
第1のタップ係数更新器72は、原則的には適応等化フィルタ71と同様のフィルタを有することができる。その場合、適応等化フィルタ71の出力Xout(n)及びYout(n)の供給を受けずに、当該第1のタップ係数更新器72の回路内でCMAの計算が可能である。勿論、適応等化フィルタ71の出力Xout(n)及びYout(n)の供給を受けることで、適応等化フィルタ71と同様のフィルタを有しない構成も可能である。
なお、適応等化フィルタ71のタップ係数を求めるアルゴリズムは、CMAに限定されず、ステップサイズμを使用する各種の逐次更新アルゴリズム、例えば、その他のブラインド等化方式(RDE(Radius directed equalization)等)も使用可能である。更に、RLS(Recursive Least-Squares)やLMS(Least Mean Square)等のように、送信側でトレーニング信号又はパイロット信号等の既知信号を光信号に挿入し、伝送されてきた既知信号とこの既知信号の真値との比較誤差を最小化するようにタップ係数をステップサイズ毎に更新して求める逐次更新アルゴリズムを用いることもできる。第1のタップ係数更新器72の可変のステップサイズμは、後述するステップサイズ更新器75から供給される。
第2のタップ係数更新器73は、固定のステップサイズμ0を用いて受信信号61(Xin(n)、Yin(n))の偏波変動に応じたタップ係数を計算する回路である。ここで計算されたタップ係数は、後述する偏波状態の推定に使用されるものであり、適応等化フィルタ71には供給されない。
第2のタップ係数更新器73は、適応等化フィルタ71のタップ係数を計算するように第1のタップ係数更新器72と同じように構成されるが、高速な偏波変動を観測できるようにステップサイズμとして比較的大きめの固定ステップサイズμ0を用いている。この固定ステップサイズμ0は、使用される逐次更新アルゴリズムにおいて種々の偏波変動において最適な収束が得られるように経験的に求められたステップサイズである。特許文献3で求められるステップサイズを使用することも可能である。固定ステップサイズμ0は、例えば、後述するステップサイズ更新器75から第1のタップ係数更新器72に設定される可変のステップサイズμの想定される値の平均値よりも大きい値、或いは可変のステップサイズμの想定される最大値と最小値の中間値よりも大きい値に設定することができる。
第2のタップ係数更新器73の計算にCMAを用いた場合、第1のタップ係数更新器72と同様に、以下の計算式によって、タップ係数が計算される。
xx(n+1)=hxx(n)+μ0εx(n)Xout(n)・Xin *(n)
xy(n+1)=hxy(n)+μ0εx(n)Xout(n)・Yin *(n)
yx(n+1)=hyx(n)+μ0εy(n)Yout(n)・Xin *(n)
yy(n+1)=hyy(n)+μ0εy(n)Yout(n)・Yin *(n)
このようにして計算されたタップ係数が、偏波状態推定器74に供給される。第2のタップ係数更新器73の計算では、固定ステップサイズμ0が用いられるので、後述するステップサイズ更新器75からのステップサイズμに拠らず、タップ係数演算における収束性を安定的に確保できる。なお、第2のタップ係数更新器73でも、第1のタップ係数更新器72と同様に、タップ係数を求めるアルゴリズムは、CMAに限定されず、ステップサイズを使用する各種の逐次更新アルゴリズムが使用可能である。
偏波状態推定器74は、第2のタップ係数更新器73で計算されたタップ係数を用いて、第2のタップ係数更新器73の入力信号の偏波状態SOPを推定する回路である。第2のタップ係数更新器73の入力信号は、適応等化フィルタ71及び第1のタップ係数更新器72の入力信号と同じく、受信信号61である。よって、偏波状態推定器74で推定された偏波状態SOPは、第1のタップ係数更新器に入力される受信信号61の偏波状態SOPでもある。
偏波状態SOPの推定方法には複数の方法がある。例えば特許文献1に開示された方法を使用することができる。すなわち、適応等化フィルタ71のタップ係数から偏波状態SOPを示すストークスパラメータを求め、そのストークスパラメータから偏波状態変動量及びSOP変動速度を求める。光信号の偏波状態SOPを表現するストークスパラメータは、以下の式により計算できる。
1=cos(2φ)cos(2ψ)
2=sin(2φ)cos(2ψ)
3=sin(2φ)
ここで、φ及びψは、X偏波について、次のように表される。
φx=tan-1(|hxy(f)|/|hxx(f)|),
ψx=tan-1(Im{hxy(f)/hxx(f)}/Re{hxy(f)/hxx(f)})
但し、hxx(f)及びhxy(f)は、周波数領域でのf=0における適応等化器70の係数のセントラルタップである。これはY偏波についても同様に適用できる。
ある時刻t1からt2におけるΔt=t2−t1秒後の偏波状態変動量Δθは、ストークスベクトルS(t1)=[S1(t1),S2(t1),S3(t1)]とS(t2)=[S1(t2),S2(t2),S3(t2)]の内積で求められる。このため、SOP変動速度ωは、次のように計算できる。
ω=Δθ/Δt=cos-1(S(t1)・S(t2)/(|S(t1)||S(t2)|))/Δt
偏波状態推定器74は、偏波状態SOPとしてSOP変動速度ωをステップサイズ更新器75に出力する。なお、偏波状態推定器74は、偏波状態SOPとして偏波状態変動量Δθを出力してもよい。Δtは設定値であるから、本明細書では、偏波状態変動量ΔθとSOP変動速度ωを同義として考える。
また、偏波状態SOPの推定方法として、特許文献2に開示された方法を使用することもできる。すなわち、適応等化フィルタ71のタップ係数からSOP変動速度ωを直接求める。
適応等化によって更新される適応等化器70のFIRフィルタ91〜94のタップ係数がストークスパラメータと関係性を有していることを利用する点は、特許文献1と同じである。特許文献2では更に、タップ係数の時間変化が、偏波状態の変動を示すことを利用している。偏波状態の変動S(n)は、タップ係数hxx,i、hxy,i、hyx,i、hyy,iを用いて、下式で示される。
Figure 0006876099
但し、MはFIRフィルタ91〜94のタップ数である。iはタップの順番を示す。Δhxx,i(n)、Δhxy,i、Δhyx,i、Δhyy,iは、タップ係数の時間変化量であり、それぞれ次のように表される。
Δhxx,i(n)=hxx,i(n)−hxx,i(n-1)
Δhxy,i(n)=hxy,i(n)−hxy,i(n-1)
Δhyx,i(n)=hyx,i(n)−hyx,i(n-1)
Δhyy,i(n)=hyy,i(n)−hyy,i(n-1)
xx(n)は、タップ係数ベクトルhxxの各要素の時間変化量の絶対値を加算した値である。同様に、Sxy(n)、Syx(n)、Syy(n)もタップ係数ベクトルhxy、hyx、hyyの各要素の時間変化量の絶対値を加算した値である。上式により変動S(n)を算出し、算出された変動S(n)の周期を算出する。変動S(n)の周期の逆数(周波数)を算出することにより、偏波状態のSOP変動速度ωを推定できる。あるいは、変動S(n)に対してフーリエ変換を行い、偏波状態のSOP変動速度ωを推定してもよい。推定されたSOP変動速度ωは、偏波状態推定器74の出力として、ステップサイズ更新器75に供給される。
偏波状態推定器74の出力のイメージを図5示す。これは、所定の偏波変動速度で偏波状態がランダムに変化するようにコントロールされた受信信号61(偏波スクランブラと称する)を生成して適応等化装置70に供給した場合に、偏波状態推定器74から出力される偏波状態変動量Δθを時間に対して描いたイメージである。この図によると、偏波状態変動量Δθは、所定の偏波変動速度の周期T毎にピークを発生することが分る。従って、このピークの周期Tを検出することでも、SOP変動速度ωは推定可能である。なお、SOP変動速度ωが速い場合、SOP変動速度ωが遅い場合と比較して、所定の時間幅における偏波状態変動量Δθのピークの数が多くなる傾向にあることが実験的に把握されている。
偏波状態推定器74における偏波状態SOPの推定方法は、適応等化フィルタ71のタップ係数を用いた方法であればどの方法も利用可能であり、上記の方法に限定されない。
ステップサイズ更新器75は、偏波状態推定器74で推定された偏波状態に対応したステップサイズμを求め、そのステップサイズμで第1のタップ係数更新器72の可変のステップサイズμを更新する回路である。要するに、ステップサイズ更新器75は、偏波状態推定器74から出力されるSOP変動速度ωに応じて、第1のタップ係数更新器72のCMAのステップサイズμを適応制御する。
特許文献1によると、あるSOP変動速度ωでSOP変動が生じている受信信号61に適応等化処理を適用した際に、符号誤り率BERが最小となる最適なステップサイズは、以下の関係式で近似できる。但し、A及びBは0でない定数である。
μ=(Aω)1/2+B
この関係式を用いることで、適応等化フィルタ71に入力される受信信号61のSOP変動速度ωに対応するステップサイズμを決定することが可能となる。なお、定数A、Bは、特許文献1では、デジタル信号処理によって適応的に求められることが示されている。上記の関係式によると、小さなSOP変動速度ωに対しては、比較的小さなステップサイズμが得られ、大きなSOP変動速度ωに対しては、比較的大きなステップサイズμが得られる。なお、上記の関係式の代わりに、近似式を用いることも可能である。多少の誤差は、許容可能である。
事前の測定により、偏波状態変動量Δθ及びSOP変動速度ωの少なくとも一方に対し、最適な誤り率特性が得られるステップサイズμの値を求め、その対応関係を図6に示すメモリ97に記憶させておいてもよい。これにより、ステップサイズ更新器75は、偏波状態推定器74の推定結果に基づいて、適切なステップサイズμを即座に選択することができる。
この場合、偏波状態変動量Δθ及び/又はSOP変動速度ωの所定範囲毎に、ステップサイズμを対応づけてもよい。言い換えると、偏波状態変動量Δθ及び/又はSOP変動速度ωに対してある程度幅を持たせてステップサイズμを記憶させてもよい。これにより、メモリ量を低減することができる。また、そのような幅を設けることで、第2のタップ係数更新器73のタップ係数の計算や偏波状態推定器74の偏波状態変動量Δθ及びSOP変動速度ωの計算精度を吸収(容認)することができる。
偏波状態推定器74から出力されるSOP変動速度ωが低速の場合(例えば0.1kHz程度以下)、タップ係数の更新毎に偏波状態推定器74の出力にバタつきが見られる。このため、偏波状態推定器74で推定される偏波状態変動量Δθ及びSOP変動速度ωの各々を平均化してステップサイズ更新器75に供給することが好ましい。この際、偏波状態変動量Δθ及びSOP変動速度ωの各々を、移動平均(例えば32ビット分)するとより好ましい。これにより、ステップサイズ更新器75によるステップサイズμの設定のバタつきを防ぐことができる。このような平均化を行うには、偏波状態推定器74に図7に示すような移動平均回路98を設ければよい。この移動平均回路98は、受信信号61の偏波状態の推定値の移動平均を求め、求められた平均値をステップサイズ更新器75に供給する回路である。偏波状態推定器74とステップサイズ更新器75との間に、移動平均回路98と同様の機能を有する回路を設けてもよい。
先に、偏波状態推定器74における偏波状態SOPの推定方法に関し、図5を参照して、偏波状態変動量Δθのピークの周期TからSOP変動速度ωを推定する方法について説明した。この偏波状態変動量Δθのピークの検出に際しては、ヒステリシス特性を持たせることが好ましい。例えば、図5において、偏波状態変動量Δθが指定期間内に設定閾値を指定回数超えたらピークとして検出し、設定閾値以下を指定回数下回るとピークの検出待ちの状態にする。もちろん、設定閾値及び指定回数はピーク検出状態とピーク検出待ち状態で必ずしも同じである必要はない。これにより、ピーク検出によるSOP変動速度ωの推定結果においても、ステップサイズ更新器75によるステップサイズμの設定においても、バタつきを防ぐことができる。このようなピーク検出を行うには、偏波状態推定器74に図7に示すようなピーク検出回路99を設ければよい。このピーク検出回路99は、推定された偏波状態変動量Δθからヒステリシス特性を利用してピークを検出する回路である。
[適応等化器の動作]
次に、本発明の実施の形態に係る適応等化器70の動作を比較例と比較して説明する。図8に示す適応等化器の比較例は、特許文献1に開示されているデジタル信号処理部を本実施の形態に係る適応等化器70の機能と整合させて書き直したものである。比較例では、偏波状態推定器174の推定計算に使用されるタップ係数が、ステップサイズ更新器175によってステップサイズμが適応制御されるタップ係数更新器175から供給される。それ以外の構成は、適応等化器70の構成と同じである。
比較例では、適応等化フィルタ171のタップ係数がタップ係数更新器172によって計算される。受信信号161の偏波状態の変動が小さくなると、タップ係数更新器172で計算されるタップ係数の変化も小さくなる。このタップ係数が偏波状態推定器174に供給され、偏波状態が推定される。この場合、タップ係数の変化が小さいため、偏波状態変動量ΔθやSOP変動速度ωは小さな値として推定される。ステップサイズ更新器175には小さな値として推定されたSOP変動速度ωが供給されるので、比較的小さなステップサイズが算出され、タップ係数更新器172のステップサイズμとして設定される。
結果的に、受信信号161の偏波状態変動量ΔθやSOP変動速度ωが小さいと推定された場合、比較的小さなステップサイズμが設定される。偏波変動が小さい場合は、高速の追従性を必要としないため、ステップサイズμは小さくてもよい。
次に、その後、受信信号161の偏波状態の変動が大きくなった場合について考える。上述したように比較的小さなステップサイズμがタップ係数更新器172に一旦設定されると、その後に偏波状態の変動が大きくなっても、ステップサイズμが小さいため偏波状態の変動に追従できず、適正なタップ係数が計算できない。偏波状態推定器174には適正なタップ係数が供給されないため、適正な偏波状態が推定できず、ステップサイズも適正な値が設定されない。その結果、タップ係数の更新の追従性が遅くなるという問題が発生する。
これに対し、本実施の形態では、適応等化フィルタ71のタップ係数が、まず第1タップ係数更新器72によって計算される。受信信号61の偏波状態の変動が小さくなると、第1のタップ係数更新器72で計算されるタップ係数の変化も小さくなる。また、同時に第2のタップ係数更新器73によってもタップ係数が計算され、受信信号61の偏波状態の変動が小さくなると、第2のタップ係数更新器73で計算されるタップ係数の変化も小さくなる。
この第2のタップ係数更新器73で計算されたタップ係数が偏波状態推定器74に供給され、偏波状態が推定される。この場合、タップ係数の変化が小さいため、偏波状態変動量ΔθやSOP変動速度ωは小さな値として推定される。ステップサイズ更新器75には小さな値として推定されたSOP変動速度ωが供給されるので、比較的小さなステップサイズが算出され、第1のタップ係数更新器72のステップサイズμとして設定される。
結果的に、受信信号61の偏波状態変動量ΔθやSOP変動速度ωが小さいと推定された場合、比較的小さなステップサイズμが設定される。偏波変動が小さい場合は、高速の追従性を必要としないため、ステップサイズμは小さくてもよい。
次に、その後、受信信号61の偏波状態の変動が大きくなった場合について考える。本実施の形態では、偏波状態は、第2のタップ係数更新器73で計算されたタップ係数から推定される。第2のタップ係数更新器73には予め比較的大きな固定のステップサイズμ0が設定されているため、偏波状態の変動の大小に関わらずタップ係数は比較的大きなステップサイズμ0で更新される。従って、第2のタップ係数更新器73は、タップ係数を高速に求めることができる。高速に求められたタップ係数が偏波状態推定器74に供給され、偏波状態も高速に推定でき、偏波状態の大きな変動に対応したステップサイズμが第1のタップ係数更新器72に設定できる。その結果、タップ係数の更新の安定的な追従性が確保できる。
なお、ステップサイズμ0が比較的大きな固定値のため、タップ係数の計算精度の多少の劣化が予想される。しかし、上述したとおり、ステップサイズ更新器75のメモリ97において、各ステップサイズに対応づけられるSOP変動速度ωに幅を持たせることで、タップ係数の計算精度の劣化を吸収することが可能となる。
以上示したように、本実施の形態では、第2のタップ係数更新器73が常に比較的大きい値に固定されたステップサイズを用いてタップ係数を計算し、偏波状態推定器74が第2のタップ係数更新器73によって更新されたタップ係数を用いて偏波状態を推定する。従って、偏波状態が大きく変動した場合でも、常に偏波状態を高速に追従でき、偏波状態の変動量に応じたステップサイズμをタップ係数更新の計算に使用できる。その結果、偏波状態の変動に追従した安定な適応等化を実現できる。
なお、図3に示した適応等化器70を構成するエレメント71〜75は、ハードウェアでディスクリートに構成することができる。また、エレメント71〜75の一部又は全部の機能を、コンピュータをプログラムに従って動作させることによって実現することも可能である。
10…送信信号処理器、11a…水平偏波用データ、11b…垂直偏波用データ、20…光送信機、21…信号光源(信号LD)、22a,22b…90゜合成器、23…偏波合成器、30…光ファイバ伝送路、40…光受信機、41…偏波分離器、42…局部発振光源(局発LD)、43a,43b…90゜ハイブリッド回路、44…光電変換器、45a…X偏波信号、45b…Y偏波信号、50…AD変換器、60…波長分散補償器、61…受信信号、70…適応等化器、71…適応等化フィルタ、72…第1のタップ係数更新器、73…第2のタップ係数更新器、74…偏波状態推定器、75…ステップサイズ更新器、80…復号器、91〜94…FIRフィルタ、95,96…加算回路、97…メモリ、98…移動平均回路、99…ピーク検出器。

Claims (10)

  1. タップ係数を変更することで受信信号の偏波変動を適応的に補償する適応等化フィルタと、
    可変のステップサイズを用いて前記受信信号の偏波変動に応じたタップ係数を計算し、計算された前記タップ係数で前記適応等化フィルタの前記タップ係数を自身で更新する第1のタップ係数更新器と、
    前記受信信号を入力とし前記適応等化フィルタの出力を受けることなく自身で、固定のステップサイズを用いて前記受信信号の偏波変動に応じたタップ係数を計算する第2のタップ係数更新器と、
    前記第2のタップ係数更新器で計算されたタップ係数を用いて前記受信信号の偏波状態を推定する偏波状態推定器と、
    前記偏波状態推定器で推定された偏波状態に対応したステップサイズを求め、求められた前記ステップサイズで前記可変のステップサイズを更新するステップサイズ更新器と
    を備えることを特徴とする適応等化器。
  2. 前記第1のタップ係数更新器及び前記第2のタップ係数更新器の各々は、伝送されてきた既知信号と前記既知信号の真値との比較誤差を最小化するように前記タップ係数をステップサイズ毎に更新して求める逐次更新アルゴリズムを用いて前記タップ係数を計算することを特徴とする請求項1に記載された適応等化器。
  3. 前記ステップサイズ更新器は、偏波状態変動量及び偏波状態変動速度の少なくとも一方の所定範囲毎に最適な誤り率特性が得られるステップサイズが予め記憶されているメモリを含むことを特徴とする請求項1又は2に記載された適応等化器。
  4. 前記偏波状態推定器は、
    前記受信信号の偏波状態の推定値の移動平均を求め、求められた平均値を前記ステップサイズ更新器に供給する移動平均回路、及び、
    推定された偏波状態変動量からヒステリシス特性を利用してピークを検出するピーク検出回路
    の少なくとも一方を含むことを特徴とする請求項1〜3の何れか1項に記載された適応等化器。
  5. タップ係数を変更することで受信信号の偏波変動を適応的に補償する適応等化フィルタリングを実行するステップと、
    可変のステップサイズを用いて前記受信信号の偏波変動に応じたタップ係数を計算し、計算された前記タップ係数で前記適応等化フィルタリングの前記タップ係数を自身で更新する第1のタップ係数更新ステップと、
    前記受信信号を入力とし前記適応等化フィルタの出力を受けることなく自身で、固定のステップサイズを用いて前記受信信号の偏波変動に応じたタップ係数を計算する第2のタップ係数更新ステップと、
    前記第2のタップ係数更新ステップで計算されたタップ係数を用いて前記受信信号の偏波状態を推定する偏波状態推定ステップと、
    前記偏波状態推定ステップで推定された偏波状態に対応したステップサイズを求め、求められた前記ステップサイズで前記可変のステップサイズを更新するステップサイズ更新ステップと
    を備えることを特徴とする適応等化方法。
  6. 前記第1のタップ係数更新ステップ及び前記第2のタップ係数更新ステップの各々は、伝送されてきた既知信号と前記既知信号の真値との比較誤差を最小化するように前記タップ係数をステップサイズ毎に更新して求める逐次更新アルゴリズムを用いて前記タップ係数を計算するステップを含むことを特徴とする請求項5に記載された適応等化方法。
  7. 前記ステップサイズ更新ステップは、偏波状態変動量及び偏波状態変動速度の少なくとも一方の所定範囲毎に最適な誤り率特性が得られるステップサイズが予め記憶されているメモリを使用して、ステップサイズを選択するステップを含むことを特徴とする請求項5又は6に記載された適応等化方法。
  8. 前記偏波状態推定ステップは、
    前記受信信号の偏波状態の推定値の移動平均を求め、求められた平均値を前記ステップサイズ更新ステップに供給するステップ、及び、
    推定された偏波状態変動量からヒステリシス特性を利用してピークを検出するステップ
    の少なくとも一方を含むことを特徴とする請求項5〜7の何れか1項に記載された適応等化方法。
  9. 受信信号の波長分散による歪を補償する波長分散補償器と、
    前記波長分散補償器の出力信号の偏波変動を補償する請求項1〜4の何れか1項に記載された適応等化器と、
    前記適応等化器の出力信号を復号する復号器と
    を備えることを特徴とする光通信システム。
  10. 光信号を受信し、受信された光信号を前記受信信号に変換する光受信機
    を更に備えることを特徴とする請求項9に記載された光通信システム。
JP2019115330A 2019-06-21 2019-06-21 適応等化器、適応等化方法及び光通信システム Active JP6876099B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2019115330A JP6876099B2 (ja) 2019-06-21 2019-06-21 適応等化器、適応等化方法及び光通信システム
PCT/JP2020/022780 WO2020255809A1 (ja) 2019-06-21 2020-06-10 適応等化器、適応等化方法及び光通信システム
US17/596,617 US12040834B2 (en) 2019-06-21 2020-06-10 Adaptive equalizer, adaptive equalization method, and optical communication system
CN202080044772.9A CN114026800B (zh) 2019-06-21 2020-06-10 自适应均衡器、自适应均衡方法和光通信系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019115330A JP6876099B2 (ja) 2019-06-21 2019-06-21 適応等化器、適応等化方法及び光通信システム

Publications (2)

Publication Number Publication Date
JP2021002743A JP2021002743A (ja) 2021-01-07
JP6876099B2 true JP6876099B2 (ja) 2021-05-26

Family

ID=73995204

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019115330A Active JP6876099B2 (ja) 2019-06-21 2019-06-21 適応等化器、適応等化方法及び光通信システム

Country Status (4)

Country Link
US (1) US12040834B2 (ja)
JP (1) JP6876099B2 (ja)
CN (1) CN114026800B (ja)
WO (1) WO2020255809A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230119766A1 (en) * 2020-04-07 2023-04-20 Nec Corporation Polarization-fluctuation estimating device and polarization-fluctuation estimating method
JP7055268B2 (ja) * 2020-05-28 2022-04-18 Nttエレクトロニクス株式会社 適応等化器、適応等化方法及び光通信システム
CN114884580A (zh) * 2021-02-05 2022-08-09 富士通株式会社 偏振变化的跟踪装置、接收信号的处理装置及方法
CN115811367B (zh) * 2022-12-06 2024-07-05 华中科技大学 一种基于硬件和软件结合的水下光通信系统的均衡方法
US20230214669A1 (en) * 2023-02-22 2023-07-06 Intel Corporation Hybrid search to train memory and high-speed input/output interfaces

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101425851B (zh) * 2008-12-06 2011-11-30 华中科技大学 光通讯用电子色散补偿均衡器及其抽头调整方法
US8737847B2 (en) 2012-03-19 2014-05-27 Futurewei Technologies, Inc. Method and apparatus of using joint timing recovery for a coherent optical system
JP5826702B2 (ja) * 2012-04-17 2015-12-02 日本電信電話株式会社 デジタルコヒーレント受信機およびデジタルコヒーレント受信方法
JP5628957B2 (ja) 2013-04-08 2014-11-19 日本電信電話株式会社 光受信器及び光伝送システム
CN103812806B (zh) * 2014-02-24 2017-10-10 北京大学 一种基于时域信道估计的信道均衡方法和系统
CN105530050B (zh) * 2014-10-21 2019-05-21 中兴通讯股份有限公司 均衡与偏振解复用和相偏估计与补偿的联合处理方法及装置
CN105245478B (zh) * 2015-10-14 2018-07-17 成都泰格微电子研究所有限责任公司 一种基于qam调制方式的自适应均衡算法
JP6659613B2 (ja) * 2017-03-31 2020-03-04 日本電信電話株式会社 偏波状態推定方法及び偏波状態推定装置
CN107248965B (zh) * 2017-05-26 2019-11-22 华中科技大学 一种数据处理方法及通信设备
JP7057506B2 (ja) * 2018-09-11 2022-04-20 日本電信電話株式会社 デジタルコヒーレント受信器及びデジタルコヒーレント受信方法
CN109450552B (zh) * 2018-11-15 2020-04-21 武汉邮电科学研究院有限公司 一种基于cma的信号均衡方法及系统

Also Published As

Publication number Publication date
US12040834B2 (en) 2024-07-16
JP2021002743A (ja) 2021-01-07
WO2020255809A1 (ja) 2020-12-24
CN114026800A (zh) 2022-02-08
US20220303012A1 (en) 2022-09-22
CN114026800B (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
JP6876099B2 (ja) 適応等化器、適応等化方法及び光通信システム
JP5826702B2 (ja) デジタルコヒーレント受信機およびデジタルコヒーレント受信方法
JP7057506B2 (ja) デジタルコヒーレント受信器及びデジタルコヒーレント受信方法
JP5943001B2 (ja) 等化信号処理装置、それを用いた光受信装置および等化信号処理方法
US20190074903A1 (en) Transmission characteristic compensation apparatus, transmission characteristic compensation method, and communication apparatus
US7187873B2 (en) Compensation of polarization mode dispersion in single mode fiber for maximum-likelihood sequence estimation
KR970009687B1 (ko) 디지탈통신을 위한 적응 등화시스템
JP5157536B2 (ja) トランバーサルフィルタのタップ係数制御装置及びタップ係数の制御方法、光受信装置並びに光通信システム
US9917712B2 (en) Frequency domain equalizer for optical communications
CN105530050B (zh) 均衡与偏振解复用和相偏估计与补偿的联合处理方法及装置
JP2015122632A (ja) 光通信受信装置
JP5965356B2 (ja) 信号処理システム、及び信号処理方法
CN109891842A (zh) 缓解自适应均衡与定时恢复之间的相互作用
JP6355465B2 (ja) 光受信器、送受信装置、光通信システムおよび波形歪補償方法
JP6865900B2 (ja) 光受信機、光信号受信方法及びデータ再生装置
JP5447381B2 (ja) 信号波形歪み補償器、及び信号波形歪み補償方法
CN114303329A (zh) 信号处理装置、信号处理方法和程序
JP7055268B2 (ja) 適応等化器、適応等化方法及び光通信システム
JP6464727B2 (ja) 信号処理装置、光受信装置、光通信システム、信号処理方法および光受信方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190621

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200908

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201102

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210423

R150 Certificate of patent or registration of utility model

Ref document number: 6876099

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE

Ref document number: 6876099

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350