Pierre-Simon Laplace · 1749–1827

拉普拉斯

他把牛顿的天空算完了,又给了概率论一副骨架

2V=0\nabla^2 V = 0
每一点都等于邻居的平均 · 位势论、调和函数与整个数学物理的底座

拉普拉斯 1749 年生于诺曼底的博蒙昂诺日,家里是殷实农户。十九岁带着一封推荐信到巴黎找达朗贝尔,信没打动人,他连夜写了一篇力学论文送去,达朗贝尔当天就给他安排了军事学校的教职。此后他历经王朝、革命、帝国与复辟,做过拿破仑六个星期的内政部长(拿破仑事后刻薄地说他「把无穷小的精神带进了行政」),最后被波旁王朝封为侯爵。

他做的是两件极大的事。一是《天体力学》五卷:把牛顿之后一百年里散落的引力问题收拢成一个体系,并证明行星轨道的长期摄动是周期的而不是累积的——太阳系不会自己散架。二是《概率的分析理论》:贝叶斯推断的一般形式、中心极限定理、最小二乘的概率解释,都在这里成形。今天以他命名的东西横跨两个世界:拉普拉斯算子、拉普拉斯方程、拉普拉斯变换,以及那个著名的、后来被混沌与量子力学各驳倒一半的拉普拉斯妖

拉普拉斯肖像
让–巴蒂斯特·波兰·盖兰绘 · 凡尔赛宫藏,法国国立博物馆联合会(RMN)· Wikimedia Commons,公有领域。画的是他任帝国参议院掌玺大臣时的装束。

生平

  1. 1749
    生于博蒙昂诺日

    诺曼底的殷实农户之家。

  2. 1768
    带一封信去巴黎

    信没打动达朗贝尔,他连夜写的力学论文打动了。

  3. 1774
    「逆概率」的一般写法

    他独立地给出了后来被称为贝叶斯推断的一般形式,并用它做天文数据的归算。

  4. 1796
    《宇宙体系论》与星云假说

    太阳系由一团旋转的气体云塌缩而成——这个想法沿用至今。

  5. 1799–1825
    《天体力学》五卷

    牛顿引力的集大成之作。行星的长期摄动是周期的,太阳系不会散架。

  6. 1812
    《概率的分析理论》

    中心极限定理的一般形式、生成函数、最小二乘的概率解释。

  7. 1814
    《概率哲学随笔》

    「拉普拉斯妖」出自这本书的开头几页——它其实是用来说明「概率是无知的度量」的。

  8. 1816
    改对了声速

    牛顿算的声速偏低约 16%,他指出压缩是绝热而非等温的,一改就对上了。

  9. 1827
    卒于巴黎

    同年牛顿式的天空已被他算完,而他留下的概率与位势论还要长两百年。

展品厅

绝大多数展品都可以亲手把玩——这是本馆的立馆之本;少数以叙述为主的,做成故事展签。

镇馆之宝 · 亲手玩

不管从哪里出发,加起来都是同一条钟形曲线

这就是正态分布在现实里无处不在的原因。

SnnμσnN(0,1)\frac{S_n - n\mu}{\sigma\sqrt{n}} \to \mathcal{N}(0,1)

棣莫弗 1733 年对二项分布做出了这件事——这座馆里高斯那块高尔顿板演的就是它。拉普拉斯把它推广到了一般情形:底下那个分布长什么样根本不重要,只要方差有限,独立地加足够多个,和的分布就趋于正态。这才是这条定理真正的分量,它解释了为什么正态分布无处不在:测量误差、身高、噪声,都是一大堆互不相干的小影响加在一起的结果。演示里四个源分布都特意挑得一点不像正态,其中「双峰」中间根本是空的。也留意一件事:偏得越厉害收敛得越慢,「极偏」那一栏推到 60 还差着一点——定理说的是极限,没说多快。

换四个一点不像正态的分布,看它们被碾成同一个形状

棣莫弗 1733 年对二项分布证明了这件事——这座馆里高斯那块高尔顿板演的就是它。 拉普拉斯把它推广到了一般情形: 底下那个分布长什么样根本不重要,只要方差有限, 独立地加足够多个,和的分布就趋于同一条钟形曲线。 这才是这条定理真正的分量:它解释了为什么正态分布在现实里无处不在—— 测量误差、身高、噪声,都是一大堆互不相干的小影响加在一起的结果。 四个源分布都特意挑得一点不像正态,其中「双峰」中间根本是空的。 把滑块往右推,看它们怎么被碾成同一个形状。 也留意一件事:偏得越厉害收敛得越慢, 「极偏」那一栏推到 60 还差着一点。定理说的是极限,没说多快。

把微分方程变成一道代数题

求导变成乘 s;然后看一眼极点,就知道系统长什么样。

L{y}=sY(s)y(0)\mathcal{L}\{y'\} = sY(s) - y(0)

拉普拉斯变换之后,求导就变成了乘 s,于是一整条微分方程塌成一个分式。解微分方程变成了解代数方程——这是第一层。更要紧的是第二层:分母的根(极点)直接读出响应长什么样。实部管衰减快慢,虚部管振荡频率;极点全在左半平面就稳定,越靠左衰得越快,虚部越大晃得越密。今天每一个做控制的人看系统,看的就是这张极点图;信号处理里的传递函数、电路分析里的阻抗,走的都是这条路。演示里把「按极点写出的解析解」和「直接数值积分那条微分方程」画在一起,两条必须重合——这才是「变换是对的」的意思。

拉普拉斯变换之后,求导就变成了乘 s。 于是一整条微分方程塌成一个分式:Y(s) = ω² / [s(s² + 2ζωs + ω²)]。 解微分方程变成了解代数方程——这就是全部的把戏。 更要紧的是第二件事:分母的根(极点)直接读出响应长什么样。 实部管衰减快慢,虚部管振荡频率;极点全在左半平面就稳定,越靠左衰得越快, 虚部越大晃得越密。今天每一个做控制的人看系统,看的就是这张极点图。 右边画了两条线:一条是按极点写出的解析解,一条是直接数值积分那条微分方程。两条必须重合——这才是「变换是对的」的意思,不是背下来的口诀。

每一点都等于邻居的平均

所以一块平衡的板子内部,不会凭空冒出一个热点。

u()=14u()u(\text{中}) = \tfrac{1}{4}\sum u(\text{邻})

ΔV = 0 听着抽象,可把拉普拉斯算子按五点差分写出来就很朴素:Δu = 上 + 下 + 左 + 右 − 4×中,于是 Δu = 0 就是「中间那点等于四个邻居的平均」。连续情形对应的是调和函数的均值性质:任一点的值等于它周围任一小圆上的平均。由此立刻得到最大值原理——既然每一点都不比邻居都大,内部就冒不出热点,最大值和最小值只能出现在边界上。而这个「反复用邻居平均替换」的动作,同时也是今天真正的数值解法。这个算子后来长到了哪里?热传导、静电、流体、图像处理里的边缘检测、图论里的图拉普拉斯矩阵、机器学习里的谱聚类,都是它。

拉普拉斯方程 ΔV = 0 听着抽象,可它的意思其实很朴素。 把拉普拉斯算子按五点差分写出来:Δu = 上 + 下 + 左 + 右 − 4×中。 于是 Δu = 0 就是「中间那点等于四个邻居的平均」。 连续情形下对应的是调和函数的均值性质:任一点的值等于它周围任一小圆上的平均。 由此立刻得到最大值原理:既然每一点都不可能比邻居都大, 内部就冒不出热点——最大值和最小值只能出现在边界上。 拖迭代滑块,看内部最大值怎么也越不过边界最大值。 而这个「反复用邻居平均替换」的动作,同时也是真正的数值解法: 今天解拉普拉斯方程,底子就是这一招。 换到「解析边界」那一栏,迭代会收敛到 e^(πx)sin(πy) 本身—— 因为它本来就是调和的。

拉普拉斯妖

一个被后世驳倒了两次的思想实验——但它原本不是用来吹牛的。

1814 年《概率哲学随笔》开头几页:假设有一个智者,某一瞬间知道宇宙中每一个粒子的位置与速度,又有足够的算力,那么对它来说未来与过去一样确定,没有任何事是不可知的。这段话后来被叫做「拉普拉斯妖」,成了机械决定论的招牌。可它在原书里的作用恰恰相反——他写这段是为了给出一个参照系,好说明接下来那句话:我们不是这样的智者,所以概率论是我们「无知的度量」,是弥补这个差距的工具。整本书讲的是概率,不是宿命。这只妖后来被驳倒了两次,而且是从两个方向:庞加莱发现的混沌说,初值只要差一丁点,长期就完全不可预测(这座馆里有他的三体演示);量子力学则说,那个「某一瞬间的精确位置与速度」本身就不存在。有意思的是,两次反驳并没有削弱拉普拉斯真正在做的事——反而正好证明了他为什么需要概率。

故事展签故事展签:他造这只妖,本来是为了说明我们不是它。

「我不需要那个假设」

太阳系不会自己散架——他把这件事算了出来。

牛顿留下过一个隐忧:行星互相牵扯,长期下去轨道会不会越偏越远,最后散掉?他倾向于认为需要上帝时不时出手校正。拉普拉斯在《天体力学》里把这件事算清楚了:那些看起来在单向漂移的长期摄动,其实是**周期极长的振荡**,会自己回来。木星和土星那个让哈雷困惑了半个世纪的「大不等式」,周期约九百年——在一代人的观测跨度里看着像单向漂移,实际上不是。太阳系不需要谁来扶。据阿拉果转述,拿破仑翻了他的书之后问:为什么书里一次也没提到上帝?他答:「陛下,我不需要那个假设。」这句话的转述链条不短,细节各版本有出入,但它精准地概括了这部书做成的事——把一个此前要靠信念填的缺口,用计算填上了。

故事展签故事展签:一个此前要靠信念填的缺口,被计算填上了。

传承

他的工作没有留在十八世纪——每条链的终点,都是你今天正在使用的东西。

中心极限定理统计推断与误差理论今天每一次抽样、每一个置信区间
拉普拉斯变换传递函数与极点分析自动控制、信号处理、电路设计
拉普拉斯算子 Δ位势论与偏微分方程图拉普拉斯、谱聚类、图像处理

语录

陛下,我不需要那个假设。

—— 据阿拉果转述,答拿破仑「为什么你的书里一次也没提到上帝」;转述链条不短,各版本细节有出入

概率论归根到底不过是把常识化成了算术。

—— 《概率哲学随笔》,1814(la théorie des probabilités n'est que le bon sens réduit au calcul)