|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应
. R2 A- h# H/ }% c7 P! ]+ D) M3 ~
求LTI系统的零状态响应的函数lsim，其调用格式为
" u" [) }- _' G7 i- X. t
. }3 l% Q+ W5 z/ y! W" w式中，t表示计算系统响应的抽样点向量；f是系统输入( _2 y" D$ w! r3 Q H3 [$ R
信号，sys是LTI系统模型，用来表示微分方程。$ q- b D; W- @0 n1 B. F5 b9 S/ O
9 ?7 |( u' Z9 C4 Y系统模型sys要借助tf函数获得，其调用方式为
! D2 x. I0 f3 z
2 |1 w3 y* e/ B/ C3 d! y& psys=tf(b, a)
; U7 i0 i, H5 n P) W) Z5 s+ r" i) l
' W6 l9 }/ w9 T6 d7 r; G1 f5 c( H: Y% D u7 R5 @! j) W# k
式中，b和a分别为微分方程的右端和左端各项的系数。
6 O2 N0 l/ y0 i. i比如：: G p8 F# D/ T H a$ f4 r4 f3 S; R
; x+ r$ W& e8 k; m' i8 r8 }( q
0 \5 i/ a; j$ z; n
! Z) n! A( z' H/ u K; j$ b
& Q! E! V# T& T! _0 W
; b/ B W& e$ v+ h( w) H/ `! _6 C5 {1 o( i$ b
sys=tf([1],[1,2,77]); % tf函数获得系统模型sys
, |7 d+ H2 y8 \# h% C' Tt=0:0.01:5; % 对时间t进行离散取样3 m5 [3 }. B5 f7 O& G! T6 j
f=10*sin(2*pi*t); 7 L8 [; [2 ~6 S7 O* M
y=lsim(sys,f,t); % 求系统的零状态响应
4 e L) w) I, Q- G! K% v2 x+ Q( tplot(t,y)5 @+ l; Y8 ~! ^+ U0 k; T( Z
xlabel('Time(sec)')
7 I2 ^- b, w; P6 C- Z$ Q# Pylabel('y(t)')
* z0 d9 h5 n4 a# e; P
P% C1 X Y- A) O" x2 [
5 d3 c2 P6 g; U
5 P6 J: u: J! _2 D* P- d# a
3 z: b1 [1 o: [- I1 ]- R |
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！