文章目录
- 一、复变函数项级数
-
- 二、解析函数的孤立奇点
- 三、留数与留数定理
- 四、留数在计算实积分中的应用
-
- 1. 形如 ∫ 0 2 π R ( cos θ , sin θ ) d θ \int_0^{2\pi}R(\cos\theta,\sin\theta)\mathrm{d}\theta ∫02πR(cosθ,sinθ)dθ的积分
- 2. 形如 ∫ − ∞ + ∞ R ( x ) d x \int_{-\infty}^{+\infty} R(x)\mathrm{d}x ∫−∞+∞R(x)dx的积分
- 3. 形如 ∫ − ∞ + ∞ R ( x ) e a i x d x \int_{-\infty}^{+\infty} R(x)e^{aix}\mathrm{d}x ∫−∞+∞R(x)eaixdx( a > 0 a>0 a>0)的积分
一、复变函数项级数
泰勒级数
定理 设 f ( z ) f(z) f(z)为区域 D D D内的解析函数, z 0 ∈ D z_0\in D z0∈D, d d d为 z 0 z_0 z0到 D D D的边界上各点的最短距离,则当 ∣ z − z 0 ∣ < d |z-z_0|∣z−z0∣<d时, f ( z ) f(z) f(z)可唯一地展开为幂级数 f ( z ) = ∑ n = 0 + ∞ c n ( z − z 0 ) n f(z)=\sum\limits_{n=0}^{+\infty}c_n{(z-z_0)}^n f(z)=n=0∑+∞cn(z−z0)n其中 c n = 1 2 π i ∮ ∣ ζ − z 0 ∣ = ρ f ( ζ ) ( ζ − z 0 ) n + 1 d ζ = f ( n ) ( z 0 ) n ! c_n=\frac{1}{2\pi i}\oint_{|\zeta-z_0|=\rho}\frac{f(\zeta)}{{(\zeta-z_0)}^{n+1}}\mathrm{d}\zeta=\frac{f^{(n)}(z_0)}{n!} cn=2πi1∮∣ζ−z0∣=ρ(ζ−z0)n+1f(ζ)dζ=n!f(n)(z0)这里 ∣ ζ − z 0 ∣ = ρ |\zeta-z_0|=\rho ∣ζ−z0∣=ρ是围绕 z 0 z_0 z0的一个环域,它完全被包含在解析的区域内。
如何求出收敛环域的半径呢?
推论 设函数 f ( z ) f(z) f(z)在区域 D D D内解析, z 0 ∈ D z_0\in D z0∈D, α \alpha α为 D D D内距 z 0 z_0 z0最近的一个奇点,则使泰勒展开式成立的圆域的半径 R = ∣ α − z 0 ∣ R=|\alpha-z_0| R=∣α−z0∣。即:在 z 0 z_0 z0处的泰勒展开式的收敛半径为最近奇点到 z 0 z_0 z0的距离。
常用泰勒展开式:
- e z = 1 + z + z 2 2 ! + ⋯ + z n n ! + ⋯ ( ∣ z ∣ < + ∞ ) e^z=1+z+\frac{z^2}{2!}+\cdots+\frac{z^n}{n!}+\cdots\quad(|z|<+\infty) ez=1+z+2!z2+⋯+n!zn+⋯(∣z∣<+∞)
- sin z = e i z − e − i z 2 i = z − z 3 3 ! + z 5 5 ! − ⋯ + ( − 1 ) n z 2 n + 1 ( 2 n + 1 ) ! + ⋯ ( ∣ z ∣ < + ∞ ) \sin z=\frac{e^{iz}-e^{-iz}}{2i}=z-\frac{z^3}{3!}+\frac{z^5}{5!}-\cdots+{(-1)}^n\frac{z^{2n+1}}{(2n+1)!}+\cdots\quad(|z|<+\infty) sinz=2ieiz−e−iz=z−3!z3+5!z5−⋯+(−1)n(2n+1)!z2n+1+⋯(∣z∣<+∞)
- cos z = e i z + e − i z 2 = 1 − z 2 2 ! + z 4 4 ! − ⋯ + ( − 1 ) n z 2 n ( 2 n ) ! + ⋯ ( ∣ z ∣ < + ∞ ) \cos z=\frac{e^{iz}+e^{-iz}}{2}=1-\frac{z^2}{2!}+\frac{z^4}{4!}-\cdots+{(-1)}^n\frac{z^{2n}}{(2n)!}+\cdots\quad(|z|<+\infty) cosz=2eiz+e−iz=1−2!z2+4!z4−⋯+(−1)n(2n)!z2n+⋯(∣z∣<+∞)
- ln ( 1 + z ) = z − z 2 2 + z 3 3 − z 4 4 + ⋯ + ( − 1 ) n z n + 1 n + 1 + ⋯ ( ∣ z ∣ < 1 ) \ln(1+z)=z-\frac{z^2}{2}+\frac{z^3}{3}-\frac{z^4}{4}+\cdots+{(-1)}^n\frac{z^{n+1}}{n+1}+\cdots\quad(|z|<1) ln(1+z)=z−2z2+3z3−4z4+⋯+(−1)nn+1zn+1+⋯(∣z∣<1)
- ln ( 1 − z ) = − z − z 2 2 − z 3 3 − z 4 4 − ⋯ − z n + 1 n + 1 − ⋯ ( ∣ z ∣ < 1 ) \ln(1-z)=-z-\frac{z^2}{2}-\frac{z^3}{3}-\frac{z^4}{4}-\cdots-\frac{z^{n+1}}{n+1}-\cdots\quad(|z|<1) ln(1−z)=−z−2z2−3z3−4z4−⋯−n+1zn+1−⋯(∣z∣<1)
- ( 1 + z ) α = 1 + α z + α ( α − 1 ) 2 ! z 2 + α ( α − 1 ) ( α − 2 ) 3 ! z 3 + ⋯ + α ( α − 1 ) ⋯ ( α − n + 1 ) n ! z n + ⋯ ( ∣ z ∣ < 1 ) {(1+z)}^{\alpha}=1+\alpha z+\frac{\alpha(\alpha-1)}{2!}z^2+\frac{\alpha(\alpha-1)(\alpha-2)}{3!}z^3+\cdots+\frac{\alpha(\alpha-1)\cdots(\alpha-n+1)}{n!}z^n+\cdots\quad(|z|<1) (1+z)α=1+αz+2!α(α−1)z2+3!α(α−1)(α−2)z3+⋯+n!α(α−1)⋯(α−n+1)zn+⋯(∣z∣<1)
善用裂项;遇到平方要积分。
洛朗级数
定理 设 f ( z ) f(z) f(z)在圆环域 R 1 < ∣ z − z 0 ∣ < R 2 R_1<|z-z_0|R1<∣z−z0∣<R2内解析( R 1 R_1 R1和 R 2 R_2 R2常根据奇点确定),则在此圆环域内 f ( z ) f(z) f(z)必能唯一地展开成双边幂级数 f ( z ) = ∑ n = − ∞ + ∞ c n ( z − z 0 ) n f(z)=\sum\limits_{n=-\infty}^{+\infty}c_n{(z-z_0)}^n f(z)=n=−∞∑+∞cn(z−z0)n其中 c n = 1 2 π i ∮ C f ( z ) ( z − z 0 ) n + 1 d z ( n = 0 , ± 1 , ± 2 , ⋯ ) c_n=\frac{1}{2\pi i}\oint_C\frac{f(z)}{{(z-z_0)}^{n+1}}\mathrm{d}z\quad(n=0,\pm1,\pm2,\cdots) cn=2πi1∮C(z−z0)n+1f(z)dz(n=0,±1,±2,⋯) C C C为该圆环域内绕 z 0 z_0 z0的任一正向简单闭曲线。
常用的收敛域: 0 < ∣ z − z 0 ∣ < R 0<|z-z_0|0<∣z−z0∣<R、 R < ∣ z − z 0 ∣ < + ∞ R<|z-z_0|<+\infty R<∣z−z0∣<+∞、 0 < ∣ z − z 0 ∣ < + ∞ 0<|z-z_0|<+\infty 0<∣z−z0∣<+∞
正幂项部分称为解析部分,用于确定收敛域上界( R 2 R_2 R2);
负幂项部分称为主要部分,用于确定收敛域下界( R 1 R_1 R1)。
计算方法/注意事项:
- 善用裂项:比如 1 ( z − 1 ) ( z − 2 ) = 1 z − 2 − 1 z − 1 \frac{1}{(z-1)(z-2)}=\frac{1}{z-2}-\frac{1}{z-1} (z−1)(z−2)1=z−21−z−11;
- 注意收敛半径:比如 1 z − 1 \frac{1}{z-1} z−11,在 1 < ∣ z ∣ < 2 1<|z|<2 1<∣z∣<2的情况下, ∣ z − 1 ∣ ≤ ∣ z ∣ + ∣ 1 ∣ = 2 |z-1|\le|z|+|1|=2 ∣z−1∣≤∣z∣+∣1∣=2,不一定有 ∣ z − 1 ∣ < 1 |z-1|<1 ∣z−1∣<1,故不能直接按 1 1 − z \frac{1}{1-z} 1−z1的公式展开,而 1 z < 1 \frac{1}{z}<1 z1<1,故将其处理成 1 z 1 1 − 1 z \frac{1}{z}\frac{1}{1-\frac{1}{z}} z11−z11,展开为 1 z ∑ n = 0 ∞ 1 z n \frac{1}{z}\sum\limits_{n=0}^{\infty}\frac{1}{z^n} z1n=0∑∞zn1;在 0 < ∣ z ∣ < 1 0<|z|<1 0<∣z∣<1的情况下,就可以直接展开了, 1 z − 1 = − ∑ n = 0 ∞ z n \frac{1}{z-1}=-\sum\limits_{n=0}^{\infty}z^n z−11=−n=0∑∞zn。
- 注意积分:例如求 1 z 2 \frac{1}{z^2} z21在 z 0 = 1 z_0=1 z0=1处、在 0 < ∣ z − 1 ∣ < 1 0<|z-1|<1 0<∣z−1∣<1条件下的洛朗展开式,方法为: 1 z 2 = 1 [ ( z − 1 ) + 1 ] 2 \frac{1}{z^2}=\frac{1}{{[(z-1)+1]}^2} z21=[(z−1)+1]21,令 u = z − 1 u=z-1 u=z−1,则 0 < u < 1 00<u<1,转化为 1 ( u + 1 ) 2 \frac{1}{(u+1)^2} (u+1)21。对于 1 ( u + h ) 2 \frac{1}{(u+h)^2} (u+h)21的形式( h h h为常数, ∣ u ∣ < ∣ h ∣ |u|<|h| ∣u∣<∣h∣),处理方式如下: − 1 u + h = − 1 h 1 1 + u h = − 1 h ∑ n = 0 ∞ ( − 1 ) n u n h n -\frac{1}{u+h}=-\frac{1}{h}\frac{1}{1+\frac{u}{h}}=-\frac{1}{h}\sum\limits_{n=0}^{\infty}{(-1)}^n\frac{u^n}{h^n} −u+h1=−h11+hu1=−h1n=0∑∞(−1)nhnun两边求导得 1 ( u + h ) 2 = − 1 h ∑ n = 0 ∞ ( − 1 ) n n u n − 1 h n = − 1 h ∑ n = 1 ∞ ( − 1 ) n n u n − 1 h n \frac{1}{{(u+h)}^2}=-\frac{1}{h}\sum\limits_{n=0}^{\infty}{(-1)}^nn\frac{u^{n-1}}{h^n}=-\frac{1}{h}\sum\limits_{n=1}^{\infty}{(-1)}^n n\frac{u^{n-1}}{h^n} (u+h)21=−h1n=0∑∞(−1)nnhnun−1=−h1n=1∑∞(−1)nnhnun−1因此 1 ( u + 1 ) 2 = − ∑ n = 1 ∞ ( − 1 ) n n u n − 1 \frac{1}{{(u+1)}^2}=-\sum\limits_{n=1}^{\infty}{(-1)}^n n u^{n-1} (u+1)21=−n=1∑∞(−1)nnun−1,即 1 z 2 = ∑ n = 1 ∞ ( − 1 ) n + 1 n ( z − 1 ) n − 1 \frac{1}{z^2}=\sum\limits_{n=1}^{\infty}{(-1)}^{n+1}n {(z-1)}^{n-1} z21=n=1∑∞(−1)n+1n(z−1)n−1。
- 级数可以相乘:对于 e 1 z 1 + z 2 \frac{e^{\frac{1}{z}}}{1+z^2} 1+z2ez1这种式子,没法直接展开,但是可以用 e 1 z e^{\frac{1}{z}} ez1和 1 1 + z 2 \frac{1}{1+z^2} 1+z21的展开式相乘,从而求出前面有限项的系数。在算留数的时候有用。
- 待定系数法:对于分母比较复杂的式子,要求其洛朗展开式的某一项,可以用待定系数法。比如 f ( z ) = 1 z 2 sin z f(z)=\frac{1}{z^2\sin z} f(z)=z2sinz1令 g ( z ) = z sin z g(z)=\frac{z}{\sin z} g(z)=sinzz则 f ( z ) = 1 z 3 g ( z ) f(z)=\frac{1}{z^3}g(z) f(z)=z31g(z)其中 g ( z ) = z z − z 3 3 ! + z 5 5 ! − ⋯ = 1 1 − z 2 3 ! + z 4 5 ! − ⋯ g(z)=\frac{z}{z-\frac{z^3}{3!}+\frac{z^5}{5!}-\cdots}=\frac{1}{1-\frac{z^2}{3!}+\frac{z^4}{5!}-\cdots} g(z)=z−3!z3+5!z5−⋯z=1−3!z2+5!z4−⋯1令 h ( z ) = 1 − z 2 3 ! + z 4 5 ! − ⋯ h(z)=1-\frac{z^2}{3!}+\frac{z^4}{5!}-\cdots h(z)=1−3!z2+5!z4−⋯,则 g ( z ) h ( z ) = 1 g(z)h(z)=1 g(z)h(z)=1下面考虑 g ( z ) g(z) g(z)的洛朗展开式。显然 g ( z ) g(z) g(z)是解析函数(因为分母是解析函数且分母不为0)和偶函数,因此可令 g ( z ) g(z) g(z)的洛朗展开式(其实也是泰勒展开式)为 g ( z ) = c 0 + c 2 z 2 + c 4 z 4 + ⋯ g(z)=c_0+c_2z^2+c_4z^4+\cdots g(z)=c0+c2z2+c4z4+⋯。因此 ( c 0 + c 2 z 2 + c 4 z 4 + ⋯ ) ( 1 − z 2 3 ! + z 4 5 ! − ⋯ ) = 1 (c_0+c_2z^2+c_4z^4+\cdots)\left(1-\frac{z^2}{3!}+\frac{z^4}{5!}-\cdots\right)=1 (c0+c2z2+c4z4+⋯)(1−3!z2+5!z4−⋯)=1比较左边和右边。右边只有常数项的系数非零,其他次项的系数均为 0 0 0。由常数项为 1 1 1得 c 0 ⋅ 1 = 1 ⟹ c 0 = 1 c_0\cdot 1=1\implies c_0=1 c0⋅1=1⟹c0=1由二次项系数为 0 0 0得 c 0 ⋅ 1 + c 2 ⋅ ( − 1 3 ! ) ⟹ c 2 = 1 6 c_0\cdot 1+c_2\cdot\left(-\frac{1}{3!}\right)\implies c_2=\frac{1}{6} c0⋅1+c2⋅(−3!1)⟹c2=61以此类推可以求得 g ( z ) g(z) g(z)展开式的各项系数。因此 g ( z ) = 1 + 1 6 z 2 + ⋯ g(z)=1+\frac{1}{6}z^2+\cdots g(z)=1+61z2+⋯从而有 f ( z ) = 1 z 3 g ( z ) = 1 z 3 + 1 6 1 z + ⋯ f(z)=\frac{1}{z^3}g(z)=\frac{1}{z^3}+\frac{1}{6}\frac{1}{z}+\cdots f(z)=z31g(z)=z31+61z1+⋯由此还可以求得 Res [ f ( z ) , 0 ] = 1 6 \operatorname{Res}[f(z),0]=\frac{1}{6} Res[f(z),0]=61。
二、解析函数的孤立奇点
孤立奇点: f ( z ) f(z) f(z)在 z 0 z_0 z0处解析但在其某一去心邻域内解析,那么称 z 0 z_0 z0是 f ( z ) f(z) f(z)的孤立奇点。
注意:当存在一个趋于 z 0 z_0 z0的奇点序列时, z 0 z_0 z0就是不是孤立奇点了。比如 f ( z ) = 1 sin π z f(z)=\frac{1}{\sin\frac{\pi}{z}} f(z)=sinzπ1的奇点为 z = 1 k ( k = ± 1 , ± 2 , ⋯ ) z=\frac{1}{k}(k=\pm 1,\pm 2,\cdots) z=k1(k=±1,±2,⋯),有一个趋于 0 0 0的奇点序列 1 , 1 2 , 1 3 , ⋯ 1,\frac{1}{2},\frac{1}{3},\cdots 1,21,31,⋯,这使得 0 0 0的任意邻域内都有无穷多个奇点,因此 0 0 0不是 f ( z ) f(z) f(z)的孤立奇点。
孤立奇点分为三类:(按照洛朗级数负幂项的个数分)
- 可去奇点: f ( z ) f(z) f(z)在 z 0 z_0 z0处的洛朗展开式无负幂项,此时洛朗级数就是一个幂级数。补充定义 f ( z 0 ) = lim z → z 0 f ( z ) f(z_0)=\lim\limits_{z\to z_0}f(z) f(z0)=z→z0limf(z)后,函数就变得解析了,此时洛朗级数就是泰勒级数。例如 z = 0 z=0 z=0是 f ( z ) = e z − 1 z f(z)=\frac{e^z-1}{z} f(z)=zez−1的可去奇点。
- 极点:洛朗级数只有有限个负幂项。如果洛朗级数幂次最低的负幂项是 c − m ( z − z 0 ) − m c_{-m}{(z-z_0)}^{-m} c−m(z−z0)−m,则称 z 0 z_0 z0为 f ( z ) f(z) f(z)的 m m m级极点。此时 f ( z ) = c − m ( z − z 0 ) − m + c − m + 1 ( z − z 0 ) − m + 1 + ⋯ + c 0 + c 1 ( z − z 0 ) + ⋯ f(z)=c_{-m}{(z-z_0)}^{-m}+c_{-m+1}{(z-z_0)}^{-m+1}+\cdots+c_0+c_1(z-z_0)+\cdots f(z)=c−m(z−z0)−m+c−m+1(z−z0)−m+1+⋯+c0+c1(z−z0)+⋯令 g ( z ) = ( z − z 0 ) m f ( z ) g(z)={(z-z_0)}^mf(z) g(z)=(z−z0)mf(z),则 g ( z ) g(z) g(z)的洛朗展开式没有负幂项,那么 g ( z ) g(z) g(z)在 z 0 z_0 z0的邻域内解析,且 g ( z 0 ) ≠ 0 g(z_0)\ne 0 g(z0)=0(因为 c − m ≠ 0 c_{-m}\ne 0 c−m=0)。这其实也是 m m m级极点的充要条件。例如 f ( z ) = sin z z 3 ( z + 1 ) f(z)=\frac{\sin z}{z^3(z+1)} f(z)=z3(z+1)sinz的奇点为:二级极点 z 0 = 0 z_0=0 z0=0,一级极点 z 1 = − 1 z_1=-1 z1=−1。
- 本性奇点:洛朗展开式有无穷多个负幂项。例如 f ( z ) = e 1 z f(z)=e^{\frac{1}{z}} f(z)=ez1的本性奇点是 z 0 = 0 z_0=0 z0=0。
魏尔斯特拉斯定理 设函数在 0 < ∣ z − z 0 ∣ < δ 0<|z-z_0|<\delta 0<∣z−z0∣<δ内解析,若 z 0 z_0 z0是 f ( z ) f(z) f(z)的本性奇点,则 ∀ a ∈ C \forall a\in\mathbb{C} ∀a∈C,不论它是有限数还是无穷,都有一个收敛于 z 0 z_0 z0的复数列 { z n } \{z_n\} {zn},使得 lim z → z 0 f ( z n ) = a \lim\limits_{z\to z_0}f(z_n)=a z→z0limf(zn)=a。
形象的(不严谨的)理解:若 z 0 z_0 z0是 f ( z ) f(z) f(z)的本性奇点,则 lim z → z 0 f ( z ) \lim\limits_{z\to z_0}f(z) z→z0limf(z)可以等于任何数。这说明这个极限不存在也不为无穷。
用极限判断孤立奇点的类型(充要条件):若 z 0 z_0 z0是 f ( z ) f(z) f(z)的孤立奇点,则
- lim z → z 0 f ( z n ) \lim\limits_{z\to z_0}f(z_n) z→z0limf(zn)存在且为有限复数 ⟺ \iff ⟺ z 0 z_0 z0是 f ( z ) f(z) f(z)的可去奇点;
- lim z → z 0 f ( z n ) = ∞ ⟺ z 0 \lim\limits_{z\to z_0}f(z_n)=\infty\iff z_0 z→z0limf(zn)=∞⟺z0是 f ( z ) f(z) f(z)的可去奇点;
- lim z → z 0 f ( z n ) \lim\limits_{z\to z_0}f(z_n) z→z0limf(zn)存在且不为无穷 ⟺ \iff ⟺ z 0 z_0 z0是 f ( z ) f(z) f(z)的本性奇点。
零点:设 f ( z ) f(z) f(z)为一不恒等于 0 0 0的解析函数,若 f ( z ) = ( z − z 0 ) m φ ( z ) f(z)={(z-z_0)^m\varphi(z)} f(z)=(z−z0)mφ(z),其中 φ ( z ) \varphi(z) φ(z)在 z 0 z_0 z0处解析,并且 φ ( z 0 ) ≠ 0 \varphi(z_0)\ne 0 φ(z0)=0,则称 z 0 z_0 z0是 f ( z ) f(z) f(z)的一个 m m m级零点。
用导数确定零点:若 f ( z ) f(z) f(z)在 z 0 z_0 z0处解析,则 z 0 z_0 z0为 f ( z ) f(z) f(z)的 m m m级零点的充要条件是 f ( n ) ( z 0 ) = 0 ( n = 0 , 1 , ⋯ , m − 1 ) f^{(n)}(z_0)=0\quad(n=0,1,\cdots,m-1) f(n)(z0)=0(n=0,1,⋯,m−1)且 f ( m ) ( z 0 ) ≠ 0 f^{(m)}(z_0)\ne 0 f(m)(z0)=0。这个可以用泰勒级数来理解。 f ( z ) f(z) f(z)在 z 0 z_0 z0点泰勒级数的最低次项为 c m ( z − z 0 ) m c_m{(z-z_0)}^m cm(z−z0)m。
若 z 0 z_0 z0是 P ( z ) P(z) P(z)的 n n n级零点,也是 Q ( z ) Q(z) Q(z)的 m m m级零点,则 z 0 z_0 z0是 P ( z ) Q ( z ) P(z)Q(z) P(z)Q(z)的 n + m n+m n+m级零点。(这里可以把 l l l级极点理解为 − l -l −l级零点,公式仍然成立。)
零点与极点的关系: z 0 z_0 z0是解析函数 f ( z ) f(z) f(z)的 m m m级极点,当且仅当 z 0 z_0 z0是 1 f ( z ) \frac{1}{f(z)} f(z)1的 m m m级零点。
复变函数的洛必达法则:设函数 f ( z ) f(z) f(z), g ( z ) g(z) g(z)在 z 0 z_0 z0处解析,且 f ( z 0 ) = g ( z 0 ) = 0 f(z_0)=g(z_0)=0 f(z0)=g(z0)=0。那么, lim z → z 0 f ( z ) g ( z ) = lim z → z 0 f ′ ( z ) g ′ ( z ) \lim\limits_{z\to z_0}\frac{f(z)}{g(z)}=\lim\limits_{z\to z_0}\frac{f'(z)}{g'(z)} z→z0limg(z)f(z)=z→z0limg′(z)f′(z)这是因为,设 f ( z ) = c 1 ( z − z 0 ) + c 2 ( z − z 0 ) 2 + ⋯ f(z)=c_1(z-z_0)+c_2{(z-z_0)}^2+\cdots f(z)=c1(z−z0)+c2(z−z0)2+⋯, g ( z ) = d 1 ( z − z 0 ) + d 2 ( z − z 0 ) 2 + ⋯ g(z)=d_1(z-z_0)+d_2{(z-z_0)}^2+\cdots g(z)=d1(z−z0)+d2(z−z0)2+⋯,则 lim z → z 0 f ( z ) g ( z ) = c 1 d 1 \lim\limits_{z\to z_0}\frac{f(z)}{g(z)}=\frac{c_1}{d_1} z→z0limg(z)f(z)=d1c1(上下同时消去 z − z 0 z-z_0 z−z0即得)。而 f ′ ( z ) = c 1 + 2 c 2 ( z − z 0 ) 2 + ⋯ f'(z)=c_1+2c_2{(z-z_0)}^2+\cdots f′(z)=c1+2c2(z−z0)2+⋯, g ′ ( z ) = d 1 + 2 d 2 ( z − z 0 ) 2 + ⋯ g'(z)=d_1+2d_2{(z-z_0)}^2+\cdots g′(z)=d1+2d2(z−z0)2+⋯,所以 lim z → z 0 f ′ ( z ) g ′ ( z ) = c 1 d 1 \lim\limits_{z\to z_0}\frac{f'(z)}{g'(z)}=\frac{c_1}{d_1} z→z0limg′(z)f′(z)=d1c1因此 lim z → z 0 f ( z ) g ( z ) = lim z → z 0 f ′ ( z ) g ′ ( z ) \lim\limits_{z\to z_0}\frac{f(z)}{g(z)}=\lim\limits_{z\to z_0}\frac{f'(z)}{g'(z)} z→z0limg(z)f(z)=z→z0limg′(z)f′(z)成立。
无穷远点作为奇点:函数 f ( z ) f(z) f(z)在无穷远点没有定义,故 ∞ \infty ∞一定是 f ( z ) f(z) f(z)的奇点。若 f ( z ) f(z) f(z)在 U ˚ ( ∞ ) = { z ∣ 0 < R < ∣ z ∣ < + ∞ } \mathring{U}(\infty)=\{z|0U˚(∞)={z∣0<R<∣z∣<+∞}内解析,则称 ∞ \infty ∞为 f ( z ) f(z) f(z)的孤立奇点。研究无穷远点方法为令 t = 1 z t=\frac{1}{z} t=z1,则 f ( z ) = f ( 1 t ) f(z)=f\left(\frac{1}{t}\right) f(z)=f(t1)。令 φ ( t ) = ( 1 t ) \varphi(t)=\left(\frac{1}{t}\right) φ(t)=(t1),则该变换把 ∞ \infty ∞点映到 0 0 0点。此时 lim z → ∞ f ( z ) = lim t → 0 φ ( t ) \lim\limits_{z\to\infty}f(z)=\lim\limits_{t\to 0}\varphi(t) z→∞limf(z)=t→0limφ(t),除非这两个极限都不存在。定义:若 t = 0 t=0 t=0是 φ ( t ) \varphi(t) φ(t)的可去奇点、 m m m级极点或本性奇点,则 z = ∞ z=\infty z=∞也分别是 f ( z ) f(z) f(z)的可去奇点、 m m m级极点或本性奇点,他们分别对应 f ( z ) f(z) f(z)在 U ˚ ( ∞ ) \mathring{U}(\infty) U˚(∞)内的洛朗级数不含 z z z的正幂项、含有有限多个 z z z的正幂项且 z m z^m zm是最高次幂、含有无穷多个 z z z的正幂项的情况。
当 lim z → ∞ f ( z ) \lim\limits_{z\to\infty}f(z) z→∞limf(z)存在(即 ∞ \infty ∞是可去奇点)时,可以补充定义 f ( z ) = lim z → ∞ f ( z ) f(z)=\lim\limits_{z\to\infty}f(z) f(z)=z→∞limf(z),就可以认为 f ( z ) f(z) f(z)在 ∞ \infty ∞是解析的。如 f ( z ) = z z + 1 f(z)=\frac{z}{z+1} f(z)=z+1z。
三、留数与留数定理
函数在有限点处的留数:设 z 0 z_0 z0为 f ( z ) f(z) f(z)的有限孤立奇点,设 f ( z ) f(z) f(z)在 U ˚ ( z 0 ) \mathring{U}(z_0) U˚(z0)内洛朗展开式中 ( z − z 0 ) − 1 {(z-z_0)}^{-1} (z−z0)−1系数为 c − 1 c_{-1} c−1,定义 f ( z ) f(z) f(z)在 z 0 z_0 z0处的留数为 Res [ f ( z ) , z 0 ] = c − 1 = 1 2 π i ∮ C f ( z ) d z \operatorname{Res}[f(z),z_0]=c_{-1}=\frac{1}{2\pi i}\oint_C f(z)\mathrm{d}z Res[f(z),z0]=c−1=2πi1∮Cf(z)dz其中 C C C是围绕 z 0 z_0 z0的正向环路。
函数在无穷远点处的留数:设 ∞ \infty ∞为 f ( z ) f(z) f(z)的孤立奇点, C C C为圆环域 0 < R < ∣ z ∣ < + ∞ 00<R<∣z∣<+∞( f ( z ) f(z) f(z)在其中解析)内绕原点的任一正向简单闭曲线,则定义 f ( z ) f(z) f(z)在 ∞ \infty ∞处的留数为 Res [ f ( z ) , ∞ ] = − c − 1 = 1 2 π i ∮ C − f ( z ) d z \operatorname{Res}[f(z),\infty]=-c_{-1}=\frac{1}{2\pi i}\oint_{C^-} f(z)\mathrm{d}z Res[f(z),∞]=−c−1=2πi1∮C−f(z)dz其中 C − C^- C−是 C C C的负向,即绕 ∞ \infty ∞的正向, c − 1 c_{-1} c−1是 f ( z ) f(z) f(z)在 U ˚ ( ∞ ) \mathring{U}(\infty) U˚(∞)内洛朗展开式中 z − 1 z^{-1} z−1项的系数。
定理 设 f ( z ) f(z) f(z)在扩充复平面内包括无穷远点在内只有有限个孤立奇点,则 f ( z ) f(z) f(z)在各奇点处的留数之和为 0 0 0。
可以这么理解:假设路径 C C C包含了 f ( z ) f(z) f(z)的所有有限奇点,则 1 2 π i ∮ C f ( z ) d z \frac{1}{2\pi i}\oint_C f(z)\mathrm{d}z 2πi1∮Cf(z)dz就是 f ( z ) f(z) f(z)在各有限奇点处的留数之和。现在 Res [ f ( z ) , ∞ ] \operatorname{Res}[f(z),\infty] Res[f(z),∞]定义为它的相反数,那么我们自然知道,在无穷远点处的留数与在各有限点处的留数之和为 0 0 0。利用这个性质可以用来反向计算环路积分的值。
留数定理 设函数 f ( z ) f(z) f(z)在区域 D D D内除有限个有限孤立奇点 z 1 , z 2 , ⋯ , z l z_1,z_2,\cdots,z_l z1,z2,⋯,zl外处处解析, C C C是 D D D内包含所有奇点的任一正向简单闭曲线,则 ∮ C f ( z ) d z = 2 π i ∑ k = 1 l Res [ f ( z ) , z k ] \oint_C f(z)\mathrm{d}z=2\pi i\sum\limits_{k=1}^l \operatorname{Res}[f(z),z_k] ∮Cf(z)dz=2πik=1∑lRes[f(z),zk]
留数的计算:
- 当奇点为可去奇点——留数为 0 0 0
- 当奇点为本性奇点——只能求洛朗展开式的 − 1 -1 −1次项系数
- 当奇点为 m m m级极点:
- 若 m = 1 m=1 m=1,则 Res [ f ( z ) , z 0 ] = lim z → z 0 ( z − z 0 ) f ( z ) \operatorname{Res}[f(z),z_0]=\lim\limits_{z\to z_0}(z-z_0)f(z) Res[f(z),z0]=z→z0lim(z−z0)f(z)
- 若 m > 1 m>1 m>1,则 Res [ f ( z ) , z 0 ] = 1 ( m − 1 ) ! lim z → z 0 d m − 1 d z m − 1 { ( z − z 0 ) m f ( z ) } \operatorname{Res}[f(z),z_0]=\frac{1}{(m-1)!}\lim\limits_{z\to z_0}\frac{\mathrm{d}^{m-1}}{\mathrm{d}z^{m-1}}\{{(z-z_0)}^mf(z)\} Res[f(z),z0]=(m−1)!1z→z0limdzm−1dm−1{(z−z0)mf(z)}注意阶乘和导数阶数都是 m − 1 m-1 m−1,只有 z − z 0 z-z_0 z−z0的次数为 m m m!
理解方法:设 f ( z ) = c − m ( z − z 0 ) − m + ⋯ + c − 1 ( z − z 0 ) − 1 + c 0 + ⋯ f(z)=c_{-m}{(z-z_0)}^{-m}+\cdots+c_{-1}{(z-z_0)}^{-1}+c_0+\cdots f(z)=c−m(z−z0)−m+⋯+c−1(z−z0)−1+c0+⋯,则 ( z − z 0 ) m f ( z ) = c − m + ⋯ + c − 1 ( z − z 0 ) m − 1 + c 0 ( z − z 0 ) m + ⋯ {(z-z_0)}^m f(z)=c_{-m}+\cdots+c_{-1}{(z-z_0)}^{m-1}+c_0{(z-z_0)}^m+\cdots (z−z0)mf(z)=c−m+⋯+c−1(z−z0)m−1+c0(z−z0)m+⋯。我们知道, d n d z n ( z − z 0 ) n = n ! \frac{\mathrm{d}^n}{\mathrm{d}z^n}{(z-z_0)}^n=n! dzndn(z−z0)n=n!,所以 d m − 1 d z m − 1 ( z − z 0 ) m − 1 = ( m − 1 ) ! \frac{\mathrm{d}^{m-1}}{\mathrm{d}z^{m-1}}{(z-z_0)}^{m-1}=(m-1)! dzm−1dm−1(z−z0)m−1=(m−1)!,因此给 ( z − z 0 ) m f ( z ) {(z-z_0)}^m f(z) (z−z0)mf(z)求 m − 1 m-1 m−1阶导后,比 m − 1 m-1 m−1次数低的项变为 0 0 0了,比 m − 1 m-1 m−1次数高的项在 z → z 0 z\to z_0 z→z0时也变成 0 0 0了,只留下了 ( m − 1 ) ! c − 1 (m-1)!c_{-1} (m−1)!c−1。因此 c − 1 = 1 ( m − 1 ) ! lim z → z 0 d m − 1 d z m − 1 { ( z − z 0 ) m f ( z ) } c_{-1}=\frac{1}{(m-1)!}\lim\limits_{z\to z_0}\frac{\mathrm{d}^{m-1}}{\mathrm{d}z^{m-1}}\{{(z-z_0)}^mf(z)\} c−1=(m−1)!1z→z0limdzm−1dm−1{(z−z0)mf(z)}
- 若 m = 1 m=1 m=1,且 f ( z ) f(z) f(z)可以表示为 f ( z ) = φ ( z ) ψ ( z ) f(z)=\frac{\varphi(z)}{\psi(z)} f(z)=ψ(z)φ(z),其中 φ ( z ) \varphi(z) φ(z)、 ψ ( z ) \psi(z) ψ(z)在 z 0 z_0 z0处解析,并且 φ ( z 0 ) ≠ 0 \varphi(z_0)\ne 0 φ(z0)=0, ψ ( z 0 ) = 0 \psi(z_0)=0 ψ(z0)=0, ψ ′ ( z 0 ) ≠ 0 \psi'(z_0)\ne 0 ψ′(z0)=0,则 Res [ f ( z ) , z 0 ] = φ ( z 0 ) ψ ′ ( z 0 ) \operatorname{Res}[f(z),z_0]=\frac{\varphi(z_0)}{\psi'(z_0)} Res[f(z),z0]=ψ′(z0)φ(z0)理解方法:注意到 z 0 z_0 z0是 ψ ( z ) \psi(z) ψ(z)的一级极点。设 φ ( z ) = h 0 + h 1 ( z − z 0 ) + ⋯ \varphi(z)=h_0+h_1(z-z_0)+\cdots φ(z)=h0+h1(z−z0)+⋯, ψ ( z ) = d 1 ( z − z 0 ) + d 2 ( z − z 0 ) 2 + ⋯ \psi(z)=d_1(z-z_0)+d_2{(z-z_0)}^2+\cdots ψ(z)=d1(z−z0)+d2(z−z0)2+⋯, f ( z ) = c − 1 ( z − z 0 ) − 1 + c 0 + ⋯ f(z)=c_{-1}{(z-z_0)}^{-1}+c_0+\cdots f(z)=c−1(z−z0)−1+c0+⋯,则 f ( z ) ψ ( z ) = φ ( z ) f(z)\psi(z)=\varphi(z) f(z)ψ(z)=φ(z)。这里 f ( z ) f(z) f(z)最低次项为 − 1 -1 −1次,是因为 f ( z ) ψ ( z ) f(z)\psi(z) f(z)ψ(z)最低次项为 0 0 0次,而若 f ( z ) f(z) f(z)有更低次项,则该项与 ψ ( z ) \psi(z) ψ(z)最低次项相乘得到的是负幂次项,但 φ ( z ) \varphi(z) φ(z)中并没有负幂次项。同时我们可以得出 d 1 c − 1 = h 0 d_1 c_{-1}=h_0 d1c−1=h0,即 c − 1 = h 0 d 1 = φ ( z 0 ) ψ ′ ( z 0 ) c_{-1}=\frac{h_0}{d_1}=\frac{\varphi(z_0)}{\psi'(z_0)} c−1=d1h0=ψ′(z0)φ(z0)。
- 当奇点为无穷远点: Res [ f ( z ) , ∞ ] = − Res [ f ( 1 z ) 1 z 2 , 0 ] \operatorname{Res}[f(z),\infty]=-\operatorname{Res}\left[f\left(\frac{1}{z}\right)\frac{1}{z^2},0\right] Res[f(z),∞]=−Res[f(z1)z21,0]换言之,设 z 1 , z 2 , ⋯ , z l z_1,z_2,\cdots,z_l z1,z2,⋯,zl是 f ( z ) f(z) f(z)的所有有限孤立奇点,则 ∑ k = 1 l Res [ f ( z ) , z k ] = Res [ f ( 1 z ) 1 z 2 , 0 ] \sum\limits_{k=1}^l\operatorname{Res}[f(z),z_k]=\operatorname{Res}\left[f\left(\frac{1}{z}\right)\frac{1}{z^2},0\right] k=1∑lRes[f(z),zk]=Res[f(z1)z21,0]推导过程:设环路 C : ∣ z ∣ = ρ C:|z|=\rho C:∣z∣=ρ包含所有有限奇点,作变换 t = 1 z t=\frac{1}{z} t=z1,则 C − C^- C−被映射到正向闭曲线 L : ∣ t ∣ = 1 ρ L:|t|=\frac{1}{\rho} L:∣t∣=ρ1。于是 Res [ f ( z ) , ∞ ] = 1 2 π i ∮ C − f ( z ) d z = 1 2 π i ∮ L f ( 1 t ) ( − 1 t 2 ) d t = − Res [ f ( 1 t ) 1 t 2 , 0 ] \operatorname{Res}[f(z),\infty]=\frac{1}{2\pi i}\oint_{C^-}f(z)\mathrm{d}z=\frac{1}{2\pi i}\oint_L f\left(\frac{1}{t}\right)\left(-\frac{1}{t^2}\right)\mathrm{d}t=-\operatorname{Res}\left[f\left(\frac{1}{t}\right)\frac{1}{t^2},0\right] Res[f(z),∞]=2πi1∮C−f(z)dz=2πi1∮Lf(t1)(−t21)dt=−Res[f(t1)t21,0]
技巧:若 f ( z ) f(z) f(z)是偶函数,则 Res [ f ( z ) , 0 ] = 0 \operatorname{Res}[f(z),0]=0 Res[f(z),0]=0。这是因为在 0 0 0处的洛朗展开式没有奇次项,自然有 c − 1 = 0 c_{-1}=0 c−1=0。比如 f ( z ) = z sin 3 z f(z)=\frac{z}{\sin^3 z} f(z)=sin3zz。
四、留数在计算实积分中的应用
1. 形如 ∫ 0 2 π R ( cos θ , sin θ ) d θ \int_0^{2\pi}R(\cos\theta,\sin\theta)\mathrm{d}\theta ∫02πR(cosθ,sinθ)dθ的积分
处理方法:令 z = e i θ z=e^{i\theta} z=eiθ,则 d z = i z d θ \mathrm{d}z=iz\mathrm{d}\theta dz=izdθ, sin θ = z 2 − 1 2 i z \sin\theta=\frac{z^2-1}{2iz} sinθ=2izz2−1, cos θ = z 2 + 1 2 z \cos\theta=\frac{z^2+1}{2z} cosθ=2zz2+1。把积分表达式中的 ∫ 0 2 π \int_0^{2\pi} ∫02π换成 ∮ ∣ z ∣ = 1 \oint_{|z|=1} ∮∣z∣=1、 sin θ \sin\theta sinθ换成 z 2 − 1 2 i z \frac{z^2-1}{2iz} 2izz2−1、 cos θ \cos\theta cosθ换成 z 2 + 1 2 z \frac{z^2+1}{2z} 2zz2+1、 d θ \mathrm{d}\theta dθ换成 d z i z \frac{\mathrm{d}z}{iz} izdz。则 ∫ 0 2 π R ( cos θ , sin θ ) d θ = ∮ ∣ z ∣ = 1 R ( z 2 + 1 2 z , z 2 − 1 2 i z ) d z i z = 2 π i ∑ k = 1 l Res [ 1 i z R ( z 2 + 1 2 z , z 2 − 1 2 i z ) , z k ] = 2 π ∑ k = 1 l Res [ 1 z R ( z 2 + 1 2 z , z 2 − 1 2 i z ) , z k ] \begin{aligned} \int_0^{2\pi}R(\cos\theta,\sin\theta)\mathrm{d}\theta&=\oint_{|z|=1}R\left(\frac{z^2+1}{2z},\frac{z^2-1}{2iz}\right)\frac{\mathrm{d}z}{iz}\\ &=2\pi i\sum\limits_{k=1}^l \operatorname{Res}\left[\frac{1}{iz}R\left(\frac{z^2+1}{2z},\frac{z^2-1}{2iz}\right),z_k\right]\\ &=2\pi\sum\limits_{k=1}^l \operatorname{Res}\left[\frac{1}{z}R\left(\frac{z^2+1}{2z},\frac{z^2-1}{2iz}\right),z_k\right] \end{aligned} ∫02πR(cosθ,sinθ)dθ=∮∣z∣=1R(2zz2+1,2izz2−1)izdz=2πik=1∑lRes[iz1R(2zz2+1,2izz2−1),zk]=2πk=1∑lRes[z1R(2zz2+1,2izz2−1),zk]其中 z 1 , z 2 , ⋯ , z l z_1,z_2,\cdots,z_l z1,z2,⋯,zl是含 z z z的被积函数在 ∣ z ∣ = 1 |z|=1 ∣z∣=1以内的所有奇点。
这里积分上下限不管是 0 0 0到 2 π 2\pi 2π还是 − π -\pi −π到 π \pi π都是一样的,对应的积分环路都是 ∣ z ∣ = 1 |z|=1 ∣z∣=1。对于上下限是 0 0 0到 π \pi π的,需要考虑函数的奇偶性进行处理(一般有 ∫ 0 π = 1 2 ∫ − π π \int_0^\pi=\frac{1}{2}\int_{-\pi}^{\pi} ∫0π=21∫−ππ)。
2. 形如 ∫ − ∞ + ∞ R ( x ) d x \int_{-\infty}^{+\infty} R(x)\mathrm{d}x ∫−∞+∞R(x)dx的积分
要求: R ( x ) = P ( x ) Q ( x ) R(x)=\frac{P(x)}{Q(x)} R(x)=Q(x)P(x)是 x x x的有理函数, P ( x ) P(x) P(x)与 Q ( x ) Q(x) Q(x)为互质多项式, Q ( x ) Q(x) Q(x)的次数至少比 P ( x ) P(x) P(x)高两次
处理方法:找到 R ( z ) R(z) R(z)在上半平面的孤立奇点 z 1 , z 2 , ⋯ , z l z_1,z_2,\cdots,z_l z1,z2,⋯,zl(注意如果在上半平面有一个奇点,那么对称地在下半平面也有一个奇点),则 ∫ − ∞ + ∞ R ( x ) d x = 2 π i ∑ k = 1 l Res [ R ( z ) , z k ] \int_{-\infty}^{+\infty} R(x)\mathrm{d}x=2\pi i\sum\limits_{k=1}^l \operatorname{Res}[R(z),z_k] ∫−∞+∞R(x)dx=2πik=1∑lRes[R(z),zk]特别地,如果 R ( x ) R(x) R(x)的偶函数,则 ∫ 0 + ∞ R ( x ) d x = π i ∑ k = 1 l Res [ R ( z ) , z k ] \int_{0}^{+\infty} R(x)\mathrm{d}x=\pi i\sum\limits_{k=1}^l \operatorname{Res}[R(z),z_k] ∫0+∞R(x)dx=πik=1∑lRes[R(z),zk]
3. 形如 ∫ − ∞ + ∞ R ( x ) e a i x d x \int_{-\infty}^{+\infty} R(x)e^{aix}\mathrm{d}x ∫−∞+∞R(x)eaixdx( a > 0 a>0 a>0)的积分
要求: R ( x ) = P ( x ) Q ( x ) R(x)=\frac{P(x)}{Q(x)} R(x)=Q(x)P(x)是 x x x的有理函数, P ( x ) P(x) P(x)与 Q ( x ) Q(x) Q(x)为互质多项式, Q ( x ) Q(x) Q(x)的次数至少比 P ( x ) P(x) P(x)高一次
处理方法:和第二种情况一样,找到 R ( z ) e a i z R(z)e^{aiz} R(z)eaiz在上半平面的孤立奇点 z 1 , z 2 , ⋯ , z l z_1,z_2,\cdots,z_l z1,z2,⋯,zl(注意如果在上半平面有一个奇点,那么对称地在下半平面也有一个奇点),则 ∫ − ∞ + ∞ R ( x ) e a i x d x = 2 π i ∑ k = 1 l Res [ R ( z ) e a i z , z k ] \int_{-\infty}^{+\infty} R(x)e^{aix}\mathrm{d}x=2\pi i\sum\limits_{k=1}^l \operatorname{Res}[R(z)e^{aiz},z_k] ∫−∞+∞R(x)eaixdx=2πik=1∑lRes[R(z)eaiz,zk]同时,因为 e a i x = cos a x + i sin a x e^{aix}=\cos ax+i\sin ax eaix=cosax+isinax,也可以由此计算 ∫ − ∞ + ∞ R ( x ) cos a x d x = Re ( ∫ − ∞ + ∞ R ( x ) e a i x d x ) ∫ − ∞ + ∞ R ( x ) sin a x d x = Im ( ∫ − ∞ + ∞ R ( x ) e a i x d x ) \int_{-\infty}^{+\infty}R(x)\cos ax\mathrm{d}x=\operatorname{Re}\left(\int_{-\infty}^{+\infty} R(x)e^{aix}\mathrm{d}x\right)\\ \int_{-\infty}^{+\infty}R(x)\sin ax\mathrm{d}x=\operatorname{Im}\left(\int_{-\infty}^{+\infty} R(x)e^{aix}\mathrm{d}x\right) ∫−∞+∞R(x)cosaxdx=Re(∫−∞+∞R(x)eaixdx)∫−∞+∞R(x)sinaxdx=Im(∫−∞+∞R(x)eaixdx)对于 a < 0 a<0 a<0的情况,注意到 e z ‾ = e z ‾ \overline{e^z}=e^{\overline{z}} ez=ez, ∫ − ∞ + ∞ R ( x ) e a i x d x = ∫ − ∞ + ∞ R ( x ) e a i x ‾ d x ‾ = ∫ − ∞ + ∞ R ( x ) e − a i x d x ‾ \int_{-\infty}^{+\infty}R(x)e^{aix}\mathrm{d}x=\overline{\int_{-\infty}^{+\infty}R(x)\overline{e^{aix}}\mathrm{d}x}=\overline{\int_{-\infty}^{+\infty}R(x)e^{-aix}\mathrm{d}x} ∫−∞+∞R(x)eaixdx=∫−∞+∞R(x)eaixdx=∫−∞+∞R(x)e−aixdx这就转化为了 a > 0 a>0 a>0的情况,然后对结果取共轭即可。