如何用Matlab求连续时间系统的全响应

【例】求系统y”(t)+y(t)=cos(t)，y(0+)=y(0-)=0的全响应。
b=[1];
a=[1,0,1];
[A B C D]=tf2ss(b,a);
sys=ss(A,B,C,D);
t=0:0.1:10;
x=cos(t);
zi=[0 0];
y=lsim(sys,x,t,zi);
plot(t,y);
xlabel('时间');ylabel ('y(t)') ; title ('零状态响应');

参考里面的帮助就可以了。