[论文解读] Discrete simulation of power law noise

Neil Ashby|arXiv (Cornell University)|Mar 25, 2011

Advanced Frequency and Time Standards参考文献 7被引用 7

一句话总结

本文通过修改白相位噪声的频谱并利用傅里叶变换将其转换到时域，提出了一种用于钟和振荡器中幂律噪声的离散模拟方法。关键贡献是采用对称实矩阵的矩阵方法，其特征值可生成所有均方差分布，从而实现对任意平均时间 τ 下观测方差的精确概率预测，结果表明在一般情况下分布并非卡方分布。

ABSTRACT

A method for simulating power law noise in clocks and oscillators is presented based on modification of the spectrum of white phase noise, then Fourier transforming to the time domain. Symmetric real matrices are introduced whose traces--the sums of their eigenvalues--are equal to the Allan variances, in overlapping or non-overlapping forms, as well as for the corresponding forms of the modified Allan variance. Diagonalization of these matrices leads to expressions for the probability distributions for observing a variance at an arbitrary value of the sampling or averaging interval $τ$, and hence for estimating confidence in the measurements. A number of applications are presented for the common power-law noises.

研究动机与目标

• 开发一种可靠的方法，用于模拟原子钟和振荡器中的幂律噪声，且具有可控的谱特性。
• 推导在任意平均时间 τ 下观测到的 Allan 方差测量值的精确概率分布，超越近似方法。
• 建立一个理论框架，通过第二差分算子的矩阵对角化，将时域噪声模拟与方差统计联系起来。
• 使用相同的谱方法，推导出修正 Allan 方差和非重叠 Allan 方差的解析表达式。
• 证明观测方差的概率分布一般不符合卡方分布，仅在特殊情况下成立。

提出的方法

• 使用 N 个独立、均值为零、方差为 σ² 的复正态随机变量生成白相位噪声。
• 通过将频域振幅除以 |fₘ|ᵋ 来修改，以实现所需的幂律谱密度 S_y(f) = h_α f^α。
• 通过强制满足 w₋ₘ = (wₘ)* 条件，确保逆傅里叶变换生成实值时间序列 xₖ。
• 构建对称实矩阵，其迹等于 Allan 方差（重叠或非重叠）以及修正 Allan 方差。
• 对这些矩阵进行对角化，得到定义在时间 τ 下观测到特定方差的概率密度函数的特征值。
• 在复 ω 平面上使用围道积分计算所得概率积分，考虑分支点以及特征值构型带来的相位贡献。

实验结果

研究问题

• RQ1 如何在保持正确谱密度 S_y(f) = h_α f^α 的前提下，准确模拟离散时间序列中的幂律噪声？
• RQ2 在给定平均时间 τ 下，观测到特定 Allan 方差值的精确概率分布是什么？
• RQ3 为何观测方差的概率分布偏离卡方分布？在何种条件下它们会一致？
• RQ4 如何使用相同的谱模拟框架一致地推导出修正 Allan 方差和非重叠 Allan 方差？
• RQ5 第二差分算子的特征值在决定方差测量统计行为方面起什么作用？

主要发现

• 该方法通过 |fₘ|⁻ᵋ 缩放调制傅里叶振幅，成功生成具有精确幂律谱特性的时间序列。
• 由时间序列差分构造的对称实矩阵，其迹等于 Allan 方差，从而实现直接的统计分析。
• 通过复平面上的围道积分推导出观测方差的概率分布，仅当区间下方有奇数个特征值时才贡献。
• 所得概率表达式为实数，由于相位抵消而呈现交替符号，且对特征值的数量和位置具有非平凡依赖性。
• 当存在四个不同特征值时，概率简化为两个带交替符号的积分之和，清晰展示了分布的非卡方特性。
• 推导结果确认第二差分算子的特征值在所有条件下均保持正值，确保了该框架的数学一致性。