|
t=0:0.1:12:num=[1];
" P9 i2 K9 i# X* EZeta1=0:den1=[1 2*Zeta1 1];" j1 P0 G0 B6 K
Zeta2=0.3:den2=[1 2*Zeta2 1];
" [" n1 E! U2 \0 J3 U' e& |; GZeta3=0.5:den3=[1 2*Zeta3 1];; |8 k6 s; Z% N3 h( c
Zeta4=0.7:den4=[1 2*Zeta4 1];( f0 H7 L" `! N) \0 l; @9 Z) _9 J
Zeta5=1:den5=[1 2*Zeta5 1];4 M# ]# H7 @3 }: J& ~; T7 d+ |4 o
[y1,x,t]=step(num,den1,t);
4 p5 ^1 K" Q* U4 `/ l[y2,x,t]=step(num,den2,t);/ G+ Q m; H4 q7 N! Q. w% K7 n3 h
[y3,x,t]=step(num,den3,t);
8 z& n5 v; m7 R[y4,x,t]=step(num,den4,t);: \& U% b' W) b6 \+ ^! {- M
[y5,x,t]=step(num,den5,t);1 i( `8 @- U( q6 U
plot(t,y1,t,y2,t,y3,t,y4,t,y5);
4 S$ m. L, G. l0 f7 Ygrid on; |
|