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

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

[y1''(t); y2''(t)] + 3[y1'(t); y2'(t)] + 2[y1(t); y2(t)] = [x1'(t) + 3x1(t); 0]

则可以得到矩阵形式的微分方程：

[d^2/dt^2, 3; 3, d^2/dt^2][y1(t); y2(t)] + [2, 0; 0, 2][y1(t); y2(t)] = [1, 3; 0, 0][x1(t); 0]

现在，我们可以使用MATLAB命令来求解微分方程。

首先定义矩阵A、B和C：

A = [0, 1; -3, -2]; B = [1, 3; 0, 0]; C = [1, 0; 0, 1];

然后，定义输入信号函数x(t)和起始条件向量Y0：

t = 0:0.01:5; % 定义时间范围 x = exp(-3*t).*heaviside(t); % 定义输入信号函数 Y0 = [1; 2]; % 定义起始条件向量

接下来，使用MATLAB的ode45函数求解微分方程的零输入响应：

[t, Y] = ode45(@(t, Y) AY + B[x(t); 0], t, Y0);

然后，我们可以分别提取出零输入响应的y1(t)和y2(t)：

y1_zero_input = Y(:, 1); y2_zero_input = Y(:, 2);

接下来，我们使用MATLAB命令来求解微分方程的零状态响应：

[t, Y] = ode45(@(t, Y) A*Y, t, Y0);

然后，我们可以分别提取出零状态响应的y1(t)和y2(t)：

y1_zero_state = Y(:, 1); y2_zero_state = Y(:, 2);

最后，我们可以计算完全响应：

y1_total = y1_zero_input + y1_zero_state; y2_total = y2_zero_input + y2_zero_state;

完整的MATLAB代码如下所示：

A = [0, 1; -3, -2]; B = [1, 3; 0, 0]; C = [1, 0; 0, 1];

t = 0:0.01:5; x = exp(-3*t).*heaviside(t); Y0 = [1; 2];

[t, Y] = ode45(@(t, Y) AY + B[x(t); 0], t, Y0); y1_zero_input = Y(:, 1); y2_zero_input = Y(:, 2);

[t, Y] = ode45(@(t, Y) A*Y, t, Y0); y1_zero_state = Y(:, 1); y2_zero_state = Y(:, 2);

y1_total = y1_zero_input + y1_zero_state; y2_total = y2_zero_input + y2_zero_state;

最后，我们可以绘制y1(t)和y2(t)随时间变化的图形：

figure; subplot(2, 1, 1); plot(t, y1_total); xlabel('t'); ylabel('y1(t)'); subplot(2, 1, 2); plot(t, y2_total); xlabel('t'); ylabel('y2(t)');

通过以上步骤，我们成功使用MATLAB求解了微分方程，并得到了系统的零输入响应、零状态响应和完全响应。该方法可推广到其他微分方程的求解，为理解和分析系统的动态特性提供了重要工具。