(57)【特許請求の範囲】
【請求項１】
入力光信号の時間波形を再構成する波形再構成装置であって、
前記入力光信号のパワースペクトルを取得する入力スペクトル取得部と、
複数の強度の前記入力光信号について、非線形光学効果を有する光伝送媒体内を前記入力光信号が伝播した後に出力される出力光信号を計測することにより得られるパワースペクトルを計測パワースペクトルとしてそれぞれ取得する出力スペクトル取得部と、
（ａ）前記入力光信号の強度ごとに、前記入力光信号が所与の位相スペクトルを有すると仮定して前記入力光信号の前記光伝送媒体内の伝播シミュレーションを行うことにより、前記出力光信号の計算パワースペクトルを算出し、（ｂ）前記計算パワースペクトルと前記計測パワースペクトルとの前記入力光信号の強度ごとの差分が小さくなるように前記所与の位相スペクトルを変化させて前記伝播シミュレーションを繰り返し行うことにより、前記入力光信号の位相スペクトルを探索する位相スペクトル算出部と、
前記位相スペクトル算出部により探索された位相スペクトルと、前記入力スペクトル取得部により取得されたパワースペクトルとを周波数／時間変換することにより、前記入力光信号の時間波形を再構成する波形再構成部とを備え、
前記位相スペクトル算出部は、前記伝播シミュレーションの繰り返しの少なくとも一部において、前記光伝送媒体の非線形定数γに比例する非線形光学効果と、前記光伝送媒体の２次分散β2以上の高次分散による分散効果とのうちの非線形光学効果のみに基づいて前記伝播シミュレーションを行う
波形再構成装置。
【請求項２】
前記位相スペクトル算出部は、さらに、前記伝播シミュレーションの繰り返しの少なくとも一部において、前記入力光信号が非線形光学効果及び分散効果のうちの分散効果のみを有する媒体内を伝播して得られた信号であると仮定した場合に満たされるべき位相スペクトルの制約条件に従って、前記所与の位相スペクトルを変化させる
請求項１に記載の波形再構成装置。
【請求項３】
前記位相スペクトル算出部は、前記入力光信号の位相スペクトルを探索する際に、前記制約条件に従って前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行った後に、前記制約条件に関わらず前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行う
請求項２に記載の波形再構成装置。
【請求項４】
前記位相スペクトル算出部は、前記非線形光学効果に関するパラメータのうち、自己位相変調に関するパラメータのみを用いて、前記伝播シミュレーションを行う
請求項１に記載の波形再構成装置。
【請求項５】
前記光伝送媒体は、光ファイバである
請求項１〜４のいずれか１項に記載の波形再構成装置。
【請求項６】
入力光信号の時間波形を再構成する波形再構成システムであって、
前記入力光信号の強度を互いに異なる複数の強度に変化させる強度調節器と、
非線形光学効果を有する光伝送媒体であって、少なくとも非線形光学効果に関するパラメータが既知である光伝送媒体と、
前記強度調節器により強度が変化された前記入力光信号が前記光伝送媒体内を伝播した後に出力される出力光信号のパワースペクトルを前記入力光信号の強度ごとに計測するスペクトル計測器と、
請求項１に記載の波形再構成装置とを備え、
前記出力スペクトル取得部は、前記スペクトル計測器により計測されたパワースペクトルを計測パワースペクトルとして取得する
波形再構成システム。
【請求項７】
前記光伝送媒体は、光ファイバである
請求項６に記載の波形再構成システム。
【請求項８】
入力光信号の時間波形を再構成する波形再構成方法であって、
前記入力光信号のパワースペクトルを取得するステップと、
複数の強度の前記入力光信号について、非線形光学効果を有する光伝送媒体内を前記入力光信号が伝播した後に出力される出力光信号を計測することにより得られるパワースペクトルを計測パワースペクトルとしてそれぞれ取得するステップと、
前記入力光信号の強度ごとに、前記入力光信号が所与の位相スペクトルを有すると仮定して前記入力光信号の前記光伝送媒体内の伝播シミュレーションを繰り返し行うことにより、前記出力光信号の計算パワースペクトルを算出するステップと、
前記計算パワースペクトルと前記計測パワースペクトルとの前記入力光信号の強度ごとの差分が小さくなるように前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行うことにより、前記入力光信号の位相スペクトルを探索するステップと、
前記入力光信号の位相スペクトルを探索するステップにおいて探索された位相スペクトルと、前記入力光信号のパワースペクトルを取得するステップにおいて取得されたパワースペクトルとを周波数／時間変換することにより、前記入力光信号の時間波形を再構成するステップとを含み、
前記入力光信号の位相スペクトルを探索するステップでは、前記伝播シミュレーションの繰り返しの少なくとも一部において、前記光伝送媒体の非線形定数γに比例する非線形光学効果と、前記光伝送媒体の２次分散β2以上の高次分散による分散効果とのうちの非線形光学効果のみに基づいて前記伝播シミュレーションを行う
波形再構成方法。
【請求項９】
請求項８に記載の波形再構成方法をコンピュータに実行させるためのプログラム。
【請求項１０】
入力光信号の時間波形を再構成する波形再構成装置であって、
前記入力光信号のパワースペクトルを示す情報を取得する入力スペクトル取得部と、
複数の強度の前記入力光信号について、非線形光学効果を有する光伝送媒体内を前記入力光信号が伝播した後に出力される出力光信号を計測することにより得られるパワースペクトルを計測パワースペクトルとしてそれぞれ取得する出力スペクトル取得部と、
（ａ）前記入力光信号の強度ごとに、前記入力光信号が所与の位相スペクトルを有すると仮定して前記入力光信号の前記光伝送媒体内の伝播シミュレーションを行うことにより、前記出力光信号の計算パワースペクトルを算出し、（ｂ）前記計算パワースペクトルと前記計測パワースペクトルとの前記入力光信号の強度ごとの差異が小さくなるように前記所与の位相スペクトルを変化させて前記伝播シミュレーションを繰り返し行うことにより、前記入力光信号の位相スペクトルを探索する位相スペクトル算出部と、
前記位相スペクトル算出部により探索された位相スペクトルと、前記入力スペクトル取得部により取得された情報が示すパワースペクトルとを周波数／時間変換することにより、前記入力光信号の時間波形を再構成する波形再構成部とを備え、
前記位相スペクトル算出部は、前記伝播シミュレーションの繰り返しにおいて、非線形光学効果及び分散効果のうちの分散効果のみを有する媒体内を伝播して得られた信号であると仮定した場合に満たされるべき位相スペクトルの制約条件に従って前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行った後に、前記制約条件に関わらず前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行う
波形再構成装置。
【請求項１１】
入力光信号の時間波形を再構成する波形再構成方法であって、
前記入力光信号のパワースペクトルを示す情報を取得するステップと、
複数の強度の前記入力光信号について、非線形光学効果を有する光伝送媒体内を前記入力光信号が伝播した後に出力される出力光信号を計測することにより得られるパワースペクトルを計測パワースペクトルとしてそれぞれ取得するステップと、
前記入力光信号の強度ごとに、前記入力光信号が所与の位相スペクトルを有すると仮定して前記入力光信号の前記光伝送媒体内の伝播シミュレーションを行うことにより、前記出力光信号の計算パワースペクトルを算出するステップと、
前記計算パワースペクトルと前記計測パワースペクトルとの前記入力光信号の強度ごとの差異が小さくなるように前記所与の位相スペクトルを変化させて前記伝播シミュレーションを繰り返し行うことにより、前記入力光信号の位相スペクトルを探索するステップと、
探索された前記位相スペクトルと、取得された前記情報が示すパワースペクトルとを周波数／時間変換することにより、前記入力光信号の時間波形を再構成するステップとを含み、
前記入力光信号の位相スペクトルを探索するステップでは、前記伝播シミュレーションの繰り返しにおいて、前記伝播シミュレーションの繰り返しにおいて、非線形光学効果及び分散効果のうちの分散効果のみを有する媒体内を伝播して得られた信号であると仮定した場合に満たされるべき位相スペクトルの制約条件に従って前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行った後に、前記制約条件に関わらず前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行う
波形再構成方法。
【請求項１２】
請求項１１に記載の波形再構成方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光信号の時間波形を再構成する波形再構成装置等に関するものである。
【背景技術】
【０００２】
近年、非線形光学効果を利用した情報通信システム等の実用化に向け、光信号の正確な時間波形の情報を得ることは極めて重要となっている。そこで、光信号の正確な時間波形の情報を取得するために、光サンプリングオシロスコープ、オートコリレータなどを用いて、光信号の時間波形の強度分布が計測される。その結果、計測された強度分布から、光信号の振幅情報が取得される。しかしながら、計測された強度分布からは、光信号の位相情報は取得できない。つまり、光信号の強度分布を計測しただけでは、光信号の時間波形の情報を取得することはできない。
【０００３】
そこで、光信号の位相を取得するための様々な方法が、提案されている（例えば、非特許文献１及び２を参照）。非特許文献１及び２に記載の方法は、時間分解分光に基づいて光信号の位相を取得する方法である。具体的には、非特許文献１及び２に記載の方法では、超高速時間ゲート又は参照光源を利用することにより、光信号の位相を取得する。そして、取得した位相を用いて、光信号の時間波形が再構成される。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Ｄ．Ｊ．Ｋａｎｅ、Ｒ．Ｔｒｅｂｉｎｏ、「Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ  ｏｆ  ａｒｂｉｔｒａｒｙ  ｆｅｍｔｏｓｅｃｏｎｄ  ｐｕｌｓｅｓ  ｕｓｉｎｇ  ｆｒｅｑｕｅｎｃｙ−ｒｅｓｏｌｖｅｄ  ｏｐｔｉｃａｌ  ｇａｔｉｎｇ」、ＩＥＥＥ  Ｊ．Ｑｕａｎｔｕｍ  Ｅｌｅｃｔｒｏｎ、Ｖｏｌ．２９、１９９３、ｐｐ５７１〜ｐｐ５７９
【非特許文献２】Ｃ．Ｄｏｒｒｅｒ、Ｍ．Ｊｏｆｆｒｅ、「Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ  ｏｆ  ｔｈｅ  ｓｐｅｃｔｒａｌ  ｐｈａｓｅ  ｏｆ  ｕｌｔｒａｓｈｏｒｔ  ｌｉｇｈｔ  ｐｕｌｓｅｓ」、Ｃ．Ｒ．Ａｃａｄ．Ｓｃｉ．Ｐａｒｉｓ、Ｖｏｌ．２、２００１、ｐｐ１４１５
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記従来の方法では、超高速時間ゲート又は参照光源を必要とするために、光信号とゲートとの時間調整並びに安定性及びＳＮ（Ｓｉｇｎａｌ−Ｎｏｉｓｅ）比の確保など非常に高い技術レベルが要求される。
【０００６】
そこで、本発明は、光信号の時間波形を簡便に再構成できる波形再構成装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記目的を達成するために、本発明の一態様に係る波形再構成装置は、入力光信号の時間波形を再構成する波形再構成装置であって、前記入力光信号のパワースペクトルを取得する入力スペクトル取得部と、複数の強度の前記入力光信号について、非線形光学効果を有する光伝送媒体内を前記入力光信号が伝播した後に出力される出力光信号を計測することにより得られるパワースペクトルを計測パワースペクトルとしてそれぞれ取得する出力スペクトル取得部と、（ａ）前記入力光信号の強度ごとに、前記入力光信号が所与の位相スペクトルを有すると仮定して前記入力光信号の前記光伝送媒体内の伝播シミュレーションを行うことにより、前記出力光信号の計算パワースペクトルを算出し、（ｂ）前記計算パワースペクトルと前記計測パワースペクトルとの前記入力光信号の強度ごとの差分が小さくなるように前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行うことにより、前記入力光信号の位相スペクトルを探索する位相スペクトル算出部と、前記位相スペクトル算出部により探索された位相スペクトルと、前記入力スペクトル取得部により取得されたパワースペクトルとを周波数／時間変換することにより、前記入力光信号の時間波形を再構成する波形再構成部とを備え、前記位相スペクトル算出部は、非線形光学効果及び分散効果のうちの一方のみに基づいて、前記所与の位相スペクトルを変化させる、又は前記伝播シミュレーションを行う。
【０００８】
この構成によれば、非線形光学効果を有する光伝送媒体内を入力光信号が伝播した後に出力される出力光信号のパワースペクトルを用いて、入力光信号の時間波形を再構成することができる。したがって、パワースペクトルを計測できれば入力光信号の時間波形を再構成できるので、光信号の時間波形を簡便に再構成することが可能となる。さらに、非線形光学効果及び分散効果のうちの一方のみに基づいて、所与の位相スペクトルを変化させる、又は伝播シミュレーションを行うことできる。したがって、非線形光学効果及び分散効果の両方に基づいて又は基づかずに、伝播シミュレーションを行う、又は所与の位相スペクトルを変化させる場合よりも、位相スペクトルの探索時間を短縮することが可能となる。
【０００９】
また、前記位相スペクトル算出部は、前記入力光信号が非線形光学効果及び分散効果のうちの分散効果のみを有する媒体内を伝播して得られた信号であると仮定した場合に満たされるべき位相スペクトルの制約条件に従って、前記所与の位相スペクトルを変化させることが好ましい。
【００１０】
この構成によれば、適切な制約条件に従って所与の位相スペクトルを変化させることができる。したがって、位相スペクトルの探索において、変数の数を減少させることができので、探索時間を短縮することが可能となる。
【００１１】
また、前記位相スペクトル算出部は、前記入力光信号の位相スペクトルを探索する際に、前記制約条件に従って前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行った後に、前記制約条件に関わらず前記所与の位相スペクトルを変化させて前記伝播シミュレーションを行うことが好ましい。
【００１２】
この構成によれば、制約条件に従って所与の位相スペクトル変化させた後に、制約条件に関わらず所与の位相スペクトルを変化させることができる。したがって、位相スペクトルの探索において、探索精度の低下を抑制しつつ、探索時間を短縮することが可能となる。
【００１３】
また、前記位相スペクトル算出部は、非線形光学効果と分散効果とのうち非線形光学効果に関するパラメータのみを用いて前記伝播シミュレーションを行うことが好ましい。
【００１４】
この構成によれば、非線形光学効果と分散効果とのうち非線形光学効果に関するパラメータのみを用いて伝播シミュレーションを行うことができる。したがって、伝播シミュレーションが単純化されるので、計算速度を向上させることが可能となる。
【００１５】
また、前記位相スペクトル算出部は、前記非線形光学効果に関するパラメータのうち、自己位相変調に関するパラメータのみを用いて、前記伝播シミュレーションを行っててもよい。
【００１６】
この構成によれば、非線形光学効果に関するパラメータのうち、自己位相変調に関するパラメータのみを用いて、伝播シミュレーションを行うことができる。したがって、さらに伝播シミュレーションが単純化されるので、さらに計算速度を向上させることが可能となる。
【００１７】
また、前記光伝送媒体は、光ファイバであることが好ましい。
【００１８】
この構成によれば、非線形光学効果を有する光伝送媒体を簡易な構成とすることができる。
【００１９】
また、本発明の一態様に係る波形再構成システムは、入力光信号の時間波形を再構成する波形再構成システムであって、前記入力光信号の強度を互いに異なる複数の強度に変化させる強度調節器と、非線形光学効果を有する光伝送媒体であって、少なくとも非線形光学効果に関するパラメータが既知である光伝送媒体と、前記強度調節器により強度が変化された前記入力光信号が前記光伝送媒体内を伝播した後に出力される出力光信号のパワースペクトルを前記入力光信号の強度ごとに計測するスペクトル計測器と、上記波形再構成装置とを備え、前記出力スペクトル取得部は、前記スペクトル計測器により計測されたパワースペクトルを計測パワースペクトルとして取得する。
【００２０】
この構成によれば、非線形光学効果に関するパラメータが既知の光伝送媒体とスペクトル計測器とがあれば入力光信号の時間波形を再構成でき、光信号の時間波形を簡便に再構成することが可能となる。さらに、非線形光学効果及び分散効果のうちの一方のみに基づいて、所与の位相スペクトルを変化させる、又は伝播シミュレーションを行うことできる。したがって、非線形光学効果及び分散効果の両方に基づいて又は基づかずに、伝播シミュレーションを行う、又は所与の位相スペクトルを変化させる場合よりも、位相スペクトルの探索時間を短縮することが可能となる。
【００２１】
また、前記光伝送媒体は、光ファイバであることが好ましい。
【００２２】
この構成によれば、非線形光学効果を有する光伝送媒体を簡易な構成とすることができる。
【００２３】
なお、本発明は、このような波形再構成装置として実現することができるだけでなく、このような波形再構成装置が備える特徴的な構成部の動作をステップとする波形再構成方法として実現することができる。さらに、波形再構成方法に含まれる各ステップをコンピュータに実行させるプログラムとして実現することもできる。そして、そのようなプログラムは、ＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ−Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）等の非一時的な記録媒体又はインターネット等の伝送媒体を介して配信することができるのは言うまでもない。
【発明の効果】
【００２４】
以上の説明から明らかなように、本発明によれば、光信号の時間波形を簡便に再構成することが可能となる。
【図面の簡単な説明】
【００２５】
【図１】図１は、本発明の実施の形態１に係る波形再構成システムの全体的な構成を示す図である。
【図２】図２は、本発明の実施の形態１に係る波形再構成装置の特徴的な機能構成を示すブロック図である。
【図３】図３は、本発明の実施の形態１に係る波形再構成装置の処理動作を示すフローチャートである。
【図４】図４は、本発明の実施の形態１におけるスプリットステップフーリエ法を説明するための図である。
【図５】図５は、本発明の実施の形態１において分光器により計測された出力光信号のパワースペクトルの実験結果を示すグラフである。
【図６】図６は、本発明の実施の形態１において位相スペクトルの探索中における評価関数の変化の実験結果を示すグラフである。
【図７】図７は、本発明の実施の形態１において実験により得られたパワースペクトルと位相スペクトルとを示す図である。
【図８】図８は、本発明の実施の形態１において波形再構成装置により再構成された入力光信号の時間波形の実験結果を示す図である。
【図９】図９は、本発明の実施の形態１において７種類の入力光信号の強度に対応する出力光信号のパワースペクトルの実測値と計算値との比較結果を示すグラフである。
【図１０】図１０は、本発明の実施の形態１において分光器により計測された１０種類の入力光信号の強度に対応する出力光信号のパワースペクトルのうち、位相スペクトルの算出に利用されなかった３種類の入力光信号の強度に対応する出力光信号のパワースペクトルの実測値と計算値との比較結果を示すグラフである。
【図１１】図１１は、本発明の実施の形態２の変形例の実験で用いた入力光信号のパワースペクトルを示すグラフである。
【図１２】図１２は、本発明の実施の形態２の変形例の実験において、２つの強度の入力光信号が高非線形光ファイバを伝播した後に分光器により計測されたパワースペクトルを示すグラフである。
【図１３】図１３は、本発明の実施の形態２の変形例に係る波形再構成装置による実験結果を示す図である。
【図１４】図１４は、本発明の実施の形態２の変形例に係る波形再構成装置による実験結果を示す図である。
【図１５】図１５は、本発明の実施の形態３の変形例に係る波形再構成装置による実験結果を示す図である。
【図１６】図１６は、コンピュータのハードウェア構成の一例を示す図である。
【発明を実施するための形態】
【００２６】
本発明の一態様に係る波形再構成システム及び波形再構成装置は、ある複数のスペクトル変化の組合せを引き起こす時間波形は限定されるという非線形光学効果の特徴を利用し、入力光信号の時間波形の再構成を行う。具体的には、波形再構成システム及び波形再構成装置は、入力光信号の強度を変化させたときに得られる、非線形光学効果による複数のスペクトル変化を、実測値と計算値とで比較し、その比較結果を利用して入力光信号の時間波形の再構成を行う。時間波形の再構成に必要な位相スペクトルは、焼きなまし法などの最適化アルゴリズムに従ってシミュレーションにより計算された計算パワースペクトル（計算値）が計測パワースペクトル（実測値）に近づくように位相スペクトルを変化させていくことにより決定される。
【００２７】
以下に、本発明の実施の形態について、図面を参照しながら説明する。なお、以下で説明する実施の形態は、いずれも本発明の好ましい一具体例を示す。つまり、以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置及び接続形態、ステップ、ステップの順序などは、本発明の一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素については、より好ましい形態を構成する任意の構成要素として説明される。
【００２８】
（実施の形態１）
＜波形再構成システムの構成＞
図１は、本発明の実施の形態１に係る波形再構成システム１００の全体的な構成を示す図である。この波形再構成システム１００は、入力光信号の時間波形を再構成するシステムである。図１に示すように、波形再構成システム１００は、強度調節器１１０、高非線形光ファイバ１２０、分光器１３０及び波形再構成装置１４０を備える。
【００２９】
強度調節器１１０は、光信号生成装置２００により生成された入力光信号の強度を変化させる。
【００３０】
高非線形光ファイバ１２０は、非線形光学効果を有する光伝送媒体の一例である。強度調節器１１０により強度が変化された入力光信号は、この高非線形光ファイバ１２０内を伝播する。
【００３１】
高非線形光ファイバ１２０は、非線形光学効果に関するパラメータが既知のファイバである。本実施の形態では、高非線形光ファイバ１２０は、分散効果に関するパラメータも既知のファイバである。具体的には、例えば、２次及び３次に加え、４次分散に関するパラメータが既知である。
【００３２】
分光器１３０は、スペクトル計測器の一例であり、出力光信号を波長ごとの光に分解し、波長ごとに分解した光をＯ／Ｅ変換及びＡ／Ｄ変換することにより、デジタル値で表されたパワースペクトルを生成する。すなわち、分光器１３０は、出力光信号のパワースペクトルを計測する。ここで出力光信号とは、入力光信号が高非線形光ファイバ１２０内を伝播した後に出力される光信号である。
【００３３】
波形再構成装置１４０は、例えば、図１６に示すようなコンピュータにより実現され、入力光信号の時間波形を再構成する。波形再構成装置１４０の詳細は図２を用いて後述する。
【００３４】
光信号生成装置２００は、入力光信号を生成する。具体的には、光信号生成装置２００は、例えば、ＭＬＬＤ（Ｍｏｄｅ−Ｌｏｃｋｅｄ  Ｌａｓｅｒ  Ｄｉｏｄｅ）と、ＳＭＦ（Ｓｉｎｇｌｅ  Ｍｏｄｅ  Ｆｉｂｅｒ）と、ＥＤＦＡ（Ｅｒｂｉｕｍ  Ｄｏｐｅｄ  Ｆｉｂｅｒ  Ａｍｐｌｉｆｉｅｒ）とを備える。光信号生成装置２００は、ＭＬＬＤから出力された光パルスを、ＳＭＦにより分散補償し、ＥＤＦＡにより増幅する。
【００３５】
＜波形再構成装置の構成＞
図２は、本発明の実施の形態１に係る波形再構成装置１４０の特徴的な機能構成を示すブロック図である。図２に示すように、入力スペクトル取得部１４１、出力スペクトル取得部１４２、位相スペクトル算出部１４３及び波形再構成部１４４を備える。
【００３６】
入力スペクトル取得部１４１は、光信号生成装置２００により生成された入力光信号のパワースペクトルを取得する。例えば、光信号生成装置２００により生成された入力光信号のパワースペクトルが既知である場合、入力スペクトル取得部１４１は、記憶手段等に格納されたパワースペクトルのデータを読み出すことにより、入力光信号のパワースペクトルを取得する。一方、光信号生成装置２００により生成された入力光信号のパワースペクトルが既知でない場合、入力スペクトル取得部１４１は、図示していない分光器等を用いて計測された入力光信号のパワースペクトルを取得する。なお、入力スペクトル取得部１４１が取得するパワースペクトルとは、入力光信号の波長ごとの光の強度を示すデータである。
【００３７】
出力スペクトル取得部１４２は、入力光信号が高非線形光ファイバ１２０内を伝播した後に出力される出力光信号のパワースペクトルであって、分光器１３０により計測されたパワースペクトルを計測パワースペクトルとして複数の強度の入力光信号ごとに取得する。ここで、計測パワースペクトルとは、出力光信号の波長ごとの光の強度を示すデータである。
【００３８】
位相スペクトル算出部１４３は、入力光信号の強度ごとに、入力光信号が所与の位相スペクトルを有すると仮定して入力光信号の高非線形光ファイバ１２０内の伝播シミュレーションを行うことにより、出力光信号の計算パワースペクトルを算出する。そして、位相スペクトル算出部１４３は、計算パワースペクトルと、計測パワースペクトルとの強度ごとの差分値が小さくなるように所与の位相スペクトルを変化させて伝播シミュレーションを行うことにより、入力光信号の位相スペクトルを探索する。例えば、位相スペクトル算出部１４３は、所定のアルゴリズムに従って所与の位相スペクトルを変化させる。
【００３９】
ここで、所定のアルゴリズムとは、与えられた関数の最適解を探索するためのアルゴリズムである。例えば、所定のアルゴリズムは、焼きなまし法、共役方向法、共役勾配法、遺伝的アルゴリズムなどである。なお、所定のアルゴリズムは、これらのアルゴリズムに限定される必要はなく、どのようなアルゴリズムであってもよい。
【００４０】
なお、位相スペクトル算出部１４３が探索する位相スペクトルとは、光信号の波長ごとの位相を示すデータである。
【００４１】
また、位相スペクトル算出部１４３により実行される伝播シミュレーションは、例えば、スプリットステップフーリエ法などによるパルス伝播シミュレーションである。すなわち、位相スペクトル算出部１４３により実行される伝播シミュレーションは、高非線形光ファイバ１２０に固有の既知のパラメータを利用して、出力光信号のパワースペクトルを計算するための光信号伝播シミュレーションである。ここで用いられるパラメータは、光学非線形効果に関するパラメータを含む。本実施の形態では、パラメータは、分散効果に関するパラメータも含む。
【００４２】
また、所与の位相スペクトルとは、伝播シミュレーションにおいて与えられる位相スペクトルであって、入力光信号の仮想的な位相スペクトルである。つまり、所与の位相スペクトルとは、伝播シミュレーションで用いられる位相スペクトルである。最も評価が高い伝播シミュレーションで用いられた所与の位相スペクトルが、入力光信号の位相スペクトルとして決定される。
【００４３】
波形再構成部１４４は、位相スペクトル算出部１４３により探索された位相スペクトルと、入力スペクトル取得部１４１により取得されたパワースペクトルとを周波数／時間変換することにより、入力光信号の時間波形を再構成する。具体的には、波形再構成部１４４は、例えば、位相スペクトルとパワースペクトルとを逆フーリエ変換することにより、時間波形を再構成する。
【００４４】
＜波形再構成装置の処理動作＞
次に、以上のように構成された波形再構成装置１４０における各種動作について説明する。
【００４５】
図３は、本発明の実施の形態１に係る波形再構成装置１４０における動作を示すフローチャートである。
【００４６】
まず、入力スペクトル取得部１４１は、光信号生成装置２００により生成された入力光信号のパワースペクトルを取得する（ステップＳ１０１）。続いて、出力スペクトル取得部１４２は、分光器１３０により計測された、複数の強度の入力光信号にそれぞれ対応する出力光信号のパワースペクトルを計測パワースペクトルとして取得する（ステップＳ１０２）。
【００４７】
次に、位相スペクトル算出部１４３は、伝播シミュレーションで用いられる所与の位相スペクトルの初期値を設定する（ステップＳ１０３）。例えば、位相スペクトル算出部１４３は、任意の位相スペクトルを初期値として設定する。また、例えば、位相スペクトル算出部１４３は、入力光信号が予め定められた型のパルスであると仮定したときに得られる位相スペクトルを初期値として設定してもよい。
【００４８】
そして、位相スペクトル算出部１４３は、入力光信号の強度ごとに、入力光信号が所与の位相スペクトルを有すると仮定して入力光信号の高非線形光ファイバ１２０内の伝播シミュレーションを行うことにより、出力光信号の計算パワースペクトルを算出する（ステップＳ１０４）。伝播シミュレーションの詳細は後述する。
【００４９】
次に、位相スペクトル算出部１４３は、位相スペクトルの探索を終了するか否かを判定する（ステップＳ１０５）。例えば、位相スペクトル算出部１４３は、計算された複数の強度ごとの計算パワースペクトルと、分光器１３０により計測された（ステップＳ１０２において取得された）複数の強度ごとの計測パワースペクトルとの、各強度における差分値の合計が所定閾値以下か否かに基づいて、探索を終了するか否かを判定する。また例えば、位相スペクトル算出部１４３は、計算パワースペクトルと計測パワースペクトルとの差分の変化率に基づいて、探索を終了するか否かを判定してもよい。また例えば、位相スペクトル算出部１４３は、伝播シミュレーションの反復回数が上限回数に達したか否かに基づいて、探索を終了するか否かを判定してもよい。
【００５０】
ここで、探索を終了しないと判定された場合（ステップＳ１０５のＮｏ）、位相スペクトル算出部１４３は、計算パワースペクトルと計測パワースペクトルとの差分が小さくなるように所与の位相スペクトルを変化させる（ステップＳ１０６）。具体的には、位相スペクトル算出部１４３、伝播シミュレーションに用いる新たな所与の位相スペクトルを、例えば焼きなまし法に従って設定する。そして、再度、ステップＳ１０４からの処理が繰り返される。
【００５１】
一方、探索を終了すると判定された場合（ステップＳ１０５のＹｅｓ）、位相スペクトル算出部１４３は、最も小さな差分が得られた伝播シミュレーションで用いられた所与の位相スペクトルを入力光信号の位相スペクトルとして決定する（ステップＳ１０７）。
【００５２】
このように、位相スペクトル算出部１４３は、探索を終了すると判定されるまで、所与の位相スペクトルを変化させながら伝播シミュレーションを繰り返すことにより、入力光信号の位相スペクトルを探索する。つまり、位相スペクトル算出部１４３は、計算パワースペクトルと計測パワースペクトルとの強度ごとの差分を評価関数として用いて、位相スペクトルの最適解を探索する。
【００５３】
最後に、入力光信号のパワースペクトルと、決定された位相スペクトルとを周波数／時間変換することにより、入力光信号の時間波形を再構成し（ステップＳ１０８）、処理を終了する。
【００５４】
以上の処理により、波形再構成装置１４０は、非線形光学効果を有する光伝送媒体内を入力光信号が伝播した後に出力される出力光信号のパワースペクトルを用いて、入力光信号の時間波形を簡便に再構成することが可能となる。
【００５５】
＜伝播シミュレーション＞
次に、スプリットステップフーリエ法による伝播シミュレーションの詳細について説明する。まず、高非線形光ファイバ１２０内を伝播する光信号に影響を与える分散効果と非線形光学効果とについて説明する。
【００５６】
分散効果とは、光と物質との相互作用が光の波長によって異なることに起因する現象である。分散効果によって、入力光信号の伝播速度が周波数に応じて変化する。つまり、入力光信号に含まれる各周波数成分の位相がずれ、入力光信号の時間波形が広がる。
【００５７】
また、非線形光学効果とは、光（例えば超短光パルスのような非常に強い強度を有する光）と物質とが相互作用することにより引き起こされる非線形の多様な現象である。非線形光学効果としては、自己位相変調、自己急峻化、及びラマン応答などが挙げられる。
【００５８】
まず、自己位相変調について説明する。光ファイバなどの光伝送媒体の屈折率は、その中を伝播する光信号の強度に比例してわずかに変化するため、光信号自身に位相変調が生じる。このようにして生じる位相変調を自己位相変調という。
【００５９】
次に、自己急峻化について説明する。自己急峻化とは、時間波形が対称である入力光信号が光伝送媒体中を伝播するにつれて、時間波形が非対称となり、ピークが後ろの方へ移る現象である。自己急峻化は、群速度が強度に依存することから引き起こされる。
【００６０】
最後に、ラマン応答について説明する。物質に光が入射した場合、入射光と等しい周波数をもつ強い弾性散乱（レーリー散乱)光と、入射光の周波数からわずかにずれた周波数をもつきわめて弱い非弾性散乱光とが散乱される。非弾性散乱光は、物質中の振動する原子やイオンにより散乱されるラマン散乱光と、物質中の音波によって散乱されるブリユアン散乱光とに分けられる。ラマン散乱において入射光の強度がある閾値を超えた場合にラマン散乱光が誘導放出により強く発生する現象を誘導ラマン散乱という。この誘導ラマン散乱により光伝送媒体内では光の高周波数成分から低周波数成分へエネルギーが移され、低周波成分が強められる。このような現象をラマン応答と呼ぶ。
【００６１】
以上のような分散効果と非線形光学効果との影響を受けて高非線形光ファイバ１２０内を伝播する入力光信号の伝播方程式は、式（１）で表される。
【００６２】
【数１】
0006025013.tif 000002
【００６３】
ここで、Ｅは光の電場成分を表し、ｚは高非線形光ファイバ１２０内の距離を表す。また、Ｄは分散効果とロスとを表し、Ｎは非線形光学効果を表す。Ｄ及びＮは、次の式（２）及び式（３）で表される。
【００６４】
【数２】
0006025013.tif 000003
【数３】
0006025013.tif 000004
【００６５】
ここで、α、β、γは、高非線形光ファイバ１２０に固有のパラメータを表す。具体的には、αは、光の強度のロスに関するパラメータを表す。また、βは、各次数の分散効果に関するパラメータを表す。また、γは、非線形光学効果に関するパラメータを表す。なお、Ｔは時間を表す。
【００６６】
式（２）の右辺において、第１項は２次分散を表し、第２項は３次分散を表し、第３項は４次分散を表し、最後の項は高非線形光ファイバ１２０内を伝播する光の強度のロスを表している。また、式（３）の右辺において、第１項は自己位相変調を表し、第２項は自己急峻化を表し、第３項はラマン応答を表している。
【００６７】
この伝播方程式には、Ｅ自体に依存する項が含まれている。そのため、分散項（Ｄ）と非線形項（Ｎ）とを同時に計算することは困難である。そこで、スプリットステップフーリエ法では、図４に示すように、光信号の伝播方向に沿って高非線形光ファイバ１２０を短いステップに仮想的に区切る。そして、ステップごとに分散項と非線形項とを順に計算することにより、高非線形光ファイバ１２０中を伝播する光信号の近似解を求めることができる。
【００６８】
＜実験結果＞
次に、本実施の形態に係る波形再構成システム１００により得られた実験結果について説明する。
【００６９】
本実験では、光信号生成装置２００は、ＭＬＬＤから１０ＧＨｚの周期で出力された１．３ｐｓの光パルスを、６５ｍのＳＭＦにより分散補償した後に、ＥＤＦＡにより増幅した。このようにＥＤＦＡにより増幅された光パルスを、光信号生成装置２００は入力光信号として出力した。
【００７０】
強度調節器１１０には、可変光減衰器（ＶＯＡ：Ｖａｒｉａｂｌｅ  Ｏｐｔｉｃａｌ  Ａｔｔｅｎｕａｔｏｒ）を用いた。強度調節器１１０は、電力を１７．２ｍＷから１０．３ｍＷずつ上げていくことにより、光信号生成装置２００により生成された入力光信号の強度を１０種類の異なる強度に変更した。そして、強度調節器１１０は、強度を変更した入力光信号を高非線形光ファイバ１２０に出力した。
【００７１】
高非線形光ファイバ１２０には、２次及び３次に加え、４次分散まで考慮した、表１に示すパラメータを有する光ファイバを用いた。
【００７２】
【表１】
0006025013.tif 000005
【００７３】
分光器１３０により計測された、入力光信号の１０種類の強度に対応する出力光信号のパワースペクトルを、図５の（ａ）〜（ｊ）に示す。図において、グラフの横軸は波長（Ｗａｖｅｌｅｎｇｔｈ（ｎｍ））を示し、縦軸は強度（Ｉｎｔｅｎｓｉｔｙ（ａ．ｕ．））を示す。また、図５の（ａ）〜（ｊ）に示すパワースペクトルのそれぞれは、強度調節器１１０の電力が順に１７．２ｍＷ、２７．５ｍＷ、３７．８ｍＷ、４８．１ｍＷ、５８．４ｍＷ、６８．７ｍＷ、７９．０ｍＷ、８９．４ｍＷ、９９．７ｍＷ、及び１１０．０ｍＷであるときの出力光信号のパワースペクトルである。
【００７４】
本実験では、波形再構成装置１４０が備える位相スペクトル算出部１４３は、７種類の強度に対応する（図５の（ａ）〜（ｇ）に示す）出力光信号のパワースペクトルの実測値と、シミュレーションにより得られる出力光信号のパワースペクトルの計算値との差分値を評価関数として、焼きなまし法に従って位相スペクトルを変化させながら評価関数が最小となる位相スペクトルを探索した。その際、位相スペクトルの初期値には、オートコリレーションにより求めたパルス幅を元にＧｅｒｃｈｂｅｒｇ−Ｓａｘｔｏｎ法によりＳｅｃｈ型光パルスに近似した際の位相スペクトルを用いた。このような位相スペクトルを初期値として利用することにより、早期に評価関数を収束させることが可能となった。
【００７５】
図６は、位相スペクトルの探索中における評価関数の変化の実験結果を示すグラフである。図６において、グラフの横軸は評価関数の繰り返し計算回数（Ｎｕｍｂｅｒ  ｏｆ  ｉｔｅｒａｔｉｏｎ  ｔｉｍｅｓ）を示し、縦軸は評価関数の値（Ｅｖａｌｕａｔｉｏｎ  ｆｕｎｃｔｉｏｎ（ａ．ｕ．））を示す。
【００７６】
図６から明らかなように、比較的早期に評価関数が一定値に収束している。位相スペクトル算出部１４３は、評価関数が一定値に収束したとき（例えば、図６における２５００回目）の位相スペクトルを入力光信号の位相スペクトルと決定した。
【００７７】
図７は、実験により得られたパワースペクトルと位相スペクトルとを示す図である。図７において、横軸は波長（Ｗａｖｅｌｅｎｇｔｈ（ｎｍ））を示す。また、縦軸は左側が強度（Ｉｎｔｅｎｓｉｔｙ（ａ．ｕ．））を示し、右側が位相（Ｐｈａｓｅ（ｒａｄ））を示す。ここで、パワースペクトル６０１が左側縦軸の強度に対応し、位相スペクトル６０２が右側縦軸の位相に対応する。
【００７８】
波形再構成装置１４０が備える波形再構成部１４４は、図７のように算出された位相スペクトル６０２と、入力光信号のパワースペクトル６０１とを用いて、入力光信号の時間波形を再構成する。
【００７９】
図８は、波形再構成装置１４０により再構成された入力光信号の時間波形の実験結果を示す図である。図８において、横軸は時間（Ｄｅｌａｙ  ｔｉｍｅ（ｐｓ））を示し、縦軸は強度（Ｉｎｔｅｎｓｉｔｙ（ａ．ｕ．））を示す。
【００８０】
図８から明らかなように、再構成された入力光信号の時間波形は、一般的なパルス光として近似されるＳｅｃｈ型又はＧａｕｓｓｉａｎ型の時間波形とは異なる。すなわち、波形再構成装置１４０は、一般的な近似により得られる時間波形よりも正確な時間波形を再構成できたと考えられる。
【００８１】
次に、上記実験において位相スペクトル算出部１４３により実行されたシミュレーションと算出された位相スペクトルとの妥当性を検証する。
【００８２】
図９は、７種類の入力光信号の強度に対応する出力光信号のパワースペクトルの実測値と計算値との比較結果を示すグラフである。なお、パワースペクトルの計算値は、繰り返し計算により評価関数が収束したときに得られたパワースペクトルの値である。また、入力光信号の７種類の強度は、強度調節器１１０の電力が１７．２ｍＷ、２７．５ｍＷ、３７．８ｍＷ、４８．１ｍＷ、５８．４ｍＷ、６８．７ｍＷ、及び７９．０ｍＷのときの強度である。
【００８３】
図９の（ａ）〜（ｇ）に示すように、実測値８０１と計算値８０２とは略一致している。したがって、本実験において、入力光信号の時間波形に敏感に反応する、自己位相変調によるパワースペクトルの変化を正確にシミュレーションできていることが分かる。
【００８４】
図１０は、分光器１３０により計測された１０種類の入力光信号の強度に対応する出力光信号のパワースペクトルのうち、位相スペクトルの算出に利用されなかった３種類の入力光信号の強度に対応する出力光信号のパワースペクトルの実測値と計算値との比較結果を示すグラフである。なお、パワースペクトルの計算値は、位相スペクトル算出部１４３により算出された入力光信号の位相スペクトルを用いてシミュレーションしたときの値である。また、位相スペクトルの算出に利用されなかった３種類の入力光信号の強度は、強度調節器１１０の電力が８９．４ｍＷ、９９．７ｍＷ、及び１１０．０ｍＷのときの強度である。
【００８５】
図１０の（ａ）〜（ｃ）に示すように、実測値９０１と計算値９０２とは略一致している。つまり、焼きなまし法における評価関数の対象としなかった入力光信号の強度においても、位相スペクトル算出部１４３により算出された位相スペクトルを用いて出力光信号のパワースペクトルが正確に計算されている。したがって、位相スペクトル算出部１４３により算出された入力光信号の位相スペクトルの精度が高いことがわかる。つまり、位相スペクトル算出部１４３により算出された入力光信号の位相スペクトルを用いて再構成される入力光信号の時間波形の精度も高いということができる。
【００８６】
以上のように、本実施の形態に係る波形再構成システム１００及び波形再構成装置１４０は、自己位相変調を誘起する光ファイバ内を入力光信号が伝播した後に出力される出力光信号のパワースペクトルを用いて、入力光信号の時間波形を再構成することができる。したがって、波形再構成システム１００及び波形再構成装置１４０では、時間分解分光データが必ずしも必要ではない。すなわち、波形再構成システム１００は、光信号の時間波形を再構成するために、時間分解分光データを得るために必要な超高速時間ゲート又は参照光源が必ずしも備える必要がないので、光信号の時間波形を簡便に再構成することが可能となる。
【００８７】
さらに、波形再構成システム１００及び波形再構成装置１４０は、公知の技術である焼きなまし法に従って位相スペクトルを変化させながら、所定閾値以下となる差分値を探索することができるので、比較的短時間に精度の高い位相スペクトルを算出することが可能となる。
【００８８】
（実施の形態２）
次に、本発明の実施の形態２について説明する。
【００８９】
本実施の形態では、非線形光学効果及び分散効果のうちの分散効果のみに基づいて、所与の位相スペクトルを変化させる点が、実施の形態１と異なる。具体的には、波形再構成装置は、位相スペクトルの探索において探索時間を短縮するために、制約条件に従って所与の位相スペクトルを変化させる。以下、本実施の形態に係る波形再構成装置について、実施の形態１と異なる点を中心に説明する。
【００９０】
なお、本実施の形態に係る波形再構成システムの構成と、波形再構成装置の機能構成及び処理動作とは、実施の形態１の図１〜図３と同様であるので図示を省略する。
【００９１】
本実施の形態において、位相スペクトル算出部１４３は、非線形光学効果及び分散効果のうち分散効果のみに基づいて、所与の位相スペクトルを変化させる。具体的には、位相スペクトル算出部１４３は、制約条件に従って、所与の位相スペクトルを変化させる。この制約条件とは、入力光信号が非線形光学効果及び分散効果のうちの分散効果のみを有する媒体内を伝播して得られた信号であると仮定した場合に満たされるべき位相スペクトルの制約条件である。
【００９２】
ここで、制約条件ついてより詳細に説明する。
【００９３】
レーザーから発振した直後の超短光パルスの波形は、歪みのないフーリエ変換限界パルス（ＴＬＰ：Ｔｒａｎｓｆｏｒｍ−ＬｉｍｉｔｅｄＰｕｌｓｅ）であると仮定し、かつ、再構成の対象となる入力光信号は、ＴＬＰが媒体内を伝播した後に出力される光信号であると仮定する。この場合、ＴＬＰは、式（１）に従って媒体内を伝播する。
【００９４】
しかし、ＴＬＰが伝播する媒体が一般的な媒体（つまり、特に高い非線形光学性を有していない媒体）であれば、非線形光学効果がＴＬＰの伝播に及ぼす影響は小さい。つまり、式（１）においてＮを省略して伝播計算を行っても、計算結果は大きく変わらない。したがって、入力光信号は、分散効果のみを有する媒体内をＴＬＰが伝播した後に出力される光信号として表すことができる。ここで媒体内における分散効果は、テイラー展開により式（４）で表される。
【００９５】
【数４】
0006025013.tif 000006
【００９６】
式（４）において、０次及び１次の項は、時間波形の時間シフトに影響を与えるだけであり、位相に影響を与えない。また、一般的に光ファイバなどの媒体内の光信号の伝播シミュレーションでは、４次分散まで考慮されれば十分である。そこで、分散効果によって歪んだ入力光信号の位相は、式（５）のような多項式関数として表すことができる。
【００９７】
【数５】
0006025013.tif 000007
【００９８】
位相スペクトル算出部１４３は、この式（５）が満たされるように、所与の位相スペクトルを設定する。つまり、位相スペクトル算出部１４３は、式（５）で示す制約条件に従って、所与の位相スペクトルを変化させる。
【００９９】
式（５）によれば、位相スペクトルは、ａ、ｂ、ｃの３個の変数で表される。一方、実施の形態１では、位相スペクトルは、周波数成分ごとの変数（例えば５１２個の変数）で表される。つまり、本変形例によれば、位相スペクトルの探索において、実施の形態１よりも、変数の数を減少させることができ、探索空間を縮小することができる。したがって、本変形例によれば、位相スペクトルの探索時間を短縮することが可能となる。
【０１００】
以上のように、本実施の形態に係る波形再構成装置１４０によれば、適切な制約条件に従って所与の位相スペクトルを変化させることができる。したがって、位相スペクトルの探索において、変数の数を減少させることができるので、探索時間を短縮することが可能となる。
【０１０１】
なお、ここでは、制約条件として、２次分散から４次分散までを考慮していたが、さらに高次の分散まで考慮しても構わない。例えば、位相スペクトルを探索する際に、制約条件として考慮する分散の次数を段階的に増やしても構わない。これにより、探索空間を段階的に広げることができ、探索時間の短縮と探索精度の向上とを図ることが可能となる。
【０１０２】
（実施の形態２の変形例）
次に、本発明の実施の形態２の変形例について説明する。
【０１０３】
実施の形態２では、入力光信号が分散効果のみを有する媒体内を伝播した後に出力される信号であると仮定していた。したがって、入力光信号が分散効果以外の影響を受けていた場合、探索精度が低下する。そこで、本変形例では、位相スペクトルの探索精度を向上させるために、実施の形態１及び実施の形態２の両方の方法を組み合わせて、位相スペクトルの探索を行う。
【０１０４】
具体的には、位相スペクトル算出部１４３は、入力光信号の位相スペクトルを探索する際に、制約条件に従って所与の位相スペクトルを変化させて伝播シミュレーションを行った後に、制約条件に関わらず所与の位相スペクトルを変化させて伝播シミュレーションを行う。ここで、制約条件は、実施の形態２の制約条件と同様である。
【０１０５】
つまり、位相スペクトル算出部１４３は、まず、実施の形態２に示すように制約条件に従って所与の位相スペクトルを変化させて、位相スペクトルを探索する。その後、位相スペクトル算出部１４３は、探索結果を所与の位相スペクトルの初期値として用いて、実施の形態１に示すように制約条件に関わらず位相スペクトルを変化させて、位相スペクトルを探索する。
【０１０６】
以上のように、本変形に係る波形再構成装置１４０によれば、制約条件に従って所与の位相スペクトル変化させた後に、制約条件に関わらず所与の位相スペクトルを変化させることができる。したがって、位相スペクトルの探索において、探索精度の低下を抑制しつつ、探索時間を短縮することが可能となる。
【０１０７】
＜実験結果＞
次に、本変形例に係る波形再構成装置１４０により得られた実験結果について説明する。
【０１０８】
本実験では、光信号生成装置２００として、中心波長１５４６ｎｍ、繰り返し周波数２５ＭＨｚの超短光パルスを生成するファイバレーザを用いた。また、高非線形光ファイバ１２０として、１０ｍの長さのものを用い、実施の形態１よりも大きな非線形光学効果を有する光ファイバ（γ＝２１／Ｗ／ｋｍ）を用いた。そして、生成された超短光パルスの強度を非線形光学効果が起きないほど弱めた光信号を入力光信号として用いた。
【０１０９】
図１１は、本発明の実施の形態２の変形例の実験で用いた入力光信号のパワースペクトルを示すグラフである。図１２の（ａ）及び（ｂ）は、２つの強度の入力光信号が高非線形光ファイバ１２０を伝播した後に分光器１３０により計測されたパワースペクトルを示すグラフである。図１２に示すように、本実験では、入力光信号の強度を２段階に変化させた。
【０１１０】
評価関数としては、分光器１３０で計測された計測パワースペクトルと、所与の位相スペクトルを用いた伝播シミュレーションによって得られた計算パワースペクトルとの周波数ごとの差分の和を正規化して用いた。位相スペクトル算出部１４３は、この評価関数が小さくなるように所与の位相スペクトルを変化させながら伝播シミュレーションを繰り返すことにより、位相スペクトルを探索した。このとき、所与の位相スペクトルの初期値として、全周波数で成分値が０の位相スペクトルを用いた。また、探索時の伝播シミュレーションの反復回数を１０００回とした。
【０１１１】
このような条件で行われた実験結果を図１３及び図１４を用いて説明する。
【０１１２】
図１３及び図１４は、本発明の実施の形態２の変形例に係る波形再構成装置による実験結果を示す図である。
【０１１３】
図１３において、横軸は波長を示し、縦軸は強度（左側）又は位相（右側）を示す。そして、破線は、パワースペクトルを表し、強度の縦軸に対応する。また、実線は、位相スペクトルを表し、位相の縦軸に対応する。また、図１４において、横軸は反復回数を示し、縦軸は評価関数を示す。
【０１１４】
また、図１３及び図１４において、（ａ）は、実施の形態１に係る探索方法によって得られた実験結果を示し、（ｂ）は、本変形例に係る探索方法によって得られた実験結果を示す。
【０１１５】
なお、実施の形態１に係る探索方法では、所与の位相スペクトルの各周波数成分値は、１０００回の伝播シミュレーションのすべてで互いに独立して設定された。また、本変形例に係る探索方法では、所与の位相スペクトルの各周波数成分値は、１０００回のうちのはじめの５２回の伝播シミュレーションで、式（５）を満たすように各周波数成分値が設定された。また、残りの９４８回の伝播シミュレーションで、所与の位相スペクトルの各周波数成分値は、式（５）に関わらず互いに独立して設定された。
【０１１６】
図１３及び図１４から明らかなように、本変形例に係る探索方法では、実施の形態１に係る探索方法よりも、解が収束しやすくなり、探索速度が向上している。これは、制約条件によって、適切に探索空間を小さくすることができているからであると考えられる。
【０１１７】
なお、本実験において、実施の形態１に係る探索方法で１０００回の伝播シミュレーションによって得られた評価関数の値は、本変形例に係る方法で２９３回の伝播シミュレーションによって得られた評価関数の値と同じであった。つまり、本変形例に係る探索方法は、実施の形態１に係る探索方法よりも約３．４倍の高速化を実現することができ、探索時間を短縮することができた。
【０１１８】
また、本実験のように、所与の位相スペクトルの初期値が最適解と大きく離れている場合には、本変形例に係る探索方法の方が、実施の形態１に係る探索方法よりも探索精度が向上する。これは、制約条件によって適切に探索空間を小さくすることができるので、局所的最適解に収束することを抑制することができるからであると考えられる。
【０１１９】
以上のように、本実験により、制約条件に従って所与の位相スペクトルを変化させることにより、探索精度を向上させることができるとともに、探索時間を短縮することができることを確認できた。
【０１２０】
（実施の形態３）
次に、本発明の実施の形態３について説明する。
【０１２１】
本実施の形態では、非線形光学効果及び分散効果のうち非線形光学効果のみに基づいて伝播シミュレーションを行う点が、実施の形態１と主として異なる。以下、本実施の形態に係る波形再構成装置について、実施の形態１と異なる点を中心に説明する。
【０１２２】
なお、本実施の形態に係る波形再構成システムの構成と、波形再構成装置の機能構成及び処理動作とは、実施の形態１の図１〜図３と同様であるので図示を省略する。
【０１２３】
本実施の形態において、位相スペクトル算出部１４３は、非線形光学効果及び分散効果のうち非線形光学効果のみに基づいて、伝播シミュレーションを行う。具体的には、位相スペクトル算出部１４３は、非線形光学効果と分散効果とのうち非線形光学効果に関するパラメータのみを用いて伝播シミュレーションを行う。具体的には、位相スペクトル算出部１４３は、例えば、式（２）において２次分散、３次分散及び４次分散を表す項（第１項〜第３項）を省略して、伝播シミュレーションを行う。つまり、位相スペクトル算出部１４３は、以下の式（６）〜式（８）に基づいて、伝播シミュレーションを行う。
【０１２４】
【数６】
0006025013.tif 000008
【０１２５】
なお、非線形光学効果が大きい光伝送媒体内の光信号の伝播では、非線形光学効果が支配的であり、分散効果は限定的である。つまり、入力光信号の高非線形光ファイバ１２０内の伝播シミュレーションに対して分散効果が与える影響は小さい。したがって、非線形光学効果に関するパラメータのみを用いて、入力光信号の高非線形光ファイバ１２０内の伝播シミュレーションを行っても、シミュレーションの精度が低下することを抑制することができる。
【０１２６】
また、このように分散効果を表す項を省略することで、式（６）において、ＤとＮとを同時に計算することが可能となる。つまり、位相スペクトル算出部１４３は、スプリットステップフーリエ法を用いずに、１度の計算で伝播シミュレーションを行うことができる。したがって、伝播シミュレーションで計算速度を大幅に向上させることが可能となる。
【０１２７】
なお、高非線形光ファイバ１２０は、分散効果が「０」となる波長（以下、「ゼロ分散波長」という）を有する。入力光信号の周波数成分がゼロ分散波長に集中するほど、分散効果が光信号の伝播に与える影響が小さくなる。したがって、入力光信号の中心波長と高非線形光ファイバ１２０のゼロ分散波長との差分が小さいほど、分散効果を省略しても高い精度で伝播シミュレーションを行うことができる。
【０１２８】
したがって、入力光信号の中心波長が既知である場合には、その中心波長と近似する値のゼロ分散波長を有する高非線形光ファイバ１２０を用いることが好ましい。これにより、非線形光学効果と分散効果とのうち非線形光学効果に関するパラメータのみを用いて伝播シミュレーションを行った場合に、シミュレーション結果の精度の低下を抑制することができる。
【０１２９】
なお、実験によれば、高非線形光ファイバ１２０のゼロ分散波長が１５５０ｎｍである場合に、入力光信号の中心波長が１５２０ｎｍ〜１５８０ｎｍであれば、シミュレーションの精度の低下を効果的に抑制することができた。
【０１３０】
以上のように、本実施の形態に係る波形再構成装置１４０によれば、非線形光学効果と分散効果とのうち非線形光学効果に関するパラメータのみを用いて伝播シミュレーションを行うことができる。したがって、伝播シミュレーションが単純化されるので、計算速度を向上させることが可能となる。
【０１３１】
（実施の形態３の変形例）
次に、実施の形態３の変形例について説明する。
【０１３２】
上記実施の形態３では、自己位相変調、自己急峻化、及びラマン応答を考慮して伝播シミュレーションを行った。自己急峻化及びラマン応答を表す項（式（８）の右辺の第２項及び第３項）には、時間波形の微分が含まれている。時間波形の微分とは、ある時間における時間波形の傾きに相当する。つまり、パルス幅が狭いほど、自己急峻化及びラマン応答が伝播に及ぼす影響が大きくなる。逆に言えば、パルス幅が広い場合には、自己急峻化及びラマン応答の項を省略して伝播シミュレーションを行っても、シミュレーション結果は大きく変化しないことが考えられる。
【０１３３】
そこで、本変形例において、位相スペクトル算出部１４３は、非線形光学効果に関するパラメータのうち、自己位相変調に関するパラメータのみを用いて、伝播シミュレーションを行う。具体的には、位相スペクトル算出部１４３は、例えば、式（３）において自己急峻化及びラマン応答を表す項（第２項及び第３項）を省略した式を用いて、伝播シミュレーションを行う。つまり、位相スペクトル算出部１４３は、以下の式（９）〜式（１１）に基づいて、伝播シミュレーションを行う。
【０１３４】
【数７】
0006025013.tif 000009
【０１３５】
以上のように、本変形例に係る波形再構成装置１４０によれば、上記実施の形態３よりもさらに伝播シミュレーションを簡便化でき、伝播シミュレーションの計算速度をさらに向上させることができる。特に、入力光信号のパルス幅が広い場合には、シミュレーション精度の低下も抑制することができ、効果的である。なお、実験によれば、パルス幅が５００ｆｓ以上であれば、シミュレーションの精度の低下を効果的に抑制することができた。
【０１３６】
＜実験結果＞
次に、本変形例に係る波形再構成装置１４０により得られた実験結果について説明する。ここでは、実施の形態２の変形例における実験と同様の条件で実験を行った。なお、本実験において、高非線形光ファイバ１２０のゼロ分散波長は１５５０ｎｍであった。なお、探索方法としては、実施の形態２の変形例に係る探索方法を用いた。つまり、はじめに制約条件に従って所与の位相スペクトルを変化させた。
【０１３７】
図１５は、本発明の実施の形態３の変形例に係る波形再構成装置による実験結果を示す図である。
【０１３８】
図１５において、横軸は波長を示し、縦軸は強度（左側）又は位相（右側）を示す。そして、破線は、パワースペクトルを表し、強度の縦軸に対応する。また、実線は、位相スペクトルを表し、位相の縦軸に対応する。また、図１５の（ａ）は、スプリットステップフーリエ法に基づく伝播シミュレーションによる実験結果を示す。また、図１５の（ｂ）は、本変形例に係る伝播シミュレーションによる実験結果を示す。
【０１３９】
スプリットステップフーリエ法に基づく伝播シミュレーションでは、式（１）〜式（３）に基づいて計算を行った。また、本変形例に係る伝播シミュレーションでは、式（９）〜式（１１）に基づいて計算を行った。
【０１４０】
図１５に示すように、伝播シミュレーションを繰り返して探索された位相スペクトルは、２つの伝播シミュレーションでほぼ同じであった。つまり、分散効果などを無視して伝播シミュレーションを行っても、探索結果はあまり変わらなかった。
【０１４１】
また、本変形例に係る伝播シミュレーションでは、スプリットステップフーリエ法に基づく伝播シミュレーションよりも、大幅に計算時間を削減することができた。具体的には、１０００回の伝播シミュレーションに要した時間が、スプリットステップフーリエ法に基づく伝播シミュレーションでは９９７０ｍｓであったのに対して、本変形例に係る伝播シミュレーションではわずか５６７ｍｓであった。
【０１４２】
つまり、本実験により、非線形光学効果のうち自己位相変調のみを考慮して伝播シミュレーションを行うことにより、精度の低下を抑制しつつ、計算時間を短縮することができることを確認できた。
【０１４３】
以上、本発明の一態様に係る波形再構成装置及び波形再構成システムについて、実施の形態に基づいて説明したが、本発明は、これらの実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したもの、あるいは異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の範囲内に含まれる。
【０１４４】
例えば、上記実施の形態２又はその変形例と、上記実施の形態３又はその変形例とは、組みわせられてもよい。つまり、波形再構成装置は、制約条件に従って所与の位相スペクトルを変化させながら、非線形光学効果及び分散効果のうち非線形光学効果のパラメータのみを用いて伝播シミュレーションを繰り返してもよい。これにより、位相スペクトルの探索時間をさらに短縮することが可能となる。
【０１４５】
また、上記実施の形態１〜３において、波形再構成システムは、高非線形光ファイバを備えていたが、光ファイバとは異なる光伝送媒体を備えてもよい。つまり、波形再構成システムは、非線形光学効果を有する光伝送媒体を備えればよい。
【０１４６】
また、本発明は、このような波形再構成装置の特徴的な構成要素が行う処理を実行する波形再構成方法として実現することもできる。また、その波形再構成方法を図１６に示すようなコンピュータに実行させるためのプログラムとして実現することもできる。そして、そのようなプログラムは、ＣＤ−ＲＯＭ等の記録媒体、インターネット等の伝送媒体を介して配信することができる。
【０１４７】
図１６は、コンピュータのハードウェア構成の一例を示す図である。波形再構成方法をコンピュータに実行させるためのプログラムは、例えば、コンピュータが読取可能な媒体であるＣＤ−ＲＯＭ５１５に記憶され、ＣＤ−ＲＯＭ装置５１４を介して読み出される。また例えば、波形再構成方法をコンピュータに実行させるためのプログラムは、有線若しくは無線ネットワーク、又は放送などを介して伝送される。
【０１４８】
コンピュータ５００は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）５０１、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）５０２、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）５０３、ハードディスク５０４、通信インタフェース５０５等を備える。
【０１４９】
ＣＰＵ５０１は、ＣＤ−ＲＯＭ装置５１４を介して読み出されたプログラム、又は通信インタフェース５０５を介して受信したプログラムを実行する。具体的には、ＣＰＵ５０１は、ＣＤ−ＲＯＭ装置５１４を介して読み出されたプログラム、又は通信インタフェース５０５を介して受信したプログラムをＲＡＭ５０３に展開する。そして、ＣＰＵ５０１は、ＲＡＭ５０３に展開されたプログラム中のコード化された各命令を実行する。
【０１５０】
ＲＯＭ５０２は、コンピュータ５００の動作に必要なプログラム及びデータを記憶する読み出し専用メモリである。ＲＡＭ５０３は、ＣＰＵ５０１がプログラムを実行するときにワークエリアとして使用される。具体的には、ＲＡＭ５０３は、例えば、プログラム実行時のパラメータなどのデータを一時的に記憶する。ハードディスク５０４は、プログラム、データなどを記憶する。
【０１５１】
通信インタフェース５０５は、ネットワークを介して他のコンピュータとの通信を行なう。バス５０６は、ＣＰＵ５０１、ＲＯＭ５０２、ＲＡＭ５０３、ハードディスク５０４、通信インタフェース５０５、ディスプレイ５１１、キーボード５１２、マウス５１３及びＣＤ−ＲＯＭ装置５１４を相互に接続する。
【産業上の利用可能性】
【０１５２】
本発明の一態様に係る波形再構成装置及び波形再構成システムは、光通信等の研究分野で用いられる超高速光パルスの時間波形の測定装置として利用することができる。特に近年、位相を用いた光通信符号の利用及び光パルスの状態に敏感な非線形光信号処理が盛んに研究又は利用されるようになってきていることから、簡便に正確な時間波形を測定可能な波形再構成装置及び波形再構成システムは、産業分野だけでなく基礎科学分野にも非常に有用である。
【符号の説明】
【０１５３】
１００    波形再構成システム
１１０    強度調節器
１２０    高非線形光ファイバ
１３０    分光器
１４０    波形再構成装置
１４１    入力スペクトル取得部
１４２    出力スペクトル取得部
１４３    位相スペクトル算出部
１４４    波形再構成部
２００    光信号生成装置
【図１】
0006025013.tif 000010
【図２】
0006025013.tif 000011
【図３】
0006025013.tif 000012
【図４】
0006025013.tif 000013
【図５】
0006025013.tif 000014
【図６】
0006025013.tif 000015
【図７】
0006025013.tif 000016
【図８】
0006025013.tif 000017
【図９】
0006025013.tif 000018
【図１０】
0006025013.tif 000019
【図１１】
0006025013.tif 000020
【図１２】
0006025013.tif 000021
【図１３】
0006025013.tif 000022
【図１４】
0006025013.tif 000023
【図１５】
0006025013.tif 000024
【図１６】
0006025013.tif 000025