|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应
# u3 M( H K! t' V, [
) ^* A. v* v: b/ ?6 B: h& l求LTI系统的零状态响应的函数lsim，其调用格式为5 F5 {* k+ r* Y8 ]
/ D" r4 i7 J( g, O* k3 N! |4 E* R( B式中，t表示计算系统响应的抽样点向量；f是系统输入- }2 W z2 h. m3 W1 h
信号，sys是LTI系统模型，用来表示微分方程。
1 e4 X( S# b) _* t0 h/ t, N* n# Q1 f' J* r4 X6 c
系统模型sys要借助tf函数获得，其调用方式为
: l$ T+ r/ f: y4 m! Z/ ~7 _) P2 r& o; g/ \
sys=tf(b, a)
% J4 |7 o$ L8 F) U% E
( B/ D1 W$ h$ S' Z6 o( I- L+ g& G
* F% X' { F7 h式中，b和a分别为微分方程的右端和左端各项的系数。
0 I0 w. k3 H- l2 T- I( Q! l4 b2 Z2 ?比如：
. p+ R' ~; l/ T7 [# N& F! r: l2 e$ T- f
* T5 r' Z1 Z* H7 Z. k% l, T
2 p2 Q0 [0 t* [* R2 O
- r+ u7 H8 G1 J: S) p( Q- O0 u( Z8 J) s+ ]' y
2 I/ \, N& w' [; G0 Ssys=tf([1],[1,2,77]); % tf函数获得系统模型sys
5 A. R+ J; h# q. w( ~t=0:0.01:5; % 对时间t进行离散取样
S5 `2 i K. p6 s jf=10*sin(2*pi*t); 4 c& r! l% H& T
y=lsim(sys,f,t); % 求系统的零状态响应
9 f$ v3 Y2 g" N+ C9 J8 ^$ T# k" Dplot(t,y)/ U: @7 C' n+ }: g% V, T( g+ X
xlabel('Time(sec)')
2 P0 T. s: Z" b& P! |$ Q uylabel('y(t)')
7 j' w5 y! A& N
1 |. E4 R8 W; W- t% c. q% r2 }
k5 d; o0 b8 {/ u; ~* }
! x$ I$ d5 Q9 e* U0 e' R+ u) F6 ?
) Y+ q/ l4 c. ]: g# K0 }) @2 B |
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！