pll.txt

仿真了锁相环工作到一定时间后达到锁定状态的过程

程序默认输入信号频率fin=800kHz,压控振荡器固有振荡频率foo=1000kHz。

源代码：

function[tt1,t1,fs,fo]=My_pll(fin,foo,Ad,T,Ao,th0,dt);

global A fi fo

if nargin<1;fin=640;foo=800;Ad=1;Ao=2*pi*fo;T=80/foo;th0=[0 0];th0=[0 0];dt=1/(4*foo);end;

%--------------------------------------------------------------------------

fo=foo;fi=fin;A=Ao*Ad;

dT=1/(10*fi);[t,th]=ode45('solver',0:dT:T,th0);

t1=0:dt:T-dt;

th1=spline(t,th(:,1),t1);th2=spline(t,th(:,2),t1);

f=th2/(2*pi)+fo;

%--------------------------------------------------------------------------

mstr=['锁相环路动态仿真--绿色:瞬时频差Δw(t)||红色:相位差 φе(t)||黄色:零参考线。'];

figure('Units','normal','pos',[0.09 .12 .8 .8],'color','k');%,'color','k'

set(gcf,'NumberTitle','off','Name',mstr);

dtt=f-fi;

t2=0:dt:t1(end/2);

y=sin(0*t2);

dth=pi/2*fi*t1+th1;

for i=1:160

subplot(211);

plot(t1(1:i),dtt(1:i),'g',t1(i),dtt(i),'ro',t2,y,'y:');

axis off;legend('瞬时频差' ,[520 435 20 10]);

subplot(212); 

plot(t1(1:i),dth(1:i),'r',t1(i),dth(i),'go',t2,y,'y:'); legend('相 位 差' ,[520 420 20 10]);

axis off;

pause(0.000001);

uicontrol('pos',[20 20 60 20],'string','关闭','fontsize',12, ...

'callback','close(gcbf)'); 

end

pause(0.001);

m=40;n=length(t1);n2=n-m-250;

T=m*dt;dt=0.1*dt;fs=1/dt;

tt1=0:dt:T; tt2=t1(n2)+tt1-20*dt;

vc=th2/Ao; 

v2=spline(t1(n2:n),vc(n2:n),tt2);

y2=vco(v2,fo,fs);

y12=sin(2*pi*fi*tt2);

nstr=['环路的捕获过程动态演示---绿:vo(t) 红:vi(t) fi='];

nstr=[nstr num2str(fi) '(kHz),fo=' num2str(fo) '(kHz)'];

H2=figure('Name',nstr,'Num','off','Pos',[92 78 815 610],'color','k');

for i=1:210

plot(tt2(1:i),y2(1:i),'g',tt2(1:i),y12(1:i),'-.r',tt2(i),y2(i),'ro',tt2(i),y12(i),'go');

set(gca,'units','pix','pos',[22 150 800 350]);legend('输出信号' ,'输入信号',[520 440 20 10]);

axis off;pause(0.0000001);

uicontrol('pos',[20 20 60 20],'string','关闭','fontsize',12, ...

'callback','close(gcbf)'); 

End