從特徵值問題的動物園走向一個母形式
到此為止你已多次碰過特徵值問題：一個帶參數 lambda 的二階方程，加上兩端的邊界條件，唯有在一串特殊而離散的 lambda 值上才有非零解。上一篇指南中，對於 [0, L] 上、邊界條件 y(0) = y(L) = 0 的 y'' + lambda y = 0，你求得特徵值 lambda_n = (n pi / L)^2 與特徵函數 sin(n pi x / L)。很美——但它可能讓人覺得只是那一條整潔方程的僥倖巧合。斯圖姆－劉維爾理論的存在，正是要證明這絕非巧合。
關鍵的一步，是不再盯著每條方程那些偶然的係數，而是把每一條二階方程都逼進同一個標準外形。斯圖姆－劉維爾形式把特徵值問題寫成 (p(x) y')' + q(x) y + lambda w(x) y = 0，或等價地 -(p y')' + q y = lambda w y。慢慢讀：最高階導數那一項被包進一個外層導數 (p y')' 之內，而非裸露成 p y''。這個包裹並非裝飾——它正是後面那些優美性質得以成立的全部理由。這裡 p(x) > 0 與 w(x) > 0 是給定的正函數，而 w 稱為權重。
任何二階方程都能化成這個形式
這是讓人安心的部分：你不必被交到手上的就已經是斯圖姆－劉維爾形式。「任何」線性二階方程 a(x) y'' + b(x) y' + c(x) y + lambda d(x) y = 0，都能藉著乘上某個精心挑選的函數而被揉成這個形式。這項手法和你在一階方程裡用過的積分因子是同一個念頭，只是現在它的任務是讓前兩項塌縮成一個恰好的導數 (p y')'。
具體地說：從 y'' + (b/a) y' + ... 出發，把整條方程乘上 mu(x) = (1/a) exp(integral of (b/a) dx)。乘完之後，領頭那部分 mu y'' + mu (b/a) y' 便恰好成為 (mu y')'——乘積法則倒著跑。於是 p(x) = mu(x)，而權重 w(x) 就是經同樣乘法後坐在 lambda 前面的那個東西。只要 a(x) 在區間上不歸零，這套配方就永不失手；這正是我們能誠實地說：斯圖姆－劉維爾形式不是特例，而是每個正則二階特徵值問題都能披上的一件普世外衣。
start: a(x) y'' + b(x) y' + c(x) y + lambda d(x) y = 0
multiply by mu(x) = (1/a) * exp( integral (b/a) dx )
result: ( p y' )' + q y + lambda w y = 0
p = mu w = mu * d / a ... (the weight)
q = mu * c / a
example: Legendre (1 - x^2) y'' - 2x y' + lambda y = 0
b/a = -2x/(1-x^2), exp(integral) = (1 - x^2)
=> ( (1 - x^2) y' )' + lambda y = 0 p = 1 - x^2, w = 1對稱性，以及正交性為何被逼出來
現在來兌現。定義算子 L[y] = -(p y')' + q y，於是問題寫成 L[y] = lambda w y。整座大廈靠著一條恆等式——拉格朗日恆等式——它說：對任意兩個函數 u 與 v，組合 u L[v] - v L[u] 本身是一個完美的導數：它等於 -d/dx[ p (u v' - u' v) ]。把兩邊在區間 [a, b] 上積分：右邊便成為一個純粹的邊界項，只在兩端取值。
重點來了。對於正則斯圖姆－劉維爾問題的標準邊界條件——固定端、自由端、或週期端——那個邊界項恰好歸零。剩下的便是 integral of (u L[v] - v L[u]) dx = 0，這正說明 L 是自伴的：它能在積分內從一個位置跳到另一個位置而毫無殘餘，恰如對稱矩陣滿足 u^T A v = v^T A u。一切美好之事如今都從這唯一的方程流出，正如對稱矩陣的每個性質都源自 A = A^T。
把自伴性用到兩個特徵值不同（lambda_m 與 lambda_n）的特徵函數 y_m 與 y_n 上。兩行計算便得到 (lambda_m - lambda_n) 乘以 integral of w y_m y_n dx = 0。既然兩特徵值相異，那個積分本身就必須為零。這就是帶權重的正交性：只要 m 不等於 n，便有 integral of w(x) y_m(x) y_n(x) dx = 0。特徵函數彼此垂直——不是在尋常意義下，而是在以 w 加權的內積之下。這是整套理論的核心寶藏，也正是它使下一篇指南的特徵函數展開成為可能。
那些保證：實的、有序的、完備的
自伴性交給你一整包保證，全都是一勞永逸地證成，而非一例一例地重新發現。其一，每個特徵值都是「實」的——不會有意外的複數頻率溜進來，這正是你對一個物理振動或溫度所期望的。其二，特徵值構成一條無窮遞增的數列 lambda_1 < lambda_2 < lambda_3 < ...，一路奔向無窮，中間既不堆積、也不漏失。確有一個貨真價實的「第一個」特徵值，其上是一道乾淨的階梯。
其三——也是最深的一份禮物——是完備性。特徵函數不只是正交的；它們是一組完備的集合，意思是 [a, b] 上「任何」合理的函數都能寫成它們的無窮組合 sum c_n y_n(x)。這正是傅立葉級數底下那具抽象引擎：正弦與餘弦之所以管用，是因為它們是某一個特定斯圖姆－劉維爾問題的特徵函數，而它們之所以張滿整個空間，是因為理論保證每一個這樣的族都如此。上一節的正交性，恰恰是讓你能用單一個積分算出每個係數 c_n、而不必去解聯立方程的關鍵。
還有一個值得點名的結構性附贈：特徵函數 y_n 在內部恰有 n - 1 個零點（斯圖姆振盪定理）。第一個特徵函數在區間內部從不穿過零，第二個穿過一次，第三個兩次，依此類推——特徵值越高，形狀就越扭動。在物理上這就是泛音列：一根弦的基本模態是單一道拱，下一個諧波有一個節點，再下一個有兩個。數學早已悄悄把振動的圖像編碼進去，遠在你動手畫出它之前。
細則：正則、奇異、以及可能出錯之處
對這些條件要誠實，因為那些乾淨的保證並非無條件。完整的那一包——實的、單重的、遞增的特徵值，以及正交的、完備的特徵函數——是針對正則斯圖姆－劉維爾問題證成的：一個有限區間 [a, b]，p(x) > 0 與 w(x) > 0 在「整個」閉區間（含端點）上都成立，且兩端各有「分離型」邊界條件。一旦違反其中任何一條，你所引用的那條定理，可能就不再是適用的那一條了。
最重要的現實案例，刻意地「不」正則——它們是奇異斯圖姆－劉維爾問題。貝索方程的 p(x) = x 在端點 x = 0 處歸零；勒讓德方程的 p(x) = 1 - x^2 在 x = +/-1 兩端都歸零。當 p 在某端觸零、或區間無窮時，那條整潔的邊界條件便被一個較柔軟的要求取代——通常只要求解在那個麻煩的端點處保持「有界」即可。值得注意的是，大部分良好行為依舊存活：你仍得到正交性（配上正確的權重）與一個可用的展開，但譜的性質可能改變——有時特徵值會成為一段連續區間而非離散清單，正如全直線上的傅立葉變換。
最後一個警告，提防把理論讀過了頭。正交性與完備性對於以「封閉形式」找出特徵函數一事，什麼也沒說——對大多數斯圖姆－劉維爾問題並無初等公式，正如大多數常微分方程沒有封閉形式的解，而那些特徵函數是由問題本身所定義的特殊函數（貝索、勒讓德、厄米特……）。這套理論是關於「結構」的保證，不是通往公式的捷徑。它所承諾的是：那份結構——實的譜、帶權重的正交性、完備性——永遠在那兒任你取用，無論你寫不寫得出特徵函數。而那份結構，正是下一篇指南要化為一套可運作之展開的東西。