信號與系統(tǒng)仿真實驗報告6

上機實驗6：連續(xù)系統(tǒng)的復頻域分析驗證試驗1clearall;clc;b=[10-1];a=[1232];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',12);brighten(0.1);grid;legend('零點','極點');brighten(0.1);2clearall;clc;b=[10-1];a=[1232];zplane(b,a);legend('零點','極點');3%穩(wěn)態(tài)濾波法clearall;clc;w=8000;s=j*w;num=[01e46e7];den=[187588e6];H=polyval(num,s)/polyval(den,s);mag=abs(H);phase=angle(H)/pi*180;t=2:1e-6:2.002vg=12.5*cos(w*t);vo=12.5*mag*cos(w*t+phase*pi/180);plot(t,vg,t,vo);grid;text(0.25,0.85,'OutputVoltage','sc');text(0.07,0.35,'InputVoltage','sc');title('穩(wěn)態(tài)濾波輸出');ylabel('電壓(v)');xlabel('時間(s)');4%拉氏變換法symsst;Hs=sym('(10^4*(s+6000))/(s^2+875*s+88*10^6)');Vs=laplace(12.5*cos(8000*t));Vos=Hs*Vs;Vo=ilaplace(Vos);Vo=vpa(Vo,4);ezplot(Vo,[1,1+5e-3]);holdon;ezplot('12.5*cos(8000*t)',[1,1+5e-3]);axis([11+2e-3-5050]);5clearall;clc;num=[1e11];den=[12.5e61e120];[r,p,k]=residue(num,den);%%%%%%%%%%%%以下是運行結(jié)果%>>r%%r=%%0.0333%-0.1333%0.1000%%>>p%%p=%%-2000000%-500000%0%%>>k%%k=%%[]%%>>設(shè)計實驗1_1clearall;clc;symsst;Hs=sym('(s+2)/(s^2+4*s+3)');f1=laplace(diff(Heaviside(t)));f1s=Hs*f1;f1ss=ilaplace(f1s);f1ss=vpa(f1ss,4);subplot(2,2,1);ezplot(f1ss,[1,1+5e-3]);holdon;ezplot('(s+2)/(s^2+4*s+3)',[1,1+5e-3]);axis([11+2e-3-0.250.25]);f2=laplace(Heaviside(t));f2s=Hs*f2;f2ss=ilaplace(f2s);f2ss=vpa(f2ss,4);subplot(2,2,2);ezplot(f2ss,[1,1+5e-3]);holdon;ezplot('(s+2)/(s^2+4*s+3)',[1,1+5e-3]);axis([11+2e-3-0.250.25]);f3=laplace(cos(20*t)*Heaviside(t));f3s=Hs*f3;f3ss=ilaplace(f3s);f3ss=vpa(f3ss,4);subplot(2,2,3);ezplot(f3ss,[1,1+5e-3]);holdon;ezplot('(s+2)/(s^2+4*s+3)',[1,1+5e-3]);axis([11+2e-3-0.250.25]);f4=laplace(exp(-t)*Heaviside(t));f4s=Hs*f4;f4ss=ilaplace(f4s);f4ss=vpa(f4ss,4);subplot(2,2,4);ezplot(f4ss,[1,1+5e-3]);holdon;ezplot('(s+2)/(s^2+4*s+3)',[1,1+5e-3]);axis([11+2e-3-0.250.25]);1_3clearall;clc;b=[1.65-0.331-57690.619080];a=[10.99646397.8121318.110];zs=roots(b);ps=roots(a);subplot(2,1,1);plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',12);brighten(0.1);grid;legend('零點','極點');brighten(0.1);symsst;Hs=sym('(1.65*s^4-0.331*s^3-576*s^2+90.6*s+19080)/(s^6+0.996*s^5+463*s^4+97.8*s^3+12131*s^2+8.11*s)');f1=laplace(diff(Heaviside(t)));f1s=Hs*f1;f1ss=ilaplace(f1s);f1ss=vpa(f1ss,4);subplot(2,1,2);ezplot(f1ss,[1,1+5e-3]);holdon;ezplot('(1.65*