1. 引言：1889年对“均匀随机”直觉的致命一击

在初等概率论与日常直觉中，“随机选取（at random）”似乎是一个不言自明的概念。掷一枚均匀骰子出现点数1的概率是 \(1/6\)，抛一枚对称硬币正面向上的概率是 \(1/2\)。当样本空间是有限离散集合时，计数测度赋予了每个基本事件同等的权重。

然而，一旦我们从有限离散集合迈入连续几何空间（Continuum），这种朴素的“等可能性”直觉便会瞬间崩塌。1889年，法国数学家约瑟夫·贝特兰（Joseph Bertrand, 1822–1900）在其名著《概率演算（Calcul des probabilités）》中提出了一个看似极其初等的几何概率问题——后世称之为“贝特兰悖论（Bertrand’s Paradox）”。

设在欧几里得平面 \(\mathbb{R}^2\) 上有一个半径 \(R = 1\) 的单位圆 \(S^1 = \{(x,y) : x^2 + y^2 = 1\}\)，其内接正三角形 \(\triangle ABC\) 的边长为 \(L = 2R\sin(\pi/3) = \sqrt{3} \approx 1.732\)，圆心 \(O\) 到正三角形各边的垂直距离（即内切圆半径）为 \(r_{\text{in}} = R\cos(\pi/3) = 1/2\)。贝特兰提出如下问题：

【贝特兰问题】 在单位圆内“随机”作一条弦 \(c\)。求该弦的长度大于内接正三角形边长 \(\sqrt{3}\) 的概率 \(P(c > \sqrt{3})\)。

令人震惊的是，贝特兰给出了三种完全自然、对称且无懈可击的几何构造方式，却分别严格推导出了三个彼此矛盾的精确答案：\(P_1 = 1/3\)、\(P_2 = 1/2\) 与 \(P_3 = 1/4\)。

  flowchart TD
    Q["贝特兰问题 (1889)<br/>在单位圆内随机选一条弦，<br/>弦长 c > √3 的概率是多少？"]
    
    Q --> M1["方法1：圆周上的随机端点<br/>圆心角 θ ~ Uniform(0, π)"]
    Q --> M2["方法2：固定半径上的随机距离<br/>弦心距 r ~ Uniform(0, 1)"]
    Q --> M3["方法3：圆盘内的随机中点<br/>中点坐标 (x, y) ~ Uniform(D²)"]
    
    M1 --> D1["弦心距概率密度函数<br/>f₁(r) = 2 / (π√(1 - r²))"]
    M2 --> D2["弦心距概率密度函数<br/>f₂(r) = 1"]
    M3 --> D3["弦心距概率密度函数<br/>f₃(r) = 2r"]
    
    D1 --> P1["概率 P₁(c > √3) = 1/3 ≈ 0.3333"]
    D2 --> P2["概率 P₂(c > √3) = 1/2 = 0.5000"]
    D3 --> P3["概率 P₃(c > √3) = 1/4 = 0.2500"]

2. 三种经典解法的严格推导与概率密度函数

设弦 \(c\) 所对的圆心角为 \(\theta \in [0, \pi]\)，圆心 \(O\) 到弦的垂直距离（弦心距）为 \(r \in [0, 1]\)，法线角为 \(\phi \in [0, 2\pi)\)，弦的中点为 \(M(x, y) = (r\cos\phi, r\sin\phi) \in D^2\)。由勾股定理，弦长 \(c\) 与弦心距 \(r\) 及圆心角 \(\theta\) 满足基本恒等式：

因此，目标事件 \(E = \{c > \sqrt{3}\}\) 在三种坐标系下完全等价于：

2.1 方法1：圆周上随机选取两个端点（\(P_1 = 1/3\)）

由旋转对称性 \(SO(2)\)，固定弦的一个端点于正三角形的顶点 \(A(\theta = 0)\)，在圆周上均匀选取另一个端点 \(B\)，使其极角 \(\theta \sim \text{Uniform}(0, 2\pi)\)。当且仅当点 \(B\) 落在正三角形另外两个顶点之间的弧段 \(\theta \in (2\pi/3, 4\pi/3)\) 上时，弦长 \(c > \sqrt{3}\)。该弧长为 \(2\pi/3\)，占整个圆周长 \(2\pi\) 的三分之一：

利用 \(r = \cos(\theta/2)\)（\(\theta \in [0, \pi]\)，\(f_\Theta(\theta) = 1/\pi\)）及变量代换雅可比行列式 \(\left|\frac{d\theta}{dr}\right| = \frac{2}{\sqrt{1 - r^2}}\)，可得弦心距 \(r\) 的概率密度函数为反正弦分布：

2.2 方法2：半径上随机选取弦心距（\(P_2 = 1/2\)）

先随机选取一条半径方向 \(\phi \sim \text{Uniform}(0, 2\pi)\)，再在该半径上均匀选取弦的中点到圆心的距离 \(r \sim \text{Uniform}(0, 1)\)，并作垂直于该半径的弦。由于内接正三角形的边到圆心的距离恰为 \(r_{\text{in}} = 1/2\)，弦长 \(c > \sqrt{3}\) 当且仅当 \(0 \le r < 1/2\)：

2.3 方法3：单位圆盘内按面积均匀选取弦的中点（\(P_3 = 1/4\)）

每条非直径的弦由其中点 \(M(x, y) \in D^2\) 唯一确定。若在单位圆盘 \(D^2\) 上按二维勒贝格面积测度 \(dx\,dy\) 均匀选取中点 \(M\)（即 \(f_{X,Y}(x,y) = 1/\pi\)），则 \(c > \sqrt{3}\) 当且仅当中点 \(M\) 落在半径为 \(r_{\text{in}} = 1/2\) 的同心内切圆盘 \(D_{1/2}\) 内部。由面积之比可得：

为什么方法2与方法3结果不同？因为在极坐标变换下，面积微元包含雅可比因子 \(J = r\)，即 \(dx\,dy = r\,dr\,d\phi\)。对角度 \(\phi \in [0, 2\pi)\) 积分后，径向距离 \(r\) 的边缘密度为线性函数：

3. 拉普拉斯“不充分理由原则”的危机与凯恩斯“酒水悖论”（1921）

古典概率论奠基人皮埃尔-西蒙·拉普拉斯（Pierre-Simon Laplace）提出了著名的“不充分理由原则（Principle of Insufficient Reason）”，后被经济学家与概率哲学家约翰·梅纳德·凯恩斯（John Maynard Keynes, 1921）改称为“无差异原则（Principle of Indifference）”：若我们对若干互斥且完备的可能结果没有任何偏袒的理由，则应赋予它们相等的概率。

然而在连续空间中，若对参数 \(x \in [a, b]\) 赋予均匀密度 \(f_X(x) = \frac{1}{b-a}\)，则在任何非线性双射重参数化 \(y = g(x)\) 下，诱导密度会乘上雅可比因子（Jacobian）：

因此，“在参数 \(x\) 下均匀”与“在参数 \(y = g(x)\) 下均匀”在数学上是互斥的！凯恩斯在1921年《概率论（A Treatise on Probability）》中用“酒水悖论（Wine-Water Paradox）”将这一荒谬性推向极致：

• 设杯中酒与水的混合物满足：酒与水的比例 \(x = W_{\text{wine}}/W_{\text{water}}\) 介于 \(1/3\) 与 \(3\) 之间。若假设 \(x \sim \text{Uniform}(1/3, 3)\)，则中位数为 \(x_{\text{med}} = \frac{1/3 + 3}{2} = \frac{5}{3}\)，且 \(P_X(x \le 2) = \frac{2 - 1/3}{3 - 1/3} = \frac{5}{8} = 0.625\)。
• 但同一前提完全等价于：水与酒的比例 \(y = 1/x = W_{\text{water}}/W_{\text{wine}}\) 介于 \(1/3\) 与 \(3\) 之间！若同样按无差异原则假设 \(y \sim \text{Uniform}(1/3, 3)\)，则中位数也是 \(y_{\text{med}} = \frac{5}{3}\)（这与 \(x \cdot y = 1\) 直接矛盾！），且同一事件 \(x \le 2 \iff y \ge 1/2\) 的概率变为 \(P_Y(y \ge 1/2) = \frac{3 - 1/2}{3 - 1/3} = \frac{15}{16} = 0.9375 \neq \frac{5}{8}\)！

4. 现代数学与物理学的三大消解方案

  flowchart TD
    Crisis["拉普拉斯无差异原则在连续空间中的崩溃<br/>非线性变换 y = g(x) 产生雅可比因子 |dx/dy|"]
    
    Crisis --> R1["消解方案1：柯尔莫哥洛夫 (1933)<br/>测度论公理体系 (Ω, F, P)"]
    Crisis --> R2["消解方案2：杰恩斯 (1973) 与积分几何<br/>欧几里得运动群 E(2) 不变测度"]
    Crisis --> R3["消解方案3：杰弗里斯 (1946) 与信息几何<br/>费希尔信息度量体积元"]
    
    R1 --> R1D["未显式指定概率测度 P 的问题是不适定的；<br/>三种物理抽样机制分别对应不同的概率空间。"]
    R2 --> R2D["满足旋转、缩放与平移不变性的唯一运动学测度为<br/>dμ = dp dφ（方法2，P = 1/2，柯西-克罗夫顿公式）。"]
    R3 --> R3D["采用杰弗里斯先验 p(θ) ∝ √det I(θ) 或乘法群哈尔测度 dx/x，<br/>严格实现重参数化不变性。"]

4.1 柯尔莫哥洛夫的测度论公理化 \((\Omega, \mathcal{F}, \mathbb{P})\)（1933）

1933年，苏联数学家安德烈·柯尔莫哥洛夫（Andrei Kolmogorov）在《概率论基础概念》中确立了现代概率论的公理体系：任何概率问题必须建立在明确给定的三元组 \((\Omega, \mathcal{F}, \mathbb{P})\) 之上，其中 \(\Omega\) 为样本空间，\(\mathcal{F}\) 为事件 \(\sigma\)-代数，\(\mathbb{P}\) 为可列可加的概率测度。在柯尔莫哥洛夫框架下，“随机选一条弦”而未指定测度 \(\mathbb{P}\) 在数学上是不适定问题（Ill-posed problem）；三种解法只是在同一个可测空间 \((\mathcal{C}, \mathcal{B}(\mathcal{C}))\) 上计算了三个不同测度 \(\mathbb{P}_1, \mathbb{P}_2, \mathbb{P}_3\) 的值。

4.2 杰恩斯变换群不变性（1973）与庞加莱-桑塔洛运动学测度

1973年，物理学家埃德温·杰恩斯（Edwin T. Jaynes）在经典论文 “The Well-Posed Problem” 中指出：如果我们站在远处向画在地面上的小圆随机投掷长麦秆，由于圆的朝向、大小和位置相对于投掷者是任意的，真实的物理分布必须同时满足：

1. 旋转不变性 \(SO(2)\)：\(f(p, \phi) = g(p)/(2\pi)\)（三种方法均满足）；
2. 尺度缩放不变性：限制在半径 \(a < 1\) 的同心小圆内时分布函数形式不变（排除具有边界奇异性的方法1）；
3. 平面平移不变性 \(\mathbb{R}^2\)（欧几里得刚体运动群 \(E(2) = \mathbb{R}^2 \rtimes O(2)\)）：将原点平移 \((x_0, y_0)\) 后，直线 \(x\cos\phi + y\sin\phi = p\) 的新参数为 \(p' = p - x_0\cos\phi - y_0\sin\phi\)，\(\phi' = \phi\)。其变换雅可比行列式恒等于 \(1\)： \(J = \det\frac{\partial(p', \phi')}{\partial(p, \phi)} = \det\begin{pmatrix} 1 & x_0\sin\phi - y_0\cos\phi \\ 0 & 1 \end{pmatrix} = 1\) 要求 \(g(p') = g(p)\) 对任意平移 \((x_0, y_0)\) 成立，唯一迫使 \(g(p) = \text{const.}\)！

这正是积分几何（Integral Geometry）中由庞加莱（Poincaré）与桑塔洛（Santaló）确立的平面直线空间上的唯一 \(E(2)\)-不变运动学测度（Kinematic Measure / 哈尔测度）：

由积分几何的柯西-克罗夫顿公式（Cauchy–Crofton Formula），穿过平面上任意凸体 \(K\) 的直线集合的运动学测度等于该凸体的周长 \(\text{Perimeter}(\partial K)\)（这也是蒲丰投针问题 Buffon’s Needle 的理论基础）。因此，穿过内切圆 \(K_{1/2}\)（周长 \(\pi\)）与穿过单位圆 \(K_1\)（周长 \(2\pi\)）的测度之比立即给出方法2的唯一物理解：

4.3 杰弗里斯无信息先验（Jeffreys Prior）与重参数化不变性

在贝叶斯统计中，为了消除坐标变换雅可比因子带来的主观性，哈罗德·杰弗里斯（Harold Jeffreys, 1946）利用统计流形上的费希尔信息矩阵（Fisher Information Matrix）\(I_{ij}(\theta) = -\mathbb{E}\left[\frac{\partial^2 \ln p(X\mid\theta)}{\partial\theta_i\partial\theta_j}\right]\) 提出了杰弗里斯先验：

由于在参数变换 \(\eta = \psi(\theta)\) 下，\(\sqrt{\det I(\eta)} = \left|\det\frac{\partial\theta}{\partial\eta}\right|\sqrt{\det I(\theta)}\)，体积元 \(p_J(\theta)\,d\theta = p_J(\eta)\,d\eta\) 严格保持不变。对于凯恩斯酒水悖论中的正比例参数 \(x \in [1/3, 3]\)，乘法群 \((\mathbb{R}_{>0}, \times)\) 的不变哈尔测度为 \(d\mu(x) = \frac{dx}{x}\)，归一化密度 \(f(x) = \frac{1}{2\ln 3 \cdot x}\) 在倒数变换 \(y = 1/x\) 下完全不变，中位数恒为 \(x_{\text{med}} = y_{\text{med}} = 1\)，彻底消除了悖论。

5. Python 蒙特卡洛模拟：三种经典机制与杰恩斯偏心圆投针实验

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

import numpy as np

def simulate_bertrand(n: int = 1_000_000, seed: int = 42) -> None:
    rng = np.random.default_rng(seed)
    limit = np.sqrt(3.0)

    # 方法1：圆周上的随机端点 (P = 1/3)
    t1 = rng.uniform(0, 2 * np.pi, n)
    t2 = rng.uniform(0, 2 * np.pi, n)
    c1 = 2.0 * np.sin(np.abs(t1 - t2) / 2.0)

    # 方法2：半径上的随机距离 r ~ U(0, 1) (P = 1/2)
    r2 = rng.uniform(0.0, 1.0, n)
    c2 = 2.0 * np.sqrt(1.0 - r2**2)

    # 方法3：圆盘内的随机中点 u = r^2 ~ U(0, 1) (P = 1/4)
    r3 = np.sqrt(rng.uniform(0.0, 1.0, n))
    c3 = 2.0 * np.sqrt(1.0 - r3**2)

    # 方法4：杰恩斯物理实验（向大区域随机投掷直线，观测平移至 (3.5, -2.0) 的单位圆）
    phi = rng.uniform(0, 2 * np.pi, n * 5)
    p = rng.uniform(-10.0, 10.0, n * 5)
    r_shifted = np.abs(p - (3.5 * np.cos(phi) - 2.0 * np.sin(phi)))
    hits = r_shifted[r_shifted <= 1.0]
    c_jaynes = 2.0 * np.sqrt(1.0 - hits**2)

    for name, arr, exact in [
        ("Method 1 (Endpoints)", c1, 1/3),
        ("Method 2 (Radius)", c2, 1/2),
        ("Method 3 (Midpoint)", c3, 1/4),
        ("Jaynes Off-Center Toss", c_jaynes, 1/2),
    ]:
        print(f"{name:<24}: P(c > √3) = {np.mean(arr > limit):.5f} (Exact: {exact:.5f})")

if __name__ == "__main__":
    simulate_bertrand()

模拟结果精确收敛于 \(0.33341\)（\(1/3\)）、\(0.50012\)（\(1/2\)）、\(0.24989\)（\(1/4\)），且偏心圆物理投掷实验收敛于 \(0.49994\)（\(1/2\)），完美印证了欧几里得运动群 \(E(2)\) 不变测度的唯一性。