|
t=0:0.1:12:num=[1];
/ _% u, U0 I; e( B0 d* SZeta1=0:den1=[1 2*Zeta1 1];) ^4 v8 w0 g- N' g
Zeta2=0.3:den2=[1 2*Zeta2 1];; G; n. @3 U% D) _1 @
Zeta3=0.5:den3=[1 2*Zeta3 1];
! Z! O3 F! e o9 ?$ r0 R( }Zeta4=0.7:den4=[1 2*Zeta4 1];/ Y- T$ {* C# l5 j5 `6 D
Zeta5=1:den5=[1 2*Zeta5 1];$ P& r+ h0 E, p
[y1,x,t]=step(num,den1,t);9 j, \- ~+ ~: h8 @4 Q! |
[y2,x,t]=step(num,den2,t);
: N8 e; B5 y5 v9 b0 e1 Y[y3,x,t]=step(num,den3,t);
- f! V' E n8 |! p[y4,x,t]=step(num,den4,t);
1 T7 N _6 d+ }& Y[y5,x,t]=step(num,den5,t);( N8 T' D' D% ], f& ^& g
plot(t,y1,t,y2,t,y3,t,y4,t,y5);% F$ a8 @0 f$ P4 e
grid on; |
|