核心答案是 ode45 是 MATLAB 里面最常用的常微分方程求解器，采用自适应变步长的龙格库塔四五阶算法，也就是 Dormand-Prince 算法，适合求解非刚性常微分方程初值问题，绝大多数普通 ODE 问题优先选 ode45，它会自动调整步长，兼顾精度和计算速度，如果是刚性系统，ode45 计算会非常慢甚至失败，这个时候就要换 ode15s 这类刚性求解器，使用核心格式是 [T,Y] = ode45 (odefun,tspan,y0)，需要提前写好微分方程右端函数，设置时间区间和初始条件，这是 MATLAB 数值计算、自动控制、建模仿真课程必考内容。
基础原理
ode45 属于自适应步长 RK45，同时计算 4 阶和 5 阶解，利用两者差值估算局部截断误差，误差大就缩小步长，误差小就放大步长，默认精度适中，是 MATLAB 推荐首选求解器。
适用：非刚性 ODE，常微分方程初值问题，可以是单变量一阶 ODE，也可以是一阶微分方程组，高阶微分方程必须先降阶转换成一阶方程组，才能调用 ode45。
不适合：刚性系统，方程里面时间尺度差异非常大，变量变化快慢悬殊，ode45 迭代会非常低效。
基础调用语法对照表
| 语法形式 | 说明 |
|---|---|
| [T,Y] = ode45(odefun,tspan,y0) | 最常用，返回时间点 T 和解矩阵 Y |
| ode45(odefun,tspan,y0) | 不接收返回值，自动绘制解曲线 |
| [T,Y] = ode45(odefun,tspan,y0,options) | 带 options 选项，设置精度 事件等参数 |
参数解释
odefun：微分方程右端函数句柄，函数输入是 t 和 y，输出是 dydt 导数向量
tspan：时间区间，可以写 [t0 tf]，也可以写离散时间向量 [t0:dt:tf]
y0：初始状态向量
options：通过 odeset 创建，用来修改相对误差、绝对误差、事件函数等设置
第一步 高阶微分方程降阶（重点）
ode45 只能处理一阶微分方程组，比如二阶方程 y'' + 2y' + y = 0
设状态变量 y1 = y，y2 = y'
得到一阶方程组
dy1/dt = y2
dy2/dt = -2y2 - y1
完整示例代码
例子 求解二阶阻尼振荡方程 y'' + 0.5y' + y = 0，初始条件 y (0)=1，y'(0)=0
% 定义微分方程函数
odefun = @(t,y) [y(2); -0.5y(2)-y(1)];
tspan = [0 20];
y0 = [1;0];
[T,Y] = ode45(odefun,tspan,y0);
plot(T,Y(:,1),T,Y(:,2))
legend('y','dy/dt')
xlabel('t')
Y 矩阵每一行对应一个时间点，Y (:,1) 是 y1，Y (:,2) 是 y2
odeset 常用选项
RelTol 相对误差容限，默认 1e-3，想提高精度可以改成 1e-6
AbsTol 绝对误差容限
Events 事件函数，用来求方程到达某个阈值的时刻，比如过零点
MaxStep 最大允许步长，防止步长过大
事件函数示例作用，当状态变量达到指定值时，求解器自动停止计算。
常见误区
1 误区 ode45 可以直接求解高阶微分方程
纠正 必须手动降阶，写成一阶微分方程组，否则报错
2 误区 所有 ODE 都用 ode45
纠正 刚性系统优先 ode15s ode23s，ode45 跑刚性系统极慢
3 误区 tspan 写离散时间点就代表固定步长
纠正 即使 tspan 是离散序列，ode45 内部依旧是自适应步长，只是输出插值到指定时间点，不是固定步长 RK4
4 误区 函数输出维度和 y0 维度不匹配，会直接报错，导数向量维度必须和初始向量一致
和其他求解器对比
ode23 低阶自适应，精度更低，计算更快
ode15s 刚性系统求解器，变阶 BDF 算法
ode113 多步求解器，高精度长时仿真
总结
ode45 是 MATLAB 首选非刚性 ODE 求解器，核心调用格式 [T,Y]=ode45 (odefun,tspan,y0)，高阶方程要降阶成一阶方程组，精度可以用 odeset 调整，刚性问题不要使用 ode45，更换专门刚性求解器，课程仿真、控制系统、动力学模型大量使用这个函数。
【常见问题】
问题 1，ode45 的用法，ode45 是固定步长还是变步长
回答 1，ode45 的用法，ode45 是自适应变步长龙格库塔算法，根据误差自动调整步长，不是固定步长求解器。
问题 2，ode45 的用法，为什么求解刚性微分方程 ode45 很慢
回答 2，ode45 的用法，刚性系统变量时间尺度差距很大，ode45 为了保证稳定性被迫把步长压得极小，计算效率很低，推荐改用 ode15s。
问题 3，ode45 的用法，ode45 可以直接求解二阶微分方程吗
回答 3，ode45 的用法，ode45 不能直接求解二阶微分方程，必须做变量替换降阶，转换成一阶微分方程组之后再调用 ode45。
【免责申明】
本内容为 MATLAB 函数教学科普，不同 MATLAB 版本函数参数存在细微差异，详细说明请参考 MathWorks 官方帮助文档。
