|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
; I/ \& @3 [3 }7 }; jMatlab求解系统的响应
* J+ O! K' D& c1 p0 U% i8 ]. P! a求LTI系统的零状态响应的函数lsim,其调用格式为
) l0 `; q) _. \8 Q# i# H1 d3 [ D6 k% A0 x% N" Z- z, Z2 v/ j
式中,t表示计算系统响应的抽样点向量;f是系统输入信号,sys是LTI系统模型,用来表示微分方程。& a( m4 w; P. r, |8 V/ n+ B @. m3 [
% ?( h ?8 ^) p; T4 R; @系统模型sys要借助tf函数获得,其调用方式为
: g8 m6 Y' H9 k& `$ Q& M. m+ ?, ^
; M) @7 |7 E2 J9 {2 h$ Z% i3 i7 s- sys=tf(b, a)- f0 I+ ^- R* U/ q9 ~' [+ j
) `$ v+ D A5 y# p: M
5 }& l+ J; b' N$ u% ?* B( S5 |" J, }式中,b和a分别为微分方程的右端和左端各项的系数。
) |; i0 C! Z# `% U+ t比如:3 ~9 {! ~# E8 _3 w- F" ?* K. Q g
; C9 o% B7 d1 ]& F; z
8 S5 m: E `' a
7 N l6 v f! h8 n0 x' D3 l
. R8 B7 s3 j- k
# p9 m& v$ n0 ?3 s1 [) A% \
# y# v. t2 L5 a( J8 ]/ ?! W
- 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)')
8 `% {" M4 h" s, K( X! {
# M# ]4 C2 d" a, W. U
$ z: I' S3 i& `: U, j
/ F! h8 H. m6 `0 {- f& b: q7 S
$ N/ @; r6 [" }2 [0 r |
|