python curve_fit() 错误消息只有大小为 1 的数组可以转换为 python 标量

Question

我得到了以下等式，并要求我找出一些不同的东西，其中之一是最小二乘估计。

我有实验确定的数据并编写函数。

xData = np.array([0, 2, 5, 10, 15, 20, 25, 30]) # hours
yData = np.array([ 0.398,  0.543,  1.56,  4.34,  7.22, 9.86 , 10.5, 10.6]) # biomass, g/L

def func1(t,x_infin,x_0,k):
    'nonlinear function in a and b to fit to data'
    g = x_infin/ (1+(((x_infin-x_0)/x_0)*math.exp(-k*t)))
    
    return g

然后我设置curve_fit()

initial_guess = [10.5,0.0,16.5] # intial guesses for x_infin, x_0 and k
pars, pcov = curve_fit(func1, xData, yData, p0 = initial_guess)

得到

only size-1 arrays can be converted to Python scalars

我不确定我在哪里传递了一个不应该传递的数组。我对文档中的 curve_fit() 的理解是我传入了函数 xdata ydata 然后我对参数的猜测。

谁能指出我哪里出错了？

Answer 1

curve_fit 在做它的事情时使用 t 参数的底层数组。

math.exp 仅适用于标量输入；使用numpy.exp 处理数组