|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
, `, u+ V& `& [1 k. ~( ?2 ?! ?
Matlab求解系统的响应7 f: X( b8 B! G! J0 C
求LTI系统的零状态响应的函数lsim，其调用格式为1 ?. R- P8 U3 Y8 v8 A/ y6 p4 V
, F3 T$ y# M0 e3 C" E0 y
式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
8 Y4 G: F' e2 C W/ j) W4 m* t, Y4 c6 S
系统模型sys要借助tf函数获得，其调用方式为
/ T }, {7 P3 y m( z" l' F+ e8 {/ b
- sys=tf(b, a)' \8 o8 _! ^ A/ S! e1 a0 ^
+ f4 \4 c. ~, z4 U& |0 |4 s! Z" H
式中，b和a分别为微分方程的右端和左端各项的系数。
4 S4 e% @1 b. ]; F3 i& w比如：- K% K/ C/ L1 { }0 x# R+ {
( E" u5 F% j! N' C$ b, g( ?
. D5 F- D% P6 t) e" [+ n/ n/ c( A
' B. Z( T3 o7 i4 s6 Q I
+ l8 O. i1 Q9 y: j y$ M8 c2 ?" r3 w% z( T4 j! S- _
. D$ M: b3 C) }) j
- 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)')
A v# D" W5 q1 D- x
; @: I- O1 c, T/ U {" k$ K% b
9 u/ `; j6 }( a8 W b3 `9 L
; j4 f1 t. A: C) j
* A; S/ T. E. E5 F" f1 k |
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！