锁相环Simulink仿真模型

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

2、n()()(sin)sin()(tutututututuoooooiiiiipdlfvco当i=o，两个信号之间的相位差为一固定值，不随时间变化而变化，称两信号相位同步。当io，两个信号的相位差，不论i是否等于o，只要时间有变化， 那么相位差就会随时间变化而变化，称此时两信号不同步。若这两个信号分别为锁相环的输入和输出，则此时环路出于 失锁状态 。当环路工作时，且输入与输出信号频差在捕获带范围之内，通过环路的反馈控制，输出信号的瞬时角频率)(tv便由o向i方向变化，总会有一个时刻使得i=o，相位差等于 0 或一个非常小的常数，那么此时称为相位锁定，环路处于锁定状态 。若达到锁定状态后，输入信号

3、频率变化，通过环路控制，输出信号也继续变化并向输入信号频率靠近，相位差保持在一个固定的常数之内，则称环路此时为 跟踪状态 。锁定状态可以认为是静态的相位同步，而跟踪状态则为动态的相位同步。环路从失锁进入到锁定状态称为捕获状态。其他几个环路工作时的重要概念：快捕带 ：能使环路快捕入锁的最大频差称为环路的快捕带，记为l，两倍的快捕带为快捕范围。捕获带 ：能使环路进入锁定的最大固有频差，用p表示，两倍的捕获带为捕获范围。同步带 ：环路在所定条件下，可缓慢增加固有频差，直到环路失锁，把能够维持环路锁定的最大固有频差成为同步带，用h，oitt)()(oioittt)()()(2h为同步范围。三者关系为：

4、hpl在理想二阶环的情况下，在捕获状态下，评价捕获性能的主要指标为p、l和捕获时间pt。计算式如下：3202/2nppnlt其中，n为自然谐振角频率， 后面将介绍n在设计环路滤波器时，将与 （阻尼系数， 由于考虑到不同对多种输入信号的误差响应和输出响应的影响， 选取使响应曲线最平稳的最佳值0.707）决定滤波器两个参数的大小，仿真中可通过设定快捕带得到n。从这可以看到，pt不仅与环路参数有关，而且与初始频差有关，固有频差越大，则需捕获时间就越长。在同步状态下，重要的指标有稳态相位误差)(e和h，环路锁定后，频差等于0，但稳态相差通常会存在，它反映了环路的跟踪精度，稳态相差越小，跟踪精度越高。理

5、想二阶环条件下，h。3. 锁相环的构成及工作原理从锁相环结构图看到， 其包括鉴相器、 环路滤波器和压控振荡器。3.1. 鉴相器正弦型鉴相器即一乘法器 （有些资料后接 lpf） ，用于检测环路输入信号相位与输出信号相位间的相位误差)(te，设输入输出信号分别为：作如下变换：通过鉴相器后得到，mk为相乘系数，这里为1/2。3.2. 环路滤波器由通过检相器式子看出， 检相器输出包含了和频分量和差频分量，通过环路滤波器，由于其具有低通特性，和频分量将被滤除，输出为振荡器的控制信号)(tuc。记 f(p)为环路滤波器的传递函数，则)()()(tupftudc。)(sin)sin()(tututuiiii

6、i)(cos)(cos()(tuttutuooooo)()()()()(;)()(211ttttttttttttooiooioioioii)()(2sin)()(sin*)2/1 ()()(2121tttttuuktutukooimoim)()(tutuedd3.3. 压控振荡器压控振荡器为电压频率变换器，其瞬时频率为当)(tuc =0时，)(tv=0。瞬时相位可以表示为通过以上分析，得到模拟锁相环的相位模型为：相应的数字锁相环的模型为则 pll 的动态方程为4. 数字锁相环的设计及simulink仿真数字锁相环的设计主要在于环路滤波器和nco 的设计，而鉴相器则为一简单的数字乘法器。下面将主

7、要介绍数字环路滤波器和nco 的设计。)()()(tuktuftcoocv)(*p/)(dt)(dt)()(2t0t0tukttukttcocoov）（kdf(p)n(p)=ko/pkdf(z)n(z)=ko/(z-1)(sin)p(kf)(p)(p1tttee4.1. 数字环路滤波器设计在清楚数字环路滤波器的结构后，数字环路滤波器的系数是设计的主要部分，其结构如下图所示（simulink 仿真图） ：传输函数为：)1/()(/)()(121zcckukuzfdc由 pll 的线性化数字模型得到的传递函数，将n(z)和 f(z)代入得到12()121( )121()2(1)121k kcczk

8、kc zododh zkkcczk kczodod由 pll 的线性化模拟模型得到传递函数并代入n(s)及 f(s)，然后进行双线性变换（)1/()1)(/2(11zztss）得到22 1224() 2()()4( )22214 4() (1) 2()8ttt zttznnnnnh zttztznnn比较两式得到 c1、c2，分别为)(44/()(4*)/1()(44/(8*)/1 (22221tttkkctttkkcnnndonnndo通常dokkk取 1， 取 0.707，n可由自己设定的快捕带得到，t为抽样间隔，经计算然后可以求得两参数。4.2. nco 在介绍 nco 的设计之前先介绍

9、一下dds 算法。4.2.1. dds 算法nco 一般采用数字相位综合技术（dds） ，该技术主要是由时钟驱动读取三角函数表，基于dds 的 nco 结构如下图所示：一个 n 位字长相位累加器的dds 的基本结构图如下所示：以单频信号说明dds 的工作原理，信号为)(cos)2cos()(tutfutsooo为初始相位（即前述信号的相位初始值to）以采样频率sf对信号进行采样，得到离散相位序列ssosofftfktfk/22k2)(o即连续两次采样间的相位增量，控制可控制输出信号的频率。现将正弦函数一周期的相位2进行等分，当用 n 位字长的相相位累加三角函数发生器n 位全加器锁存器三 角 波

10、 形存储器d/a lpf 位累加器时，最小等分量为n2/2，若每次相位增量取，得到的最低频率增量为nmin2/2/ssoftf，若频率控制字为m，则可得到输出信号频率增量为n2/m2/mssft。可以预见，若m 越大，则相位累加幅度就越大，输出频率也就向目标频率变化越快，落到锁相环范围内捕获时间也就越小。4.2.2. 设计原理这样就可以清楚地得到nco 的数学模型。设 nco 的自由振荡频率为of，0)0(2，在相位累加器的字长为n、采样频率为sf确定的情况下，可确定所用 dds 频率控制字的初始值和初始相位分别为soffm/2o和sooffm/22/2n，根据环路工作原理，数字环路滤波器输出

11、的控制电压加到nco 的控制端，来调整输出频率，即当数字环路滤波器输出的数字控制电压为)(kuc时，相应的频率控制字变化量就为)(kumc，nco 输出频率和输出相位为：nn2/2/mfmffsosout)()(kk)()(kukmkkcoo式中，n2/2)(mk，定义n2/2sofk为 nco 的频率控制增益，单位为 rad/(sv)。nco 相当于一相位累加器，即一差分方程，转换到z 域，其传递方程即为)1/()()(112zzkkukoc。由此，便可以构造其仿真模型，仿真图如下所示。5. 仿真模型及所遇到的问题5.1. 无噪声模型锁相环 simulink 仿真图如下所示参数设置如表所示：

12、输入信号频率110e3hz 采样频率300e3hz n2*pi*10e3 k0 2 c1 8.8844e+004 c2 1.3159e+004 simulation time 0.002s 运行模型后得到输入与输出频谱图比较如下：动态看， nco 输出信号品率将从100e3hz快速牵引到 110e3hz，但是有杂波存在，而输出与输入有20db 的差别，所以也可接受。)(kud和)(kuc分别的波形波如下：这样可以较清楚看到捕获时间为0.0001s 左右，理论计算值为1.1256e-005 ，还是存在差别，这个问题还有待研究。总体而言，此仿真已起到了数字锁相环仿真的效果，输出信号跟上了输入信号的

13、相位，并有较好的稳定性，入锁之后能够保持同步。为了更好的看到)(kud和)(kuc的入锁稳定过程图，和更好的达到入锁效果，我们必须修改参数来达到预想效果，新参数设置如下：输入信号频率110e3hz 采样频率300e3hz n2*pi*10e3 k0 2*3e3 c1 8.4424e+004/300e3 c2 1.3159e+004/300e3 simulation time 0.125s )(kud和)(kuc的波形如下：我们这时可以非常清楚的看到在0.122s时达到稳定，此时锁相环快捕入锁。 nco 输出频谱图为：杂波较之前的仿真要轻，所以猜测因为k0 的影响。当完全稳定后 nco 的输出频谱为：问题：频偏改变与入锁时间是否有关系。5.2. 加入噪声的仿真及其结果其模型图为：分别在信噪比 -10db、0db、10db、100db 的情况下进行仿真（四种情况采用改进后的参数设置） 。5.2.1.snr = -10db 在此情况下，环路在极短时间内入锁，但很快又会失锁，nco 输出信号频率非常不稳定，锁相环无法正常工作。5.2.2.snr = 0db 仿真时间为 0.161，)(kud和)(kuc的波形及输入输出信号频谱如下：可以看到，在 0.16s的时候，环路入锁，仿真时间变为inf 时，将保持稳定。而在 -5db 时，在 0.065s左右入锁，并保持稳定。5.2.3.snr