|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
; { k# Q8 N2 oMatlab求解系统的响应2 a7 Z5 o) t" i" Q' `- D
求LTI系统的零状态响应的函数lsim，其调用格式为4 r% p4 Q' C% d: r* }7 h/ j% K/ j( i
/ K0 T4 }( i/ c式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
1 t& H2 I; X4 \5 o. D5 Y5 U; J
0 ]5 c! q, R! V系统模型sys要借助tf函数获得，其调用方式为2 h, I! ?' L; S) j2 B% _( s& n
5 r, n3 z' c7 d9 F
- sys=tf(b, a)
# U( u' [. R7 O& n% `- }& ?& r
; V: {- R% |# c, Q
) t( h" x% g5 n' @1 V式中，b和a分别为微分方程的右端和左端各项的系数。
6 @# }6 j9 B2 g- h) F比如：$ i8 y7 C5 ^( L. y3 t2 b- Q0 m1 a9 `
7 P; h0 Q1 m1 E- O( i+ D9 O7 X7 S
0 U5 l; O- f* K7 y; Z0 }0 g/ S4 o2 b1 n6 g; z" X
$ D3 A5 H2 r% B/ f0 U. T( u3 n: i! u3 e) T: q* ?
: m( U* n6 X& B' C
- 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)')
. y1 [: X( [% s- w9 |$ f, N7 B$ O
. X7 n6 E5 ~" s8 Q e/ j, v0 k L# r+ ^- Q' Z7 l9 m
+ g' H0 l' J0 d4 u+ k' y* [: k2 \; F3 X6 O+ O* M9 T# Q
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！