|
t=0:0.1:12:num=[1];( a$ F3 y. E8 S3 y* X+ D
Zeta1=0:den1=[1 2*Zeta1 1];' H( e2 j+ ~/ L* u. b
Zeta2=0.3:den2=[1 2*Zeta2 1];7 }' L n x# t* G r) E, w: ?* `
Zeta3=0.5:den3=[1 2*Zeta3 1];
# s( x& k9 g k2 ~) i& O7 r* OZeta4=0.7:den4=[1 2*Zeta4 1];2 m. a; n, b+ N/ K) ?
Zeta5=1:den5=[1 2*Zeta5 1];
2 O. l( j6 u- | q/ Y[y1,x,t]=step(num,den1,t);
( O: R$ f: L. {" E: r2 a4 x- k4 S* G[y2,x,t]=step(num,den2,t);2 E' O& E' u4 K/ G
[y3,x,t]=step(num,den3,t);
# Q! C$ x7 k- L L[y4,x,t]=step(num,den4,t);7 n7 O2 F. G |+ n! u) {% w% I7 X4 M
[y5,x,t]=step(num,den5,t);
( Y2 u# m U: w/ z; u. kplot(t,y1,t,y2,t,y3,t,y4,t,y5);
& C( `4 m, F7 c# D* h: V% ogrid on; |
|