y”(t)+4y’(t)+4y(t)=f’(t)+3f(t)，f(t)=e -tu(t)用matlab求波形

要求解该微分方程及其初始条件并绘制其波形，可以使用MATLAB的符号工具箱和常微分方程求解器。以下是可能的示例代码：

syms t y(t) f(t)
eqn = diff(y,t,2) + 4*diff(y,t) + 4*y == diff(f,t) + 3*f;
f = exp(-t)*heaviside(t);
y0 = 2;
y1 = 0;
ySol(t) = dsolve(eqn, y(0)==y0, diff(y)(0)==y1);
tVals = linspace(0, 5, 1000);
yVals = double(ySol(tVals));
plot(tVals, yVals);
xlabel('t');
ylabel('y(t)');
title('Solution of y''''(t) + 4y''(t) + 4y(t) = f''(t) + 3f(t)');

其中，`syms`关键字声明`t`、`y(t)`和`f(t)`为符号变量。`eqn`表示方程式，使用`diff`函数来表示微分。`f`表示给定的函数，使用`exp`和`heaviside`函数来表示指数和阶跃函数。

然后，我们使用`dsolve`函数来求解微分方程和初始条件，得到解析解`ySol(t)`。最后，我们使用`linspace`函数生成1000个等间隔的时间值，然后使用`double`函数将解析解转换为数值解。最后，我们绘制波形图，其中X轴是时间，Y轴是y(t)的值。

该代码的输出应该是一张波形图，显示y(t)随时间变化的曲线。注意，由于使用了符号工具箱，该代码的执行可能需要一些时间。