留数定理与实积分计算
本节侧重留数定理与积分应用。第一次学习时,先读求留数:方法选择与计算,掌握孤立奇点分类、实际极点阶数与六种计算方法,再进入专项习题和 44 题综合题集。
物理问题引入:如何用复数「绕道」解决算不出来的实积分?
在理论物理中,从经典电动力学辐射场计算、阻尼系统受迫共振响应,到量子散射截面与量子场论单圈传播子积分,物理学家无时无刻不在面对形如下式的实数反常定积分:
𝐼=∫+∞−∞cos(𝑘𝑥)𝑥2+𝑎2d𝑥,∫∞0sin𝑥𝑥d𝑥,∫∞0sin2𝑥𝑥2d𝑥I = int_-inf^+inf cos(k x)x^2 + a^2 dx, quad int_0^inf sin xx dx, quad int_0^inf sin^2 xx^2dx若使用大一微积分的初等换元或分部积分,根本无法找到初等原函数进行牛顿 - 莱布尼茨公式计算. 然而,复分析提供了一种「降维打击」式的数学策略:
- 复化延拓:将实数轴上的被积函数延拓到复平面 𝑧 =𝑥 +𝑖𝑦z = x + iy;
- 闭合围道:给实轴积分补上一段半圆弧、扇形或矩形路径,使其成为封闭围道;
- 留数定理:围道内的积分完全不依赖于整条路径的繁琐细节,仅仅由被包围奇点处的局部代数系数(留数)之和乘以 2𝜋𝑖2pi i 直接决定!
实数领域中几乎不可逾越的积分天堑,在复平面上只需做简单的代数除法与求导即可迎刃而解.
1. 留数的诞生哲学与计算工具箱
1.1 留数(Residue) 诞生的代数直觉
设 𝑧0z_0 为函数 𝑓(𝑧)f(z) 的孤立奇点.在以 𝑧0z_0 为中心的去心圆环域内将函数作洛朗级数展开:
𝑓(𝑧)=⋯+𝑏2(𝑧−𝑧0)2+𝑏1𝑧−𝑧0+𝑎0+𝑎1(𝑧−𝑧0)+𝑎2(𝑧−𝑧0)2+…f(z) = dots + b_2(z - z_0)^2 + b_1z - z_0 + a_0 + a_1(z - z_0) + a_2(z - z_0)^2 + dots现在,我们用一个包含 𝑧0z_0 的逆时针简单闭合曲线 𝐶C 对整个级数进行逐项积分.回顾基础整幂积分公式:
∮𝐶(𝑧−𝑧0)𝑛d𝑧={2𝜋𝑖,𝑛=−10,𝑛≠−1oint_C (z - z_0)^n dz = cases 2pi i, & n = -1 0, & n != -1 cases- 所有正整数次幂项 (𝑧 −𝑧0)𝑛(z-z_0)^n(𝑛 ≥0n >= 0)沿闭合回路的积分全部为零(柯西积分定理);
- 所有负二次及更高负幂次项 (𝑧 −𝑧0)−𝑘(z-z_0)^-k(𝑘 ≥2k >= 2)沿闭回路的积分也全部为零(由高阶导数推广公式);
- 全级数无穷多项积分之后,唯独只有负一次幂项 𝑏1𝑧−𝑧0b_1/z - z_0 留下了非零的数值贡献 2𝜋𝑖𝑏12pi i b_1!
因为 𝑏1b_1(习惯上也记作 𝑎−1a_-1 或 𝑐−1c_-1)是围道积分中唯一「遗留」下来的非零核心,物理学家与数学家称之为 留数(Residue),记为:
Res[𝑓(𝑧),𝑧0]=𝑏1=𝑎−1Res[f(z), z_0] = b_1 = a_-11.2 极点留数的极速计算法则
在物理实战中,无须每次都求出完整的洛朗级数,而是利用如下求导与极限法则快速求得留数:
法则 1:一阶单极点极限法与「商式法则」
若 𝑧0z_0 为一阶单极点,直接乘以 (𝑧 −𝑧0)(z - z_0) 消除分母后取极限:
Res[𝑓,𝑧0]=lim𝑧→𝑧0(𝑧−𝑧0)𝑓(𝑧)Res[f, z_0] = lim_z -> z_0 (z - z_0) f(z)商式法则(物理中应用最广的留数公式)
若被积函数为分式 𝑓(𝑧) =𝑃(𝑧)𝑄(𝑧)f(z) = P(z)Q(z),其中 𝑃(𝑧0) ≠0P(z_0) != 0,𝑧0z_0 为分母 𝑄(𝑧)Q(z) 的一阶单零点(即 𝑄(𝑧0) =0Q(z_0) = 0 且 𝑄′(𝑧0) ≠0Q'(z_0) != 0),利用洛必达法则:
Res[𝑓,𝑧0]=lim𝑧→𝑧0(𝑧−𝑧0)𝑃(𝑧)𝑄(𝑧)=lim𝑧→𝑧0𝑃(𝑧)⋅lim𝑧→𝑧0𝑧−𝑧0𝑄(𝑧)−𝑄(𝑧0)=𝑃(𝑧0)𝑄′(𝑧0)Res[f, z_0] = lim_z -> z_0 (z - z_0)P(z)Q(z) = lim_z -> z_0 P(z) * lim_z -> z_0 z - z_0Q(z) - Q(z_0) = P(z_0)Q'(z_0) 只需将极点代入分子,除以分母的导数,即可在数秒内算出留数!
法则 2:𝑚m 阶高阶极点法则
若 𝑧0z_0 为 𝑚m 阶极点(𝑚 ≥2m >= 2),两端乘 (𝑧 −𝑧0)𝑚(z-z_0)^m 并求 𝑚 −1m-1 阶导数:
Res[𝑓,𝑧0]=1(𝑚−1)!lim𝑧→𝑧0d𝑚−1d𝑧𝑚−1[(𝑧−𝑧0)𝑚𝑓(𝑧)]Res[f, z_0] = 1(m - 1)! lim_z -> z_0 d^m-1dz^m-1left[(z - z_0)^m f(z)right]法则 3:无穷远点留数与全平面留数和定理
定义扩充复平面中无穷远点 𝑧 =∞z = inf 的留数(注意围道方向关于无穷远点为 顺时针,负号来源于此外侧积分取向):
Res[𝑓(𝑧),∞]=−12𝜋𝑖∮𝐶𝑅𝑓(𝑧)d𝑧=−Res[1𝑧2𝑓(1𝑧),0]Res[f(z), inf] = -12pi ioint_C_R f(z)dz = -Resleft[1z^2fleft(1zright), 0right]全平面留数和定理:若函数 𝑓(𝑧)f(z) 在扩充复平面 ¯ℂC 上除有限个孤立奇点 𝑧1,…,𝑧𝑛z_1, dots, z_n 外处处解析,则 其在全平面所有有限奇点的留数与无穷远点留数之和严格为零:
𝑛∑𝑘=1Res[𝑓(𝑧),𝑧𝑘]+Res[𝑓(𝑧),∞]=0sum_k=1^n Res[f(z), z_k] + Res[f(z), inf] = 0当围道内部包含大量繁琐奇点,而围道外部奇点极少(或只有无穷远点)时,利用此定理从外部反求可极大简化运算.
2. 留数定理(Residue Theorem)
定理(柯西留数定理):设简单闭曲线 𝐶C 为正向闭围道,𝑓(𝑧)f(z) 在 𝐶C 上全纯解析,在 𝐶C 内部除有限个孤立奇点 𝑧1,𝑧2,…,𝑧𝑛z_1, z_2, dots, z_n 外处处解析,则:
∮𝐶𝑓(𝑧)d𝑧=2𝜋𝑖𝑛∑𝑘=1Res[𝑓(𝑧),𝑧𝑘]oint_C f(z)dz = 2pi i sum_k=1^n Res[f(z), z_k]例题:Fibonacci 生成函数在不同半径圆周上的围道积分
计算积分 𝐼𝑅 =∮|𝑧|=𝑅𝑧1−𝑧−𝑧2d𝑧I_R = oint_|z|=R z1 - z - z^2dz.
分母零点为两个单极点:𝑧1 =√5−12 ≈0.618z_1 = 5-12 ~ 0.618,𝑧2 = −√5+12 ≈ −1.618z_2 = -5+12 ~ -1.618.
由商式法则计算两奇点的留数(𝑃(𝑧) =𝑧,𝑄(𝑧) =1 −𝑧 −𝑧2,𝑄′(𝑧) = −1 −2𝑧P(z) = z, Q(z) = 1 - z - z^2, Q'(z) = -1 - 2z):
Res[𝐹,𝑧1]=𝑧−1−2𝑧∣𝑧1=√5−12−1−(√5−1)=−√5−12√5=√5−510Res[F, z_1] = left.z-1 - 2zright|_z_1 = 5-12-1 - (5-1) = -5-125 = 5 - 510 Res[𝐹,𝑧2]=𝑧−1−2𝑧∣𝑧2=−√5+12−1+(√5+1)=−√5+12√5=−5+√510Res[F, z_2] = left.z-1 - 2zright|_z_2 = -5+12-1 + (5+1) = -5+125 = -5 + 510 当 𝑅 =0.5R = 0.5 时:圆周内无任何奇点,由柯西定理:𝐼0.5 =0I_0.5 = 0;
当 𝑅 =1.0R = 1.0 时:圆周内仅包含奇点 𝑧1z_1:
𝐼1.0=2𝜋𝑖Res[𝐹,𝑧1]=2𝜋𝑖(√5−510)=𝜋𝑖√5−55I_1.0 = 2pi i Res[F, z_1] = 2pi i left(5 - 510right) = pi i 5 - 55当 𝑅 =2.0R = 2.0 时:圆周同时包含 𝑧1z_1 与 𝑧2z_2:
𝐼2.0=2𝜋𝑖(Res[𝐹,𝑧1]+Res[𝐹,𝑧2])=2𝜋𝑖(√5−510−5+√510)=−2𝜋𝑖I_2.0 = 2pi i (Res[F, z_1] + Res[F, z_2]) = 2pi i left(5 - 510 - 5 + 510right) = -2pi i
可用无穷远展开 𝐹(𝑧) = −1/𝑧 +𝑂(𝑧−2)F(z)=-1/z+O(z^-2) 核对有限留数之和为 −1-1。讲义中的半径 0.6,0.7,1.6,1.70.6,0.7,1.6,1.7 分别包围零个、一个、一个、两个极点,对应计算见综合题集。圆周经过极点时,普通留数定理不适用。
3. 理论物理实积分的五大围道计算范式
根据被积函数的代数与解析特征,求解实数定积分有五大经典围道范式:
flowchart TD
A["实积分计算需求"] --> B["三角周期积分 \int_0^{2\pi} R(\cos\theta,\sin\theta)\mathrm{d}\theta<br/>【范式一:单位圆代换】"]
A --> C["实轴无穷有理分式 \int_{-\infty}^\infty \frac{P(x)}{Q(x)}\mathrm{d}x<br/>【范式二:上半平面大半圆弧】"]
A --> D["傅里叶振荡型积分 \int_{-\infty}^\infty R(x)e^{ikx}\mathrm{d}x<br/>【范式三:若当引理半圆围道】"]
A --> E["实轴含一阶/二阶奇点(Dirichlet型)<br/>【范式四:凹陷围道法 Indented Contours】"]
A --> F["高斯傅里叶变换型积分<br/>【范式五:矩形闭合围道】"]
范式一:单位圆代换求解三角函数周期积分
适用于形如 ∫2𝜋0𝑅(cos𝜃,sin𝜃)d𝜃int_0^2pi R(costheta, sintheta)dtheta 的有理三角积分.
令 𝑧 =𝑒𝑖𝜃z = e^itheta(𝜃theta 从 00 变到 2𝜋2pi 正向绕单位圆周 𝐶 :|𝑧| =1C: |z| = 1 一周):
cos𝜃=𝑧+𝑧−12,sin𝜃=𝑧−𝑧−12𝑖,d𝜃=d𝑧𝑖𝑧costheta = z + z^-12, quad sintheta = z - z^-12i, quad dtheta = dziz将三角定积分完全化为单位圆周上的复有理分式围道积分,通过留数定理求出落在 |𝑧| <1|z| < 1 内部的极点留数之和.
例题:计算 ∫2𝜋0d𝜃𝑎+cos𝜃int_0^2pi da + costheta(𝑎 >1a > 1)
代入参数化变换:
𝐼=∮|𝑧|=11𝑎+𝑧+𝑧−12d𝑧𝑖𝑧=∮|𝑧|=12𝑖(𝑧2+2𝑎𝑧+1)d𝑧I = oint_|z|=1 1a + z + z^-12 dziz = oint_|z|=1 2i(z^2 + 2az + 1)dz 分母零点为 𝑧± = −𝑎 ±√𝑎2−1z_+/- = -a +/- a^2 - 1.由于 𝑎 >1a > 1:
- |𝑧−| =𝑎 +√𝑎2−1 >1|z_-| = a + a^2 - 1 > 1(位于圆外);
- |𝑧+| =𝑎 −√𝑎2−1 <1|z_+| = a - a^2 - 1 < 1(唯一位于单位圆内的单极点).
由商式法则计算 𝑧+z_+ 处的留数:
Res=2/𝑖2𝑧+2𝑎∣𝑧=𝑧+=2/𝑖2√𝑎2−1=1𝑖√𝑎2−1Res = left.2/i2z + 2aright|_z = z_+ = 2/i2a^2 - 1 = 1ia^2 - 1 代入留数定理:
𝐼=2𝜋𝑖⋅1𝑖√𝑎2−1=2𝜋√𝑎2−1I = 2pi i * 1ia^2 - 1 = 2a^2 - 1
范式二:无穷实轴有理分式积分
适用于 ∫+∞−∞𝑃(𝑥)𝑄(𝑥)d𝑥int_-inf^+inf P(x)/Q(x)dx,其中分母在实轴上无实根,且分母次数比分子至少高 2 次(deg𝑄 ≥deg𝑃 +2deg Q >= deg P + 2).
补上半平面大半圆弧 𝐶𝑅 :𝑧 =𝑅𝑒𝑖𝜃C_R: z = R e^itheta(𝜃 ∈[0,𝜋]theta in [0, pi]).当 𝑅 →∞R -> inf 时大圆弧积分由于 ∼1/𝑅~ 1/R 衰减为 0.
∫+∞−∞𝑃(𝑥)𝑄(𝑥)d𝑥=2𝜋𝑖∑Im𝑧𝑘>0Res[𝑃(𝑧)𝑄(𝑧),𝑧𝑘]int_-inf^+inf P(x)Q(x)dx = 2pi i sum_Im z_k > 0 Resleft[P(z)Q(z), z_kright]
范式三:傅里叶振荡积分与约当引理(Jordan's Lemma)
适用于求解频域谱响应反演积分:
𝐼=∫+∞−∞𝑅(𝑥)𝑒𝑖𝑘𝑥d𝑥(𝑘>0)I = int_-inf^+inf R(x) e^ikx ,dx quad (k > 0)引理(约当引理/Jordan's Lemma):
设 𝐶𝑅C_R 为上半平面大半圆弧 𝑧 =𝑅𝑒𝑖𝜃z = R e^itheta(𝜃 ∈[0,𝜋]theta in [0, pi]).若在此半圆弧上函数 𝑅(𝑧)R(z) 一致地满足当 𝑅 →∞R -> inf 时 𝑀𝑅 =max𝑧∈𝐶𝑅|𝑅(𝑧)| →0M_R = max_z in C_R|R(z)| -> 0(仅需一次衰减,不要求二次衰减!),且 𝑘 >0k > 0,则:
lim𝑅→∞∫𝐶𝑅𝑅(𝑧)𝑒𝑖𝑘𝑧d𝑧=0lim_R -> inf int_C_R R(z) e^ikzdz = 0补充说明
证明要点:在 𝐶𝑅C_R 上 |𝑒𝑖𝑘𝑧| =|𝑒𝑖𝑘𝑅(cos𝜃+𝑖sin𝜃)| =𝑒−𝑘𝑅sin𝜃|e^ikz| = |e^ik R(costheta + isintheta)| = e^-k Rsintheta.
∣∫𝐶𝑅𝑅(𝑧)𝑒𝑖𝑘𝑧d𝑧∣≤𝑀𝑅𝑅∫𝜋0𝑒−𝑘𝑅sin𝜃d𝜃=2𝑀𝑅𝑅∫𝜋/20𝑒−𝑘𝑅sin𝜃d𝜃left|int_C_R R(z)e^ikzdzright| <= M_R R int_0^pi e^-k Rsinthetadtheta = 2 M_R R int_0^pi/2 e^-k Rsinthetadtheta 利用约当不等式 sin𝜃 ≥2𝜋𝜃sintheta >= 2/pitheta:
≤2𝑀𝑅𝑅∫𝜋/20𝑒−𝑘𝑅2𝜃𝜋d𝜃=2𝑀𝑅𝑅⋅𝜋2𝑘𝑅(1−𝑒−𝑘𝑅)=𝜋𝑘𝑀𝑅(1−𝑒−𝑘𝑅)𝑀𝑅→0←←←←←←←←←←→0<= 2 M_R R int_0^pi/2 e^-k R2pidtheta = 2 M_R R * 2kRleft(1 - e^-kRright) = k M_R (1 - e^-kR) M_R -> 0 0例题:求解狄拉克滤波实积分 ∫∞−∞cos𝑘𝑥𝑥2+𝑎2d𝑥int_-inf^inf cos kx/x^2 + a^2dx(𝑘 >0,𝑎 >0k > 0, a > 0)
构造辅助积分 ∮𝐶𝑒𝑖𝑘𝑧𝑧2+𝑎2d𝑧oint_C e^ikzz^2 + a^2dz.围道由实轴 [ −𝑅,𝑅][-R, R] 和上半平面半圆弧 𝐶𝑅C_R 构成.
取其实部即得:
∫+∞−∞cos𝑘𝑥𝑥2+𝑎2d𝑥=𝜋𝑎𝑒−𝑘𝑎int_-inf^+inf cos kxx^2 + a^2dx = ae^-ka
范式四:实轴奇点与凹陷围道法(Indented Contours)
当被积函数的奇点恰好落在实数积分路径上(例如原点 𝑥 =0x = 0)时,我们必须构造 绕过该实数奇点的微小圆弧(顺时针或逆时针小半圆凹陷).
小圆弧引理(分数留数引理):
若 𝑧0z_0 为一阶单极点,圆弧 𝐶𝜀C_varepsilon 的半径为 𝜀varepsilon,圆心角为 𝛼alpha(顺时针方向),则当 𝜀 →0varepsilon -> 0 时:
lim𝜀→0∫𝐶𝜀𝑓(𝑧)d𝑧=−𝛼𝑖Res[𝑓(𝑧),𝑧0]lim_varepsilon -> 0 int_C_varepsilon f(z)dz = -alpha i Res[f(z), z_0]若为顺时针小半圆弧(𝛼 =𝜋alpha = pi),其积分贡献严格等于 −𝜋𝑖Res[𝑓,𝑧0]-pi i Res[f, z_0].
例题 1:Dirichlet 积分 ∫∞0sin𝑥𝑥d𝑥 =𝜋2int_0^inf sin x/xdx = pi/2
辅助复函数 𝐹(𝑧) =𝑒𝑖𝑧𝑧F(z) = e^izz,在原点具有一阶极点.
设计上半平面闭合凹陷围道 ΓGamma(由 [ −𝑅, −𝜀][-R, -varepsilon]、小半圆弧 𝐶𝜀C_varepsilon(从 −𝜀-varepsilon 顺时针绕至 𝜀varepsilon)、[𝜀,𝑅][varepsilon, R] 及大半圆弧 𝐶𝑅C_R 构成).
围道内部处处全纯,由柯西定理:
∮Γ𝑒𝑖𝑧𝑧d𝑧=0oint_Gamma e^izzdz = 0 大圆弧 𝐶𝑅C_R 由约当引理衰减趋于 0;
小半圆弧 𝐶𝜀C_varepsilon 在原点展开:𝑒𝑖𝑧𝑧 =1+𝑖𝑧+…𝑧 =1𝑧 +𝑖 +…e^izz = 1 + iz + dots/z = 1/z + i + dots.参数化 𝑧 =𝜀𝑒𝑖𝜃z = varepsilon e^itheta(𝜃theta 从 𝜋pi 变到 00):
∫𝐶𝜀𝑒𝑖𝑧𝑧d𝑧=∫0𝜋1𝜀𝑒𝑖𝜃𝑖𝜀𝑒𝑖𝜃d𝜃=𝑖∫0𝜋d𝜃=−𝑖𝜋int_C_varepsilon e^izzdz = int_pi^0 1varepsilon e^itheta ivarepsilon e^ithetadtheta = i int_pi^0 dtheta = -ipi实轴两段合并(对负半轴作代换 𝑥 → −𝑡x -> -t):
∫−𝜀−𝑅𝑒𝑖𝑥𝑥d𝑥+∫𝑅𝜀𝑒𝑖𝑥𝑥d𝑥=∫𝑅𝜀𝑒𝑖𝑥−𝑒−𝑖𝑥𝑥d𝑥=2𝑖∫𝑅𝜀sin𝑥𝑥d𝑥int_-R^-varepsilon e^ixxdx + int_varepsilon^R e^ixxdx = int_varepsilon^R e^ix - e^-ixxdx = 2iint_varepsilon^R sin xxdx
取极限 𝜀 →0,𝑅 →∞varepsilon -> 0, R -> inf 代入柯西等式:
2𝑖∫∞0sin𝑥𝑥d𝑥−𝑖𝜋=0⟹∫∞0sin𝑥𝑥d𝑥=𝜋22i int_0^inf sin xxdx - ipi = 0 implies int_0^inf sin xxdx = 2例题 2:先判实际阶数,再用凹陷围道求 ∫∞0sin2𝑥𝑥2d𝑥int_0^inf sin^2 x/x^2dx
由偶函数对称性与二倍角降幂:
𝐼=∫∞0sin2𝑥𝑥2d𝑥=14∫+∞−∞1−cos2𝑥𝑥2d𝑥I = int_0^inf sin^2 xx^2dx = 14int_-inf^+inf 1 - cos 2xx^2dx 构造复变辅助函数 𝑓(𝑧) =1−𝑒2𝑖𝑧𝑧2f(z) = 1 - e^2izz^2,由于分子在 𝑧 =0z = 0 处有一阶零点,原点为 一阶极点.
在相同的上半平面凹陷围道上积分,由柯西定理 ∮Γ𝑓(𝑧)d𝑧 =0oint_Gamma f(z)dz = 0.
- 实轴两段实部给出 4𝐼4I;
- 大圆弧衰减为 0;
- 小圆弧上展开分子:1−(1+2𝑖𝑧−2𝑧2+…)𝑧2 = −2𝑖𝑧 +2 +…1 - (1 + 2iz - 2z^2 + dots)/z^2 = -2i/z + 2 + dots.
顺时针积分为 −2𝑖( −𝑖𝜋) = −2𝜋-2i(-ipi) = -2pi.
代回柯西方程:
4𝐼−2𝜋=0⟹𝐼=∫∞0sin2𝑥𝑥2d𝑥=𝜋24I - 2pi = 0 implies I = int_0^inf sin^2 xx^2dx = 2
范式五:矩形闭合围道与高斯傅里叶变换
例题:证明高斯函数傅里叶变换的自对偶性 ∫+∞−∞𝑒−𝜋𝑥2𝑒−2𝜋𝑖𝑥𝜉d𝑥 =𝑒−𝜋𝜉2int_-inf^+inf e^-pi x^2 e^-2pi i x xidx = e^-pi xi^2
设 𝜉 >0xi > 0.考察全纯整函数 𝑓(𝑧) =𝑒−𝜋𝑧2f(z) = e^-pi z^2.
构造顶点分别为 −𝑅,𝑅,𝑅 +𝑖𝜉, −𝑅 +𝑖𝜉-R, R, R + ixi, -R + ixi 的矩形闭合围道 𝛾𝑅gamma_R.由柯西积分定理:
∮𝛾𝑅𝑒−𝜋𝑧2d𝑧=0oint_gamma_R e^-pi z^2dz = 0 - 底边(𝑦 =0y = 0):∫𝑅−𝑅𝑒−𝜋𝑥2d𝑥𝑅→∞←←←←←←←←←→1int_-R^R e^-pi x^2dx R -> inf 1(标准高斯积分);
- 左右两侧竖直边(𝑥 = ±𝑅x = +/- R):由于被积函数含有模长衰减因子 𝑒−𝜋(𝑅2−𝑦2) ≤𝑒−𝜋𝑅2𝑒𝜋𝜉2e^-pi(R^2 - y^2) <= e^-pi R^2 e^pi xi^2,当 𝑅 →∞R -> inf 时指数衰减至 0;
顶边(𝑧 =𝑥 +𝑖𝜉,𝑥z = x + ixi, x 从 𝑅R 变到 −𝑅-R):
∫−𝑅𝑅𝑒−𝜋(𝑥+𝑖𝜉)2d𝑥=−∫𝑅−𝑅𝑒−𝜋(𝑥2−𝜉2+2𝑖𝑥𝜉)d𝑥=−𝑒𝜋𝜉2∫𝑅−𝑅𝑒−𝜋𝑥2𝑒−2𝜋𝑖𝑥𝜉d𝑥int_R^-R e^-pi(x + ixi)^2dx = -int_-R^R e^-pi(x^2 - xi^2 + 2ixxi)dx = -e^pixi^2int_-R^R e^-pi x^2 e^-2pi i xxidx
四边求和令 𝑅 →∞R -> inf:
1+0−𝑒𝜋𝜉2∫+∞−∞𝑒−𝜋𝑥2𝑒−2𝜋𝑖𝑥𝜉d𝑥+0=0⟹∫+∞−∞𝑒−𝜋𝑥2𝑒−2𝜋𝑖𝑥𝜉d𝑥=𝑒−𝜋𝜉21 + 0 - e^pixi^2int_-inf^+inf e^-pi x^2 e^-2pi i xxidx + 0 = 0 implies int_-inf^+inf e^-pi x^2 e^-2pi i xxidx = e^-pixi^2
讲义 3.7 的计算链:留数、边界项、实积分
下面三个例子将求留数步骤与讲义的应用连接起来。相应训练见综合题集。
二阶极点:有理积分与余弦积分
对 𝑓(𝑧) =𝑧2/(𝑧2 +1)2f(z)=z^2/(z^2+1)^2,上半平面只有二阶极点 𝑖i,且
Res(𝑓,𝑖)=dd𝑧𝑧2(𝑧+𝑖)2∣𝑧=𝑖=−𝑖4.Res(f,i)=left.mathrm dmathrm dz z^2(z+i)^2right|_z=i=-dfrac i4.大圆上 |𝑓(𝑧)| ≤𝑅2/(𝑅2 −1)2|f(z)|<= R^2/(R^2-1)^2,弧长为 𝜋𝑅pi R,故圆弧积分为 𝑂(𝑅−1)O(R^-1)。因此
∫∞−∞𝑥2(𝑥2+1)2d𝑥=2𝜋𝑖(−𝑖4)=𝜋2.int_-inf^infx^2(x^2+1)^2 ,mathrm dx =2pi ileft(-dfrac i4right)=dfracpi2.对余弦积分,改用 𝑔(𝑧) =𝑒𝑖𝑧/(𝑧2 +1)2g(z)=e^iz/(z^2+1)^2。必须对指数与剩余分母一起求导:
Res(𝑔,𝑖)=[𝑖𝑒𝑖𝑧(𝑧+𝑖)2−2𝑒𝑖𝑧(𝑧+𝑖)3]𝑧=𝑖=−𝑖2𝑒.Res(g,i) =left[ie^iz(z+i)^2-2e^iz(z+i)^3right]_z=i =-i2e.大圆弧由约当引理消失;留数定理得到复实轴积分 𝜋/𝑒pi/e,取实部即有
∫∞−∞cos𝑥(𝑥2+1)2d𝑥=𝜋𝑒.int_-inf^infcos x(x^2+1)^2 ,mathrm dx=dfracpi e.凹陷围道:辅助函数有极点,原正弦积分可以收敛
设 𝑎,𝑚 >0a,m>0,考虑 𝐼 =∫∞0sin(𝑚𝑥)/[𝑥(𝑥2 +𝑎2)] d𝑥I=int_0^infsin(mx)/[x(x^2+a^2)] ,mathrm dx。原实被积函数在零点趋于 𝑚/𝑎2m/a^2,但辅助函数 ℎ(𝑧) =𝑒𝑖𝑚𝑧/[𝑧(𝑧2 +𝑎2)]h(z)=e^imz/[z(z^2+a^2)] 在原点有单极点,留数为 1/𝑎21/a^2。
在上半平面用从 −𝜀-varepsilon 到 𝜀varepsilon 的顺时针小半圆绕开原点,小弧贡献为 −𝑖𝜋/𝑎2-ipi/a^2。围道内只含 𝑖𝑎ia,其留数为 −𝑒−𝑚𝑎/(2𝑎2)-e^-ma/(2a^2)。实轴两段合为 2𝑖𝐼2iI,大圆弧由约当引理消失,所以
2𝑖𝐼−𝑖𝜋𝑎2=−𝜋𝑖𝑎2𝑒−𝑚𝑎,𝐼=𝜋2𝑎2(1−𝑒−𝑚𝑎).2iI-ia^2=-pi ia^2e^-ma, qquad I=2a^2(1-e^-ma).当 𝑚 →0+m->0^+ 时结果趋零;当 𝑚m 很大时趋于 𝜋/(2𝑎2)pi/(2a^2)。这些极限也可用来核对小圆弧的符号。
傅里叶积分与 Yukawa 核
对实数 𝑎 ≠0ane0、𝑝 ∈ℝpinmathbb R,应分别处理 𝑝p 的正、负、零三种情形。𝑝 >0p>0 向上闭合,𝑝 <0p<0 向下闭合并计入顺时针负号,𝑝 =0p=0 用普通有理函数的弧长估计;三者合并为
∫∞−∞𝑒𝑖𝑝𝑥𝑥2+𝑎2d𝑥=𝜋|𝑎|𝑒−|𝑎𝑝|.int_-inf^infe^ipxx^2+a^2 ,mathrm dx =dfrac|a|e^-|ap|.三维动量积分的径向核(进阶)
讲义采用没有 (2𝜋)−3(2pi)^-3 前因子的傅里叶约定。设 𝑟 =|𝐫| >0,𝜇 >0r=|mathbf r|>0,mu>0,将极轴沿 𝐫mathbf r 取向,先完成角积分,再按球形截断取径向反常积分:
𝐼(𝑟)=∫ℝ3𝑒𝑖𝐤⋅𝐫𝑘2+𝜇2d3𝑘=4𝜋𝑟∫∞0𝑘sin(𝑘𝑟)𝑘2+𝜇2d𝑘=2𝜋𝑖𝑟∫∞−∞𝑘𝑒𝑖𝑘𝑟𝑘2+𝜇2d𝑘.I(r)=int_mathbb R^3e^imathbf k*mathbf rk^2+mu^2 ,mathrm d^3k =4rint_0^infksin(kr)k^2+mu^2 ,mathrm dk =2irint_-inf^infk e^ikrk^2+mu^2 ,mathrm dk. 最后一步的全实轴积分按对称截断理解。上半平面唯一极点为 𝑖𝜇imu,留数为 𝑒−𝜇𝑟/2e^-mu r/2;有理因子 𝑘/(𝑘2 +𝜇2) =𝑂(1/𝑘)k/(k^2+mu^2)=O(1/k) 满足约当引理条件。因此
𝐼(𝑟)=2𝜋𝑖𝑟(2𝜋𝑖𝑒−𝜇𝑟2)=2𝜋2𝑟𝑒−𝜇𝑟.I(r)=2irleft(2pi ie^-mu r2right) =2pi^2re^-mu r. 原三维积分并非绝对收敛,不能任意交换积分顺序;这里采用上述球形截断解释,也可采用相容的正则化或分布解释。若傅里叶逆变换含 (2𝜋)−3(2pi)^-3,则径向核为 𝑒−𝜇𝑟/(4𝜋𝑟)e^-mu r/(4pi r)。这个指数衰减给出了 Yukawa 势的短程径向形式。
4. 理论物理中的深度连接:因果律与极点处方
在相对论与量子物理中,「结果不能先于原因发生」是物理世界的基本铁律——因果律(Causality).这一物理原理在复分析中表现为极其深刻的数学结构:
4.1 柯西主值(P.V.)与索霍茨基 - 普莱姆利(Sokhotski–Plemelj) 公式
当极点落在实轴上时,物理学家在虚轴方向加上微小的扰动 𝑥 ±𝑖𝜖x +/- iepsilon(𝜖 →0+epsilon -> 0^+)给出绕过实轴极点的处方:
lim𝜖→0+1𝑥±𝑖𝜖=P1𝑥∓𝑖𝜋𝛿(𝑥)lim_epsilon -> 0^+ 1x +/- iepsilon = P1x mp ipi delta(x)其中 PP 代表柯西主值(Cauchy Principal Value),𝛿(𝑥)delta(x) 是狄拉克 Delta 函数。等式在分布意义下成立,不能当作 𝑥 =0x=0 处的普通函数等式。用于具体响应函数时,主值项与 Delta 项的物理意义还取决于傅里叶约定和响应模型。
4.2 克拉默斯 - 克勒尼希(Kramers–Kronig) 色散关系
在光学与凝聚态物理中,物质的复介电常数 𝜀(𝜔) =𝜀1(𝜔) +𝑖𝜀2(𝜔)varepsilon(omega) = varepsilon_1(omega) + ivarepsilon_2(omega)。采用 𝜒(𝜔) =∫∞0𝜒(𝑡)𝑒𝑖𝜔𝑡 d𝑡chi(omega)=int_0^infchi(t)e^iomega t ,mathrm dt 的约定,对稳定且满足相应收敛条件的因果响应,响应函数在上半复平面解析。若高频极限为 11,且减去该极限后的大圆弧贡献消失,再结合实时间响应的对称性,可由柯西积分公式得到:
𝜀1(𝜔)−1=2𝜋P∫∞0𝜔′𝜀2(𝜔′)𝜔′2−𝜔2d𝜔′varepsilon_1(omega) - 1 = 2piPint_0^inf omega' varepsilon_2(omega')omega'^2 - omega^2domega'这说明实部与虚部由全频段的积分关系联系起来。某一频段内可以有很小的吸收而仍有色散,不能把该关系解释为每个频率上必须同时出现非零吸收。
4.3 量子场论费曼(Feynman)𝑖𝜖iepsilon 处方
在求解相对论 Klein-Gordon 波动方程自由传播子(格林函数)时,奇点落在动量实轴上(𝑝2 −𝑚2 =0p^2 - m^2 = 0).为了确保粒子向前传播、反粒子向后传播的时间因果律,费曼规定将分母微扰为:
Δ𝐹(𝑝)=1𝑝2−𝑚2+𝑖𝜖Delta_F(p) = 1p^2 - m^2 + iepsilon在度规约定 𝑝2 =(𝑝0)2 −𝐩2p^2=(p^0)^2-mathbf p^2 下,令 𝐸𝐩 =√𝐩2+𝑚2 >0E_mathbf p=mathbf p^2+m^2>0,两能量极点位于 𝑝0 ≈𝐸𝐩 −𝑖𝜖/(2𝐸𝐩)p^0~ E_mathbf p-iepsilon/(2E_mathbf p) 和 𝑝0 ≈ −𝐸𝐩 +𝑖𝜖/(2𝐸𝐩)p^0~-E_mathbf p+iepsilon/(2E_mathbf p)。这一处方定义时间排序的费曼传播子;推迟与超前格林函数采用各自的极点处方,不能与费曼传播子混同。
5. 学习衔接
知识小测与巩固练习
补充说明
课后核心技能自测点
- 用商式法则计算 𝑓(𝑧) =𝑒𝑧𝑧4−1f(z) = e^z/z^4 - 1 在 𝑧 =𝑖z = i 处的留数;
- 计算实定积分 ∫∞−∞𝑥2(𝑥2+1)(𝑥2+4)d𝑥int_-inf^inf x^2/(x^2+1)(x^2+4)dx;
- 说明为什么约当引理中圆弧半径仅需一次衰减(𝑀𝑅 →0M_R -> 0),而普通有理分式需要二次衰减(𝑀𝑅 ∼1/𝑅2M_R ~ 1/R^2)?
参考资料
本节奇点、留数定义与计算例题对照用户提供的《3.5 奇点分类、3.6 留数定理发现之旅与 3.7 留数定理》,即仓库 References/数学物理方法课件/mmplecture 3.5-3.7.md。求留数的局部方法和扩充复平面的符号约定见方法选择与计算,讲义思考题的训练与解析见综合题集。
本页面最近更新:2026/10/12 13:04:22,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用