matlab求解系统响应

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