MathLabs

分析学

留数与留数定理

通过对围道内部极点处的留数求和来计算复变周线积分的方法。

直观从孤立极点到无穷和

想象在复平面上沿一条闭合回路 γ\gamma 行走,绕过某个在若干孤立点处发散的函数 ff。在每个这样的点附近,ff 的表现就像螺旋楼梯:越靠近,数值旋转得越快、增长得越猛。留数定理揭示了一个惊人的事实:整个回路积分的值只取决于内部每个被围住的点处提取出的一个数——它的留数——而与回路的具体形状无关,只要回路不穿过任何奇点。

复变函数的彩色图,显示围绕孤立极点旋转的颜色图案。
带孤立极点函数的区域着色图;颜色围绕极点旋转,这种旋转编码了该处的留数。

大学孤立奇点与洛朗级数

定义: 孤立奇点及其留数

若 ff 在 z0z_0 周围某个去心圆盘上全纯,但在 z0z_0 处本身不全纯,则称 z0z_0 为 ff 的孤立奇点。在该去心圆盘上,ff 有唯一的洛朗展开 f(z)=∑n=−∞∞an(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} a_n (z-z_0)^n,而项 (z−z0)−1(z-z_0)^{-1} 的系数 a−1a_{-1} 称为 ff 在 z0z_0 处的留数,记作 Res(f,z0)=a−1\mathrm{Res}(f, z_0) = a_{-1}。奇点按出现的负幂个数分类:没有负幂是可去奇点(极限 lim⁡z→z0f(z)\lim_{z \to z_0} f(z) 存在),有限个负幂、最高阶为 mm 是极点,无穷多个负幂(如在 z=0z=0 处的 sin⁡(1/z)\sin(1/z))是本性奇点。

f(z)=∑n=−∞∞an(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} a_n (z-z_0)^n

实际操作中,没有人会为了求某个极点处的留数而展开整个洛朗级数。对于单极点(阶 m=1m=1),有 Res(f,z0)=lim⁡z→z0(z−z0)f(z)\mathrm{Res}(f, z_0) = \lim_{z \to z_0} (z-z_0) f(z)。对于更高阶的极点 mm,公式推广为 Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\mathrm{Res}(f, z_0) = \dfrac{1}{(m-1)!} \lim_{z \to z_0} \dfrac{d^{m-1}}{dz^{m-1}} \left[ (z-z_0)^m f(z) \right],该式在代入之前先通过求导消去奇异部分。

Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\mathrm{Res}(f, z_0) = \dfrac{1}{(m-1)!} \lim_{z \to z_0} \dfrac{d^{m-1}}{dz^{m-1}} \left[ (z-z_0)^m f(z) \right]
孤立奇点的分类及其留数计算
奇点类型洛朗级数留数公式
可去奇点无负幂项Res(f,z0)=0\mathrm{Res}(f, z_0) = 0
单极点(m=1m=1)一个负幂项Res(f,z0)=lim⁡z→z0(z−z0)f(z)\mathrm{Res}(f, z_0) = \lim_{z \to z_0} (z-z_0) f(z)
mm 阶极点有限个负幂项Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\mathrm{Res}(f, z_0) = \dfrac{1}{(m-1)!} \lim_{z \to z_0} \dfrac{d^{m-1}}{dz^{m-1}} \left[ (z-z_0)^m f(z) \right]
本性奇点无穷多个负幂项直接读取 a−1a_{-1}

大学关键定理

定理: 留数定理

设 ff 在单连通区域上除有限个孤立奇点 z1,…,znz_1, \dots, z_n 外全纯,这些奇点都在一条正向简单闭曲线 γ\gamma 内部。则 ∮γf(z) dz=2πi∑kRes(f,zk)\oint_\gamma f(z)\,dz = 2\pi i \sum_{k} \mathrm{Res}(f, z_k)。

为什么成立?

它把一个困难的几何问题(沿曲线积分)化归为一个简单的代数问题(把有限个数加起来),因为将围道在每个极点附近变形后会收缩成一个很小的圆,在那里洛朗级数完成了全部工作。

证明

第一步(变形围道)。根据柯西积分定理,对任何围住 ff 全纯区域的闭曲线都有 ∮γf(z) dz=0\oint_\gamma f(z)\,dz = 0。由于 ff 仅在 z1,…,znz_1,\dots,z_n 处不全纯,在每个 zkz_k 周围画一个半径为 ε\varepsilon 的正向小圆 CkC_k,半径足够小使这些圆彼此不相交且都在 γ\gamma 内部。从 γ\gamma 到每个 CkC_k 切开细缝,把 γ\gamma 与各 CkC_k 之间的区域变成一个 ff 全纯的单连通区域,因此该区域边界上的积分为 00;细缝两侧的贡献成对抵消,剩下 ∮γf(z) dz=∑k∮Ckf(z) dz\oint_\gamma f(z)\,dz = \sum_{k} \oint_{C_k} f(z)\,dz。

第二步(计算每个小圆)。固定 kk,在包含 CkC_k 的去心圆盘上,将 ff 在 zkz_k 附近展开为洛朗级数 f(z)=∑n=−∞∞an(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} a_n (z-z_0)^n。对满足 n≠−1n \neq -1 的每一项 an(z−zk)na_n(z-z_k)^n,它在该去心圆盘上都有单值原函数,故沿闭合小圆 CkC_k 积分为 00;只有项 a−1(z−zk)−1a_{-1}(z-z_k)^{-1} 保留下来。

第三步(计算剩余积分)。将 CkC_k 参数化为 z=zk+εeiθz = z_k + \varepsilon e^{i\theta}(θ∈[0,2π]\theta \in [0, 2\pi]),于是 dz=iεeiθ dθdz = i\varepsilon e^{i\theta}\,d\theta 且 dzz−zk=i dθ\dfrac{dz}{z-z_k} = i\,d\theta。因此 ∮Cka−1z−zk dz=a−1∫02πi dθ=2πi a−1=2πi Res(f,zk)\oint_{C_k} \dfrac{a_{-1}}{z-z_k}\,dz = a_{-1} \int_0^{2\pi} i\,d\theta = 2\pi i\, a_{-1} = 2\pi i\, \mathrm{Res}(f, z_k)。

第四步(求和)。将第三步代入第一步的等式,得到 ∮γf(z) dz=2πi∑kRes(f,zk)\oint_\gamma f(z)\,dz = 2\pi i \sum_{k} \mathrm{Res}(f, z_k),这正是留数定理。

若 ff 在 z0z_0 处有 mm 阶极点,即 g(z)=(z−z0)mf(z)g(z) = (z-z_0)^m f(z) 可全纯延拓到 z0z_0 且 g(z0)≠0g(z_0) \neq 0,则 Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\mathrm{Res}(f, z_0) = \dfrac{1}{(m-1)!} \lim_{z \to z_0} \dfrac{d^{m-1}}{dz^{m-1}} \left[ (z-z_0)^m f(z) \right]。特别地,对单极点(m=1m=1),有 Res(f,z0)=lim⁡z→z0(z−z0)f(z)\mathrm{Res}(f, z_0) = \lim_{z \to z_0} (z-z_0) f(z)。

为什么成立?

乘以 (z−z0)m(z-z_0)^m 会消去洛朗级数的奇异部分,使其变成系数恰为导数的普通泰勒级数——于是提取 a−1a_{-1} 就变成了一次普通的微积分计算。

证明

第一步(明确写出洛朗级数)。由于 ff 在 z0z_0 处有 mm 阶极点,其洛朗级数形如 f(z)=a−m(z−z0)m+⋯+a−1z−z0+a0+a1(z−z0)+⋯f(z) = \dfrac{a_{-m}}{(z-z_0)^m} + \cdots + \dfrac{a_{-1}}{z-z_0} + a_0 + a_1(z-z_0) + \cdots,其中 a−m≠0a_{-m} \neq 0。

第二步(消去极点)。两边乘以 (z−z0)m(z-z_0)^m:g(z):=(z−z0)mf(z)=a−m+a−m+1(z−z0)+⋯+a−1(z−z0)m−1+a0(z−z0)m+⋯g(z) := (z-z_0)^m f(z) = a_{-m} + a_{-m+1}(z-z_0) + \cdots + a_{-1}(z-z_0)^{m-1} + a_0(z-z_0)^m + \cdots,这是一个普通幂级数,因此 gg 在 z0z_0 处全纯。

第三步(把系数识别为导数)。在此幂级数中,a−1a_{-1} 是 (z−z0)m−1(z-z_0)^{m-1} 的系数。对任意全纯函数 g(z)=∑jcj(z−z0)jg(z) = \sum_j c_j (z-z_0)^j,泰勒公式给出 cj=g(j)(z0)j!c_j = \dfrac{g^{(j)}(z_0)}{j!},因此 a−1=g(m−1)(z0)(m−1)!a_{-1} = \dfrac{g^{(m-1)}(z_0)}{(m-1)!}。

第四步(结论)。由于 Res(f,z0)=a−1\mathrm{Res}(f,z_0) = a_{-1} 且 g(m−1)(z0)=lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]g^{(m-1)}(z_0) = \lim_{z\to z_0} \dfrac{d^{m-1}}{dz^{m-1}}[(z-z_0)^m f(z)],代入即得 Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\mathrm{Res}(f, z_0) = \dfrac{1}{(m-1)!} \lim_{z \to z_0} \dfrac{d^{m-1}}{dz^{m-1}} \left[ (z-z_0)^m f(z) \right]。令 m=1m=1 化简导数与阶乘,即得 Res(f,z0)=lim⁡z→z0(z−z0)f(z)\mathrm{Res}(f, z_0) = \lim_{z \to z_0} (z-z_0) f(z)。

大学实际应用与典型例题

留数把原本难以处理的实定积分变成了简单的算术运算,这正是它们广泛出现在物理学和工程学中的原因:在信号处理中计算振荡的傅里叶和拉普拉斯变换、在量子场论中计算散射振幅和传播子、通过拉普拉斯逆变换求解控制系统的微分方程,以及在解析数论和组合学中估计计数函数的渐近行为。

例题: 没有初等原函数捷径的实积分

用留数定理计算 ∫−∞∞dxx2+1\displaystyle\int_{-\infty}^{\infty} \frac{dx}{x^2+1}。(该积分也可用代换 x=tan⁡θx=\tan\theta 求得,但留数方法能推广到代换法无法处理的积分。)

解答

第一步(建立围道)。考虑 f(z)=1z2+1f(z) = \dfrac{1}{z^2+1} 以及由实轴上的线段 [−R,R][-R,R] 与逆时针方向的上半圆 ∣z∣=R|z|=R 组成的闭合围道 γR\gamma_R(RR 取较大值)。

第二步(估计圆弧部分)。在半圆弧上,∣f(z)∣≤1R2−1|f(z)| \le \dfrac{1}{R^2-1},而弧长为 πR\pi R,所以弧上的贡献至多为 πRR2−1→0\dfrac{\pi R}{R^2-1} \to 0(当 R→∞R \to \infty 时)。

第三步(找出围道内的极点)。f(z)=1z2+1f(z) = \dfrac{1}{z^2+1} 在 z=iz = i 和 z=−iz=-i 处有单极点;只有 z=iz = i 位于上半圆内部。由单极点公式,Res(f,i)=lim⁡z→i(z−i)1(z−i)(z+i)=12i\mathrm{Res}(f,i) = \lim_{z\to i}(z-i)\dfrac{1}{(z-i)(z+i)} = \dfrac{1}{2i}。

第四步(应用留数定理并取极限)。对任意 R>1R>1 都有 ∮γRf(z) dz=2πi⋅12i=π\oint_{\gamma_R} f(z)\,dz = 2\pi i \cdot \dfrac{1}{2i} = \pi,这等于实积分加上消失的圆弧项,令 R→∞R \to \infty 即得 ∫−∞∞dxx2+1=π\displaystyle\int_{-\infty}^{\infty} \frac{dx}{x^2+1} = \pi。

例题: 控制理论中的拉普拉斯逆变换

某控制工程师用传递函数 F(s)=1(s+1)(s+2)2F(s) = \dfrac{1}{(s+1)(s+2)^2} 对一个有阻尼的系统建模,需要求出其脉冲响应 y(t)y(t)(即拉普拉斯逆变换),以检验系统是否能不振荡地趋于稳定。

解答

第一步(回顾反演积分)。拉普拉斯逆变换为 y(t)=12πi∫c−i∞c+i∞F(s)est dsy(t) = \dfrac{1}{2\pi i}\int_{c-i\infty}^{c+i\infty} F(s) e^{st}\,ds,对于在无穷远处衰减的有理函数 F(s)F(s),该积分等于 F(s)=1(s+1)(s+2)2F(s) = \dfrac{1}{(s+1)(s+2)^2} 全部极点处 F(s)estF(s)e^{st} 的留数之和乘以 2πi2\pi i。

第二步(确定极点位置)。F(s)=1(s+1)(s+2)2F(s) = \dfrac{1}{(s+1)(s+2)^2} 在 s=−1s=-1 处有一个单极点,在 s=−2s=-2 处有一个 22 阶极点。

第三步(单极点处的留数)。Ress=−1 F(s)est=lim⁡s→−1(s+1)est(s+1)(s+2)2=e−t(−1)2=e−t\mathrm{Res}_{s=-1}\, F(s)e^{st} = \lim_{s\to -1}(s+1)\dfrac{e^{st}}{(s+1)(s+2)^2} = \dfrac{e^{-t}}{(-1)^2} = e^{-t}。

第四步(二阶极点处的留数)。Ress=−2 F(s)est=lim⁡s→−2dds[ests+1]=lim⁡s→−2test(s+1)−est(s+1)2=−te−2t−e−2t\mathrm{Res}_{s=-2}\, F(s)e^{st} = \lim_{s\to -2}\dfrac{d}{ds}\left[\dfrac{e^{st}}{s+1}\right] = \lim_{s\to -2}\dfrac{t e^{st}(s+1) - e^{st}}{(s+1)^2} = -te^{-2t} - e^{-2t}。

第五步(组合出响应)。将留数相加得到 t≥0t \ge 0 时 y(t)=e−t−(1+t)e−2ty(t) = e^{-t} - (1+t)e^{-2t};两个指数项都在衰减且都不振荡,因此工程师确认该系统是过阻尼的。

若 ff 在正向简单闭曲线 γ\gamma 内部只有一个单极点 z0z_0,且 Res(f,z0)=3\mathrm{Res}(f,z_0)=3,那么 ∮γf(z) dz\oint_\gamma f(z)\,dz 等于多少?

f(z)=1z3f(z) = \dfrac{1}{z^3} 在 z=0z=0 处极点的阶数是多少?

Resz=0 ezz\mathrm{Res}_{z=0}\, \dfrac{e^z}{z} 等于多少?

一位电气工程师需要计算 ∫−∞∞dxx2+1\displaystyle\int_{-\infty}^{\infty} \dfrac{dx}{x^2+1} 来归一化一个洛伦兹频率响应。用 z=iz=i 处的留数,应该得到什么值?