|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
* Z6 r- f, D+ n" E+ l0 G( s- q
Matlab求解系统的响应4 a4 {% w. E. J7 E$ T- ~
求LTI系统的零状态响应的函数lsim，其调用格式为
" y( k" ~+ ~( O, o/ T8 |0 t0 K3 w9 J! ^" E8 }7 o* J
式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
& |+ r; j4 D7 @" ^+ y- ], x8 S0 k n/ T2 O; e
系统模型sys要借助tf函数获得，其调用方式为
/ v6 Y- \# s0 V5 T2 n
* d5 L1 m& ~4 {4 G. y- sys=tf(b, a)! \+ A) @" @4 J/ I5 `& ]
7 y( k) H5 y H6 Q+ [2 i$ ~ `' H( a6 X. t, \2 L9 _5 R
式中，b和a分别为微分方程的右端和左端各项的系数。1 Q8 w3 _) b$ A9 k8 ]
比如：3 Z+ _! a1 G; C' h, t
8 T8 B% V0 W/ Z0 O
. ~3 u2 j' m3 ~) M, j
4 x6 S) G7 c8 e' N' }% y1 L: c
+ l; j3 q3 ?/ s$ H1 R6 c
B) y: `, j8 u- @# x9 _ g1 V( o8 c& N( }' Z# F2 i: X' W7 Q6 t! k! Q
- sys=tf([1],[1,2,77]); % tf函数获得系统模型sys
- t=0:0.01:5; % 对时间t进行离散取样
- f=10*sin(2*pi*t);
- y=lsim(sys,f,t); % 求系统的零状态响应
- plot(t,y)
- xlabel('Time(sec)')
- ylabel('y(t)')
- J1 j# r$ r1 m+ P& T8 `( b
+ E9 S& {! L# u) a; d
- V' @1 C) Y5 N9 |/ C
4 V8 W% o# I# ^. n' a7 ^
# w9 K1 ^' T! T
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！