拉普拉斯方程，又名調和方程、位勢方程，是一種偏微分方程。因為由法國數學家皮埃爾 西蒙拉普拉斯首先提出而得名。求解拉普拉斯方程是電磁學、天文學、熱力學和流體力學等領域經常遇到的一類重要的數學問題，因為這種方程以勢函式的形式描寫電場、引力場和流場等物理物件（一般統稱為 保守場 或 有勢場）的性質。 定義 三維情況下，拉普拉斯方程可由下面的形式描述，問題歸結為求解對實自變數 x、y、z 二階可微的實函式 φ： 使用笛卡爾座標， Delta f frac partial2 fpartial x2 frac partial2 fpartial y2 frac partial2 fpartial z2 0。 使用柱座標， Delta ffrac 1r frac partialpartial r left r frac partial fpartial r right frac 1r2 frac partial2 fpartial phi2 frac partial2 fpartial z2 0 使用球面座標， Delta f frac 1rho2frac partialpartial rho left rho2 frac partial fpartial rhoright frac 1rho2 sintheta frac partialpartial theta left sintheta frac partial fpartial thetaright frac 1rho2 sin2theta frac partial2 fpartial varphi2 0。 使用曲線座標， Delta f frac partialpartial xijleft frac partial fpartial xik gkiright frac partial fpartial xij gjmGamman_mn 0 或 Delta f frac 1sqrt g frac partialpartial xiileft sqrt g gij frac partial fpartial xijright 0 qquad gmathrm detg_ij。 這組方程又經常寫為 nabla2 varphi 0 或者 operatorname divoperatorname gradvarphi 0 ； 其中，div 表示向量場的散度（結果是一個標量場），grad 表示標量場的梯度（結果是一個向量場）。 這方程又可寫為 Delta varphi 0； 其中，Δ 稱為拉普拉斯運算元。 拉普拉斯方程的解稱為調和函式。 如果等號右邊是一個給定的函式 f x y z，即 Delta varphi f， 則該方程稱為泊松方程。拉普拉斯方程和泊松方程是最簡單的橢圓型微分方程。偏微分運算元 nabla2 或 Delta（可以在任意維空間中定義這樣的運算元）稱為拉普拉斯運算元。 邊界條件 拉普拉斯方程的狄利克雷問題可歸結為求解在區域 D 內定義的函式 varphi，使得 varphi 在 D 的邊界上等於某給定的函式。為方便敘述，以下采用拉普拉斯運算元應用的其中一個例子 熱傳導問題作為背景進行介紹：固定區域邊界上的溫度（是邊界上各點位置座標的函式），直到區域內部熱傳導使溫度分佈達到穩定，這個溫度分佈場就是相應的狄利克雷問題的解。 拉普拉斯方程的諾伊曼邊界條件不直接給出區域 D 邊界處的溫度函式 φ 本身，而是 φ 沿 D 的邊界法向的導數。從物理的角度看，這種邊界條件給出的是向量場的勢分佈在區域邊界處的已知效果（對熱傳導問題而言，這種效果便是邊界熱流密度）。 拉普拉斯方程的解稱為調和函式，此函式在方程成立的區域內是解析的。任意兩個函式，如果它們都滿足拉普拉斯方程（或任意線性微分方程），這兩個函式的任意線性組合同樣滿足前述方程。這種非常有用的性質稱為疊加原理。可以根據該原理將各種通解線性組合起來，以滿足所有邊界條件。 二維拉普拉斯方程 兩個自變數的拉普拉斯方程具有以下形式： frac partial2psipartial x2 frac partial2psipartial y2 equiv psi_xx psi_yy 0 解析函式 解析函式的實部和虛部均滿足拉普拉斯方程。換言之，若 z x iy，並且 f z u xy iv xy， 那麼 f z 是解析函式的充要條件是 u xy，v xy 可微，且滿足下列柯西 黎曼方程： u_x v_y quad v_x u_y。 上述方程繼續求導就得到 u_yy v_x_y v_y_x u_x_x。 所以 u 滿足拉普拉斯方程。類似的計算可推得 v 同樣滿足拉普拉斯方程。 反之，給定一個由解析函式（或至少在某點及其鄰域內解析的函式）f z 的實部確定的調和函式，若寫成下列形式： f z varphi xy i psi xy， 則等式 psi_x varphi_y quad psi_y varphi_x。 成立就可使得柯西 黎曼方程得到滿足。上述關係無法確定 ψ，只能得到它的微增量表達式： d psi varphi_y dx varphi_x dy。 φ 滿足拉普拉斯方程意味著 ψ 滿足可積條件： psi_xy psi_yx。 所以可以透過一個線積分來定義 ψ。可積條件和斯托克斯定理的滿足說明線積分的結果與積分經過的具體路徑無關，僅由起點和終點決定。於是，我們便透過複變函式方法得到了 φ 和 ψ 這一對拉普拉斯方程的解。這樣的解稱為一對共軛調和函式。這種構造解的方法只在區域性（複變函式 f z） 的解析域內）有效，或者說，建構函式的積分路徑不能圍繞有 f z 的奇點。譬如，在極座標平面 rθ 上定義函式 varphi log r， 那麼相應的解析函式為 f z log z log r itheta。 在這裡需要注意的是，極角 θ 僅在不包含原點的區域內才是單值的。 拉普拉斯方程與解析函式之間的緊密聯絡說明拉普拉斯方程的任何解都無窮階可導（這是解析函式的一個性質），因此可以展開成冪級數形式，至少在不包含奇點的圓域內是如此。這與波動方程的解形成鮮明對照，後者包含任意函式，其中一些的可微分階數是很小的。 冪級數和傅立葉級數之間存在著密切的關係。如果我們將函式 f 在複平面上以原點為中心，R 為半徑的圓域內展開成冪級數，即 f z sum_n0infty c_n zn， 將每一項係數適當地分離出實部和虛部 c_n a_n i b_n。 那麼 f z sum_n0infty left a_n rn cos n theta b_n rn sin n thetaright i sum_n1infty left a_n rn sin ntheta b_n rn cos n thetaright， 這便是 f 的傅立葉級數。這些三角函式自身也可以用倍角公式展開。 流體動力學 設 u、v 分別為滿足定常、不可壓縮和無旋條件的流體速度場的 x 和 y 方向分量（這裡僅考慮二維流場），那麼不可壓縮條件為： u_x v_y0， 無旋條件為： v_x u_y 0。 若定義一個標量函式 psi，使其微分滿足： d psi v dx u dy， 那麼不可壓縮條件便是上述微分式的可積條件。積分的結果函式 psi 稱為流函式，因為它在同一條流線上各點的值是相同的。psi 的一階偏導為： psi_x v quad psi_yu， 無旋條件即令 psi 滿足拉普拉斯方程。psi 的共軛調和函式 varphi 稱為速度勢。柯西 黎曼方程要求 varphi_xu quad varphi_yv。 所以每一個解析函式都對應著平面內的一個定常不可壓縮無旋流場。解析函式的實部為速度勢函式，虛部為流函式。 靜電學 根據麥克斯韋方程組，二維空間中不隨時間變化的電場 uv 滿足： nabla times uv v_x u_y 0， 和 nabla cdot uv rho， 其中 ρ 為電荷密度。第一個麥克斯韋方程便是下列微分式的可積條件： d varphi u dx v dy， 所以可以構造電勢函式 φ 使其滿足 varphi_x u quad varphi_y v。 第二個麥克斯韋方程即： varphi_xx varphi_yy rho， 這是一個泊松方程，當空間不包含自由電荷時，方程等號右邊變為 0，方程變為拉普拉斯方程。 三維拉普拉斯方程 基本解 拉普拉斯方程的基本解滿足 nabla cdot nabla u u_xx u_yy u_zz delta xxyyzz， 其中的三維 δ 函式代表位於 x y z 的一個點源。 由基本解的定義，若對 u 作用拉普拉斯運算元，再把結果在包含點源的任意體積內積分，那麼 iiint_V nabla cdot nabla u dV 1。 由於座標軸旋轉不改變拉普拉斯方程的形式，所以基本解必然包含在那些僅與到點源距離 r 相關的解中。如果我們選取包含點源、半徑為 a 的球形域作為積分域，那麼根據高斯散度定理 1 iiint_V nabla cdot nabla u dV iint_S u_r dS 4pi a2 u_r a。 求得在以點源為中心，半徑為 r 的球面上有 u_r r frac 14pi r2， 所以 u frac 14pi r。 經過類似的推導同樣可求得二維形式的解 u frac ln r2pi。 格林函式 格林函式是一種不但滿足前述基本解的定義，而且在體積域 V 的邊界 S 上還滿足一定的邊界條件的基本解。譬如， G xyzxyz 可以滿足 nabla cdot nabla G delta xxyyzz quad hbox in quad V， G 0 quad hbox if quad xyz quad hbox on quad S。 現設 u 為在 V 內滿足泊松方程的任意解： nabla cdot nabla u f， 且 u 在邊界 S 上取值為 g，那麼我們可以應用格林定理（是高斯散度定理的一個推論），得到 iiint_V left G nabla cdot nabla u u nabla cdot nabla G right dV iiint_V nabla cdot left G nabla u u nabla G right dV iint_S left G u_n u G_n right dS。 un 和 Gn 分別代表兩個函式在邊界 S 上的 法向 導數。考慮到 u 和 G 滿足的條件，可將這滿足狄利克雷邊界條件的公式化簡為 u xyz iiint_V G f dV iint_S G_n g dS。 所以格林函式描述了量 f 和 g 對 xyz 點函式值的影響。 圓球殼案例 格林函式在半徑為 a 的球面內的點上得值可以透過映象法求得：距球心 ρ 的源點 P 的透過球面的 反射映象 P 距球心 rho frac a2rho。 需要注意的是，如果 P 在球內，那麼 P 將在球外。於是可得格林函式為 G frac 14 pi R frac a4 pi rho R； 其中，R 表示距源點 P 的距離，R 表示距映象點 P 的距離。從格林函式上面的表示式可以推出泊松積分公式。設 ρ、θ 和 φ 為源點 P 的三個球座標分量。此處 θ 按照物理學界的通用標準定義為座標矢徑與豎直軸（z 軸）的夾角（與歐洲習慣相同，與美國習慣不同）。於是球面內拉普拉斯方程的解為： u P frac 14pi a3left 1 frac rho2a2 right iint frac g thetavarphi sin theta dtheta dvarphia2 rho2 2 a rho cos Theta32 ； 其中， cos Theta cos theta cos theta sintheta sinthetacos theta theta。 這個公式的一個顯見的結論是：若 u 是調和函式，那麼 u 在球心處的取值為其在球面上取值的平均。於是我們可以立即得出以下結論：任意一個調和函式（只要不是常函式）的最大值必然不會在其定義域的內部點取得。 參見 球面調和函式 勢流理論 電勢 參考文獻 嚴鎮軍編，《數學物理方程》，第二版，中國科學技術大學出版社，合肥，2002，ISBN 9787312007996O177 LC Evans Partial Differential Equations American Mathematical Society Providence 1998 ISBN 9780821807729 I G Petrovsky Partial Differential Equations W B Saunders Co Philadelphia 1967 A D Polyanin Handbook of Linear Partial Differential Equations for Engineers and Scientists Chapman HallCRC Press Boca Raton 2002 ISBN 9781584882992 A Sommerfeld Partial Differential Equations in Physics Academic Press New York 1949 Pijush KKundu Fluid Mechanics Academic Press 2002 外部連結 拉普拉斯方程（特解和邊值問題） 在 EqWorld 網站上的介紹（專門介紹數學方程的網站，英文）。 初始值問題舉例 exampleproblemscom 網站上使用拉普拉斯方程解題的例子（英文）。 Category 皮埃爾 西蒙拉普拉斯 Category 偏微分方程 Category 橢圓型偏微分方程 Category 調和函式