通信系统的计算机模拟第八讲

通信系统的计算机模拟第八讲1第一页，共四十一页，编辑于2023年，星期三第七讲回顾FIR滤波器：综合技术与滤波器特性根据幅度响应进行的FIR设计根据冲激响应进行的设计FIR滤波器仿真模型的实现FIR滤波器的计算机辅助设计

2第二页，共四十一页，编辑于2023年，星期三第六章案例研究：锁相环与微分方程法锁相环是实现现代通信系统中许多子系统的基本构建模块。锁相环广泛用于频率合成，以完成倍频器与分频器、载波与符号的同步以及相干接收机的实现。微分方程（从而与其对应的系统）可以是非线性的、时变的，或者既是非线性的又是时变的。因此，到本章结束时，我们将会掌握仿真很多种极其复杂的系统所需要的工具。3第三页，共四十一页，编辑于2023年，星期三

6.1基本锁相环概念捕获：非线性跟踪：用一个简单的线性模型进行分析就能给出满意的结果，而不需要仿真。正如我们将看到的，标准的环路参数是根据线性模型来定义的。(6-1)(6-2)4第四页，共四十一页，编辑于2023年，星期三6.1.1锁相环模型第一步是建立鉴相器的模型。鉴相器的特性在很大程度上决定着锁相环的工作特性。最常见的鉴相器模型称为正弦鉴相器其输出与相位差的正弦成正比实际应用中并不需要低通滤波器。稳态时，锁相环的相位差可能为零也可能不为零，这要取决于输入信号和环路滤波器的特性。5第五页，共四十一页，编辑于2023年，星期三正弦鉴相特性(t)=θ(t）PPL输入和VCO输出的相位正交（phasequadrature）鉴相器的输出是相位差的奇函数，这一点则是必需的若式（6-1）和式（6-2）都取余弦函数，而自变量不变，则鉴相器的输出与cos((t)-θ(t))

成正比，是一个关于相位差的偶函数。分辨不出负相位差与正相位差。6第六页，共四十一页，编辑于2023年，星期三锁相环模型cont.乘上环路放大器增益μ后，使用传递函数为F(s)而单位冲激响应为f(t)的环路滤波器。对鉴相器的输出ed(t)进行滤波。这样，VCO的输入为7第七页，共四十一页，编辑于2023年，星期三锁相环模型cont.下一步是将VCO相位偏差与VCO输入关联起来。由定义得，VCO输出的频率偏差与VCO输入信号成正比，这样

Kd-VCO常数Hz/V式中G定义为环路增益，由下式给出滤波器的冲激响应f(t)一直是任意的.8第八页，共四十一页，编辑于2023年，星期三参数合并的影响如果仿真的目的是确定系统输入输出特性或者系统的整体特性，比如锁相环完成相位锁定所需要的时间，合并项是一个有效措施。如果进行的仿真是为了考查系统中各功能模块输入或输出端的波形，定义的各参数不能合并在一起。9第九页，共四十一页，编辑于2023年，星期三6.1.2非线性相位模型

(t)和θ(t）之间的关系跟载波频率完全没关系，在仿真中模型中不需要考虑载波频率。我们需要一个不考虑载波情况下，描述(t)和θ(t）之间的关系的模型。10第十页，共四十一页，编辑于2023年，星期三6.1.2非线性相位模型cont锁相环的非线性相位模型。Sin-非线性它又是一种相位模型：输入相位偏差和VCO相位偏差之间的关系，而不是建立如式（6-1）和式（6-2）所分别表示的实际环路的输入信号与VCO输出信号之间的关系，其输入为带通信号的相位偏差，而非实际带通信号。11第十一页，共四十一页，编辑于2023年，星期三非线性相位模型的表达能力如果已知输入信号的相位偏差、载波频率和信号幅度，那么，就可完全确定式（6-1）。如果已知相位偏差、载波频率和相位幅度(6.2)的VCO输出信号非线性相位模型表示了输入相位和VCO相位间的关系-我们感兴趣的物理量因为环路输入和VCO输出的相位偏差均是低通信号，采样率低12第十二页，共四十一页，编辑于2023年，星期三关于低通滤波器低通滤波器用于滤除鉴相器模型中乘法器所产生的载波二次谐波。滤波器仅是概念模型中的一部分，不会出现在物理器件中。表示VCO模型只是一个积分器，由于积分器是一个低通滤波器（它在f=0时，其增益是无穷大；而在f=1/2πHz时，是单位增益）。13第十三页，共四十一页，编辑于2023年，星期三例6-1建立鉴相器模型的一种通用方法通用方法建立任意特性的鉴相器。这种通用方法通过一个傅里叶级数来表示将鉴相器输出与输入关联起来的函数。Ck-傅立叶系数，ψ(t)=θ（t）-（t)任意精度的模型。级数中仅包含奇数项，从而使ed(t)是相位差的奇函数。例如：14第十四页，共四十一页，编辑于2023年，星期三例6-1Matlab实现式中C‘是C的转置。由于B和C是完全由鉴相器模型来定义的，是固定不变的，应在仿真循环的外面定义它们。最常用的蒙特卡罗方法，通常需要很长的仿真运行时间才能得到结果因而尽可能地使用最有效的算法就变得非常重要。避免常用来求和的环路运算，因而就采用向量形式来计算pdout。15第十五页，共四十一页，编辑于2023年，星期三6.1.3具有复输入的非线性模型通常锁相环模型，使得其环路输入为相位偏差量（t）系统内部使用锁相环时，使其环路输入是表示式（6-1）的复包络低通信号Aexp[j(t)]正弦鉴相器16第十六页，共四十一页，编辑于2023年，星期三6.1.4线性模型与环传递函数若相位差很小，可做如下的线性近似（6-12）则环路方程式（6-8）变为（6-13）锁相环的线性相位模型

17第十七页，共四十一页，编辑于2023年，星期三6.1.4线性模型与环传递函数cont式(6-13)进行拉普拉斯变换，井注意到积分的变换相当于除以s，，而时域上卷积的变换相当于频域相乘（6-13）关联VCO相位和输入相位的传递函数为

（6-15）必须牢记，式（6-15）得出的传递函数是基于线性假设的，严格地说，非线性系统的传递函数是不存在的。线性模型还是有用的。环路参数几乎是按照线性模型来定义的。解析法通常是容易分析，（健全检查）仿真结果。那些被跟踪信号变化速率慢于环路带宽的应用中的跟踪特性。18第十八页，共四十一页，编辑于2023年，星期三6.2一阶环与二阶环锁相环的捕捉和跟踪特性很大程度上取决于环路的阶数。锁相环实现的阶数等于式（6-15）给出的传递函数中有限极点的个数。锁相环实现的阶数比传递函数的极点数大1，这多出来的极点是来自于VCO模型的积分器。一阶和二阶19第十九页，共四十一页，编辑于2023年，星期三6.2.1一阶锁相环对于一阶锁相环，F（s）=1，所以（6-16）将上式代入式（6-8）得（6-17）抽样特性（siftingproperty）对进行积分运算得（6-19）对t求导得微分方程（6-18）20第二十页，共四十一页，编辑于2023年，星期三一阶锁相环cont（6-19）根据相位差ψ

(t)=θ（t）-（t)重写式（6-19），就给出了关联相位差和输入相位的微分方程，即（6-20）为了研究一阶锁相环对时刻t0的一个大小为的fΔHz频率阶跃的响应，我们令（6-22）当时t>t0，式（6-20）变为（6-21）这个式子给出了t>t0时频率差与相位差之间的关系,相平面方程式。21第二十一页，共四十一页，编辑于2023年，星期三相平面方程式或简称为相平面，描述了系统的动态特性

（6-22）每一点的相位差和频率差必须在每一点满足（6-22）这些与时间相关的点称为工作点。相平面的上半部分工作点是自左向右移动而相平面的下半部分是自右向左移动。22第二十二页，共四十一页，编辑于2023年，星期三相平面方程式（6-22）Δψ和Δt分别表示相位差和时间的小增量由于时间总是递增的，有Δt>0，因此，在相平面的上半部分，dψ>0,而在相平面的下半部分，dψ<0

上半部分的相位差是递增的（从左向右移动）；而下半部分的相位差递减的（从右向左移动）稳态工作点只能位于上半相平面和下半相平面之间的边界上。B和点D是稳定的工作点，而点C是不稳定的工作点23第二十三页，共四十一页，编辑于2023年，星期三相平面方程式（6-22）假设对应于频率阶跃的初始工作点是A，由式（6-22）可以看出，若2πfΔ<G，那么，稳态工作点是稳定点B。在这个点上频率差为零，稳态相位差是在条件下式（6-22）的解。若2πfΔ>G

，则对零频率差dψ/dt，式（6-22）无解；对于fΔ

>0，工作点将一直向右边移动；对于fΔ

<0

，工作点将一直向左边移动。环路增益G变为锁定范围。对于一阶环路G也是环路带宽（单位rad/s）。24第二十四页，共四十一页，编辑于2023年，星期三例6-2fΔ=5，2πfΔ=31.42，G=30和G=40G=30,则对零频率差，相位无法锁定，系统将一直振荡25第二十五页，共四十一页，编辑于2023年，星期三例6-2fΔ=5，2πfΔ=31.42，G=30和G=40输入频率偏差d/dt和VCO频率偏差dθ/dt我们已看到，一阶锁相环的锁定范围和环路带宽都由参数G决定的。有大的锁定范围和小的环路带宽一阶环路无法满足这个要求，二阶环路是一个实用的系统。26第二十六页，共四十一页，编辑于2023年，星期三6.2.2二阶锁相环一阶环只有一个参数，所以通过调整环路参数来满足一套给定的工作指标的能力是相当有限的。通过改变环路滤波器来提高锁相环的工作特性和设计能力，这导致了二阶锁相环的出现。对于二阶锁相环，环路滤波器的一般形式是（6-25）λ=0定义为理想的二阶锁相环，它表示包含理想积分器（极点s=0）的一个环路滤波器。典型应用中，λ<<1环路线性模型传递函数，F（s）代人（6-15）27第二十七页，共四十一页，编辑于2023年，星期三二阶锁相环化简若用ξ表示系统的阻尼因子，而ωn用表示系统的固有频率----表征线性二阶系统的参数传递函数的分母通常称为特征多项式把（6-26）改为标准形式而对给定的阻尼因子和固有频率进行锁相环设计，从而就可确定达到和设计值所需的物理参数（这里指G和a）。令式（6-27）中具有相同的s幂次方的项相等，得（6-28）（6-29）假设λ是一个已知的常数，我们可以求解式（6-29）以得出a28第二十八页，共四十一页，编辑于2023年，星期三求解a&G假定λ是常数，求解式（6-29）以得出a，并带入（6-28）由于a是一个实参数，λ必须小于ξ2。ξ2的典型值在1/2左右（通常选择ξ2

=1/2）.而正如前面提及的，λ

<<1。对于理想的二阶锁相环，，而。29第二十九页，共四十一页，编辑于2023年，星期三6.3案例研究：仿真锁相环

6.3.1仿真结构把仿真软件分成几个完全独立的程序通常是有益的。预处理：全部参数（系统参数），和设置管理仿真运行的参数（固有参数）。如滤波器阶数、滤波器类型与带宽、放大器增益、码率、扩频比、载波频率、比特率和信噪比。采样率、稳定时间（确保起始暂态值已衰减到忽略不计所需的时间）要处理的采样个数。要传给后处理器做分析或绘图的数据必须在在仿真过程中保存下来。预处理器必须指定后处理器所需要的数据文件，以便在运行仿真时，保存合适的数据。是用于后处理器对信号图的时间轴进行校准的一个波形采样向量和一个采样时间向量。一旦在预处理器中指定了所有必要的信息，通常把预处理器产生的数据写入一个文件，使得仿真引擎和后处理器都可以调用这个文件。我们把预处理器数据写入MATLAB工作区。仿真引擎读出存储在文件中（或者在工作区中）由预处理器建立的数据、并执行仿真。仿真的目的是为后处理器产生数据以供随后的研究用。这可以采用数值的形式（信噪比、误比特率和编码增益等），也可以是由后处理器进行其他处理的采样值向量。后处理器提取仿真引擎产生的数据并生成用户所需的最终仿真结果。波形图、信号星座图、以Eh/NO为自变量的误比特率函数曲线、所研究系统中某一点的功率谱密度图、眼图和直方图。后处理器显然需要相当程度的图形支持。30第三十页，共四十一页，编辑于2023年，星期三6.3.2仿真除了环路滤波器外，仿真模型是简单明了的（6-25）传递函数不是一个真分数函数。应用长除法得F(s)可由两个传递函数并联来实现。第一个传递函数是常数，而第二个是F1(s)Y1(s)、X1（s）分别表示子滤波器F1(s)输出和输入的拉氏变换,交叉相乘，得(6-35)其时域表示为(6-36)31第三十一页，共四十一页，编辑于2023年，星期三6.3.2仿真--环路滤波器(6-36)32第三十二页，共四十一页，编辑于2023年，星期三环路滤波器—信号流图33第三十三页，共四十一页，编辑于2023年，星期三MatlabCodefori=1:nptss1=phin(i)-phivco; %phaseerrors2=sin(s1); %sinusoidalphasedetectors3=G*s2;s4=a1*s3;s4a=s4-a2*s5; %loopfilterintegratorinputw1b=s4a+w2b; %filterintegrator(step1)w2b=s4a+w1b; %filterintegrator(step2) s5=w1b/twofs; %generatefiteroutputs6=s3+s5; %VCOintegratorinputw1c=s6+w2c; %VCOintegrator(step1)w2c=s6+w1c; %VCOintegrator(step2)phivco=w1c/twofs; %generateVCOoutputphierror(i)=s1; %buildphaseerrorvectorfvco(i)=s6/twopi; %buildVCOinputvectorend34第三十四页，共四十一页，编辑于2023年，星期三仿真结果t0单位阶跃

35第三十五页，共四十一页，编辑于2023年，星期三Result36第三十六页，共四十一页，编辑于2023年，星期三6.3.4仿真的误差源在仿真中会出现许多误差源。物理设备映射到解析模型和从解析模型转换到仿真模型。解析模型很多假设，理想化，物理器件的高精确度仿真时，它们可能是不正确的。通常需要进行实验室