若被積函數包含有 zα 型式在內,其中 0<α<1,考慮此函數之 log0 分支(即分支切 (branch cut) 在正 x 軸),則若 z=∣z∣eiarg0z=reiθ=0, 0<θ≤2π,因
zα=eαlog0(z)=eα[lnr+iθ]=eαlnreiαθ=eαlnr(cosαθ+sinαθ), 因此 zα 在區域 {reiθ:r>0, 0<θ<2π} 內為解析。
定理8.6-1. 令 P 與 Q 為多項式,且 deg(Q)≥deg(P)+2。設 Q(x)=0, ∀x>0, Q 在原點的零點階數最多為 1,以及 f(z)=Q(z)zαP(z), 0<α<1 ,則
P.V.∫0∞Q(x)xαP(x)dx=1−eiα2π2πij=1∑kRes[f,zj],(8.6-1) 其中 z1,z2,…,zk 為 QP 之非零極點。
[證明]
令正向繞上半圓 CR={Reiθ:0≤θ≤π}、 Cr={reiθ:2π≤θ≤23π}、 Γ1={x−ir):−R≤−x≤0} 以及 Γ2={x+ir:0≤x≤R}。參考下圖,積分路徑選為 C=CR+Γ1−Cr+Γ2 ,選擇 R 使得 C 含 Q 的零點在內,因此
R→∞lim∮Cf(z)dz=2πij=1∑kRes[f,zj]. %20760b71d76c194d598b84b95be174ce65/Fig_8_6_1.png)
從定理8.3-1.的證明過程可知
R→∞lim∫CRf(z)dz=0. 由於 Q 在原點的零點接數最多為 1,亦即存在多項式 g(z) 使得
r→0+lim∫Crf(z)dz=r→0+lim∫Crzα[zRes[QP,0]+g(z)]dz=r→0+lim[irα∫2π23πeiαθ+ir1+α∫2π23πg(reiθ)eiθ]dθ=0. 設 z∈Γ1 ,對應的參數式為
z=x+ir=x2+r2ei(2π−θ), −R≤−x≤0, θ=tan−1xr,
則有
r→0+limzr→0+limzα=r→0+limx−ir=x,=r→0+lim(x−ir)α=r→0+lim(x2+r2)2αeiα(2π−θ)=xα⋅eiα2π. 因此
r→0+lim∫Γ1Q(z)zαP(z)dz=∫R0eiα2πQ(x)xαP(x)dx=−eiα2π∫0RQ(x)xαP(x)dx=−eiα2π∫Γ2Q(x)xαP(x)dx. 可得
r→0+R→∞lim∫Cf(z)dz=R→∞lim∫CRf(z)dz=R→∞lim∫Γ1f(z)dz−r→0+lim∫Crf(z)dz+R→∞lim∫Γ1f(z)dz=R→∞lim(1−eiα2π)∫0RQ(x)xαP(x)dx, 所以
P.V.∫0∞Q(x)xαP(x)dx=1−eiα2π2πij=1∑kRes[f,zj],f(z)=Q(z)zαP(z). ■
範例8.6-1. 計算 ∫0∞x(x+1)xαdx 之值,其中 0<α<1。
[暇積分存在性]
被積函數可以改寫成 x(x+1)xα=x1−α(x−1)1,可得 x1−α(x+1)=x2−α+x1−α≥x2−α≥x1−α, ∀x≥0,且 2−α>1−α>0,因此
∫01x(x+1)xαdx∫1∞x(x+1)xαdx≤∫01x1−α1=r→0+lim∫0rx1−α1dx=r→0+limαxα1b=r→0+limα1−rα=α1,≤∫1∞x2−α1=b→∞lim∫r1x2−α1dx=b→∞limα−1xα−11b=b→∞lim1−α1−b1−α1=1−α1, 整理可得
∫0∞x(x+1)xαdx≤α1+1−α1=α(1−α)1. 即此瑕積分存在。
[解]
由於被積函數分母 x(x+1) 在原點的零極點階數為 1 ,因此其值利用定理8.6-1可得;由於
Res[z(z+1)zα,−1]=z→−1lim(z+1)z(z+1)zα=z→−1limzzα=−1(−1)α=−1eiαπ=−eiαπ 則
P.V.∫0∞x(x+1)xαdx=∫0∞x(x+1)xαdx=1−eiα2π2πiRes[z(z+1)zα,−1]=−1−eiα2π2πieiαπ=2ieiαπ−e−iαππ=sinαππ.
範例8.6-2. 推導 P.V.∫0∞x2+a2lnxdx 之值,其中 a>0。
[解]
令 f(z)=z2+a2log−2πz ,一階極點 z=±ia,以及
log−2πz=ln∣z∣+iarg−2πz=lnr+iθ,z=reiθ=0, −2π<θ≤23π. 此外,
Res[f,ia]=z→ialim(z+ia)z2+a2log−2πz=z→ialim2zlog−2πz+(z+ia)z1=2ialog−2π(ia)=2i1aln∣ia∣+iarg−2πz=2i1alna+i2π=2i1(alna+i2aπ). 設 Cr:z=reiθ 與 CR:z=Reiθ,其中 0≤θ≤π,以及取如下圖所示的積分路徑 C=CR+([−R,−r]×{0})−Cr+([r,R]×{0}),則
%20760b71d76c194d598b84b95be174ce65/Fig_8_6_2.png)
∫Cf(z)dz=2πiRes[f,ia]=2πi2i1(alna+i2aπ)=aπlna+i2aπ2. 又當 R>a 時,透過 ML不等式知
R→∞lim∫CRf(z)dzr→0+lim∫Crf(z)dz=R→∞lim∫0πR2ei2θ+a2lnR+iθiReiθdθ≤R→∞limR2−a2lnR+πRπ=0,=r→0+lim∫0πr2ei2θ+a2lnr+iθireiθdθ≤r→0+lima2−r2π−lnrrπ=0. 所以
∫−R−rf(x)dx=∫−R−rx2+a2ln∣x∣+iπdx=∫rRx2+a2ln∣x∣+iπdx=∫rRx2+a2lnxdx+iπ∫rRx2+a21dx, 因此
r→0+R→∞lim∮Cf(z)dz=R→∞lim∫CRf(z)dz+r→0+R→∞lim∫−R−rf(x)dx−r→0+lim∫Crf(z)dz+r→0+R→∞lim∫rRf(x)dx=r→0+R→∞lim∫rRx2+a22lnx+iπdx=2∫0∞x2+a2lnxdx+iπ∫0∞x2+a21dx. 又
∫0∞x2+a21dx=a1∫0∞1+(ax)21d(ax)=a1tan−1ax0∞=2aπ, 整理可得
P.V.∫0∞x2+a2lnxdx=21(r→0+R→∞lim∮Cf(z)dz−iπ∫0∞x2+a21dx)=21(aπlna+i2aπ2−πi2aπ)=2aπlna.
習題
- 計算 P.V.∫0∞x32(1+x)1 dx。
- 計算 P.V.∫0∞x21(1+x)1 dx。
- 計算 P.V.∫0∞(1+x)2x21 dx。
- 計算 P.V.∫0∞x2+1ln(x2+1) dx。 註:取 f(z)=z2+1log(z+i)。
- 計算 P.V.∫0∞(x2+1)2lnx dx。