|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应
5 m" h' d$ H# S
- S6 M. S" K5 Y7 U; u4 g. b; q5 k求LTI系统的零状态响应的函数lsim，其调用格式为3 n: f- }7 ]. D% g" P2 Y; O& K
5 _+ l( q' f; F" U. q6 |: {
式中，t表示计算系统响应的抽样点向量；f是系统输入
& }: ?% V7 N- p% g8 t信号，sys是LTI系统模型，用来表示微分方程。( Q/ Z" k o; j5 \
% m5 q% K: K) H6 v5 L0 L3 V$ B+ O系统模型sys要借助tf函数获得，其调用方式为
% B1 u$ B4 y2 F8 H. S' k4 I6 V8 C: B/ h6 G; R/ S, ]0 A
sys=tf(b, a)4 J, v- t1 O2 d6 _' |
) y- b- J4 N& Q* X3 I$ _! n& f8 Y1 ~
t/ J# T' @9 Y; r' t6 l式中，b和a分别为微分方程的右端和左端各项的系数。
" S$ d7 P9 J; w/ z# L k比如：0 ~' u! V& j N
# Z5 y( J* |7 V. I) Z/ L
E& K/ X" y- G% A" h9 c! w) z5 _% \& V3 i8 @
7 Q% \7 o4 Y$ s! ^# T, S+ s! P$ e7 s& d
+ i' I7 w( [; j* a8 ~sys=tf([1],[1,2,77]); % tf函数获得系统模型sys
: t5 o9 R X3 @$ @: O, x+ {% ot=0:0.01:5; % 对时间t进行离散取样5 S* B( W, H) p! k8 C! i5 g: t# y, V
f=10*sin(2*pi*t);
, K1 K$ [9 g' d! u2 C- V8 l* O Ky=lsim(sys,f,t); % 求系统的零状态响应* G. D M E; {" ^: m
plot(t,y)
$ W; J; P$ y; I5 J. m' l0 Xxlabel('Time(sec)'): C h' {2 ^" X/ e( ]3 m
ylabel('y(t)')$ k7 p/ t( ?) Q2 C" Q" Z. L5 o# @
* u) H5 L- V4 ]9 {2 m6 _$ {
1 A" u) p7 i* N b3 a/ i- |" I% G; Y3 c0 `* E
- G8 v5 U n. U# r |
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！