1. 变量与随机变量的区别#

变量是用于表示某个数值的符号。若其取值由确定的条件决定,它可以看作普通的确定性变量。

随机变量是将随机试验的每个样本点映射为实数的函数:

X:Ω→R.X:\Omega\rightarrow\mathbb R.

随机变量的取值在试验前不能确定,但每个可能取值具有相应的概率规律。

例如,掷一次骰子时,“点数”可以定义为随机变量 XX,其可能取值为 1,2,…,61,2,\ldots,6。

2. 古典概型与几何概型#

古典概型适用于样本空间中只有有限个等可能基本事件的情形。若事件 AA 包含 mm 个基本事件,样本空间包含 nn 个基本事件,则:

P(A)=mn.P(A)=\frac{m}{n}.

例如,掷一枚均匀骰子,出现偶数的概率为 3/6=1/23/6=1/2。

几何概型适用于在某个几何区域内均匀随机取点的情形。概率由长度、面积或体积之比确定:

P(A)=事件 A 对应区域的测度整个样本空间的测度.P(A) = \frac{\text{事件 }A\text{ 对应区域的测度}} {\text{整个样本空间的测度}}.

例如,在长度为 LL 的线段上均匀取一点,落入长度为 ll 的子区间的概率为 l/Ll/L。

几何概型的关键不是“可能结果无限”,而是结果在几何区域内按照测度均匀分布。

3. 联合概率、边缘概率与条件概率#

联合概率表示多个事件同时发生的概率:

P(A∩B).P(A\cap B).

边缘概率表示单个事件发生的概率。它可以由联合概率对其他变量求和或积分得到。例如:

P(A)=∑bP(A,B=b).P(A)=\sum_b P(A,B=b).

对连续随机变量:

fX(x)=∫−∞+∞fX,Y(x,y) dy.f_X(x)=\int_{-\infty}^{+\infty}f_{X,Y}(x,y)\,dy.

条件概率表示在事件 BB 已经发生的条件下,事件 AA 发生的概率:

P(A∣B)=P(A∩B)P(B),P(B)>0.P(A\mid B)=\frac{P(A\cap B)}{P(B)},\qquad P(B)>0.

由此可得乘法公式:

P(A∩B)=P(B)P(A∣B)=P(A)P(B∣A).P(A\cap B)=P(B)P(A\mid B)=P(A)P(B\mid A).

4. 全概率公式与贝叶斯公式#

设事件 B1,B2,…,BnB_1,B_2,\ldots,B_n 两两互斥,且构成样本空间的一个完备划分:

Bi∩Bj=∅(i≠j),⋃i=1nBi=Ω.B_i\cap B_j=\varnothing\quad(i\ne j),\qquad \bigcup_{i=1}^n B_i=\Omega.

全概率公式用于“由原因推结果”。若要计算结果 AA 发生的概率,可以将所有可能原因分别考虑:

P(A)=∑i=1nP(Bi)P(A∣Bi).P(A)=\sum_{i=1}^n P(B_i)P(A\mid B_i).

贝叶斯公式用于“由结果推原因”。在已知 AA 已经发生时,原因 BiB_i 发生的后验概率为:

P(Bi∣A)=P(A∣Bi)P(Bi)∑j=1nP(A∣Bj)P(Bj).P(B_i\mid A) = \frac{P(A\mid B_i)P(B_i)} {\sum_{j=1}^n P(A\mid B_j)P(B_j)}.

其中,P(Bi)P(B_i) 是先验概率,P(A∣Bi)P(A\mid B_i) 是似然,P(Bi∣A)P(B_i\mid A) 是后验概率。

5. 先验概率、似然与后验概率#

在贝叶斯推断中,设 θ\theta 为未知参数,DD 为观测数据。

  • 先验概率 P(θ)P(\theta):观察数据之前,对参数 θ\theta 的已有认识。
  • 似然 P(D∣θ)P(D\mid\theta):在参数为 θ\theta 时,观测到数据 DD 的可能性。
  • 后验概率 P(θ∣D)P(\theta\mid D):观察数据之后,对参数概率分布的更新。

贝叶斯公式为:

P(θ∣D)=P(D∣θ)P(θ)P(D)∝P(D∣θ)P(θ).P(\theta\mid D) = \frac{P(D\mid\theta)P(\theta)}{P(D)} \propto P(D\mid\theta)P(\theta).

先验概率不是由当前数据计算得到的;后验概率则综合了先验知识与观测数据。

6. 常见概率分布#

分布概率质量函数或密度函数期望方差
伯努利分布 B(1,p)B(1,p)P(X=k)=pk(1−p)1−k, k=0,1P(X=k)=p^k(1-p)^{1-k},\ k=0,1ppp(1−p)p(1-p)
二项分布 B(n,p)B(n,p)P(X=k)=(nk)pk(1−p)n−kP(X=k)=\binom nkp^k(1-p)^{n-k}npnpnp(1−p)np(1-p)
泊松分布 Pois⁡(λ)\operatorname{Pois}(\lambda)P(X=k)=λke−λk!P(X=k)=\dfrac{\lambda^ke^{-\lambda}}{k!}λ\lambdaλ\lambda
超几何分布 H(N,M,n)H(N,M,n)P(X=k)=(Mk)(N−Mn−k)(Nn)P(X=k)=\dfrac{\binom Mk\binom{N-M}{n-k}}{\binom Nn}nMN\dfrac{nM}{N}nMN(1−MN)N−nN−1n\dfrac MN\left(1-\dfrac MN\right)\dfrac{N-n}{N-1}
几何分布 Geo⁡(p)\operatorname{Geo}(p)P(X=k)=(1−p)k−1p, k≥1P(X=k)=(1-p)^{k-1}p,\ k\ge11p\dfrac1p1−pp2\dfrac{1-p}{p^2}
负二项分布 NB⁡(r,p)\operatorname{NB}(r,p)P(X=k)=(k−1r−1)pr(1−p)k−r, k≥rP(X=k)=\binom{k-1}{r-1}p^r(1-p)^{k-r},\ k\ge rrp\dfrac rpr(1−p)p2\dfrac{r(1-p)}{p^2}
均匀分布 U(a,b)U(a,b)f(x)=1b−a, a<x<bf(x)=\dfrac1{b-a},\ a<x<ba+b2\dfrac{a+b}{2}(b−a)212\dfrac{(b-a)^2}{12}
正态分布 N(μ,σ2)N(\mu,\sigma^2)f(x)=12πσe−(x−μ)22σ2f(x)=\dfrac1{\sqrt{2\pi}\sigma}e^{-\frac{(x-\mu)^2}{2\sigma^2}}μ\muσ2\sigma^2
指数分布 Exp⁡(λ)\operatorname{Exp}(\lambda)f(x)=λe−λx, x≥0f(x)=\lambda e^{-\lambda x},\ x\ge01λ\dfrac1\lambda1λ2\dfrac1{\lambda^2}

补充说明:

  • 伯努利分布也称 0–1 分布,是两点分布的一种特殊形式。
  • 超几何分布描述不放回抽样;二项分布通常描述独立重复试验或近似放回抽样。
  • 当 nn 很大、pp 很小且 λ=np\lambda=np 适中时:
B(n,p)≈Pois⁡(np).B(n,p)\approx\operatorname{Pois}(np).
  • 泊松分布常用于描述固定时间或空间范围内随机事件的发生次数。
  • 指数分布常用于描述泊松过程中相邻事件之间的等待时间。

7. 概率密度函数#

连续型随机变量 XX 的概率密度函数记为 fX(x)f_X(x),满足:

fX(x)≥0,∫−∞+∞fX(x) dx=1.f_X(x)\ge0, \qquad \int_{-\infty}^{+\infty}f_X(x)\,dx=1.

区间概率由密度函数积分得到:

P(a<X≤b)=∫abfX(x) dx.P(a<X\le b)=\int_a^bf_X(x)\,dx.

累积分布函数为:

FX(x)=P(X≤x)=∫−∞xfX(t) dt.F_X(x)=P(X\le x) = \int_{-\infty}^xf_X(t)\,dt.

若 FXF_X 可导,则:

fX(x)=FX′(x).f_X(x)=F_X'(x).

需要注意:

  • fX(x)f_X(x) 本身不是点 xx 处的概率。
  • 对连续型随机变量,P(X=x)=0P(X=x)=0。
  • 概率密度值可以大于 11,只要总积分等于 11。

8. 期望与方差#

期望描述随机变量取值的平均水平或集中位置。

离散型随机变量:

E(X)=∑xxP(X=x).E(X)=\sum_x xP(X=x).

连续型随机变量:

E(X)=∫−∞+∞xfX(x) dx.E(X)=\int_{-\infty}^{+\infty}xf_X(x)\,dx.

期望具有线性性质:

E(aX+bY+c)=aE(X)+bE(Y)+c.E(aX+bY+c)=aE(X)+bE(Y)+c.

该性质不要求 XX 与 YY 独立。

方差描述随机变量围绕期望的离散程度:

Var⁡(X)=E[(X−E(X))2]=E(X2)−[E(X)]2.\operatorname{Var}(X) = E\left[(X-E(X))^2\right] = E(X^2)-[E(X)]^2.

常用性质:

Var⁡(aX+b)=a2Var⁡(X).\operatorname{Var}(aX+b)=a^2\operatorname{Var}(X).

标准差为:

σX=Var⁡(X).\sigma_X=\sqrt{\operatorname{Var}(X)}.

9. 协方差、协方差矩阵与相关系数#

协方差描述两个随机变量共同变化的方向:

Cov⁡(X,Y)=E[(X−E(X))(Y−E(Y))]=E(XY)−E(X)E(Y).\operatorname{Cov}(X,Y) = E\bigl[(X-E(X))(Y-E(Y))\bigr] = E(XY)-E(X)E(Y).
  • Cov⁡(X,Y)>0\operatorname{Cov}(X,Y)>0:二者总体上同向变化。
  • Cov⁡(X,Y)<0\operatorname{Cov}(X,Y)<0:二者总体上反向变化。
  • Cov⁡(X,Y)=0\operatorname{Cov}(X,Y)=0:二者不存在线性相关关系,但不一定独立。

两个随机变量之和的方差为:

Var⁡(X+Y)=Var⁡(X)+Var⁡(Y)+2Cov⁡(X,Y).\operatorname{Var}(X+Y) = \operatorname{Var}(X) + \operatorname{Var}(Y) + 2\operatorname{Cov}(X,Y).

若 XX 与 YY 独立,则协方差为 00,因此:

Var⁡(X+Y)=Var⁡(X)+Var⁡(Y).\operatorname{Var}(X+Y) = \operatorname{Var}(X)+\operatorname{Var}(Y).

对于随机向量:

X=[X1,X2,…,Xn]T,\mathbf X=[X_1,X_2,\ldots,X_n]^T,

其协方差矩阵为:

Σ=(σij)n×n,σij=Cov⁡(Xi,Xj).\Sigma=(\sigma_{ij})_{n\times n}, \qquad \sigma_{ij}=\operatorname{Cov}(X_i,X_j).

协方差矩阵是对称半正定矩阵,对角线元素是各变量的方差。

相关系数是标准化后的协方差:

ρX,Y=Cov⁡(X,Y)Var⁡(X)Var⁡(Y).\rho_{X,Y} = \frac{\operatorname{Cov}(X,Y)} {\sqrt{\operatorname{Var}(X)}\sqrt{\operatorname{Var}(Y)}}.

其取值范围为:

−1≤ρX,Y≤1.-1\le\rho_{X,Y}\le1.

10. 独立与不相关的区别#

若随机变量 XX 和 YY 独立,则对任意合适的集合 A,BA,B:

P(X∈A,Y∈B)=P(X∈A)P(Y∈B).P(X\in A,Y\in B) = P(X\in A)P(Y\in B).

在二阶矩存在时,独立能够推出:

Cov⁡(X,Y)=0.\operatorname{Cov}(X,Y)=0.

因此:

独立⟹不相关.\text{独立}\Longrightarrow\text{不相关}.

反之通常不成立。不相关仅表示二者之间没有线性关系,仍可能存在非线性依赖。

例如,设 X∼U(−1,1)X\sim U(-1,1),Y=X2Y=X^2。此时 Cov⁡(X,Y)=0\operatorname{Cov}(X,Y)=0,但 YY 完全由 XX 决定,因此二者并不独立。

对于联合正态分布的随机变量,不相关可以推出独立。

11. 独立与互斥的关系#

事件 AA 与 BB 独立是指:

P(A∩B)=P(A)P(B).P(A\cap B)=P(A)P(B).

事件 AA 与 BB 互斥是指:

A∩B=∅,A\cap B=\varnothing,

因此:

P(A∩B)=0.P(A\cap B)=0.

若 P(A)>0P(A)>0 且 P(B)>0P(B)>0,则互斥事件不可能独立,因为:

P(A)P(B)>0≠P(A∩B).P(A)P(B)>0\ne P(A\cap B).

直观上:

  • 互斥表示两个事件不能同时发生。
  • 独立表示一个事件是否发生不会改变另一个事件的概率。

12. 高斯随机变量的线性组合与乘积#

若 Xi∼N(μi,σi2)X_i\sim N(\mu_i,\sigma_i^2) 相互独立,则它们的线性组合仍服从正态分布:

Y=b+∑i=1naiXi∼N(b+∑i=1naiμi, ∑i=1nai2σi2).Y=b+\sum_{i=1}^n a_iX_i \sim N\left( b+\sum_{i=1}^na_i\mu_i,\, \sum_{i=1}^na_i^2\sigma_i^2 \right).

特别地:

X+Y∼N(μX+μY,σX2+σY2)X+Y\sim N(\mu_X+\mu_Y,\sigma_X^2+\sigma_Y^2)

成立的前提是 XX 与 YY 独立。

常数倍和平移也保持正态性:

X∼N(μ,σ2)⟹aX+b∼N(aμ+b,a2σ2).X\sim N(\mu,\sigma^2) \quad\Longrightarrow\quad aX+b\sim N(a\mu+b,a^2\sigma^2).

但两个高斯随机变量的乘积 XYXY 一般不服从高斯分布。只有线性组合具有稳定的高斯形式。

13. 大数定律#

大数定律说明:当重复试验次数足够多时,样本平均值会以较大概率接近期望。它建立了“理论期望”和“实际平均结果”之间的联系。

切比雪夫大数定律

若 X1,X2,…,XnX_1,X_2,\ldots,X_n 相互独立,期望和方差均存在,并且方差一致有界,则:

1n∑i=1nXi−1n∑i=1nE(Xi)→P0.\frac{1}{n}\sum_{i=1}^n X_i - \frac{1}{n}\sum_{i=1}^n E(X_i) \xrightarrow{P}0.

伯努利大数定律

若 Xn∼B(n,p)X_n\sim B(n,p) 表示 nn 次独立伯努利试验中成功的次数,则:

Xnn→Pp.\frac{X_n}{n}\xrightarrow{P}p.

也就是说,试验次数足够多时,事件发生的频率依概率收敛于其概率。

辛钦大数定律

若 X1,X2,…X_1,X_2,\ldots 独立同分布,且 E(Xi)=μE(X_i)=\mu 存在,则:

X‾n=1n∑i=1nXi→Pμ.\overline X_n=\frac{1}{n}\sum_{i=1}^nX_i \xrightarrow{P}\mu.

大数定律描述的是样本均值的收敛性,并不直接说明误差服从什么分布。

14. 中心极限定理#

中心极限定理说明:大量独立随机变量之和经过标准化后,通常近似服从标准正态分布。

若 X1,X2,…,XnX_1,X_2,\ldots,X_n 独立同分布,且

E(Xi)=μ,Var⁡(Xi)=σ2>0,E(X_i)=\mu,\qquad \operatorname{Var}(X_i)=\sigma^2>0,

则:

∑i=1nXi−nμσn=X‾n−μσ/n→dN(0,1).\frac{\sum_{i=1}^nX_i-n\mu}{\sigma\sqrt n} = \frac{\overline X_n-\mu}{\sigma/\sqrt n} \xrightarrow{d}N(0,1).

因此,当 nn 较大时:

∑i=1nXi≈N(nμ,nσ2),X‾n≈N(μ,σ2n).\sum_{i=1}^nX_i\approx N(n\mu,n\sigma^2), \qquad \overline X_n\approx N\left(\mu,\frac{\sigma^2}{n}\right).

棣莫弗–拉普拉斯定理

若 Xn∼B(n,p)X_n\sim B(n,p),则当 nn 较大时:

Xn−npnp(1−p)→dN(0,1).\frac{X_n-np}{\sqrt{np(1-p)}}\xrightarrow{d}N(0,1).

即:

Xn≈N(np,np(1−p)).X_n\approx N\bigl(np,np(1-p)\bigr).

使用正态分布近似离散分布时,通常可以加入 0.50.5 连续性修正。

例如,每发炮弹命中数的期望为 22、方差为 1.691.69,独立发射 100100 次,则总命中数 SS 近似满足:

S≈N(200,169).S\approx N(200,169).

于是:

P(180≤S≤220)≈Φ(220.5−20013)−Φ(179.5−20013)≈0.885.P(180\le S\le220) \approx \Phi\left(\frac{220.5-200}{13}\right) - \Phi\left(\frac{179.5-200}{13}\right) \approx0.885.

15. 极大似然估计#

极大似然估计用于根据样本估计未知参数。

设样本 x1,x2,…,xnx_1,x_2,\ldots,x_n 来自参数为 θ\theta 的分布,则似然函数为:

L(θ)=∏i=1nf(xi;θ).L(\theta) = \prod_{i=1}^n f(x_i;\theta).

极大似然估计量定义为:

θ^MLE=arg⁡max⁡θL(θ).\hat\theta_{\mathrm{MLE}} = \arg\max_\theta L(\theta).

通常使用对数似然:

ℓ(θ)=log⁡L(θ)=∑i=1nlog⁡f(xi;θ),\ell(\theta) = \log L(\theta) = \sum_{i=1}^n\log f(x_i;\theta),

因为对数可以将乘积变为求和,而且不会改变最大值的位置。

一般步骤:

  1. 写出似然函数 L(θ)L(\theta)。
  2. 取对数,得到 ℓ(θ)\ell(\theta)。
  3. 对参数求导并令导数为 00。
  4. 求出候选解,并检查其是否为最大值以及是否满足参数范围。

极大似然估计得到的是基于有限样本的估计值。在适当的正则条件下,随着样本量增大,它通常具有一致性和渐近正态性,但这些性质并非在所有模型中自动成立。

16. 假设检验#

假设检验是在总体分布或参数未知时,根据样本信息判断某个假设是否具有足够证据被拒绝。

  • 原假设 H0H_0:默认成立、准备接受检验的假设。
  • 备择假设 H1H_1:当证据足够强时所支持的假设。

两类错误:

  • 第一类错误:H0H_0 为真,却拒绝 H0H_0,概率为 α\alpha。
  • 第二类错误:H0H_0 为假,却未拒绝 H0H_0,概率为 β\beta。

检验功效为:

1−β,1-\beta,

表示原假设确实为假时,正确拒绝它的概率。

一般步骤:

  1. 提出 H0H_0 和 H1H_1。
  2. 选择检验统计量及其在 H0H_0 下的分布。
  3. 给定显著性水平 α\alpha。
  4. 根据样本计算检验统计量或 pp 值。
  5. 若 p≤αp\le\alpha,则拒绝 H0H_0;否则只能说“没有足够证据拒绝 H0H_0”,不能直接证明 H0H_0 正确。

在样本量固定时,α\alpha 与 β\beta 往往存在权衡;增加样本量通常可以同时改善两类错误。

17. 机器学习为什么使用概率#

机器学习处理的数据通常具有噪声、不确定性和随机性,因此概率论是描述数据生成规律的重要工具。

主要作用包括:

  • 用概率分布描述数据、标签、噪声和模型参数。
  • 用条件概率表达预测问题,如 P(y∣x)P(y\mid x)。
  • 用似然函数和后验概率进行参数学习与推断。
  • 用期望风险描述模型的平均预测误差。
  • 用概率输出表达模型的不确定性。
  • 许多常见损失函数可以解释为负对数似然。

例如,二分类交叉熵可以写成伯努利模型的负对数似然;均方误差则可由高斯噪声假设推导得到。

18. 用不均匀硬币产生均匀随机数#

设一枚硬币产生 00 的概率为 pp,产生 11 的概率为 1−p1-p,且每次投掷相互独立。

连续投掷两次:

P(01)=p(1−p),P(10)=(1−p)p.P(01)=p(1-p), \qquad P(10)=(1-p)p.

因此,0101 和 1010 等概率。冯·诺依曼提取法规定:

  • 出现 0101,输出 00。
  • 出现 1010,输出 11。
  • 出现 0000 或 1111,丢弃并重新投掷。

这样可以得到无偏的 00–11 随机比特,不需要知道 pp 的具体值。

要产生 1,2,…,m1,2,\ldots,m 上的均匀随机整数,可以先利用上述方法生成无偏二进制位,再采用拒绝采样:

  1. 取 k=⌈log⁡2m⌉k=\lceil\log_2m\rceil。
  2. 用 kk 个无偏比特生成 00 到 2k−12^k-1 的整数 RR。
  3. 若 R<mR<m,输出 R+1R+1;否则丢弃并重新生成。

19. 伪随机数与均匀分布#

计算机中的伪随机数由确定性算法生成。给定相同的初始种子,生成的序列完全相同,但在统计性质上可以近似随机序列。

一种经典方法是线性同余生成器:

Xn=(aXn−1+c) mod m,X_n=(aX_{n-1}+c)\bmod m,

其中:

  • mm 为模数;
  • aa 为乘数;
  • cc 为增量;
  • X0X_0 为初始种子。

将整数序列归一化可得到 [0,1)[0,1) 上的伪均匀随机数:

Un=Xnm.U_n=\frac{X_n}{m}.

参数选择会影响随机序列的周期和统计质量。现代程序通常使用质量更高的伪随机数生成算法;密码学场景则必须使用密码学安全伪随机数生成器。

20. 由均匀分布生成高斯分布#

Box–Muller 变换可以将两个相互独立的均匀随机变量转换为两个相互独立的标准正态随机变量。

设:

U1,U2∼i.i.d.U(0,1),U_1,U_2\overset{\text{i.i.d.}}{\sim}U(0,1),

则:

Z0=−2ln⁡U1cos⁡(2πU2),Z_0 = \sqrt{-2\ln U_1}\cos(2\pi U_2), Z1=−2ln⁡U1sin⁡(2πU2).Z_1 = \sqrt{-2\ln U_1}\sin(2\pi U_2).

此时:

Z0,Z1∼i.i.d.N(0,1).Z_0,Z_1\overset{\text{i.i.d.}}{\sim}N(0,1).

若需要生成一般正态分布 N(μ,σ2)N(\mu,\sigma^2),可进行线性变换:

X=μ+σZ,Z∼N(0,1).X=\mu+\sigma Z, \qquad Z\sim N(0,1).