|
t=0:0.1:12:num=[1];
9 N& }8 Z$ G2 FZeta1=0:den1=[1 2*Zeta1 1];% J7 p L3 U f( e; D4 L
Zeta2=0.3:den2=[1 2*Zeta2 1];1 V% c3 w: H6 {) `+ P- e! J$ @4 r/ S
Zeta3=0.5:den3=[1 2*Zeta3 1];
0 z/ ~% ?; G/ Y! W: h" WZeta4=0.7:den4=[1 2*Zeta4 1];; ]- H9 [/ Y& Y* |
Zeta5=1:den5=[1 2*Zeta5 1];
$ t; E5 \. r& V0 l+ Y: ~: v$ E/ l[y1,x,t]=step(num,den1,t);: W; P; C7 o: O8 A
[y2,x,t]=step(num,den2,t);% R. n/ [; g* A8 t( H. X8 \2 U2 F9 G0 c6 W
[y3,x,t]=step(num,den3,t);
4 k% T3 f9 b+ n$ L1 b6 J[y4,x,t]=step(num,den4,t);
* s' [ \1 U4 [' m[y5,x,t]=step(num,den5,t);" w9 m! G) j* G! F, w2 W( q
plot(t,y1,t,y2,t,y3,t,y4,t,y5);3 V; j8 L& f& @: f2 @
grid on; |
|