特許 7701646 ネットワークコントローラ、推定方法及びコンピュータプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2025-06-24

(45)【発行日】2025-07-02

(54)【発明の名称】ネットワークコントローラ、推定方法及びコンピュータプログラム

(51)【国際特許分類】

   H04B 10/079 20130101AFI20250625BHJP

【ＦＩ】

H04B10/079 150

【請求項の数】 6

(21)【出願番号】P 2023546633

(86)(22)【出願日】2021-09-08

(86)【国際出願番号】 JP2021033039

(87)【国際公開番号】W WO2023037454

(87)【国際公開日】2023-03-16

【審査請求日】2024-02-19

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(74)【代理人】

【識別番号】110001634

【氏名又は名称】弁理士法人志賀国際特許事務所

(72)【発明者】

【氏名】笹井 健生

【審査官】対馬 英明

(56)【参考文献】

【文献】国際公開第２０１９／０４９６１６（ＷＯ，Ａ１）

【文献】特開２０１０－２７８５６２（ＪＰ，Ａ）

【文献】SASAI T. et al.，Simultaneous Detection of Anomaly Points and Fiber types in Multi-span Transmission Links Only by Receiver-side Digital Signal Processing，2020 Optical Fiber Communications Conference and Exhibition (OFC)，米国，OSA，2020年03月，Th1F.1

【文献】SASAI T. et al.，Digital Backpropagation for Optical Path Monitoring: Loss Profile and Passband Narrowing Estimation，2020 European Conference on Optical Communications (ECOC)，IEEE，2020年12月

【文献】SASAI T. et al.，Revealing Raman-amplified Power Profile and Raman Gain Spectra with Digital Backpropagation，2021 Optical Fiber Communications Conference and Exhibition (OFC)，米国，OSA，2021年03月，M3I.5

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｈ０４Ｂ １０／００－１０／９０

Ｈ０４Ｊ １４／００－１４／０８

(57)【特許請求の範囲】

【請求項1】

光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラであって、
前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定するシミュレーション部と、
前記シミュレーション部により推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性を推定して出力する理想出力推定部と、
を備え、
前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なり、
前記シミュレーション部により推定された前記デバイスの特性と、理想的なデバイスの特性を表す理想出力とを対応付けた１以上の第１データセットと、前記実空間で事前に取得した前記デバイスの特性と、理想的なデバイスの特性を表す理想出力とを対応付けた１以上の第２データセットとを記憶する記憶部をさらに備え、
前記理想出力推定部は、前記第１データセットと前記第２データセットとを任意の割合で混ぜて、前記理想出力と、前記デバイスの特性とに基づく演算を行うことで、二乗誤差を最小化するデジタルフィルタを構築し、構築した前記デジタルフィルタに、通信時に前記実空間で得られた前記光伝送路に設けられるデバイスの特性を入力することによって、前記理想的なデバイスの特性を推定するネットワークコントローラ。

【請求項2】

光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラであって、
前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定するシミュレーション部と、
前記シミュレーション部により推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性を推定して出力する理想出力推定部と、
を備え、
前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なり、
前記理想出力推定部は、前記シミュレーション部により推定された前記デバイスの特性それぞれと、前記実空間で得られた前記光伝送路に設けられるデバイスの特性との尤度を算出し、算出した尤度が最も高い前記シミュレーション部により推定されたデバイスの特性を、前記理想的なデバイスの特性として推定するネットワークコントローラ。

【請求項3】

前記光送信機から送信された光信号を受信した前記光受信機の出力に基づいて、前記光伝送路に設けられるデバイスの特性を推定するデジタル信号処理部をさらに備える、
請求項１又は２に記載のネットワークコントローラ。

【請求項4】

光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラが行う推定方法であって、
前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定し、
推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性を推定して出力し、
前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なり、
推定された前記デバイスの特性と、理想的なデバイスの特性を表す理想出力とを対応付けた１以上の第１データセットと、前記実空間で事前に取得した前記デバイスの特性と、理想的なデバイスの特性を表す理想出力とを対応付けた１以上の第２データセットとを任意の割合で混ぜて、前記理想出力と、前記デバイスの特性とに基づく演算を行うことで、二乗誤差を最小化するデジタルフィルタを構築し、構築した前記デジタルフィルタに、通信時に前記実空間で得られた前記光伝送路に設けられるデバイスの特性を入力することによって、前記理想的なデバイスの特性を推定する、推定方法。

【請求項5】

光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラが行う推定方法であって、
前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定し、
推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性を推定して出力し、
前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なり、
推定された前記デバイスの特性それぞれと、前記実空間で得られた前記光伝送路に設けられるデバイスの特性との尤度を算出し、算出した尤度が最も高い推定されたデバイスの特性を、前記理想的なデバイスの特性として推定する推定方法。

【請求項6】

請求項１から３のいずれか一項に記載のネットワークコントローラとしてコンピュータを機能するためのコンピュータプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、ネットワークコントローラ、推定方法及びコンピュータプログラムに関する。

【背景技術】

【0002】

光伝送システムには、光ファイバ、光アンプ、光フィルタ等の様々なデバイスが組み込まれている。従来では、これらのデバイス特性をＯＴＤＲ(Optical Time Domain Reflectometer)や光スペクトルアナライザ等のアナログ測定器により測定していた。しかしながら、アナログ測定器を用いた測定は、光ノード、光ファイバ毎に測定器が必要であり、設備コスト、運用コストが大きくなってしまう。

【0003】

そこで、近年、アナログ測定器による測定に代わり、受信側のデジタル信号処理により、光伝送システム内の様々なデバイスの特性を抽出する技術が提案されている（例えば、非特許文献１～４参照）。非特許文献１～４に記載の技術では、光伝送システムの受信信号に対し、デジタル信号処理のみで伝送システム中の様々なデバイスの応答を推定することが可能になる。例えば、非特許文献１～４に記載の技術では、光ファイバ損失又は分散分布、光アンプゲインスペクトル及び光フィルタ応答等が推定可能である。

【先行技術文献】

【非特許文献】

【0004】

【文献】T. Sasai, et al., “Simultaneous detection of anomaly points and fiber types in multi-span transmission links only by receiver-side digital signal processing”, in OFC2020, Paper Th1F.1.

【文献】T. Sasai, et al., “Physics-oriented learning of nonlinear Schrodinger equation: optical fiber loss and dispersion profile identification”, arXiv:2104.05890.

【文献】T. Sasai, et al., “Digital backpropagation for optical path monitoring: loss profile and passband narrowing estimation”, ECOC2020, Paper Tu2D.1.

【文献】T. Sasai, et al., “Revealing Raman-amplified power profile and Raman gain spectra with digital backpropagation”, in OFC2021, Paper M3I.5.

【発明の概要】

【発明が解決しようとする課題】

【0005】

しかしながら、非特許文献１～４に記載の技術では、アナログ測定器に比べ、簡易に測定ができるものの、デバイスの特性の推定精度や分解能が劣るという問題があった。

【0006】

上記事情に鑑み、本発明は、光伝送システムを構成するデバイスの特性を高精度に推定することができる技術の提供を目的としている。

【課題を解決するための手段】

【0007】

本発明の一態様は、光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラであって、前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定するシミュレーション部と、前記シミュレーション部により推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性の出力を推定する理想出力推定部と、を備え、前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なるネットワークコントローラである。

【0008】

本発明の一態様は、光送信機と、光受信機とが光伝送路を介して通信する光伝送システムに備えられるネットワークコントローラが行う推定方法であって、前記光送信機と、前記光受信機と、前記光伝送路とを含む光伝送システムにおける光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に前記光伝送路に設けられるデバイスの特性を推定し、推定された前記デバイスの特性と、実空間で得られた前記光伝送路に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性の出力を推定し、前記光伝送システムにおける光伝送路の状態は、前記シミュレーションモデル毎に異なる、推定方法である。

【0009】

本発明の一態様は、上記のネットワークコントローラとしてコンピュータを機能するためのコンピュータプログラムである。

【発明の効果】

【0010】

本発明により、光伝送システムを構成するデバイスの特性を高精度に推定することが可能となる。

【図面の簡単な説明】

【0011】

【図1】第１の実施形態における光伝送システムの概要を説明するための図である。

【図2】第１の実施形態における光伝送システムの構成を示す図である。

【図3】第１の実施形態におけるデジタル信号処理部の機能構成の一例を示す図である。

【図4】第１の実施形態における非線形光学補償部の機能構成を表す概略ブロック図である。

【図5】第１の実施形態におけるネットワークコントローラの処理の流れを示すフローチャートである。

【図6】第２の実施形態における光伝送システムの概要を説明するための図である。

【図7】第２の実施形態における光伝送システムの構成を示す図である。

【図8】第２の実施形態におけるネットワークコントローラの処理の流れを示すフローチャートである。

【図9】デジタル信号処理部の機能構成の別例を示す図である。

【発明を実施するための形態】

【0012】

以下、本発明の一実施形態を、図面を参照しながら説明する。
（第１の実施形態）
図１は、第１の実施形態における光伝送システムの概要を説明するための図である。第１の実施形態では、光伝送システムを構成するデバイスの特性の推定精度を向上させる方法として、デジタルフィルタを用いる。なお、以下の説明では、デバイスの特性として、光ファイバの損失分布を例に説明するが、光ファイバの分散分布、アンプゲインスペクトル及び光フィルタ応答についても同様の手法が適用可能である。以下の説明では、非特許文献１～４に記載の技術をＤＬＭ（Digital longitudinal monitoring）と記載する。

【0013】

デバイスの特性の推定精度向上のために、ＤＬＭにより得られる出力（以下「ＤＬＭ出力」という。）を理想的な出力に変換するデジタルフィルタを設計することが考えられる。しかし、あらゆる伝送路デバイスの特性を精度良く出力するフィルタを構成するには、あらゆる光伝送路を実際に用意し、伝送路デバイス特性のデータセット(ＤＬＭ出力と理想出力のセット)を膨大な量用意することが必要であり、現実的ではない。

【0014】

そこで、第１の実施形態では、デジタルフィルタを構成するためのデータセット不足を解消するため、図１の下図に示す通り、光伝送システムのシミュレーションモデル(光伝送路のデジタルツイン)を構築する。なお、光伝送路モデルは豊富に研究されており、例えば下記参考文献１の技術が利用できる。光フィルタモデルには、使用する光フィルタに合わせて任意のフィルタ形状を使用可能であり、例えばガウシアンフィルタが用いられてもよい。
（参考文献１：E. Ip et al., “Compensation of Dispersion and Nonlinear Impairments Using Digital Backpropagation”, J. Lightw. Technol., 26(20),2008）

【0015】

第１の実施形態では、実際の光伝送システムのように、デジタルツイン上であらゆる光伝送路の状態(損失量、損失位置、分散量etc等)を仮想的に再現し、ＤＬＭ出力を得て、ＤＬＭ出力と理想出力のデータセットを複数保持する。ここで、理想出力とは、デジタルツイン上でユーザが与えた光ファイバの損失分布等(ＤＬＭの推定対象が光ファイバの分散分布、アンプゲインスペクトル及び光フィルタ応答等の場合はそれら)のことであり、ＤＬＭの推定対象である伝送路の状態のことである。そして、保持している複数のデータセットに基づいてデジタルフィルタを構築することによって、上記のデータセット不足を解消する。これにより、結果としてデバイスの特性を精度よく出力するデジタルフィルタを構築することができる。以下、具体的な構成について説明する。

【0016】

図２は、第１の実施形態における光伝送システム１００の構成を示す図である。光伝送システム１００は、光送信機１と、光受信機２と、光伝送路３と、ネットワークコントローラ６とを備える。光送信機１と、光受信機２とは、光伝送路３を介して通信可能に接続される。光伝送路３は、例えば光ファイバＦで構成される。光伝送路３には、伝搬中に減衰する光信号を増幅して中継する１以上の中継ノード３１Ａが光増幅器（以下「光アンプ」ともいう。）として伝送路の途中に挿入される。なお、光伝送路３には、光フィルタが設けられてもよい。

【0017】

光送信機１は、外部の情報源から与えられる送信情報を符号化して電気信号を生成し、生成した電気信号を光信号に変換して光伝送路３を介して光受信機２に送信する。

【0018】

光受信機２は、コヒーレント受信器２１と、デジタル信号処理部２２とを備える。コヒーレント受信器２１は、ベースバンド光信号を偏波面が直交する２つの光信号に分離する。これらの光信号と局発光源（不図示）の局発光が９０°ハイブリッド回路（不図示）に入力され、両光を互いに同相及び逆相で干渉させた１組の出力光、直交（９０°）及び逆直交（－９０°）で干渉させた１組の出力光の計４つの出力光が得られる。これらの出力光はフォトダイオード（不図示）によりそれぞれアナログ信号に変換される。コヒーレント受信器２１は、これらのアナログ信号をデジタル信号に変換する。

【0019】

光伝送路３を光信号が伝搬する際に、信号の光パワーに比例して信号の位相が回転する非線形光学効果によって信号波形が歪む。デジタル信号処理部２２は、コヒーレント受信器２１が出力するデジタル信号を受信信号として取り込み、取り込んだ受信信号に対して非線形光学補償を行う。

【0020】

デジタル信号処理部２２は、受信された光信号をもとに光伝送路３におけるデバイスの特性を推定してＤＬＭ出力としてネットワークコントローラ６に通知する。例えば、デジタル信号処理部２２は、伝搬方向における光信号の強度分布を示す信号パワープロファイルをＤＬＭ出力として生成する。デジタル信号処理部２２は、生成したＤＬＭ出力をネットワークコントローラ６に供給する。

【0021】

ネットワークコントローラ６は、実空間の光伝送システムを構成する光送信機１と光受信機２と光伝送路３とを含む複数のシミュレーションモデルを構築する。この際、ネットワークコントローラ６は、あらゆる光伝送路の状態(損失量、損失位置、分散量etc等)を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築する。ネットワークコントローラ６は、構築したシミュレーションモデルそれぞれにより得られた仮想ＤＬＭ出力と、仮想ＤＬＭ出力に対応する理想出力とを含むデータセットと、光受信機２から得られたＤＬＭ出力とに基づいて、理想的なＤＬＭ出力を推定する。

【0022】

図３は、第１の実施形態におけるデジタル信号処理部２２の機能構成の一例を示す図である。デジタル信号処理部２２は、非線形光学補償部２３、適応等化部２４、周波数オフセット補償部２５、キャリア位相雑音補償部２６、係数更新部２７、および伝送特性推定部２８を備える。

【0023】

図４は、第１の実施形態における非線形光学補償部２３の機能構成を表す概略ブロック図である。非線形光学補償部２３は、複数の線形補償部２３１－１～２３１－Ｎ及び複数の非線形補償部２３２－１～２３２－Ｎを備える。１つの線形補償部２３１及び１つの非線形補償部２３２が線形補償及び非線形補償を行う１つのセットであり、非線形光学補償部２３はこのセットによる処理をＮステップ行うためにＮ個のセットを備えている。

【0024】

線形補償部２３１－１は、フーリエ変換部２３３－１、波長分散補償部２３４－１及び逆フーリエ変換部２３５－１を備える。フーリエ変換部２３３－１は、時間領域の受信信号に対してＦＦＴを行うことによって、時間領域の受信信号から周波数領域の受信信号に変換する。

【0025】

波長分散補償部２３４－１は、周波数領域の受信信号に対して、所定の値（例えば、ｅｘｐ＾（－ｊβ_ｋω_２））を乗算することによって波長分散補償を行う。なお、“＾”の記号は、“＾”以降の値がｅｘｐの上付きであることを意味する。例えば、ｅｘｐ＾（－ｊβ_ｋω_２）の場合には、（－ｊβ_ｋω_２）がｅｘｐの上付きであることを意味する。“＾”については以降の説明においても同様である。波長分散補償部２３４－１は、処理開始時には初期値として設定された分散係数β_ｋを用いて波長分散補償を行い、係数更新部２７から分散係数β_ｋが更新される度に更新後の分散係数β_ｋを用いて波長分散補償を行う。

【0026】

逆フーリエ変換部２３５－１は、波長分散補償部２３４－１から出力された信号に対してＩＦＦＴを行うことによって、波長分散補償された受信信号を時間領域の受信信号に変換する。
非線形光学補償部２３は、逆フーリエ変換部２３５から出力された信号系列に対して、所定の値（例えば、ｅｘｐ＾（－ｊφ_ｋ））を乗算することによって非線形光学効果の補償を行う。具体的には、非線形光学補償部２３は、処理開始時には初期値として設定された位相回転量φ_ｋを用いて非線形光学効果の補償を行い、係数更新部２７から位相回転量φ_ｋが更新される度に更新後の位相回転量φ_ｋを用いて非線形光学効果の補償を行う。

【0027】

線形補償部２３１－Ｎは、線形補償部２３１－１と同様の処理を行う。また、非線形補償部２３２－Ｎは、非線形補償部２３２－１と同様の処理を行う。

【0028】

図３に戻って、光受信機２の説明を続ける。適応等化部２４は、光伝送路３において光信号の波形に生じた歪みを補償する機能部である。すなわち、適応等化部２４は、光伝送路３において符号間干渉（シンボル間干渉）によって光信号に生じた符号誤りを訂正する機能部である。適応等化部２４は、設定されたタップ係数に応じて、ＦＩＲフィルタ（有限インパルス応答フィルタ）によって適応等化処理を実行する。

【0029】

周波数オフセット補償部２５は、適応等化処理が実行された４つのデジタル信号に対して、周波数オフセットを補償する処理を実行する。

【0030】

キャリア位相雑音補償部２６は、周波数オフセットが補償された４つのデジタル信号に対して、位相オフセットを補償する処理を実行する。

【0031】

係数更新部２７は、非線形光学補償部２３において使用される全係数（例えば、分散係数β_ｋ、位相回転量φ_ｋ等）を全ステップにおいて更新する。第１実施形態では、例えば係数更新部２７は、キャリア位相雑音補償部２６からの出力信号と、トレーニング信号とに基づいて、非線形光学補償部２３において使用される全係数（例えば、分散係数β_ｋ、位相回転量φ_ｋ等）を全ステップにおいて更新する。係数更新部２７は、更新した係数を、非線形光学補償部２３の各機能部に設定する。第１実施形態において出力信号と比較されるトレーニング信号は、電気信号に変換された送信信号である。

【0032】

伝送特性推定部２８は、光伝送路３の伝送特性を推定する。例えば、伝送特性推定部２８は、最適化された位相回転量φ_ｋを用いて損失分布を推定する。伝送特性推定部２８は、推定した損失分布を示す情報をＤＬＭ出力としてネットワークコントローラ６に供給する。なお、伝送特性推定部２８は、最適化された分散係数β_ｋを用いて分散分布を推定してもよい。この場合、伝送特性推定部２８は、分散分布を示す情報をＤＬＭ出力としてネットワークコントローラ６に供給する。

【0033】

次に、図２に戻ってネットワークコントローラ６の具体的な構成について説明する。ネットワークコントローラ６は、シミュレーション部６１、記憶部６２及び理想出力推定部６３を備える。

【0034】

シミュレーション部６１は、光伝送システム１００におけるあらゆる光伝送路３の状態(損失量、損失位置、分散量etc等)を仮想的に再現するための１以上のシミュレーションモデルを構築する。シミュレーション部６１は、構築した１以上のシミュレーションモデルそれぞれにより、１以上の仮想ＤＬＭ出力と、１以上の仮想ＤＬＭ出力それぞれに対応する理想出力を取得する。

【0035】

記憶部６２には、シミュレーション部６１により得られた、１以上の仮想ＤＬＭ出力と、１以上の仮想ＤＬＭ出力それぞれに対応する理想出力とを対応付けた１以上の第１データセットが記憶される。さらに、記憶部６２には、実環境で得られた１以上のＤＬＭ出力と、１以上のＤＬＭ出力それぞれに対応する理想出力とを対応付けた１以上の第２データセットが記憶される。理想出力は、人手で対応付けられてもよいし、シミュレーション部６１により自動的に対応付けられてもよい。記憶部６２は、磁気記憶装置や半導体記憶装置などの記憶装置を用いて構成される。

【0036】

理想出力推定部６３は、フィルタ構築部６３１と、理想出力部６３２とで構成される。フィルタ構築部６３１は、記憶部６２に記憶されている第１データセット及び第２データセットを用いてデジタルフィルタを構築する。理想出力部６３２は、フィルタ構築部６３１により構築されたデジタルフィルタと、光受信機２から得られたＤＬＭ出力とに基づいて光伝送システム１００の理想的なデバイスの特性の出力を推定する。

【0037】

図５は、第１の実施形態におけるネットワークコントローラ６の処理の流れを示すフローチャートである。なお、図５の処理の開始時には、シミュレーション部６１により得られたシミュレーション結果が記憶部６２に記憶されているものとする。
フィルタ構築部６３１は、記憶部６２に記憶されている第１データセット及び第２データセットを用いてデジタルフィルタを構築する（ステップＳ１０１）。具体的には、フィルタ構築部６３１は、第１データセットと、第２データセットを任意の割合で混ぜて、例えば以下の式（１）を元に変形した式（２）に基づいて、二乗誤差｜｜ｙ－Ａｈ｜｜^２を最小化する最適なＦＩＲフィルタｈを構築する。

【0038】

【数1】

【0039】

【数2】

【0040】

なお、ここで使用するデジタルフィルタは任意である。上記の例では、ＦＩＲフィルタのような線形フィルタを構築する構成を示したが、ボルテラフィルタやニューラルネットワークのような非線形フィルタを用いてもよい。ボルテラフィルタやニューラルネットワークのような非線形フィルタを用いた方がフィルタの表現力が向上し、より精度の高いフィルタを構成できる可能性がある。

【0041】

なお、デジタルフィルタの構築の方法も任意である。上記の例のような最小二乗法を用いた最適フィルタでも良いし、勾配法を用いて最適化されたフィルタでもよい。勾配法における目的関数も二乗誤差以外ものを用いてもよい。例えば式（２）内のＡ^ＴＡに任意の正則化項λＲを加えても良い。Ｒは任意の行列である。

【0042】

デジタルフィルタが構築された後、光伝送システム１００のユーザは、光送信機１と光受信機２との通信を開始させる。これにより、光受信機２では、光送信機１から送信された光信号に対してデジタル信号処理を行うことでＤＬＭ出力を取得することができる。光受信機２は、取得したＤＬＭ出力をネットワークコントローラ６に出力する。ネットワークコントローラ６の理想出力部６３２は、光受信機２により得られたＤＬＭ出力を入力とする（ステップＳ１０２）。

【0043】

理想出力部６３２は、入力したＤＬＭ出力を、フィルタ構築部６３１により構築されたデジタルフィルタに入力することで理想的なデバイスの特性の出力を推定する（ステップＳ１０３）。例えば、デジタルフィルタが高精度なフィルタである場合、入力されたＤＬＭ出力から雑音が除去されて、デバイスの特性がより明確に表れた出力結果を得ることができる。理想出力部６３２は、推定した理想的なデバイスの特性の出力の情報を外部に出力する。

【0044】

以上のように構成された光伝送システム１００によれば、光送信機１と、光受信機２と、光伝送路３とを含む光伝送システム１００における複数の光伝送路の状態を異なる仮想空間に仮想的に再現することで複数のシミュレーションモデルを構築して、シミュレーションモデル毎に光伝送路３に設けられるデバイスの特性を推定するシミュレーション部６１と、シミュレーション部６１により推定されたデバイスの特性と、実空間で得られた光伝送路３に設けられるデバイスの特性とに基づいて、理想的なデバイスの特性の出力を推定する理想出力推定部６３と備える。これにより、あらゆる光伝送路の状態を踏まえて、理想的なデバイスの特性の出力を推定することができる。そのため、伝送システムを構成するデバイスの特性を高精度に推定することが可能になる。

【0045】

さらに、光伝送システム１００では、シミュレーション部６１によるシミュレーションによってあらゆる伝送路デバイスの特性を推定している。そして、光伝送システム１００では、第１データセットと、第２データセットとを用いてデジタルフィルタを構築し、構築したデジタルフィルタに、実空間で得られた光伝送路に設けられるデバイスの特性の出力を入力することによって、理想的なデバイスの特性の出力を推定する。上述したように、デバイスの特性の推定精度向上のためには、あらゆる伝送路デバイスの特性を精度良く出力するフィルタを構成する必要があり、あらゆる光伝送路を実際に用意し、伝送路デバイス特性のデータセットを膨大な量用意することは現実的ではない。それに対して、本発明では、シミュレーションによってあらゆる伝送路デバイスの特性を推定しているため、より多くのデータセットを用意することができる。その結果、デバイスの特性を精度よく出力するデジタルフィルタを構築することができる。そのため、光伝送システム１００では、構築したデジタルフィルタを用いることにより、光伝送システムを構成するデバイスの特性を高精度に推定することが可能になる。

【0046】

（第２の実施形態）
第２の実施形態では、光伝送システムを構成するデバイスの特性の推定精度を向上させる方法として、最尤推定の手法を用いる。

【0047】

図６は、第２の実施形態における光伝送システムの概要を説明するための図である。第２の実施形態においても、第１の実施形態と同様に、デジタルツイン上であらゆる光伝送路の状態(損失量、損失位置、分散量etc等)を仮想的に再現し、ＤＬＭ出力を得て、ＤＬＭ出力と理想出力のデータセットを複数保持する。そして、ネットワークコントローラにおいて、保持している複数のＤＬＭ出力と、実際に得られたＤＬＭ出力との尤度をそれぞれ算出し、最大の尤度を与えるＤＬＭ出力を特定する。特定したＤＬＭ出力を出力した光伝送路の状態（理想出力）が、求める伝送路の状態となる。以下、具体的な構成について説明する。

【0048】

図７は、第２の実施形態における光伝送システム１００ａの構成を示す図である。光伝送システム１００ａは、光送信機１と、光受信機２と、光伝送路３と、ネットワークコントローラ６ａとを備える。光伝送システム１００ａは、ネットワークコントローラ６ａの構成以外は光伝送システム１００と同様である。そのため、第１の実施形態との相違点について説明する。

【0049】

ネットワークコントローラ６ａは、シミュレーション部６１、記憶部６２ａ及び理想出力推定部６３ａを備える。

【0050】

記憶部６２ａには、シミュレーション部６１により得られた、１以上の仮想ＤＬＭ出力と、１以上の仮想ＤＬＭ出力それぞれに対応する理想出力が記憶される。記憶部６２ａは、磁気記憶装置や半導体記憶装置などの記憶装置を用いて構成される。

【0051】

理想出力推定部６３ａは、尤度算出部６３３と、理想出力部６３２ａとで構成される。尤度算出部６３３は、１以上の仮想ＤＬＭ出力それぞれと、光受信機２から得られたＤＬＭ出力とを尤度を算出する。理想出力部６３２ａは、尤度算出部６３３により算出された複数の尤度のうち、最大の尤度を与える仮想ＤＬＭ出力を特定する。理想出力部６３２ａは、特定した仮想ＤＬＭ出力に対応する光伝送路の状態(理想出力)を、光伝送システム１００ａの理想的なデバイスの特性の出力として推定する。

【0052】

図８は、第２の実施形態におけるネットワークコントローラ６ａの処理の流れを示すフローチャートである。なお、図８の処理の開始時には、シミュレーション部６１により得られたシミュレーション結果が記憶部６２に記憶されているものとする。
尤度算出部６３３は、光受信機２により得られたＤＬＭ出力を入力とする（ステップＳ２０１）。尤度算出部６３３は、入力したＤＬＭ出力と、記憶部６２ａに記憶されている仮想ＤＬＭ出力それぞれとの尤度を算出する（ステップＳ２０２）。例えば、尤度算出部６３３は、入力したＤＬＭ出力と、記憶部６２ａに記憶されている仮想ＤＬＭ出力それぞれとの尤度を、以下の式（３）に基づいて算出してもよい。

【0053】

【数3】

【0054】

式（３）において、ａはＤＬＭ出力を表し、ａ´は仮想ＤＬＭ出力を表す。なお、尤度算出部６３３は、式（３）に示す二乗誤差（負号つき）に限らず、別の方法で尤度を算出してもよい。

【0055】

尤度算出部６３３は、算出した各尤度の情報を理想出力部６３２ａに出力する。理想出力部６３２ａは、尤度算出部６３３から出力された各尤度の情報の中から、最大の尤度を与える仮想ＤＬＭ出力を特定する。理想出力部６３２ａは、特定した仮想ＤＬＭ出力に対応する光伝送路の状態を、光伝送システム１００ａの理想的なデバイスの特性の出力として推定する（ステップＳ２０３）。理想出力部６３２ａは、推定した理想的なデバイスの特性の出力の情報を外部に出力する。

【0056】

以上のように構成された光伝送システム１００ａでは、まずシミュレーション部６１により推定されたデバイスの特性それぞれと、実空間で得られた光伝送路に設けられるデバイスの特性との尤度を算出する。そして、光伝送システム１００ａでは、算出した複数の尤度のうち、最大の尤度を与える仮想ＤＬＭ出力を特定する。そして、光伝送システム１００ａでは、特定した仮想ＤＬＭ出力に対応する光伝送路の状態(理想出力)を、光伝送システム１００ａの理想的なデバイスの特性の出力として推定する。このように、実空間で得られたＤＬＭ出力に近い光伝送路の状態を推定することができる。そのため、伝送システムを構成するデバイスの特性を高精度に推定することが可能になる。

【0057】

（第１の実施形態及び第２の実施形態に共通する変形例）
上述した構成では、ＤＬＭ出力の計算を光受信機２で行う構成を示したが、ＤＬＭ出力の計算をネットワークコントローラ６，６ａで行うように構成されてもよい。このように構成される場合、光受信機２のデジタル信号処理部２２の機能がネットワークコントローラ６，６ａに実装される。

【0058】

各実施形態においてＤＬＭ出力を得るための構成は、上述したデジタル信号処理部２２の構成に限定される必要はない。デジタル信号処理部２２は、ＤＬＭ出力を得ることができる構成であれば他の構成であってもよい。図９は、デジタル信号処理部２２ａの機能構成の別例を示す図である。デジタル信号処理部２２ａは、波長分散補償部２９、適応等化部２４、周波数オフセット補償部２５、キャリア位相雑音補償部２６、波長分散付加部３０、非線形光学補償部２３、係数更新部２７、及び伝送特性推定部２８を備える。デジタル信号処理部２２ａは、デジタル信号処理部２２の構成に、波長分散補償部２９及び波長分散付加部３０をさらに備える。

【0059】

適応等化部２４、周波数オフセット補償部２５、キャリア位相雑音補償部２６、非線形光学補償部２３、係数更新部２７、及び伝送特性推定部２８の機能は、図３と同じであるため説明を省略する。波長分散補償部２９は、光伝送路３で生じた波長分散による歪みを補償する。具体的には、波長分散補償部２９は、周波数領域の受信信号に対して、所定の値（例えば、ｅｘｐ＾（－ｊβ_ｋω^２））を乗算することによって波長分散補償を行う。波長分散付加部３０は、波長分散補償部２９で補償した波長分散を再度、信号に付加する機能部である。具体的には、波長分散付加部３０は、周波数領域の受信信号に対して、波長分散補償部２９で乗算した値の逆数（例えば、ｅｘｐ＾（+ｊβ_ｋω^２））を乗算することによって波長分散の付加を行う。

【0060】

上述した光送信機１、光受信機２、及びネットワークコントローラ６，６ａの一部または全部の機能をコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。

【0061】

さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。

【0062】

以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。

【産業上の利用可能性】

【0063】

本発明は、光伝送システムを構成するデバイスの特性を推定する技術に適用できる。

【符号の説明】

【0064】

１００，１００ａ…光伝送システム、１…光送信機、２…光受信機、２１…コヒーレント受信器、２２…デジタル信号処理部、２３…非線形光学補償部、２３１…線形補償部、２３２…非線形補償部、２３３…フーリエ変換部、２３４…波長分散補償部、２３５…逆フーリエ変換部、２４…適応等化部、２５…周波数オフセット補償部、２６…キャリア位相雑音補償部、２７…係数更新部、２８…伝送特性推定部、２９…波長分散補償部２９、３０…波長分散付加部、３…光伝送路、３１Ａ…中継ノード、６，６ａ…ネットワークコントローラ、６１…シミュレーション部、６２，６２ａ…記憶部、６３，６３ａ…理想出力推定部、６３１…フィルタ構築部、６３２，６３２ａ…理想出力部、６３３…尤度算出部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】