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

為何「型」決定了一切

判別式給了你一個標籤;現在來看看為何這個標籤統治了整條方程。同一個 B^2 - AC 的正負號,既命名了橢圓型、拋物型或雙曲型,也決定了哪種資料才是適定的、訊號傳得多快,以及解是會變得更平滑還是更尖銳。

一個正負號,三個迥異的世界

到現在,你已經做完了這一級最硬的部分。你寫下了二階線性偏微分方程的一般式 A u_xx + 2B u_xy + C u_yy +(低階項)= 0,算出了判別式 B^2 - AC,也看著標準形在你沿著特徵線前進時浮現出來。到目前為止,「型」感覺像是一個歸檔的標籤——把方程整齊地歸到橢圓型拋物型雙曲型底下的一種方式。這篇導讀講的是回報:這個標籤不是方便歸檔用的,它是統治整條方程的法律。單一一個數字 B^2 - AC 的正負號,暗中決定了你被允許給什麼資料、效應傳播得多快,以及解究竟會平滑化還是維持鋸齒。

把三個模型方程擺在眼前當作錨點,因為底下的一切其實都是關於它們的故事。橢圓型的模型是拉普拉斯方程 Laplacian u = 0(這裡是 u_xx + u_yy = 0):完全沒有時間,是純粹的平衡。拋物型的模型是熱方程 u_t = k u_xx:時間只往前走,東西會擴散。雙曲型的模型是波動方程 u_tt = c^2 u_xx:時間可逆,訊號會行進。這三者不是隨意挑的例子——根據標準形的結果,每一條二階線性偏微分方程在局部都長得像它們其中之一。弄懂這三者,你就弄懂了整片森林。

type         sign of B^2 - AC   model equation        characteristics   time?
----------   ----------------   -------------------   ---------------   ------------
elliptic     B^2 - AC < 0       u_xx + u_yy = 0       none (complex)    no time
parabolic    B^2 - AC = 0       u_t = k u_xx         one family        forward only
hyperbolic   B^2 - AC > 0       u_tt = c^2 u_xx      two families      reversible
整級內容濃縮成一張卡片:B^2 - AC 的正負號固定了型、模型、有幾族實特徵線,以及時間扮演的角色。

「型」告訴你哪種資料是被允許的

這是第一個、也是最深的結果。偏微分方程本身有太多太多的解了——回想一下,通解帶著的是任意「函數」,而不只是任意常數,所以你需要資料來挑出唯一的那一個。但你不能隨便給資料;每一種型都想要某一種特定的資料,給錯了種類就會讓問題不適定。一個問題是適定的(阿達瑪的三項要求)需滿足:解存在、唯一,且連續地依賴於資料——資料動一點點,解只動一點點。把資料配對到正確的型,你三項全拿;配錯了,你至少丟掉一項,而且往往是災難性地丟掉。

這個配對直接從型衍生出來。橢圓型方程描述沒有時間方向的平衡,所以它們想要在一個封閉區域的整個邊界上給邊界資料——狄利克雷問題(在邊緣固定 u)或諾伊曼問題(固定它的法向斜率)。拋物型方程有單向的時間,所以它們想要一個初邊值問題:時間上的起始輪廓,加上空間上的邊界條件,並且只往前推進。雙曲型方程有可逆的時間與有限速度的訊號,所以它們想要嚴格意義下的柯西問題:在初始時刻同時給定 u 和它的時間導數 u_t,然後往前推進——或往後,它不在乎。

兩則警世故事:配錯了要付什麼代價

說「配錯的資料不適定」是一回事;親眼看著它失敗則遠遠更令人難忘。看看最著名的例子,源自阿達瑪。拉普拉斯方程 u_xx + u_yy = 0 是橢圓型的,所以它想要邊界資料——但假設我們固執地對它施加柯西資料,沿著直線 y = 0 給定 u 和 u_y,彷彿它是一個波。取資料 u(x,0) = 0、u_y(x,0) = sin(n x) / n,這對大的 n 而言極其微小——隨著 n 增大,它處處趨近於零。符合的解是 u(x,y) = sin(n x) sinh(n y) / n^2。當 n 增大時,sinh(n y) 像 e^(n y) 一樣爆炸,於是這個解從幾乎等於零的資料中變得無比巨大。連續依賴性被摧毀了:資料中一個察覺不到的改變,造成了答案中無界的改變。

第二則故事是逆向熱方程——把擴散倒著跑。熱方程往時間前方走是漂亮地適定的,但若你試著問「一秒鐘前溫度長什麼樣?」,你就是在把一個平滑化的過程倒著跑,也就是「去平滑化」。一杯已經冷卻到室溫的咖啡,可能源自無數種不同的更早狀態;從模糊的現在重建尖銳的過去是不穩定的,因為現在的測量中最微小的雜訊,往回放大後就會爆掉。這就是為什麼你能把奶油倒進咖啡卻無法把它倒出來:擴散是不可逆的,而數學在任何物理之前就這麼說了。

有限速度、無限速度,以及平滑化

「型」也決定了效應傳得多快、以及尖銳特徵能否存活——而波動方程與熱方程正處於兩個相反的極端。雙曲型的波動方程具有有限的傳播速度:一個點上的擾動,要等到一個以速度 c 行進的訊號有時間抵達之前,任何地方都感受不到。這給出了一個貨真價實的依賴域——在 (x, t) 處的值只依賴於距離 c*t 以內的初始資料——而那些承載曲線,正是你從 B^2 - AC > 0 找到的實特徵族。擾動沿著特徵線以一個嚴格的速度上限前進。

拋物型的熱方程恰恰相反。它從一個點狀尖峰出發的解是熱核,一個高斯函數,在 t > 0 的那一瞬間對每一個 x 都嚴格為正——所以加熱一個點會立刻把所有地方的溫度都抬高。這就是無限的傳播速度,而它伴隨著一種瞬間平滑化效應:即使是鋸齒狀、不連續的初始資料,在時間一開始的那一刻就變得無限平滑。波動方程則反其道而行:初始形狀中一個尖銳的轉角會永遠是尖銳的轉角,沿著特徵線被帶走卻從不被抹糊。雙曲型保留特徵;拋物型抹除它們——而這個差異,再一次,是 B^2 - AC 的正負號在說話。

當「型」自己改變——以及當它崩潰時

兩個誠實的複雜情況讓這幅圖景不至於太過整齊。第一,當係數 A、B、C 依賴於 x 和 y 時,B^2 - AC 的正負號可以從一個區域翻轉到另一個區域——所以同一條方程可以在某處是橢圓型、在另一處是雙曲型,這就是混合型方程。經典範例是特里科米方程 y u_xx + u_yy = 0,這裡 B = 0,判別式是 -y:在 y > 0 處橢圓型、在 y < 0 處雙曲型、恰在 y = 0 上拋物型。它確實模擬跨音速流動,也就是空氣從次音速切換到超音速之處。「型」是局部的,而在分界線上,適定性的規則也隨之改變。

第二,整篇導讀都活在線性方程的世界裡,那裡的型由你可以事先讀出的係數固定。一旦踏入非線性方程,方程的性格就可能依賴於解本身,而美好的行為可能在有限時間內崩潰。一個純量守恆律(如柏格斯方程 u_t + u u_x = 0)的完美平滑解,可以一直變陡,直到斜率變成垂直,形成一個震波——在那之後,根本不存在古典解。越過震波,你必須允許弱解,而若沒有一個額外的熵條件來挑出物理的那一個,這些弱解並不唯一。乾淨的「型決定一切」故事,只是一棟高得多的大樓的第一層。

把開放的地平線也點名出來是值得的,好讓你保持誠實的視角。上面的分類對二階線性方程已經完全理解,但物理中某些最重要的非線性偏微分方程,至今仍部分成謎。三維納維-斯托克斯方程的平滑解是否會永遠保持平滑、抑或可能爆破,是著名的未解千禧年問題之一——這提醒著我們,越過這乾淨的一級之後,前方仍橫亙著廣袤而真正未知的國度。分類,就是你帶進那片國度的地圖。