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

平滑化 vs 奇異性傳播

餵給一條偏微分方程粗糙的資料,可能發生兩件相反的事。橢圓型與拋物型方程會瞬間把角點磨平;雙曲型方程則把每一個尖角沿著一條特徵線原封不動地永遠帶下去。這份導讀會讓你看清為什麼——以及如何讀出一條方程「必定」屬於哪一種行為。

一個角點的兩種命運

本學習階梯的前三篇導讀,談的全是大小最大值原理限定了解能在哪裡攀上頂峰、比較原理替兩個解排出高下、而哈納克不等式則禁止一個正解在某一點巨大、緊鄰之處卻又微小。這篇導讀問的是另一個問題,關乎質地。取一個一開始就粗糙的解——一個帶尖角的溫度、一條帶折點的撥動弦、一個坐在單一點上的電荷——然後看著解演化時,方程對那份粗糙做了什麼。令人驚訝的事實是:答案由方程的所決定,而且恰好只有兩種風味,且這兩種再對立不過了。

一邊坐著橢圓型拋物型方程——拉普拉斯方程與熱方程是當家代表。它們是撫平者。只要你一進到區域內部(對橢圓型而言)、或處於任何正的時刻(對拋物型而言),無論資料當初多麼鋸齒,解都是無限可微的。角點消失了。另一邊坐著雙曲型方程——波動方程是它的範本。它們是傳播者。資料裡的一個折點絲毫不會變柔;它原封不動地存活下來,以一個有限的速度沿著一條特徵線行進,而你能準確指出它已走到了哪裡。

為什麼擴散會磨掉每一個角點

先從撫平的那一側講起,因為它有一套令人滿意的機制。回想熱方程那一階梯的熱核圖像:解 u_t = k u_xx,就是把你的資料與一個高斯卷積,而高斯是無限平滑的。卷積永遠會繼承較平滑那一因子的平滑度,所以即使是帶跳躍的資料,在 t > 0 的那一瞬間也變成 C-無窮。傅立葉的觀點把同一件事說得更鋒利:粗糙住在高頻模態裡,而因子 e^(-k n^2 t) 壓垮波數為 n 的模態。對任何 t > 0,那個指數裡的 n^2 幾乎完全殲滅高諧波,剩下的東西已沒有快速的波折可供它粗糙了。

做成定理,這就是拋物型方程的瞬間平滑性質:對任何 t > 0,無論初始資料有多糟,解在空間上都是 C-無窮,甚至是實解析的。橢圓型方程在一個凍結的單一時刻帶著相對應的陳述——根本沒有時間,只有一個區域的內部——名字叫內部正則性。它的重點台詞,就是一口氣說完的橢圓正則性定理:若 Laplacian u = f 且 f 在某個開集上光滑,那麼 u 在那裡也光滑,無論 u 或它的邊界資料在別處可能多麼粗糙。一個調和函數(f = 0)在內部自動是實解析的。方程由內而外逼出了平滑。

有一個可愛的方式可以體會為什麼橢圓型方程必須這麼做。一個調和函數服從均值性質:它在某一點的值,等於它在繞著該點的任一小球面上的平均——正是替導讀一的最大值原理供能的那種平均。平均是一種平滑運算。一個被強迫在各處都等於自身局部平均的函數,不可能有孤立的尖峰、沒有摺痕、沒有角點;平均會立刻把它們供出來。平滑在這裡不是一份額外的禮物——它早已烘焙進「等於你自己的平均」這句話本身的意思裡。

為什麼波會把折點永遠帶著走

現在來看相反的命運。最乾淨的一扇窗,是直線上波動方程 u_tt = c^2 u_xx 的達朗貝爾公式:通解是 u(x,t) = F(x - c t) + G(x + c t),一個右行波與一個左行波以速度 c 滑行、形狀不變。那條公式裡沒有任何東西會平滑任何東西。若你的初始分布 F 在 x = 0 處有一個角點,那麼在時刻 t,那個角點就坐在 x = c t 處,跟它一開始時一樣銳利。波動方程把資料運送過去;它從不平均資料。一個折點是一名乘客,而不是一個有待溶解的問題。

這就是奇異性傳播的核心:在一條雙曲型方程裡,資料中的一個奇異點既不會消失、也不會憑空出現——它會行進,而且它沿著特徵線行進,正是你在一階與分類那些階梯上學會去尋找的那些曲線。資料光滑之處,解就光滑;資料有折點之處,那個折點就被守恆並以有限速度帶著走。一個解「不光滑」的那個集合,技術上的名字叫它的奇異支集,而定理說的是:對一條波動方程而言,奇異支集就只是沿著特徵線移動。你能準確預測在任何之後的時刻,粗糙會落在哪裡。

有限速度 vs 無限速度:那條線劃在哪裡

平滑與速度是同一枚硬幣的兩面。波動方程具有有限傳播速度:一個在 t = 0 時被侷限在一處的擾動,無法在錐 |x| <= c t 之外被感受到。有一道銳利的波前、一個訊號尚未到達的精確區域,由特徵線圍起——依賴域是一個貨真價實、大小有限的三角形,而不是整條直線。相對地,熱方程具有無限傳播速度:高斯熱核對每一個 x、在每個 t > 0 都嚴格為正,所以原點處的一個熱脈衝會瞬間在各處被感受到——遠處微小到無法察覺,但永遠不會恰好為零。

留意這深刻的耦合。無限速度正是讓熱方程得以平滑的原因:瞬間感受到一個遙遠的源,就是去當資料的一個全域平均,而全域平均是平滑的。有限速度正是迫使波動方程不平滑的原因:若一個訊號只能藉由沿著一條特徵線行進才能抵達你,那麼任何搭著那條特徵線進來的東西——連同它的角點——都原封不動地抵達。所以撫平者是那個一次聽見一切、把它攪成一則謠言的長舌婦;傳播者則是那個以固定速度帶著一封密封信件、原封不動送達的跑者。你沒法兩者兼得。

             SMOOTHERS                          PROPAGATORS
          elliptic / parabolic                    hyperbolic
   Laplacian u = 0,  u_t = k u_xx              u_tt = c^2 u_xx
  --------------------------------     --------------------------------
  texture  | corners vanish; C-inf    | corners kept; ride forever
  speed    | infinite (positive K)    | finite, equals c
  carrier  | global average          | F(x - c t) + G(x + c t)
  fate of  | erased, instantly       | moves along the characteristic
  a kink   |   smoothed away         |   x = c t   (singular support)
  reverse? | no  (irreversible)      | yes (time-reversible)
兩種命運並排對照。同一枚硬幣:無限速度換來平滑;有限速度則禁止它。

誠實的邊界:平滑承諾了什麼、又沒承諾什麼

三個告誡,免得你把這條規則讀過了頭。第一,平滑是內部的。橢圓正則性撫平的是一個區域的內部;它對邊界隻字未提,在那裡解可以保持得跟你所施加的邊界資料一樣粗糙。一個帶著鋸齒狄利克雷條件的調和函數,在內部是解析的,卻在邊緣被它的角點釘住——平滑唯有在你步離邊界時才綻放。拋物型的版本也一樣:對 t > 0 是平滑的,但 t = 0 處的資料則愛多醜就多醜。

第二,平滑的代價是在時間上倒著付的。擴散之所以平滑,是因為它摧毀精細細節,而一個摧毀資訊的過程無法倒著跑。這正是為什麼逆向熱方程是不適定的——那個讓向前擴散如此乖巧的東西,正是讓過去無法復原的東西。波動方程不付這樣的代價:因為它只是運送、什麼也不丟失,所以它在時間上倒著跑同樣乾淨俐落。平滑與不可逆性綁在一起;有限速度與時間可逆性綁在一起。你拆不開它們。

第三,留意奇異性這個詞本身。一個基本解——熱核,或波動方程對一個點源的回應——是對一個狄拉克 δ的回應,而它是一個分布、不是一個函數,因此它按其構造本身就是奇異的。平滑是一句關於正當資料如何傳播的陳述,而非關於核本身:熱核會平滑任何你拿來與它卷積的東西,但這個核自己在 t = 0、x = 0 處就帶著一個奇異點。而對非線性雙曲型方程,這條規則有一個確確實實的例外——光滑的資料可以製造出一個先前並不存在的奇異點。一個守恆律的光滑解可以一路變陡,直到在有限時間內形成一道激震,而在那之後,若沒有熵條件,弱解甚至不唯一。所以「奇異性只會傳播、絕不出現」是一句線性的陳述;非線性的世界能製造出它們。

如何讀出行為

把它拿來用。面對一條二階線性方程 A u_xx + 2B u_xy + C u_yy + ... = 0,你已經知道如何替它分類了:算出判別式 B^2 - A C。當初替方程命名那個型的同一個正負號,現在就免費告訴你它的質地行為——這正是當初替方程分類的回報。

  1. 算出主部(最高階部分)的判別式 B^2 - A C——正是型那一階梯裡的分類步驟。
  2. 若 B^2 - A C < 0,它是橢圓型(沒有實特徵線):預期內部平滑——無論邊界資料多粗糙,解在內部都是 C-無窮。
  3. 若 B^2 - A C = 0,它是拋物型(一族特徵線):預期對 t > 0 瞬間平滑、無限速度、且不可逆。
  4. 若 B^2 - A C > 0,它是雙曲型(兩條實特徵線):預期不會平滑——奇異點以有限速度沿著特徵線行進,且過去可復原。

帶走一句話。橢圓型與拋物型方程瞬間且全域地平滑它們的資料;雙曲型方程則把每一個奇異點原封不動地沿著它們的特徵線、以有限速度傳播——而你得到哪一種命運,可以直接從型讀出,正是你老早算過的那個 B^2 - A C。本學習階梯的下一篇、也是最後一篇導讀,會替這個定性故事補上最後一章:不只是解是否保持粗糙,而是它如何隨長時間衰減,用能量方法看著一個振動或擴散的系統沉澱向靜止。