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

特徵方程:三種情形

猜 y = e^(rt),一個常係數方程就塌縮成一個簡單的二次式。它的兩個根——相異、重根,或複數——直接交給你一組基底,並把整個問題恰好劈成三種情形。

一個幸運的猜測,把微積分變成代數

從本層前兩篇指南中,你已經知道答案的形狀:一個齊次二階方程擁有一個二維解空間,所以只要你能找到兩個線性獨立的解——由非零的朗斯基行列式確認——它們的組合就給出每一個解。你還沒有的,是真正「找出」那兩個解的辦法。本篇指南為最重要的特殊情形提供了它,而且只用一個便宜得驚人的把戲。

把注意力限制在常係數方程 a y'' + b y' + c y = 0,其中 a、b、c 是純粹的數,不是 x 的函數。這個把戲,就是「猜」解是一個指數,y = e^(rt),其中 r 是某個尚待求出的常數。為什麼這樣猜?因為指數是唯一一種導數只是自身複本乘上常數的函數:對 e^(rt) 微分永遠不會產生新的形狀,只會多出一個因子 r。於是方程的三項都會共有同一個 e^(rt),它就能被提出來。

看它發生。把 y = e^(rt) 代入,得 y' = r e^(rt)、y'' = r^2 e^(rt)。代入:a r^2 e^(rt) + b r e^(rt) + c e^(rt) = 0。每一項都帶有因子 e^(rt),而它永遠不為零,所以我們可以把它除掉,微積分就這樣蒸發了。剩下的是 a r^2 + b r + c = 0——一個關於 r 的普通二次式。微分方程被換成了一個你從學生時代就會解的代數問題。

特徵方程與它的判別式

那個二次式 a r^2 + b r + c = 0 就是微分方程的特徵方程(有些書稱它為輔助方程——同一個東西)。它承載了關於系統無外力行為的全部資訊。這裡有一本乾淨的字典在運作:方程的係數逐項變成二次式的係數。一旦你信任這個樣式,就能用觀察法直接寫下特徵方程——把 y'' 換成 r^2、y' 換成 r、y 換成 1——完全不必代入。

a y'' + b y' + c y = 0     ->     a r^2 + b r + c = 0

  roots:  r = ( -b  +/-  sqrt(b^2 - 4ac) ) / (2a)

  discriminant  D = b^2 - 4ac  decides everything:
     D > 0   two distinct real roots r1, r2
     D = 0   one repeated real root  r
     D < 0   complex conjugate pair  alpha +/- i*beta
正是你一直在用的同一條二次公式;判別式 b^2 - 4ac 的正負號,把每個問題恰好分進三種情形之一。

二次公式給出 r = (-b +/- sqrt(b^2 - 4ac)) / (2a),而根號下那一個量,也就是判別式 D = b^2 - 4ac,主宰了整個故事。它可以是正、是零,或是負,而這三種可能不是無關緊要的細微差異——它們產生的是「種類」上真正不同的運動。正的判別式意味純粹的衰減或增長;零則處在刀刃之上;負則帶來振盪。接下來三節,逐一處理它們。

情形一——兩個相異實根

當 D > 0 時,二次式有兩個不同的實根 r1 與 r2。每一個都餵回猜測:y1 = e^(r1 t) 與 y2 = e^(r2 t) 都是解。它們獨立嗎?是的——兩個增長率不同的指數永遠不可能互為常數倍,所以這就是相異實根情形,而這一對 {e^(r1 t), e^(r2 t)} 構成一組基本解集。通解於是是 y = C1 e^(r1 t) + C2 e^(r2 t)。

用 y'' - 5 y' + 6 y = 0 把它具體化。它的特徵方程是 r^2 - 5 r + 6 = 0,因式分解為 (r - 2)(r - 3) = 0,所以 r1 = 2、r2 = 3。通解是 y = C1 e^(2t) + C2 e^(3t)。從物理上想像它:這是重度阻尼的情況,任何擾動都只是滑回靜止、從不越過它——那種緩緩闔上、而非用力甩動再彈回的過阻尼門弓。

情形二——重根,以及第二個解藏在哪裡

當 D = 0 時,兩個根碰撞合而為一,成為單一的重根值 r = -b/(2a)。現在指數把戲只交給你「一個」解,y1 = e^(rt)。但解空間是二維的——你仍然少了一個獨立的解,而且你不能只是換個常數重用 e^(rt),因為 C1 e^(rt) + C2 e^(rt) 不過是 (C1 + C2) e^(rt),一個假裝成兩個的單參數族。這就是真正棘手的重根情形

那個缺席的夥伴,原來是 y2 = t e^(rt)——同樣的指數乘上 t。你可以用代入法驗證它確實解出方程,但它之所以出現的更深理由,是本層緊接著下一篇指南的方法,降階法:把你手上僅有的一個解餵進去,它就藉由系統化的步驟 y2 = v(t) y1,解一個較簡單的關於 v 的方程,製造出第二個解。對重根而言,那個程序每次都交出 v = t,這正是為什麼那個 t 因子不是魔術,而是一套誠實方法可預測的產物。

所以當 r 是重根時,基本解集是 {e^(rt), t e^(rt)},通解是 y = (C1 + C2 t) e^(rt)。具體地,y'' - 4 y' + 4 y = 0 給出 r^2 - 4 r + 4 = (r - 2)^2 = 0,重根 r = 2,所以 y = (C1 + C2 t) e^(2t)。物理上這是臨界阻尼情形:在不過衝的前提下盡可能最快地回到靜止,正是介於過阻尼那種遲緩滑回,與情形三即將釋放的振盪之間的那道刀鋒。

情形三——複根,先睹為快

當 D < 0 時,負數的平方根迫使根成為一對複共軛,r = alpha +/- i*beta,二次式根本沒有實數解。形式上指數把戲仍然管用,產生 e^((alpha + i*beta) t)——但那是一個複數值的函數,而原方程的係數是實數,所以我們預期的是實數運動。這就是複共軛根情形,也正是振盪的來源:物理上,D < 0 是欠阻尼的情況,那座鳴響的鐘、那個擺盪的單擺。

從複指數通回誠實實數解的橋樑,是歐拉公式 e^(i*beta*t) = cos(beta t) + i sin(beta t),它把虛指數換成正弦與餘弦。其結果——在下一篇指南中完整展開——是實數基本解集 {e^(alpha t) cos(beta t), e^(alpha t) sin(beta t)}:一個以頻率 beta 振盪、被包絡 e^(alpha t) 包住的運動,當 alpha < 0 時包絡衰減。我們在這裡刻意停在門口——情形三值得擁有它自己的一篇指南。

配方、限制,與誠實的細則

  1. 用觀察法寫出特徵方程:把 y'' 換成 r^2、y' 換成 r、y 換成 1,得到 a r^2 + b r + c = 0。
  2. 計算判別式 D = b^2 - 4ac 並讀出情形:D > 0 相異實根、D = 0 重根、D < 0 複根。
  3. 建立基本解集:{e^(r1 t), e^(r2 t)}、或 {e^(rt), t e^(rt)}、或 {e^(alpha t) cos(beta t), e^(alpha t) sin(beta t)}。
  4. 用兩個常數組合出通解;若給了初始條件,解出 C1 與 C2 以釘住特解。

對這個方法的邊界要誠實。整個把戲建立在一個假設上:係數 a、b、c 是「常數」。一旦某個係數依賴於 x——例如 x^2 y'' + x y' + y = 0——代入 e^(rt) 就不再乾淨地提出公因子,特徵方程蒸發,你就需要別的工具(柯西-歐拉代換,或冪級數)。三情形的圖像之所以如此漂亮,正因為它特殊;它不是二階方程的通用求解器。

有兩個值得帶往後面的細點。其一,同樣的想法可以放大:一個 n 階常係數方程有一個 n 次的特徵多項式,而一個重複 m 次的根貢獻出族 e^(rt)、t e^(rt)、……、t^(m-1) e^(rt)——情形二的 t 因子樣式得到推廣。其二,別把「判別式為零」與「沒有阻尼」混為一談:D = 0 是臨界阻尼的刀鋒,而真正的無阻尼振盪(b = 0)則穩穩落在情形三,帶有純虛根。把情形分類的,是判別式,而非單獨的阻尼係數。