|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
: F1 o$ t+ c2 m! Q+ E) W
Matlab求解系统的响应8 t3 W3 J: T- X/ g8 S
求LTI系统的零状态响应的函数lsim，其调用格式为
" V: T n$ ~7 }1 P. h! C, Y" A- C0 n5 L! u L- E
式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
7 R6 y; K1 F) J5 C* m5 N& H! W% n2 v$ o# L
系统模型sys要借助tf函数获得，其调用方式为
3 Z7 D6 w4 u) Y* P9 ]+ r" Q( }% {
3 I: u- \3 b& u4 Z( X% ?! \- sys=tf(b, a)4 m* v6 z, K6 g3 c8 I, b
" b% }& w; a9 q0 b# T
7 L. k; p7 E) m+ n式中，b和a分别为微分方程的右端和左端各项的系数。
4 _( R- ?9 ^ u) a比如：5 b: S$ K; N# b, M
. }& f& R: i3 U; n" |/ D; W
9 S% d& h# u9 d7 j
% n) J2 Y3 X6 O
* F2 H# O; f% B# V0 j: O! o
; V# c- `4 T+ a9 w! ~
8 r. M7 i; V/ d) n) W' E
- 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)')( s- T; x5 F: B$ x1 u
+ {: v; {2 j, X, N1 q8 C! V/ ? o
, z& l0 R( C+ f. c2 [
( K: f4 I# }; X! _/ ]( Z
; z9 h5 [. P d8 ?# d( z8 Q, i6 v0 {
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！