3. 随机过程与功率谱密度基础

在量子传感实验中,我们记录到的噪声电压或光强信号本质上是随机过程的一次实现。要从中提取有物理意义的功率谱密度(PSD),必须先明确这个随机过程满足哪些数学条件,否则算出的谱可能随时间漂移、无法复现,甚至把系统漂移误判为宽带噪声。本节围绕四个核心概念展开:平稳性、各态历经性、自相关函数与维纳-辛钦定理,它们共同构成噪声PSD分析的数学地基。

平稳性:噪声统计特性不随时间平移而改变

设随机过程为 X(t),若其任意有限维联合分布对时间平移不变,则称严平稳。工程上更常用的是宽平稳(弱平稳),只要求前两阶矩满足两个条件:均值 μX = E[X(t)] 为常数,自相关函数 RX(t1, t2) 只依赖时间差 τ = t2 − t1,即 RX(τ) = E[X(t)X(t+τ)]。宽平稳是PSD存在的充分条件之一,也是维纳-辛钦定理成立的前提。

在量子传感中,判断平稳性有直接的实验含义。例如用光电探测器测量金刚石NV色心的荧光强度,若激光功率、温度或磁场在测量期间缓慢漂移,均值就不再是常数,信号属于非平稳。此时直接做FFT得到的谱会把漂移能量堆在低频端,形成看似巨大的1/f噪声,实际却是实验条件变化。实操中应先做去趋势(减去线性或多项式拟合基线),再检查分段均值是否稳定:把长记录切成若干段,若各段均值差异远大于段内标准差除以√N,说明仍存在非平稳成分,需要缩短单次测量时间或加强环境控制。

各态历经性:用一条长记录代替无穷多次实验

平稳性只保证统计特性不随时间变,但理论上PSD需要对随机过程的集合平均求极限,即无穷多个平行实验的平均。实际实验中我们只有一条时间序列,这就需要各态历经性(ergodicity):时间平均以概率1收敛于集合平均。对宽平稳过程,若 RX(τ) 满足 RX(τ) → μX2(τ → ∞),则均值具有各态历经性;自相关函数各态历经还需其协方差随时间差衰减足够快。

各态历经性是“一条长记录能代表整个随机过程”的合法性依据。工程上无法严格证明,但可通过分段一致性检验近似判断:把总时长 T 的记录分成 K 段,分别估计PSD,若各段谱在统计误差范围内一致,则可认为过程近似各态历经。若某段谱明显偏离,往往意味着存在非平稳事件(如一次振动、一次温度跳变),应剔除该段或改用更短的分析窗。

自相关函数:噪声在时间上的“记忆”

自相关函数 RX(τ) 描述同一随机过程在两个时刻取值的相关程度,是连接时域与频域的桥梁。它具有三条实用性质:RX(0) = E[X2] 等于信号平均功率;RX(−τ) = RX(τ) 为偶函数;|RX(τ)| ≤ RX(0)。白噪声的自相关是δ函数,意味着任意两个时刻不相关;而1/f噪声或带限噪声的自相关衰减缓慢,说明相邻采样点存在强关联。

在量子传感中,自相关函数还能帮助识别噪声类型。例如对磁力仪输出做自相关,若 RX(τ) 在几个采样间隔内迅速降到接近零,说明噪声带宽较宽、相关时间短;若衰减很慢甚至呈指数长尾,则提示存在低频漂移或1/f成分。相关时间 τc 可粗略取 RX(τ) 降到 RX(0)/e 的时间,它决定了独立采样所需的最小间隔,也影响PSD估计的频率分辨率。

维纳-辛钦定理:自相关与PSD互为傅里叶变换

维纳-辛钦定理指出:对宽平稳随机过程,其自相关函数 RX(τ) 与功率谱密度 SX(f) 构成一对傅里叶变换:

S_X(f) = ∫_{-∞}^{+∞} R_X(τ) e^{-j2πfτ} dτ
R_X(τ) = ∫_{-∞}^{+∞} S_X(f) e^{+j2πfτ} df

这一定理是PSD估计方法的理论源头。它说明:只要能从有限数据中估计出自相关函数,再对其做傅里叶变换,就得到PSD;反过来,PSD的逆变换给出自相关。实际计算中,常用周期图法(直接对数据做FFT取模平方)作为维纳-辛钦定理的近似实现,因为有限长数据的自相关估计与周期图在数学上等价。理解这一点,就能明白为什么加窗、分段平均会影响谱估计的偏差与方差。

下面这张图把本节四个概念串成一条从“随机过程假设”到“PSD可计算”的逻辑链,并标出每个环节在实验中的检查动作。

从随机过程到功率谱密度:四个数学条件的逻辑链 随机过程 X(t) 噪声电压 / 光强序列 一次实验 = 一条实现 宽平稳性 均值恒定 · R 只依赖 τ 检查:分段均值 / 去趋势 各态历经性 时间平均 ≈ 集合平均 检查:分段谱一致性 自相关函数 R(τ) R(0)=平均功率 衰减快慢 → 相关时间 维纳-辛钦定理 R(τ) ⇄ S(f) 傅里叶变换对 PSD S(f) 周期图 / Welch 单位:V²/Hz 实验判据速查 ① 平稳性:分段均值漂移 < 段内标准误,否则去趋势或缩短单次测量; ② 各态历经性:K 段 PSD 在统计误差内重合,异常段剔除; ③ 自相关:R(τ) 衰减到 R(0)/e 的时间给出相关时间 τc,决定采样间隔与频率分辨率; ④ 维纳-辛钦:周期图是 R(τ) 傅里叶变换的有限数据近似,加窗与分段平均用于控制偏差和方差。

把四个概念落到一次PSD估计流程中

假设我们用NV色心磁力仪采集了一段时长 T = 10 s、采样率 fs = 10 kHz 的电压序列,目标是得到 1 Hz~1 kHz 的噪声谱。可按以下顺序操作:先检查均值是否漂移,若存在缓慢趋势则减去多项式基线;再计算自相关,观察相关时间是否远小于 T,以确认各态历经近似合理;然后采用Welch方法,把数据分成若干段、加Hann窗、段间重叠50%、对各段周期图取平均。段长决定频率分辨率 Δf = fs/N,段数决定谱估计的方差。若分段谱之间差异过大,应回到平稳性检查,而不是盲目增加平均次数。

维纳-辛钦定理要求自相关函数绝对可积,理想白噪声不满足这一条件,其PSD应理解为广义函数。实际测量中探测器带宽总是有限的,带限白噪声的自相关为sinc型,定理仍严格成立。因此分析前应明确探测器带宽,避免把带外噪声混叠进谱内。

掌握这四个概念后,后续讨论周期图偏差、Welch平均、Allan方差与PSD的关系时,就有了统一的判断依据:任何谱估计方法,本质上都是在有限数据下对维纳-辛钦关系的一次近似,而平稳性与各态历经性决定了这次近似是否可信。