超球体积可以由高斯积分统一求出:同一个积分,一方面按坐标分离,另一方面按半径分层,再借助伽马函数确定未知的体积系数。这一方法把高维几何问题化为一维积分,也为机器学习中的高维空间直觉提供基础。
体积公式与对象
设 n∈N+、r≥0,欧氏空间 Rn 中半径为 r 的闭球为
Bn(r)={x∈Rn:∥x∥2≤r}.
记其 n 维体积为 Vn(r)。这里求的是包含内部的球体体积;它的边界是 (n−1) 维球面。
Vn(r)=Γ(2n+1)πn/2rn.
令 Cn=Vn(1)。对 r>0,缩放 x=ru 的雅可比行列式绝对值为 rn,所以
Vn(r)=Cnrn.
r=0 时球体退化为一点,因 n≥1,体积为零,同一公式仍成立。以下只需确定 Cn。
高斯积分的坐标分离
一维高斯积分给出
∫−∞+∞e−x2dx=π.
考虑
In=∫Rne−∥x∥22dx1⋯dxn.
被积函数非负,且可以写成各坐标函数的乘积。由 Tonelli 定理可按坐标逐次积分:
In=∫Rni=1∏ne−xi2dx1⋯dxn=i=1∏n∫−∞+∞e−xi2dxi=(π)n=πn/2.
这里用的是乘积结构与非负积分定理,不需要额外引入随机变量的独立性假设。
按半径分层积分
被积函数只依赖半径
ρ=∥x∥2=x12+⋯+xn2.
对 ρ>0,厚度为 Δρ>0 的球壳体积为
Vn(ρ+Δρ)−Vn(ρ)=Cn((ρ+Δρ)n−ρn)=nCnρn−1Δρ+o(Δρ).
因此,径向分层积分的体积权重是
dVn(ρ)=nCnρn−1dρ.
这可由球壳上的阶梯函数积分,再取极限得到;三维情形与球面坐标下先积角变量的结果相同。于是
In=nCn∫0∞e−ρ2ρn−1dρ.
整个步骤只用了体积的缩放关系,尚未假设 Cn 的具体值。
化为伽马函数积分
对实数 s>0,伽马函数定义为
Γ(s)=∫0∞ts−1e−tdt,Γ(s+1)=sΓ(s).
令 t=ρ2,则
dρ=21t−1/2dt,ρn−1=t(n−1)/2.
因为 n/2>0,下述 Gamma 积分收敛:
∫0∞e−ρ2ρn−1dρ=21∫0∞e−ttn/2−1dt=21Γ(2n).
代回径向积分,得到
In=2nCnΓ(2n)=CnΓ(2n+1).
与坐标分离所得的 In=πn/2 比较:
Cn=Γ(2n+1)πn/2,Vn(r)=Γ(2n+1)πn/2rn.
Gamma 函数出现的直接原因,是径向权重 ρn−1 在 t=ρ2 换元后变成了 tn/2−1。
偶数维与低维核验
由 Γ(k+1)=k!,当 n=2k、k≥1 时,
V2k(r)=k!πkr2k.
奇数维对应 Gamma 函数的半整数值。由
Γ(21)=2∫0∞e−u2du=π
和递推关系,例如
Γ(25)=23⋅21π=43π,
便得到
V3(r)=Γ(5/2)π3/2r3=34πr3.
| 维度 n | Vn(r) | 几何对象 |
|---|
| 1 | 2r | 区间长度 |
| 2 | πr2 | 圆盘面积 |
| 3 | 34πr3 | 球体体积 |
| 4 | 2π2r4 | 四维球体积 |
维度递推与高维极限
对 n≥3、r>0,利用 Gamma 函数递推,
Vn−2(r)Vn(r)=Γ(n/2+1)πn/2rnπ(n−2)/2rn−2Γ(n/2)=n2πr2.
所以
Vn(r)=n2πr2Vn−2(r),n≥3.
该等式在 r=0 时也成立;计算正整数维度时,以 V1(r)=2r、V2(r)=πr2 为起点即可。
对单位球,
Cn=n2πCn−2.
当 n>2π 时,Cn<Cn−2,比较的是同奇偶性的维度。初期体积并非随维度一直减小,例如 C2>C1。
进一步,当 n≥13 时,递推系数 2π/n<1/2。因此奇数维、偶数维两个子序列都最终被等比递减序列控制,从而
n→∞limVn(1)=0.
与机器学习中的高维几何的联系
体积的缩放关系直接给出球内均匀采样的半径分布。设 r>0,X 在 Bn(r) 内按体积均匀分布,则对 0≤a≤1,
P(∥X∥2≤ar)=Vn(r)Vn(ar)=an.
因此,对固定的相对厚度 0<ε<1,外层球壳包含的概率为
P((1−ε)r<∥X∥2≤r)=1−(1−ε)n⟶1.
例如在 100 维球内均匀采样,最外侧 5% 半径厚度的球壳约包含
1−0.95100≈99.41%
的样本概率。这说明高维中的体积主要分布在靠近边界的区域,有助于理解高维采样与邻域覆盖为何不同于二维、三维直觉。这里的结论以球内均匀分布为前提,应用到数据时还需明确数据分布与特征尺度。
高斯积分的按坐标分解与按半径分层,也说明高斯分布与球面几何为何经常共同出现:各向同性高斯密度只依赖半径,可以利用同一径向体积权重分析。其半径分布还包含高斯密度的权重,不能直接套用球内均匀分布的 an。