全面加强政府网站建设/网络营销策划方案ppt模板

如何使用SciPy曲线拟合函数拟合高斯曲线来关闭此错误？换句话说，如果它不符合模型的峰值，那么它就不是峰值，所以我不想返回任何东西。还有，有没有更快的方法？曲线拟合可能太慢，我的应用程序看大数据。在RuntimeError: Optimal parameters not found: Number of calls to function has reached maxfev = 800.from scipy.optimize import curve_fit

from scipy import asarray as ar,exp

import matplotlib.pyplot as plt

from numpy import sqrt, pi, exp, loadtxt

data = loadtxt('data/model1d_gauss.dat')

x = data[:, 0]

y = data[:, 1]

n = len(x) #the number of data

mean = sum(x*y)/n #note this correction

sigma = sum(y*(x-mean)**2)/n #note this correction

def gaus(x,a,x0,sigma):

return a*exp(-(x-x0)**2/(2*sigma**2))

def gaussian(x, amp, cen, wid):

"1-d gaussian: gaussian(x, amp, cen, wid)"

return (amp/(sqrt(2*pi)*wid)) * exp(-(x-cen)**2 /(2*wid**2))

popt,pcov = curve_fit(gaus,x,y,p0=[1,mean,sigma])

#popt,pcov = curve_fit(gaussian,x,y,p0=[5,1,1])

plt.plot(x,y,'bo:',label='data')

plt.plot(x,gaus(x,*popt),'ro:',label='fit')

plt.legend()

plt.show()`enter code here`