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

首先，我们需要将微分方程转化为矩阵形式。令y1(t) = y(t)，y2(t) = y'(t)，x1(t) = x(t)，x2(t) = x'(t)，则有：

y1'(t) = y2(t) y2'(t) = -2y1(t) - 3y2(t) + x2(t) + 3x1(t) x1(t) = e^(-3t) * ξ(t) x2(t) = -3e^(-3t) * ξ(t) + δ(t)

其中，δ(t)为单位冲激函数。

将上述方程转化为矩阵形式，可以得到：

Y'(t) = AY(t) + BX(t)

其中，

Y(t) = [y1(t); y2(t)] X(t) = [x1(t); x2(t)] A = [0 1; -2 -3] B = [0 1; 3 1]

接下来，我们可以使用MATLAB命令进行求解。

首先，定义相关矩阵和向量：

A = [0 1; -2 -3]; B = [0 1; 3 1]; X0 = [1; 2]; % 起始条件

然后，定义输入信号函数：

t = 0:0.01:5; % 时间范围 x = exp(-3*t); % 输入信号

接下来，使用MATLAB的ode45函数进行求解：

[t, Y] = ode45(@(t, Y) AY + B[x(t); -3exp(-3t)], t, X0);

其中，@(t, Y) AY + B[x(t); -3exp(-3t)]是一个匿名函数，用于表示微分方程。

最后，绘制零输入响应、零状态响应和完全响应的图像：

figure; subplot(3, 1, 1); plot(t, Y(:, 1)); % 零输入响应 title('Zero Input Response'); xlabel('t'); ylabel('y(t)');

subplot(3, 1, 2); plot(t, Y(:, 2)); % 零状态响应 title('Zero State Response'); xlabel('t'); ylabel('y''(t)');

subplot(3, 1, 3); plot(t, Y(:, 1) + Y(:, 2)); % 完全响应 title('Total Response'); xlabel('t'); ylabel('y(t)');

运行以上代码即可得到零输入响应、零状态响应和完全响应的图像。