分析学
留数与留数定理
通过对围道内部极点处的留数求和来计算复变周线积分的方法。
直观从孤立极点到无穷和
想象在复平面上沿一条闭合回路 γ 行走,绕过某个在若干孤立点处发散的函数 f。在每个这样的点附近,f 的表现就像螺旋楼梯:越靠近,数值旋转得越快、增长得越猛。留数定理揭示了一个惊人的事实:整个回路积分的值只取决于内部每个被围住的点处提取出的一个数——它的留数——而与回路的具体形状无关,只要回路不穿过任何奇点。
带孤立极点函数的区域着色图;颜色围绕极点旋转,这种旋转编码了该处的留数。大学孤立奇点与洛朗级数
定义: 孤立奇点及其留数
若 f 在 z0 周围某个去心圆盘上全纯,但在 z0 处本身不全纯,则称 z0 为 f 的孤立奇点。在该去心圆盘上,f 有唯一的洛朗展开 f(z)=∑n=−∞∞an(z−z0)n,而项 (z−z0)−1 的系数 a−1 称为 f 在 z0 处的留数,记作 Res(f,z0)=a−1。奇点按出现的负幂个数分类:没有负幂是可去奇点(极限 limz→z0f(z) 存在),有限个负幂、最高阶为 m 是极点,无穷多个负幂(如在 z=0 处的 sin(1/z))是本性奇点。
f(z)=n=−∞∑∞an(z−z0)n 实际操作中,没有人会为了求某个极点处的留数而展开整个洛朗级数。对于单极点(阶 m=1),有 Res(f,z0)=limz→z0(z−z0)f(z)。对于更高阶的极点 m,公式推广为 Res(f,z0)=(m−1)!1limz→z0dzm−1dm−1[(z−z0)mf(z)],该式在代入之前先通过求导消去奇异部分。
Res(f,z0)=(m−1)!1z→z0limdzm−1dm−1[(z−z0)mf(z)] 孤立奇点的分类及其留数计算| 奇点类型 | 洛朗级数 | 留数公式 |
|---|
| 可去奇点 | 无负幂项 | Res(f,z0)=0 |
| 单极点(m=1) | 一个负幂项 | Res(f,z0)=limz→z0(z−z0)f(z) |
| m 阶极点 | 有限个负幂项 | Res(f,z0)=(m−1)!1limz→z0dzm−1dm−1[(z−z0)mf(z)] |
| 本性奇点 | 无穷多个负幂项 | 直接读取 a−1 |
大学关键定理
设 f 在单连通区域上除有限个孤立奇点 z1,…,zn 外全纯,这些奇点都在一条正向简单闭曲线 γ 内部。则 ∮γf(z)dz=2πi∑kRes(f,zk)。
为什么成立?
它把一个困难的几何问题(沿曲线积分)化归为一个简单的代数问题(把有限个数加起来),因为将围道在每个极点附近变形后会收缩成一个很小的圆,在那里洛朗级数完成了全部工作。
证明
第一步(变形围道)。根据柯西积分定理,对任何围住 f 全纯区域的闭曲线都有 ∮γf(z)dz=0。由于 f 仅在 z1,…,zn 处不全纯,在每个 zk 周围画一个半径为 ε 的正向小圆 Ck,半径足够小使这些圆彼此不相交且都在 γ 内部。从 γ 到每个 Ck 切开细缝,把 γ 与各 Ck 之间的区域变成一个 f 全纯的单连通区域,因此该区域边界上的积分为 0;细缝两侧的贡献成对抵消,剩下 ∮γf(z)dz=∑k∮Ckf(z)dz。
第二步(计算每个小圆)。固定 k,在包含 Ck 的去心圆盘上,将 f 在 zk 附近展开为洛朗级数 f(z)=∑n=−∞∞an(z−z0)n。对满足 n=−1 的每一项 an(z−zk)n,它在该去心圆盘上都有单值原函数,故沿闭合小圆 Ck 积分为 0;只有项 a−1(z−zk)−1 保留下来。
第三步(计算剩余积分)。将 Ck 参数化为 z=zk+εeiθ(θ∈[0,2π]),于是 dz=iεeiθdθ 且 z−zkdz=idθ。因此 ∮Ckz−zka−1dz=a−1∫02πidθ=2πia−1=2πiRes(f,zk)。
第四步(求和)。将第三步代入第一步的等式,得到 ∮γf(z)dz=2πi∑kRes(f,zk),这正是留数定理。
若 f 在 z0 处有 m 阶极点,即 g(z)=(z−z0)mf(z) 可全纯延拓到 z0 且 g(z0)=0,则 Res(f,z0)=(m−1)!1limz→z0dzm−1dm−1[(z−z0)mf(z)]。特别地,对单极点(m=1),有 Res(f,z0)=limz→z0(z−z0)f(z)。
为什么成立?
乘以 (z−z0)m 会消去洛朗级数的奇异部分,使其变成系数恰为导数的普通泰勒级数——于是提取 a−1 就变成了一次普通的微积分计算。
证明
第一步(明确写出洛朗级数)。由于 f 在 z0 处有 m 阶极点,其洛朗级数形如 f(z)=(z−z0)ma−m+⋯+z−z0a−1+a0+a1(z−z0)+⋯,其中 a−m=0。
第二步(消去极点)。两边乘以 (z−z0)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 在 z0 处全纯。
第三步(把系数识别为导数)。在此幂级数中,a−1 是 (z−z0)m−1 的系数。对任意全纯函数 g(z)=∑jcj(z−z0)j,泰勒公式给出 cj=j!g(j)(z0),因此 a−1=(m−1)!g(m−1)(z0)。
第四步(结论)。由于 Res(f,z0)=a−1 且 g(m−1)(z0)=limz→z0dzm−1dm−1[(z−z0)mf(z)],代入即得 Res(f,z0)=(m−1)!1limz→z0dzm−1dm−1[(z−z0)mf(z)]。令 m=1 化简导数与阶乘,即得 Res(f,z0)=limz→z0(z−z0)f(z)。
大学实际应用与典型例题
留数把原本难以处理的实定积分变成了简单的算术运算,这正是它们广泛出现在物理学和工程学中的原因:在信号处理中计算振荡的傅里叶和拉普拉斯变换、在量子场论中计算散射振幅和传播子、通过拉普拉斯逆变换求解控制系统的微分方程,以及在解析数论和组合学中估计计数函数的渐近行为。
例题: 没有初等原函数捷径的实积分
用留数定理计算 ∫−∞∞x2+1dx。(该积分也可用代换 x=tanθ 求得,但留数方法能推广到代换法无法处理的积分。)
解答
第一步(建立围道)。考虑 f(z)=z2+11 以及由实轴上的线段 [−R,R] 与逆时针方向的上半圆 ∣z∣=R 组成的闭合围道 γR(R 取较大值)。
第二步(估计圆弧部分)。在半圆弧上,∣f(z)∣≤R2−11,而弧长为 πR,所以弧上的贡献至多为 R2−1πR→0(当 R→∞ 时)。
第三步(找出围道内的极点)。f(z)=z2+11 在 z=i 和 z=−i 处有单极点;只有 z=i 位于上半圆内部。由单极点公式,Res(f,i)=limz→i(z−i)(z−i)(z+i)1=2i1。
第四步(应用留数定理并取极限)。对任意 R>1 都有 ∮γRf(z)dz=2πi⋅2i1=π,这等于实积分加上消失的圆弧项,令 R→∞ 即得 ∫−∞∞x2+1dx=π。
例题: 控制理论中的拉普拉斯逆变换
某控制工程师用传递函数 F(s)=(s+1)(s+2)21 对一个有阻尼的系统建模,需要求出其脉冲响应 y(t)(即拉普拉斯逆变换),以检验系统是否能不振荡地趋于稳定。
解答
第一步(回顾反演积分)。拉普拉斯逆变换为 y(t)=2πi1∫c−i∞c+i∞F(s)estds,对于在无穷远处衰减的有理函数 F(s),该积分等于 F(s)=(s+1)(s+2)21 全部极点处 F(s)est 的留数之和乘以 2πi。
第二步(确定极点位置)。F(s)=(s+1)(s+2)21 在 s=−1 处有一个单极点,在 s=−2 处有一个 2 阶极点。
第三步(单极点处的留数)。Ress=−1F(s)est=lims→−1(s+1)(s+1)(s+2)2est=(−1)2e−t=e−t。
第四步(二阶极点处的留数)。Ress=−2F(s)est=lims→−2dsd[s+1est]=lims→−2(s+1)2test(s+1)−est=−te−2t−e−2t。
第五步(组合出响应)。将留数相加得到 t≥0 时 y(t)=e−t−(1+t)e−2t;两个指数项都在衰减且都不振荡,因此工程师确认该系统是过阻尼的。
若 f 在正向简单闭曲线 γ 内部只有一个单极点 z0,且 Res(f,z0)=3,那么 ∮γf(z)dz 等于多少?
f(z)=z31 在 z=0 处极点的阶数是多少?
Resz=0zez 等于多少?
一位电气工程师需要计算 ∫−∞∞x2+1dx 来归一化一个洛伦兹频率响应。用 z=i 处的留数,应该得到什么值?