2021-04-24 Python 最小二乘法求解线性回归模型

红太狼 2023-01-16 06:19 119阅读 0赞

## Python 最小二乘法求解线性回归模型 ##

机器学习线性回归模型

线性回归（linear regression）是一种线性模型，它假设输入变量 x 和单个输出变量 y 之间存在线性关系

具体来说，利用线性回归模型，可以从一组输入变量 x 的线性组合中，计算输出变量 y：y=ax+b

![图片][191fbafe07a3b0f508f4a6d78fbd7d09.png]

给定有d个属性（特征）描述的示例 x =（x1; x2; …; xd），其中xi是x在第i个属性（特征）上的取值，线性模型（linear model）试图学得一个通过属性（特征）的线性组合来进行预测的函数，即：

![20210424102245681.gif][]

一般用向量形式写成：

![图片][081ef61056f4f1dbde2c72240b042941.png]

其中 w = \{w1,w2,w3,,,,,,\}。

许多功能更为强大的非线性模型可在线性模型的基础上通过引入层级结构或高维映射而得。

## 最小二乘法求解 ##

基于均方误差最小化来进行模型求解的方法称为“最小二乘法”（least square method）

它的主要思想就是选择未知参数，使得理论值与观测值之差的平方和达到最小。

![图片][c8e426a5ae90af0a8abf7e6693fc75ec.png]

我们假设输入属性（特征）的数目只有一个：

![图片][642875868f4e2e0e601859c4a64d3884.png]

在线性回归中，最小二乘法就是试图找到一条直线，使所有样本到直线上的欧式距离之和最小。

![图片][dbf67bd024876d3e7137a91813d1b2bd.png]

求解线性回归

求解w和b，使得![图片][fa3f5c40d2be7e3e05f9b9940e4f8692.png]最小化的过程，称为线性回归模型的“最小二乘参数估计”。

将E(w,b)分别对w和b求导，可以得到

![图片][337c42c0f1d7df070bc0725de3a869e8.png]

![图片][5133155af301a12b13cdc47ca35520cb.png]

令偏导数都为0，可以得到

![图片][28e46fde8bc21ad209352091c13c6b8b.png]

![图片][ec9827842fcafa887bb6a62570125ec8.png]

—— 其中

![图片][b321edef82512176961c6f4ad132519b.png]

代码实现

### 0.引入依赖
    import numpy as np
    import matplotlib.pyplot as plt
    ### 1.导入数据
    points = np.genfromtxt('data.csv',delimiter=',')
    # 提取points中的两列数据，分别作为x,y
    x = points[:, 0]
    y = points[:, 1]
    # 用plt画出散点图
    plt.scatter(x, y)
    plt.show()
    
    
    ###  2.定义损失函数
    # 损失函数是系数的函数，另外还要传入数据的x，y
    def compute_cost(w, b, points):
        total_cost = 0
        M = len(points)
        
        #逐点计算平方损失，然后计算平均值
        for i in range(M):
            x = points[i, 0]
            y = points[i, 1]
            total_cost += (y - w*x -b) ** 2
        return total_cost/M
        
        
    ### 3.定义算法拟合函数
    # 先定义一个求均值的函数
    def average(data):
        sum = 0
        num = len(data)
        for i in range(num):
            sum += data[i]
            
        return sum/num
     
    # 定义核心拟合函数
    def fit(points):
        M = len(points)
        x_bar = average(points[:,0])
        
        sum_yx = 0
        sum_x2 = 0
        sum_delta = 0
        
        for i in range(M):
            x = points[i, 0]
            y = points[i, 1]
            sum_yx += y * (x - x_bar)
            sum_x2 += x ** 2
        # 根据公式计算w
        w = sum_yx / (sum_x2 - M * (x_bar ** 2))
        
        for i in range(M):
            x = points[i, 0]
            y = points[i, 1]
            sum_delta += (y - w * x)
        b = sum_delta / M
        
        return w, b
    
    
    ### 4. 测试
    w ,b = fit(points)
    print("w is: ", w)
    print("b is: ", b)
    cost = compute_cost(w, b, points)
    print("cost is : ",cost)
    ### 5.画出拟合曲线
    plt.scatter(x,y)
    #针对每一个x，计算得出预测的y值
    pred_y = w * x + b
    plt.plot(x, pred_y, c='r')
    plt.show()

Matlab版本感兴趣的读者可以自己试一下求解，![图片][84f04dd7663e83d99895f466fd7f30e8.png]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FpbmdmZW5neGQx_size_16_color_FFFFFF_t_70][]

[191fbafe07a3b0f508f4a6d78fbd7d09.png]: /images/20221022/317eff6eabba45538d287d92f9c86291.png
[20210424102245681.gif]: /images/20221022/20aa9b822463437ab5697b67236d7c94.png
[081ef61056f4f1dbde2c72240b042941.png]: /images/20221022/b3a557a428cd40fea3f0c843576a06e9.png
[c8e426a5ae90af0a8abf7e6693fc75ec.png]: /images/20221022/2297e2c93f6b4aca9017f493d5457b9e.png
[642875868f4e2e0e601859c4a64d3884.png]: /images/20221022/b355db74835d4fce9bd54c1f08146734.png
[dbf67bd024876d3e7137a91813d1b2bd.png]: /images/20221022/f740c719730b4e4ba3c5a40c8b06d412.png
[fa3f5c40d2be7e3e05f9b9940e4f8692.png]: /images/20221022/683237d7ac804c44bc2c43a0da088eb9.png
[337c42c0f1d7df070bc0725de3a869e8.png]: /images/20221022/a370f2c5b9854d05b8934d036e387a60.png
[5133155af301a12b13cdc47ca35520cb.png]: /images/20221022/50b10f511c1a4e259bffac2abf560c3f.png
[28e46fde8bc21ad209352091c13c6b8b.png]: /images/20221022/1bb50fbb44bc4230825a842c6ecc7838.png
[ec9827842fcafa887bb6a62570125ec8.png]: /images/20221022/8940cae676c9479e9e942b7103b666b0.png
[b321edef82512176961c6f4ad132519b.png]: /images/20221022/c5542f49cf3a4934a47447da15ec8fa3.png
[84f04dd7663e83d99895f466fd7f30e8.png]: /images/20221022/5e804888d24a49adae9abfde2f141183.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FpbmdmZW5neGQx_size_16_color_FFFFFF_t_70]: /images/20221022/a606b6d72fc443328cb7f7fd400ee9a7.png