MathLabs

分析学

特殊函数

扩展初等函数、遍布数学与物理各处的伽玛函数、贝塔函数、zeta函数与贝塞尔函数。

直观超越阶乘:填补整数之间的空隙

阶乘 n!=1⋅2⋯nn! = 1\cdot 2 \cdots n 本来只对整数有意义,但它增长得如此平滑——1,2,6,24,120,…1, 2, 6, 24, 120, \dots——似乎在呼唤一条穿过每一个点、并在每个实数甚至复数输入处都有定义的曲线。特殊函数正是如此:这是一小族源自自然积分或无穷级数的函数,它们把初等算术(阶乘、幂、求和)的规律插值并推广为概率论、数论与物理学共享的平滑、普适的工具。

彩色复平面图,显示一系列沿负实轴排列的极点。
包含负实轴区域的区域着色图,伽玛函数在每个非正整数处产生一列单极点。

大学伽玛函数、贝塔函数与zeta函数

定义: 伽玛函数

对 Re⁡(z)>0\operatorname{Re}(z) > 0,定义 Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt。分部积分可得 Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z),又因 Γ(1)=∫0∞e−t dt=1\Gamma(1) = \int_0^\infty e^{-t}\,dt = 1,由归纳法对任意非负整数 nn 有 Γ(n+1)=n!\Gamma(n+1)=n!:伽玛函数插值了阶乘。该递推关系也把 Γ\Gamma 亚纯延拓到除 z=0,−1,−2,…z = 0, -1, -2, \dots 外的整个 zz,在这些点处它有单极点。

Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt

贝塔函数 B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dt 可以表示为伽玛函数之比 B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)},它是贝塔概率分布背后自然的归一化常数。一个更深的恒等式把二者对称地联系起来:欧拉反射公式 Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)} 对任意非整数 zz 成立,把 Γ(z)\Gamma(z) 的极点与 sin⁡(πz)\sin(\pi z) 的零点联系起来。

Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}

定义: 黎曼zeta函数及其解析延拓

对 Re⁡(s)>1\operatorname{Re}(s) > 1,定义 ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s}。该级数在其他区域发散,但 ζ(s)\zeta(s) 通过函数方程 ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s) 延拓为整个 C\mathbb{C} 上的亚纯函数,仅在 s=1s=1 处有一个单极点;该函数方程把 ζ(s)\zeta(s) 用 ζ(1−s)\zeta(1-s) 与一个伽玛因子表示出来。正是这一点使人们能够严格地理解诸如 ζ(−1)=−112\zeta(-1) = -\dfrac{1}{12} 这样的值,它与字面上求和 1+2+3+⋯1+2+3+\cdots 毫无关系。

ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s)
四个特殊函数一览
函数定义公式关键性质
伽玛, Γ(z)\Gamma(z)Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt插值 n!n!;在 0,−1,−2,…0,-1,-2,\dots 处有单极点
贝塔, B(x,y)B(x,y)B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dtB(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}
Zeta, ζ(s)\zeta(s)ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s}亚纯延拓;在 s=1s=1 处有一个极点
贝塞尔, Jν(x)J_\nu(x)Jν(x)=∑k=0∞(−1)kk! Γ(k+ν+1)(x2)2k+νJ_\nu(x) = \displaystyle\sum_{k=0}^{\infty} \dfrac{(-1)^k}{k!\,\Gamma(k+\nu+1)} \left(\dfrac{x}{2}\right)^{2k+\nu}满足 x2y′′+xy′+(x2−ν2)y=0x^2y''+xy'+(x^2-\nu^2)y=0

进阶关键定理

对任意非负整数 nn,有 Γ(n+1)=n!\Gamma(n+1)=n!,其中 Γ\Gamma 对 Re⁡(z)>0\operatorname{Re}(z) > 0 由 Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt 定义。

为什么成立?

这表明伽玛积分并非随意的推广,而是唯一能再现阶乘乘法结构的光滑延拓,这正是它可以代入非整数的原因所在。

证明

第一步(基础情形)。Γ(1)=∫0∞t0e−t dt=∫0∞e−t dt=[−e−t]0∞=0−(−1)=1=0!\Gamma(1) = \int_0^\infty t^0 e^{-t}\,dt = \int_0^\infty e^{-t}\,dt = \big[-e^{-t}\big]_0^\infty = 0-(-1) = 1 = 0!,故公式在 n=0n=0 时成立。

第二步(用分部积分得到递推关系)。对 Re⁡(z)>0\operatorname{Re}(z) > 0,取 u=tzu=t^z、dv=e−t dtdv=e^{-t}\,dt 对 Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt 分部积分:Γ(z+1)=∫0∞tze−t dt=[−tze−t]0∞+z∫0∞tz−1e−t dt=0+z Γ(z)\Gamma(z+1) = \int_0^\infty t^z e^{-t}\,dt = \big[-t^z e^{-t}\big]_0^\infty + z\int_0^\infty t^{z-1}e^{-t}\,dt = 0 + z\,\Gamma(z),这里用到边界项在两端都为零。这就证明了 Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z)。

第三步(归纳)。假设对某个非负整数 kk 有 Γ(k+1)=k!\Gamma(k+1) = k!。由第二步取 z=k+1z=k+1,得 Γ(k+2)=(k+1) Γ(k+1)=(k+1)⋅k!=(k+1)!\Gamma(k+2) = (k+1)\,\Gamma(k+1) = (k+1)\cdot k! = (k+1)!。

第四步(结论)。由于基础情形 n=0n=0 成立(第一步),且归纳步骤把等式从 n=kn=k 传递到 n=k+1n=k+1(第三步),由数学归纳法,Γ(n+1)=n!\Gamma(n+1)=n! 对任意非负整数 nn 成立。

对满足 0<Re⁡(z)<10 < \operatorname{Re}(z) < 1 的任意 zz,有 Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}。

为什么成立?

它把一个看起来纯代数的对象(两个伽玛值之比)与三角函数联系起来,恰好是使 Γ\Gamma 在非正整数处的极点与 sin⁡(πz)\sin(\pi z) 在所有整数处的零点对齐的恒等式——正是围道积分中的留数计算产生了右边的正弦函数。

证明

第一步(化为贝塔函数)。取 y=1−zy=1-z 使用 B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)},得 Γ(z)Γ(1−z)=B(z,1−z)=∫01tz−1(1−t)−z dt\Gamma(z)\Gamma(1-z) = B(z,1-z) = \int_0^1 t^{z-1}(1-t)^{-z}\,dt。代换 u=t/(1−t)u = t/(1-t) 后,对 0<Re⁡(z)<10<\operatorname{Re}(z)<1,这变为 ∫0∞uz−11+u du\displaystyle\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du。

第二步(设置钥匙孔形围道)。为计算 I=∫0∞uz−11+u duI=\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du,取沿正实轴切开的 wz−1w^{z-1} 分支,并沿钥匙孔形围道对 wz−11+w\dfrac{w^{z-1}}{1+w} 积分:紧贴切口上方向外、半径为 RR 的大圆、紧贴切口下方返回,以及绕 00 的半径为 ε\varepsilon 的小圆。

第三步(唯一极点处的留数)。wz−11+w\dfrac{w^{z-1}}{1+w} 在围道内部唯一的奇点是 w=−1w=-1 处的单极点,在那里(用幅角为 π\pi 的分支)Res=(−1)z−1=eiπ(z−1)\mathrm{Res} = (-1)^{z-1} = e^{i\pi(z-1)},故留数定理给出钥匙孔积分等于 2πi eiπ(z−1)2\pi i\, e^{i\pi(z-1)}。

第四步(比较切口两侧,令 R→∞R\to\infty、ε→0\varepsilon\to 0)。当 0<Re⁡(z)<10<\operatorname{Re}(z)<1 时大圆和小圆的贡献消失,两条直线段的差别在于 wz−1w^{z-1} 穿过切口时获得的相位,由此得到 I(1−e2πi(z−1))=2πi eiπ(z−1)I\big(1 - e^{2\pi i(z-1)}\big) = 2\pi i\, e^{i\pi(z-1)}。解出 II 并利用 sin⁡θ=eiθ−e−iθ2i\sin\theta = \frac{e^{i\theta}-e^{-i\theta}}{2i} 化简指数,得 I=πsin⁡(πz)I = \dfrac{\pi}{\sin(\pi z)},结合第一步即得 Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}。

进阶实际应用与典型例题

伽玛函数为用于建模等待时间和贝叶斯先验的概率分布(伽玛分布、贝塔分布、卡方分布、学生t分布)提供归一化常数;贝塞尔函数是柱面几何中波动方程和热方程的自然解,因此出现在鼓面声学、光学衍射图样和天线辐射方向图中;而zeta函数是解析数论(素数分布)的支柱,同时也通过zeta函数正规化出现在物理学中,用于给卡西米尔效应和弦论中的发散级数赋予有限值。

例题: 归一化卡方分布

某统计学家用卡方密度 f(x)=C xk/2−1e−x/2f(x) = C\, x^{k/2-1}e^{-x/2}(x>0x>0)对 k=5k=5 个独立标准正态变量的平方和建模。求归一化常数 CC,使得 ∫0∞f(x) dx=1\int_0^\infty f(x)\,dx = 1。

解答

第一步(代换以匹配伽玛积分)。令 t=x/2t = x/2,则 x=2tx=2t,dx=2 dtdx = 2\,dt:∫0∞xk/2−1e−x/2 dx=∫0∞(2t)k/2−1e−t 2 dt=2k/2∫0∞tk/2−1e−t dt\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = \int_0^\infty (2t)^{k/2-1}e^{-t}\, 2\,dt = 2^{k/2}\int_0^\infty t^{k/2-1}e^{-t}\,dt。

第二步(识别出伽玛函数)。剩下的积分按定义恰好是 Γ(k/2)\Gamma(k/2),所以 ∫0∞xk/2−1e−x/2 dx=2k/2 Γ(k/2)\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = 2^{k/2}\,\Gamma(k/2)。

第三步(代入 k=5k=5)。此时 k/2=5/2k/2 = 5/2,利用递推关系 Γ(5/2)=(3/2)(1/2)Γ(1/2)=34π\Gamma(5/2) = (3/2)(1/2)\Gamma(1/2) = \dfrac{3}{4}\sqrt{\pi},化简 25/2=422^{5/2}=4\sqrt{2} 后积分等于 42⋅34π=32π4\sqrt{2}\cdot\frac{3}{4}\sqrt{\pi}=3\sqrt{2\pi}。

第四步(解出 CC)。由于 C⋅32π=1C \cdot 3\sqrt{2\pi} = 1,归一化常数为 C=132πC = \dfrac{1}{3\sqrt{2\pi}},与一般卡方公式 C=12k/2Γ(k/2)C = \dfrac{1}{2^{k/2}\Gamma(k/2)} 一致。

例题: 圆形鼓面的基频

一个半径为 aa 的圆形鼓面具有径向对称的振动模式 u(r,t)=J0(kr)cos⁡(ωt)u(r,t) = J_0(kr)\cos(\omega t),它必须在刚性边缘处满足 u(a,t)=0u(a,t)=0(对所有 tt)。已知 J0J_0 的第一个正零点为 j0,1≈2.405j_{0,1}\approx 2.405,用 aa 和波速 cc(其中 ω=ck\omega = ck)表示基本角频率 ω1\omega_1。

解答

第一步(应用边界条件)。由于 u(a,t)=J0(ka)cos⁡(ωt)u(a,t)=J_0(ka)\cos(\omega t) 必须对所有 tt 都为零,需要 J0(ka)=0J_0(ka) = 0,所以 kaka 必须是贝塞尔函数 J0J_0 的一个零点。

第二步(选取基本模式)。基本模式对应最小的正 kk,即 J0J_0 的最小正零点 j0,1≈2.405j_{0,1}\approx 2.405;因此 k1a=j0,1k_1 a = j_{0,1},故 k1=j0,1/ak_1 = j_{0,1}/a。

第三步(转换为角频率)。利用膜波动方程的色散关系 ω=ck\omega = ck,得 ω1=ck1=c j0,1a\omega_1 = c k_1 = \dfrac{c\, j_{0,1}}{a}。

第四步(解释结果)。代入数值得 ω1≈2.405 ca\omega_1 \approx \dfrac{2.405\,c}{a}:半径 aa 增大一倍,基本频率就减半,这与更大的鼓产生更低沉声音的物理直觉完全一致。

Γ(6)\Gamma(6) 等于多少?

Γ ⁣(12)\Gamma\!\left(\dfrac{1}{2}\right) 等于多少?

ζ(s)\zeta(s) 的亚纯延拓唯一的极点在哪里,是什么类型?

一位声学工程师设计一个圆形鼓面,需要在中心保持有限、在刚性边缘处为零的模式形状。哪一族特殊函数给出这些径向模式形状?

参考文献

  1. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 5: Gamma Function
  2. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 25: Zeta and Related Functions
  3. David J. Platt, Timothy S. Trudgian (2021). The Riemann Hypothesis Is True Up to 3×10^12 · arXiv:2004.09765