5.2 留数理论

5.2 留数理论

本节摘要：SOURCE 5.2：“5.2.1 留数定义与计算方法”将深入各类奇点的计算技艺；“5.2.2 留数定理”将严格证明其普适性，并探讨其在多连通区域、无穷远点等广义情形下的推广；“5.2.3 实际积分计算”则将展示其作为“万能钥匙”的磅礴力量，从三角积分、有理函数积分，到含参数的积分、反常积分，乃至椭圆积分的复变解法。

留数的定义：被精心标记的系数

5.1 节确立了一个判据：奇点类型由洛朗展开的主部形状决定。本节要提取主部中那个"最值钱"的系数。设 f(z) 在孤立奇点 z_0 的穿孔邻域 0<|z-z_0|<R 内解析，其洛朗展开为

则系数 c_{-1}（即 \frac{1}{z-z_0} 项的系数）称为 f 在 z_0 处的留数，记作 \operatorname{Res}(f,z_0)。

留数不是级数中偶然被标记的系数，而是围道积分的不变量。由柯西积分公式直接可得：对任意包围且仅包围 z_0 的正向简单闭曲线 \gamma，

这个等式揭示了一个惊人的事实：无论 \gamma 如何扭曲、如何贴近或远离 z_0，只要它同伦于一个以 z_0 为中心的小圆周而不跨过其他奇点，积分值除以 2\pi i 后恒等于同一个复数 c_{-1}。留数是奇点所携带的"拓扑电荷"——它不关心奇点内部如何狂乱，只忠实地编码围绕它的"旋转总量"。

留数的计算方法

计算方法的选择，本质上是对奇点类型的认知策略：

1. 可去奇点：主部为空，c_{-1}=0，留数恒为零。

2. 极点：若 z_0 是 m 阶极点，则有通用公式

特别地，一阶极点（m=1）时公式化简为 \operatorname{Res}(f,z_0)=\lim_{z\to z_0}(z-z_0)f(z)。若 f=p/q 而 q(z_0)=0,\ q'(z_0)\ne0，则还有更常用的形式 \operatorname{Res}(f,z_0)=p(z_0)/q'(z_0)。

3. 本性奇点：无通用求导公式，必须回归洛朗展开本身。例如 f(z)=e^{1/z} 的展开为 \sum_{n=0}^\infty\frac{z^{-n}}{n!}，故 c_{-1}=\frac{1}{1!}=1。

留数定理：复平面的积分守恒律

留数定理：设 f(z) 在单连通区域 D 内除有限个孤立奇点 z_1,\dots,z_n 外处处解析，\gamma 是 D 内一条分段光滑的正向简单闭曲线且其内部包含全部这些奇点，则

这个公式美得令人屏息：一个看似复杂的围道积分，其值竟完全由几个离散点上的局部信息线性叠加决定。路径的形状、长度、曲折程度——一切几何细节都退居次要，真正主宰结果的是被包围奇点所携带的"拓扑电荷"。

证明思路用"挖洞法"：在每个 z_k 周围作一个极小正向圆周 \gamma_k，使所有 \gamma_k 互不相交且均在 \gamma 内部。由复连通区域的柯西定理，

留数定理的伟大在于将全局积分问题降维为局部奇点分析问题。它还可推广到多连通区域（只需保证所有奇点位于 \gamma 所围成的有界区域内）以及无穷远点：若 f 在扩充复平面上除有限个奇点外解析，则全部留数（含无穷远点处的）之和为零。

用留数定理计算实积分

留数定理最震撼的应用，是为实变函数中棘手的定积分与反常积分提供统一解法。其核心思想是维度跃迁：把一维实轴上的积分嵌入二维复平面，借助围道积分将难以捉摸的实积分转化为可精确计算的离散留数和。

例 1（有理函数在无穷区间）：计算 \int_{-\infty}^{\infty}\frac{dx}{1+x^2}。

取 F(z)=\frac{1}{1+z^2}，它在全平面除 z=\pm i 外解析。对 z=i（上半平面内的单极点），

取"实轴段 [-R,R] + 上半圆弧 \Gamma_R"构成的闭合围道，当 R\to\infty 时上半圆弧上的积分趋于零（因为 |F(z)|\sim 1/|z|^2，弧长 \sim\pi R，乘积 \to0）。故

例 2（三角积分）：计算 I=\int_0^{2\pi}\frac{d\theta}{a+\cos\theta}，a>1。令 z=e^{i\theta}，则 d\theta=dz/(iz)、\cos\theta=\frac12(z+z^{-1})，转化为单位圆上的围道积分

被积函数在单位圆内的极点是 z=-a+\sqrt{a^2-1}（另一个根 z=-a-\sqrt{a^2-1} 在圆外），计算留数后得 I=\frac{2\pi}{\sqrt{a^2-1}}。

例 3（Fourier 型积分）：\int_{-\infty}^{\infty}\frac{e^{i\xi x}}{x^2+a^2}\,dx 可用 e^{i\xi z} 在上半平面的指数衰减（\operatorname{Im}z>0 时 |e^{i\xi z}|=e^{-\xi\operatorname{Im}z}\to0）控制辅助弧，留数方法同样有效。这揭示了指数衰减函数的 Fourier 变换为有理函数，极点位置直接编码原函数的衰减速率。

数值验证：围道积分与留数

import numpy as np

# 验证 Res(f, i) = 1/(2i)：数值围道积分 f(z)=1/(1+z^2) 绕 z=i
r = 0.4                                   # 小圆半径，绕 z=i 且不包含 z=-i
t = np.linspace(0, 2*np.pi, 100000, endpoint=False)
z = 1j + r*np.exp(1j*t)
dz = 1j*r*np.exp(1j*t)
f = 1.0 / (1.0 + z**2)
I = np.trapz(f*dz, t)                     # 数值积分 (1/(2πi))∮f dz
res = I / (2*np.pi*1j)
print("数值留数: %.8f %+.8fi (理论 1/(2i) = -0.5i)" % (res.real, res.imag))
​

# 验证留数定理：∮_{|z|=2} dz/(z^2-1) = 2πi * [Res(z=1)+Res(z=-1)]
# 两个单极点均在圆 |z|=2 内，留数分别为 1/2 和 -1/2，和为零
r = 2.0
t = np.linspace(0, 2*np.pi, 200000, endpoint=False)
z = r*np.exp(1j*t)
dz = 1j*r*np.exp(1j*t)
I = np.trapz((1.0/(z**2-1))*dz, t)
print("∮ dz/(z^2-1) = %.6f %+.6fi" % (I.real, I.imag))
print("2πi·(1/2 - 1/2) = 0.000000 +0.000000i")
​

# 验证例 2 的结果：∫_0^{2π} dθ/(a+cosθ) = 2π/√(a^2-1)
a = 3.0
th = np.linspace(0, 2*np.pi, 400000, endpoint=False)
numer = np.trapz(1.0/(a+np.cos(th)), th)
theory = 2*np.pi/np.sqrt(a*a-1)
print("数值积分: %.8f, 留数理论: %.8f, 误差 %.2e"
      % (numer, theory, abs(numer-theory)))
​

从局部到全局：留数的现代回响

留数理论的影响远超积分计算本身。在量子场论中，费曼传播子的计算本质是求解含 i\varepsilon 规则的复积分；在信号处理中，滤波器响应稳定性分析依赖系统函数极点是否位于左半平面；在空气动力学中，Kutta-Joukowski 定理把机翼升力归结为复速度势绕翼型的环量，其核心仍是留数。代数视角下，留数是洛朗环上的线性泛函，在代数几何中被推广为 Serre 对偶理论的核心构件——留数映射；拓扑视角下，围道积分的值只依赖曲线的缠绕数与奇点留数，即 f\,dz 诱导的上同调类在生成元上的取值。

⚠️ 常见坑：只记结论不记适用边界——留数定理要求奇点是孤立的且不在围道上；奇点位于 \gamma 上时积分可能发散，定理失效。无穷远点处的留数常被遗漏，检查"全部留数之和为零"可作自检。

💡 关键直觉：围道积分的值等于 2\pi i 乘以内部所有奇点留数之和——这是贯穿积分计算的核心杠杆；三个数值算例分别核对了留数本身、留数定理求和与三角积分公式。

本节延伸：留数求和的应用

留数定理还能计算无穷级数。若 f 在实轴上无极点且满足适当衰减条件，则

其中 \pi\cot\pi z 在全体整数处以留数 1 的极点为"采样器"。这一技巧把离散求和转化为有限个留数的计算，是解析数论中估和的标准工具。本节掌握的核心能力：识别奇点类型 → 选择留数公式 → 构造围道 → 取极限，四个步骤缺一不可，而其正确性可用数值围道积分随时核对。