3. 弱磁信号特性与噪声分析

锁相检测之所以能在强噪声背景中“捞出”微弱信号，前提是我们对信号本身和噪声环境有清晰认识。本节先建立弱磁信号的幅度与频率画像，再逐一分析四类常见噪声的来源与频谱特征，最后给出信噪比的定义、改善思路，并通过一个混合模型的估算任务把概念落到可计算的层面。

一、弱磁信号的幅度范围与频率分布

“弱磁”没有绝对统一的数值边界，工程上通常把待测磁感应强度在皮特斯拉到纳特斯拉量级、且被环境噪声淹没的信号视为弱磁信号。典型对象包括：心磁/脑磁（约 10⁻¹³～10⁻¹² T，即 0.1～1 pT）、金属缺陷漏磁的微弱分量（约 0.1～10 nT）、地磁异常中的小尺度目标（约 1～100 nT）、以及磁通门、SQUID、光泵磁力仪等传感器自身可分辨的微小场变化。需要强调，传感器灵敏度（如 1 pT/√Hz）与信号幅度是两个概念：前者决定噪声底，后者决定待测信息量。

频率分布上，弱磁信号往往集中在低频段。生物磁信号主要能量在 0.1～100 Hz；漏磁检测中，若激励频率为 f₀，缺陷响应常出现在 f₀ 及其谐波附近；地磁异常则接近直流或极低频。这一“低频集中”特性非常关键：它既意味着信号与 1/f 噪声、工频干扰在频域上容易重叠，也意味着单纯靠滤波器难以分离，必须借助调制/解调把信号搬移到噪声更低的频段——这正是锁相检测的用武之地。

二、四类常见噪声

热噪声源于导体中载流子的随机热运动，功率谱密度在很宽频带内平坦（白噪声），单边谱密度为 √(4kTR)，k 为玻尔兹曼常数、T 为绝对温度、R 为等效电阻。它对传感器前置放大电路尤为关键，降低 R、降低温度或压缩带宽都能减小其有效值。

1/f 噪声（闪烁噪声）功率谱密度与频率成反比，S(f) ∝ 1/f，在低频段迅速抬升，是弱磁低频检测的主要障碍之一。其拐点频率因器件而异，磁传感器常在几赫兹到几百赫兹之间。由于 1/f 噪声在低频强、高频弱，把信号调制到较高频率再解调，可显著避开其影响。

工频干扰来自 50 Hz 电力系统及其谐波（100 Hz、150 Hz…），表现为窄带强谱线，幅度常远大于弱磁信号。它可通过屏蔽、差分测量、陷波器抑制，但陷波器会带来相位失真和带宽损失，需谨慎使用。

环境磁噪声包括地磁日变化、城市电磁环境、邻近设备（电机、变压器、电梯）产生的磁场波动，频率覆盖极低频到几十赫兹，幅度可达纳特斯拉甚至更高，且具有非平稳性。这类噪声难以用固定模型描述，通常依靠磁屏蔽室、梯度线圈、参考通道自适应抵消等手段处理。

三、信噪比定义与改善需求

信噪比（SNR）定义为信号功率与噪声功率之比，工程上常用分贝表示：SNR(dB) = 10·lg(P_s/P_n) = 20·lg(A_s,rms/A_n,rms)。对弱磁检测而言，原始 SNR 常低至 −20 dB 甚至 −40 dB，即噪声幅度是信号的十倍到百倍。改善需求来自两方面：一是检测门限，要从噪声中可靠判决目标，通常需要解调后 SNR 达到 10 dB 以上；二是测量精度，SNR 每提高 10 dB，参数估计的方差约下降一个数量级。锁相检测通过窄带化把等效噪声带宽压缩到参考频率附近的极窄范围，从而以“牺牲时间换信噪比”的方式实现改善，理论上 SNR 改善量约为 10·lg(B_wide/B_lock)。

四、实战任务：混合模型信噪比估算

设某弱磁检测场景的混合模型为：待测信号 s(t) = 2 nT·sin(2π·10·t)，即幅度 2 nT、频率 10 Hz 的正弦；叠加白噪声（热噪声等效）功率谱密度 N₀ = 0.5 nT²/Hz（单边），1/f 噪声在 10 Hz 处谱密度为 1 nT/√Hz，另有 50 Hz 工频干扰幅度 20 nT。分析带宽取 0～100 Hz。下面按步骤估算原始 SNR 与锁相后的 SNR。

步骤1  信号有效值：A_s,rms = 2/√2 ≈ 1.414 nT
步骤2  白噪声在 100 Hz 带宽内的功率：P_white = N₀·B = 0.5×100 = 50 nT²
步骤3  1/f 噪声功率（按 S(f)=k/f，k=1 nT²/Hz 在 10 Hz 处）：
        P_flicker = ∫(1/f)df 从 f_low 到 100 Hz，取 f_low=0.1 Hz
                  = ln(100/0.1) ≈ 6.91 nT²
步骤4  工频干扰功率：P_hum = (20/√2)² = 200 nT²（落在带宽内）
步骤5  总噪声功率：P_n ≈ 50 + 6.91 + 200 ≈ 256.9 nT²
步骤6  原始 SNR = 10·lg(1.414²/256.9) = 10·lg(2/256.9) ≈ −21.1 dB
步骤7  锁相后等效噪声带宽取 B_lock = 0.1 Hz（时间常数约 1.6 s）：
        白噪声功率降为 0.5×0.1 = 0.05 nT²
        1/f 噪声在 10 Hz 附近按窄带近似，功率约 1×0.1 = 0.1 nT²
        工频干扰被参考频率选择性抑制，残余按 −40 dB 计：200×10⁻⁴ = 0.02 nT²
        P_n,lock ≈ 0.05 + 0.1 + 0.02 = 0.17 nT²
步骤8  锁相后 SNR = 10·lg(2/0.17) ≈ +10.7 dB
步骤9  SNR 改善量 ≈ 10.7 − (−21.1) = 31.8 dB

上述数值为教学示例中的合理假设，用于演示估算流程，并非实测结果。从计算可以看出三点：第一，工频干扰在原始噪声中占比最大（约 78%），若不先做前置抑制，锁相输出的动态范围会被压缩；第二，1/f 噪声在窄带内因频率低而贡献仍不可忽略，说明参考频率不宜选得过低；第三，SNR 改善主要来自带宽压缩比（100 Hz → 0.1 Hz，约 30 dB），与理论 10·lg(1000) = 30 dB 基本吻合，剩余差异来自 1/f 与工频残余。

完成估算后，应回到设计层面判断：若目标要求解调后 SNR ≥ 10 dB，本例在 B_lock = 0.1 Hz 时刚好达标，但余量仅 0.7 dB，实际工程中宜进一步降低前置工频干扰或适当增大时间常数，并预留 3～6 dB 裕量。这一“先画像、再分解噪声、后估算验证”的流程，是后续锁相参数选择的直接依据。