|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应
# d( @8 l3 `4 T% U, r9 ^. `! U
! `1 Q: j5 g- O, H! p* I. B求LTI系统的零状态响应的函数lsim，其调用格式为: g1 |0 ^# N- W
+ t }& Y+ ]! G3 v% e式中，t表示计算系统响应的抽样点向量；f是系统输入 _9 O. J. R8 d! U% I# r* R- A% |
信号，sys是LTI系统模型，用来表示微分方程。
[) U6 m0 }+ c! \! c9 f; f9 q$ P" N! |
系统模型sys要借助tf函数获得，其调用方式为4 w$ |: y+ E: R' X9 o; X
$ H' m: i8 I0 O" l% Q( a8 o0 _8 Dsys=tf(b, a)5 A3 y- o. m8 O. _$ U- P' z
" e6 Y" I. e7 ?: U2 S) w @ b; E4 m8 r! P8 ?7 }/ J
式中，b和a分别为微分方程的右端和左端各项的系数。
0 {4 q( R2 `6 `1 |, B3 u, ^ I, n) e2 D' F7 `比如：
, L; g+ I- `0 A1 c9 A; L( s( x2 t3 s0 I* U/ |
9 y& N+ H9 [1 M) S$ {5 b( J: E
[+ i/ r# T2 \7 z# l
. i* h1 A I; g0 V a' t1 v2 a# o7 x8 C
4 b! }9 J3 F( @' K6 V
' x; k2 _0 l& r$ Psys=tf([1],[1,2,77]); % tf函数获得系统模型sys
" n) m" c- A) B6 l1 f$ H, bt=0:0.01:5; % 对时间t进行离散取样
6 Y* ^& d) w* {& R7 o3 z/ B1 D, lf=10*sin(2*pi*t);
v1 y0 ]/ F2 Y( D v1 ]y=lsim(sys,f,t); % 求系统的零状态响应
- T" L5 Q2 i/ Hplot(t,y)# I5 v5 o3 O J' Y G
xlabel('Time(sec)')
% }: n) {2 i3 s9 \+ z( Lylabel('y(t)')0 M% n1 |+ f: D0 T( n/ z2 v
( ~7 y D' a8 V3 K+ M
( Q: l: P) H1 n! V& a" j! S/ u1 j0 I& X3 }
5 ~& Y% F+ V7 |) q$ p
|
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！