一個能讀出整條方程的數字
在上一篇導讀裡,你認識了兩變數的 二階線性偏微分方程通式,也學到只有它最高階的導數——主部——才決定它的性格。把那個主部寫成 A u_xx + 2 B u_xy + C u_yy,其餘所有東西(一階導數、u 那一項、右端項)都歸入不影響型別的低階雜物。本學習階梯整段的問題是:這條方程住在三個截然不同世界中的哪一個?令人驚訝的是,答案僅由 A、B、C 三者所組成的單一數字決定。
那個數字就是 判別式 B^2 - A C。它的正負號——而且只有正負號——為方程命名。若 B^2 - A C 為負,方程是 橢圓型;若恰為零,方程是 拋物型;若為正,方程是 雙曲型。這個數字的大小無關緊要,重要的只是它落在零之下、零之上,還是恰在零上。這就是本學習階梯得名的「大三分」,整個學科都繞著它組織起來。
它挑出的三條模型方程
橢圓、拋物、雙曲這些名字並非隨意選取——也不是(雖然這個說法很誘人)因為解的形狀像那些圓錐曲線。它們來自代數:把 u_xx、u_xy、u_yy 換成符號 X^2、X Y、Y^2,你就得到 A X^2 + 2 B X Y + C Y^2,而正是這個式子所對應的圓錐曲線 A X^2 + 2 B X Y + C Y^2 = 常數,會恰好按照 B^2 - A C 的正負號而成為橢圓、拋物線或雙曲線。這套分類,就是穿著偏微分方程戲服的圓錐曲線代數。每一型都有一個最乾淨的代表,就是那條 模型方程,每當你聽到型別的名字時都該把它記在心裡。
A u_xx + 2 B u_xy + C u_yy + (lower order) = 0
B^2 - A C < 0 ELLIPTIC model: u_xx + u_yy = 0 (Laplace)
B^2 - A C = 0 PARABOLIC model: u_t = k u_xx (heat)
B^2 - A C > 0 HYPERBOLIC model: u_tt = c^2 u_xx (wave)把這張表套在每個模型上。拉普拉斯方程 u_xx + u_yy = 0 有 A = 1、B = 0、C = 1,故 B^2 - A C = -1 < 0:橢圓型,穩態平衡的方程。熱方程 u_t = k u_xx,以 (x, t) 為兩個變數來讀,有一個 u_xx 項卻完全沒有 u_tt 項,故 A = k、B = 0、C = 0,B^2 - A C = 0:拋物型,不可逆擴散的方程。波動方程 u_tt = c^2 u_xx 重排成 c^2 u_xx - u_tt = 0,得 A = c^2、C = -1(符號相反),故 B^2 - A C = c^2 > 0:雙曲型,有限速度訊號的方程。三個正負號、三種行為、三條有史以來最重要的方程。
名字的由來:數一數特徵方向
還有第二種、更動態的方式來讀判別式,它為下一篇導讀作準備。去尋找那些讓主部退化的特殊方向——也就是 特徵方向。試一個只沿方向 dy/dx = m 變化的解,會導出二次式 A m^2 - 2 B m + C = 0(正負號隨慣例而異,但重點在判別式)。解出斜率 m,得 m = (B ± sqrt(B^2 - A C)) / A。這條小二次式有幾個實數解,由——你猜對了——B^2 - A C 的正負號掌管。
- 若 B^2 - A C > 0(雙曲型),平方根為實數且非零,於是有「兩條」相異的實特徵斜率——兩族讓訊號沿之傳播的實曲線。這正是雙曲型方程會傳播的原因。
- 若 B^2 - A C = 0(拋物型),平方根消失,兩個斜率併成「一條」重根的實特徵——單一的承載方向,正是夾在另外兩型之間的臨界情形。
- 若 B^2 - A C < 0(橢圓型),平方根為虛數,於是完全「沒有」實特徵——沒有任何實曲線承載資訊,這正是橢圓型方程描述「瞬時、處處耦合、毫無傳播」之平衡狀態的原因。
所以判別式實際上就是在數實特徵族的數目:兩個、一個,或零個。雙曲型方程擁有讓訊號沿之奔跑的實軌道——波動方程可因式分解為一個左行波與一個右行波,各自沿著自己的特徵前進,正如你在下一級的傳輸方程裡所見。橢圓型方程一個也沒有,所以任一處的擾動會被處處同時感受到。拋物型恰好坐在刀鋒上,只有單一條重根軌道。下一篇導讀會把這些特徵曲線化為新座標,並用它們把每一型摺疊成它的 標準型。
如何分類任一條方程——包括會變型的方程
這就是整套程序的縮影,是你將執行數百次的動作。比方說取 u_xx + 4 u_xy + u_yy = 0。讀出 A = 1,由交叉項 4 u_xy = 2 B u_xy 得 B = 2,且 C = 1。於是 B^2 - A C = 4 - 1 = 3 > 0:雙曲型。或取 u_xx + 2 u_xy + u_yy = 0:此時 B = 1,B^2 - A C = 1 - 1 = 0:拋物型。配方永不改變——僅從二階項認出 A、B、C,算出 B^2 - A C,讀正負號。低階項只是搭順風車,永遠不會改變型別。
一個關鍵的誠實:當 A、B、C 不是常數而是 x、y 的函數時,判別式 B^2 - A C 本身就是位置的函數,因此同一條方程可以在某一區域是橢圓型、在另一區域是雙曲型。這些就是 混合型方程,其中最著名的是 特里科米方程 y u_xx + u_yy = 0:此處 A = y、B = 0、C = 1,故 B^2 - A C = -y。它在 y > 0 處是橢圓型、在 y < 0 處是雙曲型,而恰好沿著直線 y = 0 是拋物型——一條同時身兼三型的方程,這正是穿音速氣流(次音速與超音速區域並肩相鄰)如此微妙的原因。
為何這單一標籤決定一切
如果一套分類只是按代數把方程歸類,那它就太貧弱了。三分之所以贏得核心地位,是因為型別掌控行為——解長什麼樣子,而最重要的是,什麼樣的資料能讓問題可解。這就是 依型別而定的適定性,是接下來兩篇導讀的核心,現在先預覽一下很值得,好讓判別式感覺起來不只是記帳。每一型想要不同種類的資料,而給錯種類不只是更難——它在數學上是壞掉的。
幾個關於這些世界有多不同的誠實預覽。橢圓型問題(拉普拉斯)想要繞著一個區域一整圈的邊界資料,並給出單一光滑的平衡;但為拉普拉斯方程指定柯西資料——在一段開弧上的值與法向導數——卻是出了名的不適定,其解會在資料最微小的抖動下爆炸。拋物型問題(熱)沿時間向前奔跑,並把任何粗糙的起始瞬間抹平,但把熱方程沿時間倒著跑卻是不適定的,實質上無法穩定地進行。雙曲型問題(波)以有限速度沿特徵傳播,保留尖銳特徵而非抹平它們,並且樂於接受柯西資料——初始位置與速度。同一套三分,三種截然不同的命運。
所以判別式不是為了整潔而貼上的標籤;它是一份預報。算出 B^2 - A C,你在解出任何東西之前就已經知道:該期待抹平還是尖銳波前、無限還是有限速度、邊界資料還是初始資料、平衡還是演化。這份預測力,正是這單一數字坐鎮整個學科門口的原因——也是本學習階梯其餘導讀要逐一拆解「每一型如何活出其判別式所指派之命運」的原因。