3. 漂移的时域与频域表征方法
漂移是量子测量系统中最难对付的误差之一：它不像白噪声那样可以通过多次平均迅速压低，而是随时间缓慢累积，表现为零点缓慢偏移、灵敏度逐渐下降、长期一致性变差。要补偿漂移，第一步不是“调参”，而是把漂移“看清楚”——用合适的统计量把它从原始时间序列中分离出来，判断它属于哪一类噪声、在哪个时间尺度上占主导、是否可以通过反馈或建模消除。本节围绕四种互补的表征工具展开：自相关函数、功率谱密度、Allan方差，以及由它们导出的漂移率与稳定度指标。
一、从时间序列到统计量：先明确参照约定
设量子测量系统输出一段长度为 N 的采样序列 x[k]，采样间隔为 τ₀，总观测时长 T = Nτ₀。本文所有频率 f 均指相对于采样率的普通频率（单位 Hz），功率谱密度 S(f) 采用单边谱定义，单位随被测量而定，例如磁场测量中为 nT²/Hz。Allan方差 σ²(τ) 中的 τ 为平均时间，与采样间隔 τ₀ 区分开。参照约定一旦确定，全文与图示保持一致，避免把“漂移率”误读为斜率或把“稳定度”误读为标准差。
二、自相关函数：漂移在时域中的“记忆”
自相关函数 R(m) = E{x[k]x[k+m]} 描述序列在不同时间间隔 mτ₀ 上的相似程度。白噪声只在 m=0 处有尖峰，随后迅速衰减到零；而漂移（如随机游走、1/f 噪声）会让 R(m) 在很长的滞后上仍保持正值，衰减极慢。实际计算时通常先去除均值，再对有限长度序列做有偏估计，并注意归一化到 R(0)。如果 R(m) 在几十秒甚至几分钟的滞后上仍显著非零，就说明系统存在慢漂移，单纯增加平均次数无法有效抑制。
三、功率谱密度：漂移在频域中的“指纹”
功率谱密度把总功率按频率分解，是识别噪声类型的核心工具。常见量子测量噪声在低频段呈现不同斜率：白噪声 S(f) ∝ f⁰；1/f 噪声 ∝ f⁻¹；随机游走（频率漂移）∝ f⁻²；线性漂移则表现为极低频处能量集中。工程上常用 Welch 平均周期图法估计 PSD：分段、加窗（如 Hann 窗）、FFT、取模平方、平均。分段长度决定频率分辨率 Δf = 1/(N_seg τ₀)，段数越多方差越小但分辨率越粗。观察 S(f) 在双对数坐标下的斜率，就能判断漂移类型，并据此选择补偿策略。
四、Allan方差：连接时域与频域的桥梁
Allan方差专为振荡器与传感器稳定度分析设计，能有效区分白噪声、闪烁噪声和随机游走。其计算过程为：把序列按平均时间 τ = mτ₀ 分组，求每组均值，再对相邻组均值之差求平方平均：
σ²(τ) = 1/(2(M-1)) · Σ_{i=1}^{M-1} (ȳ_{i+1} - ȳ_i)²
其中 ȳ_i 为第 i 个长度为 τ 的均值，M 为组数。在双对数坐标下，σ(τ) 随 τ 的斜率直接对应噪声类型：斜率 −1/2 为白噪声，斜率 0 为闪烁噪声（1/f），斜率 +1/2 为随机游走。Allan偏差 σ(τ) 的最小值点给出系统的最佳平均时间，超过该点后漂移开始主导，继续平均反而变差。这正是漂移补偿的切入点：在最佳平均时间附近设计反馈带宽。
五、漂移率与稳定度指标：把曲线变成可比较的数字
有了上述曲线，还需要可比较的标量指标。常用指标包括：
- 漂移率：单位时间内的零点偏移量，例如 nT/min 或 ppm/h。可由长时间线性拟合斜率得到，也可由 Allan 偏差在长 τ 段的上升斜率反推。
- Allan偏差 σ(τ)：在指定平均时间下的稳定度，常写作 σ(1 s)、σ(10 s)。
- 最佳平均时间 τ_opt：σ(τ) 最小值对应的 τ，是反馈带宽设计的重要依据。
- 噪声本底：白噪声段在 τ=1 s 处的 σ 值，反映系统短期极限。
以某磁强计示例（合理假设数据）为例：采样率 100 Hz，采集 2 小时。计算得 σ(1 s) ≈ 0.8 nT，σ(10 s) ≈ 0.35 nT，σ(100 s) ≈ 0.5 nT，τ_opt ≈ 20 s。这说明 20 s 以内平均有效，超过 20 s 后漂移开始主导，若反馈环路带宽设为 0.05 Hz（对应 20 s），可在漂移被抑制的同时保留足够带宽。若 PSD 在 0.01 Hz 以下呈 f⁻² 斜率，则说明存在随机游走，需引入积分反馈或定期零点校准。
总结来看，自相关函数回答“漂移有多长的记忆”，功率谱密度回答“漂移在哪些频率占主导”，Allan 方差回答“在什么平均时间下系统最稳”，而漂移率与稳定度指标则把这些信息压缩成可比较、可验收的数字。四者配合，才能为后续的漂移补偿与稳定性控制提供可靠依据。