最小二乘法(Python)
发布时间
阅读量:
阅读量
直接调用scipy.optimize.leastsq的实践与反思
下面的图片中左边是库拟合的,右边是我自己拟合的

自行体会上图
# 这里用的方程是y=k*x+b
# 输入的是列表list,分别是x和y的值
# 输出的是k和b的值
def least_square(x, y):
if len(x) != len(y):
return False
aver_x = sum(x) / len(x)
aver_y = sum(y) / len(y)
# 分母和分子初始化
denominator = 0
molecule = 0
for i in range(0, len(x)):
denominator += (x[i] - aver_x) * (y[i] - aver_y)
molecule += pow((x[i] - aver_x), 2)
k = denominator / mole
全部评论 (0)
还没有任何评论哟~
