锁相检测技术的原理
简介
锁相放大器技术于 20 世纪 30 年代问世 [1, 2, 3],并于 20 世纪中叶投入商业应用 [4]。这种电子仪器能够在噪声极强的环境中提取出信号的幅值和相位信息(见图 1)。锁相放大器采用零差检测方法和低通滤波技术,通过对照一个周期性的参考信号来测量待测信号的幅值和相位。锁相测量方法能够提取出以参考频率为中心的特定频段内的信号,有效滤除所有其他频率分量。如今,市面上性能最出色的锁相放大器具有高达 120 dB 的动态储备 [5],可以在噪声幅值高达待测信号幅值百万倍的情况下实现精准测量。
几十年来,随着这项技术的不断发展,研究人员已经将锁相放大器应用于诸多不同领域。其中最为重要的用途是作为精密交流电压表和交流相位计、噪声测量装置、阻抗谱仪、网络分析仪、频谱分析仪以及锁相环中的鉴相器。相关的科研领域几乎包括了所有空间尺度和温度条件 ,例如在全日光条件下观测日冕 [6]、测量分数量子霍尔效应 [7],或者对分子中各个原子间的键合特性进行直接成像 [8]。锁相放大器的应用范围极为广泛。它与频谱分析仪和示波器一样,是物理、工程和生命科学领域的各种实验系统中一件极为常用、必不可少的工具。与大多数功能强大的工具一样,要想充分发挥锁相放大器的价值并成功设计各种实验,用户必须充分了解其工作原理及特性。
本文将简要说明锁相放大技术的工作原理和最为重要的测量设置,相关说明将同时从时域和频域角度展开。此外,本文还将详细说明如何利用信号调制,实现在保持较短采集时间的同时提高信噪比 (SNR)。最后,我们还将讨论最近在锁相检测技术领域中出现的创新成果及该技术的最新发展水平。
锁相放大器的工作原理
锁相放大器利用信号的时间相关性将信号从嘈杂背景中提取出来。首先,锁相放大器会将输入信号与参考信号相乘(这一过程有时又称下混频或外差/零差检测),然后通过一个可调低通滤波器对相乘后的结果进行滤波。这种方法称为解调或相敏检测,可以将目标频率的信号分离出来,并滤掉所有其他频率分量。参考信号可由锁相放大器自身生成,也可由外部信号源提供给锁相放大器和实验设备。
参考信号通常为正弦波,但也可以采用其他波形。如果采用纯正弦波进行解调,就可以有选择性地在基频或其任意谐波的频率处进行测量。有些仪器使用方波 [9] 作为参考信号,但这会同时捕捉到信号的所有奇次谐波,因而可能引入系统性的测量误差。
为了帮助大家更好地理解锁相检测技术,下文将分别从时域和频域两个角度来介绍混频和滤波过程。
双相解调
如图 2 所示,在典型试验中,通常使用正弦信号来激励被测设备。锁相放大器将利用设备响应信号 Vs(t) 和参考信号 Vr(t) 来确定幅值 R 和相位 θ。这一过程是通过图 2 (b) 中所示的双相解调电路实现的。输入信号将被拆分成两路,并分别与参考信号及其 90 度相移信号相乘。混频器的输出信号经可调低通滤波器滤波,得到 X 和 Y 两个输出,分别称为同相分量和正交分量。通过下式将笛卡尔坐标 X 和 Y 转换为极坐标,即可轻松得到幅值 R 和相位 θ。
\[\begin{align*}
&R=\sqrt{X^2+Y^2},\\
&\Theta = atan2(Y,X)
\end{align*}\]
请注意,为了使相位角的输出范围覆盖全部四个象限,即 (−π, π],我们使用了 atan2 函数,而不是 atan 函数。
如图 2 (b) 所示,为了使用两个不同的相位对输入信号进行解调,锁相放大器必须将输入信号拆分。与模拟仪器相比,利用数字技术拆分信号可以避免SNR 下降和通道不匹配的问题。
信号混频在时域上的表现
我们可以运用复数,以一种简洁的数学形式表达解调过程中的计算。利用基本三角函数定律
\[cos(x)=\frac{1}{2}e^{+ix}+\frac{1}{2}e^{-ix}\]
我们可以将输入信号 Vs(t) 重新表示为复平面上两个矢量之和,这两个矢量的长度均为 R/√2,以相同角速度 ωs 旋转,但一个顺时针旋转,另一个逆时针旋转:
\[\begin{align*}
V_S(t)&=\sqrt{2}R\cdot cos(\omega_St+\Theta) \\
&= \frac{R}{\sqrt{2}}e^{+i(\omega_st+\Theta)}+\frac{R}{\sqrt{2}}e^{-i(\omega_st+\Theta)}.
\end{align*}\]
从图 3 (a) 和 (b) 可以看出,两个矢量在 x 轴上的投影之和(实部)正好是 Vs(t),而矢量之和在 y 轴上的投影(虚部)始终为零。
双相下混频过程在数学上可以表示为输入信号与复参考信号相乘:
\[V_r(t)=\sqrt{2}e^{-i\omega_rt}=\sqrt{2}cos(\omega_rt)-i\sqrt{2}sin(\omega_rt).\]
(4)
混频后的复信号如下式所示:
\[\begin{align*}
Z(t) &= X(t)+iY(t)=V_S(t)\cdot V_r(t)\\
&= R\left [ e^{i\left [ (\omega_s-\omega_r)t+\Theta \right ]} +e^{-i\left [ (\omega_s+\omega_r)t+\Theta \right ]}\right ],
\end{align*}\]
得出的信号带有两个分量,分别位于信号频率与参考频率之和与两者之差处。如图 3 (c) 所示,复混频相当于有一位身处原点,并以频率 ωr 沿逆时针方向旋转的观察者。
在这位观察者看来,有两个箭头正分别以 ωs - ωr 和 ωs + ωr这两个不同的角速度旋转,并且如果信号频率接近于参考频率,角速度为 ωs + ωr的箭头会旋转得快得多。
随后的滤波过程在数学上可表示为求出这两个时变矢量在一段时间内的平均值,并用尖角括号 ⟨…⟩ 来表示。在滤波过程中,通过设置 ⟨exp [-i (ωs + ωr) t + iθ] ⟩ = 0,将 |ωs + ωr| 处的快速旋转项滤除。解调后的信号平均值变为:
\[Z(t)=R\cdot e^{i[(\omega_s-\omega_r)t+\Theta]}.\]
如果信号同频 ωs = ωr,则该式可进一步简化为:
\[Z(t)=R\cdot e^{i\Theta}.\]
公式 7 代表解调后的信号,也就是锁相放大器的主要输出,其中,绝对值 |Z| = R 表示信号的均方根幅值,辐角 arg(Z) = θ 表示输入信号相对于参考信号的相位。
解调信号 Z(t) 的实部和虚部分别为同相分量 X 和正交分量 Y。这两个分量可通过欧拉公式 exp(iωst) ≡ cos(ωst) + i sin(ωst) 计算得出:
\[\begin{align*}
X&= Re(Z)=\left \langle V_s(t)cos(\omega_st) \right \rangle=R \; cos \; \Theta,\\
Y&= Im(Z)=-\left \langle V_s(t)sin(\omega_st) \right \rangle=R \; sin \; \Theta.
\end{align*}\]
图中所示,ωs = ωr 意味着逆时针旋转的箭头看起来静止不动。而另一个箭头则以两倍于频率的速度 (−2ωs)顺时针旋转,这通常称为 2ω 分量。低通滤波器通常会将 2ω 分量完全滤除。
图 4 为信号在混频和滤波前后通过示波器显示出的不同状态。图 4 (a) 所示为示例正弦信号 Vs 和 Vr 随时间变化的情 况,两个信号频率相同,分别为 ωs 和 ωr。图 4 (b) 中的蓝线表示混频后的信号,主要为 2ω 分量。绿线表示滤波后的信号,仅留下直流分量,其值等于 Vs 的同相幅值 X。如图 4 (c) 所示,如果信号频率与参考频率不同,则混频后的信号将不再是简单的正弦波,并且滤波后的平均值为零,如图 4 (d) 所示。这是一个典型的同步检测示例,它仅提取出与参考频率相干的信号,并滤除所有其他信号。
信号混频在频域上的表现
我们将利用傅里叶变换 [10] ,把视角从时域转向频域。傅里叶变换是线性变换,可将时域中频率为 f0 的正弦函数转换为频域中的狄拉克 δ 函数 δ(f-f0),即在频谱中频率 f0 点处的单个峰值。由于任何周期性信号均可表示为正弦信号与余弦信号的叠加[11],因此对于由少量频谱分量组成的信号,我们通常可以非常直观地理解其变换过程。
图 5 (a) 是一个有噪正弦信号在时域中的表示,图 5 (b) 是该信号经过傅里叶变换后在频域中的表示。该正弦信号在频谱中的 +fs 和 -fs频率处都出现峰值。零频率处较小的峰值是由输入信号中的直流偏移量导致。图 5 (c) 中的蓝线表示混频后的时域信号。其相关频谱如图 5 (d) 所示,该图与图 (b) 基本一致,但向低频方向平移了相当于参考频率 fr 的距离。
图 (d) 中用红色虚线表示的低通滤波过程将只允许频率低于特定滤波器带宽 fBW 的信号通过。(c) 中的红线代表输出信号,包含 (d) 中所示频谱的直流分量与滤波器带宽 |f| < fBW 内的噪声。从图中可以明显看出,滤波器带宽必须远小于信号频率 fs,才能有效抑制输入信号中的偏移量。在下面几节中,我们将进一步讨论如何根据具体的实验条件选择合适的滤波器特性。
低通滤波在频域上的表现
\[Q_{out}(\omega)=H(\omega)Q_{in}(\omega).\]
H(ω) 表示滤波器的传递函数。Qin(ω) 与 Qout(ω) 分别是时域中输入信号 Qin(t) 和输出信号 Qout(t) 的傅里叶变换。
为了完全滤除频谱中不需要的部分,我们可能会想要找到一种理想的滤波器,它能够允许所有低于 fBW 的频率通过(即通带),并彻底滤除任何其他频率(即阻带)。然而遗憾的是,这样理想的“矩形滤波器”根本不可能实现,因为这种滤波器的脉冲响应在时间上从 -∞ 延伸到了 +∞,违背了因果性。我们只能达到与之基本类似的效果,为此我们采用如图 6 所示的 RC 滤波器。这种类型的滤波器在模拟域和数字域都很容易实现。模拟 RC 滤波器的传递函数可以用下式近似表示:
\[H(\omega)=\frac{1}{1+i\omega\tau},\]
其中,τ = RC 称为滤波时间常数,R 为电阻,C 为电容。图 7 (a) 和 (b) 中的蓝线是此传递函数的波特图,表示 20log|H(2πf)| 和 arg[H(2πf)] 与 log(f) 之间的函数关系。
| 阶数 | 时间 | 滚降系数 | 带宽,单位:1/τ | 稳定时间,单位:τ | ||||||||
| n | 常数 τ | dB/oct | dB/dec | f-3dB | fNEP | fNEP/f-3dB | 63.2% | 90% | 99% | 99.9% | ||
| 1 | 1 | 6 | 20 | 0.159 | 0.250 | 1.57 | 1.00 | 2.30 | 4.61 | 6.91 | ||
| 2 | 1 | 12 | 40 | 0.102 | 0.125 | 1.23 | 2.15 | 3.89 | 6.64 | 9.23 | ||
| 3 | 1 | 18 | 60 | 0.081 | 0.094 | 1.16 | 3.26 | 5.32 | 8.41 | 11.23 | ||
| 4 | 1 | 24 | 80 | 0.069 | 0.078 | 1.13 | 4.35 | 6.68 | 10.05 | 13.06 | ||
| 5 | 1 | 30 | 100 | 0.061 | 0.069 | 1.12 | 5.43 | 7.99 | 11.60 | 14.79 | ||
| 6 | 1 | 36 | 120 | 0.056 | 0.062 | 1.11 | 6.51 | 9.27 | 13.11 | 16.45 | ||
| 7 | 1 | 42 | 140 | 0.051 | 0.057 | 1.11 | 7.58 | 10.53 | 14.57 | 18.06 | ||
| 8 | 1 | 48 | 160 | 0.048 | 0.053 | 1.10 | 8.64 | 11.77 | 16.00 | 19.62 | ||
表 1. 具有相同时间常数的 n 阶 RC 滤波器的滤波特性。动态应用通常关注 f−3dB 和稳定时间,而在噪声测量中,确保 fNEP