JP7069993B2 - 光スペクトル線幅演算方法、装置およびプログラム - Google Patents

Description

本発明は、自己遅延ヘテロダイン／ホモダイン法を用いた光スペクトル線幅演算方法、装置およびプログラムに関する。

コヒーレント光通信や光周波数領域反射測定等において、レーザのコヒーレンス性は性能を左右する重要な要素である。これら応用技術の高性能化に向けて狭線幅（高コヒーレンス）レーザの開発が進められており、正確なレーザ性能評価のためには高精度な光スペクトル線幅測定法が要求される。

線幅測定に広く用いられている技術に自己遅延ヘテロダイン／ホモダイン法がある（例えば非特許文献１）。図１は、自己遅延ヘテロダイン／ホモダイン法の一般的な光スペクトル線幅測定装置の構成例を示す図である。光スペクトル線幅測定装置１００は、マッハツェンダ干渉計１１０と、マッハツェンダ干渉計１１０から出射された光干渉信号を受光する受光器１２０と、受光器１２０から出力された信号を処理するＲＦスペクトル解析装置１３０とを備える。マッハツェンダ干渉計１１０は、一方のアーム導波路上には光遅延ファイバ１１１が設けられ、他方のアーム導波路上には光周波数シフタ１１２が設けられている。

本測定法では、マッハツェンダ干渉計１１０を用いて、被測定レーザ２００から出射された光から互いに遅延差τを持つ２つの光を生成し、それらを合波することによって得られる光干渉信号のパワースペクトル形状から光スペクトル線幅σを求める。ヘテロダイン法では、遅延差τを持つ光の一方に任意の光周波数シフトｆ_bを与えて干渉させ、ｆ_bを中心周波数とする光干渉信号のスペクトルを得る。これに対しホモダイン法では、光周波数シフタ１１２を省いて、光周波数シフトを与えずに同一の周波数で遅延差τを持つ光を干渉させて０Ｈｚを中心とする光干渉信号のスペクトルを得る。

光干渉信号のパワースペクトルＳ（ｆ）は次式のように記述できる。

ここでτは干渉計で与えられる遅延差、τ_cはレーザのコヒーレンス時間、δ（ｆ）はデルタ関数である。ｆ_bは図１の光周波数シフタで与えられる周波数シフトであり、ホモダイン法ではｆ_b＝０、ヘテロダイン法ではｆ_b≠０である。干渉計で与えられる遅延差τがコヒーレンス時間τ_cに対して十分長い場合（τ＞＞τ_c）、式（１）は次式のように近似できる。

ここでσはレーザのスペクトル線幅であり、発振スペクトルがローレンツ型のレーザにおいてσ＝１／（πτ_c）が成り立つことを用いた。式（２）に示されるように、τ＞＞τ_cにおいて光干渉信号のパワースペクトルは半値全幅２σのローレンツ関数となるため、パワースペクトルの半値半幅からレーザのスペクトル線幅σを求めることができる。

T.Okoshi et al., "Novel method for high resolution measurement of laser output spectrum," Electronics Letters,Vol. 16, No. 16, 1980, pp.630-631

前述のように、従来の自己遅延ヘテロダイン／ホモダイン法ではレーザのコヒーレンス時間に対して十分長い遅延差を干渉計に与える必要があるため、評価可能なレーザはτ＞＞τ_cを満たすレーザに限られるという制約がある。例えば線幅１ｋＨｚ以下のレーザを評価するためには長さ数１００ｋｍよりも十分長い遅延ファイバを用意する必要があり、必ずしも容易なことではない。

本発明は上記事情を鑑みてなされたものであり、その目的は、遅延線路を有する光干渉計により生成された光干渉信号から、遅延線路による遅延時間よりも長い遅延時間を持つ光干渉信号の位相に基づき被測定レーザのスペクトル線幅を算出できる光スペクトル線幅演算方法、装置およびプログラムを提供することにある。

上記の課題を解決するために、本発明の一態様は、被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、を有することを特徴とする。

本発明の別の一態様は、被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与え、かつ前記２分岐した光のいずれかの光周波数をシフトさせることで前記２分岐した光の間に光周波数差ｆ_bを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、
前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、を有することを特徴とする。

本発明の別の一態様は、被測定レーザから出射された光の周波数を掃引速度γで周波数掃引し、前記周波数掃引した光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、を有することを特徴とする。

本発明の別の態様では、前記位相Ｘ_N（ｔ）を算出するステップおよびパワースペクトルを算出するステップは、前記パワースペクトルがローレンツ関数に収束するまでＮを増加させながら繰り返し実施し、前記被測定レーザのスペクトル線幅を算出するステップは、ローレンツ関数に収束した前記パワースペクトルの半値半幅を前記被測定レーザのスペクトル線幅とすることを特徴とする。

本発明の別の一態様は、被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、を備えたことを特徴とする。

本発明の別の一態様は、被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与え、かつ前記２分岐した光のいずれかの光周波数をシフトさせることで前記２分岐した光の間に光周波数差ｆ_bを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、前記パワースペクトルの半値半幅から被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、を備えたことを特徴とする。

本発明の別の一態様は、被測定レーザから出射された光の周波数を掃引速度γで周波数掃引し、前記周波数掃引した光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、

前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、前記パワースペクトルの半値半幅から被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、を備えたことを特徴とする。

本発明の別の態様では、前記パワースペクトルがローレンツ関数に収束するまでＮを増加させて前記第２の位相演算部に前記位相Ｘ_N（ｔ）を算出させ、前記パワースペクトルがローレンツ関数に収束したとき、前記パワースペクトルを前記スペクトル線幅演算部に入力する比較演算部をさらに備えたことを特徴とする。

本発明の別の一態様は、プログラムであって、プロセッサを光スペクトル線幅演算装置として各種演算部として機能させることを特徴とする。

本発明を用いることにより、τ＜τ_cの条件下でも線幅が測定可能であるため、従来よりも高分解能な線幅測定が実現できる。また、従来は線幅分解能向上のためには干渉計の遅延ファイバを延長する必要があったのに対し、本発明は光干渉信号に対する信号処理のみで実施されるものであるため、装置構成を変更することなく分解能向上が可能である。したがって、自己遅延ヘテロダイン／ホモダイン法を用いた既存の線幅測定器に対して、本発明のアルゴリズムを信号処理部に組み込むだけで高分解能化できる。

自己遅延ヘテロダイン／ホモダイン法の一般的な光スペクトル線幅測定装置の構成例を示す図である。 式（５）の計算を視覚的に示すイメージ図である。 ｃｏｓＸ_N（ｔ）のパワースペクトル｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²の演算結果を示すイメージ図である。 パワースペクトルの半値半幅のＮ依存性の一例を示す図である。 本発明の第１の実施形態に係る光スペクトル線幅測定装置の構成例を示す図である。 本発明の光スペクトル線幅測定装置の演算処理装置で行う処理の流れを示すフローチャートである。 本発明の光スペクトル線幅測定装置の演算処理装置の機能構成図である。 本発明の第２の実施形態に係る光スペクトル線幅測定装置の構成例を示す図である。

以下、本発明の実施の形態について、詳細に説明する。

上記目的を達成するために、本発明では自己遅延ヘテロダイン／ホモダイン法において得られる光干渉信号を用いて、干渉計で与えられる遅延量を整数倍増加させた場合の光干渉信号に対応する位相関数を数値計算により導出し、導出した位相関数を位相とする正弦波信号のパワースペクトルの半値半幅から被測定レーザのスペクトル線幅を求める。

自己遅延ヘテロダイン／ホモダイン法で得られる光干渉信号Ｉ（ｔ）は次式のように記述できる。

ここでＰ₀は光強度、θ（ｔ）は位相雑音である。Ｉ（ｔ）の位相成分Ｘ₁（ｔ）は次式により求められる。

ここでＨ［Ｉ（ｔ）］はＩ（ｔ）のヒルベルト変換である。なお、ここでは位相がπ／２遅れた光干渉信号をヒルベルト変換により算出しているが、代わりに光９０度ハイブリッド回路を用いて取得してもよい。Ｘ₁（ｔ）を用いて次式の計算を行うことで、干渉計の遅延差τをＮ倍（Ｎは自然数）に増加させた場合の光干渉信号の位相Ｘ_N（ｔ）を算出できる。

図２は、式（５）の計算を視覚的に示すイメージ図である。τ間隔で遅延を長くした位相Ｘ₁（ｔ－ｎτ）を足し合せることにより、θ（ｔ－τ）、θ（ｔ－２τ）、・・・、θ（ｔ－（Ｎ－１）τ）が隣接する遅延差τの項の間で互いに打ち消され、結果的に位相雑音項はθ（ｔ）－θ（ｔ－Ｎτ）のみが残留する。パワースペクトルの半値半幅はｔに依存する項のみで決まる。そのため、式（４）と式（５）との比較からｔに依存しない項（Ｎ（Ｎ－１）πｆ_bτ）を無視すると、式（５）によって算出される位相Ｘ_N（ｔ）は、遅延差Ｎτを持つ干渉計で得られる位相とパワースペクトルの半値半幅が等しくなる。

図３は、ｃｏｓＸ_N（ｔ）のパワースペクトル｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²の演算結果を示すイメージ図である。また図４は、パワースペクトルの半値半幅のＮ依存性の一例を示す図である。ここでＦ［ｃｏｓＸ_N（ｔ）］はｃｏｓＸ_N（ｔ）のフーリエ変換である。

図３及び図４が示すように、Ｎを増加させるにしたがって、パワースペクトル｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²はレーザのスペクトル線幅σを半値半幅とするローレンツ関数に収束する。したがって、自己遅延ヘテロダイン／ホモダイン法において、マッハツェンダ干渉計が与える遅延差τがたとえτ＜τ_cの条件にあっても、Ｎτ＞＞τ_cとみなせる十分大きな値をＮに用いてｃｏｓＸ_N（ｔ）のパワースペクトル｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²を計算してパワースペクトルの半値半幅を求め、レーザのスペクトル線幅σを求めることができる。

（第１の実施形態）
図５は、本発明の第１の実施形態に係る光スペクトル線幅測定装置の構成例を示す図である。光スペクトル線幅測定装置５００は、マッハツェンダ干渉計５１０と、マッハツェンダ干渉計５１０から出射された光干渉信号を受光する受光器５２０と、受光器５２０から出力されたアナログ電気信号をデジタル電気信号に変換するＡ／Ｄ変換器５３０と、デジタル電気信号を処理する演算処理装置５４０とを備える。

マッハツェンダ干渉計５１０は、一方のアーム導波路上には光遅延ファイバ５１１が設けられ、他方のアーム導波路上には光周波数シフタ５１２が設けられている。尚、ホモダイン法を用いる場合は、光周波数シフタ５１２を省くものとする。従来の自己遅延ヘテロダイン／ホモダイン法と同様に、マッハツェンダ干渉計５１０を用いて、被測定レーザ２００から出射された光から互いに遅延差τを持つ２つの光を生成し、それらを合波して光干渉信号を生成する。

マッハツェンダ干渉計５１０から出射された光干渉信号は、受光器５２０およびＡ／Ｄ変換器５３０を介してデジタル信号として演算処理装置５４０に入力され、以下に説明するように処理される。

図６に、本発明の光スペクトル線幅測定装置の演算処理装置で行う処理の流れを示すフローチャートを示す。初めに、受光器５２０から光干渉信号Ｉ（ｔ）を取得し（Ｓ６０１）、式（４）により位相Ｘ₁（ｔ）を求める（Ｓ６０２）。次に、マッハツェンダ干渉計５１０で与えられる遅延差τをＮ倍増加させた場合の位相に対応する位相Ｘ_N（ｔ）を、Ｘ₁（ｔ）を用いて数値計算により求める（Ｓ６０３）。Ｘ_N（ｔ）の算出には、自己遅延ヘテロダイン法の場合は式（６）を用い、自己遅延ホモダイン法の場合は式（７）を用いる。

ここでｉは、ステップＳ６０３と後述するステップＳ６０４の実行済み回数を示し、Ｎ_iは、Ｎ₀が２以上の任意の整数でＮ_i＞Ｎ_i-1とする。和をとる位相Ｘ₁（ｔ）は、測定して記録された値を遡って
次にｃｏｓＸ_N0（ｔ）のフーリエ変換Ｆ［ｃｏｓＸ_N0（ｔ）］を計算し、パワースペクトル｜Ｆ［ｃｏｓＸ_N0（ｔ）］｜²を求める（Ｓ６０４）。次に、Ｎ₀よりも大きい整数Ｎ₁における位相Ｘ_N1（ｔ）を計算し、パワースペクトル｜Ｆ［ｃｏｓＸ_N1（ｔ）］｜²を求める（Ｓ６０５）。このとき、Ｎ₀がＮ₀＞＞τ_c／τを満たす場合、｜Ｆ［ｃｏｓＸ_N0（ｔ）］｜²と｜Ｆ［ｃｏｓＸ_N1（ｔ）］｜²はどちらもレーザ線幅σを半値半幅とするほぼ同一のローレンツ関数になる。

一方、２つのパワースペクトル｜Ｆ［ｃｏｓＸ_N0（ｔ）］｜²と｜Ｆ［ｃｏｓＸ_N1（ｔ）］｜²とが互いに異なるスペクトル形状になる場合、Ｎ₀はＮ₀＞＞τ_c／τを満たしていないことを意味し、より大きな遅延差が必要であるため、Ｎ₀よりも大きなＮ_iを用いてパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²を得る必要がある。

その場合、ｉを増加させてＮ_iを大きくしながら、Ｎ＝Ｎ_i、Ｎ_i+1について位相Ｘ_N（ｔ）及びパワースペクトル｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²を計算し、パワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²と｜Ｆ［ｃｏｓＸ_Ni+1（ｔ）］｜²とが、差異が十分小さく等しいとみなせるローレンツ関数として得られるまで繰り返す（Ｓ６０６）。Ｎの増加に対して｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²の形状変化が収束した時点で、｜Ｆ［ｃｏｓＸ_N（ｔ）］｜²の半値半幅から被測定レーザのスペクトル線幅σを求める（Ｓ６０７）。

図７に、本発明の光スペクトル線幅測定装置の演算処理装置の機能構成図を示す。図６のフローチャートで説明したように、演算処理装置５４０に光干渉信号Ｉ（ｔ）が入力されると、第１の位相演算部５４１において光干渉信号Ｉ（ｔ）から位相Ｘ₁（ｔ）を算出する。次に、第２の位相演算部５４２において、位相Ｘ₁（ｔ）に基づき２以上の任意の整数Ｎ_iで位相Ｘ_Ni（ｔ）を算出する。位相Ｘ_Ni（ｔ）が算出されると、パワースペクトル演算部５４３において位相Ｘ_Ni（ｔ）の正弦波のパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²を算出する。比較演算部５４４は、比較判定部５４５と記憶部５４６を備えており、比較判定部５４５が入力されたパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²を記憶部５４６に記憶し、パワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²と｜Ｆ［ｃｏｓＸ_Ni+1（ｔ）］｜²とを比較してそれらが所定の基準において同一のローレンツ関数であるとみなせるか否かを判定する。比較判定部５４５は、それらが同一であるとみなせる場合、そのパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²を出力し、それらが同一であるとみなせない場合、第２の位相演算部５４２に対してより大きいＮ_iに関して位相Ｘ_Ni（ｔ）を算出させる。比較判定部５４５からパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²が出力されると、スペクトル線幅演算部５４７においてパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²の半値半幅から被測定レーザのスペクトル線幅σを算出する。

上述のように、遅延差の倍率を変化させて複数の遅延差Ｎ_iτに関して光干渉信号のパワースペクトルを得ることは、線幅測定精度を高める上で重要な意味を持つ。従来の自己遅延ヘテロダイン／ホモダイン法では、レーザのコヒーレンス時間τ_cよりも十分長い遅延差τを干渉計に与える必要がある。しかし、被測定レーザのコヒーレンス時間τ_cは多くの場面において未知であり、干渉計の遅延差τがτ＞＞τ_cを満たしているか確認する手段がなかった。そのため、従来技術では、τ＞＞τ_cを保証できず、測定された線幅の値が十分な信頼を得られない場合があり、狭線幅レーザの仕様書では「線幅○○Ｈｚ以下」と記載し、具体的な線幅の値を記載できない場合がある。

一方、本発明では干渉計の遅延差τの任意の整数倍の遅延差Ｎ_iτについて光干渉信号のパワースペクトルを生成できるため、複数の異なる遅延差の光干渉信号のパワースペクトル形状を取得し、それらの変化状況を把握することができる。これにより、複数の遅延差Ｎ_iτ、Ｎ_i+1τの光干渉信号のパワースペクトルが同じローレンツ型に収束することを確認することで、Ｎ_iτ＞＞τ_cを満足しているか判別することができる。本発明では、確実にＮ_iτ＞＞τ_cを満たすパワースペクトル｜Ｆ［ｃｏｓＸ_Ni（ｔ）］｜²の半値半幅からレーザのスペクトル線幅を算出するため、従来よりも高精度な線幅評価が実現できる。

尚、被測定レーザのコヒーレンス時間τ_cが予め分かっている場合は、Ｎτ＞＞τ_cの条件を満たす倍率Ｎも予め分かるため、図６に示すフローチャートにおけるＳ６０５、Ｓ６０６のステップ、および図７に示す演算処理装置の機能構成図の比較演算部５４４は省くことができる。

（第２の実施形態）
図８は、本発明の第２の実施形態に係る光スペクトル線幅測定装置の構成例を示す図である。光スペクトル線幅測定装置８００は、光周波数掃引装置８１０と、光周波数掃引装置８１０で周波数掃引された光が入射されるマッハツェンダ干渉計８２０と、マッハツェンダ干渉計８２０から出射された光干渉信号を受光する受光器８３０と、受光器８３０から出力されたアナログ電気信号をデジタル電気信号に変換するＡ／Ｄ変換器８４０と、デジタル電気信号を処理する演算処理装置８５０とを備える。マッハツェンダ干渉計８２０は、一方のアーム導波路上に光遅延ファイバ８２１が設けられている。

被測定レーザ２００から出射された光は、光周波数掃引装置８１０に入射され、周波数掃引速度γで周波数掃引される。周波数掃引した光がマッハツェンダ干渉計８２０に入射されると２分岐され、遅延時間差τを与えられて合波され、光干渉信号が出射される。マッハツェンダ干渉計８２０から出射された光干渉信号は、受光器８３０で検出され、得られる光干渉信号Ｉ（ｔ）は次式のように記述できる。

式（３）と式（８）との比較から明らかなように、本実施形態で得られる光干渉信号は、自己遅延ヘテロダイン法でｆ_b＝γτの場合に得られる光干渉信号と等価である。そのため、式（６）にｆ_b＝γτを代入することにより、本実施形態においても第１の実施形態と同様の演算処理が適用できる。つまり、Ｉ（ｔ）の位相Ｘ₁（ｔ）を用いて次式を計算することにより、干渉計で与えられる遅延差をＮ倍増加させたＸ_N（ｔ）を求めることができる。

Ａ／Ｄ変換器８４０でデジタル信号に変換したＩ（ｔ）を演算処理装置８５０に転送し、演算処理装置８５０において、第１の実施形態と同様に図６に示す演算処理を施すことにより、被測定レーザ２００のスペクトル線幅σを求める。演算処理装置８５０の機能構成も、図７に示すものと同様の構成とすることができる。

本実施形態は、光周波数領域反射測定に用いるレーザの性能評価において特に有用である。光周波数領域反射測定では、レーザのコヒーレンス性だけでなく、光周波数掃引の非線形性も測定分解能に影響する。本実施形態においては、レーザの位相雑音だけでなく掃引非線形性の影響もスペクトル線幅に影響されるため（掃引非線形性の影響も式（８）中のθ（ｔ）に含まれる）、光周波数領域反射測定の測定性能を見積もる目的では本実施形態を用いることが望ましい。

尚、実施形態１、２の演算処理装置５４０、８５０は、コンピュータとそのコンピュータによって図６に示すステップを実行させるためのプログラムによっても実現できる。このプログラムは、記録媒体に記録して提供することも、ネットワークを通して提供することも可能である。

１００、５００、８００ 光スペクトル線幅測定装置
１１０、５１０、８２０ マッハツェンダ干渉計
１１１、５１１、８２１ 光遅延ファイバ
１１２、５１２ 光周波数シフタ
１２０、５２０、８３０ 受光器
１３０ ＲＦスペクトル解析装置
２００ 被測定レーザ
５３０、８４０ Ａ／Ｄ変換器
５４０、８５０ 演算処理装置
５４１ 第１の位相演算部
５４２ 第２の位相演算部
５４３ パワースペクトル演算部
５４４ 比較演算部
５４５ 比較判定部
５４６ 記憶部
５４７ スペクトル線幅演算部
８１０ 光周波数掃引装置

Claims (9)

1. 被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、
   前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、
   を有することを特徴とする光スペクトル線幅演算方法。
2. 被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与え、かつ前記２分岐した光のいずれかの光周波数をシフトさせることで前記２分岐した光の間に光周波数差ｆ_bを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、
   前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、
   を有することを特徴とする光スペクトル線幅演算方法。
3. 被測定レーザから出射された光の周波数を掃引速度γで周波数掃引し、前記周波数掃引した光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算方法であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出するステップと、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出するステップと、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するステップと、
   前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するステップと、
   を有することを特徴とする光スペクトル線幅演算方法。
4. 前記位相Ｘ_N（ｔ）を算出するステップおよびパワースペクトルを算出するステップは、前記パワースペクトルがローレンツ関数に収束するまでＮを増加させながら繰り返し実施し、
   前記被測定レーザのスペクトル線幅を算出するステップは、ローレンツ関数に収束した前記パワースペクトルの半値半幅を前記被測定レーザのスペクトル線幅とすることを特徴とする請求項１乃至３のいずれか一項に記載の光スペクトル線幅演算方法。
5. 被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、
   前記パワースペクトルの半値半幅から前記被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、
   を備えたことを特徴とする光スペクトル線幅演算装置。
6. 被測定レーザから出射された光を２分岐し、２分岐した光の間に遅延時間差τを与え、かつ前記２分岐した光のいずれかの光周波数をシフトさせることで前記２分岐した光の間に光周波数差ｆ_bを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、
   前記パワースペクトルの半値半幅から被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、
   を備えたことを特徴とする光スペクトル線幅演算装置。
7. 被測定レーザから出射された光の周波数を掃引速度γで周波数掃引し、前記周波数掃引した光を２分岐し、２分岐した光の間に遅延時間差τを与えて合波することによって得られる光干渉信号から前記被測定レーザの光スペクトル線幅を算出する光スペクトル線幅演算装置であって、
   前記光干渉信号の位相Ｘ₁（ｔ）を算出する第１の位相演算部と、
   前記位相Ｘ₁（ｔ）から、次式により位相Ｘ_N（ｔ）（Ｎは２以上の整数）を算出する第２の位相演算部と、
   

   

   前記位相Ｘ_N（ｔ）を位相とする正弦波のパワースペクトルを算出するパワースペクトル演算部と、
   前記パワースペクトルの半値半幅から被測定レーザのスペクトル線幅を算出するスペクトル線幅演算部と、
   を備えたことを特徴とする光スペクトル線幅演算装置。
8. 前記パワースペクトルがローレンツ関数に収束するまでＮを増加させて前記第２の位相演算部に前記位相Ｘ_N（ｔ）を算出させ、前記パワースペクトルがローレンツ関数に収束したとき、前記パワースペクトルを前記スペクトル線幅演算部に入力する比較演算部をさらに備えたことを特徴とする請求項５乃至７のいずれか一項に記載の光スペクトル線幅演算装置。
9. プロセッサを請求項５乃至８のいずれか一項に記載の光スペクトル線幅演算装置として各種演算部として機能させることを特徴とするプログラム。

Images