锁相环的MATLAB的仿真程序

%Thephaselockedloop(PLL),adjuststhephaseofalocaloscillator%w.r.ttheincomingmodulatedsignal.Inthisway,thephaseofthe%incomingsignalislockedandthesignalisdemodulated.Thisscheme%isusedinPMandFMaswell.%Wewillimplementitbyusingaclosedloopsystem.Controlsystems%techniquesareappliedhere.%**************************************************************%STEPRESPONSEOFTHEFIRSTORDERCLOSEDLOOPTRANSMITTANCEOFPLL%H(S)=1;%SYSTEMTYPENUMBER=1;%THETAo/THETAi(outputphase/inputphase)closeallkv=1;kd=1;dt=0.01t=0:dt:2u=ones(1,length(t))g11=[tf([2*pi*kv*kd],[12*pi*kv*kd])]%itsthetransferfunctiongiveninthehandout[y11t]=lsim(g11,u,t)figureplot(t,y11)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOF1stORDERCLOSEDLOOPTRANSMITTANCE')%***************************************************************%STEPRESPONSEOFTHEFIRSTORDERCLOSEDLOOPERRORTRANSMITTANCEOFPLL%ALLTHEOTHERFACTORSH(S)etcARESAMEHERE%THETAe/THETAi(sameinterp.asabove)g12=[tf([10],[12*pi*kv*kd])]%errortransmittancegiveninthehandout[y12t]=lsim(g12,u,t)figureplot(t,y12)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOF1stORDERCLOSEDLOOPERRORTRANSMITTANCE')%****************************************************************%STEPRESPONSEOFTHEFIRSTORDERCLOSEDLOOPTRANSMITTANCEOFPLL%BETWEENVCOANDINPUTSIGNALPHASE%H(S)=1;%SYSTEMTYPENUMBER=1;%V2/THETAiKd=1;g13=[tf([Kd0],[12*pi*kv*kd])]%vcovoltageandinputsignaltransmittance[y13t]=lsim(g13,u,t)figureplot(t,y13)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOF1stORDERCLOSEDLOOPTRANSMITTANCEB/WVCOANDINPUTPHASE')%********************************************************************%STEPRESPONSEOFTHESECONDORDERCLOSEDLOOPTRANSMITTANCEOFPLL%SYSTEMTYPENUMBER=2;%THETAo/THETAia=3.15zeta=sqrt((pi*kv*kd)/(2*a))omegan=sqrt(2*pi*kv*kd*a)g21=[tf([2*zeta*omeganomegan^2],[12*zeta*omeganomegan^2])][y21t]=lsim(g21,u,t)figureplot(t,y21)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOFSECONDORDERCLOSEDLOOPTRANSMITTANCEOFPLL')%*********************************************************************%STEPRESPONSEOFTHESECONDORDERCLOSEDLOOPERRORTRANSMITTANCEOFPLL%SYSTEMTYPENUMBER=2;%THETAe/THETAig22=[tf([100],[12*zeta*omeganomegan^2])][y22t]=lsim(g22,u,t)figureplot(t,y22)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOFSECONDORDERCLOSEDLOOPERRORTRANSMITTANCEOFPLL')%*********************************************************************%STEPRESPONSEOFTHESECONDORDERCLOSEDLOOPTRANSMITTANCEOFPLL%BETWEENVCOANDINPUTSIGNALPHASE%SYSTEMTYPENUMBER=2;%V2/THETAig23=[tf([kdkd*a0],[12*pi*kv*kd2*pi*kv*kd*a])][y23t]=lsim(g23,u,t)figureplot(t,y23)xlabel('TIMEINSECONDS')ylabel('AMPLITUDE')title('STEPRESPONSEOFSECONDORDERCLOSEDLOOPTRANSMITTANCEB/WVCOANDINPUTPHASE')%**************************************************************************%%PLL.Mfunctionpll%STUDYINGTHEPHASELOCKEDLOOP(PLL)CHARACHTERISTICS%Stepresponseof1storderclosedlooptransmittanceofpllkv=1;kd=1;dt=.01;t=0:dt:2;step=ones(1,length(t));clt_1=tf([2*pi*kv*kd],[12*pi*kv*kd]);[out_1t]=lsim(clt_1,step,t);figuresubplot(3,1,1)plot(out_1);xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof1stOrderTransmittanceofPLL')gridon;%Stepresponseof1storderclosedlooperrortransmittanceofpllclt_e1=tf([10],[12*pi*kv*kd]);[out_e1t]=lsim(clt_e1,step,t);subplot(3,1,2)plot(out_e1)xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof1stOrderErrorTransmittanceofPLL')gridon;%Stepresponseof1stordertransmittancebetweenVCO&Inputclt1=tf([kd0],[12*pi*kv*kd]);[out1t]=lsim(clt1,step,t);subplot(3,1,3)plot(out1)xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof1stOrderTransmittancebetweenVCO&Input')gridon;%Stepresonseof2ndorderclosedlooptransmittanceofplla=3.15;zeta=sqrt((pi*kv*kd)/(2*a));wn=sqrt(2*pi*kv*kd*a);clt_2=tf([2*zeta*wnwn^2],[12*zeta*wnwn^2]);[out_2t]=lsim(clt_2,step,t);figuresubplot(3,1,1)plot(out_2)xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof2ndOrdertransmittanceofPLL')gridon;%Stepresponseof2ndorderclosedlooperrortransmittanceofpllclt_e2=tf([100],[12*zeta*wnwn^2]);[out_e2t]=lsim(clt_e2,step,t);subplot(3,1,2)plot(out_e2)xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof2ndOrderErrorTransmittanceofPLL')gridon;%Stepresponseof2ndordertransmittancebetweenVCO&Inputclt2=tf([kdkd*a0],[12*pi*kv*kd2*pi*kv*kd*a]);[out2t]=lsim(clt2,step,t);subplot(3,1,3)plot(out2)xlabel('Timeinseconds')ylabel('Amplitude')TITLE('StepResponseof2ndOrderTransmittancebetweenVCO&Input')gridon;%IMPLEMENTATIONOFPHASELOCKEDLOOPINAFREQUENCYDEMODULATOR%Signalgenerationt0=.15;%signaldurationts=0.0005;%samplingintervalfc=200;%carrierfrequencykf=50;%modulationindexfs=1/ts;%samplingfrequencyt=[0:ts:t0];%timevectordf=0.25;%requiredfrequencyresolutionc=cos(2*pi*fc*t);%carriersignalm=[2*ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];%Frequencymodulationint_m(1)=0;for(i=1:length(t)-1)%integralofmint_m(i+1)=int_m(i)+m(i)*ts;endu=cos(2*pi*fc*t+2*pi*kf*int_m);%modulatedsignalfiguresubplot(3,1,1)plot(m(1:300))TITLE('ModulatingSignal')gridon;subplot(3,1,2)plot(