UA MATH524 复变函数13 补充:留数计算的例题
UA MATH524 復變函數(shù)13 補充:留數(shù)計算的例題
- 概念與公式總結(jié)
- 例題
- pole的留數(shù)公式
- Laurent級數(shù)法
概念與公式總結(jié)
孤立奇點類型與判斷
假設z0z_0z0?是復變函數(shù)fff的孤立奇點,則
- z0z_0z0?為可去奇點,如果lim?z→z0f(z)\lim_{z \to z_0}f(z)limz→z0??f(z)存在且有限(此時Laurent級數(shù)與冪級數(shù)相同);
- z0z_0z0?為mmm階pole,如果H(z)=f(z)(z?z0)mH(z)=f(z)(z-z_0)^mH(z)=f(z)(z?z0?)m在z0z_0z0?的鄰域解析并且H(z0)≠0H(z_0) \ne 0H(z0?)?=0(此時lim?z→z0f(z)\lim_{z \to z_0}f(z)limz→z0??f(z)存在且為無限、Laurent級數(shù)主部有mmm項)
- z0z_0z0?為本性奇點,如果lim?z→z0f(z)\lim_{z \to z_0}f(z)limz→z0??f(z)不存在(此時Laurent級數(shù)主部有無窮項)
不同類型孤立奇點的留數(shù)計算方法
假設fff在B°(z0,r)\overset{\circ}{B}(z_0,r)B°(z0?,r)上解析,則留數(shù)的定義為
Res(f;z0)=12πi∫∣w?z0∣=sf(w)dw,0<s<rRes(f;z_0)=\frac{1}{2 \pi i} \int_{|w-z_0|=s}f(w)dw,0<s<rRes(f;z0?)=2πi1?∫∣w?z0?∣=s?f(w)dw,0<s<r
- 當fff在z0z_0z0?處解析或者fff在z0z_0z0?處滿足lim?z→z0∣f(z)∣<∞\lim_{z \to z_0}|f(z)|<\inftylimz→z0??∣f(z)∣<∞時,Res(f;z0)=f(z0)Res(f;z_0)=f(z_0)Res(f;z0?)=f(z0?)
- 當z0z_0z0?是fff的mmm階pole時,Res(f;z0)=H(m?1)(z0)(m?1)!,H(z)=f(z)(z?z0)mRes(f;z_0)=\frac{H^{(m-1)}(z_0)}{(m-1)!},H(z)=f(z)(z-z_0)^mRes(f;z0?)=(m?1)!H(m?1)(z0?)?,H(z)=f(z)(z?z0?)m
- 當z0z_0z0?是fff的本性奇點時,Res(f;z0)=c?1Res(f;z_0)=c_{-1}Res(f;z0?)=c?1?其中c?1c_{-1}c?1?是fff在z0z_0z0?處的Laurent級數(shù)的(z?z0)?1(z-z_0)^{-1}(z?z0?)?1項的系數(shù)
例題
pole的留數(shù)公式
例1 f(z)=z2+3z?1z+2f(z)=\frac{z^2+3z-1}{z+2}f(z)=z+2z2+3z?1?
這個函數(shù)的一階pole為z0=?2z_0=-2z0?=?2,令H(z)=(z+2)f(z)=z2+3z?1H(z)=(z+2)f(z)=z^2+3z-1H(z)=(z+2)f(z)=z2+3z?1,
Res(f;?2)=H(1?1)(?2)(1?1)!=?3Res(f;-2)=\frac{H^{(1-1)}(-2)}{(1-1)!}=-3Res(f;?2)=(1?1)!H(1?1)(?2)?=?3
例2 f(z)=ez(z?1)3f(z)=\frac{e^z}{(z-1)^3}f(z)=(z?1)3ez?
這個函數(shù)的三階pole為z0=1z_0=1z0?=1,令H(z)=(z?1)3f(z)=ezH(z)=(z-1)^3f(z)=e^zH(z)=(z?1)3f(z)=ez,
Res(f;1)=H(3?1)(1)(3?1)!=e2Res(f;1)=\frac{H^{(3-1)}(1)}{(3-1)!}=\frac{e}{2}Res(f;1)=(3?1)!H(3?1)(1)?=2e?
例3
f(z)=Log(z)(1+z2)2=Log(z)(z?i)2(z+i)2f(z)=\frac{Log(z)}{(1+z^2)^2}=\frac{Log(z)}{(z-i)^2(z+i)^2}f(z)=(1+z2)2Log(z)?=(z?i)2(z+i)2Log(z)?
所以z1=i,z2=?iz_1=i,z_2=-iz1?=i,z2?=?i是這個函數(shù)的兩個二階pole,記H1(z)=Log(z)(z+i)2,H2(z)=Log(z)(z?i)2H_1(z)=\frac{Log(z)}{(z+i)^2},H_2(z)=\frac{Log(z)}{(z-i)^2}H1?(z)=(z+i)2Log(z)?,H2?(z)=(z?i)2Log(z)?,
Res(f;i)=H1′(i)1!=z+iz?2Log(z)(z+i)3∣z=i=2?iπ(2i)3=π8+i4Res(f;?i)=H2′(?i)1!=z?iz?2Log(z)(z?i)3∣z=?i=2+iπ?(2i)3=π8+i4Res(f;i)=\frac{H_1'(i)}{1!}=\left. \frac{\frac{z+i}{z}-2Log(z)}{(z+i)^3} \right|_{z=i}=\frac{2-i\pi}{(2i)^3}=\frac{\pi}{8}+\frac{i}{4} \\ Res(f;-i)=\frac{H_2'(-i)}{1!}=\left. \frac{\frac{z-i}{z}-2Log(z)}{(z-i)^3} \right|_{z=-i}=\frac{2+i\pi}{-(2i)^3}=\frac{\pi}{8}+\frac{i}{4}Res(f;i)=1!H1′?(i)?=(z+i)3zz+i??2Log(z)?∣∣∣∣∣?z=i?=(2i)32?iπ?=8π?+4i?Res(f;?i)=1!H2′?(?i)?=(z?i)3zz?i??2Log(z)?∣∣∣∣∣?z=?i?=?(2i)32+iπ?=8π?+4i?
例4 計算下面的函數(shù)在z0=0z_0=0z0?=0處的留數(shù)
f(z)=cot?(αz)=cos?(αz)sin?(αz)f(z)=\cot(\alpha z)=\frac{\cos(\alpha z)}{\sin(\alpha z)}f(z)=cot(αz)=sin(αz)cos(αz)?
z0=0z_0=0z0?=0是fff的一階pole,記H(z)=zf(z)=1αcos?(αz)sin?(αz)αzH(z)=zf(z)=\frac{1}{\alpha}\frac{\cos(\alpha z)}{\frac{\sin (\alpha z)}{\alpha z}}H(z)=zf(z)=α1?αzsin(αz)?cos(αz)?
其中lim?z→z0sin?(αz)αz=1\lim_{z \to z_0}\frac{\sin (\alpha z)}{\alpha z}=1limz→z0??αzsin(αz)?=1,所以
Res(f;0)=H(0)=1αRes(f;0)=H(0)=\frac{1}{\alpha}Res(f;0)=H(0)=α1?
例5
f(z)=eiaz?eibzz2,a≠bf(z)=\frac{e^{iaz}-e^{ibz}}{z^2},a \ne bf(z)=z2eiaz?eibz?,a?=b
z0=0z_0=0z0?=0是分母的二階零點,但是lim?z→z0(eiaz?eibz)=0\lim_{z \to z_0}(e^{iaz}-e^{ibz})=0limz→z0??(eiaz?eibz)=0,此時函數(shù)成為0/00/00/0不定型。用洛必達法則考慮
lim?z→z0eiaz?eibzz=lim?z→z0iaeiaz?ibeibz1=i(a?b)\lim_{z \to z_0} \frac{e^{iaz}-e^{ibz}}{z}=\lim_{z \to z_0} \frac{iae^{iaz}-ibe^{ibz}}{1}=i(a-b)z→z0?lim?zeiaz?eibz?=z→z0?lim?1iaeiaz?ibeibz?=i(a?b)
因此z0=0z_0=0z0?=0是fff的一階pole,記H(z)=zf(z)=eiaz?eibzzH(z)=zf(z)=\frac{e^{iaz}-e^{ibz}}{z}H(z)=zf(z)=zeiaz?eibz?,
Res(f;0)=H(0)=i(a?b)Res(f;0)=H(0)=i(a-b)Res(f;0)=H(0)=i(a?b)
Laurent級數(shù)法
例6 求下面的函數(shù)在z0=iz_0=iz0?=i處的留數(shù)
f(z)=eizz(z2+1)2=eizz(z?i)2(z+i)2f(z)=\frac{e^{iz}}{z(z^2+1)^2}=\frac{e^{iz}}{z(z-i)^2(z+i)^2}f(z)=z(z2+1)2eiz?=z(z?i)2(z+i)2eiz?
z0z_0z0?是fff的二階pole,記H(z)=(z?z0)2f(z)=eizz(z+i)2H(z)=(z-z_0)^2f(z)=\frac{e^{iz}}{z(z+i)^2}H(z)=(z?z0?)2f(z)=z(z+i)2eiz?,如果用公式,
Res(f;i)=H′(z0)=(iz2?4z?i)eizz2(z+i)3∣z=z0=?34eRes(f;i)=H'(z_0)=\left.\frac{(iz^2-4z-i)e^{iz}}{z^2(z+i)^3} \right|_{z=z_0}=-\frac{3}{4e}Res(f;i)=H′(z0?)=z2(z+i)3(iz2?4z?i)eiz?∣∣∣∣?z=z0??=?4e3?
就會需要計算一個有一點點難算的導數(shù),所以我們考慮用待定系數(shù)法求fff的Laurent級數(shù)中的c?1c_{-1}c?1?。首先,c?1c_{-1}c?1?是fff的Laurent展開中(z?z0)?1(z-z_0)^{-1}(z?z0?)?1項的系數(shù),那么它就是HHH的展開∑n=0∞an(z?z0)n\sum_{n=0}^{\infty}a_n(z-z_0)^n∑n=0∞?an?(z?z0?)n中(z?z0)(z-z_0)(z?z0?)項的系數(shù)a1a_1a1?,令t=z?z0=z?it=z-z_0=z-it=z?z0?=z?i,則
eiz=ei(t+i)=e?1∑n=0+∞inn!tnz(z+i)2H(z)=(t+i)(t+2i)2∑n=0+∞antne^{iz}=e^{i(t+i)}=e^{-1} \sum_{n=0}^{+\infty} \frac{i^n}{n!}t^n \\ z(z+i)^2H(z)=(t+i)(t+2i)^2 \sum_{n=0}^{+\infty} a_nt^neiz=ei(t+i)=e?1n=0∑+∞?n!in?tnz(z+i)2H(z)=(t+i)(t+2i)2n=0∑+∞?an?tn
根據(jù)ttt項的系數(shù)相等可以得到
e?1i1=i?(2i)2a1+[i?(4i)+(2i)2]a0e^{-1}\frac{i}{1}=i \cdot (2i)^2a_1+[i\cdot (4i)+(2i)^2]a_0e?11i?=i?(2i)2a1?+[i?(4i)+(2i)2]a0?
根據(jù)常數(shù)項相等可以得到
e?1=i?(2i)2a0e^{-1}=i \cdot (2i)^2a_0e?1=i?(2i)2a0?
聯(lián)立可得a1=?34ea_1=-\frac{3}{4e}a1?=?4e3?。
PS. 有一說一,這個問題我寧愿直接求導。。。。
例7 求下面的函數(shù)在z0=0z_0=0z0?=0處的留數(shù)
f(z)=1z2sin?zf(z)=\frac{1}{z^2 \sin z}f(z)=z2sinz1?
z0=0z_0=0z0?=0是fff的三階pole,記
H(z)=z3f(z)=zsin?zH(z)=z^3f(z)=\frac{z}{\sin z}H(z)=z3f(z)=sinzz?
這個也可以直接用公式,需要求HHH的二階導;也可以用待定系數(shù)法,找HHH的展開中z2z^2z2項的系數(shù)。答案為1/61/61/6。
總結(jié)
以上是生活随笔為你收集整理的UA MATH524 复变函数13 补充:留数计算的例题的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA OPTI544 量子光学13 场的
- 下一篇: 贝叶斯多元Logistics回归理论基础