东大matlab实践课三.docx

1、东大matlab实践课三实验六一2b=10 10;a=1 0;z,p,k=tf2zp(b,a)z = -1p = 0k =103b=10;a=1 -5;z,p,k=tf2zp(b,a)z = Empty matrix: 0-by-1p = 5k =10二1 z=0; p=-1 -4; k=1; sys1=zpk(z,p,k) Zero/pole/gain: s-(s+1) (s+4)2 Ns=1,1; Ds=1,0,-1; sys2=tf(Ns,Ds) Transfer function: s + 1-s2 13 Ns=1,6,6,0; Ds=1,6,8; sys3=tf(Ns,Ds) Tra

2、nsfer function:s3 + 6 s2 + 6 s- s2 + 6 s + 8实验七三1 z=1,2; p=-1,-2; zplane(z,p)2num=1;den=1,0;z,p,k=tf2zp(num,den)zplane(z,p)z = Empty matrix: 0-by-1p = 0k = 13num=1,0,1;den=1,2,5;z,p,k=tf2zp(num,den)zplane(z,p)z = 0 + 1.0000i 0 - 1.0000ip = -1.0000 + 2.0000i -1.0000 - 2.0000ik = 14num=1.8,1.2,1.2,3;d

3、en=1,3,2,1;z,p,k=tf2zp(num,den)zplane(z,p)z = -1.2284 0.2809 + 1.1304i 0.2809 - 1.1304ip = -2.3247 -0.3376 + 0.5623i -0.3376 - 0.5623ik =1.8000五5 b=100*conv(1,0,conv(1,2,conv(1,2,conv(1,3,2,1,3,2); a=conv(1,-1,conv(1,1,conv(1,3,5,2,conv(1,0,2,0,4,1,0,2,0,4); z,p,k=tf2zp(b,a)zplane(z,p)z = 0 -2.0004

4、+ 0.0004i -2.0004 - 0.0004i -1.9996 + 0.0004i -1.9996 - 0.0004i -1.0000 + 0.0000i -1.0000 - 0.0000ip = 1.0000 0.7071 + 1.2247i 0.7071 - 1.2247i 0.7071 + 1.2247i 0.7071 - 1.2247i -1.2267 + 1.4677i -1.2267 - 1.4677i -0.7071 + 1.2247i -0.7071 - 1.2247i -0.7071 + 1.2247i -0.7071 - 1.2247i -1.0000 -0.546

5、6 k = 100不稳定6 b=conv(1,-1.414,1,1,1); a=conv(1,0.9,0.81,1,-0.3); z,p,k=tf2zp(b,a)zplane(z,p)z = -1.0000 0.7070 + 0.7072i 0.7070 - 0.7072ip = -0.4500 + 0.7794i -0.4500 - 0.7794i 0.3000 k = 1稳定七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); nyquist(sys);3 b=1; a=1,2,2; sys=

6、tf(b,a); nyquist(sys);实验八一1 b=1,3; a=1,3,2; sys=tf(b,a); impulse(sys); B=1,3; A=1,3,2; a,b,c,d=tf2ss(B,A)a = -3 -2 1 0b = 1 0c = 1 3d = 0 sys=ss(a,b,c,d); x0=4,-1; initial(sys,x0); grid;2 b=1,3; a=1,2,2; sys=tf(b,a); impulse(sys); B=1,3; A=1,2,2; a,b,c,d=tf2ss(B,A); sys=ss(a,b,c,d); x0=1.6,-0.2; ini

7、tial(sys,x0); grid;3 b=1,3; a=1,2,1; sys=tf(b,a); impulse(sys); B=1,3; A=1,2,1; a,b,c,d=tf2ss(B,A); sys=ss(a,b,c,d); x0=7/4,-1/4; initial(sys,x0); grid;二b=1; a=1,1,1,0;sys=tf(b,a);subplot(1,2,1);impulse(sys); subplot(1,2,2);step(sys);subplot(1,2,1); impulse(sys); subplot(1,2,2); step(sys); t=0:0.08:10 e=sin(t);sys=tf(b,a);lsim(sys,e,t);lsim(sys,e,t);三1 clear; b=1,3; a=1,3,2; t=0:0.08:10; e=exp(-3.*t); sys=tf(b,a); lsim(sys,e,t);2b=1,3;a=1,2,2;t=0:0.08:10;e=exp(-0.*t);sys=tf(b,a);lsim(sys,e,t);3b=1,3;a=1,2,1;t=0:0.08:10;e=exp(-2.*t);sys=tf(b,a);lsim(sys,e,t);