[论文解读] 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}} $，以分析稳定性度量。
• 分析频率计数器的脉冲响应，表明Π型计数器使用矩形平均窗，而Λ型计数器使用三角形（积分）窗。
• 推导出两种计数器类型的有效权重函数：Π型的 $ w_{\Pi}(t) $ 与Λ型的 $ w_{\Lambda}(t) $，这些函数决定了实际测量的是AVAR还是MVAR。
• 证明Λ型计数器的重叠测量产生三角形形状的有效窗，当输入AVAR公式时，导致MVAR结果。
• 表明对 $ m $ 个连续的Λ型测量值进行平均，可得到梯形窗，其渐近趋近于矩形窗，当 $ m $ 较大时收敛至AVAR。
• 采用小波类比解释方差估计器，表明AVAR可视为以 $ 1/(2\tau) $ 为中心的半倍频程带通滤波器。

实验结果

研究问题

• RQ1 Π型与Λ型频率计数器在有效平均窗与方差估计方面有何不同？
• RQ2 为何将Λ型计数器的数据输入AVAR公式会得到MVAR而非AVAR？
• RQ3 在何种条件下，Λ型计数器的重叠平均会收敛至AVAR而非MVAR？
• RQ4 窗函数形状（矩形与三角形）对方差估计结果有何影响？
• RQ5 为何无法在不增加测量时间 $ \tau $ 的前提下，将Λ型计数器数据转换为AVAR？

主要发现

• 具有重叠测量的Λ型计数器产生三角形形状的有效平均窗，当使用标准AVAR公式时，实际测量结果为改进阿伦方差（MVAR）。
• 对于长期稳定性分析（$ m \gg 1 $），Λ型计数器的平均过程趋近于矩形窗，使得AVAR公式输出经典阿伦方差（AVAR），而非MVAR。
• 在短期稳定性测量中，必须使用系数 $ \{1,2,\ldots,\lceil m/2 \rceil,\ldots,2,1\} $ 的加权平均，以保持三角形窗形状，从而在 $ \tau = m\tau_B $ 处正确估计MVAR。
• 采用栅栏法（picket-fence method）——即零交叉点的绝对定时——可提供无歧义的数据，可用于计算AVAR、MVAR与TotVar，且无解释误差。
• 本文解决了一个关键模糊点：将原始Λ型计数器数据输入AVAR公式时，实际得到的是MVAR而非AVAR，原因在于其底层的三角形窗函数。
• 对于白相位噪声，使用重叠平均时，不确定性按 $ \sigma^2_\nu \propto 1/\tau^3 $ 改善，优于经典 $ 1/\tau^2 $ 的标度。