Gamma函数与Beta函数及其关系
本文介绍了Beta函数和Gamma函数的定义及其相互关系。Gamma函数作为阶乘函数的推广,通过欧拉积分形式定义在复数域上,满足递推关系Γ(z+1)=zΓ(z),并与阶乘有直接对应关系Γ(n+1)=n!。Beta函数定义为第一类欧拉积分,与Gamma函数存在重要关系式B(a,b)=Γ(a)Γ(b)/Γ(a+b)。通过变量代换和雅克比行列式证明了这一关系,并推导出当参数为整数时的特殊形式B(n+1,
文章目录
前言
Beta 函数 ( B ( x , y ) B(x, y) B(x,y)) 和 Gamma 函数 ( Γ ( z ) \Gamma(z) Γ(z)) 是特殊函数领域中最为基础且联系紧密的两个函数。它们在物理、工程以及量子信息科学等领域有广泛应用。
一、Gamma函数的定义
Gamma 函数阶乘函数的推广,也是连接离散与连续数学的桥梁。普通阶乘 n ! n! n! 只对整数有定义,但Gamma函数可以定义在几乎所有复数上,这让我们能:
- 把离散的阶乘运算扩展到连续变量(比如积分里的非整数次幂)
- 处理组合数、Beta函数、泊松分布等大量数学工具
Gamma函数 Γ ( z ) \Gamma(z) Γ(z) 的标准定义是欧拉积分形式:
Γ ( z ) = ∫ 0 ∞ t z − 1 e − t d t , Re ( z ) > 0 \Gamma(z) = \int_{0}^{\infty} t^{z-1} e^{-t} dt, \quad \text{Re}(z) > 0 Γ(z)=∫0∞tz−1e−tdt,Re(z)>0
其中 z z z 是复数变量,只要它的实部大于0,这个积分就收敛。
二、Gamma函数与阶乘的关系
对于非负整数 n n n,Gamma函数和阶乘直接对应:
Γ ( n + 1 ) = n ! \Gamma(n+1) = n! Γ(n+1)=n!
证明:
我们先来回顾一下分部积分公式,详情查看【分部积分法的详细解释、推导过程及应用示例】:
∫ a b u d v = u v ∣ a b − ∫ a b v d u \int_{a}^{b} u \, dv = \left. u v \right|_{a}^{b} - \int_{a}^{b} v \, du ∫abudv=uv∣ab−∫abvdu
令 $ u = t^{z} , , , dv = e^{-t}dt $ 则:
d u = d d t t z d t = z t z − 1 d t du = \frac{d}{dt} t^{z} dt = z t^{z-1} dt du=dtdtzdt=ztz−1dt
v = ∫ e − t d t = − e − t v = \int e^{-t} dt = -e^{-t} v=∫e−tdt=−e−t
那么
Γ ( z + 1 ) = ∫ 0 ∞ t z e − t d t = u v ∣ 0 ∞ − ∫ 0 ∞ v d u = t z ⋅ ( − e − t ) ∣ 0 ∞ − ∫ 0 ∞ ( − e − t ) ⋅ z t z − 1 d t = − t z e − t ∣ 0 ∞ + z ∫ 0 ∞ t z − 1 e − t d t = 0 + z ∫ 0 ∞ t z − 1 e − t d t = z Γ ( z ) \begin{aligned} \Gamma(z+1) &= \int_{0}^{\infty} t^{z} e^{-t} dt \\ &=\left. u v \right|_{0}^{\infty} - \int_{0}^{\infty} v \, du\\ &=\left. t^{z} \cdot (-e^{-t}) \right|_{0}^{\infty} - \int_{0}^{\infty} (-e^{-t}) \cdot z t^{z-1} dt\\ &= -\left. t^{z} e^{-t} \right|_{0}^{\infty}+z \int_{0}^{\infty} t^{z-1} e^{-t} dt\\ &=0+z \int_{0}^{\infty} t^{z-1} e^{-t} dt\\ &=z\Gamma(z) \end{aligned} Γ(z+1)=∫0∞tze−tdt=uv∣0∞−∫0∞vdu=tz⋅(−e−t)
0∞−∫0∞(−e−t)⋅ztz−1dt=−tze−t
0∞+z∫0∞tz−1e−tdt=0+z∫0∞tz−1e−tdt=zΓ(z)
t z t^z tz 是幂函数,增长较慢; e t e^t et 是指数函数,爆炸增长。 e t e^t et影响较大所以边界项为 0。
我们得到:
Γ ( z + 1 ) = z Γ ( z ) \boxed{\Gamma(z+1) = z \, \Gamma(z)} Γ(z+1)=zΓ(z)
接下来验证整数阶乘的特例
对于非负整数 n,我们用递推关系可以验证阶乘:
-
基础值: Γ ( 1 ) = ∫ 0 ∞ e − t d t = 1 = 0 ! \Gamma(1) = \int_{0}^{\infty} e^{-t} dt = 1 = 0! Γ(1)=∫0∞e−tdt=1=0!
-
Γ ( 2 ) = 1 ⋅ Γ ( 1 ) = 1 = 1 ! \Gamma(2) = 1 \cdot \Gamma(1) = 1 = 1! Γ(2)=1⋅Γ(1)=1=1!
-
Γ ( 3 ) = 2 ⋅ Γ ( 2 ) = 2 ⋅ 1 = 2 ! \Gamma(3) = 2 \cdot \Gamma(2) = 2 \cdot 1 = 2! Γ(3)=2⋅Γ(2)=2⋅1=2!
-
Γ ( 4 ) = 3 ⋅ Γ ( 3 ) = 3 ⋅ 2 ⋅ 1 = 3 ! \Gamma(4) = 3 \cdot \Gamma(3) = 3 \cdot 2 \cdot 1 = 3! Γ(4)=3⋅Γ(3)=3⋅2⋅1=3!
以此类推,可得: Γ ( n + 1 ) = n ! \Gamma(n+1) = n! Γ(n+1)=n!
三、Beta函数的定义
Beta 函数定义为:
B ( a , b ) = ∫ 0 1 x a − 1 ( 1 − x ) b − 1 d x , R e ( a ) > 0 , R e ( b ) > 0 \mathrm{B}(a,b)=\int_{0}^{1}x^{a-1}(1-x)^{b-1}\,\mathrm{d}x,\quad \mathrm{Re}(a)>0,\ \mathrm{Re}(b)>0 B(a,b)=∫01xa−1(1−x)b−1dx,Re(a)>0, Re(b)>0
- 定义域:(a,b>0)(实数)或复平面右半部分。
- 别称:第一类欧拉积分;与Gamma函数(第二类欧拉积分)紧密关联。
四、Beta函数与Gamma函数的关系
B ( a , b ) = ∫ 0 1 x a − 1 ( 1 − x ) b − 1 d t = Γ ( a ) Γ ( b ) Γ ( a + b ) = ( a − 1 ) ! ( b − 1 ) ! ( a + b − 1 ) ! B(a,b) = \int_{0}^{1} x^{a-1} (1-x)^{b-1} dt = \frac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)}=\frac{(a-1)!(b-1)!}{(a+b-1)!} B(a,b)=∫01xa−1(1−x)b−1dt=Γ(a+b)Γ(a)Γ(b)=(a+b−1)!(a−1)!(b−1)!
当 a = n + 1 , b = m + 1 a = n+1, b = m+1 a=n+1,b=m+1 时:
B ( n + 1 , m + 1 ) = ∫ 0 1 x m ( 1 − x ) n d t = Γ ( n + 1 ) Γ ( m + 1 ) Γ ( n + m + 2 ) = n ! m ! ( n + m + 1 ) ! B(n+1,m+1) = \int_{0}^{1} x^{m} (1-x)^{n} dt = \frac{\Gamma(n+1)\Gamma(m+1)}{\Gamma(n+m+2)} = \frac{n!m!}{(n+m+1)!} B(n+1,m+1)=∫01xm(1−x)ndt=Γ(n+m+2)Γ(n+1)Γ(m+1)=(n+m+1)!n!m!
证明
1. 从 Gamma 函数乘积出发
Γ ( p ) Γ ( q ) = ( ∫ 0 + ∞ x p − 1 e − x d x ) ( ∫ 0 + ∞ y q − 1 e − y d y ) = ∫ 0 + ∞ ∫ 0 + ∞ x p − 1 y q − 1 e − ( x + y ) d x d y \Gamma(p)\Gamma(q) = \left( \int_{0}^{+\infty} x^{p-1} e^{-x} dx \right) \left( \int_{0}^{+\infty} y^{q-1} e^{-y} dy \right) = \int_{0}^{+\infty} \int_{0}^{+\infty} x^{p-1} y^{q-1} e^{-(x+y)} \, dxdy Γ(p)Γ(q)=(∫0+∞xp−1e−xdx)(∫0+∞yq−1e−ydy)=∫0+∞∫0+∞xp−1yq−1e−(x+y)dxdy
这是一个二重积分,我们需要引入雅克比行列式。
2. 变量代换
令
{ x = u v y = ( 1 − u ) v ( 0 < u < 1 , v > 0 ) \begin{cases} x = uv \\ y = (1-u)v \end{cases} \quad (0 < u < 1,\ v > 0) {x=uvy=(1−u)v(0<u<1, v>0)
计算雅克比行列式:
∣ J ∣ = ∣ ∂ x ∂ u ∂ x ∂ v ∂ y ∂ u ∂ y ∂ v ∣ = ∣ v u − v 1 − u ∣ = v ( 1 − u ) − u ( − v ) = v |J| = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \\[1.5em] \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix}= \begin{vmatrix} v & u \\ -v & 1-u \end{vmatrix}= v(1-u) - u(-v) = v ∣J∣= ∂u∂x∂u∂y∂v∂x∂v∂y = v−vu1−u =v(1−u)−u(−v)=v
因此面积微元为:
d x d y = ∣ J ∣ d v d u = v d v d u dxdy = |J| \, dvdu = v \, dvdu dxdy=∣J∣dvdu=vdvdu
3. 代入被积函数并化简
Γ ( p ) Γ ( q ) = ∫ 0 1 ∫ 0 + ∞ ( u v ) p − 1 [ ( 1 − u ) v ] q − 1 e − v ⋅ v d v d u = ∫ 0 1 u p − 1 ( 1 − u ) q − 1 d u ⋅ ∫ 0 + ∞ v p + q − 1 e − v d v \begin{aligned} \Gamma(p)\Gamma(q) &= \int_{0}^{1} \int_{0}^{+\infty} (uv)^{p-1} \left[(1-u)v\right]^{q-1} e^{-v} \cdot v \, dvdu \\ &= \int_{0}^{1} u^{p-1} (1-u)^{q-1} du \cdot \int_{0}^{+\infty} v^{p+q-1} e^{-v} dv \end{aligned} Γ(p)Γ(q)=∫01∫0+∞(uv)p−1[(1−u)v]q−1e−v⋅vdvdu=∫01up−1(1−u)q−1du⋅∫0+∞vp+q−1e−vdv
4. 匹配 Beta 与 Gamma 定义
∫ 0 1 u p − 1 ( 1 − u ) q − 1 d u = B ( p , q ) \int_{0}^{1} u^{p-1} (1-u)^{q-1} du = \mathrm{B}(p,q) ∫01up−1(1−u)q−1du=B(p,q)
∫ 0 + ∞ v p + q − 1 e − v d v = Γ ( p + q ) \int_{0}^{+\infty} v^{p+q-1} e^{-v} dv = \Gamma(p+q) ∫0+∞vp+q−1e−vdv=Γ(p+q)
因此:
Γ ( p ) Γ ( q ) = B ( p , q ) Γ ( p + q ) \Gamma(p)\Gamma(q) = \mathrm{B}(p,q) \, \Gamma(p+q) Γ(p)Γ(q)=B(p,q)Γ(p+q)
移项即得核心结论:
B ( p , q ) = Γ ( p ) Γ ( q ) Γ ( p + q ) B ( n + 1 , m + 1 ) = Γ ( p ) Γ ( q ) Γ ( p + q ) = n ! m ! ( n + m + 1 ) ! \mathrm{B}(p,q) = \dfrac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}\\[1.5em] \mathrm{B}(n+1,m+1) = \dfrac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}=\frac{n!m!}{(n+m+1)!} B(p,q)=Γ(p+q)Γ(p)Γ(q)B(n+1,m+1)=Γ(p+q)Γ(p)Γ(q)=(n+m+1)!n!m!
openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构
更多推荐



所有评论(0)