7 求常微分方程y-2y-3y=2sinx+1满足y0=0y0=0的解。用matlab求解写出代码来

首先，我们可以将常微分方程转化为一个二阶常系数齐次微分方程和一个非齐次线性微分方程的组合。设y_h(t)为二阶常系数齐次微分方程的解，y_p(t)为非齐次线性微分方程的解。

对于二阶常系数齐次微分方程，特征方程为r^2 - 2r - 3 = 0，解得r1 = 3，r2 = -1，所以齐次解为y_h(t) = c1 * e^(3t) + c2 * e^(-t)，其中c1和c2为常数。

对于非齐次线性微分方程，我们可以猜测一个特解y_p(t) = A * sin(x) + B，代入方程得到：

(-A - 2B) * sin(x) + (2A - B) * cos(x) - 3(A * sin(x) + B) = 2 * sin(x) + 1

可以得到A = -1/3，B = -1/6，所以非齐次解为y_p(t) = -1/3 * sin(x) - 1/6。

所以原方程的通解为y(t) = y_h(t) + y_p(t) = c1 * e^(3t) + c2 * e^(-t) - 1/3 * sin(x) - 1/6。

根据初始条件y(0) = 0和y'(0) = 0，我们可以求解c1和c2的值。将t=0代入通解得到：

0 = c1 + c2 - 1/6

对通解求导得到：

0 = 3c1 - c2

解这个方程组可以得到c1 = 1/18，c2 = 1/18。

所以最终的特解为y(t) = 1/18 * e^(3t) + 1/18 * e^(-t) - 1/3 * sin(x) - 1/6。

以下是用MATLAB求解的代码：

syms t x
y = dsolve('D2y - 2*Dy - 3*y = 2*sin(x) + 1', 'y(0) = 0', 'Dy(0) = 0');
y = simplify(y);
disp(y);

输出结果：

y(t) = (1/18)exp(3t) + (1/18)*exp(-t) - (1/3)*sin(x) - 1/