|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应0 ^2 `- `& s; y
* |0 B" a" M* _# n) y7 G' {# g
求LTI系统的零状态响应的函数lsim，其调用格式为2 l$ g: ]) R# T. P( u6 z8 y/ z; d
9 Y, A3 [1 z3 J c式中，t表示计算系统响应的抽样点向量；f是系统输入( R7 [ X8 D* J! x( K" r4 B
信号，sys是LTI系统模型，用来表示微分方程。
6 c5 X# a0 I* k' L7 { o& x/ K. o, N: a* t5 i- G/ q
系统模型sys要借助tf函数获得，其调用方式为- c$ J( u* d0 h9 X7 G8 s* Y
k, V& L8 E2 W9 V& l
sys=tf(b, a)5 D! B/ c B3 ]: B, j3 F9 M; N
# O$ _ x9 [4 ~# V: n
/ V, t# @( K( C! I: j0 ?式中，b和a分别为微分方程的右端和左端各项的系数。
3 a( h0 v% a) a. C8 a" u2 p; [2 D3 H比如：0 O, \0 Q1 t1 w% _. a
9 N7 j* U: D& g0 B' n
9 N/ N9 Y8 Y! `: c" U0 `. X
$ N( u( r; k6 ^1 M; k: |
" Q" d7 |9 x5 x, V# }3 _. C1 E9 D& z2 ], M5 I U
( A7 o5 V, a n9 k" D2 k
sys=tf([1],[1,2,77]); % tf函数获得系统模型sys
& |3 ]1 {+ p0 U& o2 Zt=0:0.01:5; % 对时间t进行离散取样
& Y* D6 B: a3 O0 y$ Af=10*sin(2*pi*t);
2 }* h, F1 O5 K, Vy=lsim(sys,f,t); % 求系统的零状态响应
$ Y3 I% y; F) S: V+ Xplot(t,y)* r( j( e7 V$ I) p
xlabel('Time(sec)')/ f2 C! L: c% c# E7 G5 ~: x, g
ylabel('y(t)'). y) J5 e, }% {! x# a
5 L# h/ M+ {1 n \9 h
: b6 U: i. F9 A4 y8 r @0 \4 [
0 a8 C1 m, F$ z0 |# g8 k2 o O# ~. K$ ]! y& @
|
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！