五．實驗結果（程序設計實驗源程序和結果）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);>>