Advertisement

最小二乘法(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)

还没有任何评论哟~