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

一般的二階線性偏微分方程

一個模板——A u_xx + 2B u_xy + C u_yy + 低階項 = 0——竟同時暗藏了熱方程、波動方程與拉普拉斯方程。認識這個模板,明白為何只有它最高階導數的部分對「型」有意義,並讀出那個單一的數字 B^2 - A C,它將把每一條這樣的方程分入三個家族。

一個容納三大經典的模板

到目前為止,你已用特徵線法把輪廓沿著河流帶往下游,並一條承載曲線接一條地馴服了一階方程。現在我們提高一階。物理交給你的方程中,有極大一部分都套進同一個模子:A u_xx + 2B u_xy + C u_yy + D u_x + E u_y + F u + G = 0。這裡 u(x,y) 是未知量,而係數 A、B、C、D、E、F、G 可以依賴 x 與 y,但不能依賴 u 或它的導數——這正是它之所以是一條 線性偏微分方程的原因。你正要開始的整個學習階梯,就是這一個模板如何僅憑它頂端係數的大小,分裂成三個截然不同世界的故事。

它真的容納了這三大經典。令 A = 1、C = -1,其餘係數全為零,再把 y 換成 t:你得到 u_xx - u_tt = 0,這是喬裝過的波動方程。令 A = 1、C = 1,其餘為零:你得到 u_xx + u_yy = 0,也就是拉普拉斯方程,Laplacian u = 0。而一個略為不同的套法,A = 1 配上單一個一階時間項,給出 u_t = k u_xx,即熱方程。同一個模板,三張面孔——振動、平衡、擴散。本學習階梯的任務,就是解釋它們明明長得如此相像,行為卻為何如此不同。

A u_xx + 2B u_xy + C u_yy  +  D u_x + E u_y + F u + G = 0

  A=1, C=-1, rest 0   -->  u_xx - u_tt = 0     (wave)
  A=1, C= 1, rest 0   -->  u_xx + u_yy = 0     (Laplace)
  A=1, plus a u_t     -->  u_t = k u_xx        (heat)
三條著名方程,一個二階線性模板——只憑最高階導數的係數就被區分開來。

只有頂樓重要:主部

這是分類的第一個令人驚訝的真相:低階項 D u_x + E u_y + F u + G 完全不影響「型」。只有三個最高階導數的項,A u_xx + 2B u_xy + C u_yy,決定了一切。我們給這一塊取個名字——主部——並把它當作方程的引擎室。加上一個 u_x 項或一個 F u 項會改變解的細節,但它無法把波動方程變成熱方程,或變成拉普拉斯方程。「型」純粹是頂樓的性質。

為什麼最高階導數說了算?直覺上,二階導數描述的是 u 中最尖銳、變化最快的部分——曲面如何彎曲——而正是這份彎曲,掌控著資訊如何移動:一個折角能否傳播、一個擾動是停留在局部還是擴散到處處、甚至「未來」這個問題問起來究竟有沒有意義。低階項會推一推、抑制一下答案,但「方程如何傳播、在何處行為良好」這整套結構,完全寫在 A、B、C 之中。所以從現在起,要分類,我們就瞇起眼,直到只有主部還在焦點裡。

把一切分類的那個數字:B^2 - A C

從主部 A u_xx + 2B u_xy + C u_yy,我們構造一個單一的數字,即判別式:B^2 - A C。它的正負號——正、零或負——把方程歸入三種之一。若 B^2 - A C 為負,方程是橢圓型(拉普拉斯方程是其模型)。若它恰為零,方程是拋物型(熱方程是其模型)。若它為正,方程是雙曲型(波動方程是其模型)。三個符號,三個世界。下一篇導讀將整篇用來拆解這個三分法,所以這裡我們只是先認識它,並在我們的經典方程上驗一驗。

驗一驗這三大經典,名字便自己掙得了。對拉普拉斯方程 u_xx + u_yy = 0,我們有 A = 1、B = 0、C = 1,所以 B^2 - A C = 0 - 1 = -1 < 0:橢圓型,正如所許諾。對熱方程,把 u_t = k u_xx 放在 (x,t) 裡讀、只有 u_xx 出現,主部有 A = k 而 C = 0(沒有 u_tt 項),給出 B^2 - A C = 0 - 0 = 0:拋物型。對波動方程 u_tt - c^2 u_xx = 0,頂端係數正負相反,所以 A C < 0 且 B^2 - A C > 0:雙曲型。這些標籤不是裝飾——它們直接來自圓錐曲線,在那裡 B^2 - A C 正是憑這個符號把橢圓、拋物線與雙曲線彼此區分開來。

  1. 剝掉所有低階項;只留下 A u_xx + 2B u_xy + C u_yy——也就是主部。
  2. 讀出 A(u_xx 的係數)、B(u_xy 係數的一半)以及 C(u_yy 的係數)。
  3. 計算判別式 B^2 - A C。
  4. 為負 -> 橢圓型;為零 -> 拋物型;為正 -> 雙曲型。這一個符號就定下了「型」。

為何「型」值得這般大費周章

如果「型」只是個標籤,那回報就太單薄了。它不是——它是對行為的一項預言。判別式暗中在數方程的特徵線:雙曲型方程有兩族實的特徵線,拋物型有一族,橢圓型則一族也沒有(那些曲線是虛的)。這些曲線是資訊行走的道路,所以它們的數目,正是「沿尖銳波前傳遞清晰訊號的波」、「全無波前的熱流」與「每一點同時感受到其他每一點的穩態場」三者之間的分野。數特徵線,正是本學習階梯的第三篇導讀。

「型」也決定了哪些問題問起來才算公平——這是第四篇導讀的深層主題。像波動方程這樣的雙曲型方程,取初始資料後以有限速度向前推進;這是適定的。像熱方程這樣的拋物型方程也隨時間向前運行,邊走邊平滑,但把它倒著運行卻是出了名的不適定災難:微小的誤差會爆炸。而要求像拉普拉斯方程這樣的橢圓型方程從柯西資料向前推進,同樣是不適定的——橢圓型問題要的是繞著一個封閉區域一圈的邊界值,而不是一個初始切片。同一個模板;對於「什麼才算一個合理、適定的問題」卻有著完全不同的規則。

繼續攀登前的誠實提醒

三條要帶上階梯的誠實話。第一,當 A、B、C 依賴 x 與 y 時,B^2 - A C 的符號可能隨地點而改變,所以單一條方程可以在某一區域是橢圓型、在另一區域是雙曲型。那是真正的混合型方程,而其經典範例——特里科米方程 u_yy = y u_xx——正是流體跨越音速時的跨音速流動模型。「型」是局部的,並非總是整體的。

第二,這整套分類是為線性方程打造的,其頂端係數是 x 與 y 的固定函數。一旦係數依賴於 u 本身——一條非線性方程——「型」就可能隨著解的演化而改變,而守恆律的平滑解甚至能在有限時間內撞成震波。所以把 B^2 - A C 當作線性世界的萬能鑰匙,並預期非線性的那幾級會對它加上但書。第三,這裡的配方活在兩個變數 x 與 y 中;在三個或更多變數裡,你無法把一切化約成單一個數字,而是改從算子的符號讀出「型」——一個二次型,它特徵值的符號推廣了 B^2 - A C。那正是本學習階梯的收尾導讀。