如何在Python中用2个变量（+时间）求解2个微分方程

2024-05-20 • 问答

我在python中非常新，我一直在尝试解决2个具有2个未知数的联立微分方程组。两个变量T,X都是时间的函数，T的导数：dT/dt取决于T和X，而dX_dt仅取决于X。我的第一个想法是创建函数{{ 1}}返回fun(t,T,X)和dT/dt的值并使用dX/dt，但是效果不佳。有人可以告诉我解决这类ODE系统的正确方法是什么？提前致谢。这是脚本：

solve_ivp

ysqysq1 回答：如何在Python中用2个变量（+时间）求解2个微分方程

返回值（状态的导数）是元组（或概念上的向量）的方式相同，ODE函数还必须将状态接受为相同维的元组/向量。求解器基于通过向量值向量函数实现ODE系统。

def fun(t,u):
    T,X = u;
    # etc.

进行此更改后，您的程序应运行或至少给出更多有趣的错误。

您可能想尝试一下method="Radau"是否能以更少的内部步骤实现更快的集成，通常是这样。

ode solver

本文链接：https://www.f2er.com/2972449.html

如何在Python中用2个变量（+时间）求解2个微分方程

ysqysq1 回答：如何在Python中用2个变量（+时间）求解2个微分方程

大家都在问