JP3971137B2 - Ｍｐｅｇのｐｃｒのジッタ、周波数オフセット及びドリフト・レートの測定方法 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、一般に、圧縮ビデオ信号の品質を測定する方法に関し、特に、モーション・ピクチャ・エンジニアリング・グループ（ＭＰＥＧ：Motion Picture Engineering Group又はMoving Picture Experts Group）のプログラム時刻基準参照値（ＰＣＲ：Program Clock Reference）のジッタ、周波数オフセット及びドリフト・レートを、ＰＣＲ到着時間が不均一でＰＣＲレートが可変である場合に、選択可能な一定測定帯域幅により測定する方法に関する。
【０００２】
【従来の技術】
デジタル・ビデオ放送（ＤＶＢ：Digital Video Broadcast）システムの如きデジタル伝送システムのデコーダ側にてプログラム・システム・タイミング・クロック（ＳＴＣ：System Timing Clock）を回復させることが、ビデオ信号を再生するために必要である。モーション・ピクチャ・エンジニアリング・グループ（ＭＰＥＧ：Motion Picture Engineering Group又はMoving Picture Experts Group）圧縮技術では、このＳＴＣを等時的に回復できる。このＭＰＥＧ圧縮技術において、クロック・カウントが２７ＭＨｚであるエンコーダにて、プログラムＳＴＣカウントの明示プログラム時刻基準参照値（ＰＣＲ：Program Clock Reference）サンプルは、転送ストリーム（ＴＳ：Transport Stream：誤りの発生する環境における伝送、蓄積に使用する）の中で伝送される。このＰＣＲ値は、（ｉ）原点にて正確であり、（ii）圧縮信号をデコード処理する際に問題が生じる点への伝送チェーンにおいて歪まないことが必要である。ＰＣＲ値の到着時間の間隔（インターバル）、期待値の確度、ＰＣＲ値により累積された到着時間及びそのジッタ、並びにＴＳにおける伝送到着時間の測定が、プログラム・ストリームのデコード可能性を確実にするために必要である。ジッタとして、周波数オフセット及び周波数ドリフト・レートが全体的な処理にとって重要なパラメータであり、ＰＣＲジッタを理解するために明確な定義が必要であり、その測定方法にガイダンスが求められている。
【０００３】
ＩＳＯ／ＩＥＣ１３８１８−１に記載された規格から、公称ＳＴＣ周波数の２７ＭＨｚからの周波数偏差用の限界マスクを図２に示すように定義することが可能である。この図２は、周波数の対数に対する時間の対数で表す周波数ドリフト限界、周波数オフセット限界及び位相エラー限界をマスクにより示している。周波数オフセットは、エンコーダにおけるＳＴＣの公称周波数と実際の値との差である。２７ＭＨｚでは、限界が＋／−８１０Ｈｚに設定されており、即ち、＋／−３０ｐｐｍとして正規化されている。変動の周波数レート又はドリフト・レートは、プログラムＳＴＣの周波数が時間と共に変化する「速度」である。すなわち、これは、時間に対する周波数の第１導関数であるか、又は時間に対する位相の第２導関数である。この限界を２７ＭＨｚのプログラムＳＴＣに対して７５ｍＨｚ／秒に設定するか、又は２７ＭＨｚのＳＴＣ周波数に１０ｐｐｍ／時間として正規化する。この位相許容度を＋／−５００ナノ秒に設定するが、これは、ＴＳにおける時間位置に応じた値に対するＰＣＲ値の最大エラーを表す。
【０００４】
周波数オフセット及びドリフト・レートの限界をＳＴＣに対して設けている。これは、ＳＴＣを対応するＰＣＲフィールドの値で表せるためである。これら限界には、プログラムＳＴＣの影響や、ＰＣＲ計算における任意の可能性のあるエラーを含む。５００ナノ秒の限界を、プログラムＳＴＣには設けるが、ＴＳにおけるＰＣＲ位置に対するＰＣＲ値を表す確度には設けない。しかし、ＳＴＣを再現するためにデコード点にてＰＣＲ値を用いる際、ＰＣＲエラーは、位相及びジッタのエラーと完全に等価である。
【０００５】
プログラムＳＴＣ周波数エラー、ドリフト・レート及びＰＣＲジッタを試験する方法が、１９９６年に発行されたＩＳＯ／ＩＥＣ１３８１８−９に記載されている。この試験は、１度にＴＳにて１つのプログラムに対して実施される。この手順は、ＰＣＲ値の図を用いる。この図は、図９に示すように、Ｎ個のサンプルに適合した線形回帰線をＰＣＲ到着時間に対するＰＣＲ値（時間に変換されている）で表す。すなわち、ＰＣＲ到着時間に対して、データが伝送された際のプログラムＳＴＣ時間を表す。到達時間は、Ｘ軸上にプロットされるが、受信器での正確な時間を表す。ＰＣＲ値は、Ｙ軸上にプロットされる。この線形回帰直線は、Ｎサンプルに適合し、この直線の傾きマイナス１が周波数オフセットである。このジッタは、この直線の上／下のＰＣＲ値の距離として定義される。
【０００６】
この方法は、充分なサンプルを累積するのに大きなＮがどの程度かを決めなければならないという問題がある。Ｎが大きすぎると、低周波数タイミングの不安定さ、即ち、ワンダにより、２次方程式（quadratic）への適合（fit）よりも直線カーブへの適合を用いる際に、ジッタの測定値が過度になる。線形適合に対して、この直線の傾斜を用いて、周波数オフセットを予測する。しかし、周波数オフセットは、しばしばドリフトしたり、ワンダを起こすので、これらワンダの影響によって、ジッタ値が不当に大きくなるかもしれない。これらジッタ値は、ジッタ規格をはるかに超えるが、デコーダの位相拘束ループ（ＰＬＬ）により追跡するので、ドリフト・レートが限界値未満ならば、これらジッタの問題がない。したがって、これら値は、ジッタ測定に存在せず、ドリフト・レートに対して別に測定しなければならない。この逆もまたドリフト・レートの測定では正しく、ジッタ規格内で適当な小さなタイミング変動は、ドリフト・レート規格をしばしば遙かに超える。これは、タイミング又は位相エラーの第２導関数に比例するためである。したがって、ドリフト・レートの測定は、許容可能なジッタ限界内でピーク変動を伴うタイミング・エラーを含まないように、ロウパス・フィルタ処理しなければならない。２次方程式曲線（quadratic curve）への適合を行って、周波数ドリフト・レートに対するコンプライアンスを確認するのが好ましいが、Ｎが小さすぎる場合、ジッタ限界内のタイミング変動により、ドリフト・レートの推定値が常に規格を超えてしまう。ＰＣＲが充分な規則正しさで到着しないならば、再生したＳＴＣには、ジッタやドリフトがあるかもしれない。また、受信器／デコーダは、拘束から外れるかもしれない。
【０００７】
ＴＳインタフェース点で見えるタイミングの影響は、図１０に示すようにモデル化できる。図１０は、ＭＰＥＧ伝送システムにおけるプログラム・ストリーム用のＰＣＲ値における変動のソース（源）を示すモデルのブロック図である。これは、エンコーダ／リマルチプレクサ及びコミュニケーション・ネットワークの理想化したモデルである。エンコーダにおけるモデルの要素は、特定のプログラム・ストリームに対するＰＣＲに特定されるが、ネットワークにおける要素は、全体のＴＳに関連する。
【０００８】
このモデルは、２７ＭＨｚの如き公称周波数Ｆのシステム・クロック周波数発振器を有するが、この実際の周波数は、公称周波数から関数ｆ_dev（ｐ、ｔ）だけ離れている。この関数は、時間（ｔ）に応じて決まり、単一のプログラム・ストリーム（ｐ）に特定される。これは、公称周波数に対する周波数オフセットを表し、ＭＰＥＧ標準に応じて、特定の許容度の範囲に留まる。ドリフト・レートは、ｆ_dev（ｐ，ｔ）の時間に対する変化のレートである。このドリフト・レートは、ＭＰＥＧ標準が規定する限界未満でなければならない。
【０００９】
システム・クロック周波数発振器（基準クロック発振器）は、理想的なＰＣＲカウントＮ_p,iを発生するＰＣＲカウンタを駆動する。なお、「ｐ」は、特定のプログラム・ストリームであり、「ｉ」は、ＴＳ内のビット位置である。これにより、ＰＣＲの不正確なソースＭ_p,iからの値を付加して、ＴＳ内で見えるＰＣＲ値Ｐ_p,iを生成する。これら値を加算した結果の簡単な値は、次のようになる。
Ｐ_p,i＝Ｎ_p,i＋Ｍ_p,i
【００１０】
エンコーダ（図１０の点線より左側）の出力端において、ＴＳは、一定ビット・レートＲ_nomで伝送される。エンコーダの先のコミュニケーション・ネットワーク（図１０の点線より右側）は、可変遅延Ｊ（ｉ）が導入されることにより、ＭＰＥＧパケットの公称到着時間から、実際の到着時間が変動する。なお、Ｊ（ｉ）は、ＴＳにおけるｉ番目のビットに対する秒単位の遅延である。この公称到着時間は、Ｒ_nomと、ＴＳにおけるビット・インデックスから計算する。実際の到着時間Ａ（ｉ）の式は、次のようになる。
Ａ（ｉ）＝ｉ／Ｒ_nom ＋ Ｄ ＋ Ｊ（ｉ）
Ｄは、時間原点及びビット・インデックス原点が任意であるという事実と、コミュニケーション・ネットワークに平均遅延があるという事実を表す定数である。Ｊ（ｉ）は、ネットワーク遅延におけるジッタを表し、総ての時間にわたる平均値は、ゼロに定義される。すなわち、これは、Ｄの値を設定する。Ｊ（ｉ）＋Ｍ_p,jは、全体的なＰＣＲジッタであり、ＭＰＥＧ標準にて設定された最大エラー規格である。
【００１１】
総ての測定システムが、定まった帯域幅で、同じプログラムＰＣＲでは同じ値を測定するように、周波数オフセット、ドリフト・レート及びジッタを判定する一定帯域幅方法が好ましい。また、一定帯域幅測定は、ＭＰＥＧデコーダで回復したＳＴＣがいかにＰＣＲエラーに応答するかによって、より物理的な意味を持つ。これは、デコーダが、典型的には、一定帯域幅を有し、ＰＣＲ到着インターバル又はＰＣＲ平均レートに独立して動作するためである。
【００１２】
【発明が解決しようとする課題】
不均一なＰＣＲ到着時間及び可変ＰＣＲレートにわたって選択可能で一定の測定帯域幅により、ＭＰＥＧにおけるＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定する方法が望まれている。
【００１３】
したがって、本発明は、ＰＣＲ到着時間が不均一で、ＰＣＲレートが可変の環境にて、選択可能な一定の測定帯域幅により、ＭＰＥＧにおけるＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定できる方法を提供するものである。
【００１４】
【課題を解決するための手段】
本発明は、プログラム時刻基準参照値（ＰＣＲ）の到着時間が不均一であり、ＰＣＲレートが可変である場合に、選択可能な一定の測定帯域幅で、ＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定する方法であって；非同期精密基準発振器によりＰＣＲ到着時間の間の間隔（インターバル）を測定し；可変のＰＣＲレートから独立して、選択可能な一定測定帯域幅を、ジッタ及びワンダの間の周波数を区別する一定測定帯域幅に固定し；１組の差分方程式の形式で、ＰＣＲ到着時間の間隔、各ＰＣＲ到着時間でのＰＣＲ値、固定された選択可能な一定測定帯域幅に応じて、ＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定するものである。
【００１５】
本発明においては、選択された定帯域幅は、ジッタをワンダから区別するためであり、ＭＰＥＧ標準におけるジッタ、周波数、ドリフト・レートに対して特定した標準に基づいて限界マスクを用いて得られる。ＰＣＲ到着時間は、基準クロックにより測定され、測定ＰＣＲインターバルに変換される。測定したＰＣＲインターバル、ＰＣＲ値及び選択された定帯域幅は、一組の差分方程式に導入されて、ジッタ、周波数及びドリフト・レートの値を求める。この一組の差分方程式は、形式IIの位相拘束ループからか、又は平均最小二乗（ＬＭＳ）線形及び２次方程式推定量から求める。
【００１６】
本発明の目的、利点及び新規な特徴は、添付図を参照した以下の詳細な説明から理解できよう。
【００１７】
【発明の実施の形態】
ＭＰＥＧシステム・タイミング・クロック（ＳＴＣ）の周波数オフセット、周波数ドリフト・レート及びジッタを、伝送されたＰＣＲによって測定する２つの方法を説明する。各方法は、コンピュータによる実施に適する１組の簡単な差分方程式として記述される。これら方法は、ＰＣＲレートから独立した不均一なサンプル・データにて動作する定帯域幅周波数分析器（constant bandwidth filter）に作用する。本発明により周波数オフセット、ディザ及びジッタを測定する測定フィルタの簡略化した全体的なブロック図である図１に示すように、各方法は、２つの別個の時間入力を有する。すなわち、これら時間入力とは、（１）ＭＰＥＧ伝送ストリーム（ＴＳ）からの所定プログラム・ストリーム用のデコードされたＰＣＲ値であるＰＣＲ（ｎ）と、（２）ＰＣＲ到着時間から非同期精度基準発生器によって測定する如き所定プログラム・ストリームに対する複数のＰＣＲ値の間のＰＣＲインターバルであるＰＣＲ＿ｉｎｔ（ｎ）とである。各方法は、単一の制御パラメータである測定帯域幅（ＢＷ）を有する。この測定帯域幅ＢＷは、ジッタ及びワンダの間の境界周波数に関連する。すなわち、この測定帯域幅は、選択可能であり、平均ＰＣＲレートから独立した定数である。各方法は、３つの個別の時間出力測定結果を出力する。これら測定結果は、（１）ＳＴＣ周波数オフセットであるＦＯ（ｎ）と、（２）ＳＴＣ周波数ドリフト・レートであるＤＲ（ｎ）と、（３）ＰＣＲジッタであるＪｉｔ（ｎ）とである。
【００１８】
第１方法用の方程式は、図３に示すハイブリッド・デジタル・アナログ位相拘束ループ（ＰＬＬ）モデルから得られる。ここでは、結果としての帰還制御システムにおける種々の点から、上述の３つの出力パラメータＦＯ（ｎ）、ＤＲ（ｎ）、Ｊｉｔ（ｎ）を抽出している。ハイパス・フィルタ（ＨＰＦ）及び縮尺（拡大／縮小）（１／Ｆｏとの乗算器）をＪｉｔ（ｎ）の信号路に追加して、ワンダ成分の分離を改良すると共に、１秒間当たりの時間単位に出力を縮尺する。また、のロウパス・フィルタ（ＬＰＦ）をＦＯ（ｎ）及びＤＲ（ｎ）の信号路に追加して、ジッタ成分の除去を改善する。実際的な差分方程式をＰＬＬモデルから求めて、不均一なＰＣＲの到着時間から、選択可能な一定帯域幅により、周波数オフセット、ドリフト・レート及びジッタを個別に測定できる。
【００１９】
図３に示す形式IIのＰＬＬは、２つの理想的な積分器、即ち、デジタル・アキュムレータ（加算器＋及びラッチで構成）と電圧制御発振器（ＶＣＸＯ）とを具えている。これにより、位相引き算の出力における２次ハイパス閉ループ応答を生成する。したがって、ループ帯域幅の下において、その応答がドリフト・レートに比例し、ループ帯域幅の上において、その応答がジッタに比例する。１次ハイパス・フィルタをジッタ測定に追加して、ドリフト・レートの影響を除く。ロウパス・フィルタをドリフト・レート出力に追加して、この測定からジッタの影響を除く。生のデータにデジタル信号処理（ＤＳＰ）を使用して、これらフィルタを最良に実現する。また、ＰＣＲレートはサンプル・レートなので、平均ＰＣＲレートを求め、係数を調整して、正確な測定フィルタ帯域幅及び制御ループ・ダンピングを近似的に維持する。自走発振器及びデコードしたＰＣＲ値を用いて、ＰＬＬを、測定したＰＣＲインターバルで駆動される完全サンプル・システムに変換することが可能であるので、かかる帯域幅は、不均一ＰＣＲ値から独立であり、明示電圧制御発振器を用いずにソフトウェアにより測定値を計算する。なお、ＤＡＣは、デジタル・アナログ変換器であり、なお、ＰＬＬは、ＤＡＣ、ループ・フィルタ、電圧制御発振器（ＶＣＸＯ）、このＶＣＸＯの出力をクロックしてロード信号ＬＤでＰＣＲ（Ｏ）をロードするカウンタ、ラッチ回路、加算器＋及び積分器で構成される。
【００２０】
０．１ｐｐｍの確度の自走水晶発振器（ＯＣＸＯ）である電圧制御発振器（ＶＣＸＯ）と、＋／−３０ｐｐｍレンジでこのＶＣＸＯを拘束する位相拘束ループ（ＰＬＬ）とを有するサブシステムから、ＶＣＸＯ回路を構成する。このＶＣＸＯ回路を用いて、望ましくは、０．１ｐｐｍの確度で周波数オフセット出力を校正できる。そうでない場合には、ＶＣＸＯ回路を単独で用いてもよく、既知の高確度の周波数を供給し、その後の測定結果からエラーを減算して、周波数エラー又はオフセットを確認できる。これは、ＰＬＬの実施（インプリメンテーション）をハイブリッドＤＳＰ及びアナログ信号処理方法として記述できる。
【００２１】
閉ループ応答を複製したＰＬＬモデルから差分方程式の実際の設定を導出できるので、不均一なＰＣＲ到着時間から独立して選択可能な固定（一定）帯域幅により、周波数オフセット、ドリフト・レート及びジッタを測定できる。かかる方程式を求める等価回路を図４及び図５に示す。ＰＣＲ（ｎ）値は、微分器（図４）の入力であり、微分器からの出力は、不均一サンプル・インターバル用の２次ロウパス・フィルタ（ＬＰＦ）に入力する。また、２次ロウパス・フィルタの入力は、最後に受信したＰＣＲ（ｎ−１）からの対応ＰＣＲ＿ｉｎｔ（ｎ）である。
【００２２】
フィルタ出力ｙ（ｎ）を全体的な開ループ直流利得であるＧｏｄｃと乗算して、正規化された周波数を得る。この正規化周波数から１（unity）を減算して、バイポーラ周波数オフセットＦＯを求める。２つの連続的なサンプルｙ（ｎ）及びｙ（ｎ−１）であるフィルタ出力を第２微分器（図５）に入力して、ジッタ・エラー出力Ｊ（ｎ）を発生する。このジッタ・エラー出力Ｊ（ｎ）は、１次ロウパス・フィルタ（ＬＰＦ）に入力されてドリフト・レートＤＲ（ｎ）を生じると共に、１次ハイパス・フィルタ（ＨＰＦ）に入力されてジッタＪｉｔ（ｎ）を生じる。この図示の回路から、次の差分方程式が得られる。
y(n)=2*y(n-1)-y(n-2)-(y(n-1)-N1*y(n-2))*5.932*BW*PCR_int(n)+(PCR(n)-PCR(n-1))/FOなお、＊は、乗算を意味する。ＦＯは、２７ＭＨｚであり、ＢＷは、ヘルツである測定帯域幅であり、Ｎ１は、1/(1+pi*BW/60)に等しい定数であり、piはπである。ｙ（ｎ）の周波数オフセットＦＯ（ｎ）は、ｙ（ｎ）＊Ｇｏｄｃ−１であり、Ｇｏｄｃは、(1-N1)*5.932*BWに等しい定数である。ジッタJit(n)は、Jit(n)=A1*(Jit(n-1)+J(n)-J(n-1))となる。なお、J(n)=y(n)-y(n-1)であり、A1は、exp(-2*pi*BW*0.68/30)となる。最後に、ドリフト・レートＤＲ（ｎ）は、A2*DR(n-1)+(Godc/Tp)*J(n)である。なお、Tpは、1/(2*pi*.01) であり、A2は、exp(-1/30*Tp)である。この方法で、ＤＲ（ｎ）は、０．０１Ｈｚの固定帯域幅によって後でフィルタ処理され、ｙ（ｎ）は、１／Ｇｏｄｃに初期化されて、設定時間を最小にする。
【００２３】
第２方法は、線形フィルタ近似から、ＰＣＲ値用の一元２次時間方程式の平均最小二乗（ＬＭＳ）周波数オフセット推定量（線形項係数）及びＬＭＳドリフト・レート推定量（２次項係数）を導出する。この結果の差動方程式は、上述の方法と類似の結果を生じるが、計算がわずかに簡単である。
【００２４】
正確な基準クロックを用いて、受信したＰＣＲ（ｎ）の各々の測定到着時間としてＰＣＲ＿ＡＴ（ｎ）を得るために、次式が、プログラム・デコーダＳＴＣ周波数オフセット、ドリフト・レート及びジッタ用の項を有する期待ＰＣＲ（ｎ）値をＦｓｔｃのサイクルで定義する。
PCR(n)=a+PCR_AT(n)*Fstc+(PCR_AT(n))^2*c+∈(PCR_AT(n))
なお、^は、べき乗を意味する。この最終項は、時間の関数としてのジッタである。ＰＣＲ＿ＡＴ（ｎ）は、到着時間ジッタ及びＰＣＲ（ｎ）エンコード・エラーの組合せの影響によりＰＣＲ（ｎ）に追加されるランダム変数である。これは、真に現実的であり、ゼロ平均で、独立し、同じに分布されたランダム変数と仮定する。したがって、ＰＣＲ（ｎ）値は、それ自体がI.I.D.ランダム変数である。ａは、基準発振器に関連した初期位相エラーであり、Fstcは、システム・タイミング・クロック周波数であり、ｃは、ドリフト・レートの半分である。これは、ドリフト・レート（Ｈｚ／ｓ）が、位相の第２時間導関数ＤＲ＝２＊ｃである。
【００２５】
既知の理想的なシステム・クロック周波数であるＦｏ＝２７ＭＨｚを減算して、ＳＴＣ周波数オフセットを見つけ、Ｆｏで除算して、これらの項を正規化する。ＳＴＣドリフト・レートは、許容できる周波数オフセット値に対して典型的には非常に小さく、周波数オフセット測定内に含ませることが好ましいので、この方程式は、ドリフト・レート項ｃを用いずに、周波数オフセットｂの分析用に簡略化できる。ドリフト・レート及び高次の導関数は、高次の要素を除去した周波数オフセットの定義よりも、周波数オフセット値に影響する。ＰＣＲエラー値は、次のようになる。
PCR_err(n)=a+PCR_AT(n)*b+∈(PCR_AT(n))秒
ここで、b=(Fstc-Fo)/Foである。ジッタ項は、秒単位のジッタであり、ａは、秒単位での初期位相オフセットであり、ｂは、次元のない比としてのＳＴＣの周波数オフセットである。ドリフト・レートを無視したこの式は、ジッタによる平均平方エラーを最小にする係数ａ及びｂの推定値から決まるデータの線形回帰又は直線適合を可能にするＰＣＲ（ｎ）用の式である。これは、二重の結果を提供し、先ず、ｂ用のＬＭＳ推定値としてのＳＴＣ周波数オフセットＦＯを決定し、次に、ＭＰＥＧ標準が提案するようにｂ及び初期位相値ａによるＦＯ影響を引くことによりジッタを決定する。
【００２６】
直線適合を計算するために用いる多くのＰＣＲ（ｎ）値を決定するために、線形フィルタ処理に対する線形回帰分析の比較が必要になる。この方法は、ＰＣＲ（ｎ）データへの線形適合の平均化の実効帯域幅であり、平均ＰＣＲレート用に見つけることができる。なお、ＰＣＲ（ｎ）値は、均一平均レートで到着すると一時的に仮定する。均一なレート仮定は、非常に多くのサンプルによる前後関係にて正当化される。また、これら非常に多くのサンプルを用いて、ＭＰＥＧ標準の非常に少ないドリフト・レート許容度である１０ｐｐｍ／時間により負わされた制約の中で適切にＦＯを見積る。直線適合回帰分析に対するフィルタ応答の良好な予測が求まると、一定帯域幅フィルタ処理を不均一なサンプル・データ用に構成する。これら不均一なサンプル・データ自体が、ＰＣＲレートから独立した一定帯域幅測定を発展させるための瞬間サンプル・インターバルの関数である。
【００２７】
これは、重要であるが、ＭＰＥＧ標準の公表された限界と比較するために、ＰＣＲジッタ、ＳＴＣ周波数オフセット及びドリフト・レートの実際的で一定の見積り値を見つけるように機械化された連続的フィルタ処理が可能なためである。さらに、測定の帯域幅が一定のため、ＳＴＣドリフト・レート及びＰＣＲ精密規格に基づく標準化境界周波数を用いて、ＭＰＥＧデコーダ内の回復ＳＴＣクロック特性に直接関係した実際的で標準化されたフィルタ帯域幅を求めることができる。
【００２８】
線形回帰分析：
τＯを基準クロックにより決定した平均ＰＣＲインターバルとし、均一な到着時間が個別にτＯ秒だけ離れていると仮定する。
PCR_err(n)=a+b*τo*n+∈(n*τo) 秒
ＰＣＲ値のＮ個のサンプルを得るために、ａ及びｂの平均最小二乗（ＬＭＳ）推定値がベクトルＢで与えられる。
B=[a;b]; B= T^-1*(N*N^T)^-1*(N^T*X)
なお、
T^-1=[1,0; 0,τo^-1]; (N*N^T)^-1= [2(2N+1),-6; -6,12/(N+1)]*(1/(N(N-1)))
(N^T*X):= [Σ_{n=1 〜 N}PCR_err(n);Σ_{n=1 〜 N}PCR_err(n)*n].
である。初期位相は、次のように求まる。
a=2/(N*(N-1)*τo){Σ_{n=1 〜 N}PCR_err(n)*(2N+1-3n)}
Ｎ個の均一に離れたサンプル用の周波数オフセットを次のように見積る。
b= 6/(N*(N-1)*τo){Σ_{n=1 〜 N}PCR_err(n)*(2n/(N+1)-1)}.
【００２９】
２次方程式回帰分析：
ＳＴＣ周波数ドリフト・レートを無視することなく、次のように、ＬＭＳ２次方程式曲線適合を見つけ、係数ｃを見積る。なお、ｃ＝（ＳＴＣ周波数ドリフト・レート）／２＝ＤＲ／２である。
PCR_err(n)=a+b*τo*n+c*τo²*n²+∈(n*τo) 秒
B=[a; b; c]=T^-1(N*N^T)^-1(N^T*X); T^-1=[1,0,0; 0,τo^-1,0; 0,0,τo^-2];
(N*N^T)^-1=[3(3(N+1)+2),-18(2N+1),30;
-18(2N+1),{12(2N+1)(8N+11)}/{(N+1)(N+2)},-180/(N+2);
30,-180/(N+2), 180/{(N+1)(N+2)}]/{N(N-1)(N-2)};
(N^T*X)=[Σ_{n=1 〜 N}PCR_err(n); Σ_{n=1 〜 N}PCR_err(n)*n;Σ_{n=1 〜 N}PCR_err(n)*n²]
次に、Ｎ個の均一に離間したサンプル用の周波数ドリフト・レートに対するＬＭＳ見積りは、次のようになる。
DR=60/{N(N-1)(N-2)τo²}*Σ_{n=1 〜 N}PCR_err(n){6n²/{(N+1)(N+2)}-6n/(N+2)+1}.
【００３０】
これら２つのＬＭＳ推定値、線形回帰からの周波数オフセット、及び２次方程式回帰からの周波数ドリフト・レートを、有限インパルス応答（ＦＩＲ）フィルタとして実現できる。まず、ｚ及び平均ＰＣＲ値の数Ｎに応じたＦＩＲフィルタによって周波数オフセットを見積るＬＭＳを記述する。
Hr(z,N)=6/{N(N-1)τo}*z^-N*Σ_{k=0 〜 (N-1)}z^k{2(k+1)/(N+1)-1}
ＦＯを見積るのに用いたサンプルの数は、N=Pr/(2*BW)である。なお、Ｐｒは、平均ＰＣＲレートであり、τｏは、１／Ｆｓであり、ＢＷは、測定帯域幅である。したがって、例えば、Pr=24及びBW=0.1Hzで、N=floor(Pr/(2*BW)) = 125サンプルである。第２に、ＦＩＲフィルタをｚの関数として記述すると、周波数ドリフト・レートの見積りは次のようになる。
Hdr(z,N)={60/(N(N-1)(N-2)τo²)}z^-NΣ_{k=0 〜 N-1}{(6(k+1)²)/((N+1)(N+2))-(6(k+1))/(N+2)+1}
【００３１】
比較において、帯域幅がＢＷの２次循環バターワース・ロウパス・フィルタが微分器とカスケード接続される。循環ＬＰＦ及び微分器の応答のカスケード接続を、図７のＦＩＲ ＬＭＳ周波数オフセット推定値と比較する。ドリフト・レートＦＩＲ応答は、図８に示すようになる。なお、
N=floor(Pr/(2^1/2*Bw)=サンプル数
が、Ｐｒの平均ＰＣＲレートにてＦＩＲフィルタで処理される。比較のため、２個の微分器とカスケード接続された帯域幅ＢＷの３次循環バターワースＬＰＦを、図８のＦＩＲ ＬＭＳドリフト・レート推定値と比較する。
【００３２】
この分析より、簡単な微分器とカスケード接続された２個の線形循環フィルタ処理に対する均一なサンプル・データの線形及び２次方程式回帰の精密な近似が存在する。Ｎ個のサンプルの線形回帰分析を置換する等価循環フィルタ処理の平均ＰＣＲレートＰｒと、帯域幅ＢＷとの推定値を得ることは、ＢＷ＝Ｐｒ／（２＊Ｎ）であり、２次方程式回帰では、ＢＷ＝Ｐｒ／（２^1/2＊Ｎ）である。これら等価処理は、上述の如き微分器及びバターワースＬＰＦフィルタ用の差分方程式から構成してもよい。等価回帰フィルタは、２つの前のＰＣＲ値及び相対的に小さなかけ算／足し算を蓄積するだけでよいので、これらは、計算が簡単であり、一定帯域幅測定結果を生じる単一のプロセッサにより多くのプログラムＰＣＲに対して同時に実施できる。フィルタ処理からの推定（見積り）値を用い、ＨＰＦを適用して結果からドリフト・レート成分を除去することにより、周波数オフセットによる推定ＰＣＲ値を減算して、ＰＣＲジッタを求めてもよい。
【００３３】
この第２方法の等価回路を図６に示す。ＰＣＲ＿ｉｎｔ（ｎ）をＰＣＲ＿ＡＴ（ｎ）時間から求め、ＰＣＲ（ｎ）値の微分から減算し、正規化し、ＰＣＲ（ｎ）エラー値Ｆｏｏ（ｎ）を発生する。Ｆｏｏ（ｎ）は、第１線形ＬＭＳ推定値への入力であり、この推定値を求めるには、導関数と、２次ロウパス・フィルタ（ＬＰＦ）とを含んでおり、周波数オフセットＦＯ（ｎ）を発生する。周波数オフセットは、第２の２次方程式ＬＭＳ推定値に入力し、この推定値を求めるには、導関数、微分回路、１次ロウパス・フィルタを含んでおり、ドリフト・レートＤＲ（ｎ）を発生する。周波数オフセットをＰＣＲ＿ｉｎｔ（ｎ）により乗算して、その結果をＰＣＲ（ｎ）エラー値から減算して、瞬間ジッタ値ＪＪ（ｎ）を生成する。この瞬間ジッタ値をカスケード・フィルタで処理する。このカスケード・フィルタは、まず、１次ハイパス・フィルタが後続する１次ロウパス・フィルタを具え、ＰＣＲジッタＪｉｔ（ｎ）を発生する。その結果の差分方程式は、次のようになる。
FO(n)=PCR_int(n)*(π*BW)²*Foo(n)*E(n)+2E(n)*FO(n-1)-E(n)²*FO(n-2)
なお、Foo(n)=[PCR(n)-PCR(n-1)]/Fo-PCR_int(n);
E(n)=exp[-PCR_int(n)*π*BW] (exp(-x)〜1-x+x²/2-x³/3);
DR(n)=E(n)*DR(n-1)+π*BW*[FO(n)-FO(n-1)];
JJ(n)=Foo(n)-FO(n)*PCR_int(n);
Jit(n)=E(n)*[JJ(n)-JJ(n-1)+2*Jit(n-1)-E(n)*Jit(n-2)].
である。
【００３４】
よって、本発明の方法は、不均一なインターバルで到着するＰＣＲ値と、選択可能な一定帯域幅とからＭＰＥＧにおけるＰＣＲのジッタと、周波数オフセットと、ドリフト・レートとを求める。この帯域幅は、形式II位相拘束ループから、又はＬＭＳの１次及び２次推定値から得た１組の差分方程式を用いて、ジッタからワンダを分離する周波数である。
【００３５】
【発明の効果】
上述の如く、本発明によれば、ＰＣＲ到着時間が不均一で、ＰＣＲレートが可変である場合に、選択可能な一定測定帯域幅を利用して、ＭＰＥＧにおけるＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定できる。
【図面の簡単な説明】
【図１】本発明により周波数オフセット、ディザ及びジッタを測定する測定フィルタの簡略化した全体的なブロック図である。
【図２】本発明による周波数オフセット・ドリフト及びジッタ・マスクのグラフを示す図である。
【図３】本発明により周波数オフセット、ディザ及びジッタを測定するための位相拘束ループ（ＰＬＬ）のブロック図である。
【図４】本発明により、１組の差分方程式を得るための図３の等価回路のブロック図の一部分を示す図である。
【図５】本発明により、１組の差分方程式を得るための図３の等価回路のブロック図の他の部分を示す図である。
【図６】本発明により、１組の差分方程式を得るために平均最小二乗１次及び２次推定量を用いる等価回路のブロック図である。
【図７】本発明により、ＦＩＲ線形フィルタ及び２次回帰フィルタと、微分器及びロウパス・フィルタの組合せ及び二重微分器及びロウパス・フィルタの組合せとの比較を示す図である。
【図８】本発明により、ＦＩＲ線形フィルタ及び２次回帰フィルタと、微分器及びロウパス・フィルタの組合せ及び二重微分器及びロウパス・フィルタの組合せとの比較を示す図である。
【図９】従来技術により周波数オフセット及びジッタを見積るためのグラフを示す図である。
【図１０】ＭＰＥＧ伝送システムにおけるプログラム・ストリーム用のＰＣＲ値における変動の源を示すモデルのブロック図である。

Claims (2)

1. プログラム時刻基準参照値（ＰＣＲ）の到着時間が不均一であり、ＰＣＲレートが可変である場合に、選択可能な一定の測定帯域幅で、ＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定する方法であって、
   非同期精密基準発振器によりＰＣＲ到着時間の間の間隔を測定し、
   上記可変のＰＣＲレートから独立して、上記選択可能な一定測定帯域幅を、ジッタ及びワンダの間の周波数を区別する一定測定帯域幅に固定し、
   １組の差分方程式の形式で、上記ＰＣＲ到着時間の間隔、各ＰＣＲ到着時間でのＰＣＲ値、上記固定された選択可能な一定測定帯域幅に応じて、上記ＰＣＲのジッタ、周波数オフセット及びドリフト・レートを測定する測定方法。
2. 上記１組の差分方程式は、ＰＣＲ値用の一元２次時間方程式の平均最小二乗（ＬＭＳ）周波数オフセット及びＬＭＳドリフト・レート推定量に対する線形ろ波近似から導出したことを特徴とする請求項１の測定方法。

Images