复变积分与柯西定理 物理问题引入:解析函数的「全息摄影」特性 在实数多元微积分中,如果你仅仅知道一个二元函数 𝑢 ( 𝑥 , 𝑦 ) u ( x , y ) u(x, y) 在某个区域边界曲线 𝜕 𝐷 ∂ D d D 上的取值,你通常根本无法断言它在区域内部的数值——除非你事先完全知晓其内部的泊松方程源项分布.此外,一个实函数即使一阶可导,它的二阶导数也可能处处不存在.
然而,在复分析中,解析函数展现出了近乎「全息摄影(Holography)」的惊人特质:
闭合环路积分为零 :只要区域内部处处全纯解析无奇点,沿该区域内任意闭合回路的复线积分 恒等于零 ;边界决定内部一切 :只要知道解析函数在边界封闭围道上的值,内部任意一点的值与各阶导数都可以通过一个边界围道积分被 毫厘不差地完全重构出来 (柯西积分公式);一阶解析即无限阶可微 :只要复一阶导数存在,各阶高阶导数自动存在且全部全纯!这种「边界完全编码内部自由度」的特性,不仅在宏观上是真空静电学拉普拉斯方程解的唯一性体现,更在现代高能理论中深刻启发了引力与量子场论的 AdS/CFT 全息对应 思想.
1. 复平面围道积分定义与性质 1.1 复变线积分的定义 设 𝐶 C C 为复平面上一条连续逐段光滑的定向曲线,参数化方程为 𝑧 ( 𝑡 ) = 𝑥 ( 𝑡 ) + 𝑖 𝑦 ( 𝑡 ) z ( t ) = x ( t ) + i y ( t ) z(t) = x(t) + iy(t) (𝑡 ∈ [ 𝑎 , 𝑏 ] t ∈ [ a , b ] t in [a, b] ). 复变函数 𝑓 ( 𝑧 ) = 𝑢 ( 𝑥 , 𝑦 ) + 𝑖 𝑣 ( 𝑥 , 𝑦 ) f ( z ) = u ( x , y ) + i v ( x , y ) f(z) = u(x, y) + i v(x, y) 沿曲线 𝐶 C C 的线积分定义为:
∫ 𝐶 𝑓 ( 𝑧 ) d 𝑧 = ∫ 𝐶 ( 𝑢 + 𝑖 𝑣 ) ( d 𝑥 + 𝑖 d 𝑦 ) = ∫ 𝐶 ( 𝑢 d 𝑥 − 𝑣 d 𝑦 ) + 𝑖 ∫ 𝐶 ( 𝑣 d 𝑥 + 𝑢 d 𝑦 ) ∫ C f ( z ) d z = ∫ C ( u + i v ) ( d x + i d y ) = ∫ C ( u d x − v d y ) + i ∫ C ( v d x + u d y ) int_C f(z) dz = int_C (u + iv)(dx + idy) = int_C (u ,dx - v ,dy) + i int_C (v ,dx + u ,dy) 一个复积分在实质上完全等价于 两个定向实平面第二类线积分 的线性组合.
1.2 基础奇点整幂积分:留数理论的种子 考察以 𝑧 0 z 0 z_0 为圆心、半径为 𝑅 R R 的逆时针正向圆周围道 𝐶 C C :𝑧 − 𝑧 0 = 𝑅 𝑒 𝑖 𝜃 z − z 0 = R e i θ z - z_0 = R e^itheta (𝜃 ∈ [ 0 , 2 𝜋 ] θ ∈ [ 0 , 2 π ] theta in [0, 2pi] ),d 𝑧 = 𝑖 𝑅 𝑒 𝑖 𝜃 d 𝜃 d z = i R e i θ d θ dz = i R e^ithetadtheta . 计算整幂函数 ( 𝑧 − 𝑧 0 ) 𝑛 ( z − z 0 ) n (z - z_0)^n (𝑛 ∈ ℤ n ∈ Z n in Z )的闭回路积分:
∮ 𝐶 ( 𝑧 − 𝑧 0 ) 𝑛 d 𝑧 = ∫ 2 𝜋 0 𝑅 𝑛 𝑒 𝑖 𝑛 𝜃 ( 𝑖 𝑅 𝑒 𝑖 𝜃 ) d 𝜃 = 𝑖 𝑅 𝑛 + 1 ∫ 2 𝜋 0 𝑒 𝑖 ( 𝑛 + 1 ) 𝜃 d 𝜃 ∮ C ( z − z 0 ) n d z = ∫ 0 2 π R n e i n θ ( i R e i θ ) d θ = i R n + 1 ∫ 0 2 π e i ( n + 1 ) θ d θ oint_C (z - z_0)^n dz = int_0^2pi R^n e^intheta left(i R e^ithetaright)dtheta = i R^n+1 int_0^2pi e^i(n+1)thetadtheta 当 𝑛 = − 1 n = − 1 n = -1 时 :𝑛 + 1 = 0 ⟹ 𝑒 𝑖 0 = 1 n + 1 = 0 ⟹ e i 0 = 1 n + 1 = 0 implies e^i 0 = 1 :
∮ 𝐶 1 𝑧 − 𝑧 0 d 𝑧 = 𝑖 𝑅 0 ∫ 2 𝜋 0 1 d 𝜃 = 2 𝜋 𝑖 ∮ C 1 z − z 0 d z = i R 0 ∫ 0 2 π 1 d θ = 2 π i oint_C 1z - z_0dz = i R^0 int_0^2pi 1 ,dtheta = 2pi i 当 𝑛 ≠ − 1 n ≠ − 1 n != -1 时 :
∫ 2 𝜋 0 𝑒 𝑖 ( 𝑛 + 1 ) 𝜃 d 𝜃 = 𝑒 𝑖 ( 𝑛 + 1 ) 𝜃 𝑖 ( 𝑛 + 1 ) ∣ 2 𝜋 0 = 1 − 1 𝑖 ( 𝑛 + 1 ) = 0 ∫ 0 2 π e i ( n + 1 ) θ d θ = e i ( n + 1 ) θ i ( n + 1 ) | 0 2 π = 1 − 1 i ( n + 1 ) = 0 int_0^2pi e^i(n+1)thetadtheta = left.e^i(n+1)thetai(n+1)right|_0^2pi = 1 - 1i(n+1) = 0 核心结论 :
∮ 𝐶 ( 𝑧 − 𝑧 0 ) 𝑛 d 𝑧 = { 2 𝜋 𝑖 , 𝑛 = − 1 0 , 𝑛 ∈ ℤ , 𝑛 ≠ − 1 ∮ C ( z − z 0 ) n d z = { 2 π i , n = − 1 0 , n ∈ Z , n ≠ − 1 oint_C (z - z_0)^n dz = cases 2pi i, & n = -1 0, & n in Z, n != -1 cases 围绕奇点的幂函数闭回路积分中,所有其他整幂次全部积分归零,唯独负一次幂项 1 𝑧 − 𝑧 0 1 z − z 0 1/z-z_0 贡献了非零的拓扑荷 2 𝜋 𝑖 2 π i 2pi i !
例题:非解析函数的闭回路积分 计算共轭函数 𝑓 ( 𝑧 ) = 𝑧 ∗ f ( z ) = z ∗ f(z) = z^* 沿单位圆周 | 𝑧 | = 1 | z | = 1 |z| = 1 的积分: 参数化:𝑧 = 𝑒 𝑖 𝜃 , 𝑧 ∗ = 𝑒 − 𝑖 𝜃 , d 𝑧 = 𝑖 𝑒 𝑖 𝜃 d 𝜃 z = e i θ , z ∗ = e − i θ , d z = i e i θ d θ z = e^itheta, z^* = e^-itheta, dz = i e^ithetadtheta :
∮ | 𝑧 | = 1 𝑧 ∗ d 𝑧 = ∫ 2 𝜋 0 𝑒 − 𝑖 𝜃 ( 𝑖 𝑒 𝑖 𝜃 ) d 𝜃 = 𝑖 ∫ 2 𝜋 0 d 𝜃 = 2 𝜋 𝑖 ≠ 0 ∮ | z | = 1 z ∗ d z = ∫ 0 2 π e − i θ ( i e i θ ) d θ = i ∫ 0 2 π d θ = 2 π i ≠ 0 oint_|z|=1 z^* ,dz = int_0^2pi e^-itheta(i e^itheta)dtheta = i int_0^2pi dtheta = 2pi i != 0 闭路积分不为零直接从积分判据上证明了 𝑓 ( 𝑧 ) = 𝑧 ∗ f ( z ) = z ∗ f(z) = z^* 在全复平面上处处不解析.
2. 柯西积分定理(Cauchy's Integral Theorem) 2.1 定理陈述与格林公式证明 定理(Cauchy-Goursat 定理) :若函数 𝑓 ( 𝑧 ) f ( z ) f(z) 在单连通区域 𝐷 D D 内全纯解析,且在闭区域 ¯ 𝐷 = 𝐷 ∪ 𝜕 𝐷 D ¯ = D ∪ ∂ D D = D cup d D 上连续,则沿 𝐷 D D 内任意一条简单闭合曲线 𝐶 C C 的线积分恒等于零:
∮ 𝐶 𝑓 ( 𝑧 ) d 𝑧 = 0 ∮ C f ( z ) d z = 0 oint_C f(z) dz = 0 严密证明(基于二维格林公式) : 将复积分拆解为实部与虚部:
∮ 𝐶 𝑓 ( 𝑧 ) d 𝑧 = ∮ 𝐶 ( 𝑢 d 𝑥 − 𝑣 d 𝑦 ) + 𝑖 ∮ 𝐶 ( 𝑣 d 𝑥 + 𝑢 d 𝑦 ) ∮ C f ( z ) d z = ∮ C ( u d x − v d y ) + i ∮ C ( v d x + u d y ) oint_C f(z)dz = oint_C (u ,dx - v ,dy) + i oint_C (v ,dx + u ,dy) 应用平面 Green 公式 ∮ 𝐶 𝑃 d 𝑥 + 𝑄 d 𝑦 = ∬ 𝑆 ( 𝜕 𝑄 𝜕 𝑥 − 𝜕 𝑃 𝜕 𝑦 ) d 𝑥 d 𝑦 ∮ C P d x + Q d y = ∬ S ( ∂ Q ∂ x − ∂ P ∂ y ) d x d y oint_C Pdx + Qdy = iint_S left(d Q/d x - d P/d yright)dxdy :
实部线积分:𝑃 = 𝑢 , 𝑄 = − 𝑣 P = u , Q = − v P = u, Q = -v :
∮ 𝐶 ( 𝑢 d 𝑥 − 𝑣 d 𝑦 ) = ∬ 𝑆 ( − 𝜕 𝑣 𝜕 𝑥 − 𝜕 𝑢 𝜕 𝑦 ) d 𝑥 d 𝑦 ∮ C ( u d x − v d y ) = ∬ S ( − ∂ v ∂ x − ∂ u ∂ y ) d x d y oint_C (u ,dx - v ,dy) = iint_S left(-d vd x - d ud yright)dxdy 虚部线积分:𝑃 = 𝑣 , 𝑄 = 𝑢 P = v , Q = u P = v, Q = u :
∮ 𝐶 ( 𝑣 d 𝑥 + 𝑢 d 𝑦 ) = ∬ 𝑆 ( 𝜕 𝑢 𝜕 𝑥 − 𝜕 𝑣 𝜕 𝑦 ) d 𝑥 d 𝑦 ∮ C ( v d x + u d y ) = ∬ S ( ∂ u ∂ x − ∂ v ∂ y ) d x d y oint_C (v ,dx + u ,dy) = iint_S left(d ud x - d vd yright)dxdy 由于 𝑓 ( 𝑧 ) f ( z ) f(z) 解析,柯西 - 黎曼条件 𝑢 𝑥 = 𝑣 𝑦 u x = v y u_x = v_y 与 𝑢 𝑦 = − 𝑣 𝑥 u y = − v x u_y = -v_x 严格成立,两面积分被积表达式处处恒等于零!因此:
∮ 𝐶 𝑓 ( 𝑧 ) d 𝑧 = 0 + 𝑖 0 = 0 ∮ C f ( z ) d z = 0 + i 0 = 0 oint_C f(z)dz = 0 + i0 = 0 证毕.
2.2 闭路变形原理与多连通区域柯西定理 闭路变形原理(Contour Deformation) : 若函数 𝑓 ( 𝑧 ) f ( z ) f(z) 在两条闭合曲线 𝐶 1 C 1 C_1 与 𝐶 2 C 2 C_2 所夹的区域内处处全纯解析,则将外围道 𝐶 1 C 1 C_1 连续形变为内围道 𝐶 2 C 2 C_2 并不改变回路积分的数值:
∮ 𝐶 1 𝑓 ( 𝑧 ) d 𝑧 = ∮ 𝐶 2 𝑓 ( 𝑧 ) d 𝑧 ∮ C 1 f ( z ) d z = ∮ C 2 f ( z ) d z oint_C_1 f(z)dz = oint_C_2 f(z)dz 多连通区域柯西定理(复合闭路定理) : 若 𝑓 ( 𝑧 ) f ( z ) f(z) 在多连通区域 ¯ 𝐷 D ¯ D 内全纯解析,外边界为 𝐿 0 L 0 L_0 (逆时针正向),内部包含有限个洞,内边界依次为 𝐿 1 , 𝐿 2 , … , 𝐿 𝑘 L 1 , L 2 , … , L k L_1, L_2, dots, L_k (顺时针):
∮ 𝐿 0 𝑓 ( 𝑧 ) d 𝑧 = 𝑚 ∑ 𝑘 = 1 ∮ 𝐿 𝑘 𝑓 ( 𝑧 ) d 𝑧 ∮ L 0 f ( z ) d z = ∑ k = 1 m ∮ L k f ( z ) d z oint_L_0 f(z)dz = sum_k=1^m oint_L_k f(z)dz 3. 柯西积分公式与高阶导数推广
定理 :设区域 𝐺 G G 的边界为简单闭曲线 𝐶 C C .若 𝑓 ( 𝑧 ) f ( z ) f(z) 在闭区域 ¯ 𝐺 G ¯ G 上全纯解析,则对于 𝐺 G G 内部的任意点 𝑧 0 z 0 z_0 :
𝑓 ( 𝑧 0 ) = 1 2 𝜋 𝑖 ∮ 𝐶 𝑓 ( 𝑧 ) 𝑧 − 𝑧 0 d 𝑧 f ( z 0 ) = 1 2 π i ∮ C f ( z ) z − z 0 d z f(z_0) = 12pi i oint_C f(z)z - z_0dz 证明 : 以 𝑧 0 z 0 z_0 为中心、𝜌 ρ rho 为半径作充分小的逆时针圆周围道 𝐶 𝜌 : | 𝑧 − 𝑧 0 | = 𝜌 C ρ : | z − z 0 | = ρ C_rho: |z - z_0| = rho .由闭路变形原理:
∮ 𝐶 𝑓 ( 𝑧 ) 𝑧 − 𝑧 0 d 𝑧 = ∮ 𝐶 𝜌 𝑓 ( 𝑧 ) 𝑧 − 𝑧 0 d 𝑧 ∮ C f ( z ) z − z 0 d z = ∮ C ρ f ( z ) z − z 0 d z oint_C f(z)z - z_0dz = oint_C_rho f(z)z - z_0dz 在圆周 𝐶 𝜌 C ρ C_rho 上参数化:𝑧 = 𝑧 0 + 𝜌 𝑒 𝑖 𝜃 , d 𝑧 = 𝑖 𝜌 𝑒 𝑖 𝜃 d 𝜃 z = z 0 + ρ e i θ , d z = i ρ e i θ d θ z = z_0 + rho e^itheta,dz = irho e^ithetadtheta .将分子拆分为常数项与微小增量:𝑓 ( 𝑧 ) = 𝑓 ( 𝑧 0 ) + [ 𝑓 ( 𝑧 ) − 𝑓 ( 𝑧 0 ) ] f ( z ) = f ( z 0 ) + [ f ( z ) − f ( z 0 ) ] f(z) = f(z_0) + [f(z) - f(z_0)] :
∮ 𝐶 𝜌 𝑓 ( 𝑧 ) 𝑧 − 𝑧 0 d 𝑧 = ∫ 2 𝜋 0 𝑓 ( 𝑧 0 ) 𝜌 𝑒 𝑖 𝜃 𝑖 𝜌 𝑒 𝑖 𝜃 d 𝜃 + ∫ 2 𝜋 0 𝑓 ( 𝑧 0 + 𝜌 𝑒 𝑖 𝜃 ) − 𝑓 ( 𝑧 0 ) 𝜌 𝑒 𝑖 𝜃 𝑖 𝜌 𝑒 𝑖 𝜃 d 𝜃 ∮ C ρ f ( z ) z − z 0 d z = ∫ 0 2 π f ( z 0 ) ρ e i θ i ρ e i θ d θ + ∫ 0 2 π f ( z 0 + ρ e i θ ) − f ( z 0 ) ρ e i θ i ρ e i θ d θ oint_C_rho f(z)z - z_0dz = int_0^2pi f(z_0)rho e^itheta irho e^ithetadtheta + int_0^2pi f(z_0 + rho e^itheta) - f(z_0)rho e^itheta irho e^ithetadtheta = 2 𝜋 𝑖 𝑓 ( 𝑧 0 ) + 𝑖 ∫ 2 𝜋 0 [ 𝑓 ( 𝑧 0 + 𝜌 𝑒 𝑖 𝜃 ) − 𝑓 ( 𝑧 0 ) ] d 𝜃 = 2 π i f ( z 0 ) + i ∫ 0 2 π [ f ( z 0 + ρ e i θ ) − f ( z 0 ) ] d θ = 2pi i f(z_0) + i int_0^2pi [f(z_0 + rho e^itheta) - f(z_0)] dtheta 由于 𝑓 ( 𝑧 ) f ( z ) f(z) 在 𝑧 0 z 0 z_0 处连续,当 𝜌 → 0 ρ → 0 rho -> 0 时,最大增量 m a x | 𝑓 ( 𝑧 ) − 𝑓 ( 𝑧 0 ) | → 0 max | f ( z ) − f ( z 0 ) | → 0 max |f(z) - f(z_0)| -> 0 ,第二项积分严格收敛于 0.两边除以 2 𝜋 𝑖 2 π i 2pi i 即得定理.
3.2 柯西积分高阶导数推广公式 利用差商极限并令 𝜉 → 0 ξ → 0 xi -> 0 :
𝑓 ′ ( 𝑧 0 ) = l i m 𝜉 → 0 𝑓 ( 𝑧 0 + 𝜉 ) − 𝑓 ( 𝑧 0 ) 𝜉 = 1 2 𝜋 𝑖 ∮ 𝐶 l i m 𝜉 → 0 1 𝜉 [ 1 𝑧 − 𝑧 0 − 𝜉 − 1 𝑧 − 𝑧 0 ] 𝑓 ( 𝑧 ) d 𝑧 = 1 2 𝜋 𝑖 ∮ 𝐶 𝑓 ( 𝑧 ) ( 𝑧 − 𝑧 0 ) 2 d 𝑧 f ′ ( z 0 ) = lim ξ → 0 f ( z 0 + ξ ) − f ( z 0 ) ξ = 1 2 π i ∮ C lim ξ → 0 1 ξ [ 1 z − z 0 − ξ − 1 z − z 0 ] f ( z ) d z = 1 2 π i ∮ C f ( z ) ( z − z 0 ) 2 d z f'(z_0) = lim_xi -> 0 f(z_0 + xi) - f(z_0)xi = 12pi i oint_C lim_xi -> 0 1xileft[1z - z_0 - xi - 1z - z_0right] f(z)dz = 12pi i oint_C f(z)(z - z_0)^2dz 通过数学归纳法,对任意阶导数 𝑛 ∈ ℕ n ∈ N n in N :
𝑓 ( 𝑛 ) ( 𝑧 0 ) = 𝑛 ! 2 𝜋 𝑖 ∮ 𝐶 𝑓 ( 𝑧 ) ( 𝑧 − 𝑧 0 ) 𝑛 + 1 d 𝑧 f ( n ) ( z 0 ) = n ! 2 π i ∮ C f ( z ) ( z − z 0 ) n + 1 d z f^(n)(z_0) = n!2pi i oint_C f(z)(z - z_0)^n+1dz 物理与数学深刻推论 :
可导即无穷次可微 :若 𝑓 ( 𝑧 ) f ( z ) f(z) 一阶复可导,则它的各阶导数全部存在且依然是解析函数;莫雷拉定理(Morera's Theorem) :柯西积分定理的逆定理——若连续函数在区域内沿任意闭回路积分为零,则它必为全纯解析函数;刘维尔定理与代数基本定理 :全复平面上有界整函数必为常数;由此可直接证明一元 𝑛 n n 次代数方程必有 𝑛 n n 个复根;极大模原理与恩绍定理(Earnshaw's Theorem) :解析函数的模 | 𝑓 ( 𝑧 ) | | f ( z ) | |f(z)| 不可能在区域内部取得极大值(只能在边界取得).在物理上,静电势 𝜙 ϕ phi 满足拉普拉斯方程,因此在无电荷空间中 不存在任何静止电荷的稳定平衡悬浮点 .4. 典型闭围道积分计算实战 利用柯西积分公式及高阶导数公式,闭回路积分被完全转化为局域求导运算:
∮ 𝐶 𝑓 ( 𝑧 ) ( 𝑧 − 𝑧 0 ) 𝑛 + 1 d 𝑧 = 2 𝜋 𝑖 𝑛 ! 𝑓 ( 𝑛 ) ( 𝑧 0 ) ∮ C f ( z ) ( z − z 0 ) n + 1 d z = 2 π i n ! f ( n ) ( z 0 ) oint_C f(z)(z - z_0)^n+1dz = 2pi in! f^(n)(z_0) 例题 1:部分分式法分解多奇点 计算回路积分:
𝐼 = ∮ | 𝑧 | = 2 𝑒 𝑖 𝑧 𝑧 2 + 1 d 𝑧 I = ∮ | z | = 2 e i z z 2 + 1 d z I = oint_|z|=2 e^izz^2 + 1 ,dz 解 :被积函数奇点为 𝑧 = ± 𝑖 z = ± i z = +/- i ,模长为 1 < 2 1 < 2 1 < 2 ,两奇点均在闭围道内部.部分分式拆分:
𝑒 𝑖 𝑧 𝑧 2 + 1 = 𝑒 𝑖 𝑧 ( 𝑧 − 𝑖 ) ( 𝑧 + 𝑖 ) = 1 2 𝑖 ( 𝑒 𝑖 𝑧 𝑧 − 𝑖 − 𝑒 𝑖 𝑧 𝑧 + 𝑖 ) e i z z 2 + 1 = e i z ( z − i ) ( z + i ) = 1 2 i ( e i z z − i − e i z z + i ) e^izz^2 + 1 = e^iz(z - i)(z + i) = 12ileft(e^izz - i - e^izz + iright) 分别应用柯西积分公式(𝑓 ( 𝑧 ) = 𝑒 𝑖 𝑧 f ( z ) = e i z f(z) = e^iz 处处解析):
𝐼 = 1 2 𝑖 [ 2 𝜋 𝑖 𝑓 ( 𝑖 ) − 2 𝜋 𝑖 𝑓 ( − 𝑖 ) ] = 𝜋 ( 𝑒 𝑖 ( 𝑖 ) − 𝑒 𝑖 ( − 𝑖 ) ) = 𝜋 ( 𝑒 − 1 − 𝑒 ) I = 1 2 i [ 2 π i f ( i ) − 2 π i f ( − i ) ] = π ( e i ( i ) − e i ( − i ) ) = π ( e − 1 − e ) I = 12i left[2pi i f(i) - 2pi i f(-i)right] = pi left(e^i(i) - e^i(-i)right) = pi(e^-1 - e) 例题 2:高阶极点与高阶导数公式 计算积分:
𝐼 = ∮ | 𝑧 − 1 | = 1 s i n ( 𝜋 𝑧 4 ) ( 𝑧 2 − 1 ) 2 d 𝑧 I = ∮ | z − 1 | = 1 sin ( π z 4 ) ( z 2 − 1 ) 2 d z I = oint_|z-1|=1 sinleft(pi z/4right)(z^2 - 1)^2dz 解 :分母因式分解为 ( 𝑧 − 1 ) 2 ( 𝑧 + 1 ) 2 ( z − 1 ) 2 ( z + 1 ) 2 (z - 1)^2 (z + 1)^2 .积分圆周是以 1 1 1 为圆心、半径为 1 1 1 的圆 | 𝑧 − 1 | = 1 | z − 1 | = 1 |z - 1| = 1 .
奇点 𝑧 = 1 z = 1 z = 1 位于围道内部(二阶极点); 奇点 𝑧 = − 1 z = − 1 z = -1 位于围道外部(距离为 2 > 1 2 > 1 2 > 1 ). 将围道内全纯的部分归入分子函数: 𝑓 ( 𝑧 ) = s i n ( 𝜋 𝑧 4 ) ( 𝑧 + 1 ) 2 f ( z ) = sin ( π z 4 ) ( z + 1 ) 2 f(z) = sinleft(pi z/4right)(z + 1)^2 由高阶柯西公式(𝑛 = 1 n = 1 n = 1 ):
𝐼 = ∮ | 𝑧 − 1 | = 1 𝑓 ( 𝑧 ) ( 𝑧 − 1 ) 2 d 𝑧 = 2 𝜋 𝑖 1 ! 𝑓 ′ ( 1 ) I = ∮ | z − 1 | = 1 f ( z ) ( z − 1 ) 2 d z = 2 π i 1 ! f ′ ( 1 ) I = oint_|z-1|=1 f(z)(z - 1)^2dz = 2pi i1! f'(1) 计算导数:
𝑓 ′ ( 𝑧 ) = 𝜋 4 c o s ( 𝜋 𝑧 4 ) ( 𝑧 + 1 ) 2 − 2 ( 𝑧 + 1 ) s i n ( 𝜋 𝑧 4 ) ( 𝑧 + 1 ) 4 f ′ ( z ) = π 4 cos ( π z 4 ) ( z + 1 ) 2 − 2 ( z + 1 ) sin ( π z 4 ) ( z + 1 ) 4 f'(z) = pi/4cosleft(pi z/4right)(z + 1)^2 - 2(z + 1)sinleft(pi z/4right)(z + 1)^4 代入 𝑧 = 1 z = 1 z = 1 :
𝑓 ′ ( 1 ) = 𝜋 4 ( √ 2 2 ) ( 4 ) − 2 ( 2 ) ( √ 2 2 ) 1 6 = √ 2 𝜋 − 4 √ 2 3 2 = √ 2 8 ( 𝜋 4 − 1 ) f ′ ( 1 ) = π 4 ( 2 2 ) ( 4 ) − 2 ( 2 ) ( 2 2 ) 16 = 2 π − 4 2 32 = 2 8 ( π 4 − 1 ) f'(1) = pi/4left(22right)(4) - 2(2)left(22right)16 = 2pi - 4232 = 28left(4 - 1right) 最终得到:
𝐼 = 2 𝜋 𝑖 ⋅ √ 2 8 ( 𝜋 4 − 1 ) = √ 2 𝜋 𝑖 4 ( 𝜋 4 − 1 ) I = 2 π i ⋅ 2 8 ( π 4 − 1 ) = 2 π i 4 ( π 4 − 1 ) I = 2pi i * 28left(4 - 1right) = 2pi i4left(4 - 1right) 例题 3:常微分方程与高阶柯西公式综合题 若函数 𝑓 ( 𝑧 ) f ( z ) f(z) 在 | 𝑧 | < 1 | z | < 1 |z| < 1 内解析,且满足二阶线性常微分方程:
( 1 − 𝑧 2 ) 𝑓 ″ ( 𝑧 ) − 2 𝑧 𝑓 ′ ( 𝑧 ) + 4 2 𝑓 ( 𝑧 ) = 0 ( 1 − z 2 ) f ″ ( z ) − 2 z f ′ ( z ) + 42 f ( z ) = 0 (1 - z^2)f''(z) - 2z f'(z) + 42 f(z) = 0 同时已知初值 𝑓 ′ ( 0 ) = 3 f ′ ( 0 ) = 3 f'(0) = 3 .计算回路积分:
𝐼 = ∮ | 𝑧 | = 0 . 1 𝑓 ( 𝑧 ) 𝑧 4 d 𝑧 I = ∮ | z | = 0.1 f ( z ) z 4 d z I = oint_|z|=0.1 f(z)z^4dz 解 :
由柯西高阶导数公式,当 𝑧 0 = 0 , 𝑛 = 3 z 0 = 0 , n = 3 z_0 = 0, n = 3 时:
𝐼 = ∮ | 𝑧 | = 0 . 1 𝑓 ( 𝑧 ) 𝑧 4 d 𝑧 = 2 𝜋 𝑖 3 ! 𝑓 ‴ ( 0 ) = 𝜋 𝑖 3 𝑓 ‴ ( 0 ) I = ∮ | z | = 0.1 f ( z ) z 4 d z = 2 π i 3 ! f ‴ ( 0 ) = π i 3 f ‴ ( 0 ) I = oint_|z|=0.1 f(z)z^4dz = 2pi i3! f'''(0) = pi i3 f'''(0) 对原微分方程两边关于复变量 𝑧 z z 求导:
d d 𝑧 [ ( 1 − 𝑧 2 ) 𝑓 ″ − 2 𝑧 𝑓 ′ + 4 2 𝑓 ] = ( 1 − 𝑧 2 ) 𝑓 ‴ − 2 𝑧 𝑓 ″ − 2 𝑓 ′ − 2 𝑧 𝑓 ″ + 4 2 𝑓 ′ = 0 d d z [ ( 1 − z 2 ) f ″ − 2 z f ′ + 42 f ] = ( 1 − z 2 ) f ‴ − 2 z f ″ − 2 f ′ − 2 z f ″ + 42 f ′ = 0 ddzleft[(1 - z^2)f'' - 2z f' + 42 fright] = (1 - z^2)f''' - 2z f'' - 2f' - 2z f'' + 42 f' = 0 整理得:
( 1 − 𝑧 2 ) 𝑓 ‴ ( 𝑧 ) − 4 𝑧 𝑓 ″ ( 𝑧 ) + 4 0 𝑓 ′ ( 𝑧 ) = 0 ( 1 − z 2 ) f ‴ ( z ) − 4 z f ″ ( z ) + 40 f ′ ( z ) = 0 (1 - z^2)f'''(z) - 4z f''(z) + 40 f'(z) = 0 代入 𝑧 = 0 z = 0 z = 0 :
( 1 − 0 ) 𝑓 ‴ ( 0 ) − 0 + 4 0 𝑓 ′ ( 0 ) = 0 ⟹ 𝑓 ‴ ( 0 ) = − 4 0 𝑓 ′ ( 0 ) = − 4 0 × 3 = − 1 2 0 ( 1 − 0 ) f ‴ ( 0 ) − 0 + 40 f ′ ( 0 ) = 0 ⟹ f ‴ ( 0 ) = − 40 f ′ ( 0 ) = − 40 × 3 = − 120 (1 - 0)f'''(0) - 0 + 40 f'(0) = 0 implies f'''(0) = -40 f'(0) = -40 * 3 = -120 回代得最终积分值:
𝐼 = 𝜋 𝑖 3 ( − 1 2 0 ) = − 4 0 𝜋 𝑖 I = π i 3 ( − 120 ) = − 40 π i I = pi i3(-120) = -40pi i
5. 围道设计实战:菲涅尔(Fresnel)积分与威克旋转
5.1 菲涅尔积分(Fresnel Integrals) 扇形围道严格求解 在物理光学与衍射理论中,菲涅尔积分定义为:
𝐶 ≡ ∫ ∞ 0 c o s 𝑥 2 d 𝑥 , 𝑆 ≡ ∫ ∞ 0 s i n 𝑥 2 d 𝑥 C ≡ ∫ 0 ∞ cos x 2 d x , S ≡ ∫ 0 ∞ sin x 2 d x C equiv int_0^inf cos x^2 ,dx, quad S equiv int_0^inf sin x^2 ,dx 两者恰好是积分 ∫ ∞ 0 𝑒 𝑖 𝑥 2 d 𝑥 = 𝐶 + 𝑖 𝑆 ∫ 0 ∞ e i x 2 d x = C + i S int_0^inf e^ix^2dx = C + iS 的实部与虚部.
围道设计 构造辅助复函数 𝑓 ( 𝑧 ) = 𝑒 𝑖 𝑧 2 f ( z ) = e i z 2 f(z) = e^iz^2 .由于它在全复平面上无奇点(整函数),设计第一象限内半径为 𝑅 R R 、圆心角为 𝜋 / 4 π / 4 pi/4 的闭合扇形围道 Γ = 𝐶 1 + 𝐶 2 + 𝐶 3 Γ = C 1 + C 2 + C 3 Gamma = C_1 + C_2 + C_3 :
𝐶 1 C 1 C_1 (底边实轴) :从 0 0 0 到 𝑅 R R ,𝑧 = 𝑥 z = x z = x ;𝐶 2 C 2 C_2 (圆弧段) :𝑧 = 𝑅 𝑒 𝑖 𝜃 , 𝜃 ∈ [ 0 , 𝜋 / 4 ] z = R e i θ , θ ∈ [ 0 , π / 4 ] z = R e^itheta, theta in [0, pi/4] ;𝐶 3 C 3 C_3 (斜向返回原点射线) :𝑧 = 𝑟 𝑒 𝑖 𝜋 / 4 , 𝑟 z = r e i π / 4 , r z = r e^ipi/4, r 从 𝑅 R R 变到 0 0 0 .由柯西积分定理:
∮ Γ 𝑒 𝑖 𝑧 2 d 𝑧 = ∫ 𝐶 1 + ∫ 𝐶 2 + ∫ 𝐶 3 = 0 ∮ Γ e i z 2 d z = ∫ C 1 + ∫ C 2 + ∫ C 3 = 0 oint_Gamma e^iz^2dz = int_C_1 + int_C_2 + int_C_3 = 0 各段计算 在 𝐶 1 C 1 C_1 上 :∫ 𝐶 1 𝑒 𝑖 𝑧 2 d 𝑧 = ∫ 𝑅 0 𝑒 𝑖 𝑥 2 d 𝑥 𝑅 → ∞ ←←←←←←←←← → 𝐶 + 𝑖 𝑆 ∫ C 1 e i z 2 d z = ∫ 0 R e i x 2 d x → R → ∞ C + i S int_C_1 e^iz^2dz = int_0^R e^ix^2dx R -> inf C + iS ;
在 𝐶 3 C 3 C_3 上 :𝑧 = 𝑟 𝑒 𝑖 𝜋 / 4 ⟹ 𝑧 2 = 𝑟 2 𝑒 𝑖 𝜋 / 2 = 𝑖 𝑟 2 ⟹ 𝑒 𝑖 𝑧 2 = 𝑒 𝑖 ( 𝑖 𝑟 2 ) = 𝑒 − 𝑟 2 z = r e i π / 4 ⟹ z 2 = r 2 e i π / 2 = i r 2 ⟹ e i z 2 = e i ( i r 2 ) = e − r 2 z = r e^ipi/4 implies z^2 = r^2 e^ipi/2 = i r^2 implies e^iz^2 = e^i(ir^2) = e^-r^2 ! 神奇的几何效果 :复平面的 𝜋 / 4 π / 4 pi/4 旋转把高度振荡的因子 𝑒 𝑖 𝑥 2 e i x 2 e^ix^2 奇迹般地化成了 极速高斯指数衰减因子 𝑒 − 𝑟 2 e − r 2 e^-r^2 !
∫ 𝐶 3 𝑒 𝑖 𝑧 2 d 𝑧 = ∫ 0 𝑅 𝑒 − 𝑟 2 𝑒 𝑖 𝜋 / 4 d 𝑟 = − 𝑒 𝑖 𝜋 / 4 ∫ 𝑅 0 𝑒 − 𝑟 2 d 𝑟 𝑅 → ∞ ←←←←←←←←← → − 𝑒 𝑖 𝜋 / 4 √ 𝜋 2 ∫ C 3 e i z 2 d z = ∫ R 0 e − r 2 e i π / 4 d r = − e i π / 4 ∫ 0 R e − r 2 d r → R → ∞ − e i π / 4 π 2 int_C_3 e^iz^2dz = int_R^0 e^-r^2 e^ipi/4dr = -e^ipi/4int_0^R e^-r^2dr R -> inf -e^ipi/42 在圆弧 𝐶 2 C 2 C_2 上利用约当(Jordan)不等式估计 : 在 𝐶 2 C 2 C_2 上 𝑧 2 = 𝑅 2 ( c o s 2 𝜃 + 𝑖 s i n 2 𝜃 ) z 2 = R 2 ( cos 2 θ + i sin 2 θ ) z^2 = R^2(cos 2theta + isin 2theta) :
| 𝑒 𝑖 𝑧 2 | = | 𝑒 − 𝑅 2 s i n 2 𝜃 | = 𝑒 − 𝑅 2 s i n 2 𝜃 | e i z 2 | = | e − R 2 sin 2 θ | = e − R 2 sin 2 θ |e^iz^2| = |e^-R^2sin 2theta| = e^-R^2sin 2theta 由约当不等式:在 𝛼 ∈ [ 0 , 𝜋 / 2 ] α ∈ [ 0 , π / 2 ] alpha in [0, pi/2] 上 s i n 𝛼 ≥ 2 𝜋 𝛼 sin α ≥ 2 π α sinalpha >= 2/pialpha (正弦弧位于弦上方).令 𝛼 = 2 𝜃 α = 2 θ alpha = 2theta :
∣ ∫ 𝐶 2 𝑒 𝑖 𝑧 2 d 𝑧 ∣ ≤ ∫ 𝜋 / 4 0 𝑒 − 𝑅 2 s i n 2 𝜃 𝑅 d 𝜃 = 𝑅 2 ∫ 𝜋 / 2 0 𝑒 − 𝑅 2 s i n 𝛼 d 𝛼 ≤ 𝑅 2 ∫ 𝜋 / 2 0 𝑒 − 𝑅 2 2 𝛼 𝜋 d 𝛼 = 𝜋 4 𝑅 ( 1 − 𝑒 − 𝑅 2 ) 𝑅 → ∞ ←←←←←←←←← → 0 | ∫ C 2 e i z 2 d z | ≤ ∫ 0 π / 4 e − R 2 sin 2 θ R d θ = R 2 ∫ 0 π / 2 e − R 2 sin α d α ≤ R 2 ∫ 0 π / 2 e − R 2 2 α π d α = π 4 R ( 1 − e − R 2 ) → R → ∞ 0 left|int_C_2 e^iz^2dzright| <= int_0^pi/4 e^-R^2sin 2theta R ,dtheta = R2int_0^pi/2 e^-R^2sinalphadalpha <= R2int_0^pi/2 e^-R^22alpha/pidalpha = 4R(1 - e^-R^2) R -> inf 0 联立求解 代入柯西方程:
( 𝐶 + 𝑖 𝑆 ) + 0 − 𝑒 𝑖 𝜋 / 4 √ 𝜋 2 = 0 ⟹ ∫ ∞ 0 𝑒 𝑖 𝑥 2 d 𝑥 = 𝑒 𝑖 𝜋 / 4 √ 𝜋 2 = ( √ 2 2 + 𝑖 √ 2 2 ) √ 𝜋 2 = √ 2 𝜋 4 ( 1 + 𝑖 ) ( C + i S ) + 0 − e i π / 4 π 2 = 0 ⟹ ∫ 0 ∞ e i x 2 d x = e i π / 4 π 2 = ( 2 2 + i 2 2 ) π 2 = 2 π 4 ( 1 + i ) (C + iS) + 0 - e^ipi/42 = 0 implies int_0^inf e^ix^2dx = e^ipi/42 = left(22 + i22right)2 = 24(1 + i) 分离实部与虚部,得到著名的菲涅尔积分精确解:
∫ ∞ 0 c o s 𝑥 2 d 𝑥 = ∫ ∞ 0 s i n 𝑥 2 d 𝑥 = √ 2 𝜋 4 ∫ 0 ∞ cos x 2 d x = ∫ 0 ∞ sin x 2 d x = 2 π 4 int_0^inf cos x^2 ,dx = int_0^inf sin x^2 ,dx = 24
5.2 威克旋转(Wick Rotation) 与物理高斯积分 上述扇形围道中 𝑧 = 𝑟 𝑒 𝑖 𝜋 / 4 z = r e i π / 4 z = r e^ipi/4 的变量代换,在量子力学与量子统计场论中有着极其宏大的推广——威克旋转(Wick Rotation) .
在闵可夫斯基时空中,振荡的跃迁概率幅包含动力学相位 𝑒 𝑖 𝑆 / ℏ ∼ 𝑒 − 𝑖 𝐻 𝑡 / ℏ e i S / ℏ ∼ e − i H t / ℏ e^i S/hbar ~ e^-i H t / hbar ; 物理学家作解析延拓,将实时间轴转动到虚时间轴:
𝑡 → − 𝑖 𝜏 t → − i τ t -> -i tau 振荡的量子振幅瞬间变成统计力学中的玻尔兹曼热衰减权重 𝑒 − 𝛽 𝐻 e − β H e^-beta H (其中 𝜏 ∼ ℏ 𝛽 = ℏ 𝑘 𝐵 𝑇 τ ∼ ℏ β = ℏ k B T tau ~ hbarbeta = hbar/k_B T ),建立了量子力学路径积分与统计物理配分函数的精确对偶映射.
含源多维高斯积分(实对称正定矩阵 𝐀 A A ):
∫ + ∞ − ∞ … ∫ + ∞ − ∞ d 𝑁 𝑥 e x p ( − 1 2 𝒙 ⋅ 𝐀 ⋅ 𝒙 + 𝑱 ⋅ 𝒙 ) = ( 2 𝜋 ) 𝑁 / 2 √ d e t 𝐀 e x p ( 1 2 𝑱 ⋅ 𝐀 − 1 ⋅ 𝑱 ) ∫ − ∞ + ∞ … ∫ − ∞ + ∞ d N x exp ( − 1 2 x ⋅ A ⋅ x + J ⋅ x ) = ( 2 π ) N / 2 det A exp ( 1 2 J ⋅ A − 1 ⋅ J ) int_-inf^+infdotsint_-inf^+infd^N x , expleft(-12x* A*x + J*xright) = (2pi)^N/2det A expleft(12J* A^-1*Jright) 6. 学习衔接 知识小测与巩固练习 补充说明 课后核心技能自测点
被积函数 1 𝑧 ( 𝑧 − 1 ) 1 z ( z − 1 ) 1/z(z-1) 沿 | 𝑧 | = 2 | z | = 2 |z|=2 逆时针积分的值是多少?沿 | 𝑧 − 1 | = 0 . 5 | z − 1 | = 0.5 |z-1|=0.5 呢? 能否用柯西积分公式快速说明为什么全纯函数的实部 𝑢 ( 𝑥 , 𝑦 ) u ( x , y ) u(x,y) 满足平均值定理? 若 𝑓 ( 𝑧 ) f ( z ) f(z) 在以原点为中心的圆盘内解析,𝑓 ( 0 ) = 1 f ( 0 ) = 1 f(0)=1 ,计算 ∮ | 𝑧 | = 𝑟 𝑓 ( 𝑧 ) 𝑧 d 𝑧 ∮ | z | = r f ( z ) z d z oint_|z|=r f(z)/zdz 的值. ✏️ 本页面状态:审校中 。 欢迎参与讨论或提交补充。 本页面最近更新:2026/10/1 10:49:01 ,更新历史 发现错误?想一起完善? 在 GitHub 上编辑此页! 本页面贡献者:Leafuke , Physics Learning Wiki Team 本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用