Skip to content

留数 ​

留数理论(Residue Theory)是复分析(Complex Analysis)中的核心内容之一。它将复变函数的解析性质与积分计算紧密结合,是孤立奇点理论的自然延伸。利用留数定理,我们可以将复杂的复变函数沿闭曲线的积分转化为函数在曲线内部孤立奇点处的“留数”之和,从而极大地简化了积分的计算。此外,留数理论也是计算某些难以用初等方法求解的实定积分的强有力工具。

—

留数的定义与留数定理 ​

留数的定义 ​

设 z0z_0 是函数 f(z)f(z) 的一个孤立奇点,即 f(z)f(z) 在去心邻域 D:0<∣z−z0∣<RD: 0 < |z - z_0| < R 内解析。根据罗朗(Laurent)级数理论,f(z)f(z) 在该去心邻域内可以展开为:

f(z)=∑n=−∞+∞cn(z−z0)nf(z) = \sum_{n=-\infty}^{+\infty} c_n (z - z_0)^n

其中展开式的系数 cnc_n 由下式决定:

cn=12πi∮Cf(z)(z−z0)n+1dz(n=0,±1,±2,⋯ )c_n = \frac{1}{2\pi i} \oint_{C} \frac{f(z)}{(z - z_0)^{n+1}} dz \quad (n = 0, \pm 1, \pm 2, \cdots)

这里的 CC 是 DD 内绕 z0z_0 的任意一条正向简单闭曲线。

特别地,当 n=−1n = -1 时,我们有:

c−1=12πi∮Cf(z)dzc_{-1} = \frac{1}{2\pi i} \oint_{C} f(z) dz

这个系数 c−1c_{-1} 包含了函数在该奇点处的某种“残留信息”,因此我们将其称为 f(z)f(z) 在孤立奇点 z0z_0 处的留数(或残数),记作 Res(f,z0)\text{Res}(f, z_0) 或 Res[f(z),z0]\text{Res}[f(z), z_0]。

定义 设 z0z_0 为 f(z)f(z) 的孤立奇点,f(z)f(z) 在其去心邻域内的罗朗展开式中 (z−z0)−1(z-z_0)^{-1} 的系数 c−1c_{-1} 称为 f(z)f(z) 在 z0z_0 处的留数,即:

Res(f,z0)=c−1=12πi∮Cf(z)dz\text{Res}(f, z_0) = c_{-1} = \frac{1}{2\pi i} \oint_{C} f(z) dz

留数基本定理 ​

留数定义直接给出了单奇点包裹情形下的积分结果。当区域内包含多个孤立奇点时,通过复连通区域的柯西积分定理,可以立即导出著名的留数定理。

定理 设 DD 是一个有界区域,其边界为正向周线或闭曲线族 CC。若函数 f(z)f(z) 在 DD 内除去有限个孤立奇点 z1,z2,⋯ ,znz_1, z_2, \cdots, z_n 外处处解析,且在 Dˉ=D∪C\bar{D} = D \cup C 上除这些奇点外连续,则 f(z)f(z) 沿 CC 的正向积分等于 f(z)f(z) 在 DD 内所有奇点处的留数之和的 2πi2\pi i 倍:

∮Cf(z)dz=2πi∑k=1nRes(f,zk)\oint_{C} f(z) dz = 2\pi i \sum_{k=1}^{n} \text{Res}(f, z_k)

—

留数的计算方法 ​

计算留数的核心在于寻找罗朗级数中 c−1c_{-1} 的值。虽然可以直接展开函数,但针对不同类型的孤立奇点(主要是极点),有着更为简便的代数计算公式。

可去奇点与本性奇点处的留数 ​

  1. 可去奇点: 若 z0z_0 是 f(z)f(z) 的可去奇点,则 f(z)f(z) 在 z0z_0 的去心邻域内没有负幂项,即 c−1=0c_{-1} = 0。因此,可去奇点处的留数必为 0。

  2. 本性奇点: 若 z0z_0 是 f(z)f(z) 的本性奇点,通常没有通用的代数公式,必须将 f(z)f(z) 直接展开为罗朗级数,然后直接观察 (z−z0)−1(z-z_0)^{-1} 的系数。

极点处的留数计算公式 ​

若 z0z_0 是 f(z)f(z) 的极点,我们可以利用导数和极限来直接计算留数。

1. mm 阶极点的通用公式 ​

定理 若 z0z_0 是 f(z)f(z) 的 mm 阶极点,则 f(z)f(z) 在 z0z_0 处的留数为:

Res(f,z0)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\text{Res}(f, z_0) = \frac{1}{(m-1)!} \lim_{z \to z_0} \frac{d^{m-1}}{dz^{m-1}} \left[ (z - z_0)^m f(z) \right]

2. 1 阶极点(单极点)的特例公式 ​

当 m=1m = 1 时,上述公式简化为:

Res(f,z0)=lim⁡z→z0(z−z0)f(z)\text{Res}(f, z_0) = \lim_{z \to z_0} (z - z_0) f(z)

此外,在实际应用中,单极点经常以分式形式 f(z)=P(z)Q(z)f(z) = \frac{P(z)}{Q(z)} 出现。

定理 设 f(z)=P(z)Q(z)f(z) = \frac{P(z)}{Q(z)},其中 P(z)P(z) 和 Q(z)Q(z) 在 z0z_0 处解析。若 P(z0)≠0P(z_0) \neq 0,Q(z0)=0Q(z_0) = 0 且 Q′(z0)≠0Q'(z_0) \neq 0(即 z0z_0 是 Q(z)Q(z) 的一重零点,从而为 f(z)f(z) 的单极点),则:

Res(f,z0)=P(z0)Q′(z0)\text{Res}(f, z_0) = \frac{P(z_0)}{Q'(z_0)}

这个公式在处理三角函数或多项式分式时效率极高。

无穷远点处的留数 ​

如果 f(z)f(z) 在区域 R<∣z∣<+∞R < |z| < +\infty 内解析,则称 ∞\infty 轴为 f(z)f(z) 的孤立奇点。 定义 f(z)f(z) 在 ∞\infty 处的留数为:

Res(f,∞)=12πi∮C−f(z)dz=−12πi∮C+f(z)dz\text{Res}(f, \infty) = \frac{1}{2\pi i} \oint_{C^-} f(z) dz = - \frac{1}{2\pi i} \oint_{C^+} f(z) dz

其中 CC 是一条顺时针方向(对无穷远点而言是正向)绕原点的闭曲线。

通过作变换 ζ=1z\zeta = \frac{1}{z},可以得到无穷远点留数的计算公式:

Res(f,∞)=Res[−1ζ2f(1ζ),0]\text{Res}(f, \infty) = \text{Res}\left[ -\frac{1}{\zeta^2} f\left(\frac{1}{\zeta}\right), 0 \right]

定理 若 f(z)f(z) 在整个扩充复平面上只有有限个孤立奇点 z1,z2,⋯ ,znz_1, z_2, \cdots, z_n 以及 ∞\infty,则它们处的留数总和为零:

∑k=1nRes(f,zk)+Res(f,∞)=0\sum_{k=1}^{n} \text{Res}(f, z_k) + \text{Res}(f, \infty) = 0

—

留数的典型应用 ​

留数定理最重要的应用之一是计算实变函数中难以用常规微积分方法解决的某些实定积分。通常的方法是将实积分路线看作复平面上某条闭曲线的一部分,然后构造辅助复变函数进行闭合路径积分。

1. 计算 ∫02πR(cos⁡θ,sin⁡θ)dθ\int_{0}^{2\pi} R(\cos\theta, \sin\theta) d\theta 型积分 ​

这类积分的特点是区间为 [0,2π][0, 2\pi],被积函数是关于 sin⁡θ\sin\theta 和 cos⁡θ\cos\theta 的有理函数。

转化策略: 令 z=eiθ=cos⁡θ+isin⁡θz = e^{i\theta} = \cos\theta + i\sin\theta,当 θ\theta 从 00 变到 2π2\pi 时,zz 沿着复平面上的单位圆周 C:∣z∣=1C: |z|=1 逆时针运动。 利用欧拉公式,我们有:

dθ=dziz,cos⁡θ=z+z−12=z2+12z,sin⁡θ=z−z−12i=z2−12izd\theta = \frac{dz}{iz}, \quad \cos\theta = \frac{z + z^{-1}}{2} = \frac{z^2 + 1}{2z}, \quad \sin\theta = \frac{z - z^{-1}}{2i} = \frac{z^2 - 1}{2iz}

将其代入原积分,即可转化为沿单位圆周 CC 的复变函数有理分式积分:

∫02πR(cos⁡θ,sin⁡θ)dθ=∮∣z∣=1R(z2+12z,z2−12iz)dziz\int_{0}^{2\pi} R(\cos\theta, \sin\theta) d\theta = \oint_{|z|=1} R\left( \frac{z^2+1}{2z}, \frac{z^2-1}{2iz} \right) \frac{dz}{iz}

最后只需计算该函数在单位圆 ∣z∣<1|z| < 1 内部的所有孤立奇点的留数即可。

2. 计算 ∫−∞+∞P(x)/Q(x)dx\int_{-\infty}^{+\infty} P(x)/Q(x) dx 型广义积分 ​

其中 P(x)P(x) 和 Q(x)Q(x) 为实多项式,且满足:Q(x)Q(x) 在实轴上无零点,且 deg(Q)≥deg(P)+2\text{deg}(Q) \ge \text{deg}(P) + 2(保证积分收敛)。

转化策略: 在复平面上构造闭合回路 Γ=[−R,R]∪CR\Gamma = [-R, R] \cup C_R,其中 CRC_R 是上半平面上以原点为圆心、 RR 为半径的半圆周。 根据留数定理,当 RR 足够大以至于包含了上半平面内的所有奇点 zkz_k(Im(zk)>0\text{Im}(z_k) > 0)时:

∫−RRP(x)Q(x)dx+∫CRP(z)Q(z)dz=2πi∑Im(zk)>0Res(PQ,zk)\int_{-R}^{R} \frac{P(x)}{Q(x)} dx + \int_{C_R} \frac{P(z)}{Q(z)} dz = 2\pi i \sum_{\text{Im}(z_k)>0} \text{Res}\left( \frac{P}{Q}, z_k \right)

由大弧段引理(Jordan引理的简化版),当 R→+∞R \to +\infty 时,沿半圆周 CRC_R 的积分趋于 0。因此:

∫−∞+∞P(x)Q(x)dx=2πi∑Im(zk)>0Res(P(z)Q(z),zk)\int_{-\infty}^{+\infty} \frac{P(x)}{Q(x)} dx = 2\pi i \sum_{\text{Im}(z_k)>0} \text{Res}\left( \frac{P(z)}{Q(z)}, z_k \right)

3. 计算 ∫−∞+∞P(x)Q(x)eiαxdx\int_{-\infty}^{+\infty} \frac{P(x)}{Q(x)} e^{i\alpha x} dx 型积分(约当引理) ​

这类积分常用于傅里叶变换的求解,其中 α>0\alpha > 0,deg(Q)≥deg(P)+1\text{deg}(Q) \ge \text{deg}(P) + 1。

在同样的半圆闭回路 Γ\Gamma 下,我们需要借助**约当引理(Jordan’s Lemma)**来证明大半圆弧上的积分贡献在 R→+∞R \to +\infty 时为 0:

引理 设 CRC_R 是上半平面上的半圆周 ∣z∣=R,Im(z)>0|z|=R, \text{Im}(z) > 0。若在 CRC_R 上 lim⁡R→+∞F(z)=0\lim_{R \to +\infty} F(z) = 0 均匀成立,则对于任何 α>0\alpha > 0,有:

lim⁡R→+∞∫CRF(z)eiαzdz=0\lim_{R \to +\infty} \int_{C_R} F(z) e^{i\alpha z} dz = 0

基于此引理,这类积分的计算公式为:

∫−∞+∞P(x)Q(x)eiαxdx=2πi∑Im(zk)>0Res(P(z)Q(z)eiαz,zk)\int_{-\infty}^{+\infty} \frac{P(x)}{Q(x)} e^{i\alpha x} dx = 2\pi i \sum_{\text{Im}(z_k)>0} \text{Res}\left( \frac{P(z)}{Q(z)} e^{i\alpha z}, z_k \right)

通过取其实部或虚部,可以顺带求出形如 ∫−∞+∞P(x)Q(x)cos⁡(αx)dx\int_{-\infty}^{+\infty} \frac{P(x)}{Q(x)}\cos(\alpha x) dx 或 ∫−∞+∞P(x)Q(x)sin⁡(αx)dx\int_{-\infty}^{+\infty} \frac{P(x)}{Q(x)}\sin(\alpha x) dx 的积分。

—

典型例题解析 ​

例 计算积分 I=∫02π12+cos⁡θdθI = \int_{0}^{2\pi} \frac{1}{2 + \cos\theta} d\theta。

证明 令 z=eiθz = e^{i\theta},则 cos⁡θ=z2+12z\cos\theta = \frac{z^2+1}{2z},dθ=dzizd\theta = \frac{dz}{iz}。代入原积分得:

I=∮∣z∣=112+z2+12zdziz=∮∣z∣=124z+z2+1dzi=2i∮∣z∣=11z2+4z+1dzI = \oint_{|z|=1} \frac{1}{2 + \frac{z^2+1}{2z}} \frac{dz}{iz} = \oint_{|z|=1} \frac{2}{4z + z^2 + 1} \frac{dz}{i} = \frac{2}{i} \oint_{|z|=1} \frac{1}{z^2 + 4z + 1} dz

设 f(z)=1z2+4z+1f(z) = \frac{1}{z^2 + 4z + 1},其分母的零点(即 f(z)f(z) 的单极点)为:

z1=−2+3,z2=−2−3z_1 = -2 + \sqrt{3}, \quad z_2 = -2 - \sqrt{3}

显然,∣z1∣=2−3<1|z_1| = 2 - \sqrt{3} < 1(在单位圆内),而 ∣z2∣=2+3>1|z_2| = 2 + \sqrt{3} > 1(在单位圆外)。

利用分式型单极点留数公式,计算 z1z_1 处的留数:

Res(f,z1)=1(z2+4z+1)′∣z=z1=12z1+4=12(−2+3)+4=123\text{Res}(f, z_1) = \left. \frac{1}{(z^2 + 4z + 1)'} \right|_{z=z_1} = \frac{1}{2z_1 + 4} = \frac{1}{2(-2+\sqrt{3}) + 4} = \frac{1}{2\sqrt{3}}

由留数定理:

I=2i⋅2πi⋅Res(f,z1)=4π⋅123=2π3I = \frac{2}{i} \cdot 2\pi i \cdot \text{Res}(f, z_1) = 4\pi \cdot \frac{1}{2\sqrt{3}} = \frac{2\pi}{\sqrt{3}}

用 VitePress 搭建 · 托管在 GitHub Pages