|
t=0:0.1:12:num=[1];
( J f9 B+ D5 C) sZeta1=0:den1=[1 2*Zeta1 1];7 ^* B- r/ b' X' j# W! q4 R& ]
Zeta2=0.3:den2=[1 2*Zeta2 1];
; R8 O, P$ u. m5 ?: \4 t3 ]Zeta3=0.5:den3=[1 2*Zeta3 1];1 j; B o1 l! ]
Zeta4=0.7:den4=[1 2*Zeta4 1];: c4 |+ m, i* H$ L# i& f
Zeta5=1:den5=[1 2*Zeta5 1];% b4 R: l2 @: Z9 U* X$ X0 j
[y1,x,t]=step(num,den1,t);8 t" {" w: Q2 o( J3 ?- k/ V
[y2,x,t]=step(num,den2,t);# a4 F: C2 Y. d1 m$ H7 ^
[y3,x,t]=step(num,den3,t);/ S0 ~5 b- l5 s( V5 t6 i7 O8 \' B
[y4,x,t]=step(num,den4,t);
/ H) V8 I \8 L) r& j/ {, k7 O[y5,x,t]=step(num,den5,t);
6 r- G8 Z5 I' X; b: Xplot(t,y1,t,y2,t,y3,t,y4,t,y5);
& O8 `3 ]% u$ P5 A0 x, [grid on; |
|