我正在尝试使用 odeint 模块解决这个问题。但是第一个方程是第二个方程的函数。如果我在第一个方程中忽略 dw/dz，而第二个方程是第一个方程的函数。我可以简单地使用 odeint 来解决它。我可以通过正向微分通过 fsolve 求解这些方程。但是我可以将 Odeint 与 fsolve 一起使用吗？并同时求解前两个方程以找到。
我如何解决这些耦合微分方程？
计算科学
颂
2021-12-19 12:41:16
1个回答
我会简单地替换在第一个等式中通过其在第二个等式中给出的表达式，然后在左侧重新组合所有涉及的项.
然后，知道你的整体状态向量的值，你有一个明确的公式.
如果手动无法进行这种操作，那么您可能需要使用可以处理这种情况的隐式求解器（例如，参见 Steven Roberts 的评论），但我对这些了解不多。