奇点 Singularity,洛朗级数 Laurent series 和 残数 Residue
复变函数 f ( z ) f(z) f ( z ) 有时会存在一些不可解析点 s 0 s_0 s 0 这些点被称为 奇点 Singularity ,我们可以在这些奇点周围进行 洛朗级数 Laurent series 展开,洛朗级数和泰勒级数类似,但是它比泰勒级数多了所有的负次幂项
f ( z ) = ∑ n = 0 ∞ a n ( z − s 0 ) n ( 1 ) f ( z ) = ∑ n = − ∞ ∞ a n ( z − s 0 ) n ( 2 ) \begin{align*}
f(z) &= \sum_{n=0}^{\infty} a_n(z - s_0)^n &\quad (1) \\
f(z) &= \sum_{n=-\infty}^{\infty} a_n(z - s_0)^n &\quad (2)
\end{align*} f ( z ) f ( z ) = n = 0 ∑ ∞ a n ( z − s 0 ) n = n = − ∞ ∑ ∞ a n ( z − s 0 ) n ( 1 ) ( 2 )
( 1 ) (1) ( 1 ) 是泰勒级数, ( 2 ) (2) ( 2 ) 是洛朗级数
因为包含了负次幂的项,所以洛朗级数很有用,可以用于计算围绕奇点的闭合路径积分
残数 Residue
为什么说包含负次幂的项很有用呢?我们可以分析 ( z − s 0 ) n (z - s_0)^n ( z − s 0 ) n ,我们可以沿着围绕 s 0 s_0 s 0 的闭合曲线 C C C 进行路径积分 (假设我们从 C C C 上的点 z 0 z_0 z 0 出发又回到 z 0 z_0 z 0 )
∮ C ( z − s 0 ) n d z \oint_{C} (z - s_0)^n \ dz ∮ C ( z − s 0 ) n d z
围绕奇点的闭合曲面积分
此时结果会根据 n n n 的值不同而有所不同
如果 n ≠ − 1 n \neq -1 n = − 1 即 n = ⋯ , − 3 , − 2 , − 1 , − 0 , − 1 , − 2 , − 3 , ⋯ n = \cdots,-3,-2,-1,\phantom{-}0,\phantom{-}1,\phantom{-}2,\phantom{-}3,\cdots n = ⋯ , − 3 , − 2 , − 1 , − 0 , − 1 , − 2 , − 3 , ⋯
由柯西-古萨定理可得
∮ C ( z − s 0 ) n d z = ( z − s 0 ) n + 1 n + 1 ∣ z 0 z 0 = 0 \oint_{C} (z - s_0)^n \ dz = \left. \frac{(z-s_0)^{n+1}}{n+1} \right|_{z_0}^{z_0} =0 ∮ C ( z − s 0 ) n d z = n + 1 ( z − s 0 ) n + 1 z 0 z 0 = 0
如果 n = − 1 n = -1 n = − 1
由于复对数函数的性质,我们实际上没有回到 z 0 z_0 z 0 ,而是走到了 z 1 z_1 z 1
∮ C ( z − s 0 ) n d z = log ( z − s 0 ) ∣ z 1 z 0 = 2 π i log ( z ) = log ( ∣ z ∣ ) + i ( θ + 2 π n ) ⏟ v \begin{align*}
\oint_{C} (z - s_0)^n \ dz &= \left. \log{(z - s_0)} \right|_{z_1}^{z_0} = 2\pi i \\
\log(z) &= \log(|z|) + i\underbrace{\left(\theta + 2\pi n \right)}_v
\end{align*} ∮ C ( z − s 0 ) n d z log ( z ) = log ( z − s 0 ) ∣ z 1 z 0 = 2 π i = log ( ∣ z ∣ ) + i v ( θ + 2 π n )
走到了另一个点
要对此进行证明也很容易
z − s 0 = R e i θ ⇒ d z = i R e i θ d θ ∫ θ = 0 2 π ( z − s 0 ) − 1 d z = ∮ C i R e i θ R e i θ d θ = ∫ θ = 0 2 π i d θ = 2 π i \begin{align*}
z - s_0 = Re^i\theta \Rightarrow dz &= iRe^i\theta \ d\theta \\
\int_{\theta = 0}^{2\pi} (z - s_0)^{-1} \ dz &= \oint_{C} \frac{iRe^{i\theta}}{Re^{i\theta}} \ d\theta \\
&= \int_{\theta = 0}^{2\pi} i \ d\theta \\
&= 2\pi i
\end{align*} z − s 0 = R e i θ ⇒ d z ∫ θ = 0 2 π ( z − s 0 ) − 1 d z = i R e i θ d θ = ∮ C R e i θ i R e i θ d θ = ∫ θ = 0 2 π i d θ = 2 π i
柯西积分公式 Cauchy Integral Formula
柯西积分公式描述了复变函数在一个简单闭合曲线内的值与沿该曲线积分之间的关系,只要知道复变函数在曲线上的积分,就可以知道复变函数在曲线围成的区域内任意一点的值
具体来说,如果 f ( z ) f(z) f ( z ) 在闭合曲线 C C C 上以及其内部是解析的,并且点 a a a 是曲线围成区域内的一点,则可以定义 f ( z ) z − a \dfrac{f(z)}{z-a} z − a f ( z ) ,如此一来 a a a 就是 f ( z ) z − a \dfrac{f(z)}{z-a} z − a f ( z ) 的奇异点,有了这些信息,柯西积分公式是
∮ C f ( z ) z − a d z = 2 π i f ( a ) \oint_{C} \frac{f(z)}{z-a} \ dz = 2\pi i \ f(a) ∮ C z − a f ( z ) d z = 2 π i f ( a )
柯西积分公式示意图
如此一来,我们就可以间接的求出 f ( a ) f(a) f ( a ) 的值,从而规避可能存在问题的 f ( a ) f(a) f ( a )
我们可以对柯西积分公式进行不严谨的证明
∮ C f ( z ) z − a d z = ∮ C f ( a ) z − a d z + ∮ C f ( z ) − f ( a ) z − a d z \oint_{C} \frac{f(z)}{z-a} \ dz = \oint_{C} \frac{f(a)}{z-a} \ dz + \oint_{C} \frac{f(z) - f(a)}{z-a} \ dz ∮ C z − a f ( z ) d z = ∮ C z − a f ( a ) d z + ∮ C z − a f ( z ) − f ( a ) d z
我们首先来看 ∮ C f ( a ) z − a d z \oint_{C} \frac{f(a)}{z-a} \ dz ∮ C z − a f ( a ) d z 这一部分,根据柯西-古萨定理,我们有
∮ C f ( a ) z − a d z = f ( a ) ∮ C 1 z − a d z = f ( a ) 2 π i \oint_{C} \frac{f(a)}{z-a} \ dz = f(a) \oint_{C} \frac{1}{z-a} \ dz = f(a) \ 2\pi i ∮ C z − a f ( a ) d z = f ( a ) ∮ C z − a 1 d z = f ( a ) 2 π i
我们再来看 ∮ C f ( z ) − f ( a ) z − a d z \oint_{C} \frac{f(z) - f(a)}{z-a} \ dz ∮ C z − a f ( z ) − f ( a ) d z
由于 f ( z ) f(z) f ( z ) 在 a a a 周围是解析的,所以我们可以找到一个 ϵ > 0 \epsilon > 0 ϵ > 0 使得 ∣ f ( z ) − f ( a ) ∣ < ϵ |f(z) - f(a)| < \epsilon ∣ f ( z ) − f ( a ) ∣ < ϵ 而不管 ϵ \epsilon ϵ 多么的小(但是大于零)都可以成立,再一次根据柯西-古萨定理我们有
∮ C ϵ z − a d z = ϵ 2 π i \oint_{C} \frac{\epsilon}{z-a} \ dz = \epsilon \ 2\pi i ∮ C z − a ϵ d z = ϵ 2 π i
因为我们定义 ∣ f ( z ) − f ( a ) ∣ < ϵ |f(z) - f(a)| < \epsilon ∣ f ( z ) − f ( a ) ∣ < ϵ ,所以显然
∮ C f ( z ) − f ( a ) z − a ≤ ∮ C ϵ z − a d z = ϵ 2 π i \oint_{C} \frac{f(z) - f(a)}{z-a} \leq \oint_{C} \frac{\epsilon}{z-a} \ dz = \epsilon \ 2\pi i ∮ C z − a f ( z ) − f ( a ) ≤ ∮ C z − a ϵ d z = ϵ 2 π i
我们可以使得 ϵ → 0 \epsilon \to 0 ϵ → 0 那么 ∮ C f ( z ) − f ( a ) z − a → 0 \oint_{C} \frac{f(z) - f(a)}{z-a} \to 0 ∮ C z − a f ( z ) − f ( a ) → 0
所以
∮ C f ( z ) z − a d z = ∮ C f ( a ) z − a d z + ∮ C f ( z ) − f ( a ) z − a d z = f ( a ) 2 π i + 0 = f ( a ) 2 π i \begin{align*}
\oint_{C} \frac{f(z)}{z-a} \ dz &= \oint_{C} \frac{f(a)}{z-a} \ dz + \oint_{C} \frac{f(z) - f(a)}{z-a} \ dz \\
&= f(a) \ 2\pi i + 0 \\
&= f(a) \ 2\pi i
\end{align*} ∮ C z − a f ( z ) d z = ∮ C z − a f ( a ) d z + ∮ C z − a f ( z ) − f ( a ) d z = f ( a ) 2 π i + 0 = f ( a ) 2 π i