變係數汪洋中一座可解的島
本階段至今的一切都倚賴一項奢侈:常係數。整部特徵多項式機器——猜 y = e^(rx)、讀出一個 n 次多項式、求它的根——只有在 y、y'、y''、y^(n) 前面的數字不動時才管用。一旦某個係數依賴於 x,那個幸運的猜測 e^(rx) 就不再乾淨地重現自己,整套方法隨之瓦解。這是壞消息,而且是常態:一個一般的變係數線性方程根本沒有初等的封閉形式解,這正是後面階段要動用冪級數與數值方法的原因。
但有一個華麗的例外——一個*同樣可解*、堪比常係數的變係數方程。它就是[[ode-cauchy-euler-equation|柯西-尤拉方程]](又稱等維方程或尤拉方程)。二階時它寫作 a x^2 y'' + b x y' + c y = 0,其中 a、b、c 為常數。仔細看這個模式:x^2 配 y''、x^1 配 y'、x^0 配 y。每個係數中 x 的次方,*恰好等於它所乘的導數的階數*。這唯一的巧合就是全部的祕密。
為什麼這種精確的對應如此重要?因為它使方程等維——當你把 x 乘以一個常數重新縮放時保持不變。這種縮放對稱性並非數學上的奇趣:凡是物理問題沒有內建長度尺度之處它就出現,如點電荷周圍的電場、轉角附近的溫度、楔形中的應力。在極座標或球座標下解拉普拉斯方程,會直接吐出柯西-尤拉方程。所以這並非牽強的特例——它是一整個真實問題家族的封閉形式核心。
正確的幸運猜測:y = x^r
對常係數方程,幸運的猜測是 y = e^(rx),因為指數在微分下會重現自己。對柯西-尤拉方程,正確的猜測不同——一個樸素的冪,也就是[[trial-solution-x-to-the-r|試解 y = x^r]]。原因是一場完美的配對:對 x^r 微分每次把 x 的次方降一,而方程內建的 x 因子又恰好把它推回去。看這場抵消。若 y = x^r 則 y' = r x^(r-1)、y'' = r(r-1) x^(r-2)。相乘:x^2 y'' 成為 r(r-1) x^r,x y' 成為 r x^r,而 y 就是 x^r。現在每一項都是一個數乘以*同一個* x^r。
把共同的 x^r 提出來,它就被約掉,留下一個純粹關於 r 的代數方程,處處無 x:a r(r-1) + b r + c = 0。這就是[[ode-indicial-equation|指標方程]]——柯西-尤拉方程版本的特徵多項式。它與前幾篇是同一個想法,只是把指數換成了冪。解出這個(一般為 n 次的)方程求 r,每個根都遞給你一塊基底解 x^r。最常見的失誤就住在這裡:y'' 項貢獻的是 r(r-1),而非 r^2,因為把 x^r 微分兩次先帶下 r 再帶下 r-1。漏掉這個因子,正是柯西-尤拉方程的經典錯誤。
a x^2 y'' + b x y' + c y = 0 (Cauchy-Euler form) y = x^r -> y' = r x^(r-1) , y'' = r(r-1) x^(r-2) x^2 y'' -> r(r-1) x^r x y' -> r x^r y -> x^r ---------------------------------------------------- indicial equation: a r(r-1) + b r + c = 0
根的三種情形——帶一個對數的轉折
指標方程有三種情形,精確地映照你已熟悉的常係數故事——只是帶一個意外。相異實根 r1、r2 是最容易的:兩個解 x^(r1) 與 x^(r2),通解為 y = C1 x^(r1) + C2 x^(r2)。重根 r 才是意外所在。你也許會猜第二個解是另一個冪,但它不是——它是 x^r ln x。所以通解為 y = (C1 + C2 ln x) x^r。這裡的對數扮演的角色,恰如重指數根的 t e^(rt) 中那個多出來的 t;別寫成 x^r 與 x^(r+1),那是經典的重複計數錯誤。
複共軛根 alpha ± i beta 是第三種情形,它們同樣會振盪——但是在 ln x 裡,而非在 x 裡。兩個實解是 x^alpha cos(beta ln x) 與 x^alpha sin(beta ln x),給出 y = x^alpha (C1 cos(beta ln x) + C2 sin(beta ln x))。實部 alpha 控制振幅隨 x 的某個冪增長或縮小;虛部 beta 設定解相對於那個*對數*時鐘 ln x 擺動得多快。所以當 x 從 1 走到 e 再到 e^2 時,這個波完成相同次數的擺動——它以幾何而非算術的步伐振盪。這就是二階階段那同一個歐拉公式的訣竅,只是改在對數軸上閱讀。
更深的緣由:用 x = e^t 把軸拉直
對數為何會出現,x^r 又為何根本管用?有一個單一的變數代換能解釋一切:[[substitution-x-equals-e-to-the-t|代換 x = e^t]],等價於 t = ln x。它把 x 軸作對數式拉伸,並藉此把礙事的 x 次方完美燙平。依連鎖律,x 乘 d/dx 變成 d/dt,而 x^2 乘 d^2/dx^2 變成 d^2/dt^2 - d/dt(留意那個交叉項——丟掉 - d/dt 是常見的失誤)。把這些代入 a x^2 y'' + b x y' + c y = 0,*每個 x 都抵消*,留下一個關於 t 的樸素常係數方程:a y_tt + (b - a) y_t + c y = 0。
現在你回到了熟悉的地盤。用普通的特徵方程解這個 t 方程,求出 y 作為 t 的函數,再代回 t = ln x。看那魔法落地:解 e^(rt) 變成 e^(r ln x) = x^r——這*正是* x^r 試解,得到了解釋。重根的 t e^(rt) 變成 (ln x) x^r——你的對數在此。而 e^(alpha t) cos(beta t) 變成 x^alpha cos(beta ln x)——你在 ln x 中的振盪在此。所以這三種情形並不是三條要背的規則;它們就是*同樣*的三種常係數情形,只是透過一扇對數的窗回看。
兩個並排的完整範例
讓我們把這套機器跑兩次——一次用直接的 x^r 試解,一次經由代換——好讓你相信兩者一致。取 x^2 y'' - 2x y' + 2y = 0,在 x > 0 上。這裡 a = 1、b = -2、c = 2。
- 組出指標方程:r(r-1) - 2r + 2 = r^2 - 3r + 2 = 0。注意 y'' 項給的是 r(r-1),而非 r^2。
- 分解:(r - 1)(r - 2) = 0,所以 r = 1 與 r = 2——兩個相異實根。
- 寫出通解:y = C1 x + C2 x^2。完成——沒有對數,因為兩根相異。
- 用 x = e^t 交叉核對:方程變成 y_tt - 3 y_t + 2 y = 0,其特徵方程 r^2 - 3r + 2 = 0 有相同的根。它的 y = C1 e^t + C2 e^(2t),而 e^t = x、e^(2t) = x^2——一模一樣的答案。