跳转至

一维随机变量及其分布

在上一节中,我们用集合和事件来描述随机现象.但在物理学中,我们关心的往往是具体的数值——分子的速率是多少?测量误差有多大?为了用微积分等数学工具来分析这些问题,我们需要将随机事件「映射」为实数.这就是随机变量的概念.

学习目标

读完本页后,你应该能够:

  • 理解随机变量的概念,区分离散随机变量和连续随机变量
  • 掌握分布函数和概率密度函数的定义与性质
  • 熟悉均匀分布、指数分布和正态分布(高斯分布)
  • 计算高斯积分,理解其在物理推导中的重要性

随机变量

随机变量

随机变量(Random Variable)是一个函数 X:Ω→RX: Omega -> R,它将样本空间中的每个结果映射到一个实数.

  • 离散随机变量:只能取有限个或可数个值.例如掷骰子的点数 X∈{1,2,3,4,5,6}X in 1,2,3,4,5,6.
  • 连续随机变量:可以取某个区间内的任意值.例如分子的速率 X∈[0,+∞)X in [0, +inf).

引入随机变量后,概率问题就转化为对实数(或实数区间)的概率计算,从而可以使用微积分工具.

分布函数

分布函数

随机变量 XX 的 分布函数(Cumulative Distribution Function, CDF)定义为:

F(x)=P(X≤x)F(x) = P(X <= x)

即 XX 取值不超过 xx 的概率.

分布函数具有以下基本性质:

  1. 单调不减:若 x1<x2x_1 < x_2,则 F(x1)≤F(x2)F(x_1) <= F(x_2)
  2. 右连续:limε→0+F(x+ε)=F(x)lim_varepsilon -> 0^+ F(x + varepsilon) = F(x)
  3. 边界值:F(−∞)=0F(-inf) = 0,F(+∞)=1F(+inf) = 1

用分布函数计算概率:

P(a<X≤b)=F(b)−F(a)P(a < X <= b) = F(b) - F(a)

概率密度函数

对于连续随机变量,分布函数 F(x)F(x) 通常是连续可导的.我们定义其导数为 概率密度函数(Probability Density Function, PDF):

概率密度函数 f(x)=dF(x)dx=F′(x)f(x) = dF(x)dx = F'(x)

反过来,分布函数是概率密度函数的积分:

F(x)=∫−∞xf(t)dtF(x) = int_-inf^x f(t) ,dt

概率密度函数的两条基本性质:

  1. 非负性:f(x)≥0f(x) >= 0(因为 F(x)F(x) 单调不减)
  2. 归一化:
∫−∞+∞f(x)dx=F(+∞)−F(−∞)=1int_-inf^+inf f(x) ,dx = F(+inf) - F(-inf) = 1

用概率密度函数计算概率:

P(a≤X≤b)=∫abf(x)dxP(a <= X <= b) = int_a^b f(x) ,dx
概率密度不是概率

f(x)f(x) 本身 不是概率——它可以大于 1.f(x)f(x) 表示的是概率在 xx 轴上的分布 密度.概率是密度在某个区间上的积分(即曲线下面积).类比:线密度 ρ(x)rho(x) 不是质量,∫ρ(x)dxint rho(x) ,dx 才是质量.

常见一维分布

均匀分布

在区间 [a,b][a, b] 内,取任何值的概率密度都相等.

f(x)={1b−a,a≤x≤b0,其他f(x) = cases 1b-a, & a <= x <= b 0, & 其他 cases

记作 X∼U(a,b)X ~ U(a, b).均匀分布是最简单的连续分布,常用于描述「无偏好」的随机取值.

指数分布

描述独立随机事件发生的时间间隔的概率.若事件以恒定速率 λlambda 发生,则等待时间 TT 服从指数分布:

f(t)={λe−λt,t≥00,t<0f(t) = cases lambda e^-lambda t, & t >= 0 0, & t < 0 cases

指数分布有一个重要性质——无记忆性:已知等待了时间 ss 后,还需再等待时间 tt 的条件概率与 ss 无关:

P(T>s+t|T>s)=P(T>t)P(T > s + t ,| ,T > s) = P(T > t)

在物理学中,放射性衰变的等待时间服从指数分布.

正态分布(高斯分布)

正态分布是自然界和科学研究中最重要的连续分布.

正态分布(高斯分布)

若连续随机变量 XX 的概率密度函数为:

f(x)=1σ2πexp(−(x−μ)22σ2)f(x) = 1sigma2pi ,exp !left(-(x-mu)^22sigma^2right)

则称 XX 服从参数为 μmu(均值)和 σ2sigma^2(方差)的 正态分布(Normal Distribution),记作 X∼N(μ,σ2)X ~ N(mu, sigma^2).

正态分布的特征:

  • 曲线关于 x=μx = mu 对称,呈钟形
  • 在 x=μ±σx = mu +/- sigma 处有拐点
  • μmu 决定分布的中心位置,σsigma 决定分布的宽度
  • 当 μ=0mu = 0,σ=1sigma = 1 时,称为 标准正态分布 N(0,1)N(0, 1)

正态分布之所以在物理学中极为重要,是因为 中心极限定理:大量独立随机变量之和(或均值)近似服从正态分布,无论每个变量本身服从什么分布.

高斯积分

在推导麦克斯韦速率分布律等物理问题中,我们需要反复计算以下形式的积分,统称为 高斯积分:

高斯积分公式

设 α>0alpha > 0,定义:

In(α)=∫0+∞xne−αx2dxI_n(alpha) = int_0^+inf x^n ,e^-alpha x^2 ,dx

常用结果:

I0(α)=∫0+∞e−αx2dx=12παI_0(alpha) = int_0^+inf e^-alpha x^2 ,dx = 12alpha I2(α)=∫0+∞x2e−αx2dx=14πα3I_2(alpha) = int_0^+inf x^2 ,e^-alpha x^2 ,dx = 14alpha^3 I4(α)=∫0+∞x4e−αx2dx=38πα5I_4(alpha) = int_0^+inf x^4 ,e^-alpha x^2 ,dx = 38alpha^5

一般地,对于偶数 n=2kn = 2k:

I2k(α)=(2k−1)!!2k+1πα2k+1I_2k(alpha) = (2k-1)!!2^k+1alpha^2k+1

对于奇数 n=2k+1n = 2k+1:

I2k+1(α)=k!2αk+1I_2k+1(alpha) = k!2alpha^k+1

其中 (2k−1)!!=1⋅3⋅5⋯(2k−1)(2k-1)!! = 1 * 3 * 5 cdots (2k-1) 为双阶乘.

I0I_0 的推导

I0I_0 是最基本的高斯积分,其推导利用了一个巧妙的技巧——将一维积分转化为二维极坐标积分:

I02=(∫−∞+∞e−αx2dx)(∫−∞+∞e−αy2dy)=∫−∞+∞∫−∞+∞e−α(x2+y2)dxdyI_0^2 = left(int_-inf^+inf e^-alpha x^2 ,dxright) !left(int_-inf^+inf e^-alpha y^2 ,dyright) = int_-inf^+inf ! !int_-inf^+inf e^-alpha(x^2+y^2) ,dx ,dy

转化为极坐标 (r,θ)(r, theta):

I02=∫02πdθ∫0+∞e−αr2rdr=2π⋅12α=παI_0^2 = int_0^2pi !dthetaint_0^+inf e^-alpha r^2 ,r ,dr = 2pi * 12alpha = alpha

因此:

I0=∫−∞+∞e−αx2dx=παI_0 = int_-inf^+inf e^-alpha x^2 ,dx = alpha

递推关系

其他高斯积分可以通过对 αalpha 求导得到:

In(α)=−ddαIn−2(α)I_n(alpha) = -ddalpha ,I_n-2(alpha)

这使得所有高斯积分都可以从 I0I_0 出发递推得到.

例题:验证正态分布的归一化

题目:验证 f(x)=1σ2πexp(−(x−μ)22σ2)f(x) = 1sigma2pi ,exp !left(-(x-mu)^22sigma^2right) 满足归一化条件 ∫−∞+∞f(x)dx=1int_-inf^+inf f(x) ,dx = 1.

解答:令 t=x−μσ2t = x - sigma2,则 dx=σ2dtdx = sigma2 ,dt,代入得:

∫−∞+∞f(x)dx=1σ2π∫−∞+∞e−t2σ2dt=22π∫−∞+∞e−t2dt=1π⋅π=1int_-inf^+inf f(x) ,dx = 1sigma2pi int_-inf^+inf e^-t^2 ,sigma2 ,dt = 22pi int_-inf^+inf e^-t^2 ,dt = 1pi * pi = 1

归一化条件得到验证.这里用到了高斯积分 ∫−∞+∞e−t2dt=πint_-inf^+inf e^-t^2 ,dt = pi(即 I0(1)I_0(1) 的全实数轴形式).

高斯积分在物理学中的重要性

高斯积分是统计物理和量子力学中最常用的数学工具之一.在麦克斯韦速率分布律的推导中,我们需要计算形如 ∫0+∞vne−mv2/2kTdvint_0^+inf v^n e^-mv^2/2kT ,dv 的积分,这正是高斯积分 In(m/2kT)I_n(m/2kT).熟练掌握高斯积分的公式和递推方法,是理解统计物理推导的关键.

常见分布的参数汇总

分布记号PDF定义域
均匀分布U(a,b)U(a,b)1b−a1b-a[a,b][a,b]
指数分布Exp(λ)Exp(lambda)λe−λxlambda e^-lambda x[0,+∞)[0,+inf)
正态分布N(μ,σ2)N(mu,sigma^2)1σ2πe−(x−μ)2/2σ21sigma2pie^-(x-mu)^2/2sigma^2(−∞,+∞)(-inf,+inf)

各分布的期望和方差将在 随机变量的数字特征 页面中汇总.

学习衔接



评论