把整階兌現
這一篇正是前四篇為之而建的。你現在手上有了試驗函數、有了作為線性規則的分布、有了作為最乾淨那種規則的 狄拉克 δ、有了能對任何東西求導的 分布導數,以及被推廣到 緩增分布 上的傅立葉轉換。每一件工具登場時都帶著一個承諾:它會在這裡兌現。兌現的成果,是一句五篇之前我們無法老實寫下的話:一個線性算子 L 的 基本解,是一個滿足 L E = delta 的分布 E。我們用這一篇拆解:為什麼那句話恰好是對的,以及它換來了什麼。
把那些字一個個慢慢拆開。L 是一個常係數的微分算子——熱算子 d/dt - k d^2/dx^2、波算子 d^2/dt^2 - c^2 d^2/dx^2,或拉普拉斯算子 d^2/dx^2 + d^2/dy^2。右邊的 delta 是一個點源:一個理想化的單位輸入,在某一瞬、某一處被打開。未知的 E 就是 基本解,系統對那一次單一、無限集中的戳動所給的響應。在舊語言裡,這一切都解析不通——delta 不是函數,所以 L E = delta 不是方程。這一整階存在的理由,就是要讓那一行上的每個符號都有意思。
「L E = delta」對 E 究竟要求什麼
這裡是讓方程變嚴謹的那一步,而它正是前面每一篇裡的同一個反射動作:絕不直接碰 E,永遠拿它去對一支探針檢驗。要在分布意義下說 L E = delta,意思是:對每一個試驗函數 phi,把 L E 與 phi 配對所得的數,等於把 delta 與 phi 配對所得的數。但 L 是由導數搭成的,而我們早已知道 分布導數 怎麼運作——它藉分部積分把每個導數轉交給平滑的探針。於是 < L E, phi > 被改寫成 < E, L* phi >,其中 L* 是形式伴隨算子(同一個算子,但帶上分部積分留下的符號),方程就成了:對每一個 phi,< E, L* phi > = phi(0)。
把那個條件大聲念出來,它就不再嚇人了。它說:不管 E 是什麼,當你把算子的導數推到一支平滑探針上、再拿 E 對結果求平均時,你必須永遠還原出探針在原點的值。這是個乾淨、可檢驗的要求——沒有無限、沒有除以零,兩邊各是一個數。一個 分布解 恰恰就是一個對每一支探針都通過這項檢驗的 E,而基本解就是源為 delta 的那個特例。前面幾篇並不只是替這個構想做裝飾;它們是這條方程根本有意義可供檢驗的唯一理由。
Want (as distributions): L E = delta Test against every probe: < L E, phi > = < delta, phi > Move derivatives onto phi: < E, L* phi > = phi(0) for all phi (L* = formal adjoint: integration by parts flips the signs) So E is a fundamental solution <=> < E, L* phi > = phi(0) for every test phi
一個響應,所有的源:卷積這把戲
為什麼要把這麼多心思全傾注在對單一點源的響應上?因為有一個事實,讓這份功夫換來巨額的回報:對一個常係數的 L,知道對一個 delta 的響應,就等於拿到了對 每一個 源的響應。理由是:一般的源 f 不過是點源的連續堆疊——f 是對所有位置 y、把 f(y) 乘上一個坐落在 y 的 delta 加總起來。因為 L 是線性的、又不在乎你身在何處,對那堆疊的響應,就是同一堆 E 的平移副本的疊加。那一群平移後基本解的疊加,恰恰就是 卷積 所算的東西。
於是這帖配方短得令人屏息:要解 L u = f,就取 u = E * f,也就是基本解與源的 卷積。施以 L,導數滑到 E 身上(卷積與算子可交換),得到 L u = (L E) * f = delta * f = f——因為與 delta 卷積,會原封不動地交還你一開始的東西。delta 是卷積的單位元,是這個世界裡那個乘法上的「1」,這 正是 為什麼基本解配得上它的名字:它是組裝出每一個解所用的積木。
三個模型算子,用嚴謹的眼光看
看著這套機械把你已經認識的方程一一點亮。對熱算子,基本解是 熱核——那個擴散開來、在 t = 0 之前為零、在此後每一刻積分都等於 1 的高斯。用分布的眼光讀,「熱核解 u_t = k u_xx,且在原點有一點熱量」這句話,就成了精確的方程 (d/dt - k d^2/dx^2) E = delta。它在任何 t > 0 都無限平滑的形狀,正是熱方程瞬間磨平這件事嚴謹的面孔:用一個 delta 戳它,下一刻響應就已經是解析的了——而且,老實說,這也意味著熱有 無限的傳播速度,因為高斯瞬間就在每一處都非零。
拉普拉斯算子說的是另一個、更銳利的故事。在平面上,Laplacian E = delta 的基本解是對數位勢 (1/2pi) log|x|;在三維裡,它是牛頓位勢 -1/(4pi|x|),也就是重力與靜電學裡那條 1/r 律。所以庫侖定律並不是另外栓上去的一塊物理——它根本就是一個橢圓型算子的基本解,而把它與電荷密度 f 卷積,就藉 u = E * f 重建出整個位勢。注意它與熱的對照:這個 E 處處平滑,唯獨 在源點除外,它在那裡爆掉。那一個壞點就是它的 奇異支撐,離開它之後,響應要多平滑有多平滑——這是橢圓型算子的指紋。
波算子為這三部曲收尾,並且拒絕表現得像前兩者中的任何一個。它的基本解 根本不 磨平任何東西——在三維裡,對原點一道閃光的響應,恰好就活在那擴張的光錐 |x| = c t 上,在錐內與錐外都為零,這正是有限速度與惠更斯原理的嚴謹陳述(你聽見一聲拍掌,俐落地響一次,而不是拖泥帶水地糊成一片)。這個基本解是一個支撐在一張曲面上的分布,根本不是任何函數。三個算子、三個截然不同的 E——而每一個 E 的奇異性都把物理編了進去:熱瞬間磨平並擴散開、拉普拉斯算子在一點奇異、波則永遠保持它的邊緣銳利。
怎麼找出 E,以及誠實的附帶細則
你究竟要怎麼算出這些 E 之一?乾淨的路子要用到上一篇。對 L E = delta 取傅立葉轉換;因為轉換把導數變成乘法、又把 delta 送成常數 1,那條微分方程就塌縮成一條代數方程 P(k) * E-hat(k) = 1,其中 P 是算子的符號。除過去就好:E-hat = 1/P(k)。然後反轉換把 E 換回來。解一條常係數偏微分方程的整個問題,就被化約為「除以一個多項式、再反轉一個 緩增 傅立葉轉換」——這是「為什麼第四篇值得花那番功夫」最乾淨不過的例證。
- 寫下你要 E 滿足的方程:L E = delta,其中 L 是常係數算子。
- 對兩邊取傅立葉轉換;導數變成乘上符號 P(k),delta 變成常數 1,得到 P(k) * E-hat(k) = 1。
- 解那道代數:E-hat(k) = 1/P(k),並決定如何處理 P 的零點(你正是在這裡挑出某一個特定的 E)。
- 反轉換把 E 找回來,再把 L u = f 的解構造成卷積 u = E * f。
現在來看誠實的界限,因為這台漂亮的機器有一座圍起來的院子。第一,E 一開始竟然存在,這本身是一條深刻的定理,而非理所當然:馬爾格朗日–埃倫普賴斯定理 保證 每一個 非零的常係數算子都有基本解——這是個確實艱深的結果,而且一旦係數開始變動、或算子變成非線性,它就悄悄失效。第二,整個卷積奇蹟靠的是 常係數 與 線性;係數一變動,你就得退回到一個綁定具體幾何的 格林函數,與單一個 E 卷積便不再夠用。
第三,也是整階最該記住的一點:這套框架骨子裡是 線性 的。E * f 之所以管用,靠的是疊加原理,而疊加原理恰恰是非線性偏微分方程所摧毀的東西。對那些「平滑的起步能在有限時間內爆發成震波」的方程,你造不出一個基本解;在那裡一個 弱解 未必唯一,你必須伸手到分布理論之外去找一個附加條件——熵條件——來挑出物理的答案。所以,珍視你在這裡所建的,但要把它擺在它真正的尺度上:分布讓線性、常係數的世界變得徹底嚴謹、乾淨得驚人,那是一場真實而有界的勝利——而不是一劑萬用溶劑。