SciPy curve_fit 显示直线并且不适合数据

Question

我正在尝试将一组数据拟合为 CDF 指数函数。但是，我不确定我的代码或初始参数猜测出了什么问题，但它只会创建一条直线。数据是从 CSV 文件导入的。

#Plot Data
plt.figure(1,dpi=120)
plt.title("Cell A3")
plt.xlabel(rawdata[0][0])
plt.ylabel(rawdata[0][1])
plt.scatter(xdata,ydata,label="A3 Cell 1")



#Define Function
def func(t,lam):
    return 1 - (np.exp(-lam * t))

funcdata = func(xdata,1.17)
plt.plot(xdata,funcdata,label="Model")
plt.legend() 

#CurveFit data to model
popt, pcov = curve_fit(func,xdata,ydata,p0=(-0.64))
perr = np.sqrt(np.diag(pcov))

Image of the graph I get with the initial data and the straight line that the curve_fit gives

Answer 1

想想你的数据和你的功能。 ydata 是一个相当大的值。

的最大值是多少

def func(t,lam):
    return 1 - (np.exp(-lam * t))

我想你会发现函数的最大值出现在 lam 接近无穷大时，函数接近 1。最大值 == 1 的函数如何拟合 1000 秒内的数据？如果您希望能够扩展到 1 以上，则您的函数中需要更多参数。试试

def func(t,lam,scale):
    return ( 1 - (np.exp(-lam * t)) ) * scale

看看 scipy 是否能够更好地拟合数据。

编辑：我设法使它起作用，但是，您甚至没有绘制最佳参数。为此，请使用模拟的 xdata 和 ydata 查看我的代码：

#Plot Data
import numpy as np
from scipy.optimize import curve_fit
from matplotlib import pyplot as plt

def func(t,lam,scale):
    return ( 1 - (np.exp(-lam * t)) ) * scale

xdata = np.arange(25.)
ydata = func(xdata, 1.12, 2000.)

plt.figure(1,dpi=120)
plt.title("Cell A3")
plt.xlabel(rawdata[0][0])
plt.ylabel(rawdata[0][1])
plt.scatter(xdata,ydata,label="A3 Cell 1")

#CurveFit data to model
popt, pcov = curve_fit(func,xdata,ydata,p0=[0.5, 1000.1])
plt.plot(np.arange(25),func(np.arange(25), *popt),label="Model")
plt.legend()

输出：

Answer 2

您无法正确拟合这种简单的指数函数：

y=( 1 - (np.exp(-lam * t)) ) * 比例尺

数据，因为这个函数的形状在0

最好考虑逻辑类的函数，例如：

SciPy curve_fit 显示直线并且不适合数据

SciPy curve_fit displays straight line and is not fit to the data

python

curve-fitting

scipy