那個讓偏微分方程變得可解的大膽猜測
你已經學會了給一個偏微分方程分類,能看出它是拋物型、雙曲型還是橢圓型。現在登場的,是第一個真正的求解技術,也是你日後會一次次伸手去拿的那一個:**分離變量法**。這個策略大膽得近乎放肆。一個偏微分方程把兩個獨立變量纏在一起——對於桿中的熱方程,是位置 x 和時間 t,被那些偏導數打成一個結。我們不去解這個結,而是乾脆假設它從來就不存在:我們去找形如「一個只含 x 的函數」乘以「一個只含 t 的函數」的解,u(x, t) = X(x) T(t)。這要求似乎太過分了。可了不起的是,對於合適的問題,它真的管用。
把畫面落到實處。取一根長度為 L 的細金屬桿,沿 x 軸從 0 擺到 L,兩端都埋在冰裡,因而被固定在溫度零。在最初那一刻,你交給它某個溫度分布 f(x)——也許中間燙、兩端冷。熱方程 du/dt = k d^2u/dx^2(其中 k > 0 是熱擴散率)支配著此後的每一刻:每一點處溫度隨時間的變化率,正比於那裡空間上的彎曲程度。分布向上鼓起(上凸,偏導數的二階為負)的地方溫度下降;向下凹陷(下凹)的地方溫度上升。熱從熱處流向冷處,凸包逐漸鬆弛。我們的任務,是從這一個初始形狀 f(x),預測出整個未來 u(x, t)。
把猜測代回去,再把變量撬開
把猜測 u = X(x) T(t) 餵進 du/dt = k d^2u/dx^2。時間導數只碰 T,於是 du/dt = X(x) T'(t)。兩次空間導數只碰 X,於是 d^2u/dx^2 = X''(x) T(t)。熱方程變成 X T' = k X'' T。現在每一項都是「一個 x 的東西」乘「一個 t 的東西」,我們可以把它們分類:兩邊同除以 k X T,把所有含 t 的項歸到一邊,所有含 x 的項歸到另一邊。結果是 T'/(k T) = X''/X。盯著這個方程,因為下一句話就是整個訣竅。
左邊只依賴 t;右邊只依賴 x。然而它們對所有 x 和所有 t 同時相等。一個純粹的 t 的函數要處處等於一個純粹的 x 的函數,唯一的可能是兩者都等於同一個常數——你撥動 t,左邊不能動(一動就破壞了等式,而凍結在 x 裡的右邊拒絕跟隨),反之亦然。所以每一邊都等於一個固定的數。我們把它記作負 lambda,寫成 T'/(k T) = X''/X = -lambda,那個負號是有遠見地選的,因為(我們將會看到)邊界條件會逼著 lambda 取正。一舉之間,一個含兩個變量的偏微分方程,裂變成了兩個分立的常微分方程:X'' + lambda X = 0 與 T' + k lambda T = 0。
邊界條件挑選出被允許的模式
現在來解空間方程 X'' + lambda X = 0——但不是在真空裡解。桿的兩端被固定在零,u(0, t) = 0 與 u(L, t) = 0 對所有 t 成立,而既然 u = X T,這就逼出 X(0) = 0 與 X(L) = 0(我們丟棄平凡的 T = 0)。這是 X 上的齊次狄利克雷邊界條件。常微分方程 X'' + lambda X = 0 有一個你在二階線性常微分方程裡熟知的特徵方程 r^2 + lambda = 0。它的解依賴於 lambda 的符號,而把這三種情形都查一遍,正是把除一支特殊的離散族之外的一切都殺掉的過程。
X'' + lambda X = 0, X(0) = 0, X(L) = 0 Case lambda < 0 (write lambda = -mu^2): X = A cosh(mu x) + B sinh(mu x) X(0)=0 => A = 0; X(L)=0 => B sinh(mu L) = 0 => B = 0. Only X == 0. Rejected. Case lambda = 0: X = A + B x X(0)=0 => A = 0; X(L)=0 => B L = 0 => B = 0. Only X == 0. Rejected. Case lambda > 0 (write lambda = mu^2): X = A cos(mu x) + B sin(mu x) X(0)=0 => A = 0; X(L)=0 => B sin(mu L) = 0. Nontrivial (B != 0) needs sin(mu L) = 0 => mu L = n pi, n = 1, 2, 3, ... ==> lambda_n = (n pi / L)^2, X_n(x) = sin(n pi x / L)
讀懂剛剛發生的事。邊界條件起了過濾器的作用:在連續多個可能的 lambda 之中,只有離散的階梯 lambda_n = (n pi / L)^2(n = 1, 2, 3, ...)通過,每一個都帶著自己的空間形狀 X_n(x) = sin(n pi x / L)。這些特殊的數就是這個問題的本徵值,那些正弦形狀就是它的本徵函數——桿的自然模式,完全由它的幾何與兩端被固定的方式所確定。它們恰恰是長度為 L 的弦能容納的駐波形狀:一個鼓包、兩個、三個。手握每一個 X_n,時間方程 T' + k lambda_n T = 0 就是一行就能解完的可分離一階常微分方程,它的解是純粹的衰減,T_n(t) = e^(-k lambda_n t)。於是一個分離解就是 u_n(x, t) = sin(n pi x / L) e^(-k (n pi / L)^2 t)。
疊加:用這些模式搭出任意形狀
單獨的每一個 u_n 都是貨真價實的解,但單獨一個正弦幾乎永遠配不上你拿到的那個雜亂初始分布 f(x)。這裡,熱方程的線性救了我們。因為方程是線性且齊次的,疊加原理成立:任意多個解之和仍是解。於是我們組成最一般的組合,u(x, t) = 對 n 求和的 b_n sin(n pi x / L) e^(-k (n pi / L)^2 t),其中常數 b_n 仍然自由。無論 b_n 怎麼選,這個和都自動滿足偏微分方程與兩個邊界條件。只剩一個條件還沒用上——初始形狀——而它恰好足以把每一個 b_n 都定下來。
令 t = 0。每個衰減因子都變成 e^0 = 1,於是要求 u(x, 0) = f(x) 讀作 f(x) = b_n sin(n pi x / L) 之和。這就是 f 在區間 [0, L] 上的一個傅立葉正弦級數——傅立葉那幾篇裡的半區間展開,如今帶著用途登場。這些係數不是猜出來的;它們靠本徵函數的正交性算出來。不同 n 的正弦 sin(n pi x / L) 兩兩正交:sin(n pi x / L) sin(m pi x / L) 在 [0, L] 上的定積分,當 n 不等於 m 時為零,相等時等於 L/2。把兩邊同乘 sin(m pi x / L),從 0 到 L 積分,整個和就坍縮了——除第 m 項外每一項都積成零——剩下 b_m = (2/L) 乘以 f(x) sin(m pi x / L) 從 0 到 L 的積分。
- 分離:代入 u = X(x) T(t),除以 k X T,令兩邊都等於常數 -lambda。偏微分方程裂成 X'' + lambda X = 0 與 T' + k lambda T = 0。
- 解本徵值問題:在空間常微分方程上施加 X(0) = X(L) = 0;只有 lambda_n = (n pi / L)^2 倖存,本徵函數為 X_n = sin(n pi x / L)。
- 對每個模式解時間常微分方程:T_n(t) = e^(-k lambda_n t),一個衰減率由 lambda_n 決定的純指數衰減。
- 疊加並擬合初始數據:求和 b_n X_n(x) T_n(t),再令 t = 0,挑選 b_n 使這個和等於 f(x)——它們就是傅立葉正弦係數 b_n = (2/L) 乘以 f(x) sin(n pi x / L) 從 0 到 L 的積分。
讀懂這個答案,也讀懂它的小字
完成的解是 u(x, t) = b_n sin(n pi x / L) e^(-k (n pi / L)^2 t) 之和,它講了一個鮮活的物理故事。每個模式都指數衰減,但衰減率是 k (n pi / L)^2——正比於 n 的平方。那些曲率尖銳、n 很高的扭動模式,兇猛地飛快消融;而平緩的第一模式 sin(pi x / L) 褪色最慢。所以無論你以多麼參差的分布起步,桿幾乎一瞬間就抹平成那個最低的單峰形狀,隨後悄悄地塌向零。這就是一個你早已刻在骨子裡的事實的數學面孔:尖銳的溫差最先被抹平,而熱永遠模糊細節。這個本徵函數展開,已把單獨一個偏微分方程,變成了一支由各自獨立衰減的音符組成的無窮樂隊。
對這個級數承諾了什麼、沒承諾什麼,要誠實。b_n 恰是傅立葉正弦係數,所以若 f 帶有跳躍,級數最初的那一片就攜帶著你先前見過的同一個吉布斯過衝——在 t = 0 處,截斷和會在拐角上過衝。但這裡物理是寬容的:對任意 t > 0,因子 e^(-k (n pi / L)^2 t) 把高 n 項壓得如此猛烈,以至於熱一開始流動,級數就在那一瞬間變得無窮光滑。吉布斯振鈴只活在初始那一刻,此後所有時間裡都已消失。這種「瞬時抹平」是拋物型熱方程獨有的;用一模一樣的分離方法求解的波動方程,並沒有這樣的阻尼因子——它的模式不衰減而是永遠振盪,所以初始數據裡的一個拐角會原封不動地一路同行。同樣的方法,截然不同的物理。
在你把這套帶走之前,還有兩條誠實的提醒。其一,這個配方依賴於齊次邊界條件——兩端都為零。如果兩端反而被固定在非零的溫度上,乘積形式直接失效,因為常數無法分解成 X(x)T(t);標準的補救是先減去穩態解(桿最終安定到的那條直線),讓剩下的問題重新擁有零邊界,再去分離它。其二,這套「先乘積、再疊加」的構造能重現每一個合理的 f、並且是唯一的解,並不是揮揮手就顯然的——它靠的是一條貨真價實的定理。正弦本徵函數的完備性保證了展開存在,極大值原理保證了熱方程是適定的,所以你搭出來的那個答案,就是那個答案。這套方法之所以是一名忠實的僕人,正因為有那套理論在它背後撐著。