假设我有以下方程组。我试图弄清楚将其作为 ODE 系统或 DAE 系统来解决对我来说是否更好。真正的代码可以有多达十几个方程，并且系统是僵硬的（我正在处理化学动力学，其中速率常数可以变化多个数量级）。

管理 ODE 是：

$$\frac{d\theta_{A}}{dt}=2R_{1}(\theta_{A},\theta_{*})-R_{2}(\theta_{A},\theta_{*})$$

$$\frac{d\theta{*}}{dt}=-2R_{1}(\theta_{A},\theta_{*})+R_{1}(\theta_{A},\theta_{*})$$

我可以解决$\theta_{A}$和$\theta_{*}$通过解决这个系统。然而，也有一个守恒定律使得

$$\theta_{A}+\theta_{*}=1$$

这意味着我可以求解方程组，例如，第一个 ODE 和代数守恒定律得到$\theta_{A}$和$\theta_{*}$也是。

在计算上，哪个更容易用数值求解？我更关心准确性而不是速度，但无论哪种方式都很好奇。解决两个 ODE 的系统或 ODE 和代数守恒定律之一的系统更容易吗？

如果它是相关的，我目前正在考虑ODE15s在 MATLAB 中使用，它可以解决 ODE 或 DAE 的僵硬系统。