|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
Matlab求解系统的响应9 `# S, _& I( j$ O) Q# k
- `) [- O2 ~" M" T求LTI系统的零状态响应的函数lsim，其调用格式为. {, b1 ]7 M2 U$ ~. ?5 b
/ d( H! @' t( k5 Z
式中，t表示计算系统响应的抽样点向量；f是系统输入7 v d4 H1 Q' b, e5 r2 O$ `
信号，sys是LTI系统模型，用来表示微分方程。* ` d5 A- A0 O; z) R
0 {0 C- m) ^, y系统模型sys要借助tf函数获得，其调用方式为8 S$ k( h6 ]4 F* {& h
- ^9 M& R5 r* c1 @3 d9 F$ ?" W
sys=tf(b, a)9 Z8 U6 C/ E7 t, b' l! ?
% o. u' f6 V+ B
* {9 S9 M5 z9 h1 ?" \; \
式中，b和a分别为微分方程的右端和左端各项的系数。
9 \. [0 w$ d7 X- ^% |; f比如：
+ j; e" b# g' i+ o x" A( [; e4 _- u& k) V2 o, k
2 Y0 i) V2 p ~/ s% ?+ |
! q( L* l r, j( N$ h( E
: o6 s- \" v' P, v8 A
" h7 v2 N; m; p
1 P4 \, n6 @8 Q" m/ Msys=tf([1],[1,2,77]); % tf函数获得系统模型sys/ H5 H) o0 d8 Y* n7 R$ H9 J
t=0:0.01:5; % 对时间t进行离散取样
# a8 h: \ J; @7 T* F8 cf=10*sin(2*pi*t);
. w* h" j. e5 hy=lsim(sys,f,t); % 求系统的零状态响应
& @' W2 A! u z" C9 R8 }plot(t,y)
+ G! ^6 V8 K9 E Z8 O; U( kxlabel('Time(sec)')
7 ^; x: |( A W. S0 T7 I4 m+ `3 l' wylabel('y(t)')4 H4 A! F6 @2 e0 N% @" W3 D
0 w- a+ l; @) U' u
) }; }- W( L$ f
8 H, f2 S9 a3 @4 A y' \$ M; w. i3 _$ D# K
|
|
/1 
发表于 2022-4-1 16:33
收藏
淘帖
支持！
反对！