scipy读书报告

import numpy as np

from scipy.optimize import curve_fit

#创建函数f(x) = ax + b

def func(x,a,b):

    return a*x+b

 

#创建干净数据

x = np.linspace(0,10,100)

y = func(x,1,2)

 

#添加噪声

yn = y + 0.9* np.random.normal(size=len(x))

 

#拟合噪声数据

popt,pcov = curve_fit(func,x,yn)

 

#输出最优参数

print(popt)