东南大学4系短学期matlab第二次作业.doc

实验六一 用matlab语言描述下列系统，并求出极零点。 Ns=1;Ds=1,1;sys1=tf(Ns,Ds) Transfer function: 1-s + 1 pole(sys1)ans = -1 zero(sys1)ans = Empty matrix: 0-by-1 Ns=10; Ds=1,-5,0; sys1=tf(Ns,Ds) Transfer function: 10-s2 - 5 spole(sys1)ans = 0 5 zero(sys1)ans = Empty matrix: 0-by-1二 已知系统的系统函数如下，用matlab描述下列系统。 Ns=1,0; Ds=1,5,4; sys1=tf(Ns,Ds) Transfer function: s-s2 + 5 s + 4 Ns=1,1;Ds=1,0,-1; sys1=tf(Ns,Ds) Transfer function: s + 1-s2 1 Ns=1,6,6,0;Ds=1,6,8; sys1=tf(Ns,Ds) Transfer function:s3 + 6 s2 + 6 s- s2 + 6 s + 8 六画频率特性曲线和波特图。1. num=1; den=1 1; 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=1 2 1; den=1 0 0.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,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=tf2zp(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. 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;