如果用另一种方法如何解？

一毛

已知描述连续时间系统的微分方程和激励信号f(t)分别如下;
y''(t)+4*y'(t)+2*y(t)=f''(t)+3*f(t)
f(t)=exp(-t)*u(t)
要求出系统在0~10秒时间范围内零状态响应y(t)的样值。
我用lsim函数很容易求出来，命令如下：
a=[1 4 2];
b=[1 0 3];
sys=tf(b,a);
t=0:0.1:10;
f=exp(-t);
y=lsim(sys,f,t)
这样回车就可以显示各个零状态响应样值
但是现在要求用求解微分方程的方法来实现，即求解系统的解析解，又是如何算的呢？
请各位论坛大虾指点一下迷津，谢谢啦！

一毛

没人知道方法吗？我还是不知该如何解啊！请各位帮帮忙啦。谢谢！