求微分方程中的参数

Question

我有一组数据点，我想用这组数据找出一个微分模型的3个参数。我试图从here改编代码，但有以下例外：

File "/home/djipey/Thesis/DGL/fit/25_09/etalon.py", line 13, in f
return Vmax * S / (Km + S (1 + S / Ki))
TypeError: 'numpy.ndarray' object is not callable

模型是：

-dS/dt=Vmax*S/（Km+S（1+S/Ki））

这是我使用的代码。它是完全开箱即用的。在

^{pr2}$

你知道我错了吗？在

编辑：当我说丑的时候，我的意思是：

但我不得不提最后一件事。我在科学上不确定我能用我的数据来拟合这个模型。我只是想试试这个选择。在

Answer 1

改变

def f(S, t, Vmax, Km, Ki):
    return Vmax * S / (Km + S (1 + S / Ki)) 

到

^{pr2}$

注意乘法，你所做的是用一些参数调用S，这会导致错误，而你想用这些参数来乘以S

编辑：回复评论中的问题：

http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.curve_fit.html

popt:数组

Optimal values for the parameters so that the sum of the squared error of f(xdata, *popt) - ydata is minimized