光ファイバ伝送システム及び光受信装置
【課題】本発明は、マルチモード光ファイバ伝送における長距離伝送を実現する技術を提供することを目的とする。
【解決手段】本発明は、光信号を送信する１個の光送信機１と、光信号を受信するＭ個（Ｍは２以上の整数）の光受信機２と、１個の光送信機１が送信した光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバ３と、光ファイバ３が伝搬した光信号をＭ個の光受信機２に分波する分波器４と、各光受信機２が受信した光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償するＦＩＲフィルタ装置５と、を備えることを特徴とする光ファイバ伝送システムである。
【解決手段】本発明は、光信号を送信する１個の光送信機１と、光信号を受信するＭ個（Ｍは２以上の整数）の光受信機２と、１個の光送信機１が送信した光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバ３と、光ファイバ３が伝搬した光信号をＭ個の光受信機２に分波する分波器４と、各光受信機２が受信した光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償するＦＩＲフィルタ装置５と、を備えることを特徴とする光ファイバ伝送システムである。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光ファイバ伝送の群遅延差を補償する技術に関する。
【背景技術】
【０００２】
光ファイバ伝送システムでは、光ファイバ中で発生する非線形効果やファイバヒューズが問題となり、伝送の大容量化が制限されている。これらの制限を緩和するためには、光ファイバに導波する光の密度を低減する必要があり、非特許文献１、２に示すように大コアファイバが検討されている。
【０００３】
しかし、曲げ損失低減、単一モード動作領域の拡大、実効断面積の拡大は互いにトレードオフの関係にあり、所定の条件下における実効断面積の拡大量には限界があるという課題があった。そこで、伝送路をマルチモード光ファイバとし、Ｆｉｎｉｔｅ ｉｎｐｕｌｓｅ ｒｅｓｐｏｎｓｅ（ＦＩＲ）フィルタと適応等化技術を用いて、光ファイバ中で発生するモード分散を補償する試みが行われている（例えば、非特許文献３参照）。
【０００４】
これにより、先に述べた大コア光ファイバで制限要因であった単一モード動作条件が不要になるため、さらなる大コア化及び大容量化が可能である。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Ｔ．Ｍａｔｓｕｉ，Ｋ．Ｎａｋａｊｉｍａ ａｎｄ Ｃ．Ｆｕｋａｉ，“Ａｐｐｌｉｃａｂｉｌｉｔｙ ｏｆ Ｐｈｏｔｏｎｉｃ Ｃｒｙｓｔａｌ Ｆｉｂｅｒ Ｗｉｔｈ Ｕｎｉｆｏｒｍ Ａｉｒ−Ｈｏｌｅ Ｓｔｒｕｃｔｕｒｅ ｔｏ Ｈｉｇｈ−Ｓｐｅｅｄ ａｎｄ Ｗｉｄｅ−Ｂａｎｄ Ｔｒａｎｓｍｉｓｓｉｏｎ Ｏｖｅｒ Ｃｏｎｖｅｎｔｉｏｎａｌ Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ Ｂａｎｄｓ”，Ｊ．Ｌｉｇｈｔｗａｖｅ Ｔｅｃｈｎｏｌ．２７，５４１０−５４１６，２００９．
【非特許文献２】Ｋ．Ｍｕｋａｓａ，Ｋ．Ｉｍａｍｕｒａ，Ｒ．Ｓｕｇｉｚａｋｉ ａｎｄ Ｔ．Ｙａｇｉ，“Ｃｏｍｐａｒｉｓｏｎｓ ｏｆ ｍｅｒｉｔｓ ｏｎ ｗｉｄｅ−ｂａｎｄ ｔｒａｎｓｍｉｓｓｉｏｎ ｓｙｓｔｅｍｓ ｂｅｔｗｅｅｎ ｕｓｉｎｇ ｅｘｔｒｅｍｅｌｙ ｉｍｐｒｏｖｅｄ ｓｏｌｉｄ ＳＭＦｓ ｗｉｔｈ Ａｅｆｆ ｏｆ １６０μｍ２ ａｎｄ ｌｏｓｓ ｏｆ ０．１７５ｄＢ／ｋｍ ａｎｄ ｕｓｉｎｇ ｌａｒｇｅ−Ａｅｆｆ ｈｏｌｅｙ ｆｉｂｅｒｓ ｅｎａｂｌｉｎｇ ｔｒａｎｓｍｉｓｓｉｏｎ ｏｖｅｒ ６００ｎｍ ｂａｎｄｗｉｄｔｈ”，ｔｈｅ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＯＦＣ２００８，ＯｔｈＲ１，Ｆｅｂ．２００８．
【非特許文献３】Ｘ．Ｚｈａｏ ａｎｄ Ｆ．Ｓ．Ｃｈｏａ，“Ｄｅｍｏｎｓｔｒａｔｉｏｎ ｏｆ １０−Ｇｂ／ｓ ｔｒａｎｓｍｉｓｓｉｏｎｓ ｏｖｅｒ ａ １．５−ｋｍ−ｌｏｎｇ ｍｕｌｔｉｍｏｄｅ ｆｉｂｅｒ ｕｓｉｎｇ ｅｑｕａｌｉｚａｔｉｏｎ ｔｅｃｈｎｉｑｕｅｓ”，ＩＥＥＥ Ｐｈｏｔｏｎｉｃｓ Ｔｅｃｈｎｏｌｏｇｙ Ｌｅｔｔｅｒｓ，ｖｏｌ．１４，ｐｐ．１１８７−１１８９，２００２．
【非特許文献４】Ｍ．Ｔａｙｌｏｒ，“Ｃｏｈｅｒｅｎｔ Ｄｅｔｅｃｔｉｏｎ ｆｏｒ Ｆｉｂｅｒ Ｏｐｔｉｃ Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ Ｕｓｉｎｇ Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ”，ｉｎ Ｏｐｔｉｃａｌ Ａｍｐｌｉｆｉｅｒｓ ａｎｄ Ｔｈｅｉｒ Ａｐｐｌｉｃａｔｉｏｎｓ／Ｃｏｈｅｒｅｎｔ Ｏｐｔｉｃａｌ Ｔｅｃｈｎｏｌｏｇｉｅｓ ａｎｄ Ａｐｐｌｉｃａｔｉｏｎｓ，Ｔｅｃｈｎｉｃａｌ Ｄｉｇｅｓｔ（ＣＤ）（Ｏｐｔｉｃａｌ Ｓｏｃｉｅｔｙ ｏｆ Ａｍｅｒｉｃａ，２００６），ｐａｐｅｒ ＣＴｈＢ１．
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、伝送距離の長距離化に伴ってモード間の群遅延差が増大し、ＦＩＲフィルタによる信号の復元が困難になり、伝送距離が数ｋｍに制限されることが問題となっている。
【０００７】
そこで、前記課題を解決するために、本発明は、マルチモード光ファイバ伝送における長距離伝送を実現する技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、光送信機の個数が１個であるときに、光受信機の個数をマルチモード光ファイバの伝搬モード数以上となるようにし、ＦＩＲフィルタを利用してマルチモード光ファイバ伝送の群遅延差を補償することとした。これにより、ＦＩＲフィルタにおいて、タップ数を削減したとしても、マルチモード光ファイバ伝送において、群遅延差を補償することができ、長距離伝送を実現することができる。
【０００９】
具体的には、本発明は、光信号をＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバに送出する１個の光送信機と、前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光ファイバ伝送システムである。
【００１０】
また、本発明は、光信号を送信する１個の光送信機と、前記光信号を受信するＭ個（Ｍは２以上の整数）の光受信機と、前記１個の光送信機が送信した前記光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバと、前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光ファイバ伝送システムである。
【００１１】
また、本発明は、光信号を送信する１個の光送信機から、送信後の前記光信号を伝搬しＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバを介して、伝搬後及び分波後の前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光受信装置である。
【００１２】
この構成によれば、光送信機の個数が１個であるときに、光受信機の個数をマルチモード光ファイバの伝搬モード数以上とするため、ＦＩＲフィルタにおいて、タップ数を削減したとしても、マルチモード光ファイバ伝送において、群遅延差を補償することができ、長距離伝送を実現することができる。
【００１３】
また、本発明は、各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光ファイバ伝送システムである。
【００１４】
また、本発明は、各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光受信装置である。
【００１５】
この構成によれば、ＦＩＲフィルタにおいて、タップ数を削減することができる。
【００１６】
また、本発明は、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光ファイバ伝送システムである。
【００１７】
また、本発明は、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光受信装置である。
【００１８】
この構成によれば、ＦＩＲフィルタにおいて、タップ数をより削減することができる。
【発明の効果】
【００１９】
本発明は、マルチモード光ファイバ伝送における長距離伝送を実現する技術を提供することができる。
【図面の簡単な説明】
【００２０】
【図１】ＦＩＲフィルタ装置を利用して、複数のモードの間の群遅延差を補償する、光ファイバ伝送システム及び光受信装置の構成を示す図である。
【図２】ＦＩＲフィルタ装置の構成を示す図である。
【図３】タップ係数を更新するための適応等化アルゴリズムを示す図である。
【図４】タップ係数を更新するためのトレーニングシンボルを示す図である。
【図５】ＦＩＲフィルタ装置を利用して、光送信機が１個であり光受信機が１個であるときに、基本モード及び高次モードの間の群遅延差０．５ｎｓを補償できるかどうかを調べた、シミュレーション結果を示す図である。
【図６】ＦＩＲフィルタ装置を利用して、光送信機が１個であり光受信機が１個であるときに、基本モード及び高次モードの間の群遅延差１ｎｓを補償できるかどうかを調べた、シミュレーション結果を示す図である。
【図７】ＦＩＲフィルタ装置を利用して、光送信機が１個、光受信機が１個、光ファイバの伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図８】光送信機が１個、光受信機が２個又は３個、光ファイバの伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図９】光送信機が１個、光受信機が２個、光ファイバの伝搬モード数が２個、タップ数が１個であるときに、それらのモードの間の群遅延差及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図１０】光送信機が１個、光受信機が２個又は３個、光ファイバの伝搬モード数が３個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図１１】光送信機が１個、光受信機が３個又は４個、光ファイバの伝搬モード数が４個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓ、３ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【発明を実施するための形態】
【００２１】
添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施の例であり、本発明は以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００２２】
（ＦＩＲフィルタ装置を利用する群遅延差補償）
ＦＩＲフィルタ装置を利用して、複数のモードの間の群遅延差を補償する、光ファイバ伝送システム及び光受信装置の構成を図１に示す。光ファイバ伝送システムは、光送信機１、光受信機２−１、２−２、・・・、２−Ｍ、光ファイバ３、分波器４及びＦＩＲフィルタ装置５から構成される。光受信装置Ｒは、光受信機２−１、２−２、・・・、２−Ｍ及びＦＩＲフィルタ装置５から構成される。
【００２３】
光送信機１は、光信号ｘ（ｎ）を送信する。ここで、光信号ｘ（ｎ）はｎ番目のシンボルとして送信される光信号である。
【００２４】
光受信機２−１、２−２、・・・、２−Ｍは、それぞれ光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）を受信する。ここで、Ｍは２以上の整数であり、光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）はｎ番目のシンボルとして受信される光信号である。
【００２５】
光ファイバ３は、光送信機１が送信した光信号を伝搬し、Ｐ個の伝搬モードを有する。ここで、Ｐは２以上の整数かつＭ以下の整数である。分波器４は、光ファイバ３が伝搬した光信号を光受信機２−１、２−２、・・・、２−Ｍに分波する。
【００２６】
ＦＩＲフィルタ装置５は、群遅延差補償装置として、光受信機２−１、２−２、・・・、２−Ｍがそれぞれ受信した光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償し、光信号ｕ（ｎ）を復元する。ここで、光信号ｕ（ｎ）は、ｎ番目のシンボルとして復元される光信号である。
【００２７】
受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）は、複数のモードの間の群遅延差を包含する可能性がある。復元信号ｕ（ｎ）は、複数のモードの間の群遅延差を補償されていることが望ましく、送信信号ｘ（ｎ）と一致していることが望ましい。
【００２８】
ＦＩＲフィルタ装置の構成を図２に示す。ＦＩＲフィルタ装置５は、ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍ及び合波器５２から構成される。
【００２９】
ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍは、それぞれ受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）を入力し、それぞれ光信号ｚ１（ｎ）、ｚ２（ｎ）、・・・、ｚＭ（ｎ）を生成する。合波器５２は、光信号ｚ１（ｎ）、ｚ２（ｎ）、・・・、ｚＭ（ｎ）を合波し、復元信号ｕ（ｎ）を生成する。
【００３０】
ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍは、それぞれ１番目からＬ番目までのタップから構成される。１番目からＬ番目までのタップは、それぞれ遅延素子（図２ではτで示す）及び振幅位相調整器（図２ではｗで示す）から構成される。遅延素子は、遅延時間を調整し、振幅位相調整器は、振幅及び位相を調整する。図２では、直接形のＦＩＲフィルタを利用しているが、転置形のＦＩＲフィルタを利用してもよい。
【００３１】
受信信号ｙＢ（ｎ）についてのｉ番目のタップにおいて、遅延素子での遅延時間をτｉとし、振幅位相調整器でのタップ係数をｗＢ（ｉ）とする。遅延時間τｉ及びタップ係数ｗＢ（ｉ）を調整することにより、光ファイバ３中で発生するモード分散、波長分散、偏波モード分散などによる信号劣化を補償することができる。局発光源、９０°ハイブリッド、バランスレシーバ、アナログデジタルコンバータを利用することにより、受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）の電界の振幅及び位相の情報を取得することができる（例えば、非特許文献４参照）。
【００３２】
タップ係数を更新するための適応等化アルゴリズムを図３に示す。タップ係数を更新するためのトレーニングシンボルを図４に示す。送信信号ｘは、前段にトレーニングシンボルｘ（１）、・・・、ｘ（Ｎｔｒａｉｎｉｎｇ）を有し、後段にデータ部ｘ（Ｎｔｒａｉｎｉｎｇ＋１）、・・・、ｘ（Ｎａｌｌ）を有する。ここで、Ｎｔｒａｉｎｉｎｇ＞Ｌである必要がある。光受信装置Ｒは、データ部ｘ（Ｎｔｒａｉｎｉｎｇ＋１）、・・・、ｘ（Ｎａｌｌ）を未知とするが、トレーニングシンボルｘ（１）、・・・、ｘ（Ｎｔｒａｉｎｉｎｇ）は既知とするため、以下のようにタップ係数を更新することができる。
【００３３】
まず、トレーニングシンボルを受信中であり、ｎ≦Ｎｔｒａｉｎｉｎｇであるときについて説明する。ＦＩＲフィルタ装置５は、既に設定されたタップ係数を利用することにより、受信信号ｙ（ｎ）を復元信号ｕ（ｎ）に変換する。減算器７は、復元信号ｕ（ｎ）及び既知のトレーニング信号ｘ（ｎ）の差分を生成することにより、誤差信号ｅ（ｎ）を生成する。適応等化アルゴリズム６は、誤差信号ｅ（ｎ）を小さくするように、タップ係数を更新する。この処理を、すべてのトレーニングシンボルについて繰り返す。
【００３４】
次に、データ部を受信中であり、ｎ≧Ｎｔｒａｉｎｉｎｇ＋１であるときについて説明する。タップ係数の設定は、トレーニングシンボルの受信中に完了している。ＦＩＲフィルタ装置５は、既に設定されたタップ係数を利用することにより、受信信号ｙ（ｎ）を復元信号ｕ（ｎ）に変換する。この処理を、すべてのデータ部について繰り返す。
【００３５】
なお、適応等化アルゴリズムには、非特許文献２に記載のＬｅａｓｔ ｍｅａｎ ｓｑｕａｒｅ（ＬＭＳ）アルゴリズムやＲｅｃｕｒｓｉｖｅ ｌｅａｓｔ ｓｑｕａｒｅ（ＲＬＳ）アルゴリズムを利用することができる。
【００３６】
ＦＩＲフィルタ装置５を利用して、光送信機１が１個であり光受信機２が１個であるときに、基本モード及び高次モードの間の群遅延差Ｄ＝０．５ｎｓを補償できるかどうかを調べた、シミュレーション結果を図５に示す。ＦＩＲフィルタ装置５を利用して、光送信機１が１個であり光受信機２が１個であるときに、基本モード及び高次モードの間の群遅延差Ｄ＝１ｎｓを補償できるかどうかを調べた、シミュレーション結果を図６に示す。
【００３７】
図５及び図６において、送受信信号は、１０ＧｂｐｓのＢｉｎａｒｙ ｐｈａｓｅ ｓｈｉｆｔ ｋｅｙｉｎｇ（ＢＰＳＫ）変調されたＮｏｎ ｒｅｔｕｒｎ ｔｏ ｚｅｒｏ（ＮＲＺ）信号の４０９６ｂｉｔ分である。光ファイバ３の伝搬モード数は２個であり、基本モード及び高次モードがそれぞれ２：１の割合で励振される。Ｎｔｒａｉｎｉｎｇは４００とし、Ｌは２０とする。τｉは１０Ｇｂｐｓのシンボル長Ｔ＝０．１ｎｓを単位とし、τ１＝０、τ２＝０．１ｎｓ、τ３＝０．１×２ｎｓ、・・・、τＬ＝０．１×（Ｌ−１）ｎｓとする。適応等化アルゴリズムにはＲＬＳを用いる。ＲＬＳにおけるステップサイズ及び忘却係数は、それぞれ０．００１及び１とする。
【００３８】
基本モード及び高次モードの間の群遅延差がＤ＝０．５ｎｓである図５の場合について説明する。図５の左上には、送信信号ｘ（ｎ）のコンスタレーションマップを示す。図５の右上には、受信信号ｙ（ｎ）のコンスタレーションマップを示す。受信信号ｙ（ｎ）は、モード分散に起因して、送信信号ｘ（ｎ）を正しく反映していない。図５の左下には、復元信号ｕ（ｎ）のコンスタレーションマップを示す。復元信号ｕ（ｎ）は、ＦＩＲフィルタ装置５を利用することにより、送信信号ｘ（ｎ）を正しく復元している。図５の右下には、タップ番号及びタップ係数の絶対値の間の関係を示す。タップ係数の絶対値は、周期的に大きな値を取っており、遅延時間が群遅延差Ｄ＝０．５ｎｓの整数倍に相当するタップ番号のタップ係数が、復元に大きく寄与していることが分かる。
【００３９】
基本モード及び高次モードの間の群遅延差がＤ＝１ｎｓである図６の場合について説明する。図６の左上には、送信信号ｘ（ｎ）のコンスタレーションマップを示す。図６の右上には、受信信号ｙ（ｎ）のコンスタレーションマップを示す。受信信号ｙ（ｎ）は、モード分散に起因して、送信信号ｘ（ｎ）を正しく反映していない。図６の下段には、復元信号ｕ（ｎ）のコンスタレーションマップを示す。復元信号ｕ（ｎ）は、ＦＩＲフィルタ装置５を利用しているとしても、送信信号ｘ（ｎ）を正しく復元していない。復元の精度を向上させるためには、Ｄ／Ｔの数倍程度にＬを設定する必要があることが分かる。
【００４０】
ＦＩＲフィルタ装置５を利用して、光送信機１が１個、光受信機２が１個、光ファイバ３の伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図７に示す。シミュレーション条件は、タップ数以外は、図６及び図７において同様である。
【００４１】
ＥＶＭ値とは、コンスタレーション上での送信信号ｘ（ｎ）と復元信号ｕ（ｎ）との距離の絶対値を、コンスタレーション上での送信信号ｘ（ｎ）の原点からの距離の絶対値で除したものである。つまり、ＥＶＭ値が小さいならば、精度良く復元できていることになる。ここで、ＬをＤ／Ｔと比較して十分大きくすることにより、精度良く復元できる。しかし、Ｌを大きくすることにより、適応等化アルゴリズム６における計算量が増大するため、マルチモード光ファイバを用いた長距離伝送の実現が困難となっている。
【００４２】
（光受信機数及び伝搬モード数の大小関係）
一般的には、ＦＩＲフィルタ及び適応等化を用いた復元では、光受信機２の数Ｍを増やすと、復元の精度を向上させることができる。本発明では、光受信機２の数Ｍを伝搬モードの数Ｐ以上に増やすと、タップ数を削減しても復元の精度を向上させることができる。
【００４３】
光送信機１が１個、光受信機２が２個又は３個、光ファイバ３の伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図８に示す。
【００４４】
図８の左側は、光受信機２が２個である場合を示す。基本モードが光受信機２−１、２−２にそれぞれ１：１の割合で分波され、高次モードが光受信機２−１、２−２にそれぞれ３：２の割合で分波される。他の条件は、図７と同様である。タップ数がＤ／Ｔ＝１０以下であっても１であっても、復元の精度が向上している。
【００４５】
図８の右側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波される。他の条件は、図７と同様である。タップ数がＤ／Ｔ＝１０以下であっても１であっても、復元の精度が向上している。
【００４６】
光送信機１が１個、光受信機２が２個、光ファイバ３の伝搬モード数が２個、タップ数が１個であるときに、それらのモードの間の群遅延差及び復元精度の間の関係を調べた、シミュレーション結果を図９に示す。他の条件は、図８の左側と同様である。群遅延差Ｄの増加に対して復元信号のＥＶＭ値は低い値を維持しており、マルチモード光ファイバの伝送距離を拡大した場合においても、精度良く復元できることを示している。
【００４７】
光送信機１が１個、光受信機２が２個又は３個、光ファイバ３の伝搬モード数が３個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図１０に示す。ここで、基本モード及び第１高次モードの間の群遅延差がＤ１＝１ｎｓであり、基本モード及び第２高次モードの間の群遅延差がＤ２＝２ｎｓである。そして、基本モード、第１高次モード及び第２高次モードが、それぞれ１０：５：４の割合で励振される。
【００４８】
図１０の左側は、光受信機２が２個である場合を示す。基本モードが光受信機２−１、２−２にそれぞれ１：１の割合で分波され、第１高次モードが光受信機２−１、２−２にそれぞれ３：２の割合で分波され、第２高次モードが光受信機２−１、２−２にそれぞれ４：１の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数より小さく、タップ数がＤ１／Ｔ＝１０と同等でなければ、復元の精度が向上しない。
【００４９】
図１０の右側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３にそれぞれ６：２：２の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数と等しいため、タップ数がＤ１／Ｔ＝１０やＤ２／Ｔ＝２０以下であっても１であっても、復元の精度が向上している。
【００５０】
光送信機１が１個、光受信機２が３個又は４個、光ファイバ３の伝搬モード数が４個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓ、３ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図１１に示す。ここで、基本モード及び第１高次モードの間の群遅延差がＤ１＝１ｎｓであり、基本モード及び第２高次モードの間の群遅延差がＤ２＝２ｎｓであり、基本モード及び第３高次モードの間の群遅延差がＤ３＝３ｎｓである。そして、基本モード、第１高次モード、第２高次モード及び第３高次モードが、それぞれ１０：５：４：３の割合で励振される。
【００５１】
図１１の左側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３にそれぞれ６：２：２の割合で分波され、第３高次モードが光受信機２−１、２−２、２−３にそれぞれ７：２：１の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数より小さく、タップ数がＤ１／Ｔ＝１０と同等でなければ、復元の精度が向上しない。
【００５２】
図１１の右側は、光受信機２が４個である場合を示す。基本モードが光受信機２−１、２−２、２−３、２−４にそれぞれ４：２：２：２の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ５：２：１：１の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ３：２：４：１の割合で分波され、第３高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ２：２：３：３の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数と等しいため、タップ数がＤ１／Ｔ＝１０やＤ２／Ｔ＝２０やＤ３／Ｔ＝３０以下であっても１であっても、復元の精度が向上している。
【００５３】
光受信機数が伝搬モード数以上であるときに、タップ数がたとえ１であっても、復元の精度が向上する理由について、以下に説明する。
【００５４】
まず、光送信機１が１個、光受信機２が１個、伝搬モード数が２個であるときについて説明する。光受信機２における受信信号Ｒは、基本モードの光信号Ｅ１及び高次モードの光信号Ｅ２を含んでおり、受信信号Ｒは数式１のように表される。
【数１】
Ｅ１及びＥ２についての方程式は１つしか存在しないため、所望の基本モードの光信号Ｅ１を求めることは不可能である。よって、タップ数を十分に確保する必要がある。
【００５５】
次に、光送信機１が１個、光受信機２が２個、伝搬モード数が２個であるときについて説明する。基本モードが光受信機２−１、２−２にそれぞれａ：ｂの割合で分波されるとし、高次モードが光受信機２−１、２−２にそれぞれａ’：ｂ’の割合で分波されるとすると、光受信機２−１、２−２における受信信号Ｒ１、Ｒ２は数式２のようになる。
【数２】
ここで、一般に用いられる分波器４であるマルチモードカプラは、伝搬モード毎に異なる分岐比を与える。そして、Ｅ１及びＥ２についての方程式は２つ存在するため、所望の基本モードの光信号Ｅ１を求めることが可能である。よって、タップ数が１個であっても、復元の精度が向上する。なお、マルチモードカプラの分岐比が既知であるならば、数式２の逆行列も既知であるが、マルチモードカプラの分岐比が既知でなくても、適応等化アルゴリズム６を利用して数式２の逆行列を推定することができる。
【００５６】
次に、光送信機１が１個、光受信機２が３個、伝搬モード数が３個であるときについて説明する。基本モードが光受信機２−１、２−２、２−３にそれぞれａ：ｂ：ｃの割合で分波されるとし、第１高次モードが光受信機２−１、２−２、２−３にそれぞれａ’：ｂ’：ｃ’の割合で分波されるとし、第２高次モードが光受信機２−１、２−２、２−３にそれぞれａ’’：ｂ’’：ｃ’’の割合で分波されるとすると、光受信機２−１、２−２、２−３における受信信号Ｒ１、Ｒ２、Ｒ３は数式３のようになる。
【数３】
この場合も上の場合と同様に、所望の基本モードの光信号Ｅ１を求めることが可能である。よって、タップ数が１個であっても、復元の精度が向上する。
【産業上の利用可能性】
【００５７】
本発明に係る光ファイバ伝送システム及び光受信装置は、光ファイバ中で発生する非線形効果やファイバヒューズを低減しつつ、伝送容量を大容量化し伝送距離を長距離化する、光ファイバ伝送において適用することができる。
【符号の説明】
【００５８】
Ｒ：光受信装置
１：光送信機
２、２−１、２−２、２−Ｍ：光受信機
３：光ファイバ
４：分波器
５：ＦＩＲフィルタ装置
６：適応等化アルゴリズム
７：減算器
５１、５１−１、５１−２、５１−Ｍ：ＦＩＲフィルタ
５２：合波器
【技術分野】
【０００１】
本発明は、光ファイバ伝送の群遅延差を補償する技術に関する。
【背景技術】
【０００２】
光ファイバ伝送システムでは、光ファイバ中で発生する非線形効果やファイバヒューズが問題となり、伝送の大容量化が制限されている。これらの制限を緩和するためには、光ファイバに導波する光の密度を低減する必要があり、非特許文献１、２に示すように大コアファイバが検討されている。
【０００３】
しかし、曲げ損失低減、単一モード動作領域の拡大、実効断面積の拡大は互いにトレードオフの関係にあり、所定の条件下における実効断面積の拡大量には限界があるという課題があった。そこで、伝送路をマルチモード光ファイバとし、Ｆｉｎｉｔｅ ｉｎｐｕｌｓｅ ｒｅｓｐｏｎｓｅ（ＦＩＲ）フィルタと適応等化技術を用いて、光ファイバ中で発生するモード分散を補償する試みが行われている（例えば、非特許文献３参照）。
【０００４】
これにより、先に述べた大コア光ファイバで制限要因であった単一モード動作条件が不要になるため、さらなる大コア化及び大容量化が可能である。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Ｔ．Ｍａｔｓｕｉ，Ｋ．Ｎａｋａｊｉｍａ ａｎｄ Ｃ．Ｆｕｋａｉ，“Ａｐｐｌｉｃａｂｉｌｉｔｙ ｏｆ Ｐｈｏｔｏｎｉｃ Ｃｒｙｓｔａｌ Ｆｉｂｅｒ Ｗｉｔｈ Ｕｎｉｆｏｒｍ Ａｉｒ−Ｈｏｌｅ Ｓｔｒｕｃｔｕｒｅ ｔｏ Ｈｉｇｈ−Ｓｐｅｅｄ ａｎｄ Ｗｉｄｅ−Ｂａｎｄ Ｔｒａｎｓｍｉｓｓｉｏｎ Ｏｖｅｒ Ｃｏｎｖｅｎｔｉｏｎａｌ Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ Ｂａｎｄｓ”，Ｊ．Ｌｉｇｈｔｗａｖｅ Ｔｅｃｈｎｏｌ．２７，５４１０−５４１６，２００９．
【非特許文献２】Ｋ．Ｍｕｋａｓａ，Ｋ．Ｉｍａｍｕｒａ，Ｒ．Ｓｕｇｉｚａｋｉ ａｎｄ Ｔ．Ｙａｇｉ，“Ｃｏｍｐａｒｉｓｏｎｓ ｏｆ ｍｅｒｉｔｓ ｏｎ ｗｉｄｅ−ｂａｎｄ ｔｒａｎｓｍｉｓｓｉｏｎ ｓｙｓｔｅｍｓ ｂｅｔｗｅｅｎ ｕｓｉｎｇ ｅｘｔｒｅｍｅｌｙ ｉｍｐｒｏｖｅｄ ｓｏｌｉｄ ＳＭＦｓ ｗｉｔｈ Ａｅｆｆ ｏｆ １６０μｍ２ ａｎｄ ｌｏｓｓ ｏｆ ０．１７５ｄＢ／ｋｍ ａｎｄ ｕｓｉｎｇ ｌａｒｇｅ−Ａｅｆｆ ｈｏｌｅｙ ｆｉｂｅｒｓ ｅｎａｂｌｉｎｇ ｔｒａｎｓｍｉｓｓｉｏｎ ｏｖｅｒ ６００ｎｍ ｂａｎｄｗｉｄｔｈ”，ｔｈｅ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＯＦＣ２００８，ＯｔｈＲ１，Ｆｅｂ．２００８．
【非特許文献３】Ｘ．Ｚｈａｏ ａｎｄ Ｆ．Ｓ．Ｃｈｏａ，“Ｄｅｍｏｎｓｔｒａｔｉｏｎ ｏｆ １０−Ｇｂ／ｓ ｔｒａｎｓｍｉｓｓｉｏｎｓ ｏｖｅｒ ａ １．５−ｋｍ−ｌｏｎｇ ｍｕｌｔｉｍｏｄｅ ｆｉｂｅｒ ｕｓｉｎｇ ｅｑｕａｌｉｚａｔｉｏｎ ｔｅｃｈｎｉｑｕｅｓ”，ＩＥＥＥ Ｐｈｏｔｏｎｉｃｓ Ｔｅｃｈｎｏｌｏｇｙ Ｌｅｔｔｅｒｓ，ｖｏｌ．１４，ｐｐ．１１８７−１１８９，２００２．
【非特許文献４】Ｍ．Ｔａｙｌｏｒ，“Ｃｏｈｅｒｅｎｔ Ｄｅｔｅｃｔｉｏｎ ｆｏｒ Ｆｉｂｅｒ Ｏｐｔｉｃ Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ Ｕｓｉｎｇ Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ”，ｉｎ Ｏｐｔｉｃａｌ Ａｍｐｌｉｆｉｅｒｓ ａｎｄ Ｔｈｅｉｒ Ａｐｐｌｉｃａｔｉｏｎｓ／Ｃｏｈｅｒｅｎｔ Ｏｐｔｉｃａｌ Ｔｅｃｈｎｏｌｏｇｉｅｓ ａｎｄ Ａｐｐｌｉｃａｔｉｏｎｓ，Ｔｅｃｈｎｉｃａｌ Ｄｉｇｅｓｔ（ＣＤ）（Ｏｐｔｉｃａｌ Ｓｏｃｉｅｔｙ ｏｆ Ａｍｅｒｉｃａ，２００６），ｐａｐｅｒ ＣＴｈＢ１．
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、伝送距離の長距離化に伴ってモード間の群遅延差が増大し、ＦＩＲフィルタによる信号の復元が困難になり、伝送距離が数ｋｍに制限されることが問題となっている。
【０００７】
そこで、前記課題を解決するために、本発明は、マルチモード光ファイバ伝送における長距離伝送を実現する技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、光送信機の個数が１個であるときに、光受信機の個数をマルチモード光ファイバの伝搬モード数以上となるようにし、ＦＩＲフィルタを利用してマルチモード光ファイバ伝送の群遅延差を補償することとした。これにより、ＦＩＲフィルタにおいて、タップ数を削減したとしても、マルチモード光ファイバ伝送において、群遅延差を補償することができ、長距離伝送を実現することができる。
【０００９】
具体的には、本発明は、光信号をＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバに送出する１個の光送信機と、前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光ファイバ伝送システムである。
【００１０】
また、本発明は、光信号を送信する１個の光送信機と、前記光信号を受信するＭ個（Ｍは２以上の整数）の光受信機と、前記１個の光送信機が送信した前記光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバと、前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光ファイバ伝送システムである。
【００１１】
また、本発明は、光信号を送信する１個の光送信機から、送信後の前記光信号を伝搬しＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバを介して、伝搬後及び分波後の前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、を備えることを特徴とする光受信装置である。
【００１２】
この構成によれば、光送信機の個数が１個であるときに、光受信機の個数をマルチモード光ファイバの伝搬モード数以上とするため、ＦＩＲフィルタにおいて、タップ数を削減したとしても、マルチモード光ファイバ伝送において、群遅延差を補償することができ、長距離伝送を実現することができる。
【００１３】
また、本発明は、各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光ファイバ伝送システムである。
【００１４】
また、本発明は、各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光受信装置である。
【００１５】
この構成によれば、ＦＩＲフィルタにおいて、タップ数を削減することができる。
【００１６】
また、本発明は、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光ファイバ伝送システムである。
【００１７】
また、本発明は、前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする光受信装置である。
【００１８】
この構成によれば、ＦＩＲフィルタにおいて、タップ数をより削減することができる。
【発明の効果】
【００１９】
本発明は、マルチモード光ファイバ伝送における長距離伝送を実現する技術を提供することができる。
【図面の簡単な説明】
【００２０】
【図１】ＦＩＲフィルタ装置を利用して、複数のモードの間の群遅延差を補償する、光ファイバ伝送システム及び光受信装置の構成を示す図である。
【図２】ＦＩＲフィルタ装置の構成を示す図である。
【図３】タップ係数を更新するための適応等化アルゴリズムを示す図である。
【図４】タップ係数を更新するためのトレーニングシンボルを示す図である。
【図５】ＦＩＲフィルタ装置を利用して、光送信機が１個であり光受信機が１個であるときに、基本モード及び高次モードの間の群遅延差０．５ｎｓを補償できるかどうかを調べた、シミュレーション結果を示す図である。
【図６】ＦＩＲフィルタ装置を利用して、光送信機が１個であり光受信機が１個であるときに、基本モード及び高次モードの間の群遅延差１ｎｓを補償できるかどうかを調べた、シミュレーション結果を示す図である。
【図７】ＦＩＲフィルタ装置を利用して、光送信機が１個、光受信機が１個、光ファイバの伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図８】光送信機が１個、光受信機が２個又は３個、光ファイバの伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図９】光送信機が１個、光受信機が２個、光ファイバの伝搬モード数が２個、タップ数が１個であるときに、それらのモードの間の群遅延差及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図１０】光送信機が１個、光受信機が２個又は３個、光ファイバの伝搬モード数が３個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【図１１】光送信機が１個、光受信機が３個又は４個、光ファイバの伝搬モード数が４個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓ、３ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を示す図である。
【発明を実施するための形態】
【００２１】
添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施の例であり、本発明は以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００２２】
（ＦＩＲフィルタ装置を利用する群遅延差補償）
ＦＩＲフィルタ装置を利用して、複数のモードの間の群遅延差を補償する、光ファイバ伝送システム及び光受信装置の構成を図１に示す。光ファイバ伝送システムは、光送信機１、光受信機２−１、２−２、・・・、２−Ｍ、光ファイバ３、分波器４及びＦＩＲフィルタ装置５から構成される。光受信装置Ｒは、光受信機２−１、２−２、・・・、２−Ｍ及びＦＩＲフィルタ装置５から構成される。
【００２３】
光送信機１は、光信号ｘ（ｎ）を送信する。ここで、光信号ｘ（ｎ）はｎ番目のシンボルとして送信される光信号である。
【００２４】
光受信機２−１、２−２、・・・、２−Ｍは、それぞれ光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）を受信する。ここで、Ｍは２以上の整数であり、光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）はｎ番目のシンボルとして受信される光信号である。
【００２５】
光ファイバ３は、光送信機１が送信した光信号を伝搬し、Ｐ個の伝搬モードを有する。ここで、Ｐは２以上の整数かつＭ以下の整数である。分波器４は、光ファイバ３が伝搬した光信号を光受信機２−１、２−２、・・・、２−Ｍに分波する。
【００２６】
ＦＩＲフィルタ装置５は、群遅延差補償装置として、光受信機２−１、２−２、・・・、２−Ｍがそれぞれ受信した光信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償し、光信号ｕ（ｎ）を復元する。ここで、光信号ｕ（ｎ）は、ｎ番目のシンボルとして復元される光信号である。
【００２７】
受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）は、複数のモードの間の群遅延差を包含する可能性がある。復元信号ｕ（ｎ）は、複数のモードの間の群遅延差を補償されていることが望ましく、送信信号ｘ（ｎ）と一致していることが望ましい。
【００２８】
ＦＩＲフィルタ装置の構成を図２に示す。ＦＩＲフィルタ装置５は、ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍ及び合波器５２から構成される。
【００２９】
ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍは、それぞれ受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）を入力し、それぞれ光信号ｚ１（ｎ）、ｚ２（ｎ）、・・・、ｚＭ（ｎ）を生成する。合波器５２は、光信号ｚ１（ｎ）、ｚ２（ｎ）、・・・、ｚＭ（ｎ）を合波し、復元信号ｕ（ｎ）を生成する。
【００３０】
ＦＩＲフィルタ５１−１、５１−２、・・・、５１−Ｍは、それぞれ１番目からＬ番目までのタップから構成される。１番目からＬ番目までのタップは、それぞれ遅延素子（図２ではτで示す）及び振幅位相調整器（図２ではｗで示す）から構成される。遅延素子は、遅延時間を調整し、振幅位相調整器は、振幅及び位相を調整する。図２では、直接形のＦＩＲフィルタを利用しているが、転置形のＦＩＲフィルタを利用してもよい。
【００３１】
受信信号ｙＢ（ｎ）についてのｉ番目のタップにおいて、遅延素子での遅延時間をτｉとし、振幅位相調整器でのタップ係数をｗＢ（ｉ）とする。遅延時間τｉ及びタップ係数ｗＢ（ｉ）を調整することにより、光ファイバ３中で発生するモード分散、波長分散、偏波モード分散などによる信号劣化を補償することができる。局発光源、９０°ハイブリッド、バランスレシーバ、アナログデジタルコンバータを利用することにより、受信信号ｙ１（ｎ）、ｙ２（ｎ）、・・・、ｙＭ（ｎ）の電界の振幅及び位相の情報を取得することができる（例えば、非特許文献４参照）。
【００３２】
タップ係数を更新するための適応等化アルゴリズムを図３に示す。タップ係数を更新するためのトレーニングシンボルを図４に示す。送信信号ｘは、前段にトレーニングシンボルｘ（１）、・・・、ｘ（Ｎｔｒａｉｎｉｎｇ）を有し、後段にデータ部ｘ（Ｎｔｒａｉｎｉｎｇ＋１）、・・・、ｘ（Ｎａｌｌ）を有する。ここで、Ｎｔｒａｉｎｉｎｇ＞Ｌである必要がある。光受信装置Ｒは、データ部ｘ（Ｎｔｒａｉｎｉｎｇ＋１）、・・・、ｘ（Ｎａｌｌ）を未知とするが、トレーニングシンボルｘ（１）、・・・、ｘ（Ｎｔｒａｉｎｉｎｇ）は既知とするため、以下のようにタップ係数を更新することができる。
【００３３】
まず、トレーニングシンボルを受信中であり、ｎ≦Ｎｔｒａｉｎｉｎｇであるときについて説明する。ＦＩＲフィルタ装置５は、既に設定されたタップ係数を利用することにより、受信信号ｙ（ｎ）を復元信号ｕ（ｎ）に変換する。減算器７は、復元信号ｕ（ｎ）及び既知のトレーニング信号ｘ（ｎ）の差分を生成することにより、誤差信号ｅ（ｎ）を生成する。適応等化アルゴリズム６は、誤差信号ｅ（ｎ）を小さくするように、タップ係数を更新する。この処理を、すべてのトレーニングシンボルについて繰り返す。
【００３４】
次に、データ部を受信中であり、ｎ≧Ｎｔｒａｉｎｉｎｇ＋１であるときについて説明する。タップ係数の設定は、トレーニングシンボルの受信中に完了している。ＦＩＲフィルタ装置５は、既に設定されたタップ係数を利用することにより、受信信号ｙ（ｎ）を復元信号ｕ（ｎ）に変換する。この処理を、すべてのデータ部について繰り返す。
【００３５】
なお、適応等化アルゴリズムには、非特許文献２に記載のＬｅａｓｔ ｍｅａｎ ｓｑｕａｒｅ（ＬＭＳ）アルゴリズムやＲｅｃｕｒｓｉｖｅ ｌｅａｓｔ ｓｑｕａｒｅ（ＲＬＳ）アルゴリズムを利用することができる。
【００３６】
ＦＩＲフィルタ装置５を利用して、光送信機１が１個であり光受信機２が１個であるときに、基本モード及び高次モードの間の群遅延差Ｄ＝０．５ｎｓを補償できるかどうかを調べた、シミュレーション結果を図５に示す。ＦＩＲフィルタ装置５を利用して、光送信機１が１個であり光受信機２が１個であるときに、基本モード及び高次モードの間の群遅延差Ｄ＝１ｎｓを補償できるかどうかを調べた、シミュレーション結果を図６に示す。
【００３７】
図５及び図６において、送受信信号は、１０ＧｂｐｓのＢｉｎａｒｙ ｐｈａｓｅ ｓｈｉｆｔ ｋｅｙｉｎｇ（ＢＰＳＫ）変調されたＮｏｎ ｒｅｔｕｒｎ ｔｏ ｚｅｒｏ（ＮＲＺ）信号の４０９６ｂｉｔ分である。光ファイバ３の伝搬モード数は２個であり、基本モード及び高次モードがそれぞれ２：１の割合で励振される。Ｎｔｒａｉｎｉｎｇは４００とし、Ｌは２０とする。τｉは１０Ｇｂｐｓのシンボル長Ｔ＝０．１ｎｓを単位とし、τ１＝０、τ２＝０．１ｎｓ、τ３＝０．１×２ｎｓ、・・・、τＬ＝０．１×（Ｌ−１）ｎｓとする。適応等化アルゴリズムにはＲＬＳを用いる。ＲＬＳにおけるステップサイズ及び忘却係数は、それぞれ０．００１及び１とする。
【００３８】
基本モード及び高次モードの間の群遅延差がＤ＝０．５ｎｓである図５の場合について説明する。図５の左上には、送信信号ｘ（ｎ）のコンスタレーションマップを示す。図５の右上には、受信信号ｙ（ｎ）のコンスタレーションマップを示す。受信信号ｙ（ｎ）は、モード分散に起因して、送信信号ｘ（ｎ）を正しく反映していない。図５の左下には、復元信号ｕ（ｎ）のコンスタレーションマップを示す。復元信号ｕ（ｎ）は、ＦＩＲフィルタ装置５を利用することにより、送信信号ｘ（ｎ）を正しく復元している。図５の右下には、タップ番号及びタップ係数の絶対値の間の関係を示す。タップ係数の絶対値は、周期的に大きな値を取っており、遅延時間が群遅延差Ｄ＝０．５ｎｓの整数倍に相当するタップ番号のタップ係数が、復元に大きく寄与していることが分かる。
【００３９】
基本モード及び高次モードの間の群遅延差がＤ＝１ｎｓである図６の場合について説明する。図６の左上には、送信信号ｘ（ｎ）のコンスタレーションマップを示す。図６の右上には、受信信号ｙ（ｎ）のコンスタレーションマップを示す。受信信号ｙ（ｎ）は、モード分散に起因して、送信信号ｘ（ｎ）を正しく反映していない。図６の下段には、復元信号ｕ（ｎ）のコンスタレーションマップを示す。復元信号ｕ（ｎ）は、ＦＩＲフィルタ装置５を利用しているとしても、送信信号ｘ（ｎ）を正しく復元していない。復元の精度を向上させるためには、Ｄ／Ｔの数倍程度にＬを設定する必要があることが分かる。
【００４０】
ＦＩＲフィルタ装置５を利用して、光送信機１が１個、光受信機２が１個、光ファイバ３の伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図７に示す。シミュレーション条件は、タップ数以外は、図６及び図７において同様である。
【００４１】
ＥＶＭ値とは、コンスタレーション上での送信信号ｘ（ｎ）と復元信号ｕ（ｎ）との距離の絶対値を、コンスタレーション上での送信信号ｘ（ｎ）の原点からの距離の絶対値で除したものである。つまり、ＥＶＭ値が小さいならば、精度良く復元できていることになる。ここで、ＬをＤ／Ｔと比較して十分大きくすることにより、精度良く復元できる。しかし、Ｌを大きくすることにより、適応等化アルゴリズム６における計算量が増大するため、マルチモード光ファイバを用いた長距離伝送の実現が困難となっている。
【００４２】
（光受信機数及び伝搬モード数の大小関係）
一般的には、ＦＩＲフィルタ及び適応等化を用いた復元では、光受信機２の数Ｍを増やすと、復元の精度を向上させることができる。本発明では、光受信機２の数Ｍを伝搬モードの数Ｐ以上に増やすと、タップ数を削減しても復元の精度を向上させることができる。
【００４３】
光送信機１が１個、光受信機２が２個又は３個、光ファイバ３の伝搬モード数が２個、それらのモードの間の群遅延差が１ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図８に示す。
【００４４】
図８の左側は、光受信機２が２個である場合を示す。基本モードが光受信機２−１、２−２にそれぞれ１：１の割合で分波され、高次モードが光受信機２−１、２−２にそれぞれ３：２の割合で分波される。他の条件は、図７と同様である。タップ数がＤ／Ｔ＝１０以下であっても１であっても、復元の精度が向上している。
【００４５】
図８の右側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波される。他の条件は、図７と同様である。タップ数がＤ／Ｔ＝１０以下であっても１であっても、復元の精度が向上している。
【００４６】
光送信機１が１個、光受信機２が２個、光ファイバ３の伝搬モード数が２個、タップ数が１個であるときに、それらのモードの間の群遅延差及び復元精度の間の関係を調べた、シミュレーション結果を図９に示す。他の条件は、図８の左側と同様である。群遅延差Ｄの増加に対して復元信号のＥＶＭ値は低い値を維持しており、マルチモード光ファイバの伝送距離を拡大した場合においても、精度良く復元できることを示している。
【００４７】
光送信機１が１個、光受信機２が２個又は３個、光ファイバ３の伝搬モード数が３個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図１０に示す。ここで、基本モード及び第１高次モードの間の群遅延差がＤ１＝１ｎｓであり、基本モード及び第２高次モードの間の群遅延差がＤ２＝２ｎｓである。そして、基本モード、第１高次モード及び第２高次モードが、それぞれ１０：５：４の割合で励振される。
【００４８】
図１０の左側は、光受信機２が２個である場合を示す。基本モードが光受信機２−１、２−２にそれぞれ１：１の割合で分波され、第１高次モードが光受信機２−１、２−２にそれぞれ３：２の割合で分波され、第２高次モードが光受信機２−１、２−２にそれぞれ４：１の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数より小さく、タップ数がＤ１／Ｔ＝１０と同等でなければ、復元の精度が向上しない。
【００４９】
図１０の右側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３にそれぞれ６：２：２の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数と等しいため、タップ数がＤ１／Ｔ＝１０やＤ２／Ｔ＝２０以下であっても１であっても、復元の精度が向上している。
【００５０】
光送信機１が１個、光受信機２が３個又は４個、光ファイバ３の伝搬モード数が４個、それらのモードの間の群遅延差が１ｎｓ、２ｎｓ、３ｎｓであるときに、タップ数及び復元精度の間の関係を調べた、シミュレーション結果を図１１に示す。ここで、基本モード及び第１高次モードの間の群遅延差がＤ１＝１ｎｓであり、基本モード及び第２高次モードの間の群遅延差がＤ２＝２ｎｓであり、基本モード及び第３高次モードの間の群遅延差がＤ３＝３ｎｓである。そして、基本モード、第１高次モード、第２高次モード及び第３高次モードが、それぞれ１０：５：４：３の割合で励振される。
【００５１】
図１１の左側は、光受信機２が３個である場合を示す。基本モードが光受信機２−１、２−２、２−３にそれぞれ４：３：３の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３にそれぞれ５：３：２の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３にそれぞれ６：２：２の割合で分波され、第３高次モードが光受信機２−１、２−２、２−３にそれぞれ７：２：１の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数より小さく、タップ数がＤ１／Ｔ＝１０と同等でなければ、復元の精度が向上しない。
【００５２】
図１１の右側は、光受信機２が４個である場合を示す。基本モードが光受信機２−１、２−２、２−３、２−４にそれぞれ４：２：２：２の割合で分波され、第１高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ５：２：１：１の割合で分波され、第２高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ３：２：４：１の割合で分波され、第３高次モードが光受信機２−１、２−２、２−３、２−４にそれぞれ２：２：３：３の割合で分波される。他の条件は、図９と同様である。光受信機数が伝搬モード数と等しいため、タップ数がＤ１／Ｔ＝１０やＤ２／Ｔ＝２０やＤ３／Ｔ＝３０以下であっても１であっても、復元の精度が向上している。
【００５３】
光受信機数が伝搬モード数以上であるときに、タップ数がたとえ１であっても、復元の精度が向上する理由について、以下に説明する。
【００５４】
まず、光送信機１が１個、光受信機２が１個、伝搬モード数が２個であるときについて説明する。光受信機２における受信信号Ｒは、基本モードの光信号Ｅ１及び高次モードの光信号Ｅ２を含んでおり、受信信号Ｒは数式１のように表される。
【数１】
Ｅ１及びＥ２についての方程式は１つしか存在しないため、所望の基本モードの光信号Ｅ１を求めることは不可能である。よって、タップ数を十分に確保する必要がある。
【００５５】
次に、光送信機１が１個、光受信機２が２個、伝搬モード数が２個であるときについて説明する。基本モードが光受信機２−１、２−２にそれぞれａ：ｂの割合で分波されるとし、高次モードが光受信機２−１、２−２にそれぞれａ’：ｂ’の割合で分波されるとすると、光受信機２−１、２−２における受信信号Ｒ１、Ｒ２は数式２のようになる。
【数２】
ここで、一般に用いられる分波器４であるマルチモードカプラは、伝搬モード毎に異なる分岐比を与える。そして、Ｅ１及びＥ２についての方程式は２つ存在するため、所望の基本モードの光信号Ｅ１を求めることが可能である。よって、タップ数が１個であっても、復元の精度が向上する。なお、マルチモードカプラの分岐比が既知であるならば、数式２の逆行列も既知であるが、マルチモードカプラの分岐比が既知でなくても、適応等化アルゴリズム６を利用して数式２の逆行列を推定することができる。
【００５６】
次に、光送信機１が１個、光受信機２が３個、伝搬モード数が３個であるときについて説明する。基本モードが光受信機２−１、２−２、２−３にそれぞれａ：ｂ：ｃの割合で分波されるとし、第１高次モードが光受信機２−１、２−２、２−３にそれぞれａ’：ｂ’：ｃ’の割合で分波されるとし、第２高次モードが光受信機２−１、２−２、２−３にそれぞれａ’’：ｂ’’：ｃ’’の割合で分波されるとすると、光受信機２−１、２−２、２−３における受信信号Ｒ１、Ｒ２、Ｒ３は数式３のようになる。
【数３】
この場合も上の場合と同様に、所望の基本モードの光信号Ｅ１を求めることが可能である。よって、タップ数が１個であっても、復元の精度が向上する。
【産業上の利用可能性】
【００５７】
本発明に係る光ファイバ伝送システム及び光受信装置は、光ファイバ中で発生する非線形効果やファイバヒューズを低減しつつ、伝送容量を大容量化し伝送距離を長距離化する、光ファイバ伝送において適用することができる。
【符号の説明】
【００５８】
Ｒ：光受信装置
１：光送信機
２、２−１、２−２、２−Ｍ：光受信機
３：光ファイバ
４：分波器
５：ＦＩＲフィルタ装置
６：適応等化アルゴリズム
７：減算器
５１、５１−１、５１−２、５１−Ｍ：ＦＩＲフィルタ
５２：合波器
【特許請求の範囲】
【請求項１】
光信号をＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバに送出する１個の光送信機と、
前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、
前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光ファイバ伝送システム。
【請求項２】
光信号を送信する１個の光送信機と、
前記光信号を受信するＭ個（Ｍは２以上の整数）の光受信機と、
前記１個の光送信機が送信した前記光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバと、
前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光ファイバ伝送システム。
【請求項３】
各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項１又は請求項２に記載の光ファイバ伝送システム。
【請求項４】
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項１から請求項３のいずれかに記載の光ファイバ伝送システム。
【請求項５】
光信号を送信する１個の光送信機から、送信後の前記光信号を伝搬しＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバを介して、伝搬後及び分波後の前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光受信装置。
【請求項６】
各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項５に記載の光受信装置。
【請求項７】
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項５又は請求項６に記載の光受信装置。
【請求項１】
光信号をＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバに送出する１個の光送信機と、
前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、
前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光ファイバ伝送システム。
【請求項２】
光信号を送信する１個の光送信機と、
前記光信号を受信するＭ個（Ｍは２以上の整数）の光受信機と、
前記１個の光送信機が送信した前記光信号を伝搬し、Ｐ個（Ｐは２以上の整数かつＭ以下の整数）の伝搬モードを有する光ファイバと、
前記光ファイバが伝搬した前記光信号を前記Ｍ個の光受信機に分波する分波器と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光ファイバ伝送システム。
【請求項３】
各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項１又は請求項２に記載の光ファイバ伝送システム。
【請求項４】
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項１から請求項３のいずれかに記載の光ファイバ伝送システム。
【請求項５】
光信号を送信する１個の光送信機から、送信後の前記光信号を伝搬しＰ個（Ｐは２以上の整数）の伝搬モードを有する光ファイバを介して、伝搬後及び分波後の前記光信号を受信するＭ個（Ｍは２以上の整数かつＰ以上の整数）の光受信機と、
各光受信機が受信した前記光信号について、遅延時間を調整する遅延素子並びに振幅及び位相を調整する振幅位相調整器からなるＦＩＲフィルタを利用して、群遅延差を補償する群遅延差補償装置と、
を備えることを特徴とする光受信装置。
【請求項６】
各光受信機が受信した前記光信号の群遅延差をＤとし、前記１個の光送信機が送信した前記光信号の１シンボル時間をＴとするとき、
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、Ｄ／Ｔ個以下の前記遅延素子及びＤ／Ｔ個以下の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項５に記載の光受信装置。
【請求項７】
前記群遅延差補償装置は、各光受信機が受信した前記光信号について、１個の前記遅延素子及び１個の前記振幅位相調整器を利用して、群遅延差を補償することを特徴とする、請求項５又は請求項６に記載の光受信装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１２−２２７７６４（Ｐ２０１２−２２７７６４Ａ）
【公開日】平成２４年１１月１５日（２０１２．１１．１５）
【国際特許分類】
【出願番号】特願２０１１−９４０５２（Ｐ２０１１−９４０５２）
【出願日】平成２３年４月２０日（２０１１．４．２０）
【出願人】（０００００４２２６）日本電信電話株式会社 (13,992)
【Ｆターム（参考）】
【公開日】平成２４年１１月１５日（２０１２．１１．１５）
【国際特許分類】
【出願日】平成２３年４月２０日（２０１１．４．２０）
【出願人】（０００００４２２６）日本電信電話株式会社 (13,992)
【Ｆターム（参考）】
[ Back to top ]