|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
; h: T( }% P3 w2 C1 f9 B3 R# ]
Matlab求解系统的响应
9 }6 ~: \! r6 ]求LTI系统的零状态响应的函数lsim，其调用格式为8 b8 c& x* e4 ], N
' U/ U) U) T, b% ?, L2 \! J式中，t表示计算系统响应的抽样点向量；f是系统输入信号，sys是LTI系统模型，用来表示微分方程。
+ s W6 ^/ J' D
/ ]3 w# O$ A; m系统模型sys要借助tf函数获得，其调用方式为 g& N; \, n( e- @" ~
% N- [- n7 Z" t G
- sys=tf(b, a)4 |1 Z# K* Q1 n; W4 P; d
* L+ Z3 S c$ s8 i0 i7 e. p$ A& m; t
$ t" l1 w8 {$ b+ ?9 |. R; d1 B式中，b和a分别为微分方程的右端和左端各项的系数。" [' t, Q: s& m: W/ X" d
比如：6 W ]8 g: [ Y4 Y2 e) b) |
' c# @! E+ z2 v' ^0 z
0 s& P1 ^) j5 P8 N
% a9 w$ C( k2 a9 L; A) c
5 M' k$ s' n/ n1 _. ?# @7 g! |1 u+ A& n2 M8 ?, {1 o6 p
6 p# ~* H- V. o, B
- 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)')1 f8 h3 M+ s" I4 n5 j. U
# n4 f2 ~, k2 Y' ?( t3 ?8 E; Q
, u4 z) y' n4 v: g, Q2 C
; K! ~/ E2 Y- _2 y9 N
9 h9 d$ s0 F) K |
|
/1 
发表于 2021-7-27 10:09
收藏
淘帖
支持！
反对！