主題
Search

拉普拉斯方程——球座標系

在 球座標系 中，尺度因子 為 h_r=1, h_theta=rsinphi, h_phi=r，分離函式為 f_1(r)=r^2, f_2(theta)=1, f_3(phi)=sinphi，得到 Stäckel 行列式 為 S=1。

拉普拉斯算符 是

del ^2=1/(r^2)partial/(partialr)(r^2partial/(partialr))+1/(r^2sin^2phi)(partial^2)/(partialtheta^2)+1/(r^2sinphi)partial/(partialphi)(sinphipartial/(partialphi)).
(1)

為了在 球座標系 中求解 拉普拉斯方程，嘗試透過寫入 變數分離法

F(r,theta,phi)=R(r)Theta(theta)Phi(phi).
(2)

那麼 亥姆霍茲微分方程 變為

(d^2R)/(dr^2)PhiTheta+2/r(dR)/(dr)PhiTheta+1/(r^2sin^2phi)(d^2Theta)/(dtheta^2)PhiR+(cosphi)/(r^2sinphi)(dPhi)/(dphi)ThetaR+1/(r^2)(d^2Phi)/(dphi^2)ThetaR=0.
(3)

現在除以 RThetaPhi,

(r^2sin^2phi)/(PhiRTheta)PhiTheta(d^2R)/(dr^2)+2/r(r^2sin^2phi)/(PhiRTheta)PhiTheta(dR)/(dr)+1/(r^2sin^2phi)(r^2sin^2phi)/(PhiRTheta)PhiR(d^2Theta)/(dtheta^2)+(cosphi)/(r^2sinphi)(r^2sin^2phi)/(PhiThetaR)(dPhi)/(dphi)ThetaR+1/(r^2)(r^2sin^2phi)/(PhiRTheta)(d^2Phi)/(dphi^2)ThetaR=0
(4)
((r^2sin^2phi)/R(d^2R)/(dr^2)+(2rsin^2phi)/R(dR)/(dr))+(1/Theta(d^2Theta)/(dtheta^2)) +((cosphisinphi)/Phi(dPhi)/(dphi)+(sin^2phi)/Phi(d^2Phi)/(dphi^2))=0.
(5)

方程 (5) 第二部分的解必須是正弦型的，所以微分方程是

(d^2Theta)/(dtheta^2)1/Theta=-m^2,
(6)

其解可以定義為 復 函式，其中 m=-infty, ..., infty

Theta(theta)=A_me^(imtheta),
(7)

或者定義為 實 數正弦和餘弦函式的和，其中 m=-infty, ..., infty

Theta(theta)=S_msin(mtheta)+C_mcos(mtheta).
(8)

將 (6) 代入 (7),

(r^2)/R(d^2R)/(dr^2)+(2r)/R(dR)/(dr)+1/(sin^2phi)((cosphisinphi)/Phi(dPhi)/(dphi)+(sin^2phi)/Phi(d^2Phi)/(dphi^2)-m^2)=0.
(9)

徑向部分必須等於一個常數

(r^2)/R(d^2R)/(dr^2)+(2r)/R(dR)/(dr)=l(l+1)
(10)
r^2(d^2R)/(dr^2)+2r(dR)/(dr)=l(l+1)R.
(11)

但這是 尤拉微分方程，所以我們嘗試 形如

R=sum_(n=0)^inftya_nr^(n+c).
(12)

那麼

r^2sum_(n=0)^infty(n+c)(n+c-1)a_nr^(n+c-2)+2rsum_(n=0)^infty(n+c)a_nr^(n+c-1) -l(l+1)sum_(n=0)^inftya_nr^(n+c)=0
(13)
sum_(n=0)^infty(n+c)(n+c-1)a_nr^(n+c)+2sum_(n=0)^infty(n+c)a_nr^(n+c) -l(l+1)sum_(n=0)^inftya_nr^(n+c)=0
(14)
sum_(n=0)^infty[(n+c)(n+c+1)-l(l+1)]a_nr^(n+c)=0.
(15)

這對於 冪 r 的所有次冪都必須成立。對於 r^c 項（其中 n=0),

c(c+1)=l(l+1),
(16)

這僅當 c=l,-l-1 且所有其他項消失時才成立。所以對於 n!=l, -l-1，a_n=0。因此，R 分量的解由下式給出

R_l(r)=A_lr^l+B_lr^(-l-1).
(17)

將 (17) 代回 (◇),

l(l+1)-(m^2)/(sin^2phi)+(cosphi)/(sinphi)1/Phi(dPhi)/(dphi)+1/Phi(d^2Phi)/(dphi^2)=0
(18)
Phi^('')+(cosphi)/(sinphi)Phi^'+[l(l+1)-(m^2)/(sin^2phi)]Phi=0,
(19)

這是關於 x=cosphi 和 m=0, ..., l 的 連帶勒讓德微分方程。因此，一般 復 數解是

sum_(l=0)^inftysum_(m=-l)^l(A_lr^l+B_lr^(-l-1))P_l^m(cosphi)e^(-imtheta) =sum_(l=0)^inftysum_(m=-1)^l(A_lr^l+B_lr^(-l-1))Y_l^m(theta,phi),
(20)

其中

Y_l^m(theta,phi)=P_l^m(cosphi)e^(-imtheta)
(21)

是（復 數）球諧函式。一般 實 數解是

sum_(l=0)^inftysum_(m=0)^l(A_lr^l+B_lr^(-l-1))P_l^m(cosphi)[S_msin(mtheta)+C_mcos(mtheta)].
(22)

P_l^m 的一些歸一化常數可以被 S_m 和 C_m 吸收，所以這個方程可能以以下形式出現

sum_(l=0)^inftysum_(m=0)^l(A_lr^l+B_lr^(-l-1))P_l^m(cosphi)[S_l^msin(mtheta)+C_l^mcos(mtheta)] =sum_(l=0)^inftysum_(m=0)^l(A_lr^l+B_lr^(-l-1))×[S_l^mY_l^(m(o))(theta,phi)+C_l^mY_l^(m(e))(theta,phi)],
(23)

其中

Y_l^(m(o))(theta,phi)=P_l^m(cosphi)sin(mtheta)
(24)
Y_l^(m(e))(theta,phi)=P_l^m(cosphi)cos(mtheta)
(25)

是 偶 和 奇 （實數）球諧函式。如果存在方位角對稱性，則 Theta(theta) 是常數，並且 Phi 分量的解是 勒讓德多項式 P_l(cosphi)。那麼，一般解是

F(r,phi)=sum_(l=0)^infty(A_lr^l+B_lr^(-l-1))P_l(cosphi).
(26)

另請參閱

亥姆霍茲微分方程——球座標系

使用 探索

參考文獻

Byerly, W. E. An Elementary Treatise on Fourier's Series, and Spherical, Cylindrical, and Ellipsoidal Harmonics, with Applications to Problems in Mathematical Physics. New York: Dover, p. 244, 1959.Moon, P. and Spencer, D. E. Field Theory Handbook, Including Coordinate Systems, Differential Equations, and Their Solutions, 2nd ed. New York: Springer-Verlag, p. 27, 1988.Morse, P. M. and Feshbach, H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, p. 514 and 658, 1953.

請引用本文獻為

Weisstein, Eric W. “拉普拉斯方程——球座標系。” 來自 —— 資源。 https://mathworld.tw/LaplacesEquationSphericalCoordinates.html

主題分類