求解系统响应
求解系统零状态响应函数lsim
y=lsim(sys,f,t)
t:表示坐标范围及单位长度。
f:表示输入信号的函数。
sys:表示系统的微分方程。
获取系统的微分方程tf
sys=tf(b,a)
b:微分方程等号右端的系数向量。
a:表示微分方程等号左边的系数向量。
例子:
sys=tf([1],[1,2,77]);%tf函数获得系统模型sys
t=0:0.01:5;%对t进行离散抽样(0~5单位长度为0.01)
f=10*sin(2*pi*t);%输入函数
y=lsim(sys,f,t);%求系统的零状态响应
plot(t,y);%画图
xlabel('Time(sec)');
ylabel('y(t)');
求解冲激响应和阶跃响应
冲激响应函数impulse
impulse(b,a)
阶跃响应函数step
step(b,a)
a,b表示方程左右两端的系数向量
a=[7 4 6];%构造系数向量
b=[1 1];
subplot(2,1,1)
impulse(b,a);%求系统的冲激响应并做图
subplot(2,1,2)
step(b,a);%求系统的阶跃响应并做图