Matlab非负 ODE 解

Posted 2023-03-15 code_kd

将 ODE 的解约束为非负解。施加非负性并不总是微不足道的，但有时由于方程的物理解释或由于解决方案的性质，这是必要的。您应该只在必要时对解决方案施加此约束，例如在没有它的情况下集成失败的情况，或者解决方案不适用的情况。

如果解决方案的某些组件必须是非负的，则使用odeset设置NonNegative这些组件的索引选项。此选项不适用于ode23s、ode15i或应用于质量矩阵问题的隐式求解器 ( ode15s、ode23t、 )。ode23tb特别是，您不能对 DAE 问题施加非负性约束，该问题必然具有奇异质量矩阵。