跳轉到

9.3 Laplace 方程與極值原理#

本章導言#

\[ \Delta u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0. \]

Laplace 方程描述一切「達到平衡」的場：穩態溫度、靜電位勢、不可壓無旋流、引力位。它沒有時間——是熱方程（9.1）與波方程（9.2）時間箭頭盡頭的雕像。方程以 Pierre-Simon Laplace 命名（《天體力學》，1785 起），但首先深入研究的是 Laplace 的同代人與論敵 Legendre（5.2 節的球面位勢）與年輕的 Carl Friedrich Gauss（1777–1855）。本節講述調和函數的三大性質——平均值性質、極值原理、解析性——它們是 PDE 理論中最優雅的一章。

平衡的方程：從熱到引力#

熱傳導穩態： \(u_t = 0\) ⇒ \(\Delta u = 0\) （9.1 節時間因子的靜態極限）。 靜電位勢：真空中無電荷處 \(\Delta\phi = 0\) （有電荷處 \(\Delta\phi = -\rho/\varepsilon_0\) ，Poisson 方程 1813）。 流體：不可壓無旋流的速度位勢滿足同一方程——三個物理學科共用一個方程（這種「方程的普適性」是數學物理的常態：1.1 節 \(F=ma\) 、9.1 熱方程亦然）。

邊值問題（Dirichlet 問題）：給定區域 \(\Omega\) 與邊界溫度 \(u|_{\partial\Omega} = f\) ，求內部解。「邊界決定內部」——與 ODE 的「初值決定未來」平行，但性質截然不同（見極值原理）。

平均值性質：調和函數的靈魂#

定理（平均值性質）： \(u\) 在 \(\Omega\) 調和（ \(\Delta u = 0\) ），則 \(u\) 在任意圓心 \(x_0\) 的值等於其圓周上的平均值：

\[ u(x_0) = \frac{1}{2\pi} \int_0^{2\pi} u(x_0 + r\cos\theta,\ x_0 + r\sin\theta)\, d\theta \qquad (\text{任意 } r \text{ 使球盤在 } \Omega \text{ 內}). \]

證明思想（Gauss 的均值定理）：考察

\[ M(r) = \frac{1}{2\pi} \int_0^{2\pi} u(x_0 + r\cos\theta, \ldots)\, d\theta. \]

對 \(r\) 求導並用散度定理（Green 恆等式）：

\[ M'(r) = \frac{1}{2\pi r} \oint_{\partial B_r} \frac{\partial u}{\partial n}\, ds = \frac{1}{2\pi r} \iint_{B_r} \Delta u\, dA = 0. \]

** \(M\) 恆定** ⇒ \(M(r) = M(0) = u(x_0)\) 。∎

散度定理把 PDE 化為邊界積分——這一步是本節所有證明的引擎（呼應微積分系列的 8.2–9.3 節向量分析）。

極值原理：內部沒有峰谷#

定理（最大值原理）：調和函數在內部不能取得最大值（除非恆常數）；最大值必在邊界上取得。

證明（反證 + 平均值性質）：若 \(u\) 在內點 \(x_0\) 取得最大值 \(M\) ，則 \(u(x_0) \ge u(x)\) 處處；但平均值性質給 \(u(x_0)\) = 圓周平均 ≤ \(M\) ，等號成立要求圓周上處處 \(= M\) ——對每個 \(r\) 成立， \(u\) 在整個盤上恆為 \(M\) ，延拓得 \(u \equiv M\) 。∎

立即推論：

  1. Dirichlet 問題解的唯一性：兩解之差邊界為零，內部最大值 ≤ 0、最小值 ≥ 0（對 \(-u\) 用同一原理）⇒ 恆零。邊界溫度唯一決定內部——與波方程（可以「埋伏」初值於內部）形成尖銳對照。
  2. 穩定性：解連續依賴於邊界數據（最大偏差被邊界最大偏差控制）——適定性（Hadamard 1923 的三條件：存在、唯一、穩定——Laplace 方程全部滿足；而「熱方程反演過去」不滿足第三條，9.1 節）。
  3. 沒有局部極值：調和函數的曲面無峰無谷（像繃緊的肥皂膜）—— \(\Delta u\) = 0 即「平均曲率張力平衡」的近似（極小曲面方程是其非線性精確版，Plateau 實驗 1873、Douglas–Radó 1931 菲爾茲級工作）。

解析性：無限光滑的奢侈品#

定理：調和函數解析（實解析：局部可展開為收斂冪級數）——且可無窮次微分。

證明路線：平均值性質 + Poisson 積分公式（圓盤上的顯式解）

\[ u(r, \theta) = \frac{1 - r^2}{2\pi} \int_0^{2\pi} \frac{f(\phi)}{1 - 2r\cos(\theta - \phi) + r^2}\, d\phi \qquad \text{（Poisson 1820s）} \]

——邊界值 \(f\) 經Poisson 核（光滑權重）積分出內部，分母的分母（ \(1 - 2r\cos + r^2 > 0\) ）可對 \(r, \theta\) 任意求導 ⇒ \(u\) 解析。∎

對照熱方程（ \(t > 0\) 瞬時平滑）與波方程（奇異性永久傳播——位移的二階導可以跳）：三兄弟的正則性階梯：波「保真」（初值多野、解多野）→ 熱「瞬時漂白」→ 調和「全時解析」。PDE 的類型決定正則性的命運——這是 9.1、9.2 節對照表的終極形態。

複分析的統一：調和函數 = 解析函數 \(w(z)\) 的實部（或虛部）（ \(w = u + iv\) 滿足 Cauchy–Riemann 方程 ⇒ \(u, v\) 調和；反之單連通域上調和 \(u\) 有調和共軛 \(v\) ）。二維 Laplace 理論就是複分析——保角映射（Riemann 映射定理 1851）解 Dirichlet 問題：把區域映到圓盤、用 Poisson 公式。Gauss、Riemann、Weierstrass 的複分析傳統與位勢理論是同一棵樹。

一個具體例題#

例題 9.3.1：單位圓盤上的 Dirichlet 問題：邊界溫度 \(f(\theta) = 1 + \cos 2\theta\) （上「半熱」下「半冷」的角向分布—— \(f > 0\) 全域但其角向不均）。求內部解 \(u(r, \theta)\) ，並用極值原理與平均值性質檢驗。

解：分離變數（極座標）： \(u = R(r)\Theta(\theta)\) ， \(\Theta\) 需 \(2\pi\) 週期 ⇒ \(\Theta = \cos n\theta, \sin n\theta\) ； \(R\) 滿足 Euler 型（5.1 節場景）

\[ r^2 R'' + r R' - n^2 R = 0 \quad\Longrightarrow\quad R = Ar^n + Br^{-n}. \]

內部正則（ \(r = 0\) 有限，3.3 節剔除 \(r^{-n}\) ）⇒ \(R = r^n\) 。疊加：

\[ u(r, \theta) = \sum_{n \ge 0} r^n (a_n \cos n\theta + b_n \sin n\theta). \]

匹配邊界： \(u(1, \theta) = 1 + \cos 2\theta\) 已是級數形式 ⇒ \(a_0 = 1\) 、 \(a_2 = 1\) 、其餘零：

\[ \boxed{\;u(r, \theta) = 1 + r^2 \cos 2\theta\;} \]

檢驗 1（調和）： \(\Delta(r^2\cos 2\theta) = (4r\cos 2\theta \cdot \frac{?}{})\cdots\) ——用 \(\Delta = \partial_{rr} + \frac1r \partial_r + \frac{1}{r^2}\partial_{\theta\theta}\) ：

\[ (r^2 \cos 2\theta)_{rr} = 2\cos 2\theta,\quad \frac1r (r^2\cos 2\theta)_r = 2\cos 2\theta,\quad \frac{1}{r^2}(r^2\cos 2\theta)_{\theta\theta} = \frac{-4r^2 \cos 2\theta}{r^2} = -4\cos 2\theta, \]

和 \(= 2 + 2 - 4 = 0\) ✓。

檢驗 2（邊界）： \(r = 1\) ： \(1 + \cos 2\theta\) ✓。

檢驗 3（平均值性質）： \(u(0, \cdot) = 1\) ；邊界平均 \(\frac{1}{2\pi}\int (1 + \cos 2\theta)\, d\theta = 1 + 0 = 1\) ✓（ \(\cos\) 的平均為零）。

檢驗 4（極值原理）： \(u \in [1 - r^2, 1 + r^2]\) ——最大值 \(2\) 在邊界 \(\theta = 0\) 取得、最小值 \(0\) 在邊界 \(\theta = \pi/2\) 取得；內部嚴格夾在 \((0, 2)\) 中（ \(r < 1\) 時 \(|r^2 \cos 2\theta| < 1\) ）——極值全在邊界 ✓。∎

解讀：內部溫度 \(u = 1 + r^2\cos 2\theta\) ：中心恆為平均值 1；角向不均隨深入而指數衰減（ \(r^n\) 因子）——邊界的高頻擾動（ \(\cos n\theta\) ）深入距離按 \(r^n\) 衰減，與 9.1 節高頻快死同構（熱方程 \(e^{-n^2 t}\) ⇔ Laplace 方程 \(r^n\) ：時間之箭 ⇔ 半徑之箭——調和函數是「沿半徑的熱方程」）。

歷史啟示#

  • 邊值思維的誕生：ODE 的世界由初值統治；Laplace 方程的世界由邊界統治——「環境塑造平衡」。這一思想在靜電學（導體邊界決定電場）、水壩設計（滲流）、成像（反問題）中全面開花。
  • 極值原理 = 適定性的憲法：Hadamard 1923 年定義「適定」時的範本正是 Dirichlet 問題；不適定問題（反演、去模糊）的求解理論（Tikhonov 正則化 1963）則是對「違憲」的補救。
  • 三型方程的三種命運：波（記憶）、熱（遺忘）、Laplace（平衡）——PDE 分類是物理時間哲學的數學地圖；三節連讀，是本書 PDE 部分的總譜。

與其他章節及領域的關聯#

  • 9.1、9.2 節：三型方程三分天下（本節是合流的第三極）。
  • 5.2 節：球面位勢的 Legendre 展開（本節極座標解的球面版）。
  • 8.3 節：Green 函數（ \(\Delta G = \delta\) ）——位勢理論的工具箱； \(G = -\frac{1}{2\pi}\ln|x|\) （2D）、 \(\frac{1}{4\pi|x|}\) （3D）——牛頓引力與庫倫定律是 Laplace 算子的基本解。
  • 4.1 節：Sturm–Liouville 的圓盤版（ \(n\) 離散譜由週期條件決定）。
  • 物理：靜電學、磁靜力學、引力位、理想流、彈性膜——五科共用。
  • 複分析：保角映射、Cauchy–Riemann、Riemann 映射定理——二維位勢理論的代數引擎。
  • 現代應用：圖像修復（Laplace 補洞）、電腦圖形（Poisson 編輯 2003）、機器學習的圖神經網絡（圖 Laplacian，2017–）——「邊界決定內部」的平衡觀念正進入 AI。

練習#

  1. 驗證 \(u = x^2 - y^2\) 、 \(u = xy\) 、 \(u = e^x\cos y\) 都是調和函數；找出它們的調和共軛（複勢 \(w\) ）。
  2. 極座標 Dirichlet：環 \(1 < r < 2\) ，內界 \(u = 1\) 、外界 \(u = 0\) 。求解（提示：徑向解 \(A\ln r + B\) ）。
  3. 用極值原理證明 Dirichlet 問題解的唯一性（完整寫出）。
  4. 對例題 9.3.1 的 \(u = 1 + r^2\cos 2\theta\) ：計算其在 \(r = \frac12\) 圓周上的平均值，並與 \(u(0)\) 比對（平均值性質的實操檢驗）。
  5. 思考題：最大值原理說調和函數「內部無峰」。但肥皂膜（極小曲面）可以有馬鞍形（內部不是極值也不是最大）。請解釋：馬鞍點與「無局部極值」是否矛盾？並說明 Laplace 方程的解可以有哪些「形狀」而仍無峰谷。

進階研讀：位勢理論的五個深化場景#

場景一：Green 函數與 Dirichlet 原理。 \(\Delta_x G(x, \xi) = -\delta(x - \xi)\) 、 \(G|_{\partial} = 0\) 的解使 Dirichlet 問題積分化： \(u(x) = -\oint_{\partial\Omega} \frac{\partial G}{\partial n} f\, ds\) 。Dirichlet 原理（Dirichlet/Riemann 1850s，Weierstrass 1870 指出漏洞，Hilbert 1900 修復）：解 = 「Dirichlet 積分」 \(\int |\nabla u|^2\) 的最小化——變分法與 PDE 的第一次大聯姻，直接孕育 Sobolev 空間與現代變分法（弱解，1930s）。

場景二：Perron 方法與次調和函數。Perron（1923）用次調和函數（ \(\Delta u \ge 0\) ，即「下凸的位勢」）的上包絡構造 Dirichlet 解——不依賴變分、純「上確界」式構造。最大值原理的反向使用（次調和函數內部無最大值）是這套方法的引擎——位勢論（potential theory）作為獨立學科由此成形。

場景三：正則性與 Wiener 判據。邊界點何時「正則」（邊界值可達）？Wiener 判據（1924）：用區域補集的電容量（容量的級數收斂性）判定——尖刺（刺入很細）處邊值失效、內點「看不見」尖刺尖端。「區域的毛刺會遮蔽邊界信息」——幾何的精細尺度決定 PDE 的行為，這是分形與位勢論的交匯（Makarov 1985 的普適性定理接續）。

場景四：隨機行走與離散 Laplacian。圖上的離散調和函數（值 = 鄰居平均）滿足離散平均值性質，其解 = 隨機行走的期望： \(u(x) = \mathbb{E}_x[u(\text{首次觸達邊界})]\) 。Dirichlet 問題 = 賭徒的期望贏額（邊界值 = 各出口的賠付）。連續極限 = 布朗運動（9.1 節熱核的表親）——調和函數、隨機行走、電阻網絡（Kirchhoff 定律）三位一體（Doyle–Snell 1984 教科書化）。機器學習的 label propagation（圖半監督學習）正是這套數學。

場景五：從 Laplace 到 mean-field。非線性位勢 \(-\Delta u = u^p\) （Yamabe 問題 1960、解決 1984）、 \(-\Delta u = e^u\) （Chern–Gaussian 幾何）、平均場博弈（Lasry–Lions 2007）——「平衡 = 場方程 + 邊界」的框架被反覆升級。圖神經網絡（GNN）的「消息傳遞收斂態」滿足圖 Laplacian 方程——調和函數的抽象在 AI 的圖結構上復活。

補充：Gauss 與 Laplace 的位勢之爭#

「位勢」（potential）一詞由 George Green（1828，自費出版的《數學分析嘗試》——磨坊主出身的自學者，此書沉寂至 1846 年由 Kelvin 挖掘）與 Gauss（1840 位勢理論論文）確立。Gauss 的《一般定理》以極值思想證明靜電平衡的存在唯一；Laplace 早三十年（1782–85）的球體引力級數展開（ \(\nabla^2 V = 0\) 的球諧解）是先驅。方法之別：Laplace 級數展開、Gauss 積分不等式、Green 積分公式——三條路至今仍是位勢理論的三大工具箱。「同一方程、三種數學人格」。

補充：為什麼「內部無峰」物理上顯然#

穩態溫度場中若內部有「熱點」（局部最大），該點向四周散熱而無源補充 ⇒ 溫度下降——與穩態矛盾。「平衡場無內源則無內峰」——極值原理是「守恆 + 平衡」的直接推論。這一物理直覺經數學化後成為強大工具：它甚至適用於非線性（ \(p\) -Laplacian）、隨機（超過程）、與圖結構——直覺一旦被定理化，就獲得跨越學科的通行證。

名詞速查表（本節出現的核心術語）#

  • Laplace 方程 / 調和函數： \(\Delta u = 0\) 及其解。
  • Dirichlet 問題：邊界值給定的內部求解。
  • 平均值性質： \(u(x_0)\) = 圓周平均；調和性的等價特徵。
  • 最大值原理：極值在邊界；唯一性與穩定性的引擎。
  • Poisson 積分公式：圓盤 Dirichlet 問題的顯式解。
  • 調和共軛： \(u \pm iv\) 構成解析函數；複分析接口。
  • Green 函數： \(\Delta G = -\delta\) ；點源響應（8.3 節的位勢版）。
  • Poisson 方程： \(\Delta u = -\rho\) ；有源位勢。
  • Wiener 判據：邊界正則性的容量準則。
  • 離散調和函數：圖上「值 = 鄰居平均」；隨機行走的期望。

年表速查（本節涉及的關鍵年份）#

  • 1782–1785：Laplace 球體引力與球諧級數。
  • 1813：Poisson 方程。
  • 1828：Green《數學分析嘗試》（位勢與 Green 函數）。
  • 1840：Gauss 位勢理論。
  • 1846：Kelvin 重提 Green（Green 函數命名）。
  • 1850s：Riemann Dirichlet 原理。
  • 1870：Weierstrass 批評 Dirichlet 原理。
  • 1900：Hilbert 修復 Dirichlet 原理。
  • 1923：Perron 方法；Hadamard 適定性；Wiener 判據（1924）。
  • 1963：Tikhonov 正則化。
  • 1984：Yamabe 猜想解決。

深化問答三則#

問：Laplace 方程與熱方程是什麼精確關係？ 答：時間方向的極限與逐點的等價。穩態： \(u_t = 0\) 把熱方程變成 Laplace（時間箭頭的終點）；特徵視角：Laplace 方程的本徵函數 \(\sin n\theta\) 深入衰減 \(r^n\) ，恰似熱方程 \(\sin nx\) 時間衰減 \(e^{-n^2 t}\) ——把 \(r\) （半徑）當作「向內流逝的時間」（記 \(r = e^{-s}\) ，則 \(\Delta = e^{-2s}(\partial_{ss} + \partial_{\theta\theta})\) 變回熱方程型）。平衡態是「已燒盡的時間」——三型方程經這種變換互通（抛物↔橢圓的字典）。

問：極值原理為什麼在波方程上失效？ 答：波方程二階時間導數使其「能量守恆」而非「值平均」——初值可以在內部「埋」一個高峰（如尖峰初位移），它將分裂為兩個行波傳播：內部極值合法存在。數學上：波的解不滿足平均值性質（它滿足的是時空光錐上的 d'Alembert 公式）。平均值的幾何決定極值原理的成立範圍——這就是為什麼「橢圓型方程有極值原理、雙曲型沒有」是定理而非經驗。

問：Green 函數的物理身份是什麼？ 答：單位點荷（點源）產生的位勢，扣除邊界感應的修正。3D 自由空間 \(G = \frac{1}{4\pi r}\) ——就是庫倫/牛頓位勢；放在接地導體盒中， \(G\) 還要加上「鏡像電荷」修正（Kelvin 鏡像法，1845–48）。「基本解 + 邊界修正 = Green 函數」——8.3 節的 \(\delta\) 響應思想在空間中的完整形態；所有靜電、引力、擴散問題的求解都以此為積分核。

（本節完）