Попытка вывести параметр модели Лотки-Вольтерры ⇐ Python

[Anonymous]

def производная(X, t, A, B, C, D): х, у = Х dotx = x * (A - B * y) doty = y * (-D + C * x) вернуть np.array([dotx, doty]) интеграция по определению (t,A,B,C,D,X0): res = odeint(производная, X0, t, args = (A,B,C,D)) вернуть разрешение Х0 = [30, 4] X = массив([[30., 4.], [47.2, 6.1], [70.2, 9.8], [77.4, 35.2], [36.3, 59.4], [20.6, 41.7], [18.1, 19. ], [21.4, 13.], [22. , 8.3], [25.4, 9.1], [27.1, 7.4], [40.3, 8. ], [57. , 12.3], [76.6, 19.5], [52.3, 45.7], [19.5, 51.1], [11.2, 29.7], [7.6, 15.8], [14.6, 9.7], [16.2, 10.1], [24.7, 8.6]]) t = [0,0, 1,0, 2,0, 3,0, 4,0, 5,0, 6,0, 7,0, 8,0, 9,0, 10,0, 11,0, 12,0, 13,0, 14,0, 15,0, 16,0, 17,0, 18,0, 19,0, 20,0] XData = т YДанные = X Curve_fit (интеграция, XData, YData) Итак, X — это мои данные, первый столбец — это вид x, а второй столбец — вид y. Я попытался определить параметры этой модели Лотки-Вольтерры, используя оду и аппроксимацию кривой. Ошибка говорит о том, что для распаковки недостаточно значений (ожидалось 2, получено 1). На самом деле я даже не уверен, следует ли мне таким образом выводить параметр.

Может ли кто-нибудь помочь мне с этим, есть ли лучшие методы вывода параметров. Заранее спасибо!