浮点相对精度或公差应用到“选项”字段可能是这种差异的原因。
典型的工作流解决代数拓扑的限制是使用“ode15s”,“ode23t”,或“ode15i。“解决这些让你指定拓扑代数约束通过指定一个单一质量矩阵和可能会更适合这些方程。
以下链接解决解决拓扑使用这些详细信息
然而,即使DAE解决者,不能保证完全满足的约束。这可能是不可能的因为浮点数的限制。
如果有约束满意到底是要求您的应用程序,一个象征性的解决者可能是一个更好的选择。这里有一个例子如何计算一个象征性的解决方案的一个常微分方程(ODE)使用MATLAB的“dsolve”功能: