资源描述:
《东南大学信号与系统MATLAB实践第二次作业》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、练习二实验六一.用MATLAB语言描述下列系统,并求出极零点、1.>>Ns=[1];Ds=[1,1];sys1=tf(Ns,Ds)实验结果:sys1=1-----s+1>>[z,p,k]=tf2zp([1],[1,1])z=Emptymatrix:0-by-1p=-1k=12.>>Ns=[10]Ds=[1,-5,0]sys2=tf(Ns,Ds)实验结果:Ns=10Ds=1-50sys2=10---------s^2-5s>>[z,p,k]=tf2zp([10],[1,-5,0])z=Emptymatrix:0-by-1p=05k=10二.已知系统的系统函数如下,用MATLAB描述
2、下列系统。1.>>z=[0];p=[-1,-4];k=1;sys1=zpk(z,p,k)实验结果:sys1=s-----------(s+1)(s+4)Continuous-timezero/pole/gainmodel.2.>>Ns=[1,1]Ds=[1,0,-1]sys2=tf(Ns,Ds)实验结果:Ns=11Ds=10-1sys2=s+1-------s^2-1Continuous-timetransferfunction.3.>>Ns=[1,6,6,0];Ds=[1,6,8];sys3=tf(Ns,Ds)实验结果:Ns=1660Ds=168sys3=s^3+6s^2+6s
3、-----------------s^2+6s+8Continuous-timetransferfunction.六.已知下列H(s)或H(z),请分别画出其直角坐标系下的频率特性曲线。1.>>clear;forn=1:400w(n)=(n-1)*0.05;H(n)=(1j*w(n))/(1j*w(n)+1);endmag=abs(H);phase=angle(H);subplot(2,1,1)plot(w,mag);title('幅频特性')subplot(2,1,2)plot(w,phase);title('相频特性')实验结果:2.>>clear;forn=1:400w(n
4、)=(n-1)*0.05;H(n)=(2*j*w(n))/((1j*w(n))^2+sqrt(2)*j*w(n)+1);endmag=abs(H);phase=angle(H);subplot(2,1,1)plot(w,mag);title('幅频特性')subplot(2,1,2)plot(w,phase);title('相频特性')实验结果:3.>>clear;forn=1:400w(n)=(n-1)*0.05;H(n)=(1j*w(n)+1)^2/((1j*w(n))^2+0.61);endmag=abs(H);phase=angle(H);subplot(2,1,1)pl
5、ot(w,mag);title('幅频特性')subplot(2,1,2)plot(w,phase);title('相频特性')实验结果:4.>>clear;forn=1:400w(n)=(n-1)*0.05;H(n)=3*(1j*w(n)-1)*(1j*w(n)-2)/(1j*w(n)+1)*(1j*w(n)+2);endmag=abs(H);phase=angle(H);subplot(2,1,1)plot(w,mag);title('幅频特性')subplot(2,1,2)plot(w,phase);title('相频特性')实验结果:实验七三.已知下列传递函数H(s)或H
6、(z),求其极零点,并画出极零图。1.>>z=[1,2]';p=[-1,-2]';zplane(z,p)实验结果:2.>>z=[1,2];p=[-1,-2];zplane(z,p)>>num=[1];den=[1,0];[z,p,k]=tf2zp(num,den);zplane(z,p)>>num=[1];den=[1,0];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果:z=Emptymatrix:0-by-1p=0k=13.>>num=[1,0,1];den=[1,2,5];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果
7、:z=0+1.0000i0-1.0000ip=-1.0000+2.0000i-1.0000-2.0000ik=14.>>num=[1.8,1.2,1.2,3];den=[1,3,2,1];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果:z=-1.22840.2809+1.1304i0.2809-1.1304ip=-2.3247-0.3376+0.5623i-0.3376-0.5623ik=1.80005.>>clear;A=[0,1,0;0,0,1;-