ARTICLE

digamma函数

digamma函数 digamma函数(digamma function),又称Ψ函数(psi function),是伽马函数(gamma function)的对数导数,定义为: 它是数学分析、数论和概率论中一类至关重要的特殊函数。由于伽马函数本身是阶乘向复平面的解析延拓,digamma函数自然成为研究级数求和、积分表示和差分方程的核心工具。它也是更广泛的多

浏览 0 更新 2025-11-07

digamma函数

digamma函数(digamma function),又称Ψ函数(psi function),是伽马函数(gamma function)的对数导数,定义为:

ψ(z)=ddzlnΓ(z)=Γ(z)Γ(z)\psi(z) = \frac{d}{dz} \ln \Gamma(z) = \frac{\Gamma'(z)}{\Gamma(z)}

它是数学分析、数论和概率论中一类至关重要的特殊函数。由于伽马函数本身是阶乘向复平面的解析延拓,digamma函数自然成为研究级数求和、积分表示和差分方程的核心工具。它也是更广泛的多伽马函数(polygamma function)家族中的第一个成员。多伽马函数是伽马函数对数导数的各阶导数,记作 ψ(m)(z)\psi^{(m)}(z),其中 digamma 对应 m=0m = 0,一阶导数 ψ(1)(z)\psi^{(1)}(z) 称为三伽马函数(trigamma function),二阶及以上统称多伽马函数。

历史背景

digamma函数的研究可以追溯到18世纪。欧拉(Leonhard Euler)在研究伽马函数时已接触了其对数导数,但系统的研究由勒让德(Adrien-Marie Legendre)和高斯(Carl Friedrich Gauss)完成。高斯在1812年给出了digamma函数在有理点处的闭式表达式,即著名的高斯digamma定理。此后,维尔斯特拉斯(Karl Weierstrass)用无穷乘积给出了digamma函数的级数表示,奠定了其在复分析中的基础地位。20世纪以来,随着数值计算和统计学的蓬勃发展,digamma函数在实际应用中的价值日益凸显。

基本性质

digamma函数具有以下几组核心性质。

递推关系(recurrence relation):

ψ(z+1)=ψ(z)+1z\psi(z+1) = \psi(z) + \frac{1}{z}

该关系直接来源于伽马函数的函数方程 Γ(z+1)=zΓ(z)\Gamma(z+1) = z\Gamma(z)。对正整数 nn 反复应用可得:

ψ(n+1)=γ+k=1n1k\psi(n+1) = -\gamma + \sum_{k=1}^{n} \frac{1}{k}

其中 γ0.57721566\gamma \approx 0.57721566欧拉-马歇罗尼常数(Euler-Mascheroni constant)。这一关系表明,digamma 在正整数处的取值直接与调和数 HnH_n 挂钩。事实上,ψ(n+1)=Hnγ\psi(n+1) = H_n - \gamma

反射公式(reflection formula):

ψ(1z)ψ(z)=πcot(πz)\psi(1-z) - \psi(z) = \pi \cot(\pi z)

该式由伽马函数的反射公式 Γ(z)Γ(1z)=π/sin(πz)\Gamma(z)\Gamma(1-z) = \pi / \sin(\pi z) 两端取对数导数得出,巧妙地将 digamma 函数与三角函数联系起来。

加倍公式(duplication formula):

ψ(2z)=12ψ(z)+12ψ(z+12)+ln2\psi(2z) = \frac{1}{2} \psi(z) + \frac{1}{2} \psi\left(z + \frac{1}{2}\right) + \ln 2

这是更一般的乘法公式(multiplication formula)当 n=2n = 2 时的特例。对任意正整数 nn,乘法公式为:

ψ(nz)=1nk=0n1ψ(z+kn)+lnn\psi(nz) = \frac{1}{n} \sum_{k=0}^{n-1} \psi\left(z + \frac{k}{n}\right) + \ln n

级数表示与积分表示

digamma函数有多种等价的解析表示形式,不同形式适用于不同场景。

维尔斯特拉斯无穷乘积形式

ψ(z)=γ1z+n=1(1n1z+n)\psi(z) = -\gamma - \frac{1}{z} + \sum_{n=1}^{\infty} \left(\frac{1}{n} - \frac{1}{z + n}\right)

该形式直接揭示了 digamma 函数在负整数处的单极点结构,每个极点的留数均为 1-1

高斯积分表示

ψ(z)=0(ettezt1et)dt,Re(z)>0\psi(z) = \int_{0}^{\infty} \left( \frac{e^{-t}}{t} - \frac{e^{-zt}}{1 - e^{-t}} \right) dt, \quad \operatorname{Re}(z) > 0

该积分形式在解析延拓和渐近分析中十分有用。

渐近展开(当 z|z| \to \inftyargz<π|\arg z| < \pi):

ψ(z)=lnz12z112z2+1120z41252z6+O(1z8)\psi(z) = \ln z - \frac{1}{2z} - \frac{1}{12z^2} + \frac{1}{120z^4} - \frac{1}{252z^6} + O\left(\frac{1}{z^8}\right)

这一斯特林级数(Stirling series)使 digamma 函数在大参数下的快速计算成为可能。实际数值计算中,通常先用递推关系将自变量降到较小范围,再结合渐近展开达到高精度。

高斯digamma定理

高斯证明了在有理点处 digamma 函数可用初等函数和常数显式表示。设 0<p<q0 < p < q 为正整数,则:

ψ(pq)=γln(2q)π2cot(πpq)+2k=1(q1)/2cos(2πkpq)ln(sinπkq)\psi\left(\frac{p}{q}\right) = -\gamma - \ln(2q) - \frac{\pi}{2} \cot\left(\frac{\pi p}{q}\right) + 2\sum_{k=1}^{\lfloor (q-1)/2 \rfloor} \cos\left(\frac{2\pi k p}{q}\right) \ln\left(\sin\frac{\pi k}{q}\right)

这一结果在数论中的狄利克雷L函数(Dirichlet L-function)和类数公式(class number formula)研究中有重要应用。它表明尽管 digamma 一般不是初等函数,但在有理点上却可化为初等形式。

特殊值

digamma函数在几个特殊点处有简洁的闭式值:

ψ(1)=γ,ψ(12)=γ2ln2\psi(1) = -\gamma,\quad \psi\left(\frac{1}{2}\right) = -\gamma - 2\ln 2
ψ(13)=γ32ln3π23\psi\left(\frac{1}{3}\right) = -\gamma - \frac{3}{2}\ln 3 - \frac{\pi}{2\sqrt{3}}
ψ(23)=γ32ln3+π23\psi\left(\frac{2}{3}\right) = -\gamma - \frac{3}{2}\ln 3 + \frac{\pi}{2\sqrt{3}}
ψ(14)=γπ23ln2,ψ(34)=γ+π23ln2\psi\left(\frac{1}{4}\right) = -\gamma - \frac{\pi}{2} - 3\ln 2,\quad \psi\left(\frac{3}{4}\right) = -\gamma + \frac{\pi}{2} - 3\ln 2
ψ(16)=γ3π22ln232ln3\psi\left(\frac{1}{6}\right) = -\gamma - \frac{\sqrt{3}\pi}{2} - 2\ln 2 - \frac{3}{2}\ln 3

这些特殊值在理论物理和统计力学的精确计算中经常出现。

与调和数的关系

digamma函数与调和数 HnH_n 之间存在简洁而深刻的直接对应关系:

ψ(n+1)=Hnγ\psi(n+1) = H_n - \gamma

其中 Hn=1+12+13++1nH_n = 1 + \frac{1}{2} + \frac{1}{3} + \cdots + \frac{1}{n}。这一关系将连续的 digamma 函数与离散的调和数无缝衔接。利用该关系,可将涉及调和数的求和转化为 digamma 函数的积分或极限求解。例如,求解级数 n=1Hnn2\sum_{n=1}^{\infty} \frac{H_n}{n^2} 时,借助 digamma 函数的积分表示可得到闭式结果 2ζ(3)2\zeta(3)

应用领域

数论与解析数论:digamma函数出现在黎曼ζ函数(Riemann zeta function)的函数方程推导中,也在素数定理的渐近分析中发挥作用。高斯digamma定理为计算狄利克雷L函数在整数点处的值提供了直接工具。

统计与概率论:在狄利克雷分布(Dirichlet distribution)和维希特分布(Wishart distribution)的矩估计中,digamma函数及其导数被用于计算分布的期望和方差。在最大熵模型(maximum entropy model)的参数估计中,常常需要求解包含 digamma 函数的非线性方程组。此外,贝塔分布(beta distribution)的对数似然函数求导也会自然引出 digamma 函数。

物理与工程:在量子力学中,digamma函数出现在球对称势阱的能级计算和卡西米尔效应(Casimir effect)的数学表述中。在信号处理领域,它用于分析阿尔法稳定分布(alpha-stable distribution)的似然函数和参数估计。

数值计算:由于 digamma 函数没有初等闭式表达式(除有理点外),实际计算依赖渐近展开、递推关系和连分数算法。现代数学软件如 MATLAB、Mathematica、SciPy 均内置了 digamma 函数的高精度实现,底层通常结合斯特林级数与递推降阶策略,在保证精度的同时提升运算速度。

与其他特殊函数的联系

digamma函数与多个重要特殊函数有着深刻的联系。它是多伽马函数的基函数,通过求导可得到任意阶的多伽马函数。在广义欧拉常数的定义中,digamma 函数提供了统一的框架。它还与双伽马函数(harmonic number function)直接等价,调和数不过是 digamma 在整数点上的采样。在复分析中,digamma 函数的拉普拉斯变换(Laplace transform)连接了伽马函数与黎曼ζ函数,揭示了特殊函数之间的深层结构。

总结

digamma函数作为伽马函数的对数导数,是连接特殊函数论、组合数学和解析数论的桥梁。其递推关系、反射公式和高斯定理构成了完整的理论体系,而丰富的跨学科应用则使其在纯粹数学与应用科学中都占有不可替代的地位。无论是理论推导还是数值计算,深刻理解 digamma 函数的性质对于从事数学、统计和物理相关研究的工作者而言都是一项基本素养。它不仅是理论数学中的精妙构造,更是解决实际问题的有力工具,其理论深度与应用广度使其成为特殊函数论中一颗璀璨的明珠。