为何diffeqpy的隐式DAE求解器总是跳回t=0？

1. 跳回t=0现象的根源 #

这种往复跳回初始点的现象并非隐式求解方式本身导致，更可能由以下具体问题引发：

• DAE初值不满足约束：DAE系统要求初始时刻y0和dy_dt_0必须满足代数约束daefun(0, y0, dy_dt_0) = 0（允许微小数值误差）。若初值不兼容，DImplicitEuler会触发回溯修正逻辑，表现为跳回t=0重新调整初值。
• 迭代收敛失败：隐式求解器每一步都需要解非线性方程组，若迭代过程中残差无法收敛到阈值内，求解器会重置到之前的可靠计算点（如t=0），尝试调整步长或迭代策略后重新计算。
• DAE函数实现错误：如果daefun中存在逻辑漏洞（比如条件分支导致t的取值异常、代数约束表述错误），会直接触发求解器的异常回溯行为。

2. 依赖t的参数/变量对求解的影响 #

若DAE系统中存在显式依赖t的参数或变量，会从以下维度影响求解过程和结果：

• 约束的时变特性提升计算复杂度：每一步迭代都需要重新计算时变的代数约束，会增加计算量；若t的变化导致约束非线性程度突变，会降低迭代收敛性，甚至加剧步长剧烈调整（比如你遇到的跳点现象可能更频繁）。
• 初值兼容性要求更高：时变DAE不仅要求初始时刻满足daefun(0, y0, dy_dt_0)=0，还需确保初始约束的导数也符合系统要求，否则求解器的初值修正过程会反复触发。
• 精度与稳定性受影响：时变项可能改变系统的稳定性区域，若DImplicitEuler的步长调整未适配时变项的变化速率，会导致求解精度下降，甚至出现数值振荡、状态跳转错误等不稳定现象。

3. 排查与优化建议 #

• 验证初始条件：检查y0和dy_dt_0是否严格满足daefun(0, y0, dy_dt_0) ≈ 0，若不满足，可使用DifferentialEquations.jl的初值修正功能（比如在初始化问题时添加initializealg=BrownFullBasic()参数）。
• 调整求解器参数：尝试开启自动微分（autodiff=True），或手动设置初始步长dt、最大迭代次数maxiters，例如：

  sol = de.solve(dae_prob, de.DImplicitEuler(autodiff=True), dt=1e-3, maxiters=1000)
  

• 检查DAE索引：高索引DAE（索引2及以上）对隐式求解器兼容性较差，可先将问题降阶为索引1的DAE后再求解。

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