案例研究：锁相环与微分方程方法 .pptVIP

第6章 案例研究：锁相环与微分方程方法 　　由信号流图可以直接得到仿真代码。在仿真循环中计算出的每个采样值所使用的标号与在信号流图中定义的标号是一致的。例如，代码行s2=sin(s1)定义了表征鉴相器的正弦非线性特性，是s1和s2在流图中进行的定义，仿真循环中所有其它的代码行也采用类似的方法。 附录A给出了二阶锁相环完整的仿真代码。 附录B和附录C分别给出了预处理器和后处理器的MATLAB代码。 第6章 案例研究：锁相环与微分方程方法 二阶锁相环完整的仿真代码： % File: c6_PLLsim.m w2b = 0; w2c = 0; s5 = 0; phivco = 0; % initialize twopi = 2*pi; % define 2*pi twofs = 2*fs; % define 2*fs G = 2*pi*fn*(zeta+sqrt(zeta*zeta-lambda)); % set loop gain a = 2*pi*fn/(zeta+sqrt(zeta*zeta-lambda)); % set filter parameter a1 = a*(1-lambda);a2 = a*lambda; % define constants phierror = zeros(1,npts); % initialize vector fvco = zeros(1,npts); % initialize vector % beginning of simulation loop for i = 1:npts s1 = phin(i) - phivco; % phase error s2 = sin(s1); % sinusoidal phase detector 第6章 案例研究：锁相环与微分方程方法 s3 = G*s2; s4 = a1*s3; s4a = s4-a2*s5; % loop filter integrator input w1b = s4a+w2b; % filter integrator (step 1) w2b = s4a+w1b; % filter integrator (step 2) s5 = w1b/twofs; % generate filter output s6 = s3+s5; % VCO integrator input w1c = s6+w2c; % VCO integrator (step 1) w2c = s6+w1c; % VCO integrator (step 2) phivco = w1c/twofs; % generate VCO output phierror(i) = s1; % build phase error vector fvco(i) = s6/twopi; % build VCO input vector end % end of simulation loop freqerror = fin-fvco; % build frequency error vector % End of script file. 第6章 案例研