前言

Beta 函数 ( B ( x , y ) B(x, y) B(x,y)) 和 Gamma 函数 ( Γ ( z ) \Gamma(z) Γ(z)) 是特殊函数领域中最为基础且联系紧密的两个函数。它们在物理、工程以及量子信息科学等领域有广泛应用。

一、Gamma函数的定义

Gamma 函数阶乘函数的推广,也是连接离散与连续数学的桥梁。普通阶乘 n ! n! n! 只对整数有定义,但Gamma函数可以定义在几乎所有复数上,这让我们能:

  1. 把离散的阶乘运算扩展到连续变量(比如积分里的非整数次幂)
  2. 处理组合数、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)=0tz1etdt,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=uvababvdu

令 $ 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=ztz1dt
v = ∫ e − t d t = − e − t v = \int e^{-t} dt = -e^{-t} v=etdt=et

那么
Γ ( 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)=0tzetdt=uv00vdu=tz(et) 00(et)ztz1dt=tzet 0+z0tz1etdt=0+z0tz1etdt=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)=0etdt=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)=21=2!

  • Γ ( 4 ) = 3 ⋅ Γ ( 3 ) = 3 ⋅ 2 ⋅ 1 = 3 ! \Gamma(4) = 3 \cdot \Gamma(3) = 3 \cdot 2 \cdot 1 = 3! Γ(4)=3Γ(3)=321=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)=01xa1(1x)b1dx,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)=01xa1(1x)b1dt=Γ(a+b)Γ(a)Γ(b)=(a+b1)!(a1)!(b1)!
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(1x)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+xp1exdx)(0+yq1eydy)=0+0+xp1yq1e(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=(1u)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= uxuyvxvy = vvu1u =v(1u)u(v)=v

因此面积微元为:

d x d y = ∣ J ∣   d v d u = v   d v d u dxdy = |J| \, dvdu = v \, dvdu dxdy=Jdvdu=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)=010+(uv)p1[(1u)v]q1evvdvdu=01up1(1u)q1du0+vp+q1evdv

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) 01up1(1u)q1du=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+q1evdv=Γ(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!

Logo

openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构

更多推荐