MATLAB 微分方程求解：零输入响应、零状态响应和完全响应

首先，我们需要求解微分方程的齐次解和特解。

齐次方程为：y''(t) + 3y'(t) + 2y(t) = 0

特征方程为：r^2 + 3r + 2 = 0

解特征方程得到特征根：r1 = -1, r2 = -2

齐次解为：y_h(t) = c1e^(-t) + c2e^(-2t)

其中，c1和c2为待定系数。

接下来，我们求解非齐次方程。

对于特解，我们可以猜测形式为：x_p(t) = A*e^(-3t)*ξ(t)

其中，A为待定系数。

将特解代入原方程得到：-9A*e^(-3t)ξ(t) + 3Ae^(-3t)δ(t) + 2Ae^(-3t)*ξ(t) = e^(-3t)*δ(t) + 3e^(-3t)*ξ(t)

整理得到：-9A*e^(-3t)ξ(t) + 2Ae^(-3t)*ξ(t) = e^(-3t)*δ(t) + 3e^(-3t)*ξ(t)

化简得到：-7A*e^(-3t)*ξ(t) = e^(-3t)*δ(t)

两边同时除以e^(-3t)得到：-7A*ξ(t) = δ(t)

由于ξ(t)在t=0时为1，其他时间为0，所以上式可以写为：-7A = δ(t)

因此，A = -1/7

特解为：x_p(t) = (-1/7)*e^(-3t)*ξ(t)

接下来，我们求解零输入响应、零状态响应和完全响应。

零输入响应为齐次方程的解：y_zi(t) = c1e^(-t) + c2e^(-2t)

由起始条件y(0) = 1和y'(0) = 2，可以得到：

c1 + c2 = 1 -c1 - 2c2 = 2

解这个方程组得到：c1 = 4/3，c2 = -1/3

零输入响应为：y_zi(t) = (4/3)*e^(-t) - (1/3)*e^(-2t)

零状态响应为特解：y_zs(t) = x_p(t) = (-1/7)*e^(-3t)*ξ(t)

完全响应为零输入响应和零状态响应之和：y(t) = y_zi(t) + y_zs(t) = (4/3)*e^(-t) - (1/3)*e^(-2t) - (1/7)*e^(-3t)*ξ(t)