(57)【特許請求の範囲】
【請求項１】
離散化された被測定信号のジッタを測定するジッタ測定装置であって、
前記被測定信号の立ち上がりエッジまたは立ち下がりエッジのいずれかのエッジの標本点から、所定の値に最も近い振幅値を有する標本点を、前記被測定信号のサイクル毎に抽出し、前記所定の値となる前記被測定信号の点と、前記標本点との位相誤差を前記被測定信号のサイクル毎に算出して位相誤差系列を生成し、抽出した前記標本点において、隣り合う前記標本点の間の時間間隔系列を生成する位相誤差推定器と、
前記位相誤差系列及び標本点間の前記時間間隔系列から、前記被測定信号の１サイクル毎の瞬時周期を算出し、算出した１サイクル毎の前記瞬時周期から前記被測定信号の基本周期を減算することで、前記被測定信号の周期ジッタ系列を算出する周期ジッタ推定器と
を備えることを特徴とするジッタ測定装置。
【請求項２】
前記位相誤差推定器は、前記被測定信号における複数のゼロクロス点と、前記ゼロクロス点に最も近い振幅値を有する近似ゼロクロス点との位相誤差をそれぞれ算出し、前記位相誤差系列を生成することを特徴とする請求項１に記載のジッタ測定装置。
【請求項３】
前記被測定信号の基本周波数を測定する基本周期測定器を更に備える請求項１に記載のジッタ測定装置。
【請求項４】
前記位相誤差推定器は、
前記近似ゼロクロス点を抽出し、隣り合う前記近似ゼロクロス点間の前記時間間隔系列を算出するゼロクロスサンプラと、
前記近似ゼロクロス点と、対応する前記ゼロクロス点との前記位相誤差系列を算出する位相誤差算出器と
を有することを特徴とする請求項２に記載のジッタ測定装置。
【請求項５】
前記周期ジッタ推定器は、
前記時間間隔系列及び前記位相誤差系列に基づいて、前記被測定信号の１サイクル毎の瞬時周期を算出する瞬時周期推定器と、
前記瞬時周期推定器が算出したそれぞれの瞬時周期値から、前記被測定信号の基本周期を減算して、前記周期ジッタ系列を生成する基本周期減算器と
を有することを特徴とする請求項４に記載のジッタ測定装置。
【請求項６】
前記被測定信号の基本周波数近傍の周波数成分を抽出して、前記位相誤差推定器に供給する帯域通過フィルタを更に備えることを特徴とする請求項５に記載のジッタ測定装置。
【請求項７】
前記位相誤差算出器は、
前記被測定信号の振幅値と、
それぞれの前記近似ゼロクロス点における振幅値と
に基づいて、それぞれの前記近似ゼロクロス点と、対応する前記ゼロクロス点との位相誤差をそれぞれ算出し、前記位相誤差系列を算出することを特徴とする請求項６に記載のジッタ測定装置。
【請求項８】
前記位相誤差推定器は、前記被測定信号を複素数の解析信号に変換し、前記ゼロクロスサンプラに供給する解析信号変換器を更に有し、
前記ゼロクロスサンプラは、前記解析信号の実数部において複数の前記近似ゼロクロス点を抽出し、
前記位相誤差算出器は、それぞれの前記近似ゼロクロス点における前記解析信号の値に基づいて、それぞれの前記近似ゼロクロス点の瞬時位相を算出し、算出したそれぞれの前記瞬時位相に基づいて、前記ゼロクロス点と前記近似ゼロクロス点との位相誤差系列を算出することを特徴とする請求項６に記載のジッタ測定装置。
【請求項９】
前記ゼロクロスサンプラは、前記被測定信号の前記ゼロクロス点近傍をサンプリングし、前記近似ゼロクロス点を抽出することを特徴とする請求項４に記載のジッタ測定装置。
【請求項１０】
前記周期ジッタ推定器は、基本周期減算器が生成した前記周期ジッタ系列におけるそれぞれの周期ジッタ値に、対応する前記近似ゼロクロス点間の時間間隔と前記被測定信号の基本周期との比を乗算し、前記周期ジッタ系列を補正する周期ジッタ系列補正器を更に有することを特徴とする請求項５に記載のジッタ測定装置。
【請求項１１】
前記周期ジッタ推定器は、
前記時間間隔系列と前記位相誤差系列とに基づいて、それぞれの前記近似ゼロクロス点間における前記被測定信号の瞬時角周波数を示す瞬時角周波数系列を算出する瞬時角周波数推定器と、
前記瞬時角周波数系列と、前記被測定信号の基本周期とに基づいて、前記周期ジッタ系列を算出する周期ジッタ算出器と
を有することを特徴とする請求項４に記載のジッタ測定装置。
【請求項１２】
前記周期ジッタ系列に基づいて、前記被測定信号の周期ジッタのピーク・ツゥ・ピーク値、実効値、ヒストグラムのうちのいずれか一つ又は複数を算出するジッタ検出器を更に備えることを特徴とする請求項１に記載のジッタ測定装置。
【請求項１３】
前記周期ジッタ系列において、隣り合う周期ジッタ値の差分をそれぞれ算出したサイクル・ツゥ・サイクル周期ジッタ系列を生成するサイクル・ツゥ・サイクル周期ジッタ推定器を更に備えることを特徴とする請求項１に記載のジッタ測定装置。
【請求項１４】
前記帯域通過フィルタは、
前記被測定信号を周波数領域の信号に変換する周波数領域変換器と、
前記周波数領域の信号に変換された前記被測定信号の、基本周波数近傍の周波数成分を抽出する帯域制限器と、
前記帯域制限器が抽出した前記周波数成分を、時間領域の信号に変換する時間領域変換器と
を有することを特徴とする請求項６に記載のジッタ測定装置。
【請求項１５】
前記被測定信号の振幅変調成分を除去し、前記位相誤差推定器に供給する波形クリッパを更に備えることを特徴とする請求項１に記載のジッタ測定装置。
【請求項１６】
前記被測定信号の前記基本周期を算出する基本周期推定器を更に備えることを特徴とする請求項５又は１１に記載のジッタ測定装置。
【請求項１７】
前記周期ジッタ推定器は、前記周期ジッタ系列における高周波成分を除去する低域通過フィルタを更に有することを特徴とする請求項５又は１１に記載のジッタ測定装置。
【請求項１８】
離散化された被測定信号のジッタを測定するジッタ測定方法であって、
前記被測定信号の立ち上がりエッジまたは立ち下がりエッジのいずれかのエッジの標本点から、所定の値に最も近い振幅値を有する標本点を、前記被測定信号のサイクル毎に抽出し、前記所定の値となる前記被測定信号の点と、前記標本点との位相誤差を前記被測定信号のサイクル毎に算出して位相誤差系列を生成し、抽出した前記標本点において、隣り合う前記標本点の間の時間間隔系列を生成する位相誤差推定ステップと、
前記位相誤差系列及び標本点間の前記時間間隔系列から、前記被測定信号の１サイクル毎の瞬時周期を算出し、算出した１サイクル毎の前記瞬時周期から前記被測定信号の基本周期を減算することで、前記被測定信号の周期ジッタ系列を算出する周期ジッタ推定ステップと
を備えることを特徴とするジッタ測定方法。
【請求項１９】
前記位相誤差推定ステップは、前記被測定信号における複数のゼロクロス点と、前記ゼロクロス点に最も近い振幅値を有する近似ゼロクロス点との位相誤差をそれぞれ算出し、前記位相誤差系列を生成することを特徴とする請求項１８に記載のジッタ測定方法。
【請求項２０】
前記被測定信号の基本周波数を測定する基本周期測定ステップを更に備える請求項１８に記載のジッタ測定方法。
【請求項２１】
前記位相誤差推定ステップは、
複数の前記近似ゼロクロス点を抽出し、隣り合う前記近似ゼロクロス点間の前記時間間隔系列を算出するゼロクロスサンプリングステップと、
前記近似ゼロクロス点と、対応する前記ゼロクロス点との前記位相誤差系列を算出する位相誤差算出ステップと
を有することを特徴とする請求項１９に記載のジッタ測定方法。
【請求項２２】
前記周期ジッタ推定ステップは、
前記時間間隔系列及び前記位相誤差系列に基づいて、前記瞬時周期の系列である瞬時周期系列を算出する瞬時周期推定ステップと、
前記瞬時周期系列のそれぞれの周期値から、前記被測定信号の基本周期を減算して、前記周期ジッタ系列を生成する基本周期減算ステップと
を有することを特徴とする請求項２１に記載のジッタ測定方法。
【請求項２３】
前記被測定信号の基本周波数近傍の周波数成分を抽出する帯域通過ステップを更に備えることを特徴とする請求項２２に記載のジッタ測定方法。
【請求項２４】
前記位相誤差算出ステップは、
前記被測定信号の振幅値と、
それぞれの前記近似ゼロクロス点における振幅値と
に基づいて、それぞれの前記近似ゼロクロス点と、対応する前記ゼロクロス点との位相誤差をそれぞれ算出し、前記位相誤差系列を算出することを特徴とする請求項２３に記載のジッタ測定方法。
【請求項２５】
前記位相誤差推定ステップは、前記被測定信号を複素数の解析信号に変換する解析信号変換ステップを更に有し、
前記ゼロクロスサンプリングステップは、前記解析信号の実数部において複数の前記近似ゼロクロス点を抽出し、
前記位相誤差算出ステップは、それぞれの前記近似ゼロクロス点における前記解析信号の値に基づいて、それぞれの前記近似ゼロクロス点の瞬時位相を算出し、算出したそれぞれの前記瞬時位相に基づいて、前記ゼロクロス点と前記近似ゼロクロス点との位相誤差系列を算出することを特徴とする請求項２３に記載のジッタ測定方法。
【請求項２６】
前記ゼロクロスサンプリングステップは、前記被測定信号の前記ゼロクロス点近傍をサンプリングし、前記近似ゼロクロス点を抽出することを特徴とする請求項２１に記載のジッタ測定方法。
【請求項２７】
前記周期ジッタ推定ステップは、基本周期減算ステップにおいて生成した前記周期ジッタ系列におけるそれぞれの周期ジッタ値に、対応する前記近似ゼロクロス点間の時間間隔と前記被測定信号の基本周期との比を乗算し、前記周期ジッタ系列を補正する周期ジッタ系列補正ステップを更に有することを特徴とする請求項２２に記載のジッタ測定方法。
【請求項２８】
前記周期ジッタ推定ステップは、
前記時間間隔系列と前記位相誤差系列とに基づいて、それぞれの前記近似ゼロクロス点間における前記被測定信号の瞬時角周波数を示す瞬時角周波数系列を算出する瞬時角周波数推定ステップと、
前記瞬時角周波数系列と、前記被測定信号の基本周期とに基づいて、前記周期ジッタ系列を算出する周期ジッタ算出ステップと
を有することを特徴とする請求項２１に記載のジッタ測定方法。
【請求項２９】
前記周期ジッタ系列に基づいて、前記被測定信号の周期ジッタのピーク・ツゥ・ピーク値、実効値、ヒストグラムのうちのいずれか一つ又は複数を算出するジッタ検出ステップを更に備えることを特徴とする請求項１８に記載のジッタ測定方法。
【請求項３０】
前記周期ジッタ系列において、隣り合う周期ジッタ値の差分をそれぞれ算出したサイクル・ツゥ・サイクル周期ジッタ系列を生成するサイクル・ツゥ・サイクル周期ジッタ推定ステップを更に備えることを特徴とする請求項１８に記載のジッタ測定方法。
【請求項３１】
前記帯域通過ステップは、
前記被測定信号を周波数領域の信号に変換する周波数領域変換ステップと、
前記周波数領域の信号に変換された前記被測定信号の、基本周波数近傍の周波数成分を抽出する帯域制限ステップと、
前記帯域制限ステップにおいて抽出した前記周波数成分を、時間領域の信号に変換する時間領域変換ステップと
を有することを特徴とする請求項２３に記載のジッタ測定方法。
【請求項３２】
前記被測定信号の振幅変調成分を除去する波形クリップステップを更に備えることを特徴とする請求項１８に記載のジッタ測定方法。
【請求項３３】
前記被測定信号の前記基本周期を算出する基本周期推定ステップを更に備えることを特徴とする請求項２２又は２８に記載のジッタ測定方法。
【請求項３４】
前記周期ジッタ推定ステップは、前記周期ジッタ系列における高周波成分を除去する低域通過ステップを更に有することを特徴とする請求項２２又は２８に記載のジッタ測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電子回路を測定するジッタ測定装置及びジッタ測定方法に関する。特に、本発明は、電子回路のジッタを測定するジッタ測定装置に関する。
【背景技術】
【０００２】
従来、周期ジッタの測定には、タイムインターバルアナライザやオシロスコープが用いられている。これらは、いわゆるゼロクロス方式を用いて周期ジッタの測定を行う。まず、図１に示すように、被試験ＰＬＬ１１からの被測定信号（クロック信号）ｘ（ｔ）がタイムインターバルアナライザ１２に供給される。被測定信号ｘ（ｔ）は、１つの立ち上がりエッジに対し、次の立ち上がりエッジが点線のように揺らぎ、両エッジの間隔Ｔｐ、すなわち被測定信号ｘ（ｔ）の周期が揺らぐ。
【０００３】
ゼロクロス方式は、被測定信号のゼロクロス間の時間間隔（周期）を測定し、周期ジッタをヒストグラム解析により測定する。図２に、タイムインターバルアナライザ１２が表示するヒストグラムの一例を示す。タイムインターバルアナライザについては、例えばD.Chu,"Phase Digitizing Sharpens Timing Measurements"IEEE Spectrum,pp.28-32,1988.J.Wilstrup,"A Method of Serial Data Jitter Analysis Using One-Shot Time Interval Measurements" , Proceeding of IEEE International Test Conference ,pp.819-823,1988.に記載されている。
【０００４】
また、補間法を用いてジッタ測定を行うデジタルオシロスコープがある。この補間法を用いたジッタ測定方法（補間ベース・ジッタ測定方法）は、サンプリングされた被測定信号の測定データのうち、信号値がゼロクロスに近いデータ間を補完し、ゼロクロスのタイミングを推定する。即ち、データ補間によってゼロクロス間の時間間隔（周期）の測定誤差を低減し、被測定信号の周期の揺らぎを推定する。
【０００５】
つまり、図３に示すように、被試験ＰＬＬ１１からの被測定信号ｘ（ｔ）は、デジタルオシロスコープ１４に供給される。図４は、デジタルオシロスコープ１４の構成を示す。アナログデジタル変換器１５は、被測定信号をデジタルデータ列に変換し、補間器１６は、デジタルデータ列中の信号値がゼロクロスに近いデータ間をデータ補完し、周期推定器１７は、データ補完されたデジタルデータ列について、ゼロクロス間の周期を測定し、ヒストグラム推定器１８は、周期推定器１７が測定した測定値のヒストグラムを表示し、実効値及びピーク・ツゥ・ピーク値算出器１９は、測定された周期の揺らぎの自乗平均値及びピークツゥピーク値を算出する。例えば、被測定信号ｘ（ｔ）が図５（Ａ）に示す波形である場合、デジタルオシロスコープ１４は、図５（Ｂ）に示すような周期ジッタを測定する。
【０００６】
しかしながら、タイムインターバルアナライザ方式のジッタ測定法は、精度よくジッタを測定することができるが、１回の周期測定の後、測定を行うことができないデッドタイムがあるため、ヒストグラム解析に必要なデータ数を測定するのに時間がかかるという問題がある。
【０００７】
また、広帯域のデジタルオシロスコープと補間法を組み合わせたジッタ測定法は、ジッタ値を過大評価するという問題がある。つまり、ジッタ測定値がタイムインターバルアナライザによる測定値と互換性を持たない。図６（Ａ）及び図６（Ｂ）は、４００ＭＨｚの被測定信号に対するタイムインターバルアナライザによるジッタ測定結果と、補間法によるジッタ測定結果の比較結果を示す。これにより、タイムインターバルアナライザによるジッタ測定値は７．７２ｐｓ（ＲＭＳ）であるのに対して、補間法による測定値は８．４７ｐｓ（ＲＭＳ）であり、ジッタ測定値を過大評価していることがわかる。また、補間法によるジッタ測定法も、長い測定時間が必要となる。
【０００８】
これに対し、本願の発明者らは、T.J.Tamaguchi ,M.Soma,M.Ishida,
T.Watanabe, and T.Ohmi,"Extraction of Peak-to-Peak and RMS Sinusoidal Jitter Using an Analytic Signal Method." Procssdings of 18th IEEE VLSI Test Symposium,pp.395-402,2000.で以下に述べるジッタ測定方法を提案した。
【０００９】
即ち、図７に示すように、被試験ＰＬＬ回路１１からのアナログのクロック波形がＡ／Ｄ変換器２２でデジタルのクロック信号に変換され、解析信号変換手段としてのヒルベルト変換対発生器２３へ供給されて複素数の解析信号に変換される。
ヒルベルト変換対発生器２３において、クロック信号の基本周波数近傍の信号成分が、帯域通過フィルタにより帯域制限され、ヒルベルト変換器２５に供給されヒルベルト変換される。ヒルベルト変換対発生器２３は、解析信号の実数部としてヒルベルト変換されない信号成分を出力し、解析信号の虚数部としてヒルベルト変換された信号成分を出力する。
【００１０】
この解析信号から瞬時位相推定器２６でクロック信号の瞬時位相が推定される。この瞬時位相はリニア位相除去器２７でリニア位相が除去されて、瞬時位相雑音波形が得られる。
【００１１】
この瞬時位相雑音がゼロクロスサンプラ４３でサンプリングされ、タイミングジッタ系列としてピーク・ツゥ・ピーク検出器３２に入力され、タイミングジッタのピーク値が求められる。
【００１２】
また、タイミングジッタ系列が二乗平均検出器３３にも入力され、タイミングジッタの二乗平均値が求められる。
【００１３】
このようにして、タイミングジッタのピーク値、タイミングジッタの二乗平均値を、瞬時位相雑音から求める方法をΔφ法と呼ぶ。
【００１４】
このΔφ法によれば、タイミングジッタを高速に測定することができる。しかしながら、被測定信号のゼロクロス点を、各ゼロクロス点に最も近いサンプリング点で近似するため、信号離散化における過剰標本化比が小さい場合に、ジッタ測定の誤差が大きくなる場合があった。
【００１５】
これらのため、精度よく且つ短時間でジッタを測定できるジッタ測定装置及び方法が望まれていた。
【００１６】
本発明は、上記課題を解決することのできる測定装置及び測定方法を提供することを目的とする。この目的は、請求の範囲における独立項に記載の特徴の組み合わせにより達成される。また従属項は本発明の更なる有利な具体例を規定する。
【発明の開示】
【００１７】
本発明の第１の形態においては、離散化された被測定信号のジッタを測定するジッタ測定装置であって、被測定信号の立ち上がりエッジまたは立ち下がりエッジのいずれかのエッジの標本点から、所定の値に最も近い振幅値を有する標本点を、被測定信号のサイクル毎に抽出し、所定の値となる被測定信号の点と、標本点との位相誤差を被測定信号のサイクル毎に算出して位相誤差系列を生成し、抽出した標本点において、隣り合う標本点の間の時間間隔系列を生成する位相誤差推定器と、位相誤差系列及び標本点間の時間間隔系列から、被測定信号の１サイクル毎の瞬時周期を算出し、算出した１サイクル毎の瞬時周期から被測定信号の基本周期を減算することで、被測定信号の周期ジッタ系列を算出する周期ジッタ推定器とを備えることを特徴とするジッタ測定装置を提供する。
【００１８】
位相誤差推定器は、被測定信号における複数のゼロクロス点と、標本点のうち、複数のゼロクロス点のいずれかに対して、予め定められた条件を満たす複数の近似ゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を生成してよい。
【００１９】
位相誤差推定器は、被測定信号の各サイクルの立ち上がり、又は各サイクルの立ち下がりにおけるゼロクロス点と、標本点のうち、各サイクルの立ち上がり又は各サイクルの立ち下がりにおいて、ゼロクロス点に最も近い振幅値を有する標本点である複数の近似ゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を生成してよい。
【００２０】
位相誤差推定器は、複数の近似ゼロクロス点を抽出するゼロクロスサンプラと、近似ゼロクロス点と、対応するゼロクロス点との位相誤差系列を算出する位相誤差算出器とを有してよい。
【００２１】
ゼロクロスサンプラは、隣り合う近似ゼロクロス点間の時間間隔系列を算出し、周期ジッタ推定器は、１サイクル毎の瞬時周期を算出する瞬時周期推定器と、瞬時周期推定器が算出したそれぞれの瞬時周期値から、被測定信号の基本周期を減算して、周期ジッタ系列を生成する基本周期減算器とを有してよい。
【００２２】
ジッタ測定装置は、被測定信号の基本周波数近傍の周波数成分を抽出して、位相誤差推定器に供給する帯域通過フィルタを更に備えてよい。
【００２３】
位相誤差算出器は、被測定信号の振幅値と、それぞれの近似ゼロクロス点における振幅値とに基づいて、それぞれの近似ゼロクロス点と、対応するゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を算出してよい。
【００２４】
位相誤差推定器は、被測定信号を複素数の解析信号に変換し、ゼロクロスサンプラに供給する解析信号変換器を更に有し、ゼロクロスサンプラは、解析信号の実数部において複数の近似ゼロクロス点を抽出し、位相誤差算出器は、それぞれの近似ゼロクロス点における解析信号の値に基づいて、それぞれの近似ゼロクロス点の瞬時位相を算出し、算出したそれぞれの瞬時位相に基づいて、ゼロクロス点と近似ゼロクロス点との位相誤差系列を算出してよい。
【００２５】
ゼロクロスサンプラは、被測定信号のゼロクロス点近傍をサンプリングし、近似ゼロクロス点を抽出してよい。
【００２６】
周期ジッタ推定器は、基本周期減算器が生成した周期ジッタ系列におけるそれぞれの周期ジッタ値に、対応する近似ゼロクロス点間の時間間隔と被測定信号の基本周期との比を乗算し、周期ジッタ系列を補正する周期ジッタ系列補正器を更に有してよい。
【００２７】
ゼロクロスサンプラは、隣り合う近似ゼロクロス点間の時間間隔系列を算出し、周期ジッタ推定器は、時間間隔系列と位相誤差系列とに基づいて、それぞれの近似ゼロクロス点間における被測定信号の瞬時角周波数を示す瞬時角周波数系列を算出する瞬時角周波数推定器と、瞬時角周波数系列と、被測定信号の基本周期とに基づいて、周期ジッタ系列を算出する周期ジッタ算出器とを有してよい。
【００２８】
周期ジッタ系列に基づいて、被測定信号の周期ジッタのピーク・ツゥ・ピーク値、実効値、ヒストグラムのうちのいずれか一つ又は複数を算出するジッタ検出器を更に備えてよい。
【００２９】
周期ジッタ系列において、隣り合う周期ジッタ値の差分をそれぞれ算出したサイクル・ツゥ・サイクル周期ジッタ系列を生成するサイクル・ツゥ・サイクル周期ジッタ推定器を更に備えてよい。
【００３０】
帯域通過フィルタは、被測定信号を周波数領域の信号に変換する周波数領域変換器と、周波数領域の信号に変換された被測定信号の、基本周波数近傍の周波数成分を抽出する帯域制限器と、帯域制限器が抽出した周波数成分を、時間領域の信号に変換する時間領域変換器とを有してよい。
【００３１】
ジッタ測定装置は、被測定信号の振幅変調成分を除去し、位相誤差推定器に供給する波形クリッパを更に備えてよい。
【００３２】
ジッタ測定装置は、被測定信号の基本周期を算出する基本周期推定器を更に備えてよい。
【００３３】
周期ジッタ推定器は、周期ジッタ系列における高周波成分を除去する低域通過フィルタを更に有してよい。
【００３４】
本発明の第２の形態においては、離散化された被測定信号のジッタを測定するジッタ測定方法であって、被測定信号の立ち上がりエッジまたは立ち下がりエッジのいずれかのエッジの標本点から、所定の値に最も近い振幅値を有する標本点を、被測定信号のサイクル毎に抽出し、所定の値となる被測定信号の点と、標本点との位相誤差を被測定信号のサイクル毎に算出して位相誤差系列を生成し、抽出した標本点において、隣り合う標本点の間の時間間隔系列を生成する位相誤差推定ステップと、位相誤差系列及び標本点間の時間間隔系列から、被測定信号の１サイクル毎の瞬時周期を算出し、算出した１サイクル毎の瞬時周期から被測定信号の基本周期を減算することで、被測定信号の周期ジッタ系列を算出する周期ジッタ推定ステップとを備えることを特徴とするジッタ測定方法を提供する。
【００３５】
位相誤差推定ステップは、被測定信号における複数のゼロクロス点と、標本点のうち、複数のゼロクロス点のいずれかに対して、予め定められた条件を満たす複数の近似ゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を生成してよい。
【００３６】
位相誤差推定ステップは、被測定信号の各サイクルの立ち上がり、又は各サイクルの立ち下がりにおけるゼロクロス点と、標本点のうち、各サイクルの立ち上がり又は各サイクルの立ち下がりにおいて、ゼロクロス点に最も近い振幅値を有する標本点である複数の近似ゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を生成してよい。
【００３７】
位相誤差推定ステップは、複数の近似ゼロクロス点を抽出するゼロクロスサンプリングステップと、近似ゼロクロス点と、対応するゼロクロス点との位相誤差系列を算出する位相誤差算出ステップとを有してよい。
【００３８】
ゼロクロスサンプリングステップは、隣り合う近似ゼロクロス点間の時間間隔系列を算出し、周期ジッタ推定ステップは、時間間隔系列及び位相誤差系列に基づいて、１サイクル毎の瞬時周期を算出する瞬時周期推定ステップと、瞬時周期推定器が算出したそれぞれの瞬時周期値から、被測定信号の基本周期を減算して、周期ジッタ系列を生成する基本周期減算ステップとを有してよい。
【００３９】
被測定信号の基本周波数近傍の周波数成分を抽出する帯域通過ステップを更に備えてよい。
【００４０】
位相誤差算出ステップは、被測定信号の振幅値と、それぞれの近似ゼロクロス点における振幅値とに基づいて、それぞれの近似ゼロクロス点と、対応するゼロクロス点との位相誤差をそれぞれ算出し、位相誤差系列を算出してよい。
【００４１】
位相誤差推定ステップは、被測定信号を複素数の解析信号に変換する解析信号変換ステップを更に有し、ゼロクロスサンプリングステップは、解析信号の実数部において複数の近似ゼロクロス点を抽出し、位相誤差算出ステップは、それぞれの近似ゼロクロス点における解析信号の値に基づいて、それぞれの近似ゼロクロス点の瞬時位相を算出し、算出したそれぞれの瞬時位相に基づいて、ゼロクロス点と近似ゼロクロス点との位相誤差系列を算出してよい。
【００４２】
ゼロクロスサンプリングステップは、被測定信号のゼロクロス点近傍をサンプリングし、近似ゼロクロス点を抽出してよい。
【００４３】
周期ジッタ推定ステップは、基本周期減算ステップにおいて生成した周期ジッタ系列におけるそれぞれの周期ジッタ値に、対応する近似ゼロクロス点間の時間間隔と被測定信号の基本周期との比を乗算し、周期ジッタ系列を補正する周期ジッタ系列補正ステップを更に有してよい。
【００４４】
ゼロクロスサンプリングステップは、隣り合う近似ゼロクロス点間の時間間隔系列を算出し、周期ジッタ推定ステップは、時間間隔系列と位相誤差系列とに基づいて、それぞれの近似ゼロクロス点間における被測定信号の瞬時角周波数を示す瞬時角周波数系列を算出する瞬時角周波数推定ステップと、瞬時角周波数系列と、被測定信号の基本周期とに基づいて、周期ジッタ系列を算出する周期ジッタ算出ステップとを有してよい。
【００４５】
周期ジッタ系列に基づいて、被測定信号の周期ジッタのピーク・ツゥ・ピーク値、実効値、ヒストグラムのうちのいずれか一つ又は複数を算出するジッタ検出ステップを更に備えてよい。
【００４６】
周期ジッタ系列において、隣り合う周期ジッタ値の差分をそれぞれ算出したサイクル・ツゥ・サイクル周期ジッタ系列を生成するサイクル・ツゥ・サイクル周期ジッタ推定ステップを更に備えてよい。
【００４７】
帯域通過ステップは、被測定信号を周波数領域の信号に変換する周波数領域変換ステップと、周波数領域の信号に変換された被測定信号の、基本周波数近傍の周波数成分を抽出する帯域制限ステップと、帯域制限ステップにおいて抽出した周波数成分を、時間領域の信号に変換する時間領域変換ステップとを有してよい。
【００４８】
被測定信号の振幅変調成分を除去する波形クリップステップを更に備えてよい。
【００４９】
被測定信号の基本周期を算出する基本周期推定ステップを更に備えてよい。
【００５０】
周期ジッタ推定ステップは、周期ジッタ系列における高周波成分を除去する低域通過ステップを更に有してよい。
【００５１】
尚、上記の発明の概要は、本発明の必要な特徴の全てを列挙したものではなく、これらの特徴群のサブコンビネーションも又、発明となりうる。
【図面の簡単な説明】
【００５２】
【図１】従来の、タイムインターバルアナライザを用いた周期ジッタの測定を説明する。
【図２】タイムインターバルアナライザ１２が表示するヒストグラムの一例を示す。
【図３】従来の、デジタルオシロスコープを用いた周期ジッタの測定を説明する。
【図４】デジタルオシロスコープ１４の構成を示す。
【図５】図５（Ａ）は、被測定信号ｘ（ｔ）を示す。図５（Ｂ）は、周期ジッタを示す。
【図６】図６（Ａ）は、タイムインターバルアナライザによるジッタ測定結果を示す。図６（Ｂ）は、補間法によるジッタ測定結果を示す。
【図７】従来のΔφ法による周期ジッタの測定を説明する。
【図８】本発明の実施形態に係るジッタ測定装置１００の構成の一例を示す。
【図９】ゼロクロス点及び近似ゼロクロス点を示す。
【図１０】周期ジッタ系列を補正した場合と補正しない場合のジッタ測定結果を示す。
【図１１】図１１（Ａ）は、測定した周期ジッタのＲＭＳ値の一例を示す。図１１（Ｂ）は、測定した周期ジッタのピーク・ツゥ・ピーク値の一例を示す。
【図１２】１周期毎の立ち上がりの時間間隔を周期ジッタとして算出した場合の、被測定信号の周期ジッタ波形Ｊ［ｎ］を示す。
【図１３】図１３（Ａ）は、ジッタ測定装置１００が測定した周期ジッタのヒストグラムの一例を示す。図１３（Ｂ）は、従来のタイムインターバルアナライザで測定した周期ジッタのヒストグラムを示す。
【図１４】ジッタ測定装置１００が測定した周期ジッタのＲＭＳ値及びピーク値、並びに従来のタイムインターバルアナライザで測定したこれらの値とを示す。
【図１５】サイクル・ツゥ・サイクル周期ジッタ系列の一例を示す。
【図１６】ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８が算出したサイクル・ツゥ・サイクル周期ジッタのピーク値Ｊｃｃ、ｐｐ及びＲＭＳ値Ｊｃｃ、ＲＭＳの一例を示す。
【図１７】離散化された被測定信号ｘ（ｔ）の一例を示す。
【図１８】周波数領域の信号に変換された被測定信号のパワースペクトルの一例を示す。
【図１９】帯域制限された周波数領域の信号ＸＢＰ（ｆ）の一例を示す。
【図２０】帯域制限された時間領域の信号波形ｘＢＰ（ｔ）の一例を示す。
【図２１】解析信号の実数部の波形の一例を示す。
【図２２】ジッタ測定装置１００の動作の一例を示すフローチャートである。
【図２３】ジッタ測定装置１００の構成の他の例を示す。
【図２４】図２３において説明したジッタ測定装置１００の動作を説明するフローチャートである。
【図２５】位相誤差推定器１０４の構成の一例を示す。
【図２６】図２５において説明した位相誤差推定器１０４の動作の一例を示すフローチャートである。
【図２７】位相誤差推定器１０４の構成の他の例を示す。
【図２８】図２７において説明した位相誤差推定器１４０の動作の一例を示すフローチャートである。
【図２９】周期ジッタ推定器１０５の構成の一例を示す。
【図３０】図２９において説明した周期ジッタ推定器１０５の動作の一例を示すフローチャートである。
【図３１】周期ジッタ系列推定器１０５の構成の他の例を示す。
【図３２】図３１において説明した周期ジッタ系列推定器１０５の動作の一例を示すフローチャートである。
【図３３】帯域通過フィルタ１０３の構成の一例を示す。
【図３４】図３３において説明した帯域通過フィルタ１０３の動作の一例を示す。
【図３５】帯域通過フィルタ１０３の構成の他の例を示す。
【図３６】図３５において説明した帯域通過フィルタ１０３の動作の一例を示すフローチャートである。
【図３７】図３７（Ａ）は、ジッタ測定装置１００の構成例の一部を示す。図３７（Ｂ）は、ジッタ測定装置１００の構成例の一部を示す。図３７（Ｃ）は、ジッタ測定装置１００の構成例の一部を示す。
【図３８】振幅変調成分を有する被測定信号の一例を示す。
【図３９】波形クリッパ１９０１が波形クリップした被測定信号の一例を示す。
【図４０】図４０（Ａ）は、ＡＤコンバータ１７０１及び周期推定器１０１の動作の一例を示す。図４０（Ｂ）は、波形クリッパ１９０１及び周期推定器１０１の動作の一例を示す。
【図４１】周期ジッタ推定器１０５の構成の他の例を示す。
【図４２】図２９において説明した周期ジッタ推定器１０５の動作の他の例を示すフローチャートである。
【発明を実施するための最良の形態】
【００５３】
以下、発明の実施の形態を通じて本発明を説明するが、以下の実施形態は特許請求の範囲にかかる発明を限定するものではなく、又実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
【００５４】
図８は、本発明の実施形態に係るジッタ測定装置１００の構成の一例を示す。ジッタ測定装置１００は、周期性を有する被測定信号（入力信号）のジッタを測定する。ジッタ測定装置１００は、基本周期推定器１０１と、メモリ１０２と、帯域通過フィルタ１０３と、位相誤差推定器１０４と、周期ジッタ推定器１０５と、ジッタ検出器１０６とを備える。
【００５５】
基本周期推定器１０１は、被測定信号の基本周波数の周期を算出し、メモリ１０２は、基本周期推定器１０１が算出した被測定信号の基本周期を格納する。帯域通過フィルタ１０３は、被測定信号の所望の周波数成分を抽出し、位相誤差推定器１０４に供給する。本例において、帯域通過フィルタ１０３は、被測定信号の基本周波数近傍の周波数成分を抽出する。
【００５６】
位相誤差推定器１０４は、受け取った被測定信号の周波数成分に基づいて、被測定信号において所定の値をとる点に対して、予め定められた条件を満たす複数の標本点を抽出する。つまり、被測定信号は離散化されたデジタル信号であって、位相誤差推定器１０４は、離散化された複数の標本点から、予め定められた条件を満たす複数の標本点を抽出する。例えば、位相誤差推定器１０４は、被測定信号における複数のゼロクロス点に対して、最も近い振幅値を有する標本点をそれぞれ抽出する。このとき、位相誤差推定器１０４は、被測定信号の各サイクルの立ち上がり、又は各サイクルの立ち下がりのいずれかから、複数の標本点を抽出する。
【００５７】
位相誤差推定器１０４は、抽出した複数の標本点の位相と、被測定信号が当該所定の値をとる位相との位相誤差を算出する。このとき、位相誤差推定器１０４は、被測定信号の、同一の立ち上がり波形部分又は立ち下がり波形部分における、標本点の位相と、被測定信号が所定の値をとる位相との位相誤差を算出する。位相誤差推定器１０４は、被測定信号の各サイクル毎に、当該位相誤差を算出し、算出した位相誤差系列を周期ジッタ推定器１０５に供給する。また、位相誤差推定器１０４は、抽出したそれぞれの標本点間の時間間隔系列を周期ジッタ推定器１０５に供給する。
【００５８】
周期ジッタ推定器１０５は、位相誤差系列、時間間隔系列、及びメモリ１２が格納した基本周期に基づいて、被測定信号の周期ジッタ系列を算出する。まず、周期ジッタ推定器１０５は、位相誤差系列及び標本点間の時間間隔系列から、被測定信号の１サイクル毎の瞬時周期を算出し、算出した１サイクル毎の瞬時周期から基本周期を減算することにより、被測定信号の周期ジッタ系列を算出する。周期ジッタ推定器１０５における周期ジッタ系列の算出方法については後述する。
【００５９】
ジッタ検出器１０６は、周期ジッタ系列に基づいて、被測定信号のジッタを算出する。本例において、ジッタ検出器１０６は、周期ジッタ系列の最大値と最小値との差を算出するピーク・ツゥ・ピークジッタ検出器１０７と、周期ジッタ系列の二乗平均値を算出し、被測定信号のジッタの実効値（ＲＭＳ値）を算出するＲＭＳジッタ検出器１０８と、周期ジッタ系列のヒストグラムを生成するヒストグラム推定器１０９とを有する。
【００６０】
また、帯域通過フィルタ１０３は、アナログフィルタであってよく、またデジタルフィルタであってもよい。また、ＦＦＴ等を用いてソフトウェアで構成されてもよい。
【００６１】
次に、周期ジッタ推定器１０５における周期ジッタ系列の算出方法の一例について説明する。本例においては、被測定信号としてクロック信号を用いて説明するが、被測定信号はクロック信号に限定されない。つまり、ジッタ測定装置１００は、周期性を有する信号であれば、その信号のジッタを算出することができる。また、本例においては、周期ジッタ推定器１０５は、被測定信号のゼロクロス点の時間間隔の揺らぎに基づいて、被測定信号のジッタを算出する場合について説明する。
【００６２】
ジッタのないクロック信号は、基本周波数ｆ０を有する方形波である。この信号は、フーリエ変換によって周波数ｆ０、３ｆ０、５ｆ０、・・・からなる高調波に分解できる。ジッタは、被測定信号の基本周波数の揺らぎに対応するため、基本周波数近傍の周波数成分から、被測定信号のジッタを算出することができる。ジッタ測定装置１００は、帯域通過フィルタ１０３によって、被測定信号の基本周波数近傍の周波数成分を抽出するため、効率よくジッタを算出することができる。
【００６３】
ジッタをもつクロック信号の基本サイン波成分ｘ（ｔ）は、振幅をＡ、基本周期をＴ０とすると、
【数１】0004090989.tif 000002
で表される。ここで、φ（ｔ）は、被測定信号の瞬時位相を示し、基本周期Ｔ０を含むリニア瞬時位相成分２πｔ／Ｔ０と、初期位相成分θと、瞬時位相雑音成分Δφ（ｔ）との和で表される。
【００６４】
実信号ｘ（ｔ）の解析信号ｚ（ｔ）は、次式の複素信号で定義される。
【数２】0004090989.tif 000003
ここで、ｊは虚数単位であり、複素信号ｚ（ｔ）の虚数部
【数３】0004090989.tif 000004
は、実数部ｘ（ｔ）のヒルベルト変換である。一方、時間波形ｘ（ｔ）のヒルベルト変換は、次式で定義される。
【数４】0004090989.tif 000005
【００６５】
ここで、
【数３】0004090989.tif 000006
は式（３）に示すように関数ｘ（ｔ）と、１／πｔとの畳み込み演算で与えられる。すなわち、ヒルベルト変換は、ｘ（ｔ）を全帯域通過フィルタを通過させた場合のフィルタの出力と等価である。
【００６６】
ただし、このときの出力は、スペクトル成分の大きさは変わらないが、その位相はπ／２だけシフトする。解析信号及びヒルベルト変換については、例えば、A.Papoulis,Probability,Ramdom Variables,and Stochastic Processes,2nd edition,McGraw-Hill Book Company,1984.に記載されている。
【００６７】
従って、被測定信号の基本コサイン波ｘ（ｔ）の解析信号は、
【数５】0004090989.tif 000007
で与えられる。解析信号ｚ（ｔ）の瞬時周波数は、下式で与えられる。
【数６】0004090989.tif 000008
すなわち、
【数７】0004090989.tif 000009
より、周期ジッタＪは次式で算出できる。
【数８】0004090989.tif 000010
【００６８】
一方、ｘ（ｔ）の瞬時角周波数ω（ｔ）は、
【数９】0004090989.tif 000011
で与えられる。従って、式（７）及び式（８）より、周期ジッタＪ（ｔ）は、瞬時角周波数ω（ｔ）を用いて以下のように表される。
【数１０】0004090989.tif 000012
【００６９】
位相誤差推定器１０４は、被測定信号ｘ（ｔ）が所定の値をとる点と、当該所定の値をとる点と同一の立ち上がり波形部分又は立ち下がり波形部分における標本点との位相誤差系列を算出する。本例においては、図９に示すように、位相誤差推定器１０４は、被測定信号ｘ（ｔ）の立ち上がりの各ゼロクロス点ｔｋ、ｔｋ＋１に、最も近い標本点ｔ［ｋ］、ｔ［ｋ＋１］を抽出し、各ゼロクロス点と対応する標本点との位相誤差をそれぞれ算出する。
【００７０】
以下、ゼロクロス点を例にして説明するが、位相誤差推定器１０４が被測定信号ｘ（ｔ）の所定の値をとる点について位相誤差系列を算出する場合についても、同様に位相誤差系列は算出できる。
【００７１】
近似ゼロクロス点ｔ［ｋ］と、被測定信号におけるゼロクロス点ｔｋとの間の位相誤差をδ［ｋ］とすると、隣り合う二つの近似ゼロクロス点間［ｔ［ｋ］、ｔ［ｋ＋１］］における瞬時角周波数ω［ｋ］は、下式に示すように［ｔ［ｋ］、ｔ［ｋ＋１］］間にすすむ位相角と経過時間との比で与えられる。
【数１１】0004090989.tif 000013
ここで、Ｔｋ、ｋ＋１は隣り合う二つの近似ゼロクロス点ｔ［ｋ］とｔ［ｋ＋１］との間の時間間隔であり、近似ゼロクロス点のタイミングは、標本点のサンプリングタイミングであるため、近似ゼロクロス点のタイミング系列ｔ［ｋ］は既知であるか、又は容易に算出できる。よって、Ｔｋ。ｋ＋１は、近似ゼロクロス点のタイミング系列において、対応するタイミングの差分を算出することにより、下式のように算出できる。位相誤差推定器１０４は、下式に基づいて隣り合う近似ゼロクロス点間の時間間隔系列を算出する。
【数１２】0004090989.tif 000014
【００７２】
また、位相誤差δ［ｋ］は、近似ゼロクロス点ｔ［ｋ］における位相角φ（ｔ［ｋ］）とゼロクロス点ｔｋにおける位相角φ（ｔｋ）により、下式のように表される。
【数１３】0004090989.tif 000015
ここで、帯域制限された被測定信号ｘ（ｔ）は、式（１）に示すように振幅Ａを有する基本コサイン波であり、位相角φ（ｔ）と振幅値ｘ（ｔ）は次のような関係を有する。
【数１４】0004090989.tif 000016
【００７３】
このため、被測定信号の時刻ｔにおける位相角φ（ｔ）は、時刻ｔの振幅値ｘ（ｔ）と振幅Ａから、逆正弦関数を用いて、
【数１５】0004090989.tif 000017
で表される。ただし、逆正弦関数は、−π／２から＋π／２までの範囲の位相を用いて表されるため、式（１４）はコサイン波の位相角が−πから０までの範囲、即ち、コサイン波の立ち上がり領域において成り立つ。従って、基本コサイン波ｘ（ｔ）の立ち上がりゼロクロス点に対して、近似ゼロクロス点ｔ［ｋ］における振幅値を
【数１６】0004090989.tif 000018
とすると、ｔ［ｋ］における位相角φ（ｔ［ｋ］）は、式（１４）より、
【数１７】0004090989.tif 000019
で与えられる。一方、ゼロクロス点ｔｋにおける振幅値は０であるから、ゼロクロス点ｔｋにおける位相角φ（ｔｋ）は、
【数１８】0004090989.tif 000020
である。従って、立ち上がりゼロクロス点ｔｋとその近似ゼロクロス点ｔ［ｋ］との位相誤差δ［ｋ］は、式（１２）、式（１６）、及び式（１７）より、
【数１９】0004090989.tif 000021
で算出することができる。、コサイン波の位相角が０からπまでの範囲、即ち立ち下がり領域については、時刻ｔにおける位相角φ（ｔ）は、時刻ｔの振幅値ｘ（ｔ）と振幅Ａから、逆正弦関数を用いて、
【数２０】0004090989.tif 000022
で表される。従って、基本コサイン波ｘ（ｔ）の立ち下がりにおいて、近似ゼロクロス点とゼロクロス点との位相誤差δ［ｋ］は、次式で算出できる。
【数２１】0004090989.tif 000023
【００７４】
位相誤差推定器１０４は、ゼロクロス点が立ち上がり波形部分にあるか、又は立ち下がり波形部分にあるかを判定し、対応する式（１８）又は式（１９）のいずれかを用いて、ゼロクロス点と近似ゼロクロス点との位相誤差を算出してよい。以上から明らかなように、位相誤差δ［ｋ］は、近似ゼロクロス点における被測定信号の瞬時値と、被測定信号の振幅値とから算出することができる。
【００７５】
また、位相誤差推定器１０４は、被測定信号ｘ（ｔ）を解析信号に変換し、位相誤差δ［ｋ］を算出してもよい。すなわち、位相誤差推定器１０４は、下式のように、被測定信号を解析信号に変換する。
【数２２】0004090989.tif 000024
そして、近似ゼロクロス点ｔ［ｋ］における解析信号データ
【数２３】0004090989.tif 000025
から、下式のように近似ゼロクロス点における瞬時位相を算出する。
【数２４】0004090989.tif 000026
位相誤差推定器１０４は、算出した近似ゼロクロス点における瞬時位相と、ゼロクロス点における瞬時位相との差分を算出することにより、位相誤差δ［ｋ］を算出してもよい。例えば、立ち上がりゼロクロス点におけるゼロクロス点の瞬時位相は−π／２であるため、下式によって位相誤差δ［ｋ］を算出してよい。
【数２５】0004090989.tif 000027
【００７６】
また、立ち下がりゼロクロス点におけるゼロクロス点の瞬時位相はπ／２であるため、下式によって位相誤差δ［ｋ］を算出してもよい。
【数２６】0004090989.tif 000028
【００７７】
以上により、周期ジッタＪ［ｋ］は、位相誤差δ［ｋ］を用いて、式（９）、及び式（１０）から
【数２７】0004090989.tif 000029
で算出することができる。ここで、右辺の（  ）で囲まれた項は、図９に示すように、被測定信号の瞬時周期
【数２８】0004090989.tif 000030
と基本周期Ｔ０との差、即ち周期ジッタを示す。周期ジッタ推定器１０５は、式（２７）を用いて周期ジッタを算出する。
【００７８】
また、（Ｔ０／Ｔｋ、ｋ＋１）の項は、上記周期ジッタに対する補正項である。この補正により、より高い精度で周期ジッタを算出することができる。つまり、周期ジッタ推定器１０５は、近似ゼロクロス点のタイミングと被測定信号の基本周期とに基づいて、算出した周期ジッタを補正してよい。
【００７９】
上記補正項は、図１０に示すように、周期ジッタのＲＭＳ値とピーク・ツゥ・ピーク値に対して、ジッタ推定値と理論値との誤差を小さくできる。つまり、補正項を用いない場合は、ＲＭＳ値及びピーク・ツゥ・ピーク値はそれぞれ＋０．５３％及び１５．９％の誤差があったが、補正項を用いた例においては、これらの誤差をそれぞれ−０．００４％及び＋０．０４％に低減できた。尚、図１０は、サイン波ジッタを有する信号に対するシミュレーションにより算出した。
【００８０】
また、実波形を用いた実験結果を図１１（Ａ）に示す。図１１（Ａ）において横軸は１周期当たりのサンプリング数を示し、縦軸は周期ジッタのＲＭＳ値を示す。図１１（Ａ）に示すように、周期ジッタのＲＭＳ値はサンプリング数に関係なくほぼ一定の値を算出することができる。また、補正値を用いる方法も、補正値を用いない方法も、同様の結果が得られる。
【００８１】
また、図１１（Ｂ）に周期ジッタのピーク・ツゥ・ピーク値の実験結果を示す。補正項を用いない場合、周期ジッタのピーク・ツゥ・ピーク値は、過剰標本化比が小さくなるにつれて過大評価される。ここで過剰標本化比は、サンプリング定理により波形を再現するために必要な標本数（再現するべき波形の１周期に対して標本数２）に対するサンプリング数の比を指す。図１１（Ｂ）からわかるように、ジッタ測定装置１００によれば、過剰標本化比が小さい場合に生じる誤差を低減することができる。本例の実験結果によれば、過剰標本化比が４のとき約８％、過剰標本化比が１．５のとき約１８％の誤差が生じるが、上述した補正項を用いることにより、この誤差を低減することができる。
【００８２】
また、ジッタ測定装置１００は、はじめに被測定信号の基本周波数成分を取り出すために、帯域通過フィルタ１０３を用いて被測定信号を帯域制限する。この帯域制限された被測定信号の各ゼロクロス点に最も近い標本点を近似ゼロクロス点として抽出し、この近似ゼロクロスサンプル列から近似ゼロクロス点と被測定信号のゼロクロス点との位相誤差系列（δ［ｎ］）を算出する。次に位相誤差系列から式（２６）によって被測定信号の周期ジッタ系列を算出する。
【００８３】
ここで、ジッタ測定装置１００は、所望の立ち上がり波形又は立ち下がり波形において、周期ジッタを算出してよい。例えば、被測定信号の立ち上がりと、隣接する立ち下がりとの間の時間間隔の変動を周期ジッタとして算出してよく、また被測定信号の所望の周期間隔における立ち上がり（又は立ち下がり）の周期ジッタを算出してもよい。例えば、２サイクル毎における立ち上がりの時間間隔の変動を周期ジッタとして算出してもよい。
【００８４】
次に、周期ジッタ推定器１０５は、算出した周期ジッタ系列に前述した補正項を乗算し、周期ジッタ系列の値を補正する。図１２に、１サイクル毎の立ち上がりの時間間隔を周期ジッタとして算出した場合の、被測定信号の周期ジッタ波形Ｊ［ｎ］を示す。ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８は、周期ジッタ推定器１０５がこのようにして算出した周期ジッタ系列の二乗平均、及び最大値と最小値との差をそれぞれ算出する。ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８は、それぞれ式（２８）、及び式（２９）により周期ジッタ系列のＲＭＳ値及びピーク・ツゥ・ピーク値を算出する。
【数２９】0004090989.tif 000031
【数３０】0004090989.tif 000032
【００８５】
ここで、Ｍは測定された周期ジッタのデータ数である。図１３（Ａ）に、ジッタ測定装置１００が測定した周期ジッタのヒストグラムの一例を示す。ヒストグラム推定器１０９は、図１３（Ａ）に示すような周期ジッタのヒストグラムを生成する。また、図１３（Ｂ）に、従来のタイムインターバルアナライザで測定した周期ジッタのヒストグラムを示す。図１３（Ａ）及び図１３（Ｂ）に示すように、本例におけるジッタ測定装置１００によれば、従来のタイムインターバルアナライザと互換性の有る測定を行うことができる。
【００８６】
また、図１４にジッタ測定装置１００が測定した周期ジッタのＲＭＳ値及びピーク値、並びに従来のタイムインターバルアナライザで測定したこれらの値とを示す。ここで、測定された周期ジッタのピーク値Ｊｐｐは、イベント数（ゼロクロス数）の対数の平方根にほぼ比例する。例えば、５０００イベント程度においては、Ｊｐｐ＝４５ｐｓとなる。図１４に、Ｊｐｐ＝４５ｐｓを正値とした場合における、測定装置１００の測定結果の誤差を示す。図１３及び図１４に示すように、測定装置１００によれば、ジッタ測定の誤差が小さく、従来のタイムインターバルアナライザと互換性のある測定を行うことができる。
【００８７】
更に、測定装置１００によれば、周期ジッタと、サイクル・ツゥ・サイクル周期ジッタとを同時に測定することができる。サイクル・ツゥ・サイクル周期ジッタは、被測定信号において連続するサイクル間の周期変動であり、下式により算出することができる。
【数３１】0004090989.tif 000033
【００８８】
ジッタ測定装置１００は、上式に基づいて、サイクル・ツゥ・サイクル周期ジッタを算出してよい。また、ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８は、サイクル・ツゥ・サイクル周期ジッタに基づいて、サイクル・ツゥ・サイクル周期ジッタのピーク・ツゥ・ピーク値及びＲＭＳ値を算出してもよい。この場合、、ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８は下式に基づいてサイクル・ツゥ・サイクル周期ジッタのピーク・ツゥ・ピーク値Ｊｃｃ、ｐｐ及びＲＭＳ値Ｊｃｃ、ＲＭＳを算出する。
【数３２】0004090989.tif 000034
【数３３】0004090989.tif 000035
ここで、Ｌは測定されたサイクル・ツゥ・サイクル周期ジッタの標本数を示す。図１５に、サイクル・ツゥ・サイクル周期ジッタ系列の一例を示す。また、ピーク・ツゥ・ピークジッタ検出器１０７及びＲＭＳジッタ検出器１０８が算出したサイクル・ツゥ・サイクル周期ジッタのピーク値Ｊｃｃ、ｐｐ及びＲＭＳ値Ｊｃｃ、ＲＭＳの一例を図１６に示す。
【００８９】
また、ジッタ測定装置１００は、アナログの被測定信号を離散化した後で帯域制限フィルタ処理を適用してよく、アナログの被測定信号に帯域制限フィルタ処理を適用してから波形を離散化してもよい。帯域通過フィルタ１０３は、アナログフィルタであってよく、またデジタルフィルタであってもよい。また、フーリエ変換を用いてソフトウェアによりフィルタ処理を行ってもよい。アナログ信号の離散化には、高速なＡＤ変換器、デジタイザ、デジタルサンプリングオシロスコープ等を用いることが好ましい。
【００９０】
また、ジッタ測定装置１００は、波形クリップ手段を有し、被測定信号のジッタに対応する位相変調成分を保持した状態で、振幅変調成分を取り除いて、ジッタの測定を行ってもよい。振幅変調成分を除去することにより、精度よくジッタの測定を行うことができる。
【００９１】
次に、帯域通過フィルタ１０３について説明する。本例において、帯域通過フィルタ１０３は、デジタル信号の帯域制限を行う。離散化されたデジタル信号の帯域制限は、デジタルフィルタでも実現してよく、またフーリエ変換によって実現してもよい。ここでは、フーリエ変換を用いて帯域制限を行う場合について説明する。
【００９２】
図１７に、離散化された被測定信号ｘ（ｔ）の一例を示す。帯域通過フィルタ１０３は、離散化された被測定信号を、高速フーリエ変換（ＦＦＴ）により周波数領域の信号Ｘ（ｆ）に変換する。図１８に周波数領域の信号に変換された被測定信号のパワースペクトルの一例を示す。次に、被測定信号の基本周波数（本例においては４００ＭＨｚ）近傍のデータのみを残し、他の周波数成分を零とすることにより、帯域制限を行う。
【００９３】
図１９に、帯域制限された周波数領域の信号ＸＢＰ（ｆ）の一例を示す。最後に、帯域制限された信号ＸＢＰ（ｆ）に逆ＦＦＴを施し、帯域制限された時間領域の信号波形ｘＢＰ（ｔ）を算出する。以上の動作により、帯域通過フィルタ１０３は、被測定信号の基本周波数近傍の成分のみを抽出した信号を生成することができる。図２０に、帯域制限された時間領域の信号波形ｘＢＰ（ｔ）の一例を示す。
【００９４】
次に、位相誤差推定器１０４における、近似ゼロクロス点の検出法について説明する。はじめに、位相誤差推定器１０４は、入力された被測定信号の解析信号の実数部ｘ（ｔ）の最大値を１００％レベル、最小値を０％レベルとして、検出するべきゼロクロス点の信号値レベルＶｚ（５０％）を算出する。次に、被測定信号の隣り合うそれぞれのサンプル値（ｘ（ｊ−１）、ｘ（ｊ）、ただしｊは整数）と、算出した信号値Ｖｚとの差分（ｘ（ｊ−１）−Ｖｚ、ｘ（ｊ）−Ｖｚ）を算出する。次に、算出した差分の積（（ｘ（ｊ−１）−Ｖｚ）×（ｘ（ｊ）−Ｖｚ））を算出する。被測定信号のゼロクロス点は、算出した積の値の符号が変化する点である。
【００９５】
ここで、位相誤差推定器１０４は、被測定信号のそれぞれの標本点に対して、上述した演算を行ってよい。演算した結果に基づいて、前述した積の値の符号が変化する標本点の組を抽出する。この場合、位相誤差推定器１０４は、前述した積の値が、例えば正から負に変化する標本点の組を検出してよく、また負から正に変化する標本点の組を検出してもよい。
【００９６】
位相誤差推定器１０４は、抽出した標本点の組に対して、それぞれの標本点のうち絶対値の小さい標本点を近似ゼロクロス点として抽出する。また、抽出した近似ゼロクロス点におけるタイミングを、近似ゼロクロスタイミングとして、近似ゼロクロス点間の時間間隔系列を算出する。図２１に、解析信号の実数部の波形の一例を示す。図２１において、丸印は、検出された近似ゼロクロス点を示す。本例においては、位相誤差推定器１０４は、波形の立ち上がりにおける近似ゼロクロス点を検出する。
【００９７】
また、位相誤差推定器１０４において近似ゼロクロス点を算出する被測定信号は、近似ゼロクロス点を精度よく算出するために、オーバーサンプリングされたデジタル信号であることが好ましい。例えば、基本周波数の少なくとも３倍以上の周波数で、離散化されていることが好ましい。例えば、帯域通過フィルタ１０３がアナログ信号を位相誤差推定器１０４に供給する場合、位相誤差推定器１０４は、当該アナログ信号の基本周波数の少なくとも３倍以上の周波数で、当該アナログ信号をサンプリングした標本点を算出する。
【００９８】
図２２は、ジッタ測定装置１００の動作の一例を示すフローチャートである。まず、ステップ２０１において基本周期推定器１０１が、被測定信号の基本周期を算出し、メモリ１０２に格納する。次に、ステップ２０２において帯域通過フィルタ１０３が、被測定信号の基本周波数近傍の周波数成分を抽出し、位相誤差推定器１０４に供給する。
【００９９】
次に、ステップ２０３において位相誤算推定器１０４が、帯域通過フィルタ１０３が抽出した被測定信号の周波数成分において、前述した近似ゼロクロス点の系列を算出し、それぞれの近似ゼロクロス点と、対応するゼロクロス点との位相誤差δ［ｋ］を算出する。また、位相誤差推定器１０４は、前述したように隣り合う近似ゼロクロス点の間の時間間隔系列を算出する。
【０１００】
次に、ステップ２０４において周期ジッタ推定器１０５が、位相誤差推定器１０４が算出した近似ゼロクロス点の系列、及び近似ゼロクロス点の時間間隔系列と、メモリ１０２が格納した被測定信号の基本周期とに基づいて、被測定信号の周期ジッタ系列を式（２６）に関連して説明したように算出する。そして、ステップ２０５においてジッタ検出器１０６が、位相誤差推定器１０５が算出した周期ジッタ系列に基づいて、被測定信号の周期ジッタ値を算出する。
【０１０１】
また、被測定信号の基本周期が、メモリ１０２に予め格納されている場合、ステップ２０１は省略できる。
【０１０２】
図２３は、ジッタ測定装置１００の構成の他の例を示す。図２３において、図８と同一の符号を付した構成要素は、図８に関連して説明した構成要素と同一又は同様の機能及び構成を有する。本例におけるジッタ測定装置１００は、被測定信号のサイクル・ツゥ・サイクル周期ジッタ値を更に算出する。本例におけるジッタ測定装置１００は、図８に関連して説明したジッタ測定装置１００の構成に加え、サイクル・ツゥ・サイクル周期ジッタ推定器３０１とスイッチ３０２とを更に備える。
【０１０３】
サイクル・ツゥ・サイクル周期ジッタ推定器３０１は、周期ジッタ推定器１０５が算出した周期ジッタ系列を受け取り、式（３０）に関連して説明したように、周期ジッタ系列の隣り合うデータの差分を算出し、被測定信号のサイクル・ツゥ・サイクル周期ジッタ系列を算出する。
【０１０４】
スイッチ３０２は、周期ジッタ推定器１０５が算出した周期ジッタ系列、又はサイクル・ツゥ・サイクル周期ジッタ推定器３０１が算出したサイクル・ツゥ・サイクル周期ジッタ系列のいずれかを、ジッタ検出器１０６に供給する。ジッタ検出器１０６は、受け取った周期ジッタ系列、又はサイクル・ツゥ・サイクル周期ジッタ系列に基づいて、被測定信号の周期ジッタ値、又はサイクル・ツゥ・サイクル周期ジッタ値を算出する。
【０１０５】
また、ジッタ測定装置１００は、スイッチ３０２を備えなくともよい。この場合、ジッタ測定装置１００は、サイクル・ツゥ・サイクル周期ジッタ値を算出し、周期ジッタ値を算出しない。
【０１０６】
図２４は、図２３において説明したジッタ測定装置１００の動作を説明するフローチャートである。図２４において、ステップ２０１〜ステップ２０５においては、ジッタ測定装置１００は図２２において説明したジッタ測定装置１００と同一又は同様の動作を行う。
【０１０７】
ステップ４０１において、サイクル・ツゥ・サイクル周期ジッタ推定器３０１は、サイクル・ツゥ・サイクル周期ジッタ系列を算出する。そして、スイッチ３０２は、サイクル・ツゥ・サイクル周期ジッタ系列をジッタ検出器１０６に供給する。
【０１０８】
次にステップ４０２において、ジッタ検出器１０６は、サイクル・ツゥ・サイクル周期ジッタ系列に基づいて、被測定信号のサイクル・ツゥ・サイクル周期ジッタ値を算出する。また、ジッタ測定装置１００が被測定信号の周期ジッタ値を算出しない場合、ステップ２０５は省略できる。
【０１０９】
図２５は、位相誤差推定器１０４の構成の一例を示す。位相誤差推定器１０４は、ゼロクロスサンプラ５０１と、位相誤差算出器５０２とを有する。
【０１１０】
ゼロクロスサンプラ５０１は、帯域通過フィルタ１０３から帯域制限された被測定信号を受け取り、上述した被測定信号の複数の近似ゼロクロス点を抽出する。例えば、ゼロクロスサンプラ５０１は、被測定信号を所定のサンプリング周波数でサンプリングし、サンプリングした標本点から、近似ゼロクロス点を抽出する。また、ゼロクロスサンプラ５０１は抽出した近似ゼロクロス点の時間間隔系列を算出する。ゼロクロスサンプラ５０１は、抽出した近似ゼロクロス点の系列を位相誤差算出器５０２に供給し、算出した時間間隔系列を周期ジッタ推定器１０５に供給する。
【０１１１】
位相誤差算出器５０２は、上述したように、それぞれの近似ゼロクロス点と、対応するゼロクロス点との位相誤差系列を算出し、周期ジッタ推定器１０５に供給する。
【０１１２】
図２６は、図２５において説明した位相誤差推定器１０４の動作の一例を示すフローチャートである。まず、ステップ６０１において、ゼロクロスサンプラ５０１が、帯域通過フィルタ１０３から受け取った信号に基づいて、近似ゼロクロス点の系列を算出する。
【０１１３】
次に、ステップ６０２において、ゼロクロスサンプラ５０１は、近似ゼロクロス点の系列における時間間隔系列を算出する。ここで、ゼロクロスサンプラ５０１は、式（１１）を用いて近似ゼロクロス点の系列における時間間隔系列を算出する。
【０１１４】
次に、ステップ６０３において、位相誤差算出器５０２が、近似ゼロクロス点の系列に基づいて、近似ゼロクロス点の系列と、被測定信号のゼロクロス点との位相誤差の系列を算出する。ここで、位相誤差算出器５０２は、式（１８）又は式（２０）を用いて当該位相誤差の系列を算出する。帯域通過フィルタ１０３によって、被測定信号を基本周波数近傍の周波数成分を有する正弦波とするため、式（１８）又は式（２０）を用いて位相誤差を算出することができる。また、式（１８）及び式（２０）における振幅Ａは、帯域通過フィルタ１０３が抽出する各サンプルの二乗平均値又は最大値と最小値との差から算出することができる。
【０１１５】
図２７は、位相誤差推定器１０４の構成の他の例を示す。本例における位相誤差推定器１０４は、図２５において説明した位相誤差推定器１０４の構成に加え、解析信号変換器７０１を更に有する。解析信号変換器７０１は、帯域通過フィルタ１０３から受け取った信号を、複素数の解析信号に変換し、ゼロクロスサンプラ５０１に供給する。例えば、解析信号変換器７０１は、受け取った信号を解析信号の実数部として出力し、受け取った信号のヒルベルト変換を解析信号の虚数部として出力する。また、解析信号変換器７０１は、フーリエ変換によって解析信号を生成してもよい。
【０１１６】
ゼロクロスサンプラ５０１は、解析信号の実数部における近似ゼロクロス点のタイミングにおいて、解析信号をサンプリングし、式（２２）に示すような近似ゼロクロス複素データ系列を算出する。また、位相誤差算出器５０２は、近似ゼロクロス複素データ系列から、式（２３）を用いて近似ゼロクロス点の位相を算出し、式（２４）、及び式（２５）を用いて、位相誤差系列δ［ｋ］を算出する。また、ゼロクロスサンプラ５０１は、式（１１）を用いて近似ゼロクロス点間の時間間隔系列を算出する。
【０１１７】
図２８は、図２７において説明した位相誤差推定器１４０の動作の一例を示すフローチャートである。まず、ステップ８０１において、解析信号変換器７０１が、帯域通過フィルタ１０３から受け取った信号を解析信号に変換する。次に、ステップ８０２において、ゼロクロスサンプラ５０１が、解析信号の実数部における近似ゼロクロスタイミングにおいて、解析信号をサンプリングし、近似ゼロクロス複素データ系列を算出する。
【０１１８】
次に、ステップ８０３において、ゼロクロスサンプラ５０１は、近似ゼロクロス複素データ系列における時間間隔系列を算出する。次に、ステップ８０４において、位相誤差算出器５０２が、近似ゼロクロス複素データ系列から、近似ゼロクロス複素データと、被測定信号のゼロクロス点との位相誤差系列を算出する。また、ゼロクロスサンプラ５０１は、近似ゼロクロス点間の時間間隔系列を算出する。
【０１１９】
図２９は、周期ジッタ推定器１０５の構成の一例を示す。周期ジッタ推定器１０５は、瞬時周期推定器９０１と、基本周期減算器９０２と、周期ジッタ系列補正器９０３とを有する。
【０１２０】
瞬時周期推定器９０１は、位相誤差推定器１０４が算出した位相誤差系列、及び近似ゼロクロス点間の時間間隔系列、並びにメモリ１０２が格納した被測定信号の基本周期に基づいて、式（２７）に関連して説明したように、被測定信号の瞬時周期系列を算出する。
【０１２１】
基本周期減算器９０２は、瞬時周期推定器９０１が算出した瞬時周期系列のそれぞれの周期値から、被測定信号の基本周期値を減算した周期ジッタ系列を算出する。
【０１２２】
周期ジッタ系列補正器９０３は、基本周期減算器９０２が算出した周期ジッタ系列を、式（２６）に関連して説明した補正項を用いて補正して出力する。
【０１２３】
図３０は、図２９において説明した周期ジッタ推定器１０５の動作の一例を示すフローチャートである。まず、ステップ１００１において、上述したように瞬時周期推定器９０１が、被測定信号の瞬時周期系列を算出する。次に、ステップ１００２において、基本周期減算器９０２が、瞬時周期系列のそれぞれの瞬時周期値から、被測定信号の基本周期を減算し、周期ジッタ系列を生成する。次に、ステップ１００３において、周期ジッタ系列補正器９０３が、上述したように周期ジッタ系列を補正して出力する。
【０１２４】
図３１は、周期ジッタ系列推定器１０５の構成の他の例を示す。周期ジッタ推定器１０５は、瞬時角周波数推定器１１０１と、周期ジッタ算出器１１０２とを有する。本例における周期ジッタ系列推定器１０５は、それぞれの近似ゼロクロスタイミング間における、被測定信号の瞬時角周波数に基づいて、被測定信号の周期ジッタ系列を算出する。
【０１２５】
瞬時角周波数推定器１１０１は、位相誤差推定器１０４から受け取った位相誤差系列と、近似ゼロクロス点間の時間間隔系列とから、式（１０）に関連して説明したように、それぞれの近似ゼロクロスタイミング間における、被測定信号の瞬時角周波数を算出する。
【０１２６】
周期ジッタ算出器１１０２は、それぞれの近似ゼロクロスタイミング間における、被測定信号の瞬時角周波数と、被測定信号の基本周期とに基づいて、式（９）を用いて周期ジッタ系列を算出する。
【０１２７】
また、周期ジッタ推定器１０５は、周期ジッタ算出器１１０２が算出した周期ジッタ系列において、低周波成分を通過させる低域通過フィルタ２１０１を更に有することが好ましい。
【０１２８】
図３２は、図３１において説明した周期ジッタ系列推定器１０５の動作の一例を示すフローチャートである。まず、ステップ１２０１において、瞬時角周波数推定器１１０１が、上述したように被測定信号の瞬時角周波数系列を算出する。次に、ステップ１２０２において、周期ジッタ算出器１１０２は、瞬時角周波数系列及び被測定信号の基本周期に基づいて、被測定信号の周期ジッタ系列を算出する。
【０１２９】
図３３は、帯域通過フィルタ１０３の構成の一例を示す。帯域通過フィルタ１０３は、周波数領域変換器１３０１と、帯域制限器１３０２と、時間領域変換器１３０３とを有する。
【０１３０】
周波数領域変換器１３０１は、被測定信号を例えばＦＦＴによって、周波数領域の信号に変換する。帯域制限器１３０２は、周波数領域の信号に変換された被測定信号の、基本周波数近傍の周波数成分を抽出して、時間領域変換器１３０３に供給する。時間領域変換器１３０３は、受けとった周波数領域の信号を、例えば逆ＦＦＴによって、時間領域の信号に変換し、位相誤差推定器１０４に供給する。
【０１３１】
図３４は、図３３において説明した帯域通過フィルタ１０３の動作の一例を示す。まず、ステップ１４０１において周波数領域変換器１３０１が被測定信号を周波数領域の信号に変換する。次に、ステップ１４０２において帯域制限器１３０２が、受け取った周波数領域の信号のうち、被測定信号の基本周波数近傍の周波数成分を抽出する。次に、ステップ１４０３において時間領域変換器１３０３が、帯域通過フィルタが抽出した周波数成分を、時間領域の信号に変換する。
【０１３２】
図３５は、帯域通過フィルタ１０３の構成の他の例を示す。帯域通過フィルタ１０３は、被測定信号を蓄積するバッファメモリ１５０１と、バッファメモリ１５０１が蓄積した信号を順次選択する信号選択器１５０２と、信号選択器１５０２が選択した部分信号に、予め定められた窓関数を乗算する窓関数乗算器１５０３と、窓関数が乗算された部分信号を、周波数領域の信号に変換する周波数領域変換器１５０４と、周波数領域に変換された信号のうち、被測定信号の正の基本周波数近傍の周波数成分を抽出する帯域制限器１５０５と、帯域制限器１５０５が抽出した周波数成分を、時間領域の信号に変換する時間領域変換器１５０６と、時間領域の信号に、窓関数の逆数を乗算し、解析信号を生成する振幅補正器１５０７とを有する。また、信号選択器１５０２は、直前に選択した部分信号と、一部が重複する部分信号を選択する。
【０１３３】
被測定信号ｘ（ｔ）の部分信号に、窓関数を乗算した場合、被測定信号ｘ（ｔ）は振幅変調される。しかし、本例における帯域通過フィルタ１０３によれば、振幅補正器１５０８において窓関数の逆数を乗算することにより、被測定信号ｘ（ｔ）の振幅変調を補正することができる。
【０１３４】
図３６は、図３５において説明した帯域通過フィルタ１０３の動作の一例を示すフローチャートである。はじめに、ステップ１６０１において、バッファメモリ１５０１は、被測定信号を蓄積する。つぎに、ステップ１６０２において、信号選択器１５０２は、バッファメモリ１５０１が蓄積した信号の一部を取り出す。つぎに、ステップ１６０３において、窓関数乗算器１５０３は、信号選択器１５０２が選択した部分信号に窓関数を乗算する。つぎに、ステップ１６０４において、周波数領域変換器１５０４は、窓関数を乗算された部分信号にFFTを施し、時間領域の信号を周波数領域の両側スペクトルに変換する。つぎに、ステップ１６０５において、帯域制限器１５０６は、変換された周波数領域の両側スペクトルに対し、負の周波数成分をゼロに置き換える。つぎに、帯域制限器１５０６は、負の周波数成分をゼロに置き換えられた片側スペクトルに対し、被測定信号の基本周波数付近の成分のみを残しその他の周波数成分をゼロに置き換え、周波数領域の信号を帯域制限する。つぎに、ステップ１６０６において、時間領域変換器１５０７は、帯域制限された周波数領域の片側スペクトルに逆FFTを施し、周波数領域の信号を時間領域の信号に変換する。つぎに、ステップ１６０７において、振幅補正器１５０８は、逆変換された時間領域の信号に窓関数の逆数を乗算し、帯域制限された解析信号を生成する。次に、ステップ１６０８において、帯域通過フィルタ１０３は、バッファメモリ１５０１に処理されていないデータが存在するか否かを確認し、処理されていないデータが存在する場合、ステップ１６０９において、信号選択器１５０２は、次の部分信号を選択する。信号選択器１０２が、前回取り出した分と一部重複させながら部分信号を順次取り出した後、帯域通過フィルタ１０３は、上述した処理を繰り返す。
【０１３５】
図３７（Ａ）は、ジッタ測定装置１００の構成例の一部を示す。図３７（Ａ）に示すように、ジッタ測定装置１００は、アナログの被測定信号をサンプリングして、デジタルの被測定信号を生成し、基本周期推定器１０１及び帯域通過フィルタ１０３に供給するＡ／Ｄ変換器１７０１を更に備えてもよい。Ａ／Ｄ変換器１７０１は、高速なサンプリングが行える、Ａ／Ｄ変換器、デジタイザ、デジタルサンプリングオシロスコープであることが好ましい。
【０１３６】
図３７（Ｂ）は、ジッタ測定装置１００の構成例の一部を示す。図３７（Ｂ）に示すように、ジッタ測定装置１００は、アナログの被測定信号を波形クリップし、振幅変調成分を除去する波形クリッパ１９０１を更に備えてもよい。また、図３７Ｃに示すように、ジッタ測定装置１００は、波形クリッパ１９０１の前段又は後段にＡ／Ｄ変換器１７０１を更に備えてもよい。
【０１３７】
次に、波形クリッパ１９０１における波形クリップについて説明する。波形クリッパ１９０１は、被測定信号から振幅変調成分を除去し、周期ジッタに対応する位相変調成分を残す。例えば、波形クリッパ１９０１は、アナログ又はデジタルの被測定信号の値を所定の定数倍し、予め定められた第１の閾値より大きい信号値を第１の閾値に置き換え、予め定められた第２の閾値より小さい信号値を第２の閾値に置き換えることにより、被測定信号の振幅変調成分を除去する。
【０１３８】
図３８に振幅変調成分を有する被測定信号の一例を示す。被測定信号が振幅変調しているため、波形の包絡線が変動している。図３９に、波形クリッパ１９０１が波形クリップした被測定信号の一例を示す。図３９においては、振幅成分が除去されているため、波形の包絡線が一定の値となっている。
【０１３９】
図４０（Ａ）は、図３７（Ａ）において説明したように、ジッタ測定装置１００がＡＤコンバータ１７０１を備える場合の、ＡＤコンバータ１７０１及び周期推定器１０１の動作の一例を示す。まず、ステップ１８０１においてＡＤコンバータ１７０１がアナログの被測定信号を離散化して、デジタルの被測定信号に変換する。次に、ステップ２０１で、周期推定器１０１が離散化された被測定信号に基づいて、被測定信号の基本周期を算出する。
【０１４０】
図４０（Ｂ）は、図３７（Ｂ）において説明したように、ジッタ測定装置１００が波形クリッパ１９０１を備える場合の、波形クリッパ１９０１及び周期推定器１０１の動作の一例を示す。まず、ステップ２００１において、波形クリッパ１９０１が被測定信号の振幅変調成分を除去する。次に、ステップ２０１において周期推定器１０１が被測定信号の基本周期を算出する。
【０１４１】
図４１は、周期ジッタ推定器１０５の構成の他の例を示す。本例における周期ジッタ推定器１０５は、図２９において説明した周期ジッタ推定器１０５の構成に加え、低域通過フィルタ２１０１を更に有する。低域通過フィルタ２１０１は、周期ジッタ系列補正器９０３が補正した周期ジッタ系列における高周波成分を除去する。低域通過フィルタ２１０１は、被測定信号の基本周波数程度の雑音成分を除去することにより、周期ジッタ系列の誤差を低減することができる。例えば、低域通過フィルタ２１０１は、被測定信号の基本周波数の１／２以上の周波数成分を除去する。本例における低域通過フィルタ２１０１と、図３１において説明した低域通過フィルタ２１０１は同様の機能を有してよい。
【０１４２】
図４２は、図２９において説明した周期ジッタ推定器１０５の動作の他の例を示すフローチャートである。本例において周期ジッタ推定器１０５は、図３０において説明した動作に加え、ステップ２２０１において、周期ジッタ系列の高調波成分を除去する。これにより、より精度よく周期ジッタ値を測定することができる。
【０１４３】
また、ジッタ検出器１０６は、ピーク・ツゥ・ピークジッタ検出器１０７、ＲＭＳジッタ検出器１０８、及びヒストグラム推定器１０９のいずれか１つ又は複数を備えてよい。
【０１４４】
また、ジッタ測定装置１００は、例えばコンピュータ等であってよい。この場合、コンピュータは、コンピュータをジッタ測定装置１００として機能させるプログラムを格納する。
【０１４５】
以上説明したように、ジッタ測定装置１００によれば、被測定信号の周期ジッタを精度よく測定することができる。特に、被測定信号を低い過剰標本化比で離散化した場合の、周期ジッタの測定精度を向上することができる。このため、従来のΔφ法に比べ、大幅に周期ジッタの測定精度を向上することができる。
【０１４６】
以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施形態に記載の範囲には限定されない。上記実施形態に、多様な変更または改良を加えることができる。そのような変更または改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲の記載から明らかである。
【０１４７】
上記説明から明らかなように、本発明のジッタ測定装置によれば、被測定信号の周期ジッタを精度よく測定することができる。
【産業上の利用可能性】
【０１４８】
上述説明から明らかなように、本発明に係るジッタ測定装置、及びジッタ測定方法によれば、被測定信号の周期ジッタを、高速且つ精度よく測定することができる。
【図１】
0004090989.tif 000036
【図２】
0004090989.tif 000037
【図３】
0004090989.tif 000038
【図４】
0004090989.tif 000039
【図５】
0004090989.tif 000040
【図６】
0004090989.tif 000041
【図７】
0004090989.tif 000042
【図８】
0004090989.tif 000043
【図９】
0004090989.tif 000044
【図１０】
0004090989.tif 000045
【図１１】
0004090989.tif 000046
【図１２】
0004090989.tif 000047
【図１３】
0004090989.tif 000048
【図１４】
0004090989.tif 000049
【図１５】
0004090989.tif 000050
【図１６】
0004090989.tif 000051
【図１７】
0004090989.tif 000052
【図１８】
0004090989.tif 000053
【図１９】
0004090989.tif 000054
【図２０】
0004090989.tif 000055
【図２１】
0004090989.tif 000056
【図２２】
0004090989.tif 000057
【図２３】
0004090989.tif 000058
【図２４】
0004090989.tif 000059
【図２５】
0004090989.tif 000060
【図２６】
0004090989.tif 000061
【図２７】
0004090989.tif 000062
【図２８】
0004090989.tif 000063
【図２９】
0004090989.tif 000064
【図３０】
0004090989.tif 000065
【図３１】
0004090989.tif 000066
【図３２】
0004090989.tif 000067
【図３３】
0004090989.tif 000068
【図３４】
0004090989.tif 000069
【図３５】
0004090989.tif 000070
【図３６】
0004090989.tif 000071
【図３７】
0004090989.tif 000072
【図３８】
0004090989.tif 000073
【図３９】
0004090989.tif 000074
【図４０】
0004090989.tif 000075
【図４１】
0004090989.tif 000076
【図４２】
0004090989.tif 000077