JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

表示公式與杜阿梅原理

本梯級的壓軸:基本解與格林函數如何組裝成單一的表示公式,把任一點的解寫成對資料的積分;而杜阿梅原理又如何在時間方向做同一件事——把驅動項重塑為一串脈衝之流,每個脈衝都在你已信賴的求解器下演化。

本梯級一路在搭建什麼

你如今已經把本梯級的整套工具備齊。你認識了「對點源的回應」這個想法;你寫下了拉普拉斯、熱與波三種算子的基本解——即對狄拉克 δ 的回應;你又把它升級成邊界值問題的格林函數,方法是加上一個能消去邊界資料的校正項;你還學了鏡像法,能在半空間與球上把那些校正項明確寫出。每一塊都在為單一的源回答同一個問題。這篇收尾的指南把這一切兌現:只要你知道對一個點源的回應,你就一次知道了對每一種源分布的回應。

完成兌現的機制有兩張臉,一張屬於空間、一張屬於時間,而它們暗地裡是同一個想法。在空間上,格林恆等式把基本解變成一個表示公式,把任一內部點的解值印成對資料的單一積分。在時間上,杜阿梅原理把一個無源的演化求解器變成對任意驅動都管用的求解器,方法是把對一串瞬間脈衝的回應積分起來。兩者都是把疊加原理用到極致:把資料切成點,逐點回應,再加總。這就是貫穿整個梯級的單一脈絡。

格林恆等式:長大成熟的分部積分

空間表示公式背後的引擎,是一對記帳用的恆等式,你可以從散度定理一行推出。它們不過是多變數的分部積分,讓你把一個拉普拉斯算子從一個函數搬到另一個函數上,代價是一個邊界項。格林第一恆等式搬一階導數;格林第二恆等式——對稱的那一個——搬兩階,而它正是這裡真正幹活的那一個。

Green's 2nd identity (over a region D with boundary S):

  integral over D of ( u (Laplacian v) - v (Laplacian u) ) dV
      =  integral over S of ( u v_n - v u_n ) dS

  where  _n  is the outward normal derivative on the boundary S.

Notice it is SYMMETRIC in u and v: swapping them flips a sign.
That built-in symmetry is the seed of Green's-function reciprocity.
格林第二恆等式:拉普拉斯算子可以從 u 傳到 v,只留下一個邊界積分。它的對稱性正是讓格林函數對稱的根源。

接下來是把它變成 u 本身公式的把戲。把第二個函數 v 選成以某固定內部點為中心的基本解——正是你先前打造的「對 δ 的回應」這個物件。那麼 Laplacian v 就是坐在那一點的狄拉克 δ,於是 v(Laplacian u) 這一項對著方程 Laplacian u = 源 積分,而 u(Laplacian v) 這一項對著 δ 積分,便直接把 u 在中心的值挑出來。整個積分坍縮成關於一個數的陳述:u 在所選點的值。那次坍縮——一個 δ 伸進來取樣 u——就是全部的魔法。

表示公式:解=資料對核的積分

把那次坍縮算到底,你就得到表示公式:u 在某內部點的值,等於源對基本解的體積分,加上資料對同一個核及其法向導數的邊界積分。慢慢讀,因為它是一切的妙語結局:某點上未知的解,被完全用已知的東西寫出——內部的源與邊界上的資料——每一項都依「那裡的單位源對這一點影響多強」加權。核就是那份影響;公式不過是把影響加起來。

光禿禿的基本解仍會留下一個彆扭的邊界項,涉及 u 的法向導數,而那通常是你不知道的。這正是格林函數被造出來要補的缺口。把裸基本解換成完整的格林函數——基本解加上讓它在邊界上歸零的校正項——那個彆扭的項就被關掉。對一個狄利克雷問題,你剩下一個乾淨的公式,只用到指定的邊界值與格林函數的法向導數。在球上,那個法向導數有閉形式可算,結果恰好是你也許曾把它當成獨立奇蹟見過的帕松積分公式——如今揭示為一台通用機器的特例。

對稱、互易,以及核的意義

回頭看格林第二恆等式,注意它對 u 與 v 是對稱的。把以兩個不同點為中心的兩個格林函數餵進去,那對稱性便免費交給你一個漂亮的事實:格林函數對它的兩個變元是對稱的——G(x, y) = G(y, x)。在物理上這就是互易性:位於 y 的單位源在 x 處施加的影響,等於位於 x 的單位源在 y 處施加的影響。在這裡擺動一個電荷、量那裡的電位,再把角色對調,你會讀到同一個數。這不是公式的巧合;它是拉普拉斯算子的對稱性(它的自伴性)被看見了。

互易性不只是優雅——它是一個承重的合理性檢查,也是一條計算捷徑。它告訴你這個核真的是微分算子的反算子:先施加算子、再對格林函數積分,會還原成恆等運算,恰如先乘一個矩陣、再乘它的反矩陣一樣。它也把你的工作量減半,因為你能用單一個精選的源算出許多點上的回應,而不必反過來做。它還是你早已從有限維線性代數信賴的事實的離散迴響:一個對稱矩陣有一個對稱的反矩陣。

杜阿梅原理:同一個想法,但在時間裡

到目前為止解的都是穩態問題——拉普拉斯與帕松,那裡沒有時間。像熱方程波動方程這樣的演化方程加進了時間,也隨之帶來一種新的源:一個在每一瞬間都持續注入的驅動項,而不只是一個初始條件。杜阿梅原理就是表示公式在這個情境下的孿生兄弟,而它倚賴完全相同的一招。把運轉中的源切成一串短暫的瞬間脈衝。在時刻 s 點火的脈衝表現得像一個在那一刻開啟的全新初始條件,而從那以後它就只在無源方程下演化剩下的時間 t - s——用的正是你已信賴的求解器,熱核或達朗貝爾公式。

  1. 在某一瞬間 s 把驅動凍結。讀出它在那裡注入的量 f(x, s) ds,把那一片重新詮釋為一個在時刻 s 開啟的全新初始條件。
  2. 讓那個凍結的脈衝在無源方程下向前演化剩下的時間 t - s,用的正是你對這個算子已經精通的齊次求解器。
  3. 疊加:把這些各自獨立演化的脈衝對從 0 到 t 的每個點火時刻 s 積分。用 S(tau) 表示對時間 tau 的無源演化,受迫的那部分就是 u(x, t) = integral from 0 to t of S(t - s) f(., s) ds。
  4. 若初始資料也不為零,再用疊加加上它自己單獨演化的回應——驅動由杜阿梅處理,起點由普通演化處理,兩者直接相加。

用核寫出來,熱方程對源的回應是 u(x, t) = integral from 0 to t of [ integral of G(x - y, t - s) f(y, s) dy ] ds。這正是 ODE 的參數變異法在 PDE 上的對應,而內層的物件 G(x - y, t - s) 是時間方向的因果格林函數——有時就把杜阿梅公式當作它的定義。因果這個詞帶著真實的內容:核只活在 t - s > 0 之處,編碼了「現在的溫度只能依賴過去加入的熱、絕不依賴未來」。杜阿梅原理說白了就是把時間當作其中一個變數、並把因果性內建到「核被允許非零的位置」之中的格林函數思維。

誠實的限度,以及貫穿一切的那一個想法

表示公式與杜阿梅從頭到腳都是線性陳述——它們之所以奏效,全靠對個別點源或個別脈衝的回應能夠疊加。方程一旦變成非線性,乾淨的求和就崩解:對一個依賴 u 自身的源 f,杜阿梅公式不再是最終答案,而變成一個你必須朝不動點迭代的積分方程,而對完整非線性解的真正表示公式,一般並不存在。這正是那道牆,它讓非線性守恆律的光滑解能在有限時間內形成激波,而在那之後,弱解若沒有額外的熵條件甚至不唯一。疊加是線性問題的特權;用它時要心裡有數。

因果性與適定性也約束了你能寫什麼。因果格林函數在時間上只指一個方向,正是因為向前的熱與波問題適定、而它們向後的版本不適定——把擴散倒著跑是不適定的,核的 t - s > 0 支集就是這個事實化成的幾何。而時間核的形狀本身就改變了物理:對奇數空間維的波動方程,對一個脈衝的回應是一層銳利的殼,掠過即逝,這就是惠更斯原理——一個乾淨、無回聲的訊號——而在偶數維,回應會留下一條拖長的尾巴。表示公式不只記下答案;它的核的支集就是傳播律。