|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
" H6 G( ]" Y3 C; n
Matlab求解系统的响应
8 U1 Y6 z0 d% N Y. b求LTI系统的零状态响应的函数lsim，其调用格式为$ a/ e! h% @2 k2 |' J3 f, n2 O
4 H) ~* t7 w8 D9 {. }, v
式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。/ G2 k: b( x$ W* p$ W
; s; w+ ]& r7 d
系统模型sys要借助tf函数获得，其调用方式为
9 t) ?+ e3 S4 Q( @2 u3 f
( t2 u9 ~! [ _$ M# @- sys=tf(b, a)
. e# V* F- p* S6 A9 \5 W
# D. m/ R- ?3 }6 f8 G/ l9 y
5 T* Z; J4 J1 f4 \- b" Y式中，b和a分别为微分方程的右端和左端各项的系数。, p0 ^' Q- W/ ]. d5 t/ ?
比如：
0 z4 r5 R1 `6 C8 ~. X. ^: O4 L) N7 k+ ]
+ `7 P- `* m5 e5 @& z
. {( B# y5 ?$ G" \3 ^ h
6 V, J$ F% s- F" B) I# q4 ^3 S6 u1 {1 ^1 ~" S" c
; G/ m7 Z" {% w( ?
- 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)')# x* y# z$ w9 B, w7 C( y8 V2 `
`+ M9 y( Y2 i
# n# n' H4 j+ }7 F2 r
6 a0 L5 U8 K) e+ V$ r: x/ S" ~+ ~! a, U1 \, M# o3 v
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！