外观
Chapter 5 留数定理
约 5831 字大约 19 分钟
2025-11-28
Part 1 留数定理
· 留数定理
定理 1 (留数定理)
设有界区域 G 的边界 C 为分段光滑的简单闭合曲线。若除有限个孤立奇点 bk, k=1,2,3,⋯,n 外,函数 f(z) 在 G 内单值解析,在 G 中连续,并且 f(z) 在边界 C 上连续,则沿区域 G 边界正向的积分
∮Cf(z)dz=2πik=1∑nresf(bk).
resf(bk) 称为 f(z) 在 bk 处的留数,它等于 f(z) 在孤立奇点 bk 的空心邻域内 Laurent 展开
f(z)=l=−∞∑∞al(k)(z−bk)l,0<∣z−bk∣<r
中 (z−bk)−1 的系数 a−1(k)。
/proof/
围绕每个孤立奇点 bk 作简单闭合曲线 γk,使 γk 均在 G 内,且互不交叠,则根据(多连通有界区域)Cauchy 定理及函数作 Laurent 展开的系数公式,就有
∮Cf(z)dz=k=1∑n∮γkf(z)dz=2πik=1∑na−1(k)=2πik=1∑nresf(bk).
所以,留数定理就是(多连通有界区域的)Cauchy 定理和 Laurent 展开系数公式的直接推论。它告诉我们,为了计算解析函数的围道积分值,只需计算出函数在孤立奇点处的留数即可。
求出 f(z) 在孤立奇点 b 处的留数,就是要求出 f(z) 在 z=b 的空心邻域内 Laurent 展开中 (z−b)−1 项的系数。在极点的情况下,通过微商的计算就可以比较容易地求出来。
设 z=b 是 f(z) 的 m 阶极点,则
f(z)=a−m(z−b)−m+a−m+1(z−b)−m+1+⋯+a−1(z−b)−1+a0+a1(z−b)+⋯,0<∣z−bk∣<r.
两端乘上 (z−b)m,有
(z−b)mf(z)=a−m+a−m+1(z−b)+⋯+a−1(z−b)m−1+a0(z−b)m+a1(z−b)m+1+⋯.
这时 a−1 是 (z−b)mf(z) 的展开式中 (z−b)m−1 项的系数,故
a−1=(m−1)!1dzm−1dm−1(z−b)mf(z)z=b.
特别是 z=b 为一阶极点的情形,
f(z)=a−1(z−b)−1+a0+a1(z−b)+a2(z−b)2+⋯,0<∣z−bk∣<r.
所以函数 f(z) 在该点的留数就是
a−1=z→blim(z−b)f(z).
可见,函数 f(z) 在复平面内一阶极点 z=b 处的留数一定不为零。
常见的情况是 f(z) 可以表示为 Q(z)P(z) 的形式,其中 P(z) 和 Q(z) 都在 b 点及其邻域内解析,且 P(b)=0,z=b 是 Q(z) 的一阶零点,即 Q(b)=0,Q′(b)=0,则
a−1=z→blim(z−b)f(z)=z→blim(z−b)Q(z)P(z)=Q′(b)P(b).
应用留数概念,可以方便地讨论有理函数的部分分式。
例如,要求将函数
f(z)=(z−1)(z−2)(z−3)1
部分分式,
(z−1)(z−2)(z−3)1=z−1A+z−2B+z−3C.
那么,三个待定常数,A, B 和 C,正好就是函数 f(z) 在一阶极点 z=1, z=2 和 z=3 点处的留数。因此
A=res(z−1)(z−2)(z−3)1z=1=21,
B=res(z−1)(z−2)(z−3)1z=2=−1,
C=res(z−1)(z−2)(z−3)1z=3=21.
如果函数 f(z) 具有高阶极点,也可以类似地处理。例如,
(z−1)2(z−2)(z−3)1=(z−1)2A+z−1B+z−2C+z−3D.
容易看出
A=res(z−1)2(z−2)(z−3)1z=1=21,
B=res(z−1)2(z−2)(z−3)1z=1=43,
C=res(z−1)2(z−2)(z−3)1z=2=−1,
D=res(z−1)2(z−2)(z−3)1z=3=41.
· 例题
/example/
求 z2+11 在复平面 C 内孤立奇点处的留数。
/proof/
z=±i 是它的一阶极点,
resf(±i)=2z1z=±i=±2i.
/example/
求 z2eiaz−eibz 在复平面 C 内孤立奇点处的留数。
/proof/
z=0 是它的一阶极点,
resf(0)=z→0limz⋅z2eiaz−eibz=z→0limzeiaz−eibz=i(a−b).
/example/
求 1/(z2+1)3 在复平面 C 内孤立奇点处的留数。
/proof/
z=±i 是它的三阶极点,
resf(±i)=2!1dz2d2[(z∓i)3⋅(z2+1)31]z=±i=±163i.
· 无穷远点的留数
以上的讨论均局限于复平面 C 内的孤立奇点。如果 ∞ 点不是 f(z) 的非孤立奇点,那么可以定义
resf(∞)=2πi1∮C′f(z)dz,
其中 C′ 为绕 ∞ 点正向(即顺时针方向)一周的简单封闭曲线,在围道内除 ∞ 点外的其他点均解析,∞ 点是 f(z) 的唯一可能的孤立奇点。需要注意,resf(∞) 并不是 f(z) 在 ∞ 点邻域内 Laurent 展开中 z1 项的系数。
这是因为,作变换 t=1/z,则
resf(∞)=2πi1∮C′f(z)dz=−2πi1∮Ct2f(1/t)dt
=−t2f(1/t) 在 t=0 \text{ 点邻域内幂级数展开中 } t−1 项的系数
=−f(1/t) 在 t=0 点邻域内幂级数展开中 t1 项的系数
=−f(z) 在 z=∞ 点邻域内幂级数展开中 z−1 项的系数.
在这个结果中,与有限远处不同之处在于:
(1) 从结果上说,函数 f(z) 在 ∞ 点的留数,等于 −f(z) 在 ∞ 点邻域内幂级数展开中 z−1 项的系数,这里多了一个负号。
(2) 从概念上说,由于 z−1 项属于 f(z) 在 ∞ 点邻域内幂级数展开式的正则部分,因此,即使 f(z) 在 ∞ 点解析,resf(∞) 也可能不为 0。反之,即使 ∞ 点是 f(z) 的孤立奇点,甚至是一阶极点,resf(∞) 也可能为 0。
Part 2 留数定理应用
· 有理三角函数积分
作为应用留数定理计算定积分的第一类例子,研究有理三角函数的积分
I=∫02πR(sinθ,cosθ)dθ,
其中 R 是 sinθ, cosθ 的有理函数,即 R 是由 sinθ, cosθ 的多项式相除得到的函数,R 在积分区间 [0,2π] 上是连续的。作变换 z=eiθ,则
sinθ=2izz2−1,cosθ=2zz2+1,dθ=izdz,
相应的积分路径则变为 z 复平面内的一条封闭围道——单位圆的圆周 ∣z∣=1。于是,
I=∮∣z∣=1R(2izz2−1,2zz2+1)izdz.
有理三角函数 R(sinθ,cosθ) 在积分区间 [0,2π] 上连续,保证了有理函数 R(2izz2−1,2zz2+1) 在积分围道(单位圆的圆周)上无奇点。
如果被积函数 R(2izz2−1,2zz2+1)iz1 在单位圆内部只有有限个孤立奇点,那么可以应用留数定理得
∮∣z∣=1R(2izz2−1,2zz2+1)izdz=2π∣z∣<1∑res{z1R(2izz2−1,2zz2+1)}.
· 无穷积分
第二类可以用留数定理计算的定积分是无穷积分
I=∫−∞∞f(x)dx.
在复平面上看,积分的路径是实轴,并不构成闭合围道。为了应用留数定理计算这种类型的积分,其基本原则是:
(1) 将实函数 f(x) 的定义域延拓到复平面成为复函数 f(z),
(2) 补上适当的积分路径从而形成复平面内的闭合围道。为此可以在上半平面补上以原点 O 为圆心、R 为半径的半圆弧 CR。由实轴上的线段 [−R,R] 及 CR 就构成闭合围道。如果函数 f(z) 在闭合围道上没有奇点,而且在围道内部只有有限个孤立奇点,则可以应用留数定理计算围道积分,而后令 R→∞,如果极限存在,即可算出积分
I=∫−∞∞f(x)dx.
/example/
计算积分 ∫−∞∞(1+x2)3dx。
/proof/
考虑复变积分 ∮C(1+z2)3dz,根据留数定理,有
∮C(1+z2)3dz=∫−RR(1+x2)3dx+∫CR(1+z2)31dz=2πi⋅res(1+z2)31z=i.
因为
z→∞lim[z⋅(1+z2)31]=0,
根据大圆弧引理即得
R→∞lim∫CR(1+z2)31dz=0.
取极限 R→∞,并代入
resf(±i)=2!1dz2d2[(z∓i)3⋅(z2+1)31]z=±i=±163i.
最后就求得
∫−∞∞(1+x2)3dx=83π.
不难看出,应用留数定理计算这种类型的无穷积分时,函数 f(z) 应满足下列条件:
(1) f(z) 在上半平面除有限个孤立奇点外处处解析,在实轴上没有奇点;
(2) 在 0≤argz≤π 范围内,当 ∣z∣→∞ 时,zf(z) 一致地趋于 0,即 ∀ε>0, ∃M(ε)>0,使当 ∣z∣≥M,并且 0≤argz≤π 时,∣zf(z)∣<ε。
上面关于留数定理的应用条件还可以放宽:
如果函数 f(z) 在上半平面内没有非孤立奇点,但是有无穷多个孤立奇点 bn, n=1,2,⋯,只要存在曲线序列 {Cm},每一条 Cm 都与实轴上从 −Rm 到 Rm 的直线段构成闭合围道,在围道上没有 f(z) 的奇点,
且当 m→∞ 时,Cm 上的点 z 的模 ∣z∣ 和 Rm 都趋于 ∞,同时 m→∞lim∫Cmf(z)dz=0,则
∫−∞∞f(x)dx=2πin=1∑∞resf(bn).
在上述基本原则下,围道的选取具有一定灵活性。
· 含三角函数无穷积分
第三类可以用留数定理计算的定积分是
I=∫−∞∞f(x)cospxdxorI=∫−∞∞f(x)sinpxdx.
这里不妨假设 p>0。
处理这种类型的积分,仍可采用无穷积分的半圆形围道。至于被积函数,通常并不取为 f(z)cospz 或 f(z)sinpz。这是因为 z=∞ 是 cospz, sinpz 的本性奇点,当 ∣z∣=R→∞ 时,函数 cospz, sinpz 的行为略显复杂,在计算
R→∞lim∫CRf(z)cospzdzorR→∞lim∫CRf(z)sinpzdz
时似乎会遇到一点困难。避开这一困难的方法是将被积函数取为 f(z)eipz。如果 f(z)eipz 在上半平面内只有有限个孤立奇点,则可以利用留数定理计算沿闭合围道的积分,有
∮Cf(z)eipzdz=∫−RRf(x)eipxdx+∫CRf(z)eipzdz=∫−RRf(x)(cospx+isinpx)dx+∫CRf(z)eipzdz.
这样,只要能够计算出 R→∞lim∫CRf(z)eipzdz,然后分别比较实部和虚部,就可以求得积分 ∫−∞∞f(x)cospxdx 和 ∫−∞∞f(x)sinpxdx。
为此,介绍一个引理。
引理 1 (Jordan 引理)
设在 0≤argz≤π 范围内,当 ∣z∣→∞ 时 Q(z) 一致地趋于 0,则
R→∞lim∫CRQ(z)eipzdz=0,
其中 p>0,CR 是以原点为圆心、R 为半径的上半圆弧。
/proof/
当 z 在 CR 上时,z=Reiθ,在 0≤argz≤π 范围内,当 ∣z∣→∞ 时 Q(z) 一致地趋于 0,意味着 ∀ε>0,存在与 argz 无关的 M(ε)>0,使得当 ∣z∣=R>M,且 0≤argz≤π 时,∣Q(z)∣<ε,
∫CRQ(z)eipzdz=∫0πQ(Reiθ)eipR(cosθ+isinθ)Reiθidθ≤∫0π∣Q(Reiθ)∣e−pRsinθRdθ<εR∫0πe−pRsinθdθ=2εR∫0π/2e−pRsinθdθ.
证明的关键在于精确估计 sinθ 值。当 0≤θ≤π/2 时有 sinθ≥2θ/π,所以
∫CRQ(z)eipzdz<2εR∫0π/2e−pR⋅2θ/πdθ=pεπ(1−e−pR).
这样,就证明了
R→∞lim∫CRQ(z)eipzdz=0.
上面介绍了应用留数定理计算含三角函数的无穷积分
∫−∞∞Q(x)sinpxdxand∫−∞∞Q(x)cospxdx,p>0
的一种做法,就是采用半圆形围道计算复变积分 ∮CQ(z)eipzdz。这种做法的优点是,只要 Q(z) 在上半平面范围内,当 ∣z∣→∞ 时一致地趋于 0,根据 Jordan 引理就能判断
R→∞lim∫CRQ(z)eipzdz=0.
我们可以列举出种种理由,说明在构造复变积分时,为什么不是简单地将实函数 Q(x)sinpx, Q(x)cospx 延拓为 Q(z)sinpz, Q(z)cospz,核心的理由是 z=∞ 是 sinpz, cospz 的本性奇点,或者说,sinpz, cospz 中含有 e−ipz,因而给处理沿 CR 的积分带来一些困难。
应该说,这种分析与讨论有助于我们理解如何选择复变积分(包括被积函数与积分围道两个方面),但绝不可以将上面提到的困难绝对化,更不应该引申出不正确的结论。
例如,因为 z=∞ 是 sinpz, cospz 的本性奇点,当 z 按不同方式逼近 ∞ 时,sinpz, cospz 可以逼近不同的值,或者说,z→∞ 时函数 sinpz, cospz 的极限均不存在,这些说法无疑都是正确的,但是,我们并不能由此就得出
R→∞lim∫CRQ(z)sinpzdzorR→∞lim∫CRQ(z)cospzdz
不存在”的臆断。
正是基于这一思想,我们现在就来探讨应用留数定理直接计算围道积分 ∮CQ(z)sinpzdz 与 ∮CQ(z)cospzdz 的可行性。应用这种新方法,在计算某些积分时可能更加简单。
为此目的,只需要建立一个新的引理,先称之为“补充引理”,它是留数定理与 Jordan 引理相结合的产物。
补充引理:
设函数 Q(z) 只有有限个奇点,且在下半平面的范围内,当 ∣z∣→∞ 时一致地趋近于 0,则
R→∞lim∫CRQ(z)e−ipzdz=2πi×Full plane∑res{Q(z)e−ipz}
=−2πi×res{Q(z)e−ipz}z=∞,
其中 p>0,CR 是以原点为圆心、R 为半径的半圆弧,位于上半平面
· 瑕积分
如果实变积分是瑕积分(例如瑕点为 x=a),则在处理相应的复变积分 ∮Cf(z)dz 时,实轴上的 z=a 点也是被积函数的奇点,那么我们选择积分路径时必须绕开奇点而构成闭合的积分围道。
/example/
计算主值积分 v.p.∫−∞∞x(1+x+x2)dx。
/proof?
这是一个反常积分,其反常性既表现在积分区间为无穷区间,又表现为 x=0 是被积函数的瑕点。此积分在主值意义下存在。因此,在应用留数定理计算此积分时,应考虑复变积分
∮Cz(1+z+z2)dz,
积分围道 C 由以原点为圆心、δ 为半径的小半圆弧 Cδ 和以原点为圆心、R 为半径的大半圆弧 CR 以及直线段 −R→−δ 和 δ→R 构成。于是,根据留数定理,有
∮Cz(1+z+z2)dz=∫−R−δx(1+x+x2)dx+∫Cδz(1+z+z2)dz+∫δRx(1+x+x2)dx+∫CRz(1+z+z2)dz=2πi⋅resz(1+z+z2)1z=e2πi/3=−3π−iπ.
由大圆弧引理,可以判断
R→∞lim∫CRz(1+z+z2)dz=0.
为了计算沿小圆弧 Cδ 积分的极限值,则可以应用小圆弧引理。因为
z→0lim[z⋅z(1+z+z2)1]=1,
所以有
δ→0lim∫Cδz(1+z+z2)dz=−πi.
取极限 R→∞, δ→0,就得到
v.p.∫−∞∞x(1+x+x2)dx=−3π.
在有些情况(例如,含三角函数的无穷积分)下,由于在计算复变积分中,并不是简单地将被积函数 f(x) 换成 f(z),因而本来实积分并不是瑕积分,但在复变函数的围道积分中,积分路径上却可能出现奇点。
如果考虑复变积分 ∮Czeizdz。虽然 x=0 不是原来实积分中被积函数 xsinx 的瑕点,但是 z=0 却是新的复变积分中被积函数 zeiz 的奇点。因此积分路径需要绕过奇点 z=0 。
∮Czeizdz=∫−R−δxeixdx+∫Cδzeizdz+∫δRxeixdx+∫CRzeizdz.
在积分围道包围的区域内,被积函数解析,故围道积分为 0。根据 Jordan 引理和小圆弧引理,又有
R→∞lim∫CRzeizdz=0,δ→0lim∫Cδzeizdz=−πi.
因此
∫−∞∞xeixdx=πi.
比较两端的实部和虚部,即得
v.p.∫−∞∞xcosxdx=0,∫−∞∞xsinxdx=π.
就复变积分而言,在积分路径上可以有奇点,但这种奇点一般说来只能是一阶极点。如果是二阶或二阶以上的极点,或是本性奇点,沿小圆弧 Cδ 的积分极限值就可能不存在。
· 多值函数复变积分
如果把定义域扩展到复平面,这些积分的被积函数是多值函数。我们运用围道积分计算这类积分时,就需要明确规定被积函数的函数值,例如,适当规定被积函数的单值分支。一种常见的类型就是
I=∫0∞xs−1Q(x)dx,
其中 s 为实数。从复平面上看,Q(z) 单值,在正实轴上没有奇点。被积函数中的 zs−1,当 s 不等于整数时,就是一个多值函数。原来积分中的积分变量应该理解为 argz=0。
为了计算这种类型的积分,通常考虑相应的复变积分为 ∮Czs−1Q(z)dz。
由于 z=0 及 ∞ 是 zs−1 的分支点,所以需要将复平面沿正实轴割开,并规定割线上岸 argz=0。
这时的积分路径由割开的大、小同心圆弧(圆心为 z=0 点,半径分别为 R 和 δ)及割线上、下岸组成。割线成为(玦形)积分围道的一部分。只是需要注意,在计算留数时,要遵守上面对于 zs−1 所作的限制,即 0≤argz≤2π。
· 计算无穷级数
设 f(z) 是已知函数,在 C 内除有限个极点外解析,而且 f(z) 的极点都不是整数。如果存在另一个函数 G(z),在 C∖{0,±1,±2,⋯} 内解析,而 z=0,±1,±2,⋯ 是 G(z) 的一阶极点,且在这些极点处的留数均为 1。于是,作闭合围道 CN,将 n=0,±1,±2,⋯,±N 包围在内,根据留数定理,就有
∮CNG(z)f(z)dz=2πi⎩⎨⎧n=−N∑Nf(n)+f(z)’s pole∑res[G(z)f(z)]⎭⎬⎫.
如果当 N→∞ 时,能求出 ∮CNG(z)f(z)dz 的极限值(例如,在一定条件下为零),就可以算出无穷级数 n=−∞∑∞f(n) 的和。
这里存在两个问题。一是要找到这样的函数 G(z),二是如何求出 N→∞lim∮CNG(z)f(z)dz 的极限值。前一个问题的答案是 G(z) 可取为 πcotπz。后一个问题,则用引理 2 解决。
引理 2
设 f(z) 在 C 内除了有限个孤立奇点外处处解析,若存在常数 R>0 和 M>0,使当 ∣z∣>R 时,∣f(z)∣≤M,则
N→∞lim∮CNπcotπzf(z)dz=0,
其中积分围道 CN 为正方形,四个顶点位于 (N+1/2)(1±i) 和 −(N+1/2)(1±i)。
根据这个引理,立即可以证明下面的定理。
定理 2
若函数 f(z) 在 C 上除有限个非整数的极点外处处解析,且存在常数 R>0 和 M>0,使当 ∣z∣>R 时,∣zf(z)∣≤M,则
n=−∞∑∞f(n)=N→∞limn=−N∑Nf(n)=−f(z)’s pole∑res{πcotπzf(z)}.
更新日志
2026/3/21 08:55
查看所有更新日志
0abe4-complex5于3fbb7-complex于67c9f-ts于