all_lessons/数学的逻辑/29第 30 课 / 共 44 课

第九部分 · 概率、统计与信息

随机变量与分布:给不确定性建模

优化教会了机器"把损失降到最低",但它一直假装数据是确定的。真实数据从骨子里就是随机的。这一课,我们把第 10 课那套讲不确定的道理,升级成机器学习真正用的工具——把每一个随机的量,变成能写公式、能算、能学的对象。

线性回顾
上一课:我们用凸性给优化收了尾——凸问题里局部最优即全局最优,梯度下降"下到底"就保证是最好的;约束优化交给拉格朗日乘子。整个第八部分,我们都在最小化一个确定、已知的损失
留下的问题:可数据本身充满随机与不确定。它是从某个我们看不见的"分布"里抽出来的样本,带着噪声、偶然、抽样误差。优化的语言只会说"把这个数降到最低",根本没法描述"这个量有一半概率是 3、有一半概率是 5"。第 10 课我们为单次随机讲过道理(样本空间 Ω、期望、贝叶斯),但那还停在"掷骰子"的层面。
本课新增:随机变量——把随机结果翻译成数;概率分布——pmf(离散)与 pdf(连续)两种"概率的画像";以及机器学习离不开的几张常见分布脸谱:伯努利、类别、均匀、高斯。这是给"不确定"本身建模的第一块砖。
本课路线
(1) 从样本空间到随机变量:为什么要给随机结果"贴个数字"。(2) 离散分布与概率质量函数 pmf:每个值分多少概率。(3) 连续分布与概率密度函数 pdf:为什么单点概率是 0、面积才是概率。(4) 期望与方差:分布的"重心"和"胖瘦",承接第 10 课。(5) 四张必备脸谱:伯努利、类别、均匀、高斯,各自在机器学习里管什么。(6) 钟形曲线为什么无处不在:中心极限定理一句话。

从样本空间到随机变量:给随机贴上数字

第 10 课我们把"所有可能结果"摆成样本空间 Ω(读作"Omega",全部可能的清单)。但 Ω 里装的常常不是数——"正面 / 反面""红 / 绿 / 蓝""猫 / 狗"。要想算、想优化、想喂给一台机器,我们得先把这些结果翻译成数。这个翻译器,就是随机变量

随机变量
一个随机变量(记作大写 X)就是一个把每个随机结果映射成一个实数的函数。比如掷硬币,规定 X(正面)=1X(反面)=0;掷两骰,可令 X = 两数之和。它本身不是某个固定的值,而是"一台还没摇出结果的、会吐出数字的机器"。摇一次,它落定成某个具体数,叫一次取值

有了这层翻译,所有关于随机的提问都变成关于数的提问:P(X=1)(读作"X 取值为 1 的概率")、P(X ≤ 3)("X 不超过 3 的概率")、𝔼[X]("X 的期望",读作"E X")。机器学习里,一张图片的标签、一次预测的误差、一个像素的亮度——统统被建模成随机变量。把世界的不确定性"数字化",是这一切的起点。

随机变量分两大类,对应两种描述"概率怎么分配"的方式:能一个一个列出取值的叫离散,取值填满一整段连续区间的叫连续。下面分别看它们的"概率画像"。

离散:概率质量函数 pmf

当随机变量只取有限个(或可数个)值——骰子的 1 到 6、类别标签 0/1/2——我们就能给每一个可能取值明明白白地标上一份概率。这张"哪个值分到多少概率"的表,叫概率质量函数(pmf,probability mass function):

p(x) = P(X = x),且 ∑ p(x) = 1

读作"p 在 x 处的值,就是 X 恰好取 x 的概率;把所有取值的概率加起来,必须等于 1"。这里 (读作"求和")那条规矩,正是第 10 课柯尔莫哥洛夫公理"全空间概率为 1"的化身——概率不会凭空多出来,也不会漏掉。每一份概率像一块有重量的"砝码"压在某个取值上,"质量"二字由此而来。公平骰子的 pmf 就是六根等高的柱子,每根 1/6

连续:概率密度函数 pdf,与一个反直觉

可很多量是连续的:一个人的身高、测量的噪声、神经网络某个输出。它能取 1.701.7000001……区间里有无穷多个值,没法给每个值都分一块正概率——真要分,无穷多份正数加起来会爆掉 1。这里冒出一个著名的反直觉:

连续世界的怪事
对一个连续随机变量,任何单独一个值的概率都是 0P(身高 = 1.70000…米) = 0。这不是说"不可能",而是说在连续统里,一个点细得没有任何"宽度",分不到正的概率。可"接近 1.70 米"的人明明存在——那概率藏在哪?

答案是:连续世界里,概率不属于,而属于区间,由一条概率密度函数(pdf,probability density function)f(x)面积来表示:

P(a ≤ X ≤ b) = ∫ₐᵇ f(x) dx

读作"X 落在 a 到 b 之间的概率,等于密度曲线 f 从 a 到 b 围出的面积"。那个 (读作"积分",求曲线下面积)正是第 08 课积分在这里的兑现——没有积分,连续概率根本写不出来。密度 f(x) 本身不是概率(它可以大于 1),而是"概率的浓度":哪里曲线高,哪里附近就更密集地聚着可能性。整条曲线下的总面积,仍然必须是 1

这就解释了那个怪事:单点 P(X=1.70) 是一条"宽度为零的竖线",面积当然是 0;但 P(1.69 ≤ X ≤ 1.71) 是一小片有宽度的面积,就有正概率了。pmf 用柱子的高度表示概率,pdf 用曲线下的面积——这是离散与连续最本质的分野。

期望与方差:分布的重心与胖瘦

一张分布图信息很全,但我们常想用一两个数把它的"性格"捏住。第 10 课已经介绍过第一个:期望 𝔼[X](读作"X 的期望"),它是分布的重心——把每个取值按概率加权平均:

𝔼[X] = ∑ x·p(x) (连续时换成 ∫ x·f(x) dx)

读作"每个取值乘上它的概率(密度),再全部加(积)起来"。期望回答"长期平均落在哪"。但只有重心还不够——两堆数据可以重心一样,却一个挤成一团、一个散得满地都是。刻画这个"散不散"的,是方差 σ²(读作"sigma 平方"):

σ² = Var(X) = 𝔼[(X − μ)²]

读作"方差等于:取值与均值 μ(读作'mu',即期望)之差,平方之后再求期望"。平方有两个用意:让正负偏差不互相抵消,并放大大的偏离。方差越大,分布越"胖"、越分散。它的平方根 σ(读作"sigma",标准差)和取值同一个单位,更好读——它就是"典型偏离均值多远"。期望管中心、方差管散布,这两个数是后面统计与机器学习里反复出现的主角。

四张必备脸谱

机器学习不会每次都从零设计分布,而是从一个"常用脸谱库"里挑一张去套数据。认识这四张脸,你就能读懂半数模型的假设。

① 伯努利分布(Bernoulli)——一次是非题。只有两个结果:成功(1)/失败(0)、正/反、点击/不点击。一个参数 p = 成功概率,于是 P(X=1)=pP(X=0)=1−p。它是二分类的灵魂——逻辑回归输出的就是一个伯努利的 p

② 类别分布(Categorical)——一次多选题。伯努利的多面版:K 个互斥结果(猫/狗/鸟……),每个分一份概率,加起来为 1。它是多分类的灵魂——神经网络末端的 softmax 输出,就是一个类别分布。

③ 均匀分布(Uniform)——一视同仁。在一段区间内每个值同样可能,pdf 是一条水平线。它表示"我对这个量一无所知,先假设处处等可能"——常用来做随机初始化、随机采样的起点。

④ 高斯分布(Gaussian / 正态)——那条钟形曲线。机器学习里出场率最高的一张脸,由两个参数完全决定:均值 μ(钟的中心位置)和方差 σ²(钟的胖瘦),记作 N(μ, σ²)(读作"均值 mu、方差 sigma 平方的正态分布")。它的密度是那条对称、中间高两边迅速衰减的钟形曲线。噪声、测量误差、初始化权重、扩散模型里的隐变量——到处都假设成高斯。下一课整堂都在讲它的多维版本。

从分布采样:直方图如何逼近它的"真身"
选一张脸谱,按按钮往里"扔样本"。每个样本是从这个分布里随机摇出来的;柱状直方图统计样本落在各处的频率。样本越多,直方图越贴近那条理论曲线(红色)——这正是第 10 课大数定律在分布上的回响:随机的样本,群体上服从它背后的分布。同时看样本均值如何收敛到理论均值。
样本数
0
样本均值
理论均值

钟形曲线为什么无处不在

你也许会奇怪:凭什么高斯分布的脸出现得那么频繁?身高、误差、噪声、考试分数,怎么都凑成那条同一形状的钟?这不是巧合,背后有一条概率论里最深刻、最有"宇宙感"的定理——中心极限定理(CLT):

中心极限定理(一句话)
大量相互独立的小随机因素加起来,无论每个小因素自己长什么样(甚至很奇形怪状),它们的总和(或平均)的分布都会趋近一条高斯钟形曲线

一个人的身高,是无数基因与营养的小影响叠加;一次测量误差,是无数微小扰动的总和。凡是"很多独立小因素相加"造出来的量,结果几乎必然呈高斯——这就是钟形曲线统治世界的原因。它也解释了为什么机器学习默认把噪声当成高斯:噪声往往正是无数看不清的小因素之和。中心极限定理把"无数随机的细节"煮成了"一条确定的曲线",是随机世界里少有的、几乎免费的秩序。

常见误解

一句话带走
随机变量把随机结果翻译成数,分布给这些数分配概率——离散用 pmf(柱子高度是概率),连续用 pdf(曲线下面积是概率,单点概率为 0)。期望是重心、方差是胖瘦;伯努利/类别/均匀/高斯四张脸覆盖了机器学习的大半假设,而中心极限定理解释了高斯钟形为何无处不在。不确定性,从此可以写公式、能计算。
下一步
我们学会了给一个随机变量建模。可现实里的不确定量从不孤立——一张图片的相邻像素、一只股票和大盘、身高与体重,它们彼此相关:知道了一个,就改变了对另一个的预期。怎么同时刻画好几个随机变量、并精确地描述它们之间"一起怎么动"?这就要把方差升级成一整张协方差矩阵,把钟形曲线升级成多元高斯。更妙的是,那张协方差矩阵的几何形状是一个椭圆,而它的主轴恰恰就是第 22 课的特征向量、第 23 课的主成分——相关性的几何,和我们早就埋下的线性代数,将在这里重逢。→ 第 30 课《多元高斯与协方差:相关性的几何》。