分类课程智能体AI
文章
订阅
分类课程AI导师
文章
价格
课程进度
2 / 12
上一节绪论与热的基本概念下一节温度与玻尔兹曼因子
自在学

© 2025 - 2026 株洲市自在学教育科技有限公司 版权所有

公网安备湘公网安备43020302000292号 | 湘ICP备2025148919号-1

关于我们隐私政策使用条款

© 2025 - 2026 株洲市自在学教育科技有限公司 版权所有

公网安备湘公网安备43020302000292号湘ICP备2025148919号-1

物理热力学与统计物理概率论基础

概率论基础

热力学与统计物理的核心,在于用概率的语言描述大量粒子的集体行为。当我们面对一个具体的气体分子时,无法准确预测它在某一时刻的速度、位置和动能,因为微观世界充满了随机性和不确定性。然而,若我们观察的是数以亿万计的分子整体,它们的速度分布、能量分布却表现出惊人的规律性和可预见性。这种从“随机个体”到“规律整体”的转变,是统计物理最本质、最奇妙的特点之一。

正是由于个体的不确定、系统的巨大数目,以及复杂相互作用的影响,单靠经典力学无法直接刻画宏观物理量的表现。我们需要概率论和统计方法,来描述粒子群体的行为、预测物理量的平均值与涨落规律。例如,通过概率分布函数,可以得知某一速度区间内分子的比例,进而推算出气体的压强、温度等宏观性质。换言之,统计物理的强大之处,在于用抽象的概率工具,将微观世界的无序和复杂,转化为宏观世界的简洁与确定。

理解概率论的基本原理,是打开统计物理大门的钥匙,也是后续玻尔兹曼分布、熵、涨落乃至相变等核心内容的理论基础。


随机现象与概率

掷一枚骰子,事先无法知道结果是几;但掷出足够多次,每个点数出现的频率都趋近于 1/61/61/6。这种在个体层面无法预测、在大量重复中却涌现出稳定规律的现象,称为随机现象。

概率是对随机事件发生可能性的定量描述。对某一事件 AAA,其概率 P(A)P(A)P(A) 满足:

0≤P(A)≤10 \leq P(A) \leq 10≤P(A)≤1

P(A)=0P(A) = 0P(A)=0 表示事件不可能发生,P(A)=1P(A) = 1P(A)=1 表示事件必然发生。若事件 AAA 与事件 BBB 互斥(不能同时发生),则:

P(A∪B)=P(A)+P(B)P(A \cup B) = P(A) + P(B)P(A∪B)=P(A)+P(B)

事件 AAA 的对立事件(即 AAA 不发生)记作 Aˉ\bar{A}Aˉ,满足:

P(Aˉ)=1−P(A)P(\bar{A}) = 1 - P(A)P(Aˉ)=1−P(A)

几个典型随机实验的概率结构对比如下:

welearn-96798044.png

概率反映的是大量重复实验中事件出现的频率极限,而不是单次实验的结论。单次掷骰子得到“6”完全正常,但掷 6000 次后“6”出现的次数应非常接近 1000。


离散概率分布

当随机变量 XXX 只能取有限个或可数个离散值时,描述它的工具是概率分布列。对每个可能取值 xix_ixi​,给出对应的概率 P(X=xi)=piP(X = x_i) = p_iP(X=xi​)=pi​,所有概率之和为 111:

∑ipi=1\sum_i p_i = 1i∑​pi​=1

以掷一枚均匀骰子为例,点数 XXX 的概率分布如下:

XXX111222333444555666
PPP16\dfrac{1}{6}61​16\dfrac{1}{6}61​16\dfrac{1}{6}61​16\dfrac{1}{6}61​16\dfrac{1}{6}61​16\dfrac{1}{6}61​

这是最简单的均匀分布——每个结果等概率出现。实际的物理问题中,概率分布往往不均匀。以一枚重心偏移的“不均匀硬币”为例,若正面出现的概率 p=0.6p = 0.6p=0.6,反面的概率 1−p=0.41-p = 0.41−p=0.4:

结果正面反面
概率0.60.60.60.40.40.4

在统计物理中,能量低的微观状态被占据的概率更高——这与后续的玻尔兹曼因子 e−E/kBTe^{-E/k_BT}e−E/kB​T 直接对应,是热力学统计描述的核心。


期望值与方差

期望值(均值)⟨X⟩\langle X \rangle⟨X⟩ 是随机变量取值的概率加权平均,代表大量实验的平均结果:

⟨X⟩=∑ixi pi\langle X \rangle = \sum_i x_i \, p_i⟨X⟩=i∑​xi​pi​

对掷均匀骰子,期望值为:

⟨X⟩=1×16+2×16+⋯+6×16=1+2+3+4+5+66=216=3.5\langle X \rangle = 1 \times \frac{1}{6} + 2 \times \frac{1}{6} + \cdots + 6 \times \frac{1}{6} = \frac{1+2+3+4+5+6}{6} = \frac{21}{6} = 3.5⟨X⟩=1×61​+2×61​+⋯+6×61​=61+2+3+4+5+6​=621​=3.5

注意 3.53.53.5 不是骰子能出现的点数,它是大量实验结果趋向的平均值。

期望值只告诉了分布的“中心位置”,还需要一个量描述数据的分散程度。方差 Var(X)\text{Var}(X)Var(X) 定义为偏差平方的期望:

Var(X)=⟨(X−⟨X⟩)2⟩=⟨X2⟩−⟨X⟩2\text{Var}(X) = \langle (X - \langle X \rangle)^2 \rangle = \langle X^2 \rangle - \langle X \rangle^2Var(X)=⟨(X−⟨X⟩)2⟩=⟨X2⟩−⟨X⟩2

标准差 σ\sigmaσ 是方差的平方根,与 XXX 具有相同的量纲,物理意义更为直观:

σ=Var(X)\sigma = \sqrt{\text{Var}(X)}σ=Var(X)​

对掷骰子,先计算 ⟨X2⟩\langle X^2 \rangle⟨X2⟩:

⟨X2⟩=12+22+32+42+52+626=1+4+9+16+25+366=916≈15.17\langle X^2 \rangle = \frac{1^2+2^2+3^2+4^2+5^2+6^2}{6} = \frac{1+4+9+16+25+36}{6} = \frac{91}{6} \approx 15.17⟨X2⟩=612+22+32+42+52+62​=61+4+9+16+25+36​=691​≈15.17 Var(X)=15.17−(3.5)2=15.17−12.25≈2.92,σ≈1.71\text{Var}(X) = 15.17 - (3.5)^2 = 15.17 - 12.25 \approx 2.92, \quad \sigma \approx 1.71Var(X)=15.17−(3.5)2=15.17−12.25≈2.92,σ≈1.71

期望值描述分布的“中心位置”,方差(或标准差)描述分布的“宽窄”。在统计物理中,期望值对应可测量的宏观量,标准差对应该量的涨落幅度。

三个统计量的物理对照:

welearn-28358500.png


连续概率分布

许多物理量(如分子速度、粒子位置)可以取连续的实数值,需要用概率密度函数 f(x)f(x)f(x) 来描述。f(x)f(x)f(x) 的含义是:随机变量落在 [x, x+dx][x,\, x+\mathrm{d}x][x,x+dx] 内的概率为 f(x) dxf(x)\,\mathrm{d}xf(x)dx。对有限区间积分得到概率:

P(a≤X≤b)=∫abf(x) dxP(a \leq X \leq b) = \int_a^b f(x) \, \mathrm{d}xP(a≤X≤b)=∫ab​f(x)dx

归一化条件要求:

∫−∞+∞f(x) dx=1\int_{-\infty}^{+\infty} f(x) \, \mathrm{d}x = 1∫−∞+∞​f(x)dx=1

连续分布的期望值与方差计算,将求和替换为积分:

⟨X⟩=∫−∞+∞x f(x) dx\langle X \rangle = \int_{-\infty}^{+\infty} x \, f(x) \, \mathrm{d}x⟨X⟩=∫−∞+∞​xf(x)dx Var(X)=∫−∞+∞(x−⟨X⟩)2 f(x) dx\text{Var}(X) = \int_{-\infty}^{+\infty} (x - \langle X \rangle)^2 \, f(x) \, \mathrm{d}xVar(X)=∫−∞+∞​(x−⟨X⟩)2f(x)dx

以均匀分布为例:在区间 [0,L][0, L][0,L] 上取 f(x)=1/Lf(x) = 1/Lf(x)=1/L,在区间外 f(x)=0f(x) = 0f(x)=0。直接积分可得 ⟨X⟩=L/2\langle X \rangle = L/2⟨X⟩=L/2(恰好在区间中点),方差 Var(X)=L2/12\text{Var}(X) = L^2/12Var(X)=L2/12。这与气体分子在容器中均匀分布的物理图像完全吻合——每个位置等概率出现,平均位置在容器中心。

welearn-68819360.png

概率密度函数 f(x)f(x)f(x) 本身不是概率,它的量纲是概率除以 XXX 的量纲。只有将 f(x)f(x)f(x) 对一段区间积分,才能得到该区间内的概率。f(x)f(x)f(x) 在某点处的值可以大于 111,但对全空间的积分恒等于 111。


独立随机变量的叠加

设 X1,X2,…,XNX_1, X_2, \ldots, X_NX1​,X2​,…,XN​ 是 NNN 个相互独立的随机变量,各自的期望值和方差分别为 μi\mu_iμi​ 和 σi2\sigma_i^2σi2​。定义它们的总和 SN=X1+X2+⋯+XNS_N = X_1 + X_2 + \cdots + X_NSN​=X1​+X2​+⋯+XN​,则:

⟨SN⟩=∑i=1Nμi\langle S_N \rangle = \sum_{i=1}^N \mu_i⟨SN​⟩=i=1∑N​μi​ Var(SN)=∑i=1Nσi2\text{Var}(S_N) = \sum_{i=1}^N \sigma_i^2Var(SN​)=i=1∑N​σi2​

方差的可加性要求随机变量相互独立,这是关键假设。若所有 XiX_iXi​ 同分布(均值 μ\muμ、方差 σ2\sigma^2σ2),则:

⟨SN⟩=Nμ,σSN=N σ\langle S_N \rangle = N\mu, \quad \sigma_{S_N} = \sqrt{N}\,\sigma⟨SN​⟩=Nμ,σSN​​=N​σ

相对涨落(标准差与均值之比)为:

σSN⟨SN⟩=N σNμ=σμN\frac{\sigma_{S_N}}{\langle S_N \rangle} = \frac{\sqrt{N}\,\sigma}{N\mu} = \frac{\sigma}{\mu\sqrt{N}}⟨SN​⟩σSN​​​=NμN​σ​=μN​σ​

随着 NNN 增大,相对涨落按 1/N1/\sqrt{N}1/N​ 衰减,绝对涨落在增大,而相对于均值的涨落在缩小。下表展示了这种趋势:

当 N=1023N = 10^{23}N=1023 时,相对涨落约为 10−1210^{-12}10−12,远低于任何仪器的分辨率。宏观量之所以高度稳定,原因不是粒子运动变得“整齐”,而是大量随机量叠加后相对涨落极小——这正是热力学规律得以精确成立的统计学根基。


二项分布

抛一枚硬币 NNN 次,每次正面出现的概率为 ppp,反面的概率为 1−p1-p1−p。NNN 次独立实验中恰好出现 nnn 次正面的概率,由二项分布给出:

P(n; N, p)=(Nn)pn(1−p)N−nP(n;\, N,\, p) = \binom{N}{n} p^n (1-p)^{N-n}P(n;N,p)=(nN​)pn(1−p)N−n

其中组合数 (Nn)=N!n! (N−n)!\dbinom{N}{n} = \dfrac{N!}{n!\,(N-n)!}(nN​)=n!(N−n)!N!​ 计算了从 NNN 次实验中选出 nnn 次“成功”的方案总数。二项分布的期望值与方差:

⟨n⟩=Np,Var(n)=Np(1−p),σn=Np(1−p)\langle n \rangle = Np, \quad \text{Var}(n) = Np(1-p), \quad \sigma_n = \sqrt{Np(1-p)}⟨n⟩=Np,Var(n)=Np(1−p),σn​=Np(1−p)​

以 N=4N = 4N=4、p=1/2p = 1/2p=1/2 为例,列出完整分布:

期望值 ⟨n⟩=4×0.5=2\langle n \rangle = 4 \times 0.5 = 2⟨n⟩=4×0.5=2,标准差 σ=4×0.5×0.5=1\sigma = \sqrt{4 \times 0.5 \times 0.5} = 1σ=4×0.5×0.5​=1。出现 2 次正面的概率最大(37.5%),偏离越远概率越小,分布关于均值对称。

二项分布在统计物理中有直接的物理类比:将容器分为体积相等的左右两半,每个气体分子以概率 p=1/2p = 1/2p=1/2 出现在左半。NNN 个分子中恰好有 nnn 个在左半的概率,精确地由二项分布描述。当 NNN 非常大时,几乎所有分子几乎均等地分布在两侧——偏离均值的概率按指数规律衰减,这正是气体自发均匀分布的统计学根源。

二项分布的物理意义在于:宏观上“理所当然”的均匀分布,在微观上是统计规律的必然结果,而非某种神秘的“斥力”。绝大多数微观状态对应于接近均匀的宏观分布,偏离均匀的状态数量极少——热力学第二定律的根基正在于此。


高斯近似

当 NNN 很大时,直接计算组合数 (Nn)\dbinom{N}{n}(nN​) 极为困难。利用斯特林近似 ln⁡(N!)≈Nln⁡N−N\ln(N!) \approx N\ln N - Nln(N!)≈NlnN−N,可以证明:当 N≫1N \gg 1N≫1 且 nnn 在均值附近时,二项分布收敛到高斯(正态)分布:

P(n)≈12πNp(1−p)exp⁡ ⁣(−(n−Np)22Np(1−p))P(n) \approx \frac{1}{\sqrt{2\pi N p(1-p)}} \exp\!\left(-\frac{(n - Np)^2}{2Np(1-p)}\right)P(n)≈2πNp(1−p)​1​exp(−2Np(1−p)(n−Np)2​)

这是均值为 μ=Np\mu = Npμ=Np、标准差为 σ=Np(1−p)\sigma = \sqrt{Np(1-p)}σ=Np(1−p)​ 的高斯函数,呈现出标志性的“钟形”曲线,关于均值左右对称,在 μ±σ\mu \pm \sigmaμ±σ 处有拐点。

高斯分布的几个重要性质:

welearn-18093723.png

这被称为“68–95–99.7 规则”,在实验误差估计和数据分析中广泛使用。

以 N=100N = 100N=100 次抛硬币(p=0.5p = 0.5p=0.5)为例:均值 μ=50\mu = 50μ=50,标准差 σ=5\sigma = 5σ=5。正面次数落在 [45,55][45, 55][45,55] 的概率约为 68.3%68.3\%68.3%,落在 [40,60][40, 60][40,60] 的概率约为 95.4%95.4\%95.4%。对气体中的粒子分布,N=1023N = 10^{23}N=1023 时,σ=N/2≈1.6×1011\sigma = \sqrt{N}/2 \approx 1.6 \times 10^{11}σ=N​/2≈1.6×1011,而相对涨落 σ/μ≈6.3×10−12\sigma/\mu \approx 6.3 \times 10^{-12}σ/μ≈6.3×10−12,远低于任何实验仪器的分辨能力——宏观上气体密度看起来完全均匀。

从二项分布到高斯分布的过渡,是统计物理中最重要的数学结论之一。无论原始分布的形状如何,大量独立随机变量的总和总趋向于高斯分布——这就是中心极限定理的核心思想。热力学的稳定性与确定性,最终都来源于这一深刻的数学规律。


练习

选择题

  • 题一(基本概率计算)

掷一枚均匀骰子,点数大于 444 的概率为

A. 16\dfrac{1}{6}61​   B. 13\dfrac{1}{3}31​   C. 12\dfrac{1}{2}21​   D. 23\dfrac{2}{3}32​

答案:B

点数大于 444 的结果为 {5, 6}\{5,\,6\}{5,6},共 222 个,样本空间共 666 个等可能结果,故:

P=26=13P = \frac{2}{6} = \frac{1}{3}P=62​=31​

  • 题二(期望值的计算)

某随机变量 XXX 的分布为:P(X=1)=0.2P(X=1)=0.2P(X=1)=0.2,P(X=2)=0.5P(X=2)=0.5P(X=2)=0.5,P(X=3)=0.3P(X=3)=0.3P(X=3)=0.3,期望值 ⟨X⟩\langle X \rangle⟨X⟩ 为

A. 1.81.81.8   B. 2.02.02.0   C. 2.12.12.1   D. 2.52.52.5

答案:C

⟨X⟩=1×0.2+2×0.5+3×0.3=0.2+1.0+0.9=2.1\langle X \rangle = 1 \times 0.2 + 2 \times 0.5 + 3 \times 0.3 = 0.2 + 1.0 + 0.9 = 2.1⟨X⟩=1×0.2+2×0.5+3×0.3=0.2+1.0+0.9=2.1

  • 题三(二项分布)

抛一枚均匀硬币 333 次,恰好出现 222 次正面的概率为

A. 18\dfrac{1}{8}81​   B. 28\dfrac{2}{8}82​   C. 38\dfrac{3}{8}83​   D. 48\dfrac{4}{8}84​

答案:C

由二项分布公式,N=3N=3N=3,n=2n=2n=2,p=1/2p=1/2p=1/2:

P ⁣(2; 3, 12)=(32)(12) ⁣2(12) ⁣1=3×14×12=38P\!\left(2;\,3,\,\tfrac{1}{2}\right) = \binom{3}{2}\left(\frac{1}{2}\right)^{\!2}\left(\frac{1}{2}\right)^{\!1} = 3 \times \frac{1}{4} \times \frac{1}{2} = \frac{3}{8}P(2;3,21​)=(23​)(21​)2(21​)1=3×41​×21​=83​

  • 题四(高斯分布性质)

某测量值服从正态分布,均值 μ=50\mu = 50μ=50,标准差 σ=2\sigma = 2σ=2。测量值落在区间 [46, 54][46,\,54][46,54] 内的概率约为

A. 68.3%68.3\%68.3%   B. 95.4%95.4\%95.4%   C. 99.7%99.7\%99.7%   D. 50.0%50.0\%50.0%

答案:B

区间 [46, 54]=[μ−2σ,  μ+2σ][46,\,54] = [\mu - 2\sigma,\; \mu + 2\sigma][46,54]=[μ−2σ,μ+2σ],根据高斯分布的“68–95–99.7 规则”,落在均值 ±2σ\pm 2\sigma±2σ 内的概率约为 95.4%95.4\%95.4%。


计算题

  • 题五(二项分布的期望与标准差)

一个容器被隔板分为体积相等的左右两半,其中共有 N=1000N = 1000N=1000 个气体分子,每个分子独立地以概率 p=0.5p = 0.5p=0.5 出现在左半。设左半的分子数为 nnn。求:① nnn 的期望值 ⟨n⟩\langle n \rangle⟨n⟩;② nnn 的标准差 σn\sigma_nσn​;③ nnn 相对于期望值的相对涨落 σn/⟨n⟩\sigma_n / \langle n \rangleσn​/⟨n⟩。

  • ① 期望值
⟨n⟩=Np=1000×0.5=500\langle n \rangle = Np = 1000 \times 0.5 = 500⟨n⟩=Np=1000×0.5=500
  • ② 标准差
σn=Np(1−p)=1000×0.5×0.5=250≈15.8\sigma_n = \sqrt{Np(1-p)} = \sqrt{1000 \times 0.5 \times 0.5} = \sqrt{250} \approx 15.8σn​=Np(1−p)​=1000×0.5×0.5​=250​≈15.8
  • ③ 相对涨落
σn⟨n⟩=15.8500≈3.16%=1N=11000≈3.16%\frac{\sigma_n}{\langle n \rangle} = \frac{15.8}{500} \approx 3.16\% = \frac{1}{\sqrt{N}} = \frac{1}{\sqrt{1000}} \approx 3.16\%⟨n⟩σn​​=50015.8​≈3.16%=N​1​=1000​1​≈3.16%

若粒子数增至 N=1023N = 10^{23}N=1023,相对涨落降至约 10−1210^{-12}10−12,宏观上完全不可察觉,气体表现为完美的均匀分布。


  • 题六(连续均匀分布的期望与方差)

一个随机变量 XXX 在区间 [0, 4][0,\,4][0,4] 上均匀分布,概率密度函数为 f(x)=1/4f(x) = 1/4f(x)=1/4(区间外为零)。求:① 期望值 ⟨X⟩\langle X \rangle⟨X⟩;② 方差 Var(X)\text{Var}(X)Var(X);③ 标准差 σ\sigmaσ。

  • ① 期望值
⟨X⟩=∫04x⋅14 dx=14⋅x22 ∣04=14×8=2\langle X \rangle = \int_0^4 x \cdot \frac{1}{4} \, \mathrm{d}x = \frac{1}{4} \cdot \frac{x^2}{2}\,\Bigg|_0^4 = \frac{1}{4} \times 8 = 2⟨X⟩=∫04​x⋅41​dx=41​⋅2x2​​04​=41​×8=2
  • ② 方差
⟨X2⟩=∫04x2⋅14 dx=14⋅x33 ∣04=14×643=163\langle X^2 \rangle = \int_0^4 x^2 \cdot \frac{1}{4} \, \mathrm{d}x = \frac{1}{4} \cdot \frac{x^3}{3}\,\Bigg|_0^4 = \frac{1}{4} \times \frac{64}{3} = \frac{16}{3}⟨X2⟩=∫04​x2⋅41​dx=41​⋅3x3​​04​=41​×364​=316​Var(X)=⟨X2⟩−⟨X⟩2=163−4=43≈1.33\text{Var}(X) = \langle X^2 \rangle - \langle X \rangle^2 = \frac{16}{3} - 4 = \frac{4}{3} \approx 1.33Var(X)=⟨X2⟩−⟨X⟩2=316​−4=34​≈1.33
  • ③ 标准差
σ=43=23≈1.15\sigma = \sqrt{\frac{4}{3}} = \frac{2}{\sqrt{3}} \approx 1.15σ=34​​=3​2​≈1.15

均匀分布 [0, L][0,\,L][0,L] 的通用结论:⟨X⟩=L/2\langle X \rangle = L/2⟨X⟩=L/2,Var(X)=L2/12\text{Var}(X) = L^2/12Var(X)=L2/12。本题 L=4L=4L=4,代入验证:L/2=2L/2 = 2L/2=2,L2/12=16/12=4/3L^2/12 = 16/12 = 4/3L2/12=16/12=4/3,与积分结果完全吻合。

  • 随机现象与概率
  • 离散概率分布
  • 期望值与方差
  • 连续概率分布
  • 独立随机变量的叠加
  • 二项分布
  • 高斯近似
  • 练习
    • 选择题
    • 计算题

目录

  • 随机现象与概率
  • 离散概率分布
  • 期望值与方差
  • 连续概率分布
  • 独立随机变量的叠加
  • 二项分布
  • 高斯近似
  • 练习
    • 选择题
    • 计算题