概述
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
%File:c6_nltvde.m
w2b=0;w2c=0; %?initialize?integrators
yd=0;y=0; %?initialize?differential?equation
tfinal=50; %?simulation?time
fs=100; %?sampling?frequency
delt=1/fs; %?sampling?period
npts=1+fs*tfinal; %?number?of?samples?simulated
ydv=zeros(1,npts); %?vector?of?dy/dt?samples?
yv=zeros(1,npts); %?vector?of?y(t)?samples
%
%?beginning?of?simulation?loop
for i=1:npts
t=(i-1)*delt; %?time
if t<20
ydd=4*exp(-t/2)-3*yd*abs(y)-9*y; %?de?for?t<20
else
ydd=4*exp(-t/2)-3*yd-9*y; %?de?for?t>=20
end
w1b=ydd+w2b; %?first?integrator?-?step?1
w2b=ydd+w1b; %?first?integrator?-?step?2?
yd=w1b/(2*fs); %?first?integrator?output
w1c=yd+w2c; %?second?integrator?-?step?1
w2c=yd+w1c; %?second?integrator?-?step?2?
y=w1c/(2*fs); %?second?integrator?output
ydv(1,i)=yd; %?build?dy/dt?vector?
yv(1,i)=y; %?build?y(t)?vector
end %?end?of?simulation?loop
plot(yv,ydv) %?plot?phase?plane?
xlabel('y(t)') %?label?x?axis?
ylabel('dy/dt') %?label?y?zxis
%?End?of?script?file.?
%?File:?pllpost.m
%
kk=0;
while kk==0
k=menu('Phase Lock Loop Postprocessor',...
'Input Frequency and VCO Frequency',...?
'Input Phase and VCO Phase',...
'Frequency Error','Phase Error','Phase Plane Plot',...?
'Phase Plane and Time Domain Plots','Exit Program');
if k==1
plot(t,fin,'k',t,fvco,'k')
title('Input Frequency and VCO Freqeuncy')
xlabel('Time-Seconds');ylabel('Frequency-Hertz');pause
elseif k==2
pvco=phin-phierror;plot(t,phin,t,pvco)
title('Input?Phase?and?VCO?Phase')
xlabel('Time-Seconds');ylabel('Phase-Radians');pause
elseif k==3
plot(t,freqerror);title('Frequency Error')
xlabel('Time-Seconds');ylabel('Frequency Error-Hertz');pause
elseif k==4
plot(t,phierror);title('Phase Error')
xlabel('Time-Seconds');ylabel('Phase Error-Radians');pause
elseif k==5
ppplot
elseif k==6
subplot(211);phierrn=phierror/pi;
plot(phierrn,freqerror,'k');grid;
title('Phase Plane Plot');xlabel('Phase Error/Pi');
ylabel('Frequency Error-Hertz');subplot(212)
plot(t,fin,'k',t,fvco,'k');grid
title('Input Frequency and VCO Freqeuncy')
xlabel('Time-Seconds');ylabel('Frequency-Hertz');subplot(111)
elseif k==7
kk=1;
end
end %?End?of?script?file.?
%?File:?pllpre.m
%
clear all %?be?safe
disp(' ') %?insert?blank?line
fdel=input('Enter the size of the frequency step in Hertz > ');
fn=input('Enter the loop natural frequency in Hertz > ');
lambda=input('Enter lambda,the relative pole offset > ');
disp(' ')
disp('Accept default values:')
disp('zeta=1/sqrt(2)=0.707,')
disp('fs=200*fn,and')
disp('tstop=1')
dtype=input('Enter y for yes or n for no > ','s');
if dtype=='y'
zeta=1/sqrt(2);
fs=200*fn;
tstop=1;
else
zeta=input('Enter zeta, the loop damping factor > ');
Fs=input('Enter the sampling frequency in Hertz > ');
tstop=input('Enter tstop,the simulation runtime > ');
end %
npts=fs*tstop+1; %?number?of?simulation?points
t=(0:(npts-1))/fs; %?default?time?vector
nsettle=fix(npts/10); %?set?nsettle?time?as?0.1*npts
tsettle=nsettle/fs; %?set?tsettle
%?The?next?two?lines?establish?the?loop?input?frequency?and?phase
%?deviations.?
fin=[zeros(1,nsettle),fdel*ones(1,npts-nsettle)];
phin=[zeros(1,nsettle),2*pi*fdel*t(1:(npts-nsettle))];
disp(' ') %?insert?blank?line??
%?end?of?script?file?pllpre.m
%?File:?pll2sin.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
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?fiter?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.
function []=pplane(x,y,nsettle)
%?Plots?the?phase?plane?with?phase?in?the?range?(-pi,pi)
ln=length(x);
maxfreq=max(y);
minfreq=min(y);
close %?Old?figure?discarded
axis([-1 1 1.1*minfreq 1.1*maxfreq]); %?Establish?scale
hold on %?Collect?info?for?new?fig
j=nsettle;
while j
i=1;
while x(j)
a(i) = x(j)/pi;
b(i) = y(j);
j = j+1;
i = i+1;
end
plot(a,b,'k')
a=[];
b=[];
x=x-2*pi;
end
hold off
title('Phase-Plane?Plot')
xlabel('Phase?Error/Pi')
ylabel('Frequency?Error?in?Hertz')
grid %?End?of?script?file.
%?File:?ppplot.m
%?ppplot.m?is?the?script?file?for?plotting?phase?plane?plots.?If?the
%?phase?plane?is?constrained?to?(-pi,pi)?ppplot.m?calls?pplane.m.
kz=0;
while kz==0
k=menu('Phase?Plane?Options',...
'Extended?Phase?Plane',...?
'Phase?Plane?mod(2pi)',...
'Exit?Phase?Plane?Menu');
if k==1
phierrn=phierrn/pi;
plot(phierrn,freqerror,'k')
title('Phase?Plane?Plot')
xlabel('Phase?Error/Pi')
ylabel('Frequency?Error-Hertz')
grid
pause
elseif k==2
pplane(phierrn,freqerror,nsettle+1)
pause
elseif k==3
kz=1;
end
end %?End?of?script?file.
运行时遇到如下问题
请问该如何解决?如果要改动程序,请详细指点一下,谢谢!
最后
以上就是开心麦片为你收集整理的matlab 锁相环仿真,MATLAB锁相环仿真程序求解的全部内容,希望文章能够帮你解决matlab 锁相环仿真,MATLAB锁相环仿真程序求解所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复