(57)【要約】【課題】低コストおよび短時間で伝播装置の特性を評価する。【解決手段】信号を伝播する伝播装置の特性を評価する評価装置であって、伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、分離部が分離したジッタ成分に基づいて、伝播装置の特性を評価する評価部と、を備える評価装置を提供する。【選択図】図１
2010164554.tif 000002
【特許請求の範囲】
【請求項１】
信号を伝播する伝播装置の特性を評価する評価装置であって、
前記伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、
前記出力信号のジッタの確率密度関数から、前記出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、
前記分離部が分離した前記ジッタ成分に基づいて、前記伝播装置の特性を評価する評価部と、
を備える評価装置。
【請求項２】
前記評価部は、確定成分が基準値以下であることを条件として、前記伝播装置が良品であると判定する
請求項１に記載の評価装置。
【請求項３】
前記伝播装置に与えられる入力信号のジッタの確率密度分布を表わす確率密度関数を測定する入力信号測定部を更に備え、
前記分離部は、更に、前記入力信号のジッタの確率密度関数から、前記入力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離し、
前記評価部は、前記分離部が分離した前記入力信号に含まれるジッタ成分に対する、前記分離部が分離した前記出力信号に含まれるジッタ成分の減衰率に基づいて、前記伝播装置の特性を評価する
請求項１または２に記載の評価装置。
【請求項４】
前記評価部は、前記入力信号のジッタに含まれる確定成分に対する、前記出力信号のジッタに含まれる確定成分の減衰率に基づいて、前記伝播装置の特性を評価する
請求項３に記載の評価装置。
【請求項５】
前記評価部は、更に、前記入力信号のジッタに含まれるランダム成分に対する前記出力信号のジッタに含まれるランダム成分の減衰率に基づいて、前記伝播装置の特性を評価する
請求項４に記載の評価装置。
【請求項６】
前記分離部は、
前記確率密度関数の標準偏差を算出する標準偏差算出部と、
前記確率密度関数のスペクトルを算出するスペクトル算出部と、
前記スペクトルのヌル周波数を検出するヌル周波数検出部と、
前記ヌル周波数検出部が検出した前記ヌル周波数に基づいて、予め定められた複数種類の確定成分のそれぞれに対して、スペクトルの理論値を算出する理論値算出部と、
前記標準偏差算出部が算出した前記標準偏差、および、前記スペクトル算出部が算出した前記スペクトルに基づいて、前記確率密度関数に含まれる確定成分に対して、スペクトルの測定値を算出する測定値算出部と、
前記理論値算出部が前記確定成分の種類毎に算出した前記理論値のうち、前記測定値算出部が算出した前記測定値に最も近い前記理論値に対応する前記確定成分の種類を、前記確率密度関数に含まれる確定成分の種類と判定するモデル判定部と、
前記モデル判定部が判定した前記確定成分の種類と、前記ヌル周波数検出部が検出した前記ヌル周波数とに基づいて、前記確率密度関数に含まれる前記確定成分を算出する確定成分算出部と
を有する
請求項１から５の何れかに記載の評価装置。
【請求項７】
前記分離部は、前記確定成分算出部が算出した前記確定成分を、前記確率密度関数から除去することで、前記確率密度関数に含まれるランダム成分を算出するランダム成分算出部を更に有する
請求項６に記載の評価装置。
【請求項８】
前記評価部は、前記出力信号のジッタに含まれるランダム成分および確定成分に基づいて、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを計算して、前記伝播装置の特性を評価する
請求項１から７の何れかに記載の評価装置。
【請求項９】
前記評価部は、
ビットエラーレートまたはサンプルタイミングの何れか一方の指定を受ける指定部と、
前記出力信号のジッタに含まれるランダム成分の標準偏差および前記出力信号のジッタに含まれる確定成分のピークツゥピーク値をパラメータとして代入することにより得られる、ランダム成分および所定の確率密度分布の確定成分を含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を用いて、指定された前記ビットエラーレート以下となるサンプルタイミングの範囲、または、指定された前記サンプルタイミングにおけるビットエラーレートを計算する計算部と
を有する
請求項８に記載の評価装置。
【請求項１０】
前記伝播装置は、周波数特性が外部から変更可能であり、
当該評価装置は、前記伝播装置の特性の評価結果に基づいて、前記伝播装置の周波数特性を調整する調整部を更に備える
請求項１から９の何れかに記載の評価装置。
【請求項１１】
信号を伝播する伝播装置の特性を評価する評価方法であって、
前記伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定ステップと、
前記出力信号のジッタの確率密度関数から、前記出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離ステップと、
前記分離ステップにおいて分離した前記ジッタ成分に基づいて、前記伝播装置の特性を評価する評価ステップと、
を備える評価方法。
【請求項１２】
コンピュータを、信号を伝播する伝播装置の特性を評価する評価装置として機能させるプログラムであって、
前記コンピュータを、
前記伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、
前記出力信号のジッタの確率密度関数から、前記出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、
前記分離部が分離した前記ジッタ成分に基づいて、前記伝播装置の特性を評価する評価部と、
して機能させるプログラム。
【請求項１３】
請求項１２に記載したプログラムを記憶した記憶媒体。
【請求項１４】
所定の信号を生成する電子デバイスであって、
前記所定の信号を生成して出力する動作回路と、
前記動作回路の特性を評価する評価装置と、
を備え、
前記評価装置は、
前記動作回路を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、
前記出力信号のジッタの確率密度関数から、前記出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、
前記分離部が分離した前記ジッタ成分に基づいて、前記動作回路の特性を評価する評価部と、
を有する電子デバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、評価装置、評価方法、プログラム、記憶媒体および電子デバイスに関する。
【背景技術】
【０００２】
従来より、イコライザが知られている。イコライザの特性を評価する場合、イコライザを通過した出力信号のアイ開口をオシロスコープまたはビットエラーレート測定器等により検出し、検出したアイ開口を基準値と比較していた。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかし、このような評価方法は、コストが高く、評価時間も長くなってしまう。従って、量産時における出荷試験等においては、イコライザの評価が行われない場合が多かった。また、このような評価方法では、受信回路等の他の回路と同一のチップ上に形成されたイコライザを評価することは困難であった。
【０００４】
そこで本発明の１つの側面においては、上記の課題を解決することのできる評価装置、評価方法、プログラム、記憶媒体および電子デバイスを提供することを目的とする。この目的は特許請求の範囲における独立項に記載の特徴の組み合わせにより達成される。また従属項は本発明の更なる有利な具体例を規定する。
【課題を解決するための手段】
【０００５】
即ち、本発明の第１の態様によると、信号を伝播する伝播装置の特性を評価する評価装置であって、伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、分離部が分離したジッタ成分に基づいて、伝播装置の特性を評価する評価部と、を備える評価装置を提供する。
【０００６】
本発明の第２の態様によると、信号を伝播する伝播装置の特性を評価する評価方法であって、伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定ステップと、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離ステップと、分離ステップにおいて分離したジッタ成分に基づいて、伝播装置の特性を評価する評価ステップと、を備える評価方法を提供する。
【０００７】
本発明の第３の態様によると、コンピュータを、信号を伝播する伝播装置の特性を評価する評価装置として機能させるプログラムであって、コンピュータを、伝播装置を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、分離部が分離したジッタ成分に基づいて、伝播装置の特性を評価する評価部と、して機能させるプログラム、および、このようなプログラムを記憶した記憶媒体を提供する。
【０００８】
本発明の第４の態様によると、所定の信号を生成する電子デバイスであって、所定の信号を生成して出力する動作回路と、動作回路の特性を評価する評価装置と、を備え、評価装置は、動作回路を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する出力信号測定部と、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する分離部と、分離部が分離したジッタ成分に基づいて、動作回路の特性を評価する評価部と、を有する電子デバイスを提供する。
【０００９】
なお上記の発明の概要は、本発明の必要な特徴の全てを列挙したものではなく、これらの特徴群のサブコンビネーションも又発明となりうる。
【図面の簡単な説明】
【００１０】
【図１】図１は、一つの実施形態に係る評価装置２０の機能構成例を示す図である。
【図２】図２は、評価装置２０による処理の概要を示すフローチャートである。
【図３】図３は、第１変形例に係る評価部３０の機能構成例を示す図である。
【図４】図４は、第２変形例に係る評価装置２０の機能構成例を示す図である。
【図５】図５は、一つの実施形態に係る算出装置１０の機能構成例を示す図である。
【図６】図６は、ガウス分布の確率密度関数の一例を示す図である。
【図７】図７は、デュアルディラック分布の確率密度関数の一例を示す図である。
【図８】図８は、ランダム成分とデュアルディラック分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。
【図９】図９は、ランダム成分とデュアルディラック分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。
【図１０】図１０は、一様分布の確率密度関数の一例を示す図である。
【図１１】図１１は、ランダム成分と一様分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。
【図１２】図１２は、ランダム成分と一様分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。
【図１３】図１３は、台形分布の確率密度関数の一例を示す図である。
【図１４】図１４は、ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。
【図１５】図１５は、ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。
【図１６】図１６は、算出装置１０により、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを算出する処理の概要を示すフローチャートである。
【図１７】図１７は、一つの実施形態に係る確定成分識別装置１００の機能構成例を示す図である。
【図１８Ａ】図１８Ａは、一様分布の確定成分の確率密度関数を示す図である。
【図１８Ｂ】図１８Ｂは、台形分布の確定成分の確率密度関数を示す図である。
【図１９Ａ】図１９Ａは、デュアルディラック分布の確定成分の確率密度関数を示す図である。
【図１９Ｂ】図１９Ｂは、シングルディラック分布の確定成分の確率密度関数を示す図である。
【図２０Ａ】図２０Ａは、サイン波分布の確定成分のモデルについて、確定成分の確率密度関数と、そのスペクトルとを示す図である。
【図２０Ｂ】図２０Ｂは、一様分布の確定成分のモデルについて、確定成分の確率密度関数と、そのスペクトルとを示す図である。
【図２１】図２１は、スペクトル算出部１２０が算出するスペクトルの一例を示す図である。
【図２２】図２２は、確定成分の種類毎の、時間ドメインでのモデル式、周波数ドメインでのモデル式、第１ヌル周波数ｆｚｅｒｏとピークツゥピーク値ＤＪＰ−Ｐとの関係、および、ピークツゥピーク値ＤＪＰ-Ｐと二乗平均値ＤＪＲＭＳとの関係を示す表である。
【図２３】図２３は、理論値算出部１４０が算出した、各種類の確定成分のスペクトルＹ（ｘ）の理論値と、測定値算出部１５０が算出した、確定成分のスペクトルＹ（ｘ）の測定値との一例を示す図である。
【図２４】図２４は、ヌル周波数検出部１３０の動作例を示す図である。
【図２５】図２５は、確定成分識別装置１００の他の構成例を示す図である。
【図２６】図２６は、確定成分識別装置１００の動作の概要を示すフローチャートである。
【図２７】図２７は、確定成分識別装置１００の他の構成例を示す図である。
【図２８】図２８は、図２７に関連して説明した確定成分識別装置１００の動作例を説明するフローチャートである。
【図２９】図２９は、確定成分識別装置１００の他の構成例を示す図である。
【図３０】図３０は、図２９に示した確定成分識別装置１００の動作例を説明するフローチャートである。
【図３１】図３１は、一つの実施形態に係る電子デバイス４００の構成例を示す図である。
【図３２】図３２は、一つの実施形態に係るコンピュータ１９００のハードウェア構成の一例を示す。
【発明を実施するための形態】
【００１１】
以下、発明の実施の形態を通じて本発明の（一）側面を説明するが、以下の実施形態は特許請求の範囲にかかる発明を限定するものではなく、又実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
【００１２】
図１は、一つの実施形態に係る評価装置２０の機能構成例を示す図である。評価装置２０は、信号を伝播する伝播装置４０の特性を評価する。より具体的には、評価装置２０は、伝播装置４０のジッタ特性を評価することができる。
【００１３】
伝播装置４０は、一例として、信号のジッタに含まれる確定成分を低減するイコライザであってよい。また、伝播装置４０は、入力した信号の特定の周波数成分を低減するフィルタ等であってもよい。また、伝播装置４０は、電気信号に限らず、光信号等を伝播する装置等であってもよい。
【００１４】
信号発生部２２は、伝播装置４０に信号を与える。信号発生部２２は、一例として、予めジッタを含んだ擬似ランダム信号等を与えてよい。また、信号発生部２２は、伝播装置４０がイコライザである場合、予め定められた分布の確定ジッタを含む信号を、伝播装置４０に与えてもよい。これにより、信号発生部２２は、伝播装置４０が良品であるか否かを確実に判定させることができる。
【００１５】
入力信号測定部２４は、信号発生部２２から伝播装置４０に与えられる入力信号のジッタの確率密度分布を表わす確率密度関数を測定する。出力信号測定部２６は、伝播装置４０を通過した出力信号のジッタの確率密度分布を表わす確率密度関数を測定する。
【００１６】
入力信号測定部２４および出力信号測定部２６は、一例として、ジッタを複数回測定して、ジッタの確率密度分布を表す確率密度関数を測定する。入力信号測定部２４および出力信号測定部２６は、ジッタの発生パターンの繰返し周期よりも長い期間、ジッタを測定する。なお、信号発生部２２から伝播装置４０に与えられる入力信号のジッタの確率密度関数が既知であれば、入力信号測定部２４は、入力信号のジッタの確率密度関数を予め記憶していてもよい。
【００１７】
分離部２８は、入力信号のジッタの確率密度関数から、入力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する。更に、分離部２８は、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する。
【００１８】
本実施形態においては、分離部２８は、出力信号および入力信号のジッタに含まれる確定成分のピークツゥピーク値を算出する。信号のジッタに含まれる確定成分のピークツゥピーク値は、確定成分の確率密度分布におけるジッタの範囲を表わす。例えば、ジッタの確定成分が理想的なエッジタイミングを基準に正負に分布する場合には、確定成分のピークツゥピーク値は、理想的なエッジタイミングから負側へ向かうジッタの最大距離と、理想的なエッジタイミングから正側へ向かうジッタの最大距離とを加算した距離（時間幅）であってよい。
【００１９】
更に、本実施形態において、分離部２８は、出力信号または入力信号のジッタに含まれるランダム成分の標準偏差を更に算出してもよい。ジッタに含まれるランダム成分はガウス分布により表されるので、当該標準偏差は、ランダム成分の大きさを表すことができる。
【００２０】
このような分離部２８は、詳細を図１７〜図３０において説明する確定成分識別装置１００により実現される。例えば、分離部２８は、図２５に示される確定成分識別装置１００と同様の構成であってよい。なお、分離部２８は、信号発生部２２から伝播装置４０に与えられる入力信号のジッタの確率密度関数が既知であれば、入力信号のジッタに含まれる確定成分のピークツゥピーク値およびランダム成分の標準偏差を予め記憶していてもよい。
【００２１】
評価部３０は、分離部２８が分離したジッタ成分に基づいて、伝播装置４０の特性を評価する。評価部３０は、一例として、出力信号のジッタに含まれる確定成分が基準値以下であることを条件として、伝播装置４０が良品であると判定する。より具体的には、評価部３０は、出力信号のジッタに含まれる確定成分のピークツゥピーク値が、仕様により定められた値以下であることを条件として、伝播装置４０が良品であると評価する。
【００２２】
また、評価部３０は、分離部２８が分離した入力信号に含まれるジッタ成分に対する、分離部２８が分離した出力信号に含まれるジッタ成分の減衰率に基づいて、伝播装置４０の特性を評価してもよい。評価部３０は、一例として、減衰率が仕様により定められた値以下であることを条件として、伝播装置４０が良品であると判定してよい。
【００２３】
本実施形態においては、評価部３０は、入力信号のジッタに含まれる確定成分に対する、出力信号のジッタに含まれる確定成分の減衰率に基づいて、伝播装置４０の特性を評価する。評価部３０は、一例として、入力信号のジッタに含まれる確定成分のピークツゥピーク値に対する、出力信号のジッタに含まれる確定成分のピークツゥピーク値が、仕様により定められた値以下であることを条件として、伝播装置４０が良品であると判定する。
【００２４】
また、評価部３０は、更に、入力信号のジッタに含まれるランダム成分に対する出力信号のジッタに含まれるランダム成分の減衰率に基づいて、伝播装置４０の特性を評価してもよい。例えば、評価部３０は、更に、入力信号のジッタに含まれるランダム成分の標準偏差に対する、出力信号のジッタに含まれるランダム成分の標準偏差が、仕様により定められた値以下であることを条件として、伝播装置４０が良品であると判定する。
【００２５】
図２は、評価装置２０による処理の概要を示すフローチャートである。まず、信号発生部２２は、伝播装置４０に対して信号を与える（Ｓ２１）。続いて、入力信号測定部２４は、伝播装置４０に与えられる入力信号のジッタの確率密度関数を測定する（Ｓ２２）。
【００２６】
続いて、出力信号測定部２６は、伝播装置４０を通過した出力信号のジッタの確率密度関数を測定する（Ｓ２３）。なお、ステップＳ２２およびステップＳ２３の処理は、並行して実行されても、異なる期間において実行されてもよい。
【００２７】
続いて、分離部２８は、入力信号のジッタの確率密度関数から、入力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する（Ｓ２４）。続いて、分離部２８は、出力信号のジッタの確率密度関数から、出力信号のジッタに含まれるランダム成分および確定成分の少なくとも一方のジッタ成分を分離する（Ｓ２５）。続いて、評価部３０は、分離部２８が分離したジッタ成分に基づいて、伝播装置４０の特性を評価する（Ｓ２６）。
【００２８】
以上のような評価装置２０は、オシロスコープおよびビットエラーレート測定器等を用いずにジッタを測定して、伝播装置４０の特性を評価することができる。これにより、評価装置２０によれば、低いコストおよび短い測定時間で伝播装置４０の特性を評価することができる。従って、評価装置２０によれば、量産された伝播装置４０の出荷試験において、ジッタ特性の良否を判定することができる。
【００２９】
図３は、第１変形例に係る評価部３０の機能構成例を示す図である。本変形例に係る評価部３０は、指定部１６と、計算部１８と、良否判定部３４とを有する。
【００３０】
指定部１６は、ビットエラーレートまたはサンプルタイミングの何れか一方の指定を受ける。指定部１６は、一例として、ユーザ等から指定を受けてよい。
【００３１】
計算部１８は、詳細を図５〜図１６において説明する算出装置１０に備えられる計算部１８と同一構成である。即ち、計算部１８は、分離部２８から、出力信号のジッタに含まれるランダム成分の標準偏差（σＲＪ，ＲＭＳ）および出力信号のジッタに含まれる確定成分のピークツゥピーク値（ＤＪＰ−Ｐ）を、パラメータとして受け取る。さらに、計算部１８は、分離部２８から、出力信号のジッタに含まれる確定成分の確率密度分布の種類を受け取る。計算部１８は、一例として、確定成分の確率密度分布の種類として、デュアルディラック分布、一様分布または台形分布を受け取ってよい。計算部１８は、台形分布の場合には、更に、底辺比を受け取ってよい。なお、本変形例において、分離部２８は、詳細を図１７〜図３０において説明する確定成分識別装置１００により実現される。
【００３２】
そして、計算部１８は、出力信号のジッタに含まれるランダム成分の標準偏差および出力信号のジッタに含まれる確定成分のピークツゥピーク値をパラメータとして代入することにより得られる、サンプルタイミングとビットエラーレートとの関係式を取得する。この場合において、計算部１８は、ランダム成分および所定の確率密度分布の確定成分を含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を取得する。より詳しくは、計算部１８は、ランダム成分および受け取った種類の確率密度分布の確定成分を含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を取得する。なお、計算部１８のより詳細な構成および機能については、図５〜図１６において説明する。
【００３３】
良否判定部３４は、指定されたビットエラーレート以下となるサンプルタイミングの範囲が基準範囲以上であることを条件として、伝播装置４０が良品であると評価する。または、良否判定部３４は、指定されたサンプルタイミングにおけるビットエラーレートが基準レート以下であることを条件として、伝播装置４０が良品であると評価する。
【００３４】
このように本変形例に係る評価部３０は、出力信号のジッタに含まれるランダム成分および確定成分に基づいて、指定されたビットエラーレート以下となるサンプルタイミングの範囲（即ち、アイ開口）、または、指定されたサンプルタイミングにおけるビットエラーレートを計算して、伝播装置４０の特性を評価することができる。これにより、本変形例に係る評価装置２０によれば、低いコストおよび短い測定時間で、伝播装置４０のアイ開口またはビットエラーレートを測定して、伝播装置４０の良否を判定することができる。
【００３５】
図４は、第２変形例に係る評価装置２０の機能構成例を示す図である。本変形例に係る評価装置２０は、周波数特性が外部から変更できる伝播装置４０の特性を測定して、この伝播装置４０の周波数特性を調整する。周波数特性が外部から変更できる伝播装置４０は、一例として、外部からフィルタ係数が変更可能なフィルタ回路、および、周波数帯域毎のゲインを外部から変更可能なイコライザ等であってよい。
【００３６】
本変形例に係る評価装置２０は、調整部３６を更に備える。調整部３６は、伝播装置４０の特性の評価結果に基づいて、伝播装置４０の周波数特性を調整する。
【００３７】
調整部３６は、伝播装置４０の周波数特性を順次に変更しながらジッタ特性を測定して、最も特性が良い（例えば、ジッタの減衰率が最も低い、アイ開口が最も広いまたはビットエラーレートが最も低い）周波数特性を検出する。また、調整部３６は、少なくとも仕様により定められた特性値より良い周波数特性を検出してもよい。そして、調整部３６は、このような周波数特性となるように伝播装置４０を設定する。このように本変形例に係る評価装置２０によれば、伝播装置４０を、ジッタ特性の良い周波数特性に調整することができる。
【００３８】
図５は、一つの実施形態に係る算出装置１０の機能構成例を示す図である。算出装置１０は、電気回路に伝送される電気信号または光回路に伝送される光信号等の対象信号の特性を算出する。本例においては、算出装置１０は、対象信号が有するジッタの分布を表わす確率分布関数が与えられ、指定されたビットエラーレート以下となるサンプルタイミングの範囲（即ち、時間方向のアイ開口）、または、指定されたサンプルタイミングにおけるビットエラーレートを算出する。
【００３９】
ジッタは、対象信号の位相雑音であり、より具体的には、対象信号のエッジタイミングと、理想的なエッジタイミングとの時間差であってよい。確率密度関数は、ジッタを複数回測定することにより得られる、測定値の分布を示す関数であってよい。例えば、確率密度関数は、対象信号のエッジ毎のジッタを測定したときの、測定値の分布を表わす関数であってよい。
【００４０】
ここで、対象信号のジッタには、ランダムに発生するランダム成分、および、対象信号の伝送路の特性等により確定的に生じる確定成分が含まれる。ジッタに含まれるランダム成分は、ガウス分布により表わされる。ジッタに含まれる確定成分は、その発生原因等に応じた種類の分布モデルにより表わされる。確定成分の分布モデルの種類として、例えば、デュアルディラック分布、一様分布および台形分布等が存在する。
【００４１】
算出装置１０は、モデル判別部１２と、パラメータ算出部１４と、指定部１６と、計算部１８とを備える。モデル判別部１２は、対象信号のジッタの確定密度関数に基づき、対象信号のジッタに含まれる確定成分の確率密度分布の種類を判別する。
【００４２】
モデル判別部１２は、一例として、対象信号のジッタに含まれる確定成分の確率密度分布の種類として、デュアルディラック分布、一様分布または台形分布の何れかを判別してよい。さらに、モデル判別部１２は、台形分布の場合には、台形の下辺の長さに対する上への長さの比（即ち、底辺比）を更に判別してよい。なお、モデル判別部１２は、一例として、詳細を図１７〜図３０において説明する確定成分識別装置１００により実現されてよい。
【００４３】
パラメータ算出部１４は、対象信号のジッタの確定密度関数に基づき、対象信号のジッタに含まれるランダム成分の標準偏差および対象信号のジッタに含まれる確定成分のピークツゥピーク値を算出する。なお、対象信号のジッタに含まれる確定成分のピークツゥピーク値とは、確定成分の確率密度分布におけるジッタの範囲を表わす。例えば、ジッタの確定成分が理想的なエッジタイミングを基準に正負に分布する場合には、確定成分のピークツゥピーク値は、理想的なエッジタイミングから負側へ向かうジッタの最大距離と、理想的なエッジタイミングから正側へ向かうジッタの最大距離とを加算した距離（時間幅）であってよい。なお、パラメータ算出部１４は、一例として、詳細を図１７〜図３０において説明する確定成分識別装置１００により実現されてよい。
【００４４】
指定部１６は、ビットエラーレートまたはサンプルタイミングの何れか一方の指定を受ける。指定部１６は、一例として、ユーザ等から指定を受けてよい。また、指定部１６は、メモリ等に予め記憶されたビットエラーレートまたはサンプルタイミングを読み出すことにより、指定を受けてもよい。
【００４５】
計算部１８は、モデル判別部１２から対象信号のジッタに含まれる確定成分の確率密度分布の種類を受け取る。計算部１８は、一例として、確定成分の確率密度分布の種類として、デュアルディラック分布、一様分布または台形分布を受け取ってよい。計算部１８は、台形分布の場合には、更に、底辺比を受け取ってよい。更に、計算部１８は、パラメータ算出部１４から対象信号のジッタに含まれるランダム成分の標準偏差および対象信号のジッタに含まれる確定成分のピークツゥピーク値をパラメータとして受け取る。
【００４６】
更に、計算部１８は、対象信号のジッタに含まれるランダム成分の標準偏差および対象信号のジッタに含まれる確定成分のピークツゥピーク値をパラメータとして代入することにより得られる、サンプルタイミングとビットエラーレートとの関係式を取得する。この場合において、計算部１８は、ランダム成分および所定の確率密度分布の確定成分を含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を取得する。より詳しくは、計算部１８は、ランダム成分および受け取った種類の確率密度分布の確定成分を含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を取得する。なお、関係式の具体例については詳細を後述する。
【００４７】
そして、計算部１８は、取得した関係式を用いて、指定されたビットエラーレート以下となるサンプルタイミングの範囲（即ち、時間方向のアイ開口）、または、指定されたサンプルタイミングにおけるビットエラーレートを計算する。計算部１８は、ビットエラーレートが指定された場合には、取得した関係式を用いて、対象信号における指定されたビットエラーレート以下となるサンプルタイミングの範囲を計算してよい。また、計算部１８は、サンプルタイミングが指定された場合には、取得した関係式を用いて対象信号の指定されたタイミングにおけるビットエラーレートを計算してよい。このような算出装置１０によれば、指定されたビットエラーレート以下となるサンプルタイミングの範囲、および、指定されたサンプルタイミングにおけるビットエラーレートを、精度良く算出することができる。
【００４８】
なお、算出装置１０は、他の例として、対象信号の振幅雑音の確率密度分布を表す確率密度関数が与えられ、指定されたビットエラーレート以下となる振幅の範囲（即ち、振幅方向のアイ開口）、または、指定された振幅値におけるビットエラーレートを計算する構成であってもよい。このような他の例に係る算出装置１０は、次の点において各構成要素の機能が相違する。
【００４９】
即ち、モデル判別部１２は、対象信号の振幅雑音の確定密度関数に基づき、対象信号の振幅雑音に含まれる確定成分の確率密度分布の種類を判別する。振幅雑音は、対象信号の振幅と、ビット値を識別するための理想的な振幅のスレショルドとの振幅差であってよい。この場合、確率密度関数は、振幅雑音を複数回測定することにより得られる、測定値の分布を示す関数であってよい。例えば、確率密度関数は、対象信号のビット毎の振幅雑音を測定したときの、測定値の分布を表わす関数であってよい。
【００５０】
パラメータ算出部１４は、対象信号の振幅雑音の確定密度関数に基づき、対象信号の振幅雑音に含まれるランダム成分の標準偏差および対象信号の振幅雑音に含まれる確定成分のピークツゥピーク値を算出する。なお、対象信号の振幅雑音に含まれる確定成分のピークツゥピーク値とは、確定成分の確率密度分布における振幅雑音の範囲を表わす。例えば、振幅雑音の確定成分が理想的なスレッショルドを基準に正負に分布する場合には、確定成分のピークツゥピーク値は、理想的なスレッショルドから負側へ向かう振幅雑音の最大距離と、理想的なスレッショルドから正側へ向かう振幅雑音の最大距離とを加算した距離（振幅量）であってよい。
【００５１】
指定部１６は、ビットエラーレートまたは振幅値の何れか一方の指定を受ける。計算部１８は、対象信号の振幅雑音に含まれる確定成分の確率密度分布の種類、および、対象信号の振幅雑音に含まれるランダム成分の標準偏差および対象信号の振幅雑音に含まれる確定成分のピークツゥピーク値をパラメータとして受け取る。そして、計算部１８は、対象信号の振幅雑音に含まれるランダム成分の標準偏差および対象信号の振幅雑音に含まれる確定成分のピークツゥピーク値をパラメータとして代入することにより得られる、振幅値とビットエラーレートとの関係式を取得する。他の例に係る計算部１８が取得する関係式は、サンプルタイミングを表わすパラメータが、振幅値に代わっている点において異なるが、他の点については同一である。
【００５２】
そして、計算部１８は、取得した関係式を用いて、指定されたビットエラーレート以下となる振幅値の範囲（即ち、振幅方向のアイ開口）、または、指定された振幅値におけるビットエラーレートを計算する。これにより、算出装置１０によれば、指定されたビットエラーレート以下となる振幅値の範囲、および、指定された振幅値におけるビットエラーレートを、精度良く算出することができる。
【００５３】
図６は、ガウス分布の確率密度関数の一例を示す図である。図７は、デュアルディラック分布の確率密度関数の一例を示す図である。図８は、ランダム成分とデュアルディラック分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。図６〜図８は、横軸にサンプルタイミングｔ、縦軸に確率密度をそれぞれ表わし、分布の平均をサンプルタイミングｔ＝０とした確率密度関数を表わしている。
【００５４】
ジッタに含まれるランダム成分は、図６に示されるようなガウス分布となる。ランダム成分の分布は、下記式（１０１）のガウス分布の確率密度関数ｒ（ｔ）により表わされる。
【００５５】
【数１】
2010164554.tif 000003
【００５６】
なお、ｔはサンプルタイミングを表わす。また、σＲＪ，ＲＭＳは、ランダム成分の標準偏差を表わす。以降の式においても同様である。
【００５７】
また、デュアルディラック分布は、図７に示されるような確率密度関数で表わされる。即ち、デュアルディラック分布は、下記式（１０２）により表わされる。
【００５８】
【数２】
2010164554.tif 000004
【００５９】
なお、ＤＪＰ−Ｐは、確定成分のピークツゥピーク値を表わす。以降の式においても同様である。
【００６０】
また、ランダム成分とデュアルディラック分布の確定成分とを含むジッタの分布は、ガウス分布の確率密度関数とデュアルディラック分布の確率密度関数とをコンボリューションした、図８に示されるような関数により表わされる。即ち、ランダム成分とデュアルディラック分布の確定成分とを含むジッタの確率分布関数は、下記式（１０３）により表わされる。
【００６１】
【数３】
2010164554.tif 000005
【００６２】
図９は、ランダム成分とデュアルディラック分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。ランダム成分とデュアルディラック分布の確定成分とを含むジッタを有する信号における、任意のサンプルタイミング（ｔ）のビットエラーレートは、ガウス分布の確率密度関数とデュアルディラック分布の確率密度関数とをコンボリューションした関数を、ｔから∞まで積分した関数により表わされる。即ち、このようなビットエラーレート（ＢＥＲＤＵＡＬ（ｔ））は、下記式（１０４）により表わされる。
【００６３】
【数４】
2010164554.tif 000006
【００６４】
ここで、下式に示されるような誤算関数ｅｒｆ（ｘ）を定義する。ｅｒｆ（ｘ）は、以降の式においても全て同様である。
【００６５】
【数５】
2010164554.tif 000007
【００６６】
誤差関数ｅｒｆ（ｘ）を用いて、式（１０４）から、下記式（１１）を導き出すことができる。
【００６７】
【数６】
2010164554.tif 000008
【００６８】
このような式（１１）は、ランダム成分とデュアルディラック分布の確定成分とを含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を表す。従って、計算部１８は、対象信号のジッタに含まれる確定成分の確率密度分布の種類がデュアルディラック分布の場合、式（１１）に示される関係式を用いて、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを計算してよい。
【００６９】
より具体的には、対象信号のジッタに含まれる確定成分の確率密度分布の種類がデュアルディラック分布の場合であって、サンプルタイミングが指定された場合には、計算部１８は、式（１１）に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値およびサンプルタイミングを代入して、ビットエラーレートを算出してよい。一方、ビットエラーレートが指定された場合には、計算部１８は、式（１１）の逆関数の近似式に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値およびビットエラーレートを代入して、指定されたビットエラーレートとなる条件下におけるジッタの最大値を算出してよい。
【００７０】
そして、計算部１８は、対象信号のビット間隔の平均値から算出したジッタの最大値の２倍（即ち、ジッタのピークツゥピーク値）を減算することで、指定されたビットエラーレート以下となるサンプルタイミングの範囲を算出してよい。これにより、算出装置１０によれば、対象信号のジッタに含まれる確定成分の確率密度分布の種類がデュアルディラック分布の場合における、指定されたビットエラーレート以下となる振幅の範囲、または、指定された振幅値におけるビットエラーレートを精度良く算出することができる。
【００７１】
図１０は、一様分布の確率密度関数の一例を示す図である。図１１は、ランダム成分と一様分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。図１０〜図１１は、横軸にサンプルタイミングｔ、縦軸に確率密度をそれぞれ表わし、分布の平均をサンプルタイミングｔ＝０とした確率密度関数を表わしている。
【００７２】
一様分布は、図１０に示されるような確率密度関数で表わされる。即ち、一様分布は、下記式（１０６）により表わされる。
【００７３】
【数７】
2010164554.tif 000009
【００７４】
また、ランダム成分と一様分布の確定成分とを含むジッタの分布は、ガウス分布の確率密度関数と一様分布の確率密度関数とをコンボリューションした図１１に示されるような関数により表わされる。即ち、ランダム成分と一様分布の確定成分とを含むジッタの分布は、下記式（１０７）により表わされる。
【００７５】
【数８】
2010164554.tif 000010
【００７６】
図１２は、ランダム成分と一様分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。ランダム成分と一様分布の確定成分とを含むジッタを有する信号における、任意のサンプルタイミング（ｔ）のビットエラーレート（ＢＥＲＵｎｉｆｏｒｍ（ｔ））は、ガウス分布の確率密度関数と一様分布の確率密度関数とをコンボリューションした関数を、ｔから∞まで積分した関数により表わされる。即ち、このようなビットエラーレート（ＢＥＲＵｎｉｆｏｒｍ（ｔ））は、下記式（１２）により表わされる。
【００７７】
【数９】
2010164554.tif 000011
【００７８】
このような式（１２）は、ランダム成分と一様分布の確定成分とを含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を表す。従って、計算部１８は、対象信号のジッタに含まれる確定成分の確率密度分布の種類が一様分布の場合、式（１２）に示される関係式を用いて、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを計算してよい。
【００７９】
より具体的には、対象信号のジッタに含まれる確定成分の確率密度分布の種類が一様分布の場合であって、サンプルタイミングが指定された場合には、計算部１８は、式（１２）に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値およびサンプルタイミングを代入して、ビットエラーレートを算出してよい。一方、ビットエラーレートが指定された場合には、計算部１８は、式（１２）の逆関数の近似式に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値およびビットエラーレートを代入して、指定されたビットエラーレートとなる条件下におけるジッタの最大値を算出してよい。
【００８０】
そして、計算部１８は、対象信号のビット間隔の平均値から算出したジッタの最大値の２倍（即ち、ジッタのピークツゥピーク値）を減算することで、指定されたビットエラーレート以下となるサンプルタイミングの範囲を算出してよい。これにより、算出装置１０によれば、対象信号のジッタに含まれる確定成分の確率密度分布の種類が一様分布の場合における、指定されたビットエラーレート以下となる振幅の範囲、または、指定された振幅値におけるビットエラーレートを精度良く算出することができる。
【００８１】
図１３は、台形分布の確率密度関数の一例を示す図である。図１４は、ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタの確率密度関数の一例を示す図である。図１３〜図１４は、横軸にサンプルタイミングｔ、縦軸に確率密度をそれぞれ表わし、分布の平均をサンプルタイミングｔ＝０とした確率密度関数を表わしている。
【００８２】
台形分布は、図１３に示されるような確率密度関数で表わされる。即ち、台形分布は、下記式（１０８）により表わされる。なお、式（１０８）において、αは、底辺比（台形の下辺の長さに対する上への長さの比）を表わす。以降の式においても同様である。
【００８３】
【数１０】
2010164554.tif 000012
【００８４】
また、ランダム成分と台形分布の確定成分とを含むジッタの分布は、ガウス分布の確率密度関数と台形分布の確率密度関数とをコンボリューションした図１４に示されるような関数により表わされる。即ち、ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタの分布は、下記式（１０９）により表わされる。
【００８５】
【数１１】
2010164554.tif 000013
【００８６】
図１５は、ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタを有する信号のビットエラーレートの一例を示す図である。ランダム成分と底辺比がαの台形分布の確定成分とを含むジッタを有する信号における、任意のサンプルタイミング（ｔ）のビットエラーレート（ＢＥＲＴｒａｐｅｚ（ｔ））は、ガウス分布の確率密度関数と台形分布の確率密度関数とをコンボリューションした関数を、ｔから∞まで積分した関数により表わされる。即ち、このようなビットエラーレート（ＢＥＲＴｒａｐｅｚ（ｔ））は、下記式（１３）により表わされる。
【００８７】
【数１２】
2010164554.tif 000014
【００８８】
なお、式（１３）において、ｆ１（ｔ）、ｆ２（ｔ）、ｆ３（ｔ）、ｆ４（ｔ）は、下式により表される。
【００８９】
【数１３】
2010164554.tif 000015
【数１４】
2010164554.tif 000016
【数１５】
2010164554.tif 000017
【数１６】
2010164554.tif 000018
【００９０】
式（１３）は、ランダム成分と、底辺比がαの台形分布の確定成分とを含むジッタを有する信号を伝送する伝送モデルにおける、サンプルタイミングとビットエラーレートとの関係式を表す。従って、計算部１８は、対象信号のジッタに含まれる確定成分の確率密度分布の種類が台形分布の場合、式（１３）に示される関係式を用いて、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを計算してよい。
【００９１】
より具体的には、対象信号のジッタに含まれる確定成分の確率密度分布の種類が台形分布の場合であって、サンプルタイミングが指定された場合には、計算部１８は、式（１３）に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値、底辺比およびサンプルタイミングを代入して、ビットエラーレートを算出してよい。一方、ビットエラーレートが指定された場合には、計算部１８は、式（１３）の逆関数の近似式に対して、ランダム成分の標準偏差、確定成分のピークツゥピーク値、底辺比およびビットエラーレートを代入して、指定されたビットエラーレートとなる条件下におけるジッタの最大値を算出してよい。
【００９２】
そして、計算部１８は、対象信号のビット間隔の平均値から算出したジッタの最大値の２倍（即ち、ジッタのピークツゥピーク値）を減算することで、指定されたビットエラーレート以下となるサンプルタイミングの範囲を算出してよい。これにより、算出装置１０によれば、対象信号のジッタに含まれる確定成分の確率密度分布の種類が底辺比αの台形分布の場合における、指定されたビットエラーレート以下となる振幅の範囲、または、指定された振幅値におけるビットエラーレートを精度良く算出することができる。
【００９３】
なお、算出装置１０は、指定されたビットエラーレート以下となる振幅の範囲（即ち、振幅方向のアイ開口）、または、指定された振幅値におけるビットエラーレートを算出する場合には、上記の各関係式（式（１１）、式（１２）、式（１３））のサンプルタイミングを表わす変数ｔを、振幅値を表わす変数ｇに置き換えた関係式を用いる。これにより、算出装置１０によれば、対象信号のジッタに含まれる確定成分の確率密度分布の種類が、デュアルディラック分布、一様分布または底辺比αの台形分布の場合における、指定されたビットエラーレート以下となる振幅の範囲、または、指定された振幅値におけるビットエラーレートを精度良く算出することができる。
【００９４】
図１６は、算出装置１０により、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを算出する処理の概要を示すフローチャートである。まず、算出装置１０は、電気回路または光回路等に伝送される対象信号のジッタを複数回測定して、対象信号のジッタの分布を表わす確率密度関数を生成する（Ｓ１０１）。なお、他の装置から確率密度関数が与えられる場合には、算出装置１０は、当該処理を実行しなくてよい。
【００９５】
次に、算出装置１０は、対象信号のジッタの確定密度関数に基づき、対象信号のジッタに含まれる確定成分の確率密度分布の種類を判別する（Ｓ１０２）。次に、算出装置１０は、対象信号のジッタの確定密度関数に基づき、対象信号のジッタに含まれるランダム成分の標準偏差および対象信号のジッタに含まれる確定成分のピークツゥピーク値を算出する（Ｓ１０３）。この場合において、算出装置１０は、ジッタに含まれる確定成分の確率密度分布の種別が台形分布の場合には、更に、底辺比を算出してよい。
【００９６】
次に、算出装置１０は、対象信号のジッタに含まれる確定成分の確率密度分布の種別に応じた関係式にパラメータを代入して、サンプルタイミングとビットエラーレートとの関係式を取得する（Ｓ１０４）。具体的には、算出装置１０は、ランダム成分の標準偏差、確定成分のピークツゥピーク値および底辺比をパラメータとして代入した関係式を取得する。
【００９７】
次に、算出装置１０は、取得した関係式に指定された値（ビットエラーレートまたはサンプルタイミング）を代入して、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを算出する（Ｓ１０５）。このような処理により、算出装置１０によれば、指定されたビットエラーレート以下となるサンプルタイミングの範囲、または、指定されたサンプルタイミングにおけるビットエラーレートを精度良く算出することできる。
【００９８】
図１７は、一つの実施形態に係る確定成分識別装置１００の機能構成例を示す図である。確定成分識別装置１００は、与えられる確率密度関数に含まれる確定成分の分布形状を判定する。本例の確定成分識別装置１００は、特に確定成分の種類（またはモデル）を判定する。確定成分識別装置１００は、標準偏差算出部１１０、スペクトル算出部１２０、ヌル周波数検出部１３０、理論値算出部１４０、測定値算出部１５０、および、モデル判定部１６０を備える。
【００９９】
まず、確率密度関数、確定成分、および、ランダム成分について説明する。確率密度関数に含まれる確定成分およびランダム成分が互いに独立とすると、確定成分およびランダム成分の時間軸におけるコンボリューションにより、時間軸の確率密度関数が与えられる。このため、確率密度関数のスペクトルＨ（ｆ）は、下式のように、確定成分のスペクトルＤ（ｆ）およびランダム成分のスペクトルＲ（ｆ）の乗算で与えられる。
【０１００】
【数１７】
2010164554.tif 000019
【０１０１】
また、ランダム成分は、一般にガウス分布で与えられるので、そのスペクトルＲ（ｆ）は、下式で与えられる。
【０１０２】
【数１８】
2010164554.tif 000020
【０１０３】
ただし、σＲＪ，ＲＭＳは、ランダム成分の標準偏差を示す。このように、ランダム成分の標準偏差が測定できれば、ランダム成分を決定できる。そして、与えられる確率密度関数からランダム成分を減算することで、確定成分も決定することができる。
【０１０４】
しかし、与えられる確率密度関数は、確定成分およびランダム成分が合成されているので、確率密度関数からランダム成分のみの標準偏差を測定することは困難である。ここで、確定成分およびランダム成分が互いに独立とすると、確率密度関数の標準偏差σＴＪ，ＲＭＳは、ランダム成分の標準偏差σＲＪ，ＲＭＳおよび確定成分の標準偏差ＤＪＲＭＳにより、以下のように与えられる。
【０１０５】
【数１９】
2010164554.tif 000021
【０１０６】
よって確率密度関数のスペクトルは、下式で与えられる。
【０１０７】
【数２０】
2010164554.tif 000022
【０１０８】
また、確定成分の標準偏差ＤＪＲＭＳは、確定成分のピークツゥピーク値ＤＪＰ−Ｐから求めることができる。つまり、確率密度関数の標準偏差σＴＪ，ＲＭＳ、および、確定成分のピークツゥピーク値ＤＪＰ−Ｐが測定できれば、式（１）に基づいて、確率密度関数Ｈ（ｆ）から確定成分Ｄ（ｆ）を分離することができる。
【０１０９】
標準偏差算出部１１０は、与えられる確率密度関数の標準偏差σＴＪ，ＲＭＳを算出する。標準偏差算出部１１０は、式（２）に基づいて、標準偏差算出部１１０は、時間軸の確率密度関数の標準偏差σＴＪ，ＲＭＳを算出してよい。
【０１１０】
【数２１】
2010164554.tif 000023
【０１１１】
ただし、与えられる確率密度関数のｉ番目のビンの中央値をｘｉ、ｉ番目のビンのイベント数をｙｉ、総ビン数をＮとする。また、ｐｉおよびμは以下で与えられる。
【０１１２】
【数２２】
2010164554.tif 000024
【０１１３】
スペクトル算出部１２０は、与えられる確率密度関数のスペクトルを算出する。標準偏差算出部１１０およびスペクトル算出部１２０は、時間軸の確率密度関数が並列に与えられてよい。
【０１１４】
ヌル周波数検出部１３０は、スペクトル算出部１２０が算出したスペクトルにおいて、スペクトルのパワーが略零となるヌル周波数（または、スペクトルが極小値を示すヌル周波数）を検出する。後述するように、確率密度関数のスペクトルにおけるヌル周波数から、確定成分のピークツゥピーク値ＤＪＰ−Ｐを算出することができる。
【０１１５】
図１８Ａ、図１８Ｂ、図１９Ａ、図１９Ｂは、確定成分の各モデルの確率密度関数を示す図である。図１８Ａは、一様分布の確定成分を示す。図１８Ｂは、台形分布の確定成分を示す。図１９Ａは、デュアルディラック分布の確定成分を示す。図１９Ｂは、シングルディラック分布の確定成分を示す。
【０１１６】
図１８Ａから図１９Ｂに示すように、それぞれのモデルにおいて、ピークツゥピーク値ＤＪＰ−Ｐが定まれば、確定成分の分布を一義に定めることができる。ただし、台形分布の確定成分は、上辺および下辺の比が更に与えられることが好ましい。また、シングルディラック分布の確定成分は、ピークツゥピーク値が略０の確定成分として与えられる。
【０１１７】
本例の確定成分識別装置１００は、確率密度関数のスペクトルにおけるヌル周波数に基づいて、確定成分のピークツゥピーク値ＤＪＰ−Ｐを算出する。以下では、ヌル周波数のうち、最小の周波数である第１ヌル周波数を用いる例を説明する。
【０１１８】
図２０Ａ、図２０Ｂは、所定の確定成分のモデルについて、確定成分の確率密度関数と、そのスペクトルとを示す図である。図２０Ａは、サイン波分布の確定成分のモデルについて示す。また図２０Ｂは、一様分布の確定成分のモデルについて示す。図２０Ａおよび図２０Ｂにおいて、左の波形は時間領域の確率密度関数を示し、右の波形は確率密度関数のスペクトルを示す。また、時間領域の確定成分のピークツゥピーク値をＤＪＰ−Ｐとする。
【０１１９】
図２０Ａに示すように、サイン波分布の確定成分の確率密度関数をフーリエ変換したスペクトルの第１ヌル周波数は、０．７６５／ＤＪＰ−Ｐで与えられる。即ち、当該第１ヌル周波数の逆数に、係数０．７６５を乗算することにより、確定成分のピークツゥピーク値ＤＪＰ−Ｐを算出することができる。ただし、当該係数は近似値であり、より高精度または低精度に求めた係数を用いてもよい。例えば、より有効桁数の大きい、または、小さい係数を用いてよい。
【０１２０】
また、図２０Ｂに示すように、一様分布の確定成分の確率密度関数をフーリエ変換したスペクトルの第１ヌル周波数は、１／ＤＪＰ−Ｐで与えられる。即ち、当該第１ヌル周波数の逆数を求めることにより、確定成分のピークツゥピーク値ＤＪＰ−Ｐを算出することができる。
【０１２１】
同様に、台形分布、デュアルディラック分布等の他の確定成分の種類についても、第１ヌル周波数から、ピークツゥピーク値を算出することができる。ただし、図２０Ａ、図２０Ｂに示すように、確定成分のモデルにより、第１ヌル周波数とピークツゥピーク値ＤＪＰ−Ｐとの関係が異なるので、確定成分を精度よく求めるには、確定成分のモデルを判定する必要がある。
【０１２２】
図２１は、スペクトル算出部１２０が算出するスペクトルの一例を示す図である。スペクトル算出部１２０は、与えられる確率密度関数のスペクトルを算出する。スペクトル算出部１２０は、与えられる確率密度関数をフーリエ変換することで、確率密度関数のスペクトルを算出してよい。スペクトル算出部１２０が算出するスペクトルは、図２１に示すように、確定成分のスペクトルと、ランダム成分のスペクトルとを合成したスペクトルとなる。
【０１２３】
ヌル周波数検出部１３０は、スペクトル算出部１２０が算出したスペクトルのヌル周波数を検出する。本例のヌル周波数検出部１３０は、スペクトルのヌル周波数のうち、最小の周波数となる第１ヌル周波数を検出する。
【０１２４】
上述したように、確定成分のスペクトルにおける第１ヌル周波数は、確定成分の時間軸におけるピークツゥピーク値に対応する。そして、図２１に示すように、確定成分およびランダム成分を合成したスペクトルの第１ヌル周波数は、確定成分のスペクトルにおける第１ヌル周波数とほぼ一致する。
【０１２５】
このため、与えられる確率密度関数のスペクトルにおける第１ヌル周波数を検出することで、確定成分のスペクトルにおける第１ヌル周波数を検出することができる。また、上述したように、確定成分のスペクトルにおける第１ヌル周波数から、確定成分のピークツゥピーク値ＤＪＰ−Ｐを算出することができる。ただし、算出されるピークツゥピーク値ＤＪＰ−Ｐは、確定成分の種類により異なる。
【０１２６】
図２２は、確定成分の種類毎の、時間ドメインでのモデル式、周波数ドメインでのモデル式、第１ヌル周波数ｆｚｅｒｏとピークツゥピーク値ＤＪＰ−Ｐとの関係、および、ピークツゥピーク値ＤＪＰ-Ｐと二乗平均値ＤＪＲＭＳとの関係を示す表である。なお、図２２におけるＩ０は、第１種０次のベッセル関数を示す。
【０１２７】
また、図２２におけるαは、台形分布における下辺に対する上辺の比を示す。すなわち、台形分布においてα＝１の場合が一様分布に相当し、α＝０の場合が三角分布に相当する。なお、確定成分識別装置１００が取り扱う確定成分の種類は、上述した種類に限定されない。確定成分識別装置１００は、スペクトルの第１ヌル周波数から、ピークツゥピーク値を算出できる全ての確定成分のモデルを判定してよい。
【０１２８】
ここで、それぞれの確定成分のスペクトルのモデル式を、与えられる確率密度関数から測定できる標準偏差σＴＪ、ＲＭＳと、第１ヌル周波数ｆｚｅｒｏとで表現する。サイン波分布の確定成分のスペクトルのモデル式は、式（１）および図２２から、下式で与えられる。
【０１２９】
【数２３】
2010164554.tif 000025
【０１３０】
また、台形分布の確定成分のスペクトルのモデル式は下式で与えられる。
【０１３１】
【数２４】
2010164554.tif 000026
【０１３２】
また、デュアルディラック分布の確定成分のスペクトルのモデル式は下式で与えられる。
【０１３３】
【数２５】
2010164554.tif 000027
【０１３４】
また、シングルディラック分布の確定成分のスペクトルのモデル式は下式で与えられる。
【０１３５】
【数２６】
2010164554.tif 000028
【０１３６】
理論値算出部１４０は、標準偏差算出部１１０が算出した標準偏差σＴＪ、ＲＭＳと、ヌル周波数検出部１３０が検出した第１ヌル周波数ｆｚｅｒｏとに基づいて、予め定められた複数種類の確定成分のそれぞれに対して、スペクトルの理論値を算出する。例えば理論値算出部１４０は、式（３）から式（６）に基づいて、確定成分の種類毎の理論値を算出してよい。理論値算出部１４０には、式（３）から式（６）が予め与えられてよい。
【０１３７】
また、理論値算出部１４０は、標準偏差算出部１１０が算出した標準偏差σＴＪ、ＲＭＳを有するガウス分布のスペクトルで正規化された、スペクトルの理論値を算出してもよい。例えば理論値算出部１４０は、式（３）から式（６）において、
【数２７】
2010164554.tif 000029
として、Ｙ（ｘ）を算出してよい。
【０１３８】
また、理論値算出部１４０は、第１ヌル周波数ｆｚｅｒｏで周波数ｆを正規化したスペクトルの理論値を算出してもよい。例えば理論値算出部１４０は、式（３）から式（６）において、ｆ／ｆｚｅｒｏ＝ｘとして理論値を算出してよい。
【０１３９】
これらの正規化を行う場合、式（３）から式（６）は、以下のように変形される。
【０１４０】
【数２８】
2010164554.tif 000030
【数２９】
2010164554.tif 000031
【数３０】
2010164554.tif 000032
【数３１】
2010164554.tif 000033
【０１４１】
理論値算出部１４０には、式（３'）から式（６'）が予め与えられてもよい。この場合、理論値算出部１４０は、ヌル周波数検出部１３０が検出した第１ヌル周波数を、式（３'）から式（６'）に適用することで、確定成分の各種類のスペクトルの理論値を算出してよい。
【０１４２】
測定値算出部１５０は、標準偏差算出部１１０が算出した標準偏差σＴＪ、ＲＭＳ、および、スペクトル算出部１２０が算出したスペクトルに基づいて、確率密度関数ＰＤＦに含まれる確定成分のスペクトルの測定値を算出する。測定値算出部１５０は、スペクトル算出部１２０が算出したスペクトルＨ（ｆ）の周波数を、第１ヌル周波数ｆｚｅｒｏで正規化したスペクトルＨ（ｘ）を用いてよい。
【０１４３】
測定値算出部１５０は、当該スペクトルＨ（ｘ）を、標準偏差σＴＪ、ＲＭＳを有するガウス分布のスペクトルで除算することで、確定成分のスペクトルの測定値を算出してよい。このとき、確定成分のスペクトルＹ（ｘ）の測定値は、式（７）から、下式で与えられる。
【０１４４】
【数３２】
2010164554.tif 000034
【０１４５】
図２３は、理論値算出部１４０が算出した、各種類の確定成分のスペクトルＹ（ｘ）の理論値と、測定値算出部１５０が算出した、確定成分のスペクトルＹ（ｘ）の測定値との一例を示す図である。図２３に示すように、０＜ｘ＜１の範囲において、下記の関係が成り立つ。
【０１４６】
【数３３】
2010164554.tif 000035
【０１４７】
また、台形分布における上辺／下辺の台形比αが大きいほどＹＴｒａ（ｘ，α）も大きくなる。つまり、
【０１４８】
【数３４】
2010164554.tif 000036
となる。
【０１４９】
モデル判定部１６０は、理論値算出部１４０が確定成分の種類毎に算出した理論値のうち、測定値算出部１５０が算出した測定値に最も近い理論値に対応する確定成分の種類を、確率密度関数に含まれる確定成分の種類と判定する。なお、理論値算出部１４０および測定値算出部１５０は、予め定められた同一の周波数（図２３の例では、ｆ／ｆｚｅｒｏ＝ｘａ）におけるスペクトルの値を、理論値および測定値として算出してよい。
【０１５０】
つまり、理論値算出部１４０および測定値算出部１５０は、すべての帯域に渡ってスペクトルの理論値および測定値を算出しなくともよい。ただし、理論値算出部１４０および測定値算出部１５０は、０より大きく、且つ、第１ヌル周波数より小さい周波数の範囲において、スペクトルの値を算出することが好ましい。つまり、理論値算出部１４０および測定値算出部１５０は、０＜ｘａ＜１の範囲で、スペクトルの値を算出することが好ましい。
【０１５１】
なお、図２３の例では、測定値は、台形分布（α＝０．２）の理論値に最も近いので、モデル判定部１６０は、確率密度関数に含まれる確定成分の種類を、台形分布（α＝０．２）と判定してよい。また、モデル判定部１６０は、測定値算出部１５０が算出した測定値が、一様分布の理論値と、三角分布の理論値との間にあると判定した場合、理論値が測定値により近くなるような台形比αを算出してよい。モデル判定部１６０は、算出された台形比αを有する台形分布を、確定成分のモデルとして判定してよい。
【０１５２】
このような処理により、確率密度関数に含まれる確定成分のモデルを精度よく判定することができる。このため、確率密度関数に含まれるランダム成分および確定成分を精度よく算出することができる。
【０１５３】
図２４は、ヌル周波数検出部１３０の動作例を示す図である。本例におけるヌル周波数検出部１３０は、スペクトル算出部１２０が算出したスペクトルを、周波数で２階微分した波形のピークに基づいて、スペクトルの第１ヌル周波数を検出する。
【０１５４】
本例において、スペクトルの第１ヌル周波数をｆ１とする。与えられる確率密度関数にノイズが少ない場合、スペクトルの第１ヌル周波数は精確に検出することができる。これに対し、与えられる確率密度関数にノイズが含まれる場合、図２４のスペクトルｇ（ｆ）に示すように、検出されるべき周波数ｆ１において、第１ヌル周波数を検出することができないことがある。
【０１５５】
この場合、図２４に示すように、当該スペクトルを周波数で微分することにより、第１ヌル周波数を精度よく検出することができる。図２４に示すように、当該スペクトルの２階微分波形ｇ''（ｆ）のピークが、スペクトルｇ（ｆ）のヌルに対応する。このため、ヌル周波数検出部１３０は、確率密度関数のスペクトルを２階微分し、微分波形のピーク周波数に基づいて、第１ヌル周波数を検出してよい。
【０１５６】
図２５は、確定成分識別装置１００の他の構成例を示す図である。本例における確定成分識別装置１００は、図１７に関連して説明した確定成分識別装置１００の構成に加え、確定成分算出部１７０およびランダム成分算出部１８０を更に備える。他の構成要素は、図１７に関連して説明した構成要素と同一であってよい。
【０１５７】
確定成分算出部１７０は、モデル判定部１６０が判定した確定成分の種類と、ヌル周波数検出部１３０が検出した第１ヌル周波数とに基づいて、確率密度関数に含まれる確定成分を算出する。図２２に示したように、時間ドメインまたは周波数ドメインにおける確定成分の確率密度関数は、確定成分の種類と、ピークツゥピーク値ＤＪＰ−Ｐとから決定することができる。
【０１５８】
確定成分算出部１７０は、第１ヌル周波数からピークツゥピーク値ＤＪＰ−Ｐを算出することで、確定成分の確率密度関数を算出してよい。確定成分算出部１７０には、図２２に示したような、確定成分の種類毎の、時間ドメインでのモデル式、周波数ドメインでのモデル式、第１ヌル周波数ｆｚｅｒｏとピークツゥピーク値ＤＪＰ−Ｐとの関係、および、ピークツゥピーク値ＤＪＰ-Ｐと二乗平均値ＤＪＲＭＳとの関係を示すテーブルが与えられてよい。
【０１５９】
ランダム成分算出部１８０は、確定成分算出部１７０が算出した確定成分を、確定成分識別装置１００に与えられる確率密度関数ＰＤＦから除去することで、確率密度関数ＰＤＦに含まれるランダム成分を算出する。例えばランダム成分算出部１８０は、時間ドメインの確率密度関数ＰＤＦから、時間ドメインの確定成分をデコンボリューションすることで、時間ドメインのランダム成分を算出してよい。また、ランダム成分算出部１８０は、周波数ドメインの確率密度関数ＰＤＦから、周波数ドメインの確定成分を除算することで、周波数ドメインのランダム成分を算出してよい。また、ランダム成分算出部１８０は、算出したランダム成分の標準偏差を更に算出してもよい。
【０１６０】
このような構成により、確率密度関数に含まれる確定成分およびランダム成分を、精度よく分離することができる。このため、例えばジッタを測定したときに、確定ジッタおよびランダムジッタを精度よく分離することができるので、評価対象を精度よく評価することができる。
【０１６１】
図２６は、確定成分識別装置１００により、確定成分のモデルを判定する処理の概要を示すフローチャートである。まず、標準偏差算出部１１０が、確率密度関数の標準偏差を算出する（Ｓ２００）。また、スペクトル算出部１２０が、確率密度関数のスペクトルを算出する（Ｓ２０２）。また、ヌル周波数検出部１３０は、スペクトル算出部１２０が算出したスペクトルの第１ヌル周波数を検出する（Ｓ２０４）。Ｓ２０２およびＳ２０４の処理は、Ｓ２００の処理と並行して行ってよい。
【０１６２】
理論値算出部１４０は、ヌル周波数検出部１３０が検出した第１ヌル周波数に基づいて、確定成分の各種類の理論値を算出する（Ｓ２０６）。また、測定値算出部１５０は、標準偏差算出部１１０が算出した標準偏差およびスペクトル算出部１２０が算出したスペクトルに基づいて、確定成分の測定値を算出する（Ｓ２０８）。Ｓ２０８の処理は、Ｓ２０６の処理と並行して行ってよい。
【０１６３】
そして、モデル判定部１６０は、理論値算出部１４０が算出した各理論値と、測定値算出部１５０が算出した測定値とを比較して、確率密度関数に含まれる確定成分の種類を判定する（Ｓ２１０）。このような処理により、確定成分の種類を精度よく判定することができる。
【０１６４】
図２７は、確定成分識別装置１００の他の構成例を示す図である。本例の確定成分識別装置１００は、台形分布の確定成分の台形比αを算出する。本例の確定成分識別装置１００は、図１７に関連して説明した確定成分識別装置１００の構成に対して、モデル判定部１６０に代えて比算出部１９０を備える。標準偏差算出部１１０、スペクトル算出部１２０、ヌル周波数検出部１３０、理論値算出部１４０、および、測定値算出部１５０は、図１７から図２６において同一の符号を付して説明した構成要素と同一の機能および構成を有してよい。
【０１６５】
比算出部１９０は、標準偏差算出部１１０が算出した確率密度関数の標準偏差、スペクトル算出部１２０が算出した確率密度関数のスペクトル、および、ヌル周波数検出部１３０が検出したヌル周波数に基づいて、台形比αを算出する。式（４）または式（４'）に示したように、確率密度関数に含まれる確定成分が台形分布のとき、確率密度関数のスペクトルは、第１ヌル周波数ｆｚｅｒｏ、確率密度関数の標準偏差σＴＪ、ＲＭＳ、および、台形比αにより決定される。このため、確率密度関数のスペクトル、第１ヌル周波数ｆｚｅｒｏ、および、確率密度関数の標準偏差σＴＪ、ＲＭＳから、台形比αを算出することができる。
【０１６６】
測定値算出部１５０は、標準偏差算出部１１０が算出した標準偏差σＴＪ、ＲＭＳ、および、スペクトル算出部１２０が算出したスペクトルＨ（ｘ）に基づいて、確率密度関数に含まれる確定成分に対して、スペクトルの測定値を算出する。例えば測定値算出部１５０は、式（８）に基づいて、スペクトルの測定値を算出してよい。また、測定値算出部１５０は、予め定められた周波数ｘ１（ｘ１＝ｆ１／ｆｚｅｒｏ）における、スペクトルの測定値Ｙ（ｘ１）を算出してよい。
【０１６７】
比算出部１９０は、台形比αおよび第１ヌル周波数ｆｚｅｒｏにより定まる確定成分のスペクトルの理論値が、スペクトルの測定値で近似されるように、台形比αを算出する。例えば比算出部１９０は、予め定められた複数種類の台形比αに対する確定成分のスペクトルの理論値が与えられ、測定値算出部１５０が算出した測定値に最も近い理論値を与える台形比αを選択してよい。
【０１６８】
理論値算出部１４０は、上述した複数種類の台形比αに対する、確定成分のスペクトルの理論値をそれぞれ算出してよい。このとき、理論値算出部１４０は、予め定められた周波数ｘ１における理論値を算出してよい。理論値算出部１４０は、式（４'）に基づいて、それぞれの台形比αに対する理論値を算出してよい。
【０１６９】
理論値算出部１４０は、それぞれの台形比αに対する理論値を、比算出部１９０に通知する。上述したように、比算出部１９０は、これらの理論値と、測定値算出部１５０が算出した測定値とを比較して、測定値に最も近い理論値を与える台形比αを選択してよい。また、比算出部１９０は、理論値算出部１４０が用いたそれぞれの台形比αの間における理論値を、補間法により算出することで、測定値と略同一の理論値を与える台形比αを算出してもよい。このとき、比算出部１９０は、線形補間、スプライン補間等の公知の補間法を用いてよい。
【０１７０】
また、台形比αを決定する他の方法として、比算出部１９０は、理論値算出部１４０から与えられる複数の台形比αに対する理論値に基づいて、所定の周波数ｘ１における確定成分のスペクトルの測定値Ｈ（ｘ１）を変数として台形比αを求める式を生成してもよい。例えば、比算出部１９０は、Ｎ組の台形比αおよび理論値Ｈｉｄｅａｌ（ｘ１）の組み合わせを、理論値算出部１４０から受け取ってよい。
【０１７１】
そして、比算出部１９０は、理論値Ｈｉｄｅａ（ｘ１）のＮ−１次の多項式で、αを近似する式を生成してよい。ここで、理論値Ｈｉｄｅａｌ（ｘ１）＝ｍとすると、Ｎ−１次の近似式は、下式で与えられる。ただし、ｋｉは、ｉ次項の係数を示す。
【０１７２】
【数３５】
2010164554.tif 000037
【０１７３】
比算出部１９０は、式（９）のｘに、測定値算出部１５０が算出した確定成分のスペクトルの測定値Ｈ（ｘ）を代入することで、台形比αを算出してよい。このような方法によっても、台形比αを決定することができる。
【０１７４】
また、台形比αを決定する他の方法として、比算出部１９０は、式（４）に所定の台形比αを代入して得られる理論値と、測定値算出部１５０が算出した測定値との誤差が、所定の許容範囲内となる台形比αを、２分法等の探索法で検出してもよい。このような方法によっても、台形比αを決定することができる。
【０１７５】
図２８は、図２７に関連して説明した確定成分識別装置１００の動作例を説明するフローチャートである。本例において、Ｓ２００からＳ２０８までの処理は、図２６に関連して説明したＳ２００からＳ２０８までの処理と同一であってよい。Ｓ２０８の処理の後、比算出部１９０は、確定成分の分布における上辺および底辺の比である台形比αを算出する（Ｓ２１２）。以上のような処理により、確率密度関数に含まれる台形分布の確定成分の台形比を算出することができる。
【０１７６】
図２９は、確定成分識別装置１００の他の構成例を示す図である。本例における確定成分識別装置１００は、図１７に関連して説明した確定成分識別装置１００の構成に加え、比算出部１９０を更に備える。比算出部１９０は、図２７に関連して説明した比算出部１９０と同一であってよい。
【０１７７】
本例の確定成分識別装置１００は、まずモデル判定部１６０において、確率密度関数に含まれる確定成分が台形分布か否かを判定する。モデル判定部１６０は、確率密度関数に含まれる確定成分のモデルが、サイン波分布、一様分布、デュアルディラック分布、台形分布等のいずれであるかを判定してよい。そして、モデル判定部１６０において、確率密度関数に含まれる確定成分が台形分布と判定された場合に、比算出部１９０において確定成分の台形比αを算出する。
【０１７８】
なお、図２３に関連して説明した例では、測定値が、台形比α＝０．２の台形分布の理論値に近い場合に、確定分布のモデルが台形分布と判定している。しかし、台形分布における台形比αは、０から１の間の値を取りうる。このため、モデル判定部１６０は、測定値算出部１５０が算出した測定値が、三角分布（つまり、台形比α＝０）の理論値と、一様分布（つまり、台形比α＝１）の理論値との間となる場合に、確定成分のモデルを台形分布と判定してよい。
【０１７９】
また他の例では、モデル判定部１６０は、測定値が、所定の台形比αを有する台形分布（例えば、α＝０．５）の理論値に近い場合に、確定成分のモデルを台形分布と判定してもよい。以上のような処理により、確率密度関数に含まれる確定成分のモデルを判定することができ、更に、確定成分のモデルが台形分布のときに、その台形比αを算出することができる。このため、確定成分の分布形状を精度よく識別することができる。
【０１８０】
なお、本例の確定成分識別装置１００においても、図２５に関連して説明した確定成分算出部１７０およびランダム成分算出部１８０を更に備えてよい。この場合、確定成分算出部１７０は、比算出部１９０が算出した台形比αを更に受け取ってよい。確定成分算出部１７０は、確定成分の種類、第１ヌル周波数ｆｚｅｒｏ、および、台形比αから、図２２に示した式を用いて、確定成分を算出してよい。
【０１８１】
図３０は、図２９に示した確定成分識別装置１００の動作例を説明するフローチャートである。本例において、Ｓ２００からＳ２１０までの処理は、図２６において説明したＳ２００からＳ２１０までの処理と同一であってよい。Ｓ２１０の処理の後、モデル判定部１６０は、確定成分のモデルが、台形分布であるか否かを判断する（Ｓ２１１）。確定成分のモデルが台形分布のとき、モデル判定部１６０は、その旨を比算出部１９０に通知する。また、確定成分のモデルが台形分布でない場合、確定成分識別装置１００は、処理を終了してよい。
【０１８２】
比算出部１９０は、モデル判定部１６０から、当該通知を受け取った場合に、台形比αを算出する（Ｓ２１２）。このとき、理論値算出部１４０は、図２７に関連して説明した各理論値を算出して、比算出部１９０に通知してよい。比算出部１９０は、測定値算出部１５０が算出した測定値と、理論値算出部１４０が算出した理論値とを比較して、台形比αを算出する。このような処理により、確率密度関数に含まれる確定成分のモデルを判定することができ、更に、確定成分のモデルが台形分布のときに、その台形比αを算出することができる。
【０１８３】
図３１は、一つの実施形態に係る電子デバイス４００の構成例を示す図である。本例の電子デバイス４００は、入力ピン４０２から与えられる信号に応じて動作して、生成した所定の信号を出力ピン４０４から出力する。電子デバイス４００は、動作回路４１０、および、評価装置２０を備える。
【０１８４】
動作回路４１０は、与えられる入力信号に応じて動作する。動作回路４１０は、動作結果に応じて所定の出力信号を生成して出力してよい。評価装置２０は、動作回路４１０が正常に動作するか否かを試験するＢＩＳＴ回路として機能する。評価装置２０は、図１から図４に関連して説明した評価装置２０と同一の構成である。
【０１８５】
また、評価装置２０は、良否判定結果を、試験ピン４０６を介して外部に出力してよい。以上の構成により、動作回路４１０の自己診断を精度よく行える電子デバイス４００を提供できる。
【０１８６】
図３２は、一つの実施形態に係るコンピュータ１９００のハードウェア構成の一例を示す。コンピュータ１９００は、与えられるプログラムに基づいて、図１から図４において説明した評価装置２０として機能する。プログラムは、コンピュータ１９００を、図１から図４に関連して説明した評価装置２０の各構成要素として機能させてよい。
【０１８７】
本実施形態に係るコンピュータ１９００は、ＣＰＵ周辺部、入出力部、及びレガシー入出力部を備える。ＣＰＵ周辺部は、ホスト・コントローラ２０８２により相互に接続されるＣＰＵ２０００、ＲＡＭ２０２０、グラフィック・コントローラ２０７５、及び表示装置２０８０を有する。入出力部は、入出力コントローラ２０８４によりホスト・コントローラ２０８２に接続される通信インターフェース２０３０、ハードディスクドライブ２０４０、及びＣＤ−ＲＯＭドライブ２０６０を有する。レガシー入出力部は、入出力コントローラ２０８４に接続されるＲＯＭ２０１０、フレキシブルディスク・ドライブ２０５０、及び入出力チップ２０７０を有する。
【０１８８】
ホスト・コントローラ２０８２は、ＲＡＭ２０２０と、高い転送レートでＲＡＭ２０２０をアクセスするＣＰＵ２０００及びグラフィック・コントローラ２０７５とを接続する。ＣＰＵ２０００は、ＲＯＭ２０１０及びＲＡＭ２０２０に格納されたプログラムに基づいて動作し、各部の制御を行う。グラフィック・コントローラ２０７５は、ＣＰＵ２０００等がＲＡＭ２０２０内に設けたフレーム・バッファ上に生成する画像データを取得し、表示装置２０８０上に表示させる。これに代えて、グラフィック・コントローラ２０７５は、ＣＰＵ２０００等が生成する画像データを格納するフレーム・バッファを、内部に含んでもよい。
【０１８９】
入出力コントローラ２０８４は、ホスト・コントローラ２０８２と、比較的高速な入出力装置である通信インターフェース２０３０、ハードディスクドライブ２０４０、ＣＤ−ＲＯＭドライブ２０６０を接続する。通信インターフェース２０３０は、ネットワークを介して他の装置と通信する。ハードディスクドライブ２０４０は、コンピュータ１９００内のＣＰＵ２０００が使用するプログラム及びデータを格納する。ＣＤ−ＲＯＭドライブ２０６０は、ＣＤ−ＲＯＭ２０９５からプログラム又はデータを読み取り、ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供する。
【０１９０】
また、入出力コントローラ２０８４には、ＲＯＭ２０１０と、フレキシブルディスク・ドライブ２０５０、及び入出力チップ２０７０の比較的低速な入出力装置とが接続される。ＲＯＭ２０１０は、コンピュータ１９００が起動時に実行するブート・プログラム、コンピュータ１９００のハードウェアに依存するプログラム等を格納する。フレキシブルディスク・ドライブ２０５０は、フレキシブルディスク２０９０からプログラム又はデータを読み取り、ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供する。入出力チップ２０７０は、フレキシブルディスク・ドライブ２０５０、例えばパラレル・ポート、シリアル・ポート、キーボード・ポート、マウス・ポート等を介して各種の入出力装置を接続する。
【０１９１】
ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供されるプログラムは、フレキシブルディスク２０９０、ＣＤ−ＲＯＭ２０９５、又はＩＣカード等の記録媒体に格納されて利用者によって提供される。プログラムは、記録媒体から読み出され、ＲＡＭ２０２０を介してコンピュータ１９００内のハードディスクドライブ２０４０にインストールされ、ＣＰＵ２０００において実行される。当該プログラムは、コンピュータ１９００にインストールされる。当該プログラムは、ＣＰＵ２０００等に働きかけて、コンピュータ１９００を、確定成分識別装置１００として機能させる。
【０１９２】
以上に示したプログラムは、外部の記録媒体に格納されてもよい。記録媒体としては、フレキシブルディスク２０９０、ＣＤ−ＲＯＭ２０９５の他に、ＤＶＤ、ＣＤ等の光学記録媒体、ＭＯ等の光磁気記録媒体、テープ媒体、ＩＣカード等の半導体メモリ等を用いることができる。また、専用通信ネットワークまたはインターネットに接続されたサーバシステムに設けたハードディスク又はＲＡＭ等の記憶装置を記録媒体として使用し、ネットワークを介してプログラムをコンピュータ１９００に提供してもよい。
【０１９３】
以上、本発明の（一）側面を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施の形態に記載の範囲には限定されない。上記実施の形態に、多様な変更又は改良を加えることができる。その様な変更又は改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲の記載から明らかである。
【０１９４】
特許請求の範囲、明細書、および図面中において示した装置、システム、プログラム、および方法における動作、手順、ステップ、および段階等の各処理の実行順序は、特段「より前に」、「先立って」等と明示しておらず、また、前の処理の出力を後の処理で用いるのでない限り、任意の順序で実現しうることに留意すべきである。特許請求の範囲、明細書、および図面中の動作フローに関して、便宜上「まず、」、「次に、」等を用いて説明したとしても、この順で実施することが必須であることを意味するものではない。
【０１９５】
上記説明から明らかなように、本発明の（一）実施形態によれば、低いコストおよび短い測定時間で伝播装置の特性を評価することができる評価装置を実現することができる。
【符号の説明】
【０１９６】
１０  算出装置
１２  モデル判別部
１４  パラメータ算出部
１６  指定部
１８  計算部
３４  良否判定部
２０  評価装置
２２  信号発生部
２４  入力信号測定部
２６  出力信号測定部
２８  分離部
３０  評価部
３６  調整部
４０  伝播装置
１００  確定成分識別装置
１１０  標準偏差算出部
１２０  スペクトル算出部
１３０  ヌル周波数検出部
１４０  理論値算出部
１５０  測定値算出部
１６０  モデル判定部
１７０  確定成分算出部
１８０  ランダム成分算出部
１９０  比算出部
４００  電子デバイス
４０２  入力ピン
４０４  出力ピン
４０６  試験ピン
４１０  動作回路
１９００  コンピュータ
２０００  ＣＰＵ
２０１０  ＲＯＭ
２０２０  ＲＡＭ
２０３０  通信インターフェース
２０４０  ハードディスクドライブ
２０５０  フレキシブルディスク・ドライブ
２０６０  ＣＤ−ＲＯＭドライブ
２０７０  入出力チップ
２０７５  グラフィック・コントローラ
２０８０  表示装置
２０８２  ホスト・コントローラ
２０８４  入出力コントローラ
２０９０  フレキシブルディスク
２０９５  ＣＤ−ＲＯＭ
【図１】
2010164554.tif 000038
【図２】
2010164554.tif 000039
【図３】
2010164554.tif 000040
【図４】
2010164554.tif 000041
【図５】
2010164554.tif 000042
【図６】
2010164554.tif 000043
【図７】
2010164554.tif 000044
【図８】
2010164554.tif 000045
【図９】
2010164554.tif 000046
【図１０】
2010164554.tif 000047
【図１１】
2010164554.tif 000048
【図１２】
2010164554.tif 000049
【図１３】
2010164554.tif 000050
【図１４】
2010164554.tif 000051
【図１５】
2010164554.tif 000052
【図１６】
2010164554.tif 000053
【図１７】
2010164554.tif 000054
【図１８Ａ】
2010164554.tif 000055
【図１８Ｂ】
2010164554.tif 000056
【図１９Ａ】
2010164554.tif 000057
【図１９Ｂ】
2010164554.tif 000058
【図２０Ａ】
2010164554.tif 000059
【図２０Ｂ】
2010164554.tif 000060
【図２１】
2010164554.tif 000061
【図２２】
2010164554.tif 000062
【図２３】
2010164554.tif 000063
【図２４】
2010164554.tif 000064
【図２５】
2010164554.tif 000065
【図２６】
2010164554.tif 000066
【図２７】
2010164554.tif 000067
【図２８】
2010164554.tif 000068
【図２９】
2010164554.tif 000069
【図３０】
2010164554.tif 000070
【図３１】
2010164554.tif 000071
【図３２】
2010164554.tif 000072