4. 噪声功率谱密度的定义与估计

在量子传感实验中，我们采集到的往往是一条随时间起伏的电压或计数序列，例如SQUID磁强计输出的连续电压、NV色心荧光计数率、原子钟鉴频器的误差信号。这条时域波形本身并不能直接告诉我们“噪声在哪个频段更强”，而量子传感的灵敏度、相干时间、磁噪声底等关键指标，恰恰都由噪声随频率的分布决定。功率谱密度（PSD）就是把时域随机信号转换到频域、用“单位带宽内的功率”来描述噪声强弱的核心工具。本节围绕PSD的定义约定、两种经典估计方法（周期图法与Welch平均法）以及它们在量子传感数据中的适用条件展开。

一、双边PSD与单边PSD：先约定，再计算

设采样序列为 x[n]，采样率 f_s，采样点数 N，采样总时长 T = N/f_s。离散傅里叶变换记为 X[k]，k = 0,1,…,N−1，对应频率 f_k = k·f_s/N。双边功率谱密度的定义为

S_xx(f_k) = (1 / (f_s · N)) · |X[k]|²      （双边，单位 V²/Hz）

双边谱分布在 −f_s/2 到 +f_s/2 上，正负频率各占一半功率。由于实信号的频谱共轭对称，负频率部分不携带新信息，工程上更常用单边谱：把正频率部分的功率加倍（直流与奈奎斯特频点除外），只保留 0 到 f_s/2：

S_xx^one(f_k) = 2 · S_xx(f_k)   (0 < f_k < f_s/2)
S_xx^one(0)   = S_xx(0)
S_xx^one(f_s/2) = S_xx(f_s/2)

这里必须强调一个常见错误：单边与双边相差因子2，若把两种谱混用比较，噪声底会凭空差3 dB。因此在任何报告、图表和后续灵敏度换算中，都要先声明“本文采用单边PSD”。另一个易错点是归一化：除以 f_s·N 得到的是密度（每赫兹），若只除以 N 则得到的是离散功率，量纲和数值都不同。对于量子传感中常见的“计数率”数据，PSD单位是 counts²/Hz；若再除以灵敏度系数（如 V/nT 或 counts/nT），才能换算为 nT²/Hz 或 nT/√Hz 的幅度谱密度。

二、周期图法：最直接，但方差接近100%

周期图法就是把整段数据一次性做FFT，取模平方并按上面的公式归一化。它的优点是实现简单、频率分辨率最高（Δf = f_s/N），能分辨靠得很近的谱线，例如寻找工频50 Hz及其谐波、机械泵的转动频率、或者某个已知振荡源的窄带峰。在量子传感中，周期图常被用来快速“看一眼”数据里有没有明显的干扰峰。

但周期图有一个致命弱点：它是功率谱的非一致估计。即使数据完全平稳，单个频点上的估计值也服从指数分布，相对标准差约为100%，也就是说谱线看起来“毛刺”很重，增加采样点数N并不会让噪声底的起伏变小，只会让频率网格更密。因此直接用周期图读取噪声底（例如判断磁噪声是否达到 10 pT/√Hz）会产生很大偏差，单次读数可能偏离真值数倍。

实操中若必须用周期图，可以做两点补救：一是加窗（如Hann窗）抑制频谱泄漏，避免强干扰峰的能量扩散掩盖邻近的弱噪声；二是对同一段数据做少量分段平均，哪怕只平均4~8段，也能把方差降到1/4~1/8。加窗会带来等效噪声带宽（ENBW）的变化，Hann窗的ENBW约为1.5个频率格，严格计算密度时需除以ENBW进行修正。

三、Welch平均法：量子传感噪声底估计的主力

Welch方法把长度为N的数据分成K段，每段长L，相邻段之间允许重叠（常用50%重叠），对每段加窗、做FFT、计算周期图，最后把所有段的周期图逐点平均。平均后每个频点的相对标准差降为约 1/√K，谱线明显平滑，噪声底读数稳定可重复。代价是频率分辨率从 f_s/N 降为 f_s/L，且重叠会引入段间相关性，使有效平均次数略小于K。

以一个具体的教学示例说明参数选择。假设NV色心磁强计以 f_s = 10 kS/s 采集荧光计数，记录时长 100 s，共 10⁶ 点。若关心 1 Hz 到 1 kHz 的噪声底，可取段长 L = 10⁴（对应 Δf = 1 Hz），50%重叠，则段数 K ≈ 2N/L − 1 = 199，平均后方差约 1/199，谱线非常平滑。若更关心低频（0.1 Hz 附近），则需增大L到 10⁵，Δf = 0.1 Hz，但K降到约19，方差回升到约1/19，需要在“分辨率”和“平滑度”之间权衡。这个权衡是Welch方法的核心操作决策，没有唯一正确答案，取决于你要回答的物理问题。

另一个必须注意的点是趋势项与窗函数。量子传感数据常含缓慢漂移（温度漂移、1/f噪声），若不去趋势，直流和极低频能量会泄漏到邻近频点，抬高整个低频段。标准做法是先减去每段的均值或拟合线性趋势，再加Hann窗。对于1/f噪声主导的频段，Welch平均后的谱在双对数坐标下应呈 −1 斜率，可据此判断数据是否被漂移污染。

四、在量子传感数据中的适用条件与判断依据

两种方法的选择不是偏好问题，而是由数据特性和分析目标决定。周期图法适用于：记录较短、需要最高频率分辨率找窄带干扰峰、或仅做定性筛查的场合。Welch法适用于：需要定量给出噪声底（如灵敏度曲线中的白噪声平台）、记录足够长、且信号在分段时长内近似平稳的场合。

“平稳”是Welch平均成立的前提。若量子传感器的噪声水平在测量过程中明显变化（例如磁屏蔽未稳定、激光功率漂移、或存在间歇性脉冲干扰），把不同时间段的周期图直接平均会掩盖非平稳特征。此时应先做时频分析（短时傅里叶变换或分段查看各段PSD），确认各段噪声底一致后再平均；若发现某段异常，应剔除或单独分析，而不是盲目平均。

最后给出一个可操作的检查清单，用于判断一次PSD估计是否可信：单边/双边约定是否明确并在全流程一致；归一化是否除以 f_s·N（或等效的ENBW修正）；是否去趋势、是否加窗；Welch的段长、重叠率、段数是否记录；平均后噪声底的相对起伏是否与 1/√K 量级相符；在已知白噪声频段，谱线是否平坦。满足这些条件，PSD才能作为后续灵敏度换算和噪声建模的可靠输入。