python - 有没有结合 scipy curve_fit 和 scipy solve_ivp 的好方法？

Question

我一直想知道是否有一种方法可以使用 scipys curve_fit 函数来确定可以使用 scipy 的solve_ivp 的 ode 的初始值问题的参数。

def lotkavolterra(t, z, a, b, c, d):
    x, y = z
    return [(a*x - b*x*y), (-c*y + d*x*y)]

具体来说： 例如图像，您知道您的数据遵循 Lotka-Volterra 方程（s.above），但由于噪声，您无法计算出系数（a，b，c，d）。确切地知道时间点和初始值。如何有效地获得缺失的系数？

注意在这个线程中已经讨论了关于在integrate_quad 中使用curve_fit scipy.optimize.curve_fit与scipy.integrate.quad 的定积分函数

标签： pythonscipycurve-fittingdifferential-equations