简介

在实际采集过程中,往往会遇到传感器输出信号较为微弱,并且传感器输出信号受噪声影响较大,利用普通放大滤波电路并不能将有效的待测信号直接提取出来,同时也会引入其他噪声信号对实际测量结果造成严重的偏差。针对这一现象,锁相放大器的优势变的尤为突出。锁相可以将待测信号中目标频率信号分离出来,并抑制噪声信号,减少干扰。

锁相检测的核心在于,利用一个已知的、纯净的参考信号,去检测一个未知的、通常淹没在强噪声中的目标信号,其前提是目标信号与参考信号在频率上相关(通常是同频或倍频)。锁相检测的实现可以分为相关性运算(混频),低通滤波以及直接测量或反馈控制三个步骤。

互相关函数——基本原理

微弱信号检测的基础在于被测信号在时间上具有前后相关性。所谓“相关”,是指两个函数之间存在一定的关联关系,根据概念可分为自相关与互相关。在微弱信号检测中,通常采用抗干扰能力更强的互相关检测方法。相关函数则是衡量信号之间线性相关程度的数学工具。

​ 设信号 f1(t) 为被检信号 Vs(t) 和噪声 Vn(t) 的叠加,而 f2(t) 为与被检信号同步的参考信号 Vr(t) ,二者的互相关函数为: 由于噪声 Vn(τ) 和参考信号 Vr(τ) 不相关, 故 Rnr(τ) = 0。上述分析表明,利用参考信号与有用信号的相关性,以及参考信号与噪声之间的独立性,可以通过互相关运算有效抑制噪声的影响。具体而言,将混有噪声的信号与参考信号进行相乘和积分处理,即可从任意强噪声背景中提取出深埋的微弱信号。基于该原理设计的相干检测器构成了锁相放大器的核心部分。

锁相放大器

如前文所述,锁相检测的第三步可分为直接测量或反馈控制两种方式。采用直接测量方式实现的锁相技术称为锁相放大器。其核心为相敏检波器(Phase Sensitive Detector, PSD)结合低通滤波器(Low Pass Filter, LPF),用于提取电信号的振荡幅值与相位信息。

锁相放大器原理图

相敏检波器主要用于提取信号中特定频率成分的幅值与相位。其原理是将输入信号与参考信号相乘,再通过低通滤波器滤除高频成分,从而得到所需信息。实际应用中,输入信号与参考信号通常具有相同的频率。

​ 设输入信号 Sin(t) 为: Sin(t) = Ainsin (ωint + ϕ) + Bin(t) 式中,Ainωinϕ 分别表示输入信号的幅值、角频率和相位;Bin(t) 代表高斯白噪声。

参考信号 Sref(t) 为: Sref(t) = Arefsin (ωreft) 式中,Arefωref 分别表示参考信号的幅度、角频率;参考信号的相位设为 0,且不含噪声;假设参考信号与输入信号在角频率上相等,即 ωref = ωin = ω0

​ 将两路信号输入相敏检波器,输出信号 SPSD(t) 可表示为(利用三角函数的乘积公式:): 该输出项包括三项内容:第一项是输入信号与参考信号的差频分量,近似直流信号;第二项是输入信号与参考信号的和频分量 ;第三项是 ArefBin(t)sin (ω0t) 为为高斯噪声引起的信号漂移。

SPSD(t) 通过一个窄带低通滤波器,可滤除和频分量与噪声分量,仅保留近似直流的差频信号,得到信噪比较高的输出信号(同相分量,I 通道):
​ 实际应用中,Aref 为固定常数,取决于参考信号的强度;而 Ain 的值会随噪声影响而产生变化,导致 Sout(t) 也随之波动,难以保持 ϕ 恒定。因此,发展出了正交锁相放大器,也称为双通道正交 PSD 检测技术。该方法使用相位相差 90° 的两路参考信号,构成两个正交检测通道,可有效避免单路检测中的相位干扰,提高测量的准确性与可靠性。

​ 两路正交参考信号进行相敏检波,再经低通滤波后,将其合并为复数形式,即为系统的总输出得到: 进一步对这两路输出进行幅值与相位分析,可得:

ϕ = arctan (SQ/SI)

由此可见,正交锁相放大器仅需对低通滤波器输出的两路直流信号进行运算,即可同时解析出待测信号的幅值与相位,无需预先知道其与参考信号之间的相位差。该方法能有效抑制噪声与相位波动的影响,显著提高信噪比,从而保证检测过程的高灵敏度与可靠性。

锁相环(Phase Locked Loop,PLL)

锁相环(PLL)是一种利用相位负反馈控制技术的电路系统,能够使输出信号的相位和频率与参考信号保持同步。其核心原理是通过比较输入参考信号与压控振荡器(VCO)输出信号之间的相位差,生成误差电压,经滤波后控制VCO的频率,实现两者的相位锁定。

闭环的目的:实现输出信号精确跟踪输入信号的变化。当输入信号保持稳定时,锁相环能够有效抑制内部VCO频率漂移和外部干扰引起的输出变化,维持系统稳定。

电路构成

锁相环路的组成原理

该系统由三个基本模块组成:

  • 鉴相器(Phase Detector, PD):检测输入参考信号与VCO输出信号之间的相位差,并将其转换为误差电压信号。可采用模拟乘法器实现,其输出信号 vd 表示为:

vd = Kdf(ϕr − ϕ1)

其中,f 为参考信号与VCO输出信号相位差的函数,Kd 为鉴相器的转换增益,单位为伏特/弧度(V/rad)。

  • 环路滤波器(Loop Filter, LF):通常为低通滤波器,用于滤除鉴相器输出中的高频噪声和无用谐波分量,提高系统稳定性。其输出电压记为 vf

  • 压控振荡器(Voltage Controlled Oscillator, VCO):输出信号的频率受环路滤波器输出电压 vf 的控制,其角频率控制特性为:

ω1 = ω0 + K0vf

其中,ω0 为VCO的中心角频率(对应 vf = 0 时),K0 为VCO的转换增益,单位为弧度/秒/伏(rad/s/V)或赫兹/伏(Hz/V)。

工作原理

  • 鉴相器(Phase Discriminator,PD)

​ 假设锁相环的输入信号为 vi(t) = Vimsin [ωrt + ϕi(t)] ;VCO 的输出信号为 vo(t) = Vomcos [ωrt + ϕo(t)] ,其中 ωr 为VCO的自由振荡角频率。

​ 鉴相器的作用是产生一个平均输出电压 vd,使其值正比于两输入信号的相位差,即 vd = Adsin ϕe(t)Ad 为鉴相灵敏度(单位:V/rad)。

输入与输出信号的瞬时角频率分别为: 其中,相位误差:ϕe(t) = ϕi(t) − ϕo(t) = ∫0t[ωi(t) − ωo(t)]dt ;瞬时角频差:


示例:模拟乘法器鉴相

采用“模拟乘法器 + 低通滤波器”实现鉴相,输出为: 经低通滤波后,得到 vd = Adsin ϕe(t)

鉴相器的特性图

特点:(1)输出为正弦波,频率为 Δωe = ωi(t) − ωo(t) ;(2)鉴相的特性是非线性的,当 ϕe ≤ π/6 时,正弦鉴相可以近似为线性;(3)鉴相的范围是 π/2 ∼ π/2 ;(4)鉴相的灵敏度在数值上等于其最大输出电压


  • 环路低通滤波器(Loop Filter,LF)

​ 环路滤波器的主要作用是滤除鉴相器输出中的高频成分和噪声,保留误差电压的低频或直流分量,以控制VCO。通过使用不同的电路,可选择不同的滤波特性 AF(p),从而对环路的动态性能和稳定性产生影响。常见的环路滤波器有如下的形式:简单RC滤波器、无源比例积分滤波器、有源比例积分滤波器。

​ 理想低通滤波器:对于高于截止频率的信号截止,而对于低于截止频率的信号无失真传输。其传递函数为: H(jω) = |H(jω)|ejφ(ω) 其中 j 是虚数单位,其幅频与相频特性为: |H(jω)| = u(ω + ωc) − u(ω − ωc)

φ(ω) = [u(ω + ωc) − u(ω − ωc)] • (−ωt0)

冲激响应即是传递函数的傅里叶逆变换,为:h(t) = ωc/π ⋅ Sa[ωc(t − t0)],其中 Sa(x) = sin x/x

  • 压控振荡器(Voltage Controlled Oscillator)

​ 压控振荡器的输出频率受控制电压 vc 的控制而变化,一般情况下认为 ω0 ∼ vc 特性是线性的,有: ω0(t) = ωr + Aovc(t) 其中,Ao 为压控灵敏度,其单位是 (rad/s)/V 。VCO输出信号的相位为频率的积分: φ0(t) = Ao0tvc(t)dt 若用微分算子表示 ,可得:pφ0 = Aovc(t)

环路的基本公式

PLL的传递函数

锁相环的相位控制关系可由以下基本方程描述: 其中,AF(p) 为环路滤波器的传递函数。

定义以下角频率关系:

  • 瞬时角频差Δωe(t) = ωi − ωo = pϕe(t)
  • 控制角频差Δωo(t) = ωo − ωr = UdAoAF(p)sin ϕe(t)
  • 固有角频差Δωi(t) = ωi − ωr = pϕi(t)

三者满足:Δωe(t) = Δωi(t) − Δωo(t) 。该式表明,锁相环通过自动调节控制频差 Δωo(t),使其逐渐逼近固有频差 Δωi(t),从而使瞬时频差 Δωe(t) 趋于零,实现相位锁定。

对环路基本方程的讨论

  1. 描述了输入信号与输出信号间的相位(或角频率)的关系,但环路中间是以信号幅度形式传递的。
  2. 是非线性微分方程,其非线性主要来自鉴相器 ,阶数是环路滤波器的阶数加 1。

结合锁相检测的力学传感器

以锁相环设计力传感器可以实现噪声下的信号测量。考虑在外力 fex(t) 作用下, 质量为 m 的机械振子的运动方程为: 其中 x(𝑡) 是位置随时间的函数,γ = ω0/Q 是阻尼系数, 是固有频率,𝑓ex(𝑡)​ 是外力函数。运动方程的解为: x(t) ∝ X(t)(eiω0t + ⋅ ⋅ ⋅), 其具体形式受外力决定。X(t) 是振动幅度,随时间衰减 X2(t) = X2(0)eγtX(0) 为振子的初始振幅。

将时域的运动方程进行傅立叶变化,可以得到振子运动方程在频域的表达式: 其中,ex(ω)fex(t) 的傅立叶变换。解出 (ω)ex(ω) 的响应关系: 上式就是位移对外力的响应,其中 就是机械振子的机械磁化率,记作 χm(ω),即 机械磁化率描述了机械振子位移对施加于机械振子上的外力的响应能力。

傅里叶变换推导(点击展开)

若我们在时间 t 内测得机械振子的时域位移为 x(t),通过傅里叶变换可获得其频域位移 (ω),同理可以计算出 (t) 的傅里叶变换,然后将其代入运动方程就可以得到运动方程在频域的表达方式。

  • x(t) 的傅立叶变换和傅立叶逆变换:

(ω) = ∫0te−iωtx(t)dt

  • (t) 的傅立叶变换:

  • 的傅立叶变换:

在共振条件下,即当 ω ≈ ω0 时,系统的响应会达到最大。当阻尼较小时,幅值响应 |X| 在共振附近会急剧增大,因此我们可以定义共振下的机械磁化率的幅值为: 振子响应 (ω) 与驱动力 ex(ω) 之间的相位差为: 其中,约等于公式皆处于近共振条件。锁相环通过鉴相器(本质原理是锁相放大器)测量振子响应与驱动力之间的相位差,从而跟踪锁定振子的共振频率,令 θ = ϕ − π/2 ,则在近共振条件下: dθ/dω|ω = ω0 = 2Q/ω0 ,通过上式可以得知相位噪声与频率误差之间存在如下的关系: 进一步考虑相位噪声 Δθ ,若仅考虑热噪声,单带热噪声的力学功率谱(热朗之万力PSD)为: ffth(ω) = 2Sffth(ω) = 4mγkBTen 根据前面推导的各个 PSD 之间的转化关系给出热朗之万力加速度 PSD 与 位移 PSD:

xxth(ω) = |χm(ω)|24γkBTen/m

热噪声功率谱推导(点击展开)

考虑外力只有热朗之万力 fth(t) 的情况,热朗之万力是随机力,它具有随机力的两条性质:一是相对于时间平均值为 0,二是白噪声。 fth(t)⟩ = 0;  ⟨fth(t)fth(t)⟩ = 2mγkBTenδ(t − t) 其中 Ten 是环境温度。热朗之万力的力 PSD 根据下式计算: 其中 th(ω)fth(t) 的傅里叶变换。

热噪声位移涨落与积分推导(点击展开)

在不考虑测量时间时,可以算出热噪声的位移涨落满足: 若已知测量时间 Tmea ,则热噪声位移涨落为: 上式的积分运算需要用到的积分结果:

积分结果的推导【留数定理】

f(ω) = [(ω02 − ω2)2 + ω2γ2]−1 ,函数的四个奇点是方程 (ω02 − ω2) + ω2γ2 = 0 的四个根,有: γ ≥ 2ω0 时,得出: 四个一阶奇点分别为: 选择上半复平面或下半复平面作为积分区间,则有: 其中, Resf(ω0) = limω → ω0(ω − ω0)f(ω)。最终可得出: γ ≤ 2ω0 ,结果一致,得出结论:无论噪声振子处于欠阻尼还是过阻尼区域,都有上述结果。

假设鉴相器施加振幅为 FL ,频率为 ω 的驱动信号时,锁相的 IQ 的两路振幅为(热噪声的相位紊乱,对两者的平均值并无贡献):(两路相互正交的信号,其中一路为I路(In-phase,同相位),另一路为Q路(Quadrature,正交相位)) 其中,IQ 的涨落为: ΔI = ΔQ = ⟨x2 注意:相位与 IQ 的关系为: 因此,热噪声对频率的涨落为: 或者不考虑测量时间的频率涨落: Δω = (kBTen/m)(ω0/2Q2)FL−1 这表明:随着鉴相器施加的振幅 FL 的增大,热噪声造成的频率涨落将越小。