[論文レビュー] On the measurement of frequency and of its sample variance with high-resolution counters

Enrico Rubiola|HAL (Le Centre pour la Communication Scientifique Directe)|Nov 25, 2004

Advanced Electrical Measurement Techniques被引用数 4

ひとこと要約

本稿は、高分解能カウンタを用いた周波数測定およびその標本分散の分析を通じて、異なるカウンタ型（Π型およびΛ型）がアラン分散（AVAR）および修正アラン分散（MVAR）推定に与える影響を明らかにする。重ね合わせ平均化を用いるΛ型カウンタでは、短期安定性において自然にMVARが得られ、長期平均化ではAVARに収束することが示され、周波数安定性解析における一般的な誤解を解消する。

ABSTRACT

A frequency counter measures the input frequency $\barν$ averaged over a suitable time $τ$, versus the reference clock. High resolution is achieved by interpolating the clock signal. Further increased resolution is obtained by averaging multiple frequency measurements highly overlapped. In the presence of additive white noise or white phase noise, the square uncertainty improves from $\smash{σ^2_ν\propto1/τ^2}$ to $\smash{σ^2_ν\propto1/τ^3}$. Surprisingly, when a file of contiguous data is fed into the formula of the two-sample (Allan) variance $\smash{σ^2_y(τ)=\mathbb{E}\{\frac12(\bar{y}_{k+1}-\bar{y}_k) ^2\}}$ of the fractional frequency fluctuation $y$, the result is the \emph{modified} Allan variance mod $σ^2_y(τ)$. But if a sufficient number of contiguous measures are averaged in order to get a longer $τ$ and the data are fed into the same formula, the results is the (non-modified) Allan variance. Of course interpretation mistakes are around the corner if the counter internal process is not well understood.

研究の動機と目的

• 高分解能周波数カウンタを用いた際のアラン分散（AVAR）と修正アラン分散（MVAR）の解釈に関する広範な誤解を解消すること。
• 重ね合わせ平均化およびカウンタ構造（Π型対Λ型）が有効分散推定に与える影響を明確にすること。
• 短期および長期安定性解析の両方において、AVARおよびMVARの観点から周波数カウンタ出力の正しい解釈を可能にする厳密なフレームワークを提供すること。
• 平均化および窓関数の背後にあるメカニズムを理解せずに標準分散式にカウンタデータを入力する際の重大な誤解の問題に取り組むこと。

提案手法

• 周波数および位相の揺らぎを数学的にモデル化し、分数周波数揺らぎ $ y(t) = \frac{\dot{\phi}(t)}{2\pi\nu_{00}} $ を定義することで、安定性指標を分析する。
• 周波数カウンタのインパulse応答を分析し、Π型カウンタが長方形窓平均化を用いるのに対し、Λ型カウンタが三角形（積分）窓を用いることを示す。
• 両カウンタ型の有効重み関数を導出する：Π型のための $ w_{\Pi}(t) $ およびΛ型のための $ w_{\Lambda}(t) $ であり、これらはAVARまたはMVARが実際に測定されるかを決定する。
• Λ型カウンタにおける重ね合わせ測定が、三角形形状の有効窓関数を生じさせ、AVAR式に投入するとMVARが得られることを示す。
• Λ型測定値 $ m $ 個の連続平均化が、台形窓を生成し、$ m $ が大きい場合には長方形窓に漸近的に収束することを示し、その結果、AVARに収束することが確認される。
• 分散推定器をウェーブレットに類似した解釈により分析し、AVARが中心周波数 $ 1/(2\tau) $ の半オクターブバンドパスフィルタとして機能することを示す。

実験結果

リサーチクエスチョン

• RQ1 Π型およびΛ型周波数カウンタは、有効平均化窓およびそれに伴う分散推定においてどのように異なるか？
• RQ2 なぜΛ型カウンタのデータをAVAR式に投入するとMVARが得られ、AVARではないのか？
• RQ3 Λ型カウンタにおける重ね合わせ平均化がMVARではなくAVARに収束するのはどのような条件下か？
• RQ4 窓形状（長方形対三角形）が周波数安定性分散推定に与える影響は何か？
• RQ5 なぜΛ型カウンタのデータをAVARに変換するには測定時間 $ \tau $ を延長しなければならないのか？

主な発見

• 重ね合わせ測定を伴うΛ型カウンタは、三角形形状の有効平均化窓を生成し、標準AVAR式に投入すると、修正アラン分散（MVAR）が得られる。
• 長期安定性解析（$ m \gg 1 $）において、Λ型カウンタの平均化プロセスは長方形窓に収束するため、AVAR式を用いることで古典的アラン分散（AVAR）が得られ、MVARではない。
• 短期安定性測定では、係数 $ \{1,2,\ldots,\lceil m/2 \rceil,\ldots,2,1\} $ を用いた重み付き平均化が必要であり、三角形窓を保持して $ \tau = m\tau_B $ におけるMVARを正しく推定できる。
• ピケットフェンス法（ゼロクロスの絶対的タイミング）を用いることで、解釈の誤りが生じない明確なデータが得られ、AVAR、MVAR、TotVarの計算が可能になる。
• 本稿は、重要な曖昧さを解消した：原始的なΛ型カウンタデータをAVAR式に投入すると、三角形窓関数の影響によりMVARが得られ、AVARではない。
• 白色位相雑音の場合、重ね合わせ平均化を用いることで不確実性が $ \sigma^2_\nu \propto 1/\tau^3 $ に改善され、古典的な $ 1/\tau^2 $ スケーリングを上回る。