使用matlab7.0求解线性时不变系统的响应

1. 求解系统的零状态响应和零输入响应

2y’’(t)+ 2y’(t)+ 77y(t) = f(t) 在t>=0时接入激励f(t)=10*sin(2 pi t), y’(0) =2， y(0) =1，求零状态响应，求零输入响应并图形展示
参考例句：
求解LTI系统的零状态响应函数 y = lsim(sys,f,t)
求解LTI系统的零输入响应函数y = initial(sys1,R0)
t为计算的抽样点向量 ，f是系统输入信号，sys是LTI系统模型
sys 需要借助tf函数获得 ,sys = tf(b,a) ： b,a分别是微分方程右端和左端各项的系数
sys = tf(b,a); % 获取系统模型
t = 0:0.01:5; % 对时间t进行离散抽样

零输入响应：
sys=tf([1],[2,2,77]);
t=0:0.01:5;
f=10*sin(2*pi*t);
y1=lsim(sys,f,t);
plot(t,y1);
xlabel('Time(sec)');
ylabel('y1(t)');

零状态响应：
t=0:0.01:5;
r0=[2,1];
sys=tf([1],[2 2 77]);
sys1=ss(sys);%将微分方程转换为状态方程
x=0*t;
y=lsim(sys1,x,t,r0);
plot(t,y);
xlabel(&