|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
; b( Z5 G7 H" z) V
Matlab求解系统的响应
! s7 \) O* _, R7 c- A% n4 R9 `求LTI系统的零状态响应的函数lsim,其调用格式为0 _7 q3 q- @, W& F e
2 a: @! ?; N9 S; x$ c
式中,t表示计算系统响应的抽样点向量;f是系统输入信号,sys是LTI系统模型,用来表示微分方程。; c* p0 v3 S. Q
+ G1 `) U2 P9 D1 e0 W, s系统模型sys要借助tf函数获得,其调用方式为
! C' O' A( Q; `8 t7 X$ v$ e: U
5 P7 E' g) j' J- M- sys=tf(b, a)8 a( y1 a' B9 z0 P3 ^' }
6 x7 X) `. E' U! T O; |
5 f" B% h* v2 U. q式中,b和a分别为微分方程的右端和左端各项的系数。) V# v$ Y8 i/ a) Y5 U
比如:
9 G3 W/ \* U5 Z0 R
, |8 ~; ^' f) {: C J5 Y2 F2 T
8 J- e/ H! A. {# w/ m/ b
7 J6 ^ F+ N9 C, F
0 B! y& e7 E( X c- l
! V) I8 X7 |& ]4 J/ F: R7 H/ a: c' w+ Z( D4 i7 a2 Z4 ]. b- X# 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)')2 P# D4 O) C% o6 h- [
* y8 ]$ V3 j* @% ^( ^5 N+ J+ u5 V: S0 x/ O- n# ]1 B
1 A; k8 ?, e2 y! |- Y
1 G3 C% p' ?* i6 P) Y" P/ \ |
|