4. 量子测控系统噪声基底测量与来源分析
测量对象与链路配置
噪声基底测量要求测控链路处于“静默”状态:任意波形发生器(AWG)输出关闭或幅度置零,但采样时钟、本振与采集卡保持正常工作,这样采集到的才是链路自身噪声而非信号。本例在虚拟平台上搭建一条从 AWG 输出端经衰减器、合路器到数字化仪输入的直通链路,采样率设为 1 GSa/s,采集长度 1 M 点,输入阻抗 50 Ω,耦合方式 DC,前端量程取 ±0.5 V 以避免削顶。
关键约束有三点:其一,必须确认 AWG 的 output enable 已关闭,否则残余本振泄漏会抬高底噪;其二,数字化仪不能开启内部平均或数字滤波,否则会掩盖真实噪声谱形;其三,采集时长需覆盖至少 10 个工频周期(50 Hz 下即 0.2 s),本例 1 M 点在 1 GSa/s 下仅 1 ms,因此需通过多次触发拼接或降低采样率到 10 MSa/s 采集 0.1 s 来获得工频分辨率。下面采用 10 MSa/s、1 M 点(0.1 s)配置,兼顾工频分辨与白噪声统计。
噪声数据的采集与预处理
采集得到的是时域电压序列 x[n],单位 V。预处理包括三步:去均值(消除直流偏置,否则 PSD 在 0 Hz 处出现虚假尖峰)、去趋势(一阶线性拟合扣除缓慢漂移)、加 Hann 窗(降低频谱泄漏)。去均值后的序列用于计算均方根值:
Vrms = sqrt( (1/N) * Σ (x[n] - mean(x))^2 )
本例示例数据(教学假设值,非实测):去均值后 Vrms = 42.6 µV。若前端量程 ±0.5 V、16 位 ADC,则量化步长约 15.3 µV,量化噪声理论值约 4.4 µV RMS,远低于 42.6 µV,说明量化不是主导项,噪声主要来自模拟链路与外部干扰。
功率谱密度计算与来源分解
用 Welch 法估计单边功率谱密度:分段长度 8192 点、50% 重叠、Hann 窗,频率分辨率 Δf = 10 MSa/s / 8192 ≈ 1.22 kHz。PSD 单位 V²/Hz,取 10·log10 转为 dBV/Hz 便于观察。分解思路如下:
- 工频干扰:在 50 Hz 及其谐波(100、150 Hz…)处出现窄峰,峰宽通常只有几个 Δf。本例因 Δf≈1.22 kHz 无法分辨 50 Hz 峰,故需另用低采样率(如 10 kSa/s)单独采集 1 s 做窄带分析,确认 50 Hz 峰是否存在。
- 白噪声:PSD 在宽频段平坦,其平台高度对应链路热噪声与放大器噪声之和。本例平坦段约 −142 dBV/Hz,折算到 50 Ω 系统即约 0.56 nV/√Hz 量级(示例值)。
- 仪器本底:将数字化仪输入端接 50 Ω 匹配负载(断开前级),重复同样采集与 PSD 计算,得到本底谱。若本底谱比整链谱低 10 dB 以上,说明前级贡献主导;若接近,则仪器本底不可忽略。
判断贡献占比可用频带积分:对白噪声平台在 1 kHz–1 MHz 内积分得总功率,与 Vrms² 比较;工频峰单独积分其窄带功率。本例示例:白噪声积分贡献约 38 µV RMS,工频窄带贡献约 18 µV RMS,仪器本底约 6 µV RMS,三者平方和开根约 42.7 µV,与总 Vrms 42.6 µV 吻合,说明分解自洽。
配置调整建议与验证
依据分解结果给出针对性调整,并用同一测量流程复测验证:
- 工频干扰占比高(本例 18 µV 不可忽略):检查信号线与电源线走向、单点接地、屏蔽层连接;若仍存在,可在数字化仪前加 50 Hz 陷波或改用差分输入。复测时观察 50 Hz 峰是否下降 10 dB 以上。
- 白噪声主导:白噪声来自链路热噪声与放大器,降低前端带宽(如加 1 MHz 低通)可将积分噪声按 √带宽 缩小;若需保留带宽,则更换更低噪声系数的前放。复测时观察平坦段 PSD 是否下降。
- 仪器本底接近整链:说明数字化仪自身噪声不可忽略,应提高前端增益使信号高于本底,或换用更高分辨率采集卡。复测时比较接负载与接前级两种谱的间距是否拉开。