资源描述:
《东南大学4系短学期matlab第二次作业》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验六一.用matlab语言描述下列系统,并求出极零点。>>Ns=[1];Ds=[1,1];sys1=tf(Ns,Ds)Transferfunction:1-----s+1>>pole(sys1)ans=-1>>zero(sys1)ans=Emptymatrix:0-by-1>>Ns=[10];>>Ds=[1,-5,0];>>sys1=tf(Ns,Ds)Transferfunction:10---------s^2-5spole(sys1)ans=05>>zero(sys1)ans=Emptymatrix:0-by-1一.已知系统的系统函数如下,用matlab
2、描述下列系统。>>Ns=[1,0];>>Ds=[1,5,4];>>sys1=tf(Ns,Ds)Transferfunction:s-------------s^2+5s+4>>Ns=[1,1];Ds=[1,0,-1];>>sys1=tf(Ns,Ds)Transferfunction:s+1-------s^2–1>>Ns=[1,6,6,0];Ds=[1,6,8];>>sys1=tf(Ns,Ds)Transferfunction:s^3+6s^2+6s-----------------s^2+6s+8六.画频率特性曲线和波特图。1.>>num=[1];>>den
3、=[11];>>freqz(num,den,128)>>title('频率特性')>>dbode(num,den,1);>>title('波特图')2.>>a=[1,2^(0.5),1];b=[2,0];>>sys=tf(b,a);freqs(b,a,128)>>bode(sys)3.>>num=[121];>>den=[100.61];>>freqz(num,den,128);>>title('频率特性')>>dbode(num,den,1);>>title('波特图')4.>>a=[1,3,2];b=[3,-9,6];sys=tf(b,a);freqs(b
4、,a,128)bode(sys)实验七一.求极零点,并画极零图。1.>>num=[3,-9,6];>>den=[1,3,2];>>[z,p,k]=tf2zp(num,den);>>zplane(z,p)2.>>num=[1];>>den=[1,0];>>[z,p,k]=tf2zp(num,den);>>zplane(z,p)3.>>num=[1,0,1];>>den=[1,2,5];>>[z,p,k]=tf2zp(num,den);>>zplane(z,p)4.>>num=[1.8,1.2,1.2,3];>>den=[1,3,2,1];>>[z,p,k]=tf
5、2zp(num,den);>>zplane(z,p)5.>>A=[0,1,0;0,0,1;-6,-11,-6];>>B=[0;0;1];>>C=[4,5,1];>>D=0;>>sys=ss(A,B,C,D);>>pzmap(sys)五.1.>>A=[5,2,1,0;0,4,6,0;0,-3,-6,-1;1,-2,-1,3];B=[1;2;3;4];>>A=[5,2,1,0;0,4,6,0;0,-3,-6,-1;1,-2,-1,3];>>B=[1;2;3;4];>>C=[1,2,5,2];>>D=0;>>sys=ss(A,B,C,D);>>pzmap(sys)4
6、.>>z=[-3];>>p=[-1,-5,-15];>>k=15;>>sys2=zpk(z,p,k)Zero/pole/gain:15(s+3)------------------(s+1)(s+5)(s+15)>>zplane(z,p)5.>>a=[1,3,5,2];>>r=roots(a)r=-1.2267+1.4677i-1.2267-1.4677i-0.5466七.作奈奎斯特图。1.>>b=[1];>>a=[1,3,2];>>sys=tf(b,a);>>nyquist(sys)2.>>b=[1];>>a=[1,4,4,0];>>sys=tf(b,a);
7、>>nyquist(sys)3.>>b=[1];>>a=[1,2,2];>>sys=tf(b,a);>>nyquist(sys)