scipy curve_fit 使用例子

一个使用curve_fit来模拟曲线的例子

# From https://docs.scipy.org/
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

# Define func of parameters a, b, c and variable x
def func(x, a, b, c):
    return a * np.exp(-b * x) + c

# generate some noisy data and plot it
xdata = np.linspace(0, 4, 50)
y = func(xdata, 2.5, 1.3, 0.5)
np.random.seed(1729)
y_noise = 0.2 * np.random.normal(size=xdata.size)
ydata = y + y_noise
plt.plot(xdata, ydata, 'b-', label='data')

#fit a, b, c and overplot
popt, pcov = curve_fit(func, xdata, ydata)
plt.plot(xdata, func(xdata, *popt), 'r-',
                    label='fit: a =%5.3f, b=%5.3f, c=%5.3f' % tuple(popt))

# Constrain the optimization to the region of 0 <= a <= 3.0 <= b <= 1 and 0 <= c <= 0.5
popt, pcov = curve_fit(func, xdata, ydata, bounds=(