|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
. p l: H O3 b
Matlab求解系统的响应8 V% z) M/ B0 c* d5 C. y+ b1 ~
求LTI系统的零状态响应的函数lsim，其调用格式为. ]" e# k0 G1 Q- y4 L
2 z) \5 X2 p; r) |( P式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
# M v/ w3 O0 a/ @/ {4 Y& u" k
" L) c3 I* u# M. T% h( D8 d系统模型sys要借助tf函数获得，其调用方式为
& [2 r% n. E3 t2 w# d3 E0 `8 L" n2 n9 ?! T
- sys=tf(b, a)
2 d0 t# k; S$ Y, \
' K7 w: T: d5 X1 d, l. X. }* X0 C8 |( y" t
式中，b和a分别为微分方程的右端和左端各项的系数。
3 }" ]! c6 x6 Y7 Q比如：
* n- v: v$ r4 @) E' E7 e7 a' [) m; V0 z; E
' F: N+ b6 n0 O6 g2 ?* T
/ I1 t, o* i+ R
; R/ J: \: c( t# }
) d& l" j5 I" w8 k1 \& b9 J. J) u: t$ |+ a9 @6 d% Q9 g8 P: E4 u. Z
- 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)')
9 B7 a) m' |* m- [' S3 P5 Y5 v
9 F2 ]* Q6 D7 I
* l4 r7 s o5 H% U
9 s* _; v. b: b" o+ D
5 Y* @ C/ ]5 _2 ~$ l& K3 j
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！