伽马与贝塔函数
物理问题引入:当阶乘推广到连续与非整数维度
在初等排列组合中,阶乘 𝑛! =𝑛 ×(𝑛 −1) ×⋯ ×1n! = n * (n-1) * dots * 1 只对非负整数有定义. 然而,在统计物理、高维空间与现代量子场论中,物理学家不可避免地遇到以下问题:
- 𝑛n 维空间超球体积:统计力学中由 𝑁N 个质点组成的理想气体,其动量相空间是 3𝑁3N 维欧氏空间!计算系统微态数需要求解 𝐷D 维超球体的体积,而当 𝐷D 为奇数时,公式中自然出现了「半整数阶乘」(𝐷/2)!(D/2)!;
- 高斯积分高阶矩:高斯分布的偶数矩与奇数矩平均值计算;
- 量子场论维度正则化:在消除费曼(Feynman)图紫外发散时,物理学家将时空维度解析延拓为连续变量 𝐷 =4 −𝜖D = 4 - epsilon.
莱昂哈德·欧拉(Leonhard Euler) 发现的 伽马函数(Gamma Function),正是阶乘在全复平面上的唯美延拓.
1. 伽马函数 Γ(𝑧)Gamma(z) 的定义与基本性质
1.1 欧拉第二类积分定义
对于实部为正的复数(Re(𝑧) >0Re(z) > 0),伽马函数定义为收敛的反常积分:
Γ(𝑧)=∫∞0𝑡𝑧−1𝑒−𝑡d𝑡Gamma(z) = int_0^inf t^z-1 e^-t dt1.2 基本递推性质(分部积分)
利用分部积分法:
Γ(𝑧+1)=∫∞0𝑡𝑧𝑒−𝑡d𝑡=−𝑡𝑧𝑒−𝑡∣∞0+𝑧∫∞0𝑡𝑧−1𝑒−𝑡d𝑡=𝑧Γ(𝑧)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𝑒−𝑡d𝑡 =1Gamma(1) = int_0^inf e^-tdt = 1,对任意正整数 𝑛n:
Γ(𝑛+1)=𝑛!⟺Γ(𝑛)=(𝑛−1)!Gamma(n+1) = n! quad iff quad Gamma(n) = (n-1)!伽马函数完美复刻了自然数的阶乘规律!
1.3 半整数特殊值与高斯积分
令 𝑧 =1/2z = 1/2,做变量代换 𝑡 =𝑢2t = u^2(d𝑡 =2𝑢d𝑢dt = 2udu):
Γ(1/2)=∫∞0𝑡−1/2𝑒−𝑡d𝑡=2∫∞0𝑒−𝑢2d𝑢=√𝜋Gamma(1/2) = int_0^inf t^-1/2 e^-tdt = 2int_0^inf e^-u^2du = pi借助递推公式 Γ(𝑧 +1) =𝑧Γ(𝑧)Gamma(z+1) = zGamma(z),所有半整数的伽马值均可直接算出:
Γ(3/2)=12√𝜋,Γ(5/2)=34√𝜋,Γ(𝑛+1/2)=(2𝑛−1)!!2𝑛√𝜋Gamma(3/2) = 12pi, quad Gamma(5/2) = 34pi, quad Gamma(n + 1/2) = (2n-1)!!2^npi
2. 解析延拓与欧拉余元公式
通过递推关系 Γ(𝑧) =Γ(𝑧+1)𝑧Gamma(z) = Gamma(z+1)z,伽马函数可以一步步解析延拓到全复平面,唯一的奇点是 非正整数点 𝑧 =0, −1, −2,…z = 0, -1, -2, dots(全为一阶单极点).
利用复围道积分可证明著名的反射公式:
Γ(𝑧)Γ(1−𝑧)=𝜋sin(𝜋𝑧)Gamma(z) Gamma(1 - z) = sin(pi z)令 𝑧 =1/2z = 1/2 时立刻自洽验证了 Γ(1/2)2 =𝜋Gamma(1/2)^2 = pi.
3. 贝塔函数 𝐵(𝑝,𝑞)B(p, q)
3.1 欧拉第一类积分定义
对于 Re(𝑝) >0,Re(𝑞) >0Re(p) > 0, Re(q) > 0:
𝐵(𝑝,𝑞)=∫10𝑡𝑝−1(1−𝑡)𝑞−1d𝑡B(p, q) = int_0^1 t^p-1 (1 - t)^q-1 dt通过变量代换 𝑡 =sin2𝜃t = sin^2theta,得到极具物理实用价值的三角函数形式:
𝐵(𝑝,𝑞)=2∫𝜋/20sin2𝑝−1𝜃cos2𝑞−1𝜃d𝜃B(p, q) = 2 int_0^pi/2 sin^2p-1theta cos^2q-1theta , dtheta3.2 贝塔与伽马函数的统一转换公式
将两个伽马积分相乘并转换为二维极坐标计算,可直接导出著名的联立公式:
𝐵(𝑝,𝑞)=Γ(𝑝)Γ(𝑞)Γ(𝑝+𝑞)B(p, q) = Gamma(p)Gamma(q)Gamma(p + q)计算三角积分的终极捷径
在电动力学和量子力学中计算诸如 ∫𝜋/20sin4𝜃cos2𝜃d𝜃int_0^pi/2sin^4theta cos^2thetadtheta 时,无需反复调用降幂和差化积公式,只需认出:2𝑝 −1 =4 ⟹ 𝑝 =5/22p-1 = 4 implies p = 5/2,2𝑞 −1 =2 ⟹ 𝑞 =3/22q-1 = 2 implies q = 3/2. 该积分值直接为 12𝐵(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. 经典物理建模:𝐷D 维相空间超球体积
在统计物理中,单原子理想气体系统的微观能量由动量给出:∑3𝑁𝑖=1𝑝2𝑖2𝑚 ≤𝐸sum_i=1^3N p_i^2/2m <= E. 这是一个半径为 𝑅 =√2𝑚𝐸R = 2mE 的 𝐷 =3𝑁D = 3N 维超球体.
超球体积推导
𝐷D 维球体体积可表示为:𝑉𝐷(𝑅) =𝐶𝐷𝑅𝐷V_D(R) = C_D R^D.利用高斯多重积分技巧:
∫∞−∞…∫∞−∞𝑒−(𝑥21+⋯+𝑥2𝐷)d𝑥1…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𝑒−𝑟2 ⋅(𝐷𝐶𝐷𝑟𝐷−1)d𝑟 =𝐷𝐶𝐷 ⋅12Γ(𝐷/2) =𝐶𝐷Γ(𝐷/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). 对比两边,直接得出任意 𝐷D 维超球体的精确体积公式:
𝑉𝐷(𝑅)=𝜋𝐷/2Γ(𝐷2+1)𝑅𝐷V_D(R) = pi^D/2Gammaleft(D2 + 1right) R^D将 𝐷 =3D=3 代入:Γ(5/2) =34√𝜋 ⟹ 𝑉3(𝑅) =𝜋3/234√𝜋𝑅3 =43𝜋𝑅3Gamma(5/2) = 3/4pi implies V_3(R) = pi^3/23/4pi R^3 = 43pi R^3; 将 𝐷 =2D=2 代入:Γ(2) =1 ⟹ 𝑉2(𝑅) =𝜋𝑅2Gamma(2) = 1 implies V_2(R) = pi R^2.完全精确吻合初等几何!
5. 学习衔接
- 下一节:球对称静态边值与多极展开,进入 勒让德多项式;
- 柱对称前置:在 贝塞尔函数 中,半整数阶贝塞尔函数将频繁使用伽马函数表达.
本页面最近更新:2026/10/1 10:49:01,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用