最小二乘法笔记

本文介绍了如何将线性回归问题转化为矩阵形式,通过最小二乘法求解参数a和b,以适应多元线性模型。给出了一个二次函数拟合的实际例子,展示了如何计算系数矩阵。
摘要由CSDN通过智能技术生成

资料:https://zhuanlan.zhihu.com/p/653406161

预备知识

1、矩阵的求导

在这里插入图片描述

一、问题

  • 想对一组数据进行线性回归,也就是用一条直线拟合数据(只要可以写成矩阵形式,未必是直线)。

1、变量形式

待拟合直线方程: y = a x + b y = ax + b y=ax+b
代价函数(二乘项): J ( a , b ) = 1 2 m ∑ i = 0 m [ ( a x i + b ) − y i ] 2 J(a,b) = \frac{1}{2m}\sum_ {i=0}^m[(ax_i+b)-y_i]^2 J(a,b)=2m1i=0m[(axi+b)yi]2

2、矩阵形式

待拟合方程: Y = X A = [ x 1 x 1 0 x 2 x 2 0 ] [ a b ] = [ y 1 y 2 ] Y = XA = \left[\begin{array}{c} x_1 & {x_1}^0 \\ x_2 & {x_2}^0 \end{array}\right] \left[\begin{array}{c} a \\ b \end{array}\right] = \left[\begin{array}{c} y_1 \\ y_2 \end{array}\right] Y=XA=[x1x2x10x20][ab]=[y1y2]
代价函数: J = ( X A − Y ) T ∗ ( X A − Y ) J = (XA-Y)^T*(XA-Y) J=(XAY)T(XAY)

二、最小二乘法推导

1、拆开代价函数并求导

在这里插入图片描述

2、整理上式

在这里插入图片描述

三、最小二乘法结论

在这里插入图片描述

四、最小二乘法例子

1、问题

待拟合方程: y = a x 2 + b x + c y=ax^2+bx+c y=ax2+bx+c

数据:

xy
14
26.4
38
48.8

2、计算

[ a b c ] = ( [ 1 1 1 4 2 1 9 3 1 1 4 1 ] T [ 1 1 1 4 2 1 9 3 1 1 4 1 ] ) − 1 [ 1 1 1 4 2 1 9 3 1 1 4 1 ] T [ 4 6.4 8 8.8 ] = [ − 0.4 3.6 0.8 ] \left[\begin{array}{c} a \\ b \\ c \\ \end{array}\right] ={({ \left[\begin{array}{c} 1 & 1 & 1 \\ 4 & 2 & 1 \\ 9 & 3 & 1 \\ 1 & 4 & 1 \\ \end{array}\right]^T\left[\begin{array}{c} 1 & 1 & 1 \\ 4 & 2 & 1 \\ 9 & 3 & 1 \\ 1 & 4 & 1 \\ \end{array}\right]})}^{-1}\left[\begin{array}{c} 1 & 1 & 1 \\ 4 & 2 & 1 \\ 9 & 3 & 1 \\ 1 & 4 & 1 \\ \end{array}\right]^T\left[\begin{array}{c} 4 \\ 6.4 \\ 8 \\ 8.8 \\ \end{array}\right] = \left[\begin{array}{c} -0.4 \\ 3.6 \\ 0.8 \\ \end{array}\right] abc =( 149112341111 T 149112341111 )1 149112341111 T 46.488.8 = 0.43.60.8

>> X=[1 1 1; 4 2 1;9 3 1;16 4 1]

X =

     1     1     1
     4     2     1
     9     3     1
    16     4     1

>> Y = [4;6.4;8;8.8]

Y =

    4.0000
    6.4000
    8.0000
    8.8000

>> A = inv(X.'*X)*X.'*Y

A =

   -0.4000
    3.6000
    0.8000
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值