如何理解最小二乘法?

如何理解最小二乘法?最小平方法是十九世纪统计学的主题曲 从许多方面来看 它之于统计学就相当于十八世纪的微积分之于数学 乔治 斯蒂格勒的 TheHistoryof 1 日用而不知来看一个生活中的例子 比如说 有五把尺子 用它们来分别测量一线段的长度 得到的数值分别为 颜色指不同的尺子 之所以出现不同的值可能因为 不同厂家的尺子的生产精度不同

最小平方法是十九世纪统计学的主题曲。 从许多方面来看, 它之于统计学就相当于十八世纪的微积分之于数学。

—-乔治·斯蒂格勒的《The History of Statistics》

1 日用而不知

来看一个生活中的例子。比如说,有五把尺子:

如何理解最小二乘法?

用它们来分别测量一线段的长度,得到的数值分别为(颜色指不同的尺子):

\begin{array}{c|c} \uad\uad&\uad长度\uad\\\hline \color{red}红& 10.2 \\\hline \color{blue}蓝& 10.3 \\\hline \color{orange}橙&9.8\\\hline \color{Goldenrod}黄&9.9\\\hline \color{green}绿&9.8\\\end{array}

之所以出现不同的值可能因为:

  • 不同厂家的尺子的生产精度不同
  • 尺子材质不同,热胀冷缩不一样
  • 测量的时候心情起伏不定
  • ……

总之就是有误差,这种情况下,一般取平均值来作为线段的长度:

\overline{x}=\frac{10.2+10.3+9.8+9.9+9.8}{5}=10

日常中就是这么使用的。可是作为很事’er的数学爱好者,自然要想下:

  • 这样做有道理吗?
  • 用调和平均数行不行?
  • 用中位数行不行?
  • 用几何平均数行不行?

2 最小二乘法

换一种思路来思考刚才的问题。

首先,把测试得到的值画在笛卡尔坐标系中,分别记作y_i :

如何理解最小二乘法?

其次,把要猜测的线段长度的真实值用平行于横轴的直线来表示(因为是猜测的,所以用虚线来画),记作y :

如何理解最小二乘法?

每个点都向y 做垂线,垂线的长度就是|y-y_i| ,也可以理解为测量值和真实值之间的误差:

如何理解最小二乘法?

因为误差是长度,还要取绝对值,计算起来麻烦,就干脆用平方来代表误差:

|y-y_i|\to (y-y_i)^2

总的误差的平方就是:

\epsilon=\sum (y-y_i)^2

因为y 是猜测的,所以可以不断变换:

如何理解最小二乘法?

自然,总的误差\epsilon 也是在不断变化的。

如何理解最小二乘法?

法国数学家,阿德里安-馬里·勒讓德(1752-1833,这个头像有点抽象)提出让总的误差的平方最小的y 就是真值,这是基于,如果误差是随机的,应该围绕真值上下波动(关于这点可以看下“如何理解无偏估计?”)。

这就是最小二乘法,即:

\epsilon=\sum (y-y_i)^2最小\implies 真值y

这个猜想也蛮符合直觉的,来算一下。

这是一个二次函数,对其求导,导数为0的时候取得最小值:

\begin{aligned}    \frac{d}{dy}\epsilon        &=\frac{d}{dy}\sum (y-y_i)^2=2\sum (y-y_i)\\        \quad\\        &=2((y-y_1)+(y-y_2)+(y-y_3)+(y-y_4)+(y-y_5))=0        \quad\\\end{aligned}

进而:

5y=y_1+y_2+y_3+y_4+y_5\implies y=\frac{y_1+y_2+y_3+y_4+y_5}{5}

正好是算术平均数。

原来算术平均数可以让误差最小啊,这下看来选用它显得讲道理了。

以下这种方法:

\epsilon=\sum (y-y_i)^2最小\implies 真值y

就是最小二乘法,所谓“二乘”就是平方的意思,台湾直接翻译为最小平方法。

3 推广

算术平均数只是最小二乘法的特例,适用范围比较狭窄。而最小二乘法用途就广泛。

比如温度与冰淇淋的销量:

\begin{array}{c|c} \uad\uad&\uad销量\uad\\\hline \color{red}{25^\circ}& 110 \\\hline \color{blue}{27^\circ}& 115 \\\hline \color{orange}{31^\circ}&155\\\hline \color{Goldenrod}{33^\circ}&160\\\hline \color{green}{35^\circ}&180\\\end{array}

看上去像是某种线性关系:

如何理解最小二乘法?

可以假设这种线性关系为:

f(x)=ax+b

通过最小二乘法的思想:

如何理解最小二乘法?

上图的i,x,y 分别为:

\begin{array}{c|c|c} \uad i\uad&\uad x\uad&\uad y\uad\\\hline 1&25& 110 \\\hline 2&27& 115 \\\hline 3&31&155\\\hline 4&33&160\\\hline 5&35&180\\\end{array}

总误差的平方为:

\epsilon=\sum (f(x_i)-y_i)^2=\sum (ax_i+b-y_i)^2

不同的a,b 会导致不同的\epsilon ,根据多元微积分的知识,当:

\begin{cases}    \frac{\partial}{\partial a}\epsilon=2\sum (ax_i+b-y_i)x_i=0\\    \quad\\    \frac{\partial}{\partial b}\epsilon=2\sum (ax_i+b-y_i)=0\end{cases}

这个时候\epsilon 取最小值。

对于a,b 而言,上述方程组为线性方程组,用之前的数据解出来:

\begin{cases}    a\approx 7.2\\    \quad\\    b\approx -73\end{cases}

也就是这根直线:

如何理解最小二乘法?

其实,还可以假设:

f(x)=ax^2+bx+c

在这个假设下,可以根据最小二乘法,算出a,b,c ,得到下面这根红色的二次曲线:

如何理解最小二乘法?

同一组数据,选择不同的f(x) ,通过最小二乘法可以得到不一样的拟合曲线(出处):

如何理解最小二乘法?

不同的数据,更可以选择不同的f(x) ,通过最小二乘法可以得到不一样的拟合曲线:

如何理解最小二乘法?

f(x) 也不能选择任意的函数,还是有一些讲究的,这里就不介绍了。

4 最小二乘法与正态分布

我们对勒让德的猜测,即最小二乘法,仍然抱有怀疑,万一这个猜测是错误的怎么办?

如何理解最小二乘法?

数学王子高斯(1777-1855)也像我们一样心存怀疑。

高斯换了一个思考框架,通过概率统计那一套来思考。

让我们回到最初测量线段长度的问题。高斯想,通过测量得到了这些值:

\begin{array}{c|c} \uad\uad&\uad长度\uad\\\hline \color{red}红& 10.2 \\\hline \color{blue}蓝& 10.3 \\\hline \color{orange}橙&9.8\\\hline \color{Goldenrod}黄&9.9\\\hline \color{green}绿&9.8\\\end{array}

每次的测量值x_i 都和线段长度的真值x 之间存在一个误差:

\epsilon_i=x-x_i

这些误差最终会形成一个概率分布,只是现在不知道误差的概率分布是什么。假设概率密度函数为:

p(\epsilon)

再假设一个联合概率密度函数,这样方便把所有的测量数据利用起来:

\begin{aligned}    L(x)        &=p(\epsilon_1)p(\epsilon_2)\cdots p(\epsilon_5)\\        \quad\\        &=p(x-x_i)p(x-x_2)\cdots p(x-x_5)\end{aligned}

讲到这里,有些同学可能已经看出来了上面似然函数了(关于似然函数以及马上要讲到的极大似然估计,可以参考“如何理解极大似然估计法?”)。

因为L(x) 是关于x 的函数,并且也是一个概率密度函数(下面分布图形是随便画的):

如何理解最小二乘法?

根据极大似然估计的思想,概率最大的最应该出现(既然都出现了,而我又不是“天选之才”,那么自然不会是发生了小概率事件),也就是应该取到下面这点:

如何理解最小二乘法?

当下面这个式子成立时,取得最大值:

\frac{d}{dx}L(x)=0

然后高斯想,最小二乘法给出的答案是:

x=\overline{x}=\frac{x_1+x_2+x_3+x_4+x_5}{5}

如果最小二乘法是对的,那么x=\overline{x} 时应该取得最大值,即:

\frac{d}{dx}L(x)|_{x=\overline{x}}=0

好,现在可以来解这个微分方程了。最终得到:

p(\epsilon)={1 \over \sigma\sqrt{2\pi} }\,e^{- {​{\epsilon^2 \over 2\sigma^2}}}

这是什么?这就是正态分布啊。

并且这还是一个充要条件:

x=\overline{x}\iff p(\epsilon)={1 \over \sigma\sqrt{2\pi} }\,e^{- {​{\epsilon^2 \over 2\sigma^2}}}

也就是说,如果误差的分布是正态分布,那么最小二乘法得到的就是最有可能的值。

那么误差的分布是正态分布吗?

我们相信,误差是由于随机的、无数的、独立的、多个因素造成的,比如之前提到的:

  • 不同厂家的尺子的生产精度不同
  • 尺子材质不同,热胀冷缩不一样
  • 测量的时候心情起伏不定
  • ……

那么根据中心极限定理(参考“为什么正态分布如此常见?”),误差的分布就应该是正态分布。

因为高斯的努力,才真正奠定了最小二乘法的重要地位。

文章最新版本在(有可能会有后续更新):如何理解最小二乘法?

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/214533.html原文链接:https://javaforall.net

(0)
上一篇 2026年3月18日 下午4:02
下一篇 2026年3月18日 下午4:02


相关推荐

  • 大数据面试题以及答案整理(一)

    大数据面试题以及答案整理(一)在进行大数据工作的面试中 将面试问题进行了汇总 然后对答案进行了整理

    2026年3月17日
    2
  • jQuery的下载以及使用[通俗易懂]

    jQuery的下载以及使用[通俗易懂]一、版本选择jQuery的版本有很多,大家在选择版本的时候,一般原则是越新越好,但其实不然,jQuery版本是在不断进步和发展的,最新版是当时最高技术水平,也是最先进的技术理念。如何选择jQuery版本是个值得思考的问题目前jQuery有三个大版本:1.x:兼容ie678,使用最为广泛的,官方只做BUG维护,功能不再新增。因此一般项目来说,使用1.x版本就可以了,最终版本:1.12.4(2…

    2022年5月22日
    38
  • 即梦怎么导入自己的音频

    即梦怎么导入自己的音频

    2026年3月12日
    2
  • Java学习之Spring框架入门篇

    Java学习之Spring框架入门篇0x00前言最近在学SSM的三大框架,上篇文章,已经更新了Mybatis的相关内容,那么这篇文章就来写写Spring的入门。0x01Spring概述S

    2021年12月12日
    44
  • mybatis拦截器执行顺序配置_springmvc拦截器执行顺序

    mybatis拦截器执行顺序配置_springmvc拦截器执行顺序1.原始jdbc工作流程原始jdbc工作流程以查询为例1加载驱动Class.forName(Driver.class.getName())2建立数据库连接Connectionroot=DriverManager.getConnection(“xx”,“xx”,“xx”)3预编译sql语句PreparedStatementpreparedStatement=root.prepareStatement(sql)4占位符参数赋值preparedSt

    2025年9月5日
    5
  • Java 版学生成绩管理系统,附源码[通俗易懂]

    Java 版学生成绩管理系统,附源码[通俗易懂]前言对于计算机专业的学生来讲,一定对于学生成绩管理系统课程设计十分印象深刻。设计、上机、编码、写报告、答辩等一系列流程,虽然很简单,但也可以综合运用我们所学的一些知识了。今天就来复习以下这个课题,用JavaSE来实现该课题,也算是补上当初上课的时候没有好好写的遗憾吧。虽然很简单,但是对于刚涉足编程学习的人来讲,还是有一定难度的,既要考虑界面打印,也要考虑条件判断、循环语句、输入输出控制等等技巧,所以在这里简单地实现一下,方便给初学的小伙伴们一个参考(对于我的界面比较丑的问题,就不要过于纠结了,下

    2022年7月13日
    25

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号