如何将两段ode45求解结果合并绘制在同一张Matlab图中？

我来帮你搞定这个分段求解微分方程后合并绘图的需求，其实操作起来很顺手，咱们一步步来实现：

步骤1：确保你的微分方程函数定义正确 #

首先，你得有一个定义微分方程的函数rkfunc，这里我先给个简单示例（你可以替换成自己的实际方程）：

function dydt = rkfunc(t, y)
    % 这里可以根据需要写分段的微分方程，或者保持统一
    % 示例：简单的一阶线性微分方程 dy/dt = y
    dydt = y;
end

步骤2：分段求解ODE #

按照你已经写的代码执行两段求解，注意保存好每段的时间和结果：

% 第一段求解：t ∈ [0, 2]
ainitial = 0;
arange = [0 2];
[a, A] = ode45(@rkfunc, arange, ainitial);

% 第二段求解：t ∈ [2, 4]，初始值取第一段的最后一个结果
binitial = A(end);
brange = [2 4];
[b, B] = ode45(@rkfunc, brange, binitial);

步骤3：合并时间与结果向量 #

注意两段的时间在t=2处会有重复的点，咱们可以去掉第二段的第一个点，避免绘图时出现冗余：

% 合并时间向量：第一段全部 + 第二段从第二个点开始（跳过t=2的重复点）
combined_t = [a; b(2:end)];
% 合并结果向量：同理跳过第二段的第一个结果
combined_y = [A; B(2:end)];

步骤4：绘制合并后的结果 #

现在就可以把两段结果画在同一张图里了，还可以给两段设置不同的样式方便区分：

figure;
% 绘制第一段结果，用蓝色实线
plot(a, A, 'b-', 'LineWidth', 1.5);
hold on;
% 绘制第二段结果，用红色虚线
plot(b, B, 'r--', 'LineWidth', 1.5);
% 或者直接绘制合并后的完整曲线
% plot(combined_t, combined_y, 'k-', 'LineWidth', 1.5);

% 添加图例、标题和坐标轴标签
legend('第一段 (t∈[0,2])', '第二段 (t∈[2,4])');
title('分段微分方程求解结果');
xlabel('时间 t');
ylabel('解 y(t)');
grid on;
hold off;

额外小提示 #

如果你的微分方程本身就是分段定义的（比如t<2和t≥2时方程不同），其实可以不用分两次调用ode45，直接在rkfunc里根据t的范围切换方程即可，示例如下：

function dydt = rkfunc(t, y)
    if t < 2
        % t<2时的微分方程
        dydt = y;
    else
        % t≥2时的微分方程
        dydt = -y;
    end
end

这时只需要一次调用ode45就能得到完整区间[0,4]的解，再直接绘图就行，是不是更省事？

内容的提问来源于stack exchange，提问作者obizues