如何解决该系统的IVP微分方程

Question

我有一个像这样的微分方程组

dy/dt = f(t,y,y1,y2,....,yn,x), dy1/dt = f(t,y,y1,y2,..yn,x), . . . dyn/dt = f(t,y,y1,y2,..yn,x), x(t) = f(t,y1,y2,...yn,x)

我的值是y_i（0），x（0） 如果我有dx / dt，那么只需使用scipy.integrate IVP就可以解决此问题。我无法计算x（t）的导数，因为它太复杂了。我可以在python中模拟该系统而不找到dx / dt吗？

Answer 1

没问题，定义

def derivatives(t,y):
    xt = func_x(t,y);
    dy = func_y(t,y,xt)
    return dy

其中func_y接受标量（t，xt）和向量（y）参数，并返回与y相同大小的向量。然后，它可以与solve_ivp一起很好地工作。