刚性系统应选用ode15s而非ode45：依据特征值刚性比>10³、ode45失败步数占比超15%或ode15s耗时降50%以上判别；线性负定系统、含开关逻辑非线性系统及DAE必须用ode15s；可构建try-catch自动切换流程。
如果您在MATLAB中求解常微分方程时发现ode45计算异常缓慢、步长急剧收缩、反复报出“Failure at t=... . Unable to meet integration tolerances”警告，或数值解出现非物理振荡，则极可能遭遇了刚性系统。以下是针对刚性方程选择ode45与ode15s求解器的具体操作路径：
一、依据系统刚性特征判别求解器类型
刚性系统的本质是方程中存在显著分离的时间尺度，导致显式方法（如ode45）为保持稳定性被迫采用极小步长。判别需结合数学特性与数值行为，而非仅凭经验猜测。
1、计算系统雅可比矩阵的特征值，若最大与最小特征值实部之比（刚性比）大于10³，则判定为刚性系统；
2、在ode45调用中添加选项options = odeset('Stats','on')并运行，观察输出中“Number of successful steps”与“Number of failed attempts”的比值；若失败次数占比超过15%，即提示刚性倾向；
3、将同一方程改用ode15s运行，若其成功步数远超ode45且总耗时下降50%以上，则反向验证刚性存在。
二、基于方程结构匹配求解器算法机制
ode45采用显式Runge-Kutta (4,5)法，依赖函数光滑性与缓变动力学；ode15s基于隐式NDF（后向差分公式），天然具备A-稳定性和对刚性源项的强抑制能力。二者不可互换的核心在于数值稳定性域覆盖范围不同。
1、对含快速衰减模态的线性系统（如y′ = Ay，A为负定矩阵且谱半径极大），必须选用ode15s，因ode45的稳定性域无法覆盖左半平面深远处；
2、对具有分段光滑右端函数的非线性刚性问题（如带开关逻辑的电路模型），优先启用ode15s并设置'Jacobian'选项以提供解析雅可比，避免数值微分引入误差；
3、对弱刚性问题（刚性比介于10²–10³之间）且容差宽松（RelTol ≥ 1e−3），可尝试ode23s作为过渡方案，其单步隐式特性兼顾效率与稳定性。
三、通过容错响应动态切换求解器
在未知刚性程度的批量仿真任务中，可构建自动判别流程：先以ode45启动，捕获警告信息并实时评估步长退化率，触发条件满足即切换至刚性求解器，避免人工预判失误。
1、定义匿名函数时包裹try-catch结构，捕获ode45运行时的警告ID 'MATLAB:ode45:StepSizeTooSmall'；
2、在catch块内提取当前积分点t和解y，以此为初值调用ode15s继续积分；
3、设置全局标志位isStiffDetected = true，后续相同系统直接跳过ode45试探阶段。
四、利用质量矩阵与投影约束限定求解器适用性
当ODE系统耦合代数约束（如DAE指标1系统）或存在奇异质量矩阵M(y)y′ = f(t,y)，求解器选择不再仅由刚性决定，而受制于其对M矩阵的处理能力。ode15s原生支持质量矩阵输入，而ode45完全不兼容。
1、若方程形式为M(t,y)*dy/dt = f(t,y)且M为时变或非恒等矩阵，禁止使用ode45，必须传入odeset('Mass',@massfun)配合ode15s；
2、对含守恒律的刚性系统（如能量/动量约束），启用ode15s的'MassSingular','yes'选项，并配合'MStateDependence','none'提升隐式求解效率；
3、在调用前校验M矩阵秩：若rank(full(M)) ，则表明存在冗余自由度，<strong><font color="green">必须启用ode15s的质量矩阵模式</font></strong>而非尝试降阶处理。
