|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
' k% Q+ ? o: i% y; J
Matlab求解系统的响应) c! w/ e, j3 k) E
求LTI系统的零状态响应的函数lsim，其调用格式为
0 U, v+ I9 \ I6 n
* E. @, D4 |# n9 O* N3 s式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。6 r6 X' d0 M! \8 {4 K. m
! S4 R, X* n8 [9 P
系统模型sys要借助tf函数获得，其调用方式为7 B. C8 e X+ p
9 R$ B" r" I2 p, h5 o8 U
- sys=tf(b, a): c+ X2 ?6 f8 d" s3 P) E
0 e" F3 P* a1 [1 n; N, W, j$ X
% t- t2 W/ G) a' W: S式中，b和a分别为微分方程的右端和左端各项的系数。. _: z% o9 A# F0 E5 z, I4 H* t
比如：
$ g6 h- ?4 @+ {5 u. z6 y% S
* l$ v" G1 h$ ~ F9 \, a4 P7 T* t
" Y4 s+ u) d$ h0 F1 z
) Q7 ?9 L' @0 x4 U8 q3 m6 M+ j/ L
" r; j' q- H, N |. ]/ D3 s
. G2 q( s) D; z" H0 e% V2 ]" H& F2 x* ]* G+ A
- 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)')( Z6 Z* n2 @+ `( B8 D
2 I Z% @ ]; C+ z4 C; S
" d2 p3 w h& M5 D# m5 v$ F
7 o3 V! \: u5 m. T
5 ~* H+ Y5 e5 E+ r7 S |
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！