路徑積分
討論 ∫Cf(z)dz 定義之前,須先討論 C 之分割(partition)。
設 C 為平面上之一圓滑曲線,其參數式為 z(t),t∈[a,b],函數 f 在 C 上均有定義,即若 f(z)=u(x,y)+iv(x,y) , u(x,y), v(x,y) 對 (x,y)∈C 均有定義。 在 C 上取 n 個分點 (partition points):
Pn={z0=z(a),z1=z(t1),…,zn−1=z(tn−1),zn=z(b)}, 在 zk−1 到 zk 之一小段曲線 Ck=zk−1zk⌢上任取一點 (作爲代表點,representation point) zk∗=xk∗+iyk∗,即存在有 tk∗∈[tk−1,tk] 使得 zk∗=z(tk∗)∈Ck,而且 C=∪k=1nCk。定義 Sn 為對應之黎曼和 (Riemann sum),即
Sn=k=1∑nf(zk∗)(zk−zk−1)=k=1∑nf(zk∗)Δzk(6.3-1) %20c50e44d9741b4c629f5b75d769735efb/Fig_6_2_6.png)
依據代表點 zk∗ 的選擇,常見的黎曼和分成三種:
zk∗=⎩⎨⎧zk−1,2zk−1+zk,zk,左端點,中點,右端點. 如果當 ∥Pn∥≜1≤k≤nmax∥zk−zk−1∥=1≤k≤nmax∥Δzk∥ 趨近於零時 (這時候,1≤k≤nmax∣tk−tk−1∣=1≤k≤nmax∣Δtk∣→0), Sn 趨近於一個定值,這一個定值稱為 f 沿著曲線 C 之路徑積分(path integral),記作 ∫Cf(z)dz ,此時稱 f 在 C 上可積分。
定義6.3-1. 設 C 為圓滑曲線,則延此曲線之路徑積分(path integral)定義為
∫Cf(z)dz=∥Pn∥→0limi=1∑nf(zk∗)(zk−zk−1), ∥Pn∥=kmax∣Δzk∣.(6.3-2)
範例6.3-1. 以 n=8 的中點黎曼和估算 ∫Cexp(z)dz 之值,其中 C 為連接 A=0 到 B=2+i4π 之線段。
[解]
路徑 C 之參數式為 z(t)=(2+i4π)t, 0≤t≤1。由於
tk=k81(1−0)=8k,zk=z(tk)=(2+i2π)8k=4k+i32kπ=k(41+i32π), 且 Δzk=41+i32π≜Δz,對應的 P8 則為
P8={0,41+i32π,21+i16π,43+i323π,1+i84π,45+i325π,23+i163π,27+i327π,1}. 代表點取中點,因此 zk∗=82k−1+i64(2k−1)π=22k−1Δz,對應的黎曼和則為
S(P8)=k=1∑8f(zk∗)Δz=k=1∑8exp(82k−1+i64(2k−1)π)(41+i32π)≈4.22620+5.20473i. 對應的Python程式如下:
註:正確值為
∫Cexp(z)dz=exp(z)z∈C=exp(z)02+i4π=e2+i4π−1≈4.22485+5.22485i.
除了黎曼和來求路徑積分之近似值外,是否可以算出其真正的值,而非近似值?是,可以透過參數化來達成。從式(6.3-2)來看,
∫Cf(z)dz=∥Pn∥→0limi=1∑nf(zk∗)Δzk=∥Pn∥→0limi=1∑nf(zk∗)ΔtkΔzkΔtk=∫abf(z)z′(t)dt 證明如下。
定理6.3-1. 設 C 為積分路徑,其參數化為 z(t), t∈[a,b],假設 f 在 C 上連續,則
∫Cf(z)dz=∫abf(z(t))z′(t)dt(6.3-3) [證明]
假設 f 在 C 上連續,令 f(z)=u(x,y)+iv(x,y),因 u(x,y) 與 v(x,y) 在路徑 C 上連續,因此可積分(參見微積分的 Line Integral),則 f 在 C 上必可積分。若令 zk=xk+iyk,則 zk−zk−1=(xk−xk−1)+i(yk−yk−1),由式(6.3-1)可得
Sn=k=1∑n[u(xk∗,yk∗)+iv(xk∗,yk∗)][(xk−xk−1)+i(yk−yk−1)]=k=1∑n[u(xk∗,yk∗)(xk−xk−1)−v(xk∗,yk∗)(yk−yk−1)]+ik=1∑n[v(xk∗,yk∗)(xk−xk−1)+u(xk∗,yk∗)(yk−yk−1)]=k=1∑n[u(xk∗,yk∗)Δxk−v(xk∗,yk∗)Δyk]+ik=1∑n[v(xk∗,yk∗)Δxk+u(xk∗,yk∗)Δyk] 故
∫Cf(z)dz=∫C(udx−vdy)+i∫C(vdx+udy) 因為 C 為圓滑曲線
∫C(udx−vdy)=∫ab[u(x(t),y(t))x′(t)−v(x(t),y(t))y′(t)]dt ∫C(vdx+udy)=i∫ab[v(x(t),y(t))x′(t)+u(x(t),y(t))y′(t)]dt 代入上述積分式可得
∫Cf(z)dz=∫ab(ux′−vy′)dt+i∫ab(vx′+uy′)dt 即為所求。 ■
從定理6.3-1的證明過程中,若 f=u+i v 以及 z=x+iy,可知下式(6.3-3)同義於:
∫Cf(z)dz=∫C(udx−vdy)+i∫C(vdx+udy).(6.3-4)
範例6.3-2. 以使用參數化方式計算 ∫Cexp(z)dz 之值,其中 C 為連接 A=0 到 B=2+i4π 之線段。
[解]
路徑 C 之參數式為 z(t)=(2+i4π)t, 0≤t≤1,因此 z′(t)=2+i4π,因此
∫Cexp(z)dz=∫01exp(t(2+i4π))(2+i4π)dt=exp(t(2+i4π))01=exp(2+i4π)−1.
範例6.3-3. 計算 ∫Cz2dz 之值,其中 C 之參數式為 γ(t)=t(1+i), t∈[0,1] 。
[解]
γ′(t)=1+i,因此
∫Cz2dz=∫01t2(1+i)2(1+i)dt=3t301(1+i)3=32(−1+i).
範例6.3-4. 計算 ∫C1(0)zm dz 之值,其中 m∈Z。
[解]
單位圓 C1(0) 之參數式為 z(t)=eit, t∈[−π,π],因此 z′(t)=ieit,因此
∫C1(0)zmdz=∫−ππeimtieitdt=i∫−ππei(m+1)tdt=⎩⎨⎧2πi,im+1ei(m+1)t−ππ=0,m=−1,m=−1.
範例6.3-5. 計算 ∫C1+(2)z−21dz 之值,其中 C1+(2) 表示以 2 為圓心半徑為 1 的上半圓。
[解]
以 2 為圓心的單位上半圓 C1+(2) 之參數式為 z(t)=2+eit, t∈[0,π],因此 z′(t)=ieit,因此
∫C1+(0)z−21dz=∫0πeit1ieitdt=∫0πidt=πi. 註:此例將 C1+(2) 改為整個圓,即 C1(2),則積分值變成 2πi。同時,分母等於零的點 z=2 也在 C1(2) 的內部。
設 dzdF(z)=f(z) 且 延曲線 C 之參數式為 z(t), t∈[a,b],則因
dtdF=dzdFdtdz=F′(z)z′(t)=f(z)z′(t), 可知複數微積分基本定理的另一個版本如下,參見 6.5 節:
∫Cf(z)dz=∫abf(z)z′(t)dt=F(z(t))ab=F(z(b))−F(z(a)).
這類型的積分具有下列性質:
性質:設 k∈C,以及 f 及 g 在 C 均可積分,則有
- ∫−Cf(z)dz=−∫Cf(z)dz.
- ∮C(f(z)±g(z))dz=∮Cf(z)dz±∮Cg(z)dz.
- ∫Ckf(z)dz=k∫Cf(z)dz.
- ∫C1+C2f(z)dz=∫C1f(z)dz+∫C2f(z)dz.
- 設 γ1(t) 與 γ2(t) 其中 t∈[a,b] 為路徑 C 之兩種參數表示法,則
∫Cf(z)dz=∫abf(γ1(t))γ2′(t)dt=∫abf(γ2(t))γ2′(t)dt
[證明]
只證明1的部分,其餘練習。
設
z(t), t∈[a,b] 為路徑 C 之參數式,令 τ=−t 且 γ(τ), τ∈[−b,−a] 為 −C 之參數式,即 γ(τ)=z(−τ),因此
∫−Cf(z)dz=∫−b−af(γ(τ))γ′(τ)dτ=∫−b−af(z(−τ))(−1)z′(−τ)dτ=∫−a−bf(z(−τ))z′(−τ)dτ=t=−τ∫abf(z(t))z′(t)dt=∫Cf(z)dz.
路徑與積分值相關性
範例6.3-6. 證明 ∫C1zdz=∫C2zdz=4+2i ,其中 C1 與 C2 圖示如下:
[解]
C1 之參數式為 z(t)=−1−i+t[(3+i)−(−1−i)]=−1−i+(4+2i)t, t∈[0,1],因此 z′(t)=4+2i,因此
∫C1zdz=∫01[−1−i+(4+2i)t](4+2i)dt=∫01[−(2+6i)+(12+16i)t]dt=−(2+6i)t+(12+16i)2t201=4+2i. C2 係沿著 x=(y+1)2−1 的拋物線,因此其參數式為 z(t)=x(t)+iy(t)=(t2+2t)+it=t2+(2+i)t, t∈[−1,1],因此 z′(t)=2t+2+i,因此
∫C2zdz=∫−11(t2+2t+it)(2t+2+i)dt=∫−11(2t3+6t2+3t)dt+i∫−11(3t2+4t)dt=(2t4+2t3+23t2)+i(t3+2t2)−11=4+2i. 註:令 F(z)=2z2,則 F’(z)=z,因此 ∫z(a)z(b)zdz=2z2z(a)z(b),亦即
∫−1−i3+izdz=2z2−1−i3+i=21[(3+i)2−(−1−i)2]=4+2i
範例6.3-7. 計算 ∫C1zˉdz 與 ∫C2zˉdz,其中 C1 與 C2 圖示如下:
%20c50e44d9741b4c629f5b75d769735efb/Fig_6_3_5.png)
[解]
C1 為從正向上半圓之逆向路徑,即 C1 之參數式為 z(t)=e−it, t∈[−π,0],因此 z′(t)=−ie−it,而 C2 之參數式為
⎩⎨⎧z1(t)=−1+it,z2(t)=−1+2t+i,z3(t)=1+i(1−t),0≤t≤1,0≤t≤1,0≤t≤1, 對應之導數為 z1′(t)=i, z2′(t)=2, z3′(t)=−i。因此
∫C1zˉdz=∫−π0eit(−ie−it)dt=−i∫−π0dt=−πi. 而
∫C2zˉdz=∫01(−1−it)idt+∫01(−1+2t−i)2dt+∫01[1−i(1−t)](−i)dt=∫01[(6t−3)−4i]dt=3t2−3t−4it01=−4i.
從範例6.3-6表示 ∫Czdz 之積分值與路徑無關,只與積分路徑的起點和終點座標有關;然而範例6.3-7顯示 ∫Czˉdz 之積分值則除了與積分路徑的起點與終點座標有關外,也和路徑有關。問題來了,函數具有哪些性質,可以讓 ∫Cf(z)dz 之積分值與路徑無關?
例題6.3-8. 計算 f(z)=z2 及 g(z)=z1 沿著下圖所示從 −1 到 1 之上半圓與下半圓之積分值。
%20c50e44d9741b4c629f5b75d769735efb/Fig_6_5.png)
[解]
設 C1 為從−1 到 1 的順時針繞之上半圓,其參數式為 z(t)=−e−it,t∈[−π,0],而 C2 為從−1 到 1 的逆時針繞之下半圓,其參數式為 z(t)=eit,t∈[−π,0],故
∫C1f(z)dz=∫−π0e−i2tie−itdt=−i∫−π0e−3itdt=−i−3ie−i3t−π0=31(e0−ei3π)=32, 又
∫C2f(z)dz=∫−π0ei2tieitdt=i∫−π0ei3tdt=i3iei3t−π0=31(1−e−i3π)=32. 且
∫C1g(z)dz=∫−π0−e−it1ie−itdt=−i∫−π0dt=−πi, ∫C2g(z)dz=∫−π0eit1ieitdt=i∫−π0dt=πi. 整理得
∫C1f(z)dz=∫C2f(z)dz,∫C1g(z)dz=∫C2g(z)dz. 亦即函數 f 的積分和路徑無關,但函數 g 的積分和路徑有關。
這個例子說明
f 沿著 C1 或 C2 其值均相等,事實上從 −1 到 1 點做任一個簡單圓滑曲線 C, 其值均為 32,而 32=31(13−(−1)3),即
∫Cf(z)dz=∫z=−1z=1z2dz=31z3−11=32, 其值只是起點與終點有關而與所經之路徑無關。但同樣的路徑, g(z) 則有不同之值,即 ∫Cg(z)dz 之值與路徑有關。這兩個函數的最大差異是函數 g 在單位圓內有分母為零的點(極點) z=0。又從此例中可看出單位圓 C1+(0)=−C1+C2,因此
∫C1+(0)f(z)dz=∫−C1f(z)dz+∫C2f(z)dz=−32+32=0,∫C1+(0)g(z)dz=∫C1g(z)dz+∫C2g(z)dz=πi+πi=2πi.
如果仔細檢查,若 F′(z)=f(z),則有
∫Cf(z)dz=F(z)−11=F(1)−F(−1) 為複數之微積分基本定理形式,但對於函數 g(z) 而言, 是因 (logz)′=z1, log1−log(−1)=−2kπ, k∈Z,此值不定,因為−1剛好在函數 logz 的分支切上面。回憶函數 logα(z)=ln∣z∣+iθ,其中 θ∈argαz,即 α<θ≤α+2π,因此若取 logz=log−2πz 則 log−2π1−log−2π(−1)=−πi;但若取 logz=log2πz ,則有 log2π1−log2π(−1)=2πi−πi=πi,所以我們自然會問一個問題:在什麼情況下一個複數積分只與起點、終點有關,而與所走的路徑無關呢?將在下一節討論。
以下介紹兩個重要的不等式。
定理6.3-2. (絕對值不等式)
設 f(t)=u(t)+iv(t) 為 t 之連續函數,則 ∫abf(t)dt≤∫ab∣f(t)∣dt。
[證明]
當 ∫abf(t)dt=0,此式明顯成立。假設 ∫abf(t)dt=r0eiθ0, r0=0,則有
∫abf(t)dt=r0=∫abf(t)eiθ0dt=Re∫abf(t)eiθ0dt=∫abRe(f(t)eiθ0)dt≤∫abf(t)eiθ0dt=∫ab∣f(t)∣dt. 上述最後一個等號用到性質 ∣eiθ0∣=1。 ■
定理6.3-3. (ML不等式)
設 f(z)=u(x,y)+iv(x,y) 在曲線 C 上之連續函數,則 ∫Cf(z)dz≤M L,其中 M 為 ∣f(z)∣ 在曲線 C 之上界,而 L 為曲線 C 之長度。
[證明]
設 z=γ(t), t∈[a,b] 為曲線 C 之參數化,因此
∫Cf(z)dz=∫abf(γ(t))γ′(t)dt≤定理6.3-2∫ab∣f(γ(t))∣∣γ′(t)∣dt≤M∫ab∣γ′(t)∣dt=M∫C∣dz∣=M L. ■
例題6.2-2. 估計積分之上界 ∫Cz2+11dz,其中 C 為從 2 到 2+i 之線段。
%20c50e44d9741b4c629f5b75d769735efb/Fig_6_3_5%201.png)
[解]
首先估計 ∣f(z)∣ 在 z∈C 上之上界,即
z∈Cmaxz2+11 因 ∣z2+1∣=∣z+i∣⋅∣z−i∣, 由圖形可知 ∣z−i∣≥2 (最短距離是從 2+i 到 i)且 ∣z+i∣≥5 (最短距離是從 2 到 −i), 故 z2+11≤251,而且線段 C 的長度等於1,即 L=1,所以由ML不等式知
∫Cf(z)dz≤251L=251⋅1=251.