【信号与系统】【MATLAB】二、求解系统的响应

零状态响应

% 求解系统的响应
% 求解LTI系统的零状态响应函数 y = lsim(sys,f,t)
% t为计算的抽样点向量 ，f是系统输入信号，sys是LTI系统模型
% sys 需要借助tf函数获得 ,sys = tf(b,a) ： b,a分别是微分方程右端和左端各项的系数
% 例：y''(t)+ 2y'(t)+ 3y(t) = f'(t) + 2f(t) 在t>=0时接入激励f(t)=5*sin(2 *pi *t),求零状态响应
sys=tf([1 2],[1 2 3]);
t=0:0.01:5;
f=10*sin(2*pi*t);
y=lsim(sys,f,t);
plot(t,y);
xlabel('t')
ylabel('y')

冲激响应和阶跃响应

% 冲激响应和阶跃响应 
% 以下的a、b就是各阶导数的系数，可参考上面例程的注释代码
a=[1 2 3];
b=[1 2];
subplot(2,1,1)
impulse(b,a); % 冲激响应
subplot(2,1,2)
step(b,a);    % 阶跃响应