本页目录

实验 II · 噪声、功率谱与锁相检测

对标:随机过程、频谱分析、滤波器与锁相放大器 | 前置:测量与协方差、波动与傅里叶分析 一个小信号难测,通常不是因为“噪声很大”这一句就结束了。噪声有频率分布,滤波器有等效噪声带宽,锁相放大器还把相位信息投影到同相与正交通道。本页固定一个可审计的 PSD 约定,让你先预测带宽如何改变 RMS 噪声,再判断锁相究竟保留了什么、丢掉了什么。

学习层:噪声变小了,信号真的变清楚了吗?

1. 具体谜题:一个 1 mV 信号为什么值得做频率选择?

设传感器输出

\[ v(t)=A\cos(2\pi f_0t+\phi)+n(t), \]

其中 \(A=1.000\ \mathrm{mV}\) 是峰值振幅,\(f_0=137\ \mathrm{Hz}\),噪声在相关频段用平稳白噪声近似。直接测量假定信号频率处增益为 1,噪声等效带宽为 \(B=1000\ \mathrm{Hz}\),因此会积分这一带宽内的噪声。本实验的 \(B\) 是包含信号的理想带通的单边带宽,不自动表示从 DC 到 \(B\) 的低通截止;否则把 \(B\) 调到 \(f_0\) 以下时,信号本身也会被滤掉。锁相放大器则用同频参考作乘法,再用低通保留慢变化的 \(I,Q\) 输出。

看起来“锁相后噪声变小”很诱人,但有三个容易互相打架的问题:

  1. PSD 是每 Hz 的功率,还是已经积分过的 RMS?把两者混用会差很多;
  2. 单极点滤波器的截止频率 \(f_c\) 是不是它真正的噪声带宽?
  3. 参考相位错 \(90^\circ\) 时,信号是消失了,还是只从 \(I\) 通道转移到了 \(Q\) 通道?

2. 先预测:先写答案,再打开频谱账本

本实验默认 \(S_n=1.0\times10^{-8}\ \mathrm{V^2/Hz}\)、直接带宽 \(B=1000\ \mathrm{Hz}\)、锁相时间常数 \(\tau=0.010\ \mathrm s\)、参考相位 \(\phi=0^\circ\)。这里 \(S_n\) 明确采用单边电压 PSD,单位为 \(\mathrm{V^2/Hz}\)。先预测:

  1. 白噪声带宽扩大 4 倍时,RMS 噪声会变为原来的 2 倍、4 倍,还是不变?
  2. \(\tau\) 加倍时,单极点的等效噪声带宽 \(B_{\rm ENBW}\) 会加倍还是减半?
  3. 默认参数下,直接 SNR 更接近 \(0.224\),还是锁相 SNR 更接近 \(0.224\)?
  4. \(\phi=90^\circ\) 时,锁相的 \(I\) 会为零;这是否意味着 \(Q\) 也为零?

先选择预测再揭示图。实验分别画出低通功率响应、I/Q 直流分量、噪声 RMS 与两种 SNR;各图独立标注单位和刻度。拖动相位时,\(I\) 与 \(Q\) 的分工会直接可见。

3. PSD 到 RMS:带宽是噪声积分的上下限

本页使用单边 PSD 定义:对实值噪声,只在 \(f\ge0\) 记谱密度 \(S_n(f)\),并令

\[ \sigma_n^2=\int_0^\infty S_n(f)\,df. \]

白噪声在测量带宽 \(B\) 内近似为常数 \(S_n\),所以

\[ \boxed{\sigma_{\rm direct}=\sqrt{S_nB}.} \]

默认值给出

\[ \sigma_{\rm direct} =\sqrt{(1.0\times10^{-8})(1000)} =3.16228\ \mathrm{mV}. \]

信号的 \(A\) 是峰值振幅,而正弦的 RMS 是 \(A/\sqrt2=0.707107\ \mathrm{mV}\),所以直接 SNR 为

\[ \mathrm{SNR}_{\rm direct} =\frac{A/\sqrt2}{\sqrt{S_nB}} =\frac{A}{\sqrt{2S_nB}} =0.223607. \]

这一步的关键不是记住 \(\sqrt B\),而是知道 这里的 RMS 平方是电压均方值(固定阻抗下与功率成比例),PSD 乘带宽才是这个均方值。若把 \(S_n\) 当成 \(\mathrm{V/\sqrt{Hz}}\) 的幅度谱密度,公式和单位都会错一层。

4. 单极点滤波器:截止频率不是 ENBW

锁相低通的一个教学模型是

\[ H(f)=\frac{1}{1+i2\pi f\tau}, \qquad |H(f)|^2=\frac{1}{1+(2\pi f\tau)^2}. \]

幅度下降到 \(1/\sqrt2\) 的 \(-3\ \mathrm{dB}\) 截止频率为

\[ f_c=\frac{1}{2\pi\tau}. \]

但白噪声通过滤波器后的功率要积分 \(|H|^2\),所以等效噪声带宽是

\[ B_{\rm ENBW} =\int_0^\infty |H(f)|^2\,df =\frac{1}{4\tau} =\frac{\pi}{2}f_c. \]

默认 \(\tau=0.010\ \mathrm s\) 时,\(f_c=15.9155\ \mathrm{Hz}\),而 \(B_{\rm ENBW}=25.00\ \mathrm{Hz}\)。把 \(f_c\) 直接当 ENBW 会低估经过滤波器的白噪声功率;把两者都叫“带宽”却不说明定义,是实验记录中很常见的隐患。

5. 混频与相位:锁相改变的是坐标和带宽

用参考 \(\cos(2\pi f_0t)\) 与信号相乘,三角恒等式给出

\[ \cos(2\pi f_0t+\phi)\cos(2\pi f_0t) =\frac12\cos\phi +\frac12\cos(4\pi f_0t+\phi). \]

本实验同相参考实际取 \(2\cos(2\pi f_0t)\),正交参考取 \(-2\sin(2\pi f_0t)\)。把乘法增益明确写出来后,低通输出的直流分量为

\[ I=A\cos\phi,\qquad Q=A\sin\phi. \]

因此 \(I^2+Q^2=A^2\),这是无噪声直流分量的相位旋转不变量;只看 \(I\) 时,\(\phi=90^\circ\) 会让信号暂时看起来消失,但 \(Q\) 仍为 \(A\)。参考频率偏离 \(f_0\) 时,低通输出会出现随时间旋转的拍频,不能再把它当成稳定的 DC 读数。

对 PSD 的因子也要跟校准一起声明。本实验把输入 \(S_n\) 定义为单边 PSD,并把混频后 \(I/Q\) 通道的增益校准为上式,于是白噪声的周期平均输出方差为

\[ \boxed{\sigma_{\rm LI}^2=2S_nB_{\rm ENBW},\qquad \sigma_{\rm LI}=\sqrt{2S_nB_{\rm ENBW}}.} \]

混频后的周期平均双边谱为 \(S_n^{(2)}(f-f_0)+S_n^{(2)}(f+f_0)=S_n\),故单边基带谱为 \(2S_n\)。这说明因子 \(2\) 来自两边带相加与当前增益,不是任意的校正。若改用双边 PSD \(S_n^{(2)}=S_n/2\),或者把混频器增益归一化到另一种定义,公式中的因子会随之重写;不能只换一个因子而保留其余数字。

默认锁相噪声为

\[ \sigma_{\rm LI} =\sqrt{2(1.0\times10^{-8})(25)} =0.707107\ \mathrm{mV}. \]

相位对齐时 \(I=A\),所以

\[ \mathrm{SNR}_{\rm LI} =\frac{|A\cos\phi|}{\sqrt{2S_nB_{\rm ENBW}}} =1.41421. \]

和直接 SNR 的比值是

\[ \frac{\mathrm{SNR}_{\rm LI}}{\mathrm{SNR}_{\rm direct}} =|\cos\phi|\sqrt{\frac{B}{B_{\rm ENBW}}}. \]

默认 \(\phi=0^\circ\),故这里两个已定义读数的比值为 \(\sqrt{1000/25}=6.32456\)。这个比值包含噪声带宽和读出方式的改变;它不能直接当作相同观测时间下的检测概率改善倍数。

有限低通:为什么“直流分量”不等于整个输出?

单极点滤波器不会把 \(2f_0\) 项彻底切掉。由 \(H(2f_0)\) 可知,稳态同相输出仍含纹波:

\[ I_{\rm sig}(t)=A\cos\phi+ \frac{A\cos(4\pi f_0t+\phi-\arctan(4\pi f_0\tau))} {\sqrt{1+(4\pi f_0\tau)^2}}. \]

因此纹波峰值占输入振幅的比例为 \(\eta=1/\sqrt{1+(4\pi f_0\tau)^2}\)。默认约 5.80%;若 \(f_0=20\ \mathrm{Hz},\tau=0.001\ \mathrm s\),比例约 97.0%,这时将整个输出说成稳定的 \(I=A\cos\phi\) 就会误导。图中相量只表示直流部分,旁边单列纹波;需要更长时间常数、更高阶滤波或专门的同步滤波才能进一步抑制它。

白噪声乘周期参考后也不严格平稳。以 \(h(s)=e^{-s/\tau}/\tau\)(\(s\ge0\))为低通冲激响应,用双边白噪声协方差 \(\langle n(t)n(t')\rangle=(S_n/2)\delta(t-t')\),得

\[ \operatorname{Var} I(t)=2S_n\int_0^\infty h(s)^2\cos^2[2\pi f_0(t-s)]\,ds =\frac{S_n}{2\tau}\left[1+\Re\frac{e^{i4\pi f_0t}}{1+i2\pi f_0\tau}\right]. \]

周期平均正是 \(S_n/(2\tau)=2S_nB_{\rm ENBW}\);瞬时方差的相对起伏幅度为 \(1/\sqrt{1+(2\pi f_0\tau)^2}\)。实验报告这个量,避免把周期平均 RMS 当成任意时刻都相同的高斯噪声标准差。只有载频与低通尺度充分分离时,才适合近似使用相同、互不相关的 I/Q 通道噪声。

无 JavaScript 时的静态读法:采用单边输入电压 PSD \(S_n=1.0\times10^{-8}\ \mathrm{V^2/Hz}\),信号峰值 \(A=1.000\ \mathrm{mV}\),参考频率 \(f_0=137\ \mathrm{Hz}\),直接通道带宽 \(B=1000\ \mathrm{Hz}\),单极点时间常数 \(\tau=0.010\ \mathrm s\),相位 \(\phi=0^\circ\)。白噪声假设平稳且在相关频段近似平坦,参考频率与信号对齐。

账本 默认数值 公式与解释
直接噪声 RMS \(3.16228\ \mathrm{mV}\) \(\sqrt{S_nB}\)
正弦信号 RMS \(0.707107\ \mathrm{mV}\) \(A/\sqrt2\),因为 \(A\) 是峰值
直接 SNR \(0.223607\) \((A/\sqrt2)/\sqrt{S_nB}\)
\(-3\ \mathrm{dB}\) 截止频率 \(15.9155\ \mathrm{Hz}\) \(1/(2\pi\tau)\)
锁相 ENBW \(25.00\ \mathrm{Hz}\) \(1/(4\tau)\),不是 \(f_c\)
锁相噪声 RMS \(0.707107\ \mathrm{mV}\) \(\sqrt{2S_nB_{\rm ENBW}}\)
锁相 \(I\) SNR \(1.41421\) \(\lvert A\cos\phi\rvert/\sigma_{\rm LI}\)
两种 SNR 比值 \(6.32456\) 倍 \(\sqrt{B/B_{\rm ENBW}}\),相位对齐

左上图显示实际低通功率响应 \(|H(f)|^2\) 与等面积矩形:矩形宽度是 ENBW,而 \(f_c\) 处曲线高度为 1/2;响应尾部延伸到图外,ENBW 仍积分至无穷。右上图显示无噪声直流相量;下方把噪声 RMS 与无量纲 SNR 分图比较。把 \(\phi\) 从 0° 增至 90° 时,\(I\) 下降、\(Q\) 上升,但若同时读取两通道,理想的 \(\sqrt{I^2+Q^2}\) 不应因相位旋转而改变。

6. 误区与失败边界:窄带检测不是万能滤镜

  • 混用单边与双边 PSD。实值信号的单边 PSD 把正负频率功率合并到 \(f\ge0\);双边 PSD 则在整个频轴积分。换约定时,数值密度和混频增益必须一起换,不能把 \(2\) 当成可随意调的校正因子。
  • 把 \(\mathrm{V^2/Hz}\) 与 \(\mathrm{V/\sqrt{Hz}}\) 混写。前者乘带宽后得到方差,后者乘 \(\sqrt{\text{带宽}}\) 才得到 RMS。单位检查通常比肉眼看波形更快发现错误。
  • 把截止频率当 ENBW。\(-3\ \mathrm{dB}\) 是幅度/功率响应的一个点,ENBW 是整个 \(|H|^2\) 曲线的面积。不同滤波器阶数、形状因子和数字窗会给不同 ENBW。
  • 认为锁相能救回频率漂移。参考与信号不相干、频率漂移超过低通带宽、相位随时间转动时,\(I\) 会衰减或振荡。应先测参考稳定性和时间常数,再解释低读数。
  • 不区分漂移产生的位置。输入端的低频加性漂移会被混频搬到载频附近;载频附近的边带噪声、信号振幅慢漂移,以及混频后的电子学漂移则可能留在基带。须说明是哪一种,才能讨论缩窄 ENBW 的效果。
  • 把同步检测当作“只看信号”。谐波、参考串扰、地环路、调制副产物和与参考同频的背景也会进入 \(I/Q\)。锁相选择的是相干性,不是“真实信号”的标签。
  • 忽略滤波器瞬态。时间常数越长,稳态噪声越小但建立时间越长;在刚改变频率或相位时立即读数,会把瞬态当作物理响应。

7. 迁移任务:写出能够被检验的实验决定

题 1:设 \(A=0.2\ \mathrm{mV}\),其余保持默认、相位对齐。两种 SNR 各为多少?

题 2:输入端出现 1 Hz 加性漂移,参考仍为 137 Hz;另一种情况下是信号振幅以 1 Hz 缓慢变化。这两种变化在混频后都会停在 1 Hz 吗?

题 3:把时间常数加倍,白噪声的周期平均 RMS 与建立到阶跃终值 99% 的时间各怎样变化?若取 \(f_0=20\ \mathrm{Hz},\tau=1\ \mathrm{ms}\),能否只根据直流 SNR 就认为读数稳定?

展开数值、频率搬移与实验判断

1. 直接 SNR 为 \(0.0447214\),同相直流 SNR 为 \(0.282843\),比值仍为 6.32456。这里只缩小确定性信号,噪声模型未改变。

2. 加性 \(d\cos(2\pi t)\) 乘以参考后变为 136 Hz 与 138 Hz,随后按低通响应衰减。振幅调制 \([A+d\cos(2\pi t)]\cos(2\pi137t)\) 的边带在载频附近,解调后产生真正的 1 Hz 基带变化。混频后的放大器自身漂移也不会被前面的混频移走。不能把“有 1/f 噪声”当成三者相同的理由。

3. ENBW 减半,周期平均噪声 RMS 乘以 \(1/\sqrt2\);单极点阶跃误差为 \(e^{-t/\tau}\),所以 \(t_{99}=\tau\ln100\approx4.605\tau\),等待时间加倍。20 Hz、1 ms 时二倍频纹波约为输入振幅的 97%,还要检查这种确定性波动及瞬时噪声起伏,不能只读一个直流 SNR。

实测时可扫描 \(\tau\),记录 I/Q 均值、输出 PSD、等待时间与参考漂移;用信号关闭的独立噪声数据检查 \(\sqrt{B_{\rm ENBW}}\) 标度。若偏离,先辨认漂移在信号链的位置,再修改模型。没有实际采集数据前,这只是一份实验方案。

1. 噪声的对象:随机过程、方差与频率分布

噪声不是一个固定的电压值,而是许多不可预测成分的统计集合。对零均值平稳过程 \(n(t)\),方差

\[ \sigma_n^2=\langle n^2(t)\rangle \]

描述总平均功率;功率谱密度 \(S_n(f)\) 描述这些功率按频率怎样分配。滤波器把输入变成输出 \(y(t)\),若传递函数为 \(H(f)\),则

\[ S_y(f)=|H(f)|^2S_n(f), \qquad \sigma_y^2=\int S_y(f)\,df. \]

对实值信号,最重要的是先说明积分范围。双边约定在 \(-\infty< f<\infty\) 积分;单边约定把负频率的功率折到 \(f\ge0\)。两套约定都正确,单位相同,而非零频率处的单边谱密度数值为双边的两倍。实验记录应像记录电阻单位一样记录 PSD 约定。

白噪声只是某一频段的近似:\(S_n(f)\approx S_0\)。热噪声、电阻噪声在一定频率范围内可能接近白噪声;闪烁噪声常近似 \(1/f^\alpha\);机械共振和电源纹波则会产生窄峰。把一切都叫“白噪声”会让带宽预测失去适用边界。

2. 滤波器的两种带宽

幅频图上的 \(-3\ \mathrm{dB}\) 点适合描述一个确定频率分量衰减了多少,但随机噪声需要对功率响应积分。若以 DC 增益归一化,等效噪声带宽定义为

\[ B_{\rm ENBW}=\int_0^\infty \left|\frac{H(f)}{H(0)}\right|^2df. \]

更一般地,如果频率响应不平坦,输出白噪声方差仍是 \(S_0|H(0)|^2B_{\rm ENBW}\)。数字滤波器的窗函数、采样率和归一化方式都进入 \(H\),所以“采样率是 1 kHz”不能直接当成噪声带宽。

单极点低通的一阶微分方程是

\[ \tau\frac{dy}{dt}+y=x(t). \]

它的频率响应就是学习层中的 \(H(f)\)。阶数更高的滤波器通常在过渡带衰减更快,但 ENBW 与建立时间、过冲、相位延迟之间仍有取舍。锁相放大器的 time constant 不是一个独立的“精度旋钮”,而是测量等待时间与噪声带宽的共同参数。

3. 同步检测的频域图像

乘以参考信号会把频谱平移。以余弦参考为例,乘法在时域对应频域卷积:

\[ \mathcal F\{x(t)\cos(2\pi f_0t)\} =\frac12[X(f-f_0)+X(f+f_0)]. \]

因此在 \(f_0\) 附近的窄带信号被搬到 DC,原本在 DC 附近的低频漂移则被搬到 \(\pm f_0\) 附近,随后可由低通选择。实际锁相系统还要考虑参考相位、乘法器增益、前置滤波、谐波检测和输出低通。

使用正交参考 \(-2\sin(2\pi f_0t)\) 得到 Q 通道。这个负号与输入约定 \(A\cos(2\pi f_0t+\phi)\) 配套,使低通后 \(Q=A\sin\phi\)。对一个稳定单频信号,\(I,Q\) 是复平面上的坐标:

\[ Z=I+iQ=Ae^{i\phi},\qquad A=\sqrt{I^2+Q^2},\qquad \phi=\operatorname{atan2}(Q,I). \]

这解释了锁相读数的相位依赖:改变参考相位只是在 \(I,Q\) 平面旋转向量。若硬件增益、正交度或通道噪声不匹配,\(I^2+Q^2\) 也会出现系统性变化,因而需要用已知幅度和相位的校准信号检查。

4. SNR 的定义必须指明比较对象

直接通道的正弦信号通常用 RMS 幅度比较噪声 RMS;锁相输出则常用同相峰值或复幅度比较输出噪声。两种定义都可以,但必须说明一个比较正弦 RMS,另一个比较同相直流幅度;不能把两者的比值当成同一估计器精度的直接比较。学习层的直接 SNR 用 \(A/\sqrt2\),锁相 \(I\) SNR 用 \(A\cos\phi\),所以在相位对齐且同一校准约定下比值为

\[ \sqrt{\frac{B}{B_{\rm ENBW}}}. \]

若只读 \(I\),相位失配还会带来 \(|\cos\phi|\) 因子;若同时使用 \(I,Q\) 求幅度,理想相位因子消失,但两个通道的噪声协方差和幅度偏差要进入误差预算。不要把“锁相 SNR”这个短语当成唯一标准量。

还要区分确定性相量与带噪声的幅度估计。即使真实 \(A=0\),若 I/Q 噪声近似独立同分布、零均值高斯、各标准差为 \(\sigma\),则 \(R=\sqrt{I^2+Q^2}\) 服从 Rayleigh 分布,\(E[R]=\sigma\sqrt{\pi/2}>0\)。这是取模的正偏差,不能把每次正的幅度读数都解释成发现信号;有限载频/低通分离时还须检查通道协方差。本实验的相量没有叠加随机噪声,A=0 时画零向量。

SNR 也不是检测概率。给定阈值、观测时间、先验频率范围和误报容忍度后,才可以把 SNR 映射到检测性能。扫描很多频率再挑最高峰会产生 look-elsewhere 效应;单次读数的 SNR 不能直接代表整个搜索的显著性。

5. 实验流程与残余风险

一个可重复的锁相实验可以按以下逻辑组织:先用输入端接 \(50\ \Omega\) 或明确的噪声源测 PSD;再以已知频率和幅度的参考信号校准 \(I/Q\) 增益与相位;然后在信号关闭时扫描时间常数,验证噪声 RMS 在白噪声区间的 \(\sqrt{B_{\rm ENBW}}\) 标度;最后打开信号并报告稳态等待时间、直接带宽、ENBW、参考相位和背景。

若实测曲线不服从 \(1/\sqrt{\tau}\),这不是必须“修正拟合”的坏消息,可能正是新机制:\(1/f\) 噪声在低频占主导、参考抖动把信号移出通带、输出放大器有自身噪声、滤波器存在数字截断,或者测量尚未进入稳态。保留这些失败边界,实验结论才有迁移能力。


下一页:光谱与成像会把同一套“响应函数、带宽与反演”带到空间和波长域。

单极点 ENBW、噪声带宽和仪器输出校准可对照 SRS SR860 手册的锁相基础部分。本页使用峰值校准;实际仪器可能显示输入等效 RMS,换算时须连同增益一起处理。