|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
6 m4 L# t' g# Y: c) y3 M, qMatlab求解系统的响应
& g' _8 ?6 H C V2 s6 [求LTI系统的零状态响应的函数lsim,其调用格式为9 W8 W5 O: n* k9 R! a" U5 N; O4 F% x% c
' V5 U1 \. P4 N9 q N1 Q3 k( I
式中,t表示计算系统响应的抽样点向量;f是系统输入信号,sys是LTI系统模型,用来表示微分方程。
% a1 ~+ _! y) r6 C# i8 O& I! v, o$ B8 l. l9 h& Z; e' ~8 O
系统模型sys要借助tf函数获得,其调用方式为( @! o* z4 q. o" {2 g0 ~% E
5 V; ? u) p& q1 F8 R8 ^2 R+ T
- sys=tf(b, a)
* K! i Y6 Y. d) i- b9 x* s " a6 H0 F2 _5 }/ L1 ^
1 B3 R3 ^3 l9 V) p3 W. P式中,b和a分别为微分方程的右端和左端各项的系数。9 ~! g1 j/ p8 b( d
比如:% j8 |9 e! Y" _5 N3 f' G' m
% ^& U1 l: w% @$ _! C- s
5 G# x! E/ f( E! L6 ~- Y
0 A# e+ w0 ^9 L% L. w4 f6 P6 N
2 U+ y" w- G `; X
2 j0 q- w% f+ L) I
4 N& d. h1 Z* d% O/ F
- 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)')
; @, k1 Z$ A% a0 a ; |9 @( G' J9 @% k# p# ]3 g! q% a
?# o8 [/ a& X: S& d
- y {2 ]6 j' k6 g6 e
! j7 I$ v' n/ y |
|