使用ModelingToolkit的structural_simplify将DAEs转ODE时方程丢失求助

structural_simplify()自动剔除方程通常是因为原DAE系统存在冗余（线性相关、可推导的代数约束）或非独立变量，工具会生成最小自由度的ODE系统。若你怀疑是异常丢失，可按以下步骤排查：

• 确认系统自由度与冗余性
  先计算系统理论自由度：自由度 = 微分变量数 - 独立代数约束数。用ModelingToolkit.dof(sys)查看原系统与简化后系统的自由度，若两者一致，说明丢弃的是冗余方程，属于正常行为。

• 验证简化系统的一致性
  代入初始条件求解简化后的ODE，再验证是否满足所有原方程：

  using ModelingToolkit, OrdinaryDiffEq
  # 假设sys为原DAE系统，simplified_sys为简化后系统，u0是初始条件，tspan是求解区间
  prob = ODEProblem(simplified_sys, u0, tspan)
  sol = solve(prob, Tsit5())
  # 代入原方程计算残差
  res = [substitute(eq, Dict(zip(states(sys), sol.u[1])), t => sol.t[1]) for eq in equations(sys)]
  println("原方程最大残差：", maximum(abs.(res)))
  

  若残差接近数值精度（如1e-10级别），说明简化是合理的。

• 强制保留关键方程
  若某些约束是建模必需的，可在构建系统时标记preserve=true避免被剔除：

  @variables x y z t
  D = Differential(t)
  # 标记需要保留的方程
  eq1 = D(x) ~ y + z
  eq2 = y^2 + z^2 ~ 1
  sys = ODESystem([eq1, eq2], t, [x,y,z], []; preserve=[eq2])
  simplified_sys = structural_simplify(sys)
  

• 排查建模错误
  检查原系统是否存在变量重复定义、方程线性相关（如同一约束的不同表述）、或微分变量误设为代数变量的情况，这些都可能导致过度简化。

• 查看简化细节
  启用verbose=true参数打印简化过程，明确被剔除的方程及原因：

  simplified_sys = structural_simplify(sys; verbose=true)
  

内容的提问来源于stack exchange，提问作者Thành Trần