(57)【要約】【課題】大掛かりな測定系の構築を必要とせず、簡便に光送信機に搭載される光変調器の動作条件を推定する。【解決手段】光送信機２０に搭載された光変調器２２の仮定の動作条件である仮定動作パラメータを設定入力し、光変調器２２を搭載した光送信機２０からの被測定光のスペクトラムを測定するとともに、設定された仮定動作パラメータに基づく光スペクトラムを算出する。そして、測定した被測定光のスペクトラムと各仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を算出し、このうち最小の誤差値となる仮定動作パラメータに基づく光スペクトラムを判別抽出し、この抽出した光スペクトラムと被測定光の光スペクトラムのピーク値を合せた状態のスペクトラム波形とこれらの２乗平均誤差値とを表示部１７に表示する。【選択図】図１
2011013328.tif 000002
【特許請求の範囲】
【請求項１】
光送信機（２０）に搭載される光変調器（２２）の動作条件推定方法であって、
前記光変調器から出力される被測定光のスペクトラムを測定し得られた測定値を記憶する記憶ステップと、
前記光変調器の仮定の動作条件となる仮定動作パラメータと、該仮定動作パラメータにおける光スペクトラムの算出開始時の値となる前記仮定動作パラメータの初期値と、算出終了時の値となる前記仮定動作パラメータの最終値と、前記初期値から前記最終値に達するまでの間に、光スペクトラムを算出する毎に前記仮定動作パラメータの現在値に加算するステップ値と、をそれぞれ設定する入力ステップと、
前記光変調器からの出力光のスペクトラムを、前記設定された仮定動作パラメータに基づき算出するとともに、前記記憶した被測定光のスペクトラムの測定値と、前記設定された仮定動作パラメータに基づき算出した光スペクトラムの計算値との２乗平均誤差値を算出する算出処理ステップと、
前記仮定動作パラメータの現在値が前記最終値であるか否かを判別し、前記現在値が前記最終値であった場合に、前記算出処理ステップで算出した２乗平均誤差値が最小となる前記仮定動作パラメータを判別し抽出する判別ステップと、
を含むことを特徴とする光変調器の動作条件推定方法。
【請求項２】
前記判別ステップで抽出された前記仮定動作パラメータに基づき算出した光スペクトラムの計算値と、前記記憶した被測定光のスペクトラムの測定値とのピーク値を合せたスペクトラム波形と、前記２乗平均誤差値の最小値とを表示部（１７）に表示制御する表示制御ステップを含むことを特徴とする請求項１記載の光変調器の動作条件推定方法。
【請求項３】
前記算出処理ステップは、前記光変調器（２２）から出力される被測定光の光スペクトラムの測定値と、前記設定された仮定動作パラメータに基づき算出する光スペクトラムの計算値とのピーク値を合せた状態で記憶部（１３）に記憶するステップを含むことを特徴とする請求項１又は２記載の光変調器の動作条件推定方法。
【請求項４】
前記判別ステップは、前記仮定動作パラメータの現在値が前記最終値でなかった場合に、前記現在値に前記ステップ値を加算処理して再度前記算出処理ステップで前記ステップ値が加算された前記仮定動作パラメータに基づく光スペクトラムを算出するためのステップ値加算信号を出力するステップを含むことを特徴とする請求項１〜３の何れかに記載の光変調器の動作条件推定方法。
【請求項５】
光変調器（２２）を搭載した光送信機（２０）から出射される被測定光の光スペクトラムを測定する光スペクトラムアナライザ（１）であって、
前記被測定光のスペクトラムを測定する測定部（１１）と、
前記光変調器の仮定の動作条件となる仮定動作パラメータと、該仮定動作パラメータにおける光スペクトラムの算出開始時の値となる前記仮定動作パラメータの初期値と、算出終了時の値となる前記仮定動作パラメータの最終値と、前記初期値から前記最終値に達するまでの間に、光スペクトラムを算出する毎に、前記仮定動作パラメータの現在値に加算するステップ値と、をそれぞれ設定する入力部（１２）と、
前記光送信機から出射される光のスペクトラムを前記設定された仮定動作パラメータに基づき算出するとともに、前記測定部で測定した前記被測定光のスペクトラムと、前記設定された仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を算出する算出処理部（１４）と、
前記仮定動作パラメータの現在値が前記最終値であるか否かを判別し、前記現在値が前記最終値であった場合に、前記算出処理部で算出した２乗平均誤差値が最小となる前記仮定動作パラメータを判別し抽出する判別部（１５）と、
前記判別部で抽出された前記仮定動作パラメータに基づき算出した光スペクトラムと、前記測定部で測定した被測定光のスペクトラムとのピーク値を合せたスペクトラム波形と、前記２乗平均誤差値の最小値とを表示部（１７）に表示制御する表示制御部（１６）と、
を備えたことを特徴とする光スペクトラムアナライザ。
【請求項６】
前記算出処理部（１４）は、前記光変調器（２２）から出力される被測定光のスペクトラムと、前記設定された仮定動作パラメータに基づき算出した光スペクトラムとのピーク値を合せた状態で記憶部（１３）に記憶することを特徴とする請求項５記載の光スペクトラムアナライザ。
【請求項７】
前記判別部（１５）は、前記仮定動作パラメータの現在値が前記最終値でなかった場合に、前記現在値に前記ステップ値を加算処理して再度前記ステップ値が加算された前記仮定動作パラメータに基づく光スペクトラムを算出するためのステップ値加算信号を前記算出処理部（１４）に出力することを特徴とする請求項５又は６記載の光スペクトラムアナライザ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光通信用光送信機に搭載される光変調器おける動作条件の推定方法、及びこの推定方法を備えた光スペクトラムアナライザに関するものである。
【背景技術】
【０００２】
近年、光通信におけるビットレートは高速化しており、このような高速光通信用光送信機（以下、「光送信機」と称する）では、光通信用光源としての半導体レーザダイオードモジュール（以下、「ＬＤモジュール」と称する）へ供給する電流のＯＮ／ＯＦＦによる変調ではなく、ＬＤモジュールの注入電流を一定とし、例えばニオブ酸リチウム（ＬｉＮｂＯ3 ）結晶を用いた外部光変調器（以下、「光変調器」と称する）、光変調器ドライバ、パルスパターン発生器などを用いてＬＤモジュールからの出射光を変調している。
【０００３】
この種の光送信機から出力された被測定光のスペクトラム測定をする場合、従来から多く用いられている波長分解能が５０ｐｍ（１. ５５μｍ帯の光周波数に換算して６．２ＧＨｚ）程度の低波長分解能による光スペクトラムアナライザでは、例えば高速光通信で使用される１０ＧＨｚ程度で強度変調された信号を測定しても、図４に示すような単純な単峰性のスペクトラム波形（図中実線）しか表示することができなかった。そこで、変調によって発生する側帯波の強度や発生する周波数などの光スペクトラムの詳細な情報をより正確に得るため、現在では波長分解能５ｐｍ（１. ５５μｍ帯の光周波数に換算して０．６２ＧＨｚ）以下に設定できる高波長分解能による光スペクトラムアナライザを用いることで、同図において点線で示すような被測定光の詳細なスペクトラム波形の表示が可能となっている。
【０００４】
また、光送信機の開発から製造における調整に至るまでの各段階において、光変調器の特性、光変調器の駆動電圧やバイアス電圧、光送信機としての信号光波長などの光変調器の動作条件となる下記４種類の動作パラメータの測定や調整をすることが必要となる。
【０００５】
以下、光変調器の動作条件となる４種類の動作パラメータについて示す。
・αパラメータ
光変調器の波長チャーピングの量を表すαパラメータを測定する場合、例えば非特許文献１に開示されるようにマッハ・ツェンダー干渉計と受光器とオシロスコープを用いた方法により、光送信機の開発段階において、採用した光変調器を光送信機に組み込み、実際の信号を与えた状態でαパラメータが設計値内に収まっているかを観測する方法や、特許文献１に開示されるような方法を用いている。
・駆動電圧偏差εA 
光変調器に印可する駆動電圧の駆動電圧偏差εA を測定する場合は、光送信機の製造における調整段階において、受光器とオシロスコープを用いてアイ波形を観測する方法を用いている。また、電圧調整の際には、誤り測定器を用いて誤り率の最も少なくなるように調整している。
・バイアス電圧偏差εB 
光変調器に印可するバイアス電圧のバイアス電圧偏差εB を測定する場合は、光送信機の製造における調整段階において、受光器とオシロスコープを用いてアイ波形を観測する方法を用いている。また、電圧調整の際には、誤り測定器を用いて誤り率の最も少なくなるように調整している。
・発振波長λC 
光送信機の製造における調整段階において、信号光波長が仕様に入るよう、ＬＤモジュールの発振波長を光スペクトラムアナライザを用いて観測する方法を用いている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３８６６０８２号
【非特許文献】
【０００７】
【非特許文献１】Ｒ．Ａ．Ｓａｕｎｄｅｒｓ，ｅｔ．ａｌ．，“Ｗｉｄｅｂａｎｄ  ｃｈｉｒｐ  ｍｅａｓｕｒｅｍｅｎｔ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ｈｉｇｈ  ｂｉｔ  ｒａｔｅ  ｓｏｕｒｃｅｓ”，Ｅｌｅｃｔｒｏｎｉｃｓ  Ｌｅｔｔｅｒｓ，Ｖｏｌ．３０，Ｎｏ．１６，ｐ．１３３６
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ところで、近年では光通信装置の大量生産とコストダウン要求により、時間を掛けて正確に測定するのではなく、多少精度が低くても高速、且つ安価な構成で測定できることが求められている。しかしながら、上述した４項目の動作パラメータを測定する場合、図５に示すような大掛かりな測定系を測定する度に構築する必要があり、測定系の構築の手間や測定時間が長くなるばかりでなく、多額の費用が掛かるという問題があった。
【０００９】
また、被測定光のパラメータの一つであるαパラメータを測定する方法が特許文献１に開示されているが、この方法では変調用信号源として正弦波発生器を用いるいう特殊な動作状態におけるαパラメータを測定しているため、光送信機が通信に用いられる際の動作状態におけるαパラメータの測定を行うにはパルスパターン発生器と使い分ける必要があった。
【００１０】
そこで、本発明は上記問題点に鑑みてなされたものであり、大掛かりな測定系の構築を必要とせず、簡便に光変調器を搭載した光送信機の動作条件を推定することのできる光変調器における動作条件推定方法及び光スペクトラムアナライザを提供することを目的とするものである。
【課題を解決するための手段】
【００１１】
上記した目的を達成するために、請求項１記載の光変調器の動作条件推定方法は、光送信機２０に搭載される光変調器２２の動作条件推定方法であって、
前記光変調器から出力される被測定光のスペクトラムを測定し得られた測定値を記憶する記憶ステップと、
前記光変調器の仮定の動作条件となる仮定動作パラメータと、該仮定動作パラメータにおける光スペクトラムの算出開始時の値となる前記仮定動作パラメータの初期値と、算出終了時の値となる前記仮定動作パラメータの最終値と、前記初期値から前記最終値に達するまでの間に、光スペクトラムを算出する毎に前記仮定動作パラメータの現在値に加算するステップ値と、をそれぞれ設定する入力ステップと、
前記光変調器からの出力光のスペクトラムを、前記設定された仮定動作パラメータに基づき算出するとともに、前記記憶した被測定光のスペクトラムの測定値と、前記設定された仮定動作パラメータに基づき算出した光スペクトラムの計算値との２乗平均誤差値を算出する算出処理ステップと、
前記仮定動作パラメータの現在値が前記最終値であるか否かを判別し、前記現在値が前記最終値であった場合に、前記算出処理ステップで算出した２乗平均誤差値が最小となる前記仮定動作パラメータを判別し抽出する判別ステップと、
を含むことを特徴とする。
【００１２】
請求項２記載の光変調器の動作条件推定方法は、請求項１記載の光変調器の動作条件推定方法において、前記判別ステップで抽出された前記仮定動作パラメータに基づき算出した光スペクトラムの計算値と、前記記憶した被測定光のスペクトラムの測定値とのピーク値を合せたスペクトラム波形と、前記２乗平均誤差値の最小値とを表示部１７に表示制御する表示制御ステップを含むことを特徴とする。
【００１３】
請求項３記載の光変調器の動作条件推定方法は、請求項１又は２記載の光変調器の動作条件推定方法において、前記算出処理ステップは、前記光変調器２２から出力される被測定光の光スペクトラムの測定値と、前記設定された仮定動作パラメータに基づき算出する光スペクトラムの計算値とのピーク値を合せた状態で記憶部１３に記憶するステップを含むことを特徴とする。
【００１４】
請求項４記載の光変調器の動作条件推定方法は、請求項１〜３の何れかに記載の光変調器の動作条件推定方法前記判別ステップは、前記仮定動作パラメータの現在値が前記最終値でなかった場合に、前記現在値に前記ステップ値を加算処理して再度前記算出処理ステップで前記ステップ値が加算された前記仮定動作パラメータに基づく光スペクトラムを算出するためのステップ値加算信号を出力するステップを含むことを特徴とする。
【００１５】
請求項５記載の光スペクトラムアナライザは、光変調器２２を搭載した光送信機２０から出射される被測定光の光スペクトラムを測定する光スペクトラムアナライザ１であって、
前記被測定光のスペクトラムを測定する測定部１１と、
前記光変調器の仮定の動作条件となる仮定動作パラメータと、該仮定動作パラメータにおける光スペクトラムの算出開始時の値となる前記仮定動作パラメータの初期値と、算出終了時の値となる前記仮定動作パラメータの最終値と、前記初期値から前記最終値に達するまでの間に、光スペクトラムを算出する毎に、前記仮定動作パラメータの現在値に加算するステップ値と、をそれぞれ設定する入力部１２と、
前記光送信機から出射される光のスペクトラムを前記設定された仮定動作パラメータに基づき算出するとともに、前記測定部で測定した前記被測定光のスペクトラムと、前記設定された仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を算出する算出処理部１４と、
前記仮定動作パラメータの現在値が前記最終値であるか否かを判別し、前記現在値が前記最終値であった場合に、前記算出処理部で算出した２乗平均誤差値が最小となる前記仮定動作パラメータを判別し抽出する判別部１５と、
前記判別部で抽出された前記仮定動作パラメータに基づき算出した光スペクトラムと、前記測定部で測定した被測定光のスペクトラムとのピーク値を合せたスペクトラム波形と、前記２乗平均誤差値の最小値とを表示部１７に表示制御する表示制御部１６と、
を備えたことを特徴とする。
【００１６】
請求項６記載の光スペクトラムアナライザは、請求項５記載の光スペクトラムアナライザにおいて、前記算出処理部１４は、前記光変調器２２から出力される被測定光のスペクトラムと、前記設定された仮定動作パラメータに基づき算出した光スペクトラムとのピーク値を合せた状態で記憶部１３に記憶することを特徴とする。
【００１７】
請求項７記載の光スペクトラムアナライザは、請求項５又は６記載の光スペクトラムアナライザにおいて、前記判別部１５は、前記仮定動作パラメータの現在値が前記最終値でなかった場合に、前記現在値に前記ステップ値を加算処理して再度前記ステップ値が加算された前記仮定動作パラメータに基づく光スペクトラムを算出するためのステップ値加算信号を前記算出処理部１４に出力することを特徴とする。
【発明の効果】
【００１８】
本発明によれば、光変調器の仮定動作パラメータを設定入力するだけで、その都度各仮定仮定動作パラメータを得るための測定系を構築する必要がなく、高速で、且つ簡易的に測定対象となる光送信機に搭載された光変調器における動作条件を推定することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明に係る光スペクトラムアナライザの装置構成を示す概略ブロック図である。
【図２】同装置における光スペクトラムの表示例を示す概念図である。
【図３】同装置における光変調器動作条件推定時の処理動作を説明するためのフローチャート図である。
【図４】光スペクトラムアナライザの波長分解能の違いによるスペクトラム波形を示す説明図である。
【図５】光変調器の動作パラメータを従来の方法で測定する際に使用する装置構成例を示す概略説明図である。
【発明を実施するための形態】
【００２０】
以下、本発明を実施するための形態について、添付した図面を参照しながら詳細に説明する。
【００２１】
［本発明の実施形態］
＜装置のシステム構成＞
まず、本発明に係る光スペクトラムアナライザの構成について、図１、２を参照しながら説明する。
【００２２】
本例の光スペクトラムアナライザ１は、光通信用光源としての光通信用半導体レーザダイオードモジュール２１（以下、「ＬＤモジュール２１」と称する）への注入電流を一定とし、ニオブ酸リチウム（ＬｉＮｂＯ3 ）結晶を用いた外部光変調器２２（以下、「光変調器２２」と称する）、光変調器２２用の光変調器ドライバ２３、所望のパルスパターン（擬似ランダムビット系のＰｓｅｕｄｏ−Ｒａｎｄｏｍ  Ｂｉｔ  Ｓｅｑｕｅｎｃｅ：ＰＲＢＳなど）を発生させるパルスパターン発生器２４を備えた光送信機２０と光ファイバ３０を介して接続され、光変調器２２によって変調された被測定光のスペクトラムを測定するとともに、ユーザが任意に設定する光変調器２２の仮定の動作条件である仮定動作パラメータの設定入力に基づき、被測定光のスペクトラムから光変調器２２の動作条件を推定している。
【００２３】
図１に示すように、本例の光スペクトラムアナライザ１は、測定部１１、入力部１２、記憶部１３、算出処理部１４、判別部１５、表示制御部１６、表示部１７とを備えて構成される。
【００２４】
測定部１１は、従来より公知の光スペクトラムアナライザにおける分光手段、光電変換手段、Ａ／Ｄ変換手段の機能（例えば特許第２８９２６７０号を参照のこと）を備えており、光送信機２０からの被測定光を分光、光電変換処理して電気信号に変換し、さらにこの電気信号をＡ／Ｄ処理したディジタル信号を算出処理部１４に出力している。
【００２５】
具体的構成としては、光ファイバ３０を介して光送信機２０から出力された被測定光を入力する光入力端子と、回折格子，プリズム，干渉フィルタなどの分光素子を利用した分光によって入力した被測定光から特定波長を抽出して出力する分光手段と、分光手段から出力される抽出光を受光してこれを電気信号に変換するフォトダイオード（ＰＤ）などの光電変換素子からなる光検出手段と、光検出手段で検出された微小電気信号を増幅する増幅手段、増幅手段から出力されるアナログ信号を量子化してディジタル信号に変換するＡ／Ｄ（アナログ／ディジタル）変換手段とを備えている。
【００２６】
入力部１２は、例えば各種操作ボタンやテンキー、キーボードなどの入力機器で構成され、使用する光送信機２０の装置構成や仕様に基づきユーザが任意に設定する光変調器２２の仮定の動作条件である仮定動作パラメータの入力時に操作される。なお、入力された仮定動作パラメータは、記憶部１３に記憶される。
【００２７】
ここで、入力される仮定動作パラメータについて説明する。仮定動作パラメータとは、測定対象となる光変調器２２を搭載した光送信機２０からの被測定光のスペクトラムを測定することで推定可能となる光変調器２２の仮定の動作条件（数値）を、光送信機２０の装置構成や変調方式（例えば、Ｎｏｎ−Ｒｅｔｕｒｎ  ｔｏ  Ｚｅｒｏ, ＯＮ−Ｏｆｆ  Ｋｅｙｉｎｇ、以下「ＮＲＺ−ＯＯＫ」と称する）に基づき、ユーザが任意に設定するものであり、下記（１）〜（５）がある。
（１）ビットレートＢ（例：Ｂ＝１０Ｇｂｉｔ／ｓ）
（２）光変調器２２の駆動回路の帯域Ｂe （例：Ｂe ＝８ＧＨｚ）
（３）光変調器２２の振幅オフセット値εA （例：εA ＝０）
（４）光変調器２２のバイアスオフセット値εB （例：εB ＝０）
（５）光変調器２２のαパラメータ（例：α＝０．８）
【００２８】
ここで、推定対象となる仮定動作パラメータについては、（１）は固定値とし、（２）〜（５）について光スペクトラム算出時における算出開始時の値となる初期値、算出終了時の値となる最終値、初期値から最終値に達するまでの間で光スペクトラムの算出に使用した現在のパラメータ値（以下、「現在値」と称する）に対して算出する毎に加算するステップ値を、ユーザが知得したい光変調器２２の動作条件に応じて適宜入力する。
【００２９】
数値入力例としては、εA の場合であれば初期値＝−０．１、ステップ値＝０．０１、最終値＝０．１とする。この推定対象となる仮定動作パラメータの数は１つでも複数でも良く、通常Ｂe は既知であり、上記εA 、εB 、αパラメータが推定対象となるが、試験内容によっては帯域Ｂe も推定対象となり得る。また、初期値及び最終値が同値である場合は、ステップ値を「０」として設定することで、算出時は初期値で設定した値を常時使用することになる。
【００３０】
記憶部１３は、入力部１２で設定入力された仮定動作パラメータ、算出処理部１４が算出処理時に用いる計算式等の処理プログラム、測定部１１で測定された被測定光のスペクトラム、及び算出処理部１４で算出された仮定動作パラメータに基づく光スペクトラムや、被測定光のスペクトラムと設定された仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値、判別部１５で判別された最小の２乗平均誤差値となる仮定動作パラメータの判別結果の他、光スペクトラムアナライザを構成する各部の駆動制御に係る各種制御情報を記憶している。
【００３１】
算出処理部１４は、スペクトラム算出処理手段１４ａ、２乗平均誤差値算出処理手段１４ｂとで構成され、入力部１２で設定した仮定動作パラメータに基づく光スペクトラムの算出、被測定光のスペクトラムと設定された仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を算出し、算出した光スペクトラムや２乗平均誤差値を記憶部１３に記憶している。
【００３２】
なお、算出処理部１４における各手段が算出処理時に用いる計算式等の処理プログラムは、設定入力された仮定動作パラメータ、変調方式、装置構成等に応じた算出処理が行えるよう予め記憶部１３に記憶されており、算出処理内容に応じて処理プログラムを適宜選択して算出処理を行っている。
【００３３】
スペクトラム算出処理手段１４ａは、光送信機２０からの被測定光のスペクトラムを、入力部１２から設定入力された仮定動作パラメータに基づき算出している。また、スペクトラム算出処理手段１４ａは、後述する仮定動作パラメータ判別手段１５ａからステップ値加算信号を入力すると、現在値にステップ値を加算処理した仮定動作パラメータに基づく光スペクトラムを再度算出している。さらに、スペクトラム算出処理手段１４ａは、被測定光のスペクトラムと仮定動作パラメータに基づき算出した各光スペクトラムのピーク値をそれぞれ合せた状態で記憶部１３に記憶している。
【００３４】
なお、光スペクトラムの算出に際し、光スペクトラムアナライザ１の波長分解能λRES の値も必要であるが、これは測定時における光スペクトラムアナライザ１の設定内容（装置仕様）から求めることができる。
【００３５】
２乗平均誤差値算出処理手段１４ｂは、測定部１１で測定した被測定光のスペクトラムと、スペクトラム算出処理手段１４ａにおいて、仮定動作パラメータを初期値から最終値までステップ値毎に変化させながらそれぞれ算出した各光スペクトラムとの２乗平均誤差値をそれぞれ算出し、この算出した２乗平均誤差値を記憶部１３に記憶している。
【００３６】
判別部１５は、仮定動作パラメータ判別手段１５ａと、２乗平均誤差値判別手段１５ｂとで構成され、スペクトラム算出処理手段１４ａでステップ値毎に変化させている仮定動作パラメータが最終値であるか否かの判別や、記憶部１３に記憶された被測定光のスペクトラムと仮定動作パラメータに基づく光スペクトラムとの２乗平均誤差値が最小となる仮定動作パラメータを判別抽出している。
【００３７】
仮定動作パラメータ判別手段１５ａは、入力部１２で設定された仮定動作パラメータにおける現在値がその仮定動作パラメータの最終値であるか否かを判別している。このとき、最終値であった場合は、全ての仮定動作パラメータについての算出処理が完了したと判別し、被測定光のスペクトラムと、仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値の判別を行うための判別開始信号を２乗平均誤差値判別手段１５ｂに出力している。また、仮定動作パラメータの現在値が最終値でなかった場合は、現在値にステップ値を加算処理した仮定動作パラメータに基づく光スペクトラムを算出するためのステップ値加算信号をスペクトラム算出処理手段１４ａに出力している。
【００３８】
２乗平均誤差値判別手段１５ｂは、仮定動作パラメータ判別手段１５ａから判別開始信号を入力すると、記憶部１３に記憶された被測定光のスペクトラムと仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を比較して最小の誤差値となる仮定動作パラメータに基づき算出した光スペクトラムを判別抽出し、この判別結果を表示制御部１６に出力している。なお、この判別結果は、記憶部１３に記憶される。
【００３９】
表示制御部１６は、仮定動作パラメータを入力部１２から設定入力する際の入力画面表示、光送信機２０からの被測定光のスペクトラムや仮定動作パラメータに基づく光スペクトラムの表示や各光スペクトラムのピーク値を合せた重ね表示、算出処理部１４で算出した各種算出結果など、光スペクトラムアナライザ１を駆動する上で必要な表示内容を制御するための表示制御信号を表示部１７に出力している。また、表示制御部１６は、２乗平均誤差値判別手段１５ｂから判別結果を入力すると、２乗平均誤差値の最小値と、記憶部１３に記憶された被測定光のスペクトラムと、前記２乗平均誤差値が最小となる仮定動作パラメータに基づき算出した光スペクトラムとのピーク値を合せたスペクトラム波形を表示するための表示制御信号を出力している。
【００４０】
表示部１７は、例えば液晶ディスプレイ（ＬＣＤ）などの表示機器で構成され、図２に示すような被測定光のスペクトラムと仮定動作パラメータに基づき算出した光スペクトラムとの重ね表示や２乗平均誤差値など、表示制御部１６からの表示制御信号に基づく表示内容を表示している。
【００４１】
＜装置の処理動作＞
次に、上述した光スペクトラムアナライザ１における、仮定動作パラメータ設定入力から仮定動作パラメータの推定に至るまでの一連の処理動作について、図３を参照しながら説明する。
【００４２】
図３に示すように、まず測定対象となる光変調器２２を搭載した光送信機２０からの被測定光のスペクトラムを測定し、この光スペクトラムを表示する（ＳＴ１）。次に、ユーザは、光変調器２２の仮定動作パラメータ（初期値、最終値、ステップ値）を入力部１２から設定入力する（ＳＴ２）。
【００４３】
仮定動作パラメータの設定入力が終わると、この仮定動作パラメータに基づく光スペクトラムを初期値から順に算出し（ＳＴ３）、算出した光スペクトラムと被測定光のスペクトラムとのピーク値に合せた状態で記憶部１３に記憶する（ＳＴ４）。そして、被測定光のスペクトラムと仮定動作パラメータに基づく光スペクトラムの２乗平均誤差値を算出し（ＳＴ５）、記憶部１３に記憶した後（ＳＴ６）、仮定動作パラメータの現在値が最終値であるか否かを判別する（ＳＴ７）。
【００４４】
このとき、仮定動作パラメータの現在値が最終値である場合は（ＳＴ７−Ｙｅｓ）、仮定動作パラメータ判別手段１５ａから２乗平均誤差値判別手段１５ｂに判別開始信号を出力し、測定した被測定光のスペクトラムと仮定動作パラメータに基づき算出した光スペクトラムの２乗平均誤差値のうち、どの仮定動作パラメータに基づき算出した光スペクトラムとの誤差値が最小となる光スペクトラムであるかを判別し抽出する（ＳＴ８）。
そして、誤差値が最小となる仮定動作パラメータに基づき算出した光スペクトラムが抽出されると、この光スペクトラムと被測定光のスペクトラムのピーク値を合せた状態のスペクトラム波形とこれらの２乗平均誤差値とを表示部１７に表示して（ＳＴ９）、処理を終了する。
【００４５】
一方、仮定動作パラメータの現在値が最終値でない場合は（ＳＴ７−Ｎｏ）、仮定動作パラメータの現在値にステップ値を加算処理した仮定動作パラメータに基づく光スペクトラムを算出するため、仮定動作パラメータ判別手段１５ａからスペクトラム算出処理手段１４ａにステップ値加算信号を出力し（ＳＴ１０）、再度ＳＴ３へ戻る。
【００４６】
なお、上記処理動作では、設定入力された仮定動作パラメータのうちの１つパラメータにおける処理動作であるが、通常は設定した仮定動作パラメータが複数であるため、それぞれの仮定動作パラメータの組み合わせに応じてＳＴ３〜ＳＴ７までの処理内容を入れ子構造として処理を行う。
【００４７】
［実施例］
以下、本発明に係る光スペクトラムアナライザ１における具体的な処理動作例（図３に示す処理動作におけるＳＴ３，ＳＴ４に相当）について説明する。
なお、下記実施例は本発明を限定するものではなく、前・後記の趣旨に照らし合わせて設計変更することはいずれも本発明の技術的範囲に含まれるものである。
【００４８】
上述した光スペクトラムアナライザ１による具体的な実施例について説明する。ここでは、光送信機２０に対しパルスパターン発生器２４から供給される信号として擬似ランダムビット系のＰＲＢＳを使用し、光送信機２０からの被測定光の変調方式としてＮＲＺ−ＯＯＫを使用し、光変調器２２動作条件として上述した（１）〜（５）を変調方式に応じた任意の値を入力した場合について説明する。
【００４９】
まず、デ−タ信号をａ［ｎ］，ｎ＝０，１，・・・，Ｎ−１とおく。ａ［ｎ］の値は０または１で、Ｎはデ−タ数を表し、ここではＰＲＢＳのデータ列を入れる。次に１ビット内を時間的に分割する数Ｒを適宜定める。この値は時間分解能に相当する。ＮＲ個のメモリ領域ｂ［ｍ］, ｍ＝０，１，・・・，（ＮＲ−１）を用意し、下記数１で定める値を格納する。
【００５０】
【数１】
2011013328.tif 000003
【００５１】
光変調器ドライバ２３を帯域幅Ｂe のフィルタとみたて、信号ｂ［ｍ］をこのフィルタを通過させて得られる信号をｃ［ｍ］とおく。具体的には、ｂ［ｍ］を離散Ｆｏｕｒｉｅｒ変換（以下、「ＤＦＴ」と称する）した後、Ｂe とビットレートＢ，およびフィルタの周波数特性（例えば５次ベッセルフィルタなど）に対応する関数形から定まる数値を掛け算し、その結果を逆ＤＦＴすることでｃ［ｍ］が得られる。
【００５２】
光変調器２２からの出力電場を表す信号ｅ［ｍ］＝ｘ［ｍ］＋ｊｙ［ｍ］を、下記数２〜数６の一連の式を用いて算出する。電場ｅ［ｍ］は複素信号であり、ｘ［ｍ］、ｙ［ｍ］はその実部と虚部を、ｊは虚数単位をそれぞれ表す。
【００５３】
【数２】
2011013328.tif 000004
【００５４】
【数３】
2011013328.tif 000005
【００５５】
【数４】
2011013328.tif 000006
【００５６】
【数５】
2011013328.tif 000007
【００５７】
【数６】
2011013328.tif 000008
【００５８】
なお、数５のεA は光変調器２２において最大振幅を与える駆動電圧を１とした場合のオフセット値を、εB は５０％透過時の駆動電圧をバイアス電圧とした場合のこれに対するオフセット値を、数６のαは光変調器２２のαパラメータ値を表す。
【００５９】
以上で、ＮＲＺ−ＯＯＫの場合のｅ［ｍ］の算出までが終了する。ＲＺ−ＯＯＫ信号の場合は数１におけるｂ［・］の定め方を、またＯＯＫ以外の変調方式の場合には、数３〜数６を変調方式に応じて変更すればよい。
【００６０】
下記数７は、求めたｅ［ｍ］より光スペクトラムＩ［ｍ］を求める方法であり、各変調方式に共通である。信号ｅ［ｍ］をＤＦＴして電場の周波数成分をＥ［ｍ］を求め、Ｉ［ｍ］を下記数７で算出する。
【００６１】
【数７】
2011013328.tif 000009
【００６２】
そして、数７で得られた光スペクトラムに対し、光スペクトラムアナライザ１の波長分解能に応じた丸め処理を行う。丸め処理の一実施例を説明する。デ−タ数Ｋを適宜定め、フィルタ特性Ｇ［ｋ］を下記数８で算出する。
【００６３】
【数８】
2011013328.tif 000010
【００６４】
但し、Ｃは規格化定数、ＷはλRES に応じて定める数である。またＫは数８でＧ［Ｋ］が十分小さくなるように定める。
【００６５】
次に、丸め処理した光スペクトラムＩ’［ｍ］を数９で算出する。
【００６６】
【数９】
2011013328.tif 000011
【００６７】
数９の右辺で、Ｇ［・］の引数は添字ｋの絶対値｜ｋ｜であり、和の中のＩ［・］の添字ｍ＋ｋが負となる場合はＮＲを足し、ＮＲ以上となる場合はＮＲを引くことで０以上ＮＲ−１以下となるようにする。算出された光スペクトラムＩ’［ｍ］のピークを、被測定光のスペクトラムのピーク値に合わせて再配置し、記憶部１３に記憶する。そして、この算出した値を用いて被測定光のスペクトラムと仮定動作パラメータに基づく光スペクトラムの２乗平均誤差値を算出している。
【００６８】
以上説明したように、上述した光スペクトラムアナライザ１は、光変調器２２の仮定の動作条件である仮定動作パラメータを設定入力し、光変調器２２を搭載した光送信機２０からの被測定光のスペクトラムを、設定された仮定動作パラメータに基づき算出する。そして、測定した被測定光のスペクトラムと各仮定動作パラメータに基づき算出した光スペクトラムとの２乗平均誤差値を算出し、このうち最小の誤差値となる仮定動作パラメータに基づく光スペクトラムを判別抽出し、この抽出した光スペクトラムと被測定光のスペクトラムのピーク値を合せた状態のスペクトラム波形とこれらの２乗平均誤差値とを表示部１７に表示する。
【００６９】
これにより、光変調器２２の仮定動作パラメータを設定入力するだけで、高速、且つ簡易的に測定対象となる光送信機２０に搭載された光変調器２２における動作条件を推定することができる。また、装置構成として光スペクトラムアナライザ１のみで推定することができるため、仮定動作パラメータを得るためにその都度測定系を構築する必要がなく測定作業が簡素化し、また安価に光変調器２２の動作条件の推定を行うことができる。
【００７０】
ところで、上述した形態において、変調方式としてＮｏｎ−Ｒｅｔｕｒｎ  ｔｏ  Ｚｅｒｏ, ＯＮ−Ｏｆｆ  Ｋｅｙｉｎｇを用いた例としたが、Ｒｅｔｕｒｎ  ｔｏ  Ｚｅｒｏ, Ｄｉｆｆｅｒｅｎｔｉａｌ  Ｐｈａｓｅ  Ｓｈｉｆｔ  Ｋｅｙｉｎｇなど他の変調方式の場合でも同様に光変調器２２の動作条件推定を行うことができる。
【００７１】
また、上述した光変調器２２の動作条件推定に係る各構成は、光スペクトラムアナライザ１内部にプログラムを有する他、光スペクトラムアナライザ１に接続されたパソコンに上述した処理プログラムを保有し、測定データを転送して算出することもできる。
【符号の説明】
【００７２】
１…光スペクトラムアナライザ
１１…測定部
１２…入力部
１３…記憶部
１４…算出処理部（１４ａ…スペクトラム算出処理手段、１４ｂ…２乗平均誤差値算出処理手段）
１５…判別部（１５ａ…仮定動作パラメータ判別手段、１５ｂ…２乗平均誤差値判別手段）
１６…表示制御部
１７…表示部
２０…光送信機
２１…半導体レーザダイオードモジュール（ＬＤモジュール）
２２…光変調器
２３…光変調器ドライバ
２４…パルスパターン発生器
【図１】
2011013328.tif 000012
【図２】
2011013328.tif 000013
【図３】
2011013328.tif 000014
【図４】
2011013328.tif 000015
【図５】
2011013328.tif 000016