c6_pllsim.m

该程序是二阶琐相环的防震实例 仿真结果给出了二阶琐相环的相平面

   w2b=0; w2c=0; s5=0; phivco=0;
   twopi= 2*pi;
   twofs=2*fs;
   G=2*pi*fn*(zeta+sqrt(zeta*zeta-lambda));
   a=2*pi*fn/(zeta+sqrt(zeta*zeta-lambda));
   a1=a*(1-lambda); a2=a*lambda;
   phierror=zeros(1,npts);
   fvco=zeros(1,npts);
   for i = 1:npts
       s1=phin(i)-phivco;
       s2=sin(s1);
       s3=G*s2;
       s4=a1*s3;
       s4a=s4-a2*s5;
       w1b=s4a+w2b;
       w2b=s4a+w1b;
       s5=w1b/twofs;
       s6=s3+s5;
       w1c=s6+w2c;
       w2c=s6+w1c;
       phivco=w1c/twofs;
       phierror(i)=s1;
       fvco(i)=s6/twopi;
   end 
   freqerror=fin-fvco;