我try 使用Scipy的curve_fit()函数将一个高斯函数与我的数据集进行拟合,但未能使该函数与之匹配.我使用其他一些工具,如MatLab,try 了同样的方法,该函数非常适合.有没有人能帮帮我?我不确定我做错了什么.非常感谢您的帮助:)
import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
x_data = [12, 34, 56]
y_data = [1e-10, 1e-3, 1e-10]
def func(xdata, a, b, c):
return a*np.exp(-(xdata - b)**2/(2*c**2))
popt,_ = curve_fit(func, x_data, y_data)
x_fit = np.linspace(0,100, 100)
y_fit = func(x_fit, *popt)
plt.scatter(x_data, y_data)
plt.plot(x_fit,y_fit)
plt.show()
上面是我试过的代码,我得到了一条钟形曲线,它拒绝从平均值0(钟形部分超过x=0)移动.