跳到主要内容

留数理论、辐角原理与积分应用

留数定理是复分析计算的核心工具，它将复杂的路径积分转化为孤立奇点处的代数运算。

一、留数 (Residue) 的深度导出​

1. 孤立奇点与 Laurent 展开​

若 z0z_0 是 f(z)f(z) 的孤立奇点，则在圆环域 0<∣z−z0∣<R0 < |z-z_0| < R 内有： f(z)=∑n=−∞∞cn(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} c_n (z-z_0)^n 其中 c−1c_{-1} 被定义为 留数 (Residue)，记作 Res(f,z0)\text{Res}(f, z_0)。

2. 留数计算的通用准则​

  • mm 阶极点： 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}} [ (z-z_0)^m f(z) ]
  • 一阶商极点：若 f(z)=P(z)Q(z)f(z) = \frac{P(z)}{Q(z)} 且 Q(z0)=0,Q′(z0)≠0Q(z_0)=0, Q'(z_0) \neq 0，则 Res(f,z0)=P(z0)Q′(z0)\text{Res}(f, z_0) = \frac{P(z_0)}{Q'(z_0)}。

二、留数定理与无穷远点​

1. 留数基本定理​

∮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)

2. 无穷远点留数 (Residue at Infinity)​

定义 Res(f,∞)=Res(−1z2f(1z),0)\text{Res}(f, \infty) = \text{Res}\left( -\frac{1}{z^2} f(\frac{1}{z}), 0 \right)。 全平面留数总和定律：若 f(z)f(z) 在扩充复平面内只有有限个奇点，则： ∑k=1nRes(f,zk)+Res(f,∞)=0\sum_{k=1}^n \text{Res}(f, z_k) + \text{Res}(f, \infty) = 0

三、实积分的复处理：Jordan 引理与主值​

对于含有 sin⁡(ax)\sin(ax) 或 cos⁡(ax)\cos(ax) 的积分，利用 Jordan 引理： 若 lim⁡R→∞f(Reiθ)=0\lim_{R \to \infty} f(Re^{i\theta}) = 0，则： lim⁡R→∞∫CRf(z)eiazdz=0(a>0)\lim_{R \to \infty} \int_{C_R} f(z) e^{iaz} dz = 0 \quad (a > 0)

柯西主值 (Cauchy Principal Value)​

当积分路径穿过奇点时，利用半圆路径收缩定义主值： P.V.∫−∞∞f(x)dx=lim⁡ϵ→0,R→∞[∫−Rx0−ϵ+∫x0+ϵR]f(x)dx\text{P.V.} \int_{-\infty}^{\infty} f(x) dx = \lim_{\epsilon \to 0, R \to \infty} \left[ \int_{-R}^{x_0-\epsilon} + \int_{x_0+\epsilon}^{R} \right] f(x) dx

✍️ 深度练习与 C++ 模拟​

练习 1：利用 Rouché 定理判定根分布​

确定 z5+15z+1=0z^5 + 15z + 1 = 0 在圆环 1<∣z∣<21 < |z| < 2 内的根个数。

Check Solution

解析：

  1. ∣z∣<2|z| < 2 内：令 f(z)=z5,g(z)=15z+1f(z) = z^5, g(z) = 15z + 1。在 ∣z∣=2|z|=2 上，∣f(z)∣=32,∣g(z)∣≤31|f(z)| = 32, |g(z)| \le 31。故 ∣f∣>∣g∣|f| > |g|，原方程有 5 个根。
  2. ∣z∣<1|z| < 1 内：令 f(z)=15z,g(z)=z5+1f(z) = 15z, g(z) = z^5 + 1。在 ∣z∣=1|z|=1 上，∣f(z)∣=15,∣g(z)∣≤2|f(z)| = 15, |g(z)| \le 2。故 ∣f∣>∣g∣|f| > |g|，原方程有 1 个根。
  3. 结论：在 1<∣z∣<21 < |z| < 2 内有 5−1=45 - 1 = 4 个根。 □\square

练习 2：C++ 数值模拟复路径积分​

编写 C++ 程序，利用中点法则或辛普森法则模拟沿单位圆周对 f(z)=1/zf(z) = 1/z 的路径积分，验证留数定理。

Check Solution
#include <iostream>
#include <complex>
#include <vector>
#include <cmath>

/**
 * @brief Complex Path Integrator
 * 验证 \oint_{|z|=1} (1/z) dz = 2\pi i
 */
int main() {
    using namespace std;
    typedef complex<double> cd;

    int N = 10000; // 离散步数
    double dt = 2.0 * M_PI / N;
    cd integral(0, 0);

    for (int i = 0; i < N; ++i) {
        double t = i * dt + dt/2.0; // 中点
        cd z(cos(t), sin(t));       // 路径 z(t) = e^{it}
        cd dz(-sin(t), cos(t));     // 微分 dz = i*e^{it} dt
        dz *= dt;
        
        integral += (1.0 / z) * dz;
    }

    cout << "Numerical Result: " << integral << endl;
    cout << "Expected (2*pi*i): (0, " << 2.0 * M_PI << ")" << endl;
    cout << "Absolute Error:   " << abs(integral - cd(0, 2.0*M_PI)) << endl;

    return 0;
}

模拟分析： 该模拟通过数值离散化闭曲线积分，验证了留数定理的宏观表现。即使在计算精度限制下，虚部依然高度趋近于 2π2\pi（约 6.28318），而实部由于抵消效应趋于 10−1610^{-16} 级别的零。这为理解“留数是路径积分的唯一贡献源”提供了实证。

练习 3：符号推导——无穷远点留数计算​

计算 f(z)=z2+1z2−1f(z) = \frac{z^2+1}{z^2-1} 在无穷远点的留数。

Check Solution

步骤：

  1. 令 w=1/zw = 1/z，计算 g(w)=−1w2f(1w)=−1w2(1/w)2+1(1/w)2−1=−1w21+w21−w2g(w) = -\frac{1}{w^2} f(\frac{1}{w}) = -\frac{1}{w^2} \frac{(1/w)^2+1}{(1/w)^2-1} = -\frac{1}{w^2} \frac{1+w^2}{1-w^2}。
  2. Res(f,∞)=Res(g,0)\text{Res}(f, \infty) = \text{Res}(g, 0)。
  3. g(w)=−1w2(1+w2)(1+w2+w4+… )=−1w2−2−2w2−…g(w) = -\frac{1}{w^2} (1+w^2)(1+w^2+w^4+\dots) = -\frac{1}{w^2} - 2 - 2w^2 - \dots
  4. 观察 1/w1/w 的系数，为 0。
  5. 结论：Res(f,∞)=0\text{Res}(f, \infty) = 0。 □\square

🚀 综合应用：Dirichlet 积分的证明​

利用上半平面大半圆路径 C=[−R,R]∪ΓRC = [-R, R] \cup \Gamma_R，我们可以严密证明： ∫0∞sin⁡xxdx=π2\int_0^\infty \frac{\sin x}{x} dx = \frac{\pi}{2} 这是复分析在解析数论与信号处理中不可或缺的经典应用。