【matlab有偿编程】MATLAB仿真二阶锁相环PLL源代码程序.pdf

% Lab #1 - integer-N PLL parameters % % note: this matlab m-file is used to compute the PLL parameters % for the simulink model. you have to design the loop filters % such that the PLL is stable (i.e. choose a suitable damping factor % zeta where possible) % % Author: Walt Bax % Date: November 23, 1998 % % echo off; if ~exist(gc) gc=b; end str=sprintf(graph colour %s: ,gc); buff=input(str,s); if ~isempty(buff) gc=buff; end if ~exist(samples) samples=5000; % set initial # samples for time domain end npts=800; if ~exist(fr) fr=10e6; end str=sprintf(Reference frequency %4.1fMHz: ,fr/1e6); buff=input(str); if ~isempty(buff) fr=buff; end T=1/fr; if ~exist(gmode) gmode=log; end str=sprintf(x-axis (lin|log) %s: ,gmode); buff=input(str,s); if ~isempty(buff) clf; % clear figure for new graph mode gmode=buff; end f=logspace(1,log10(fr-1),npts); % % given PLL loop parameters BW=50e3; Kphi=2; Kv=50e6; N=10; K=Kphi*Kv/N; % % PFD [rad] -- [V] Bpfd=Kphi*[1]; Apfd=[1]; Hpfd=freqs(Bpfd,Apfd,2*pi*f); % % loop filter 1 - analog integrator with compensating zero % % B(s) (Tau2)s + 1 % H(s) = = % A(s) (Tau1)s if (1) zeta=; % zeta specified for this filter Tau2=; Tau1=; Bloop=[Tau2 1]; Aloop=[Tau1 0]; Hloop1=freqs(Bloop,Aloop,2*pi*f); end % % loop filter 2 - single pole low-pass filter % % B(s) 1 % H(s) = = % A(s) (Tau3)s + 1 if (0) Tau3=; Blp=[1]; Alp=[Tau3 1]; Hloop2=freqs(Blp,Alp,2*pi*f); end % % VCO [V] -- [rad] Bvco=2*pi*Kv*[1]; Avco=[1 0]; Hvco=freqs(Bvco,Avco,2*pi*f); % % PLL frequency response Hloop=Hloop2; % choose a loop filter Hopen=Hpfd.*Hloop.*Hvco/N; % PLL open-loop response Hclose=Hpfd.*Hloop.