孤立奇点
若 f(z) 在 z0 不解析,但在某个 Bδ∗(z0) 内解析,则称 z0 为 f(z) 的孤立奇点,否则称 z9 为 f(z) 的非孤立奇点。
分类
f(z)=∑n=−∞∞cn(z−z0)n
可去奇点
洛朗级数中不含负幂项。
f(z)=∑n=0∞cn(z−z0)n
有 limz→z0f(z)=c0
若定义 f(z0)=c0,则 f(z) 在 Bδ(z0) 解析
zsinz=1−3!1z2+5!1z4−…
若定义 f(0)=1,则 f(z) 为整函数
设 z0∈C 为 f(z) 孤立奇点,则下列三条等价:
- z0 为可去奇点
- limz→z0f(z)=A∈C
- f(z) 在某个 z0 的空心领域内有界
(3)⇒(1):
∣f(z)∣≤M
f(z)=∑n=−∞∞cn(z−z0)n
cn=2πi1∮Cr(z−z0)n+1f(z)dz
∣cn∣≤2π1∮Cr∣z−z0∣n+1∣f(z)∣ds≤2π1∮Crrn+1Mds=rnM
∀n<0,cn=0
极点
洛朗级数中只有有限多个负幂项且最高系数非零的负幂次为 −m,则 z0 称为 f(z) 的 m 级极点.
若 m=1,称为简单极点.
f(z)=(z−z0)m1[c−m+c−m+1(z−z0)+…]=(z−z0)mg(z)
z5sinz=z51(z−3!z3+…)
z=0 是 4 级极点.
sin5zz=zsin5z1
z=0 是 zsin5z 的 4 级零点,因此是 sin5zz 的 4 级极点.
z0∈C 是 f(z) 的孤立奇点,则下列四条等价:
- z0 为 m 级极点
- 在某个 z0 的去心领域内 f(z)=(z−z0)mg(z),其中 g(z) 在该去心领域内解析且不取 0
- limz→z0(z−z0)mf(z)=A∈C∗
- z0 为 f(z)1 的 m 级零点
(3)⇒(4):
令 g(z)=(z−z0)mf(z)
limz→z0g(z)=A∈C∗,则 z0 为 g(z) 可去奇点
g(z0)=A
f(z)1=g(z)(z−z0)m=(z−z0)mϕ(z)
因此 z0 为 f(z)1 的 m 级零点
本性奇点
e1/z=1+z1+2!⋅z21+…
设 z0∈C 为 f(z) 本性奇点,则对任意的 A∈C,都存在 {zn}n=1∞⊂Bδ∗(z0) 使得
- limn→∞zn=z0
- limn→∞f(zn)=A
设 z0∈C 为 f(z) 本性奇点,则对任意 A∈C,至多有一个例外值 A0,都存在 {zn}n=1∞⊂Bδ∗(z0) 使得
e1/z=A=0
⇔LnA=z1
zn=lnA+i⋅2nπ1
奇点 ∞
若 f(z) 在 ∞ 点空心领域 R<∣z∣<+∞ 内解析,则称 ∞ 为 f(z) 的孤立奇点,否则称 ∞ 为 f(z) 的非孤立奇点.
令 ζ=z1,则
f(z)=n=−∞∑∞c−nζn=ϕ(ζ)
∞ 的奇点类别定义为 ζ=0 的奇点类别
设 z0∈C 为 f(z) 的孤立奇点,
- z0 为可去奇点 ⇔limz→z0f(z)=A∈C
- z0 为极点 ⇔limz→z0f(z)=∞
- z0 为本性奇点 ⇔limz→z0f(z) 不存在且不为 ∞
设 f(z)=e1/(z−1)+ez−11−z1
找出 f(z) 在 C 上的所有奇点并分类
- 0,limz→0[ez−11−z1]=−21 可去
- 1,本性奇点
- 2nπi,n=0,ez−1 的 1 阶零点,1 级奇点
- ∞,非孤立奇点
留数
z0∈C 为 f(z) 的孤立奇点
f(z)=n=−∞∑∞cn(z−z0)n
且有
c−1=2πi1∮Cf(z)dz
定义 Res[f(z),z0]=c−1 为 f(z) 在 z0 点留数。
若 z0 为 f(z) 可去奇点,则 Res[f(z),z0]=0
设 f(z) 在单连通域 D 内除优先个孤立奇点 z1,z2,…,zn 处处解析,C⊂D 为一条包含各奇点的正向简单闭曲线,则
∮Cf(z)dz=2πik=1∑nRes[f(z),zk]
由复合闭路定理可以简单证明。
留数计算规则
设 z0∈C 为 f(z) 简单极点,则 Res[f(z),z0]=limz→z0(z−z0)f(z)
设 f(z)=Q(z)P(z),P(z),Q(z) 在 z0∈C 解析,P(z0)=0,Q(z0)=0,Q′(z0)=0(即 z0 为 Q(z) 的简单零点),则 z0 为 f(z) 简单极点,有
Res[f(z),z0]=Q′(z0)P(z0)
设 z0∈C 为 f(z) 的 m(≥1) 级极点,则 ∀k≥m 有
Res[f(z),z0]=(k−1)!1z→z0limdzk−1dk−1[(z−z0)kf(z)]
用洛朗展开可以简单证明.
f(z)=z(z−1)2cos2πz
计算 I=∮∣z∣=2f(z)dz
I=2πi(Res[f(z),0]+Res[f(z),1])
z=0 为简单极点
Res[f(z),0]=z→0limzf(z)=1
z=1 为简单极点
Res[f(z),1]=z→z0lim(z−1)f(z)=−2π
也可用
Res[f(z),1]=z→1lim[(z−1)2f(z)]′=−2π
f(z)=z6z−sinz
求 I=∮∣z∣=2f(z)dz
I=2πiRes[f(z),0]
f(z)=z61[z−(z−3!z3+…)]
直接得到 c−1=-1201
z=0 为 3 级极点
也可用 Rule 3 解决
In=∮∣z∣=2zn−1zdz,n≥2
In=2πin=0∑n−1Res[f(z),zk]
Res[f(z),zk]=(zn−1)′∣z=zkzk=nzk2
In={2πi,n=20,n≥3
∞ 处的留数
若 ∞ 为孤立奇点
f(z)=n=−∞∑∞cnzn,R<∣z∣<∞
定义 ∞ 处的留数为 Res[f(z),∞]=−c−1
设 f(z) 在 Cˉ 上除去有限个孤立奇点(包括 ∞)后解析,则 f(z) 在各个奇点处(包括 ∞)的留数之和为 0
Res[f(z),∞]=Res[−z21f(z1),0]
- ! ∞ 若为可去奇点,不一定有 Res[f(z),∞]=0
f(z)=n=−∞∑∞cnzn,R<∣z∣<∞
令 ζ=z1,
f(z)=n=−∞∑∞c−nζn=ϕ(ζ),0<∣ζ∣<R1
ζ2ϕ(ζ)=n=−∞∑∞c−nζn−2
因此有
c−1=Res[ζ2ϕ(ζ),0]=Res[ζ21f(ζ1),0]=Res[z21f(z1),0]
In=∮∣z∣=2zn−1zdz,n≥2
也可用全留数定理解决:
In=−2πiRes[zn−1z,∞]=2πiRes[z21z−n−1z1,0]=2πiRes[1−znzn−3,0]
f(z)=z−1ze1/z
求 I=∮∣z∣=2f(z)dz
I=2πi(Res[f,1]+Res[f,0])=…
或用全留数定理
I=−2πiRes[f,∞]=2πiRes[z21z1−1z1ez,0]=2πiRes[z2(1−z)ez,0]
f(z)=(z−3)(z+i)2025(z−1)20261
求 I=∮∣z∣=2f(z)dz
I=−2πi(Res[f,3]+Res[f,∞])
Res[f,3] 易得,
Res[f,∞]=−Res[z21(z1−3)(z1+i)2025(z1−1)20261,0]=0
f(z)=Qm(z)Pn(z) 为有理函数,且 m−n≥2,则 Res[f,∞]=0
留数在定积分计算中的应用
三角有理函数
形如 I=∫02πR(cosθ,sinθ)dθ,R 为二元有理函数
令 z=eiθ=cosθ+isinθ
cosθ=2z+z=2zz2+1
sinθ=2iz−z=2izz2−1
dz=ieiθdθ
dθ=izdz
I=∮∣z∣=1R(2zz2+1,2izz2−1)izdz
f(z)=R(2zz2+1,2izz2−1)iz1
要求 f(z) 在 ∣z∣=1 上无奇点,
令 {zk} 为 f(z) 在 ∣z∣<1 内所有极点
则 I=2πi∑kRes[f(z),zk]
I(a)=∫02πa−cosθdθ,a>1
I(a)=∮∣z∣=1a−2zz2+1izdz=−i2∮∣z∣=1z2−2az+1dz=2i⋅2πiRes[z2−2az+11,z1]=−(z2−2az+1)′∣z=z14π=…
I(p)I~(p)J(p)I(p)I~(p)=∫02π1−2pcosθ+p2cosmθdθ=∫02π1−2pcosθ+p2sinmθdθ=I(p)+iI~(p)=∫02π1−2pcosθ+p2eimθdθ=∮∣z∣=11+p2−2p2zz2+1izzmdz=−i1∮∣z∣=1p(z2+1)−z(p2+1)zmdz=i∮∣z∣=1(z−p)(pz−1)zmdz=i⋅2πiRes[(z−p)(pz−1)zm,p]=2π⋅1−p2zm=ReJ(p)=2π⋅1−p2zm=ImJ(p)=0
有理函数
形如 I=∫−∞∞R(x)dx
R(x)=Qm(x)Pn(x),且 m−n≥2,R(z)=Qm(z)Pn(z) 在 R 上无极点.
令 {zk} 为 R(z) 在 C 上所有极点
则 I=2πi∑Imzk>0Res[f(z),zk]
I=∫−∞∞1+x4dx
zk=ei(π+2kπ)/4,k=0,1,2,3
I=2πi(Res[z4+11,z0]+Res[z4+11,z1])
有理函数乘三角函数
形如 I=∫−∞∞R(x)eiaxdx=∫−∞∞R(x)cosaxdx+i∫−∞∞R(x)sinaxdx,a∈R,a=0
其中 R(x)=Qm(x)Pn(x),m−n≥1 且 R(z) 在 R 上无极点
则
I=⎩⎨⎧2πiImzk>0∑Res[R(z)eiaz,zk],a>0−2πiImzk>0∑Res[R(z)eiaz,zk],a<0
I(b)J(b)I(b)=∫−∞∞x2+b2xsinxdx,b>0=ImJ(b)=∫−∞∞x2+b2xeixdx=2πiRes[z2+b2zeiz,bi]=2πi⋅(z2+b2)′∣z=bizeiz=πebii=−πsinb