绕原点一周,1/z 的积分是 2πi,1/z2 的积分却是零。后一个函数的奇点更强,为什么积分反而没有留下东西?上一章的逐项积分已经给了答案:闭路积分只挑出 Laurent 展开的负一次幂。这一章就沿着这条线,把局部系数变成计算积分和统计零点的方法。
留数只记录一个系数
设 a 是孤立奇点,在足够小的穿孔圆盘内写成
f(z)=n=−∞∑∞an(z−a)n.
留数定义为 Res(f,a)=a−1。沿正向圆周 Cρ:,代入 可得
∮Cρ(z−a)
Laurent 级数在这条圆周上一致收敛,允许逐项积分,因此
∮Cρf(z)dz=2πiRes(f
留数不是奇点的阶数,也不衡量函数有多大。1/z2 的留数是零,但原点仍是二阶极点;e1/z 的留数是 1,原点却是本性奇点。判断类型看整个主部,计算这个积分只取其中一项。

绕小圆积分只保留负一次幂的系数;其他负幂仍能造成极点,却不贡献这个积分。
让一条大曲线分成几条小圆
设 C 是分段光滑、逆时针的简单闭曲线,f 在曲线及其内部的一个邻域内,除有限个孤立奇点 a1,…,as 外全纯,并且这些奇点都不在 C 上。给出
∮Cf(z)dz=2πij=1∑
可以把它看成上一节小圆公式的汇总。给每个奇点挖去一个很小的圆盘,圆盘两两不交且都在大曲线内部;剩下的区域中 f 全纯。带孔区域的正向边界由逆时针的大曲线和顺时针的小圆构成。用第4章的 Cauchy 定理,边界积分为零,移项后,大曲线的积分等于各条逆时针小圆积分之和。
如果还想把区域切开来检查,可以从各个孔向外边界作切线。同一条切线在边界中出现两次,方向相反,积分抵消。最终保留下来的恰好是大边界和小圆,没有额外贡献。

挖去奇点以后,带孔区域的边界积分为零;把顺时针小圆移到另一侧,就得到留数的总贡献。
曲线不必是圆,曲线内的奇点也不必都是极点。真正需要的是上面的全纯范围与边界条件。把方向反过来,积分整体变号;曲线经过奇点时,这个普通闭路积分公式不能直接使用。
先判断局部结构,再选算法
分母只有一个简单零点
若 f=g/h,g,h 在 a 附近全纯,且 h(a)=0, h,则
h(z)=(z−a)q(z),q(a)=h′(a)
于是 (z−a)f(z)=g(z)/q(z) 在 a 附近全纯,负一次幂的系数就是它在 a 的值:
Res(f,a)=h′(a)g(a).
这里有个容易漏掉的区别:只有当 g(a)=0 时,a 才确实是简单极点。若 g(a)=0,分子的零点把分母抵消,奇点可去,公式仍然给出正确的留数 0。所以“分母有简单零点”不能单独推出“函数有简单极点”。
例题: 计算正向圆 ∣z∣=3 上
∮∣z∣=3z(z−2)4z−1d
圆内只有 0,2 两个简单极点,圆上没有奇点。用极限形式计算:
Res(f,0)=z−24z
两者相加为 4,所以积分为 8πi。若把圆半径改为 1,只剩原点的贡献,积分变成 πi。函数没变,圈住哪些奇点变了。
高阶极点,求导和展开都可以
若 a 是 m 阶极点,令 ϕ(z)=(z−a)mf(z)。它在 附近全纯且 。把 展成 Taylor 级数,再除以 ,要得到负一次幂,就应取 的 次项:
Res(f,a)=(m−1)!ϕ(m−1)(a)

简单极点取极限,高阶极点提取全纯因子的相应 Taylor 系数。分母幂次须在抵消后判断。
求导次数高时,我更愿意先看看现成的级数是否够用。例如
z3e2z=z
直接给出留数 2。用高阶公式也得到 (e2z)′′∣z=0/2!=。两种方法选择哪一种,取决于哪边能少做无用计算。
再看 1/[z2(z−1)] 在原点的留数。导数法中 ϕ(z)=1/(z−1),故 。展开法则给出
−z211−z1
负号不能省。若连同 z=1 一起圈住,另一个留数是 1,总积分反而为零。
抵消后,分母幂次会骗人
对 (ez−1−z)/z4,分母写着四次幂,但分子从 z2/2 开始。相除得到
z4ez−1−z=
原点是二阶极点,留数是 1/6。真正的阶数由抵消后的最负幂决定。若直接认定四阶,既会误报类型,也会把计算弄得更繁琐。

分子从二次项开始,抵消以后最负幂是 z−2,而留数来自另外的 z−1/6。
在下面的实验中改变分子的零点阶数和分母幂次,预测它会落在哪一种情形。观察主部哪些项消失,再改变圆的半径和方向,比较数值积分与 2πi 乘留数。留数为零时,特别检查页面显示的是“可去”还是“仍有极点”。
无穷远的负号从哪里来
在无穷远附近用 w=1/z 作坐标,微分同时变成 dz=−w−2dw。因此定义
Res(f,∞)=−Res(w2f(1/w
若外部 Laurent 展开中 f(z) 的 z−1 系数是 b−1,代入后 的 系数也为 ,所以无穷远留数是 。别只把函数里的 换掉,遗漏了微分。
当有限平面内只有有限个孤立奇点、且它们之外函数全纯时,取足够大的逆时针圆 ∣z∣=R。一方面积分等于有限点留数和的 2πi 倍;另一方面由外部 Laurent 展开,它等于 2πib−1。因此
a∈C∑Res(f,a)+Res(f,∞)=0.
这里无穷远也有一个完整穿孔邻域可用,不能把公式套给在有限处有无穷多个奇点、一直延伸到无穷远的情形。
前面的有理函数满足
z(z−2)4z−1=z
所以无穷远留数为 −4,恰与有限点的 1/2+7/2=4 相抵。这个方法特别适合只要“全部有限留数之和”、不需要各点数值的题。

有限点留数和为 4,无穷远留数为 −4。负号来自坐标替换时微分的变化。
用边界上的变化,统计内部的零点
现在让被积函数变成 f′/f。设 f 在简单闭曲线 C 及其内部的邻域内亚纯,曲线上没有零点或极点,且 C 逆时针。亚纯意味着这里只有全纯点和孤立极点,不能夹进本性奇点。
若 a 是 m 阶零点,局部写成 f(z)=(z−a)mg(z),其中 。直接求导再相除:
f(z)f′(z)=
第二项在 a 附近全纯,所以 f′/f 在这里的留数为 m。在 m 阶极点处,同样的分解用负幂,留数就变成 −m。套用刚证明的留数定理得到辐角原理:
2πi1∮Cf(z)
其中 N,P 分别是内部零点和极点的总阶数。一个二阶零点按两个计,一个三阶极点贡献负三。
这个整数为什么叫“辐角”的变化?令 w=f(z),积分变成沿像曲线 f(C) 的 ∫dw/w。沿曲线连续追踪辐角,而不是每次强行跳回主值区间;一周后 log∣w∣ 回到起点,辐角的净增量是 的整数倍。这个整数就是像曲线绕原点的绕数。因此看的是 ,并非只看原曲线 绕零几周。
例如 f(z)=(z−21)2/(z+,在 内有二阶零点,没有极点,故 ;在 内多出一个简单极点,变为 。半径恰为 或 时,边界经过零点或极点,不能继续沿用计数式。

两幅图只示意包含关系,不按比例。半径由 1 增至 3,计数从二阶零点的 2 变为 2−1。
在实验中调大圆的半径,分别观察原平面的零点、极点与像曲线。改变方向前预测符号,再看连续累积的辐角,而不要只读终点的主辐角。边界触点应使计数暂停;这不是程序少算了一个点,而是定理的条件已经失效。
不求根,也能确定个数
很多多项式没法方便地分解。此时可以在边界比较一个容易计数的主项与剩余部分。设 g,h 在 C 及内部的邻域全纯,并在 C 上满足
∣h(z)∣<∣g(z)∣.
Rouché 定理断言:g+h 与 g 在内部的零点总阶数相同。
这里的严格不等式有明确作用。把 g 连续变成 g+th,0≤t≤1,边界上始终有
∣g+th∣≥∣g∣−t∣h∣≥∣g∣−∣h∣>0.
所以整个变化过程中没有零点碰上边界。更精确地,辐角原理给出的
n(t)=2πi1∮C
一直是整数;分母在紧集 [0,1]×C 上远离零,被积函数随 t 一致连续,因而积分也连续。区间上的连续整数值函数只能是常数,这就证明两端的零点数相同。
例题: p(z)=z5+5z+2 在单位圆内和半径 2 的圆内各有几个零点?
在 ∣z∣=1 上选 g=5z,因为 ∣z5+2∣≤,内部恰有一个零点。在 上换成 ,有 ,内部恰有五个零点。两个圆周都无零点,所以环域 内有四个。主项不是固定选最高次项,而是根据边界上的大小来选。
练习
练习 1|留数相同,类型相同吗。 分别判断 z2/sinz、cosz/z2、e1/ 在原点的奇点类型和留数。
sinz=z(1−z2/6+⋯),所以第一项等于 z/(1−,原点可去,留数为零。第二项为 ,原点是二阶极点,留数仍为零。第三项有无限多个负偶次幂,是本性奇点,没有负一次幂,留数也为零。三者类型不同,留数却相同;不能由留数为零推断全纯。
练习 2|方向与抵消。 沿顺时针圆 ∣z∣=2 计算 ∮(ez−1)/z3dz。
展开为 z−2+21z−1+,原点是二阶极点、留数 。圆内没有其他奇点,顺时针再乘负号,答案为 。分子的一阶零点同时影响极点阶数,不能把分母三次幂直接当三阶。
练习 3|只求总和。 求 f(z)=(3z2+2)/(z3−8) 的全部有限点留数之和,不逐个写出三个极点。
当 ∣z∣ 足够大时,f(z)=3/z+O(z−3),所以无穷远留数为 −3。有限点留数之和是 。由此也知道任一正向大圆圈住所有极点后,积分为 。
练习 4|零点落在哪一圈。 证明 z4+4z+1 在 ∣z∣<1 内恰有一个零点,在 1<∣z 内恰有三个,均按重数计。
单位圆上 ∣z4+1∣≤2<4=∣4z∣,与 4z 比较得一个。半径 2 的圆上 ,与 比较得四个。严格不等式同时排除两条边界上的零点,相减得环域内三个。两次选择主项的依据都是边界上的估计。
1正向简单闭曲线内有一个二阶零点和一个三阶极点,边界没有零点或极点;辐角原理中 N−P 是多少?
2g/h 的分母 h 在 a 有简单零点,就能保证 a 是 g/h 的简单极点。