锁相环Simulink仿真模型

锁相环学习总结通过这段的学习，我对锁相环的一些基本概念、结构构成、工作原理、主要参数以及simulink 搭建仿真模型有了较清晰的把握与理解，同时，在仿真中也出现了一些实际问题，下面我将对这段学习中对锁相环的认识和理解、设计思路以及中间所遇到的问题作一下总结：1. 概述锁相环（PLL ）是实现两个信号相位同步的自动控制系统，组成锁相环的基本部件有检相器（PD ）、环路滤波器（LF ）、压控振荡器（VCO ），其结构图如下所示：2. 锁相环的基本概念和重要参数指标锁相是相位锁定的简称，表示两个信号之间相位同步。

若两正弦信号如下所示：相位同步是指两个信号频率相等，相差为一固定值。

)(sin )sin()()(sin )sin()('t U t U t u t U t U t u o o o o o i i i i i θθωθθω=+==+=当i ω=o ω，两个信号之间的相位差 为一固定值，不 随时间变化而变化，称两信号相位同步。

当i ω≠o ω，两个信号的相位差 ，不论i θ是否等于o θ，只要时间有变化，那么相位差就会随时间变化而变化，称此时两信号不同步。

若这两个信号分别为锁相环的输入和输出，则此时环路出于失锁状态。

当环路工作时，且输入与输出信号频差在捕获带范围之内，通过环路的反馈控制，输出信号的瞬时角频率)(t v ω便由o ω向i ω方向变化，总会有一个时刻使得i ω=o ω，相位差等于0或一个非常小的常数，那么此时称为相位锁定，环路处于锁定状态。

若达到锁定状态后，输入信号频率变化，通过环路控制，输出信号也继续变化并向输入信号频率靠近，相位差保持在一个固定的常数之内，则称环路此时为跟踪状态。

锁定状态可以认为是静态的相位同步，而跟踪状态则为动态的相位同步。

环路从失锁进入到锁定状态称为捕获状态。

其他几个环路工作时的重要概念：快捕带：能使环路快捕入锁的最大频差称为环路的快捕带，记为L ω∆，两倍的快捕带为快捕范围。

捕获带：能使环路进入锁定的最大固有频差，用P ω∆表示，两倍的捕获带为捕获范围。

同步带：环路在所定条件下，可缓慢增加固有频差，直到环路失锁，把能够维持环路锁定的最大固有频差成为同步带，用H ω∆，o i t t θθθθ-=-)()('o i o i t t t θθωωθθ-+-=-)()()('2H ω∆为同步范围。

三者关系为：H P L ωωω∆<∆<∆在理想二阶环的情况下，在捕获状态下，评价捕获性能的主要指标为P ω∆、L ω∆和捕获时间P T 。

计算式如下：3202/2nP P n L T ξωωωξωω∆=∞=∆=∆其中，n ω为自然谐振角频率，后面将介绍n ω在设计环路滤波器时，将与ξ（阻尼系数，由于考虑到不同ξ对多种输入信号的误差响应和输出响应的影响，选取使响应曲线最平稳的最佳值0.707）决定滤波器两个参数的大小，仿真中可通过设定快捕带得到n ω。

从这可以看到，P T 不仅与环路参数有关，而且与初始频差有关，固有频差越大，则需捕获时间就越长。

在同步状态下，重要的指标有稳态相位误差)(∞e θ和H ω∆，环路锁定后，频差等于0，但稳态相差通常会存在，它反映了环路的跟踪精度，稳态相差越小，跟踪精度越高。

理想二阶环条件下，∞=∆H ω。

3. 锁相环的构成及工作原理从锁相环结构图看到，其包括鉴相器、环路滤波器和压控振荡器。

3.1. 鉴相器正弦型鉴相器即一乘法器（有些资料后接LPF ），用于检测环路输入信号相位与输出信号相位间的相位误差)(t e θ，设输入输出信号分别为：作如下变换：通过鉴相器后得到，m K 为相乘系数，这里为1/2。

3.2. 环路滤波器由通过检相器式子看出，检相器输出包含了和频分量和差频分量，通过环路滤波器，由于其具有低通特性，和频分量将被滤除，输出为振荡器的控制信号)(t u c 。

记F(p)为环路滤波器的传递函数，则)()()(t u p F t u d c =。

)(sin )sin()(t U t U t u i i i i i θθω=+=)(cos ))(cos()('t U t t U t u o o o o o θθω=+=)()()()()(;)()(2'11t t t t t t t t t t t t o o i o o i o i o i o i i θωθθωθθωθωωωθωωωθωθ+=+=+∆=-=∆+-+=+=)]()(2sin[)]()({sin[***)2/1()()(2121t t t t t U U K t u t u K o o i m o i m θθωθθ+++-=)()(t U t u e d d θ=3.3. 压控振荡器压控振荡器为电压频率变换器，其瞬时频率为 当)(t u c =0时，)(t v ω=0ω。

瞬时相位可以表示为 通过以上分析，得到模拟锁相环的相位模型为：相应的数字锁相环的模型为则PLL 的动态方程为4. 数字锁相环的设计及simulink 仿真数字锁相环的设计主要在于环路滤波器和NCO 的设计，而鉴相器则为一简单的数字乘法器。

下面将主要介绍数字环路滤波器和NCO 的设计。

)()]([)(t u K t u f t c o o c v +==ωω)(*p /)(dt)(dt )()(2t0t0't u K t t u K t t c o c o o v ）（=⎰+=⎰=θωωθ)(sin )p (KF )(p )(p 1t t t e e θθθ-=4.1. 数字环路滤波器设计在清楚数字环路滤波器的结构后，数字环路滤波器的系数是设计的主要部分，其结构如下图所示（simulink 仿真图）：传输函数为：)1/()(/)()(121--+==z C C k u k u z F d c由PLL 的线性化数字模型得到的传递函数，将N(z)和F(z)代入得到12()121()121[()2](1)121K K C C z K K C z o d o d H z K K C C z K K C z o d o d --+-=--++-+- 由PLL 的线性化模拟模型得到传递函数并代入N(s)及F(s)，然后进行双线性变换（)1/()1)(/2(11--+-=z z T s s ）得到22122[4()]2()[()4]()2221[44()](1)[2()8]T T T z T T z n n n n n H z T T z T zn n n ξωωωωξωξωωω--+++-=--++++- 比较两式得到C1、C2，分别为))(44/()(4*)/1())(44/(8*)/1(22221T T T K K C T T T K K C n n n d o n n n d o ωξωωωξωξω++=++=通常d o K K K =取1，ξ取0.707，n ω可由自己设定的快捕带得到，T 为抽样间隔，经计算然后可以求得两参数。

4.2. NCO在介绍NCO 的设计之前先介绍一下DDS 算法。

4.2.1. DDS 算法NCO 一般采用数字相位综合技术（DDS ），该技术主要是由时钟驱动读取三角函数表，基于DDS 的NCO 结构如下图所示：一个N 位字长相位累加器的DDS 的基本结构图如下所示：以单频信号说明DDS 的工作原理，信号为)(cos )2cos()(t U t f U t s o o o θϕ=∆+∏=ϕ∆为初始相位（即前述信号的相位初始值t o ω）以采样频率s f 对信号进行采样，得到离散相位序列ss o s o f f T f kT f k /22k 2)(o ∏=∏=∆∆+∆=∆+∏=θϕθϕθθ∆即连续两次采样间的相位增量，控制θ∆可控制输出信号的频率。

现将正弦函数一周期的相位∏2进行等分，当用N 位字长的相位累加器时，最小等分量为N 2/2∏=δ，若每次相位增量取δ，得到的最低频率增量为N min 2/2/s s o f T f =∏=δ，若频率控制字为M ，则可得到输出信号频率增量为N 2/M 2/M s s f T =∏δ。

可以预见，若M 越大，则相位累加幅度就越大，输出频率也就向目标频率变化越快，落到锁相环范围内捕获时间也就越小。

4.2.2. 设计原理这样就可以清楚地得到NCO 的数学模型。

设NCO 的自由振荡频率为o f ，0)0(2=θ，在相位累加器的字长为N 、采样频率为s f 确定的情况下，可确定所用DDS 频率控制字的初始值和初始相位分别为s o f f M /2o ∏=和s o o f f M /22/2N ∏=∏=∆ϕ，根据环路工作原理，数字环路滤波器输出的控制电压加到NCO 的控制端，来调整输出频率，即当数字环路滤波器输出的数字控制电压为)(k u c 时，相应的频率控制字变化量就为)(k u M c =∆，NCO 输出频率和输出相位为：N N 2/2/M f M f f s o s out ∆+=)()(k k θϕθ∆+∆=)()(''k u K M K k c o o =∆=∆θ式中，N 2/2)(M k ∏∆=∆θ，定义N 2/2s o f K ∏=为NCO 的频率控制增益，单位为rad/(sV)。

NCO 相当于一相位累加器，即一差分方程，转换到Z 域，其传递方程即为)1/()()(11'2---=z z K k u k o c θ。

由此，便可以构造其仿真模型，仿真图如下所示。

5.仿真模型及所遇到的问题5.1.无噪声模型锁相环simulink仿真图如下所示参数设置如表所示：输入信号频率110e3HZ采样频率300e3HZ2*pi*10e3nK0 2C1 8.8844e+004C2 1.3159e+004 Simulation time 0.002s运行模型后得到输入与输出频谱图比较如下：动态看，NCO输出信号品率将从100e3HZ快速牵引到110e3HZ，但是有杂波存在，而输出与输入有20dB的差别，所以也可接受。

) (k ud 和)(kuc分别的波形波如下：这样可以较清楚看到捕获时间为0.0001s左右，理论计算值为1.1256e-005，还是存在差别，这个问题还有待研究。

总体而言，此仿真已起到了数字锁相环仿真的效果，输出信号跟上了输入信号的相位，并有较好的稳定性，入锁之后能够保持同步。

为了更好的看到)(kud 和)(kuc的入锁稳定过程图，和更好的达到入锁效果，我们必须修改参数来达到预想效果，新参数设置如下：输入信号频率110e3HZ采样频率300e3HZn2*pi*10e3K0 2*3e3C1 8.4424e+004/300e3C2 1.3159e+004/300e3Simulation time 0.125s) (k ud 和)(kuc的波形如下：我们这时可以非常清楚的看到在0.122s时达到稳定，此时锁相环快捕入锁。