西工大信号与系统上机实验六

五．实验结果（程序设计实验源程序和结果）1．验证性实验1）．系统零极点的求解>>clear;>>b=[1,0,-1];>>a=[1,2,3,2];>>zs=roots(b);>>ps=roots(a);>>plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',12);>>grid; >>legend('零点','极点');2）.稳态滤波法求解>>w=8000;>>s=j*w;>>num=[0,1e4,6e7];den=[1,875,88e6];>>H=polyval(num,s)/polyval(den,s);>>mag=abs(H);>>phase=angle(H)/pi*180;>>t=2:1e-6:2.002;>>vg=12.5*cos(w*t);>>vo=12.5*mag*cos(w*t+phase*pi/180);>>plot(t,vg,t,vo);gridaxis([22.002-5050]);>>text(0.25,0.85,'OutputVoltage','sc');>>text(0.07,0.37,'InputVoltage','sc');>>title('稳态滤波输出');>>ylabel('电压(v)'),xlabel('时间(s)');3）.拉氏变换法求解>>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([1,1+2e-3,-50,50]);4）.将传递函数展开为部分分式>>clear;>>num=[1e11];den=[1,2.5e6,1e12,0];>>[r,p,k]=residue(num,den)r=0.033333333333333-0.1333333333333330.100000000000000p=-2000000-5000000k=02.程序设计实验2）.已知系统传递函数，确定其零极点，画出零极点分布图，确定其阶跃响应。>>clearall>>b=[2,0,1,-3,1,4];>>a=[5,2,-1,-3,5,2,-4,2,-1];>>zs=roots(b);>>ps=roots(a);>>subplot(2,1,1);>>plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',12);>>grid;>>title('零极点分布图');>>legend('零点','极点');>>sys=tf(b,a);>>t=0:0.1:10;>>y=step(sys,t);>>subplot(2,1,2);>>plot(t,y);>>grid;>>xlabel('时间（t）');>>ylabel('y(t)');>>title('单位阶跃响应');3）.已知系统传递函数，确定其零极点，画出零极点分布图，确定其冲激响应。>>clearall>>b=[1.65,-0.33,-576,90.6,19080];>>a=[1,0.996,463,97.8,12131,8.11,0];>>zs=roots(b);>>ps=roots(a);>>subplot(2,1,1);>>plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',12);>>grid;>>title('零极点分布图');>>legend('零点','极点');>>sys=tf(b,a);>>t=0:0.1:10;>>y=step(sys,t);>>