curve_fit()使用python

时间:2016-10-26 15:11:22

标签: python numpy

def model(A, x, mu, sigma):
return A*exp(-((x-mu)**2)/(2*sigma**2))
from scipy.optimize import curve_fit
mu=np.mean(d_spacing_2)
sigma=np.std(d_spacing_2)
f=intensity_2
x=d_spacing_2
popt, pcov = curve_fit(model, A, x, mu, sigma)

TypeError:model()缺少2个必需的位置参数:' mu'和西格玛'

1 个答案:

答案 0 :(得分:1)

您使用的是完全错误的curve_fit。这是一个有效的例子,来自curve_fit和一些额外的绘图:

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit
def func(x, a, b, c):
    return a * np.exp(-b * x) + c

xdata = np.linspace(0, 4, 50)
y = func(xdata, 2.5, 1.3, 0.5)
ydata = y + 0.2 * np.random.normal(size=len(xdata))

popt, pcov = curve_fit(func, xdata, ydata,p0=[2,1,1])

plt.ion()
plt.plot(xdata,ydata,'o')
xplot = np.linspace(0,4,100)
plt.plot(xplot,func(xplot,*popt))

curve_fit的第一个输入参数是函数,第二个是数据的x值,第三个是y值。您通常也应该使用可选的输入参数p0,这是解决方案的初始猜测。