在Matlab中, 求解连续时间LTI系统冲激响应可应用控制系统工具箱提供的函数impulse, 求解阶跃响应可利用函数step. 其调用方式为
y = impulse(sys, t)
y = step(sys. t)
式中 t 表示计算系统响应的抽样点向量, sys是连续时间LTI系统模型. 下面举例说明其应用.
% 连续时间LTI系统冲激响应的求解
% 在区间[0, 10]上选取5000个点
t = linspace(0,10, 5000);
% 连续时间LTI系统模型sys要借助tf函数获得
% 其调用方式为 sys(b, a)
% 其中b和a分别为微分方程右端和左端各项的系数向量.
% yzs''(t) + 3yzs'(t) + 2yzs(t) = 2x(t)
sys = tf([2], [1 3 2]);
%
y = impulse(sys, t);
% 绘制系统冲激响应波形
plot(t, y, 'r-');
%
xlabel('time(s)');
ylabel('h(t)');
%
legend('h(t)');
% 连续时间LTI系统阶跃响应的求解
% 在区间[0, 10]上选取5000个点
t = linspace(0,10, 5000);
%