分析学
特殊函数
扩展初等函数、遍布数学与物理各处的伽玛函数、贝塔函数、zeta函数与贝塞尔函数。
直观超越阶乘:填补整数之间的空隙
阶乘 n!=1⋅2⋯n 本来只对整数有意义,但它增长得如此平滑——1,2,6,24,120,…——似乎在呼唤一条穿过每一个点、并在每个实数甚至复数输入处都有定义的曲线。特殊函数正是如此:这是一小族源自自然积分或无穷级数的函数,它们把初等算术(阶乘、幂、求和)的规律插值并推广为概率论、数论与物理学共享的平滑、普适的工具。
包含负实轴区域的区域着色图,伽玛函数在每个非正整数处产生一列单极点。大学伽玛函数、贝塔函数与zeta函数
定义: 伽玛函数
对 Re(z)>0,定义 Γ(z)=∫0∞tz−1e−tdt。分部积分可得 Γ(z+1)=zΓ(z),又因 Γ(1)=∫0∞e−tdt=1,由归纳法对任意非负整数 n 有 Γ(n+1)=n!:伽玛函数插值了阶乘。该递推关系也把 Γ 亚纯延拓到除 z=0,−1,−2,… 外的整个 z,在这些点处它有单极点。
Γ(z)=∫0∞tz−1e−tdt 贝塔函数 B(x,y)=∫01tx−1(1−t)y−1dt 可以表示为伽玛函数之比 B(x,y)=Γ(x+y)Γ(x)Γ(y),它是贝塔概率分布背后自然的归一化常数。一个更深的恒等式把二者对称地联系起来:欧拉反射公式 Γ(z)Γ(1−z)=sin(πz)π 对任意非整数 z 成立,把 Γ(z) 的极点与 sin(πz) 的零点联系起来。
Γ(z)Γ(1−z)=sin(πz)π 定义: 黎曼zeta函数及其解析延拓
对 Re(s)>1,定义 ζ(s)=n=1∑∞ns1。该级数在其他区域发散,但 ζ(s) 通过函数方程 ζ(s)=2sπs−1sin(2πs)Γ(1−s)ζ(1−s) 延拓为整个 C 上的亚纯函数,仅在 s=1 处有一个单极点;该函数方程把 ζ(s) 用 ζ(1−s) 与一个伽玛因子表示出来。正是这一点使人们能够严格地理解诸如 ζ(−1)=−121 这样的值,它与字面上求和 1+2+3+⋯ 毫无关系。
ζ(s)=2sπs−1sin(2πs)Γ(1−s)ζ(1−s) 四个特殊函数一览| 函数 | 定义公式 | 关键性质 |
|---|
| 伽玛, Γ(z) | Γ(z)=∫0∞tz−1e−tdt | 插值 n!;在 0,−1,−2,… 处有单极点 |
| 贝塔, B(x,y) | B(x,y)=∫01tx−1(1−t)y−1dt | B(x,y)=Γ(x+y)Γ(x)Γ(y) |
| Zeta, ζ(s) | ζ(s)=n=1∑∞ns1 | 亚纯延拓;在 s=1 处有一个极点 |
| 贝塞尔, Jν(x) | Jν(x)=k=0∑∞k!Γ(k+ν+1)(−1)k(2x)2k+ν | 满足 x2y′′+xy′+(x2−ν2)y=0 |
进阶关键定理
对任意非负整数 n,有 Γ(n+1)=n!,其中 Γ 对 Re(z)>0 由 Γ(z)=∫0∞tz−1e−tdt 定义。
为什么成立?
这表明伽玛积分并非随意的推广,而是唯一能再现阶乘乘法结构的光滑延拓,这正是它可以代入非整数的原因所在。
证明
第一步(基础情形)。Γ(1)=∫0∞t0e−tdt=∫0∞e−tdt=[−e−t]0∞=0−(−1)=1=0!,故公式在 n=0 时成立。
第二步(用分部积分得到递推关系)。对 Re(z)>0,取 u=tz、dv=e−tdt 对 Γ(z)=∫0∞tz−1e−tdt 分部积分:Γ(z+1)=∫0∞tze−tdt=[−tze−t]0∞+z∫0∞tz−1e−tdt=0+zΓ(z),这里用到边界项在两端都为零。这就证明了 Γ(z+1)=zΓ(z)。
第三步(归纳)。假设对某个非负整数 k 有 Γ(k+1)=k!。由第二步取 z=k+1,得 Γ(k+2)=(k+1)Γ(k+1)=(k+1)⋅k!=(k+1)!。
第四步(结论)。由于基础情形 n=0 成立(第一步),且归纳步骤把等式从 n=k 传递到 n=k+1(第三步),由数学归纳法,Γ(n+1)=n! 对任意非负整数 n 成立。
对满足 0<Re(z)<1 的任意 z,有 Γ(z)Γ(1−z)=sin(πz)π。
为什么成立?
它把一个看起来纯代数的对象(两个伽玛值之比)与三角函数联系起来,恰好是使 Γ 在非正整数处的极点与 sin(πz) 在所有整数处的零点对齐的恒等式——正是围道积分中的留数计算产生了右边的正弦函数。
证明
第一步(化为贝塔函数)。取 y=1−z 使用 B(x,y)=Γ(x+y)Γ(x)Γ(y),得 Γ(z)Γ(1−z)=B(z,1−z)=∫01tz−1(1−t)−zdt。代换 u=t/(1−t) 后,对 0<Re(z)<1,这变为 ∫0∞1+uuz−1du。
第二步(设置钥匙孔形围道)。为计算 I=∫0∞1+uuz−1du,取沿正实轴切开的 wz−1 分支,并沿钥匙孔形围道对 1+wwz−1 积分:紧贴切口上方向外、半径为 R 的大圆、紧贴切口下方返回,以及绕 0 的半径为 ε 的小圆。
第三步(唯一极点处的留数)。1+wwz−1 在围道内部唯一的奇点是 w=−1 处的单极点,在那里(用幅角为 π 的分支)Res=(−1)z−1=eiπ(z−1),故留数定理给出钥匙孔积分等于 2πieiπ(z−1)。
第四步(比较切口两侧,令 R→∞、ε→0)。当 0<Re(z)<1 时大圆和小圆的贡献消失,两条直线段的差别在于 wz−1 穿过切口时获得的相位,由此得到 I(1−e2πi(z−1))=2πieiπ(z−1)。解出 I 并利用 sinθ=2ieiθ−e−iθ 化简指数,得 I=sin(πz)π,结合第一步即得 Γ(z)Γ(1−z)=sin(πz)π。
进阶实际应用与典型例题
伽玛函数为用于建模等待时间和贝叶斯先验的概率分布(伽玛分布、贝塔分布、卡方分布、学生t分布)提供归一化常数;贝塞尔函数是柱面几何中波动方程和热方程的自然解,因此出现在鼓面声学、光学衍射图样和天线辐射方向图中;而zeta函数是解析数论(素数分布)的支柱,同时也通过zeta函数正规化出现在物理学中,用于给卡西米尔效应和弦论中的发散级数赋予有限值。
例题: 归一化卡方分布
某统计学家用卡方密度 f(x)=Cxk/2−1e−x/2(x>0)对 k=5 个独立标准正态变量的平方和建模。求归一化常数 C,使得 ∫0∞f(x)dx=1。
解答
第一步(代换以匹配伽玛积分)。令 t=x/2,则 x=2t,dx=2dt:∫0∞xk/2−1e−x/2dx=∫0∞(2t)k/2−1e−t2dt=2k/2∫0∞tk/2−1e−tdt。
第二步(识别出伽玛函数)。剩下的积分按定义恰好是 Γ(k/2),所以 ∫0∞xk/2−1e−x/2dx=2k/2Γ(k/2)。
第三步(代入 k=5)。此时 k/2=5/2,利用递推关系 Γ(5/2)=(3/2)(1/2)Γ(1/2)=43π,化简 25/2=42 后积分等于 42⋅43π=32π。
第四步(解出 C)。由于 C⋅32π=1,归一化常数为 C=32π1,与一般卡方公式 C=2k/2Γ(k/2)1 一致。
例题: 圆形鼓面的基频
一个半径为 a 的圆形鼓面具有径向对称的振动模式 u(r,t)=J0(kr)cos(ωt),它必须在刚性边缘处满足 u(a,t)=0(对所有 t)。已知 J0 的第一个正零点为 j0,1≈2.405,用 a 和波速 c(其中 ω=ck)表示基本角频率 ω1。
解答
第一步(应用边界条件)。由于 u(a,t)=J0(ka)cos(ωt) 必须对所有 t 都为零,需要 J0(ka)=0,所以 ka 必须是贝塞尔函数 J0 的一个零点。
第二步(选取基本模式)。基本模式对应最小的正 k,即 J0 的最小正零点 j0,1≈2.405;因此 k1a=j0,1,故 k1=j0,1/a。
第三步(转换为角频率)。利用膜波动方程的色散关系 ω=ck,得 ω1=ck1=acj0,1。
第四步(解释结果)。代入数值得 ω1≈a2.405c:半径 a 增大一倍,基本频率就减半,这与更大的鼓产生更低沉声音的物理直觉完全一致。
Γ(6) 等于多少?
Γ(21) 等于多少?
ζ(s) 的亚纯延拓唯一的极点在哪里,是什么类型?
一位声学工程师设计一个圆形鼓面,需要在中心保持有限、在刚性边缘处为零的模式形状。哪一族特殊函数给出这些径向模式形状?