|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
Matlab求解系统的响应
9 t7 \+ B* {' m# V; O6 g) {% E2 V G. c3 B9 D& b" h1 m/ n
求LTI系统的零状态响应的函数lsim,其调用格式为
* L9 u" T' C5 }( I6 w0 ^9 d( K% Y# o, t v5 s, R* o! V$ [; Y
式中,t表示计算系统响应的抽样点向量;f是系统输入
( h: ]$ ?) `% V/ i8 C4 r信号,sys是LTI系统模型,用来表示微分方程。" j" b! v+ K/ y9 d6 m4 e* M
6 Y5 U9 z1 v9 d系统模型sys要借助tf函数获得,其调用方式为; x$ j: {- {' {& e
7 O* h, {$ F; c. z( v7 G! y: \: a
sys=tf(b, a)1 l! r8 q. w' a; ]) ~; I2 G
7 j* L7 \. q; ?$ g$ q
5 a$ ?: O4 K2 t& |5 H. _& L9 v式中,b和a分别为微分方程的右端和左端各项的系数。
6 |* C; o q9 q3 d# O& i2 J/ l比如:
) f' L5 w. g$ q( D0 h, i6 i) x/ q. |! ~3 D0 @
2 N: u3 \) D$ O0 R8 v
" w+ l# ?7 x/ `. M5 B7 q
0 w" l' o6 w, M7 B* u9 s2 W. k
3 d/ M/ Q: L) y+ P4 s6 L/ d: b7 c5 [
sys=tf([1],[1,2,77]); % tf函数获得系统模型sys
' w4 l ^' ~3 h6 y* P! l) yt=0:0.01:5; % 对时间t进行离散取样
" p g& Q+ w6 t8 a* Ef=10*sin(2*pi*t);
) q Q" J. C! M7 K) Xy=lsim(sys,f,t); % 求系统的零状态响应# n1 {. L* R! w
plot(t,y)
; G, O& k6 X, R( R0 _xlabel('Time(sec)')
" }, H1 s/ Z$ r! @' _ylabel('y(t)')9 m F9 u( N" ^
7 ^% k. g6 f. V2 l
) _/ ~, K, ?8 z+ B" m- K
# h8 S/ m. T$ h$ s& J( Q* h( q5 P o: m
|
|