|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
7 d2 b4 v% S1 ?8 ]3 e2 z8 U3 ?' ?
Matlab求解系统的响应) V) `: |$ ?! K4 E
求LTI系统的零状态响应的函数lsim，其调用格式为% H/ W2 u* P) G- g9 O
4 |: [ C7 n( U式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。6 E$ T+ V3 u# k
5 Z1 P3 ~7 V4 M8 ?系统模型sys要借助tf函数获得，其调用方式为
6 M2 u: [! |, |+ o5 O. {
, c4 ^$ S {6 a! |% M* D- sys=tf(b, a)
& G) c# h( q w1 D, F7 Y% N
: p9 _3 |- [( u4 o- ^# y# u6 }: g: ^! }9 T5 ]# q9 h
式中，b和a分别为微分方程的右端和左端各项的系数。3 z3 Y4 G& c) L6 O) V( L0 q
比如：5 Y1 h* Q- h2 }; u. H
* ?' o: a8 ~' S5 U9 a3 d
( ?! b# Y5 ~& E) S% f* u# Z$ i
) D9 D' x. T: C0 E4 _
, i7 m* | d' B
% b) M4 g& {8 v& a8 c: S" H7 k. U# C; Q8 o" c+ Y
- 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)')
- h' q3 S p+ s8 m! X8 I: S
7 ~) H* s, ]# ~2 d; e% ~
( i/ U0 [: ?* X. _
8 @2 C8 Z) n. \/ \+ A& R' l; @ D- s4 p% y+ y' a
|
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！