跳转至

伽马与贝塔函数

物理问题引入:当阶乘推广到连续与非整数维度

在初等排列组合中,阶乘 n!=n×(n−1)×⋯×1n! = n * (n-1) * dots * 1 只对非负整数有定义. 然而,在统计物理、高维空间与现代量子场论中,物理学家不可避免地遇到以下问题:

  1. nn 维空间超球体积:统计力学中由 NN 个质点组成的理想气体,其动量相空间是 3N3N 维欧氏空间!计算系统微态数需要求解 DD 维超球体的体积,而当 DD 为奇数时,公式中自然出现了「半整数阶乘」(D/2)!(D/2)!;
  2. 高斯积分高阶矩:高斯分布的偶数矩与奇数矩平均值计算;
  3. 量子场论维度正则化:在消除费曼(Feynman)图紫外发散时,物理学家将时空维度解析延拓为连续变量 D=4−ϵD = 4 - epsilon.

莱昂哈德·欧拉(Leonhard Euler) 发现的 伽马函数(Gamma Function),正是阶乘在全复平面上的唯美延拓.


1. 伽马函数 Γ(z)Gamma(z) 的定义与基本性质

1.1 欧拉第二类积分定义

对于实部为正的复数(Re(z)>0Re(z) > 0),伽马函数定义为收敛的反常积分:

Γ(z)=∫0∞tz−1e−tdtGamma(z) = int_0^inf t^z-1 e^-t dt

1.2 基本递推性质(分部积分)

利用分部积分法:

Γ(z+1)=∫0∞tze−tdt=−tze−t|0∞+z∫0∞tz−1e−tdt=zΓ(z)Gamma(z+1) = int_0^inf t^z e^-tdt = left. -t^z e^-t right|_0^inf + z int_0^inf t^z-1 e^-tdt = z Gamma(z)

由于 Γ(1)=∫0∞e−tdt=1Gamma(1) = int_0^inf e^-tdt = 1,对任意正整数 nn:

Γ(n+1)=n!⟺Γ(n)=(n−1)!Gamma(n+1) = n! quad iff quad Gamma(n) = (n-1)!

伽马函数完美复刻了自然数的阶乘规律!

1.3 半整数特殊值与高斯积分

令 z=1/2z = 1/2,做变量代换 t=u2t = u^2(dt=2ududt = 2udu):

Γ(1/2)=∫0∞t−1/2e−tdt=2∫0∞e−u2du=πGamma(1/2) = int_0^inf t^-1/2 e^-tdt = 2int_0^inf e^-u^2du = pi

借助递推公式 Γ(z+1)=zΓ(z)Gamma(z+1) = zGamma(z),所有半整数的伽马值均可直接算出:

Γ(3/2)=12π,Γ(5/2)=34π,Γ(n+1/2)=(2n−1)!!2nπGamma(3/2) = 12pi, quad Gamma(5/2) = 34pi, quad Gamma(n + 1/2) = (2n-1)!!2^npi

2. 解析延拓与欧拉余元公式

通过递推关系 Γ(z)=Γ(z+1)zGamma(z) = Gamma(z+1)z,伽马函数可以一步步解析延拓到全复平面,唯一的奇点是 非正整数点 z=0,−1,−2,…z = 0, -1, -2, dots(全为一阶单极点).

欧拉余元公式(Reflection Formula)

利用复围道积分可证明著名的反射公式:

Γ(z)Γ(1−z)=πsin⁡(πz)Gamma(z) Gamma(1 - z) = sin(pi z)

令 z=1/2z = 1/2 时立刻自洽验证了 Γ(1/2)2=πGamma(1/2)^2 = pi.


3. 贝塔函数 B(p,q)B(p, q)

3.1 欧拉第一类积分定义

对于 Re(p)>0,Re(q)>0Re(p) > 0, Re(q) > 0:

B(p,q)=∫01tp−1(1−t)q−1dtB(p, q) = int_0^1 t^p-1 (1 - t)^q-1 dt

通过变量代换 t=sin2⁡θt = sin^2theta,得到极具物理实用价值的三角函数形式:

B(p,q)=2∫0π/2sin2p−1⁡θcos2q−1⁡θdθB(p, q) = 2 int_0^pi/2 sin^2p-1theta cos^2q-1theta , dtheta

3.2 贝塔与伽马函数的统一转换公式

将两个伽马积分相乘并转换为二维极坐标计算,可直接导出著名的联立公式:

B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p, q) = Gamma(p)Gamma(q)Gamma(p + q)
计算三角积分的终极捷径

在电动力学和量子力学中计算诸如 ∫0π/2sin4⁡θcos2⁡θdθint_0^pi/2sin^4theta cos^2thetadtheta 时,无需反复调用降幂和差化积公式,只需认出:2p−1=4⟹p=5/22p-1 = 4 implies p = 5/2,2q−1=2⟹q=3/22q-1 = 2 implies q = 3/2. 该积分值直接为 12B(5/2,3/2)=Γ(5/2)Γ(3/2)2Γ(4)=π3212B(5/2, 3/2) = Gamma(5/2)Gamma(3/2)2Gamma(4) = 32!


4. 经典物理建模:DD 维相空间超球体积

在统计物理中,单原子理想气体系统的微观能量由动量给出:∑i=13Npi22m≤Esum_i=1^3N p_i^2/2m <= E. 这是一个半径为 R=2mER = 2mE 的 D=3ND = 3N 维超球体.

超球体积推导

DD 维球体体积可表示为:VD(R)=CDRDV_D(R) = C_D R^D.利用高斯多重积分技巧:

∫−∞∞…∫−∞∞e−(x12+⋯+xD2)dx1…dxD=(π)D=πD/2int_-inf^inf dots int_-inf^inf e^-(x_1^2 + dots + x_D^2) dx_1 dots dx_D = (pi)^D = pi^D/2

转换为高维球极坐标积分:∫0∞e−r2⋅(DCDrD−1)dr=DCD⋅12Γ(D/2)=CDΓ(D/2+1)int_0^inf e^-r^2 * (D C_D r^D-1) dr = D C_D * 12Gamma(D/2) = C_D Gamma(D/2 + 1). 对比两边,直接得出任意 DD 维超球体的精确体积公式:

VD(R)=πD/2Γ(D2+1)RDV_D(R) = pi^D/2Gammaleft(D2 + 1right) R^D

将 D=3D=3 代入:Γ(5/2)=34π⟹V3(R)=π3/234πR3=43πR3Gamma(5/2) = 3/4pi implies V_3(R) = pi^3/23/4pi R^3 = 43pi R^3; 将 D=2D=2 代入:Γ(2)=1⟹V2(R)=πR2Gamma(2) = 1 implies V_2(R) = pi R^2.完全精确吻合初等几何!


5. 学习衔接

  • 下一节:球对称静态边值与多极展开,进入 勒让德多项式;
  • 柱对称前置:在 贝塞尔函数 中,半整数阶贝塞尔函数将频繁使用伽马函数表达.


评论