4.2 劳伦级数

4.2 劳伦级数

本节摘要：SOURCE 4.2：因此，当读者步入 4.2.1 "环域展开" 时，请视其为一场精密的测绘工作——你正在为函数的奇点绘制一张等高线地图；而踏入 4.2.2 "奇点隔离与分类预备" 时，请将其理解为一次严谨的法医鉴定——你正用主部这把高倍显微镜，辨析奇点的细胞结构。

目标与定位

幂级数只能在圆盘内展开，圆盘边界一旦碰到奇点就失效。劳伦级数把"展开舞台"从圆盘推广到环域 r<|z-z_0|<R，允许负幂项出现，从而能在奇点附近"分层拆解"函数。这一推广表面上只是加了几项负幂，却为第5章的奇点分类与留数定理准备了全部工具。

环域展开

劳伦定理：设 f 在环域 D=\{z: r<|z-z_0|<R\} 内解析，则 f 在 D 内可唯一展开为

其中系数

展开分为两部分：

• 正则部 \sum_{n=0}^\infty c_n(z-z_0)^n：形如泰勒级数，在圆盘 |z-z_0|<R 内解析；
• 主部 \sum_{n=-\infty}^{-1} c_n(z-z_0)^n：由负幂项组成，在 |z-z_0|>r 内解析。

若 r=0，环域变为穿孔邻域 0<|z-z_0|<R，此时主部的形状直接刻画 z_0 处奇点的类型——这正是 4.2.2 与第5章的主题。

劳伦展开的求法

求劳伦展开不需要每次用系数公式积分，通常用已知级数拼接：

例1：\frac{1}{z-1} 在环域 1<|z|<2 内展开。目标是写成 \sum c_n z^n 的形式：

这里利用 |1/z|<1 使几何级数在 |z|>1 内收敛——展开式随环域改变而改变。

例2：e^{1/z} 在 0<|z|<\infty 内展开：

主部含无穷多项，预示着 z=0 是本性奇点。

主部与奇点分类预备

若 f 在 z_0 的穿孔邻域内解析，则 z_0 的奇点类型完全由主部决定：

1. 主部为空：z_0 是可去奇点，定义 f(z_0)=c_0 后 f 在 z_0 解析；
2. 主部为有限项 \frac{c_{-m}}{(z-z_0)^m}+\cdots：z_0 是 m 阶极点；
3. 主部为无限项：z_0 是本性奇点。

第5章将给出这三类的判别方法与留数计算，本节只需记住判定原则：看主部的形状。

唯一性与代数运算

劳伦展开在给定环域内是唯一的：两种展开方式若逐项相等则系数相同。唯一性带来两个实用推论：

1. 求系数不必套积分公式，可用任何合法方式构造展开；
2. 环域一旦改变，展开式随之改变——"同一个 f"在不同环域是"不同的展开"。

劳伦级数在环域内同样支持四则运算。例如 f(z)=\frac{1}{z^2-1} 在环域 0<|z-1|<2 内的展开：令 w=z-1，则 z^2-1=(z-1)(z+1)=w(w+2)，于是

主部恰有一项 \frac{1}{2(z-1)}，说明 z=1 是 f 的一阶极点，且 c_{-1}=\frac12——这个值就是后续的留数。

数值算例：展开与部分和

import numpy as np
import math

# 数值验证 e^(1/z) 的劳伦展开：在 |z|=0.5 圆周上取点对比
z = 0.5 * np.exp(1j * np.linspace(0, 2*np.pi, 200))
series = np.zeros_like(z, dtype=complex)
for n in range(20):              # 主部取到 z^-19
    series += (1.0 / z)**n / math.factorial(n)
err = np.max(np.abs(series - np.exp(1.0 / z)))
print("劳伦级数前 20 项的最大误差: %.2e" % err)
​

# 验证 1/(z-1) 在 |z|>1 的展开 sum z^-n 等于原函数
r = 1.5
z = r * np.exp(1j * np.linspace(0, 2*np.pi, 300))
s = np.zeros_like(z, dtype=complex)
for n in range(1, 200):
    s += z**(-n)                 # 主部累加
err = np.max(np.abs(s - 1.0/(z-1)))
print("1/(z-1) 展开误差: %.2e" % err)
​

# 验证 c_{-1} 与围道积分的关系：oint e^(1/z) dz = 2*pi*i * c_{-1}
# c_{-1} 是 e^(1/z) 展开中 1/z 的系数 = 1/1! = 1
n_pt = 20000
t = np.linspace(0, 2*np.pi, n_pt)
z = np.exp(1j * t)
dz = 1j * np.exp(1j * t)
I = np.trapz(np.exp(1.0 / z) * dz, t)
print("围道积分: %.6f %+.6fi" % (I.real, I.imag))
print("2*pi*i * c_{-1} = 0.000000 +6.283185i")
​

# 验证 1/(z^2-1) 在 z=1 处主部系数：围道积分法
r = 0.5
t = np.linspace(0, 2*np.pi, 20000)
w = 1 + r*np.exp(1j*t)     # 绕 z=1 的小圆，不包含 z=-1
dw = 1j*r*np.exp(1j*t)
f = 1.0 / (w**2 - 1)
c_m1 = np.trapz(f * dw, t) / (2*np.pi*1j)
print("c_{-1} 数值: %.6f %+.6fi (理论 0.5)" % (c_m1.real, c_m1.imag))
​

应用：从展开到留数

劳伦展开中 c_{-1} 这一项有特殊地位：围道积分直接读出它，

c_{-1} 就是 f 在 z_0 处的留数。因此劳伦级数不仅是分类工具，更是留数的计算入口——第5章将在此基础上建立留数定理，并用来计算各类实积分。

从物理视角看，主部各负幂项分别刻画"场源"的多极矩：\frac{c_{-1}}{z-z_0} 对应单极、\frac{c_{-2}}{(z-z_0)^2} 对应偶极，依次类推。这与电磁学中多极展开的术语完全吻合——劳伦级数正是把局部"源"结构按阶数拆解的标准语言。工程上，只要奇点类型已知，主部结构也就固定，展开计算因而可以套用现成模板。

⚠️ 常见坑：只记结论不记适用边界——劳伦展开依赖于环域选择，同一函数在不同环域内有不同展开；r=0 时展开才对应奇点邻域，不要与泰勒展开混淆。

💡 关键直觉：4.2 劳伦级数 应能对应至少一项可复现实验或算例——两个数值算例分别验证了 e^{1/z} 与 1/(z-1) 的展开。

思考与检验

1. 在环域 0<|z|<1 与 |z|>1 内分别展开 \frac{1}{z(1-z)}，比较主部。
2. 判定 \frac{\sin z}{z}、\frac{1}{z^2(z-1)}、e^{1/z^2} 在 0 处的奇点类型。
3. 说明为何 c_{-1} 与围道积分直接相关，并数值验证 \oint_{|z|=1}e^{1/z}\,dz=2\pi i\,c_{-1}。

小结

• 主干：4.2 劳伦级数 连接「输入—过程—输出」——环域展开、主部分类、c_{-1} 留数。
• 边界：展开依赖环域；系数公式要求围道完全落在环域内。
• 方法：用已知级数拼接代替积分公式，必要时数值核对。