5.2 Legendre 方程#

本章導言#

\[ (1 - x^2)\, y'' - 2x\, y' + \lambda\, y = 0. \]

這個方程誕生於十八世紀最優雅的數學物理問題：球面上的位勢。1780 年代，Adrien-Marie Legendre（勒讓德，1752–1833） 在研究行星形狀與引力吸引時系統分析了它，其多項式解 \(P_n(x)\) 成為球面幾何的「字母表」——今日從地球重力場模型到量子角動量、從電腦圖形的球諧函數到機器學習的球面卷積， \(P_n\) 無處不在。本節講述 Legendre 方程的級數解、多項式的誕生機制、正交性與其在球面上的身份。

方程的來源：球面上的 Laplace 方程#

在球座標 \((r, \theta, \phi)\) 下解 Laplace 方程 \(\nabla^2 V = 0\) （9.3 節），分離變數 \(V = R(r)\, \Theta(\theta)\, \Phi(\phi)\) 後，極角部分給出

\[ (1 - x^2)\, \Theta'' - 2x\, \Theta' + n(n+1)\, \Theta = 0, \qquad x = \cos\theta. \]

參數的特殊形式 \(\lambda = n(n+1)\) 不是巧合，而是角動量量子化的經典前身：只有這些值使解在 \(\theta = 0, \pi\) （兩極）正則。「物理上的正則性條件，數學上表現為參數的離散化」——這一機制在 1920 年代成為量子力學的核心原理。

歷史：Legendre 1784–1787 年的行星形狀論文與 1793 年的引力研究建立了 \(P_n\) ；同代 Laplace 在天體力學（《天體力學》1799–1825）中把它們用到極致，故文獻中亦見「Laplace 係數」之名。

級數解與多項式化#

\(x_0 = 0\) 是常點（3.3 節）。設 \(y = \sum a_n x^n\) ，代入（注意 \((1 - x^2) y''\) 項）：

\[ (1 - x^2) \sum n(n-1) a_n x^{n-2} - 2x \sum n a_n x^{n-1} + \lambda \sum a_n x^n = 0. \]

對齊 \(x^k\) （記 \(\lambda = n(n+1)\) ， \(n\) 暫為參數）得遞推：

\[ (k+2)(k+1)\, a_{k+2} = \big[k(k+1) - n(n+1)\big] a_k. \]

奇蹟發生在 \(k = n\) （ \(n\) 非負整數）：係數 \(k(k+1) - n(n+1) = 0\) ，遞推斷裂——級數在 \(x^n\) 截斷，成為多項式：

\[ P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n}\big(x^2 - 1\big)^n \qquad \text{（Rodrigues 公式，1816）}. \]

前幾個：

\[ P_0 = 1, \quad P_1 = x, \quad P_2 = \frac{3x^2 - 1}{2}, \quad P_3 = \frac{5x^3 - 3x}{2}, \quad P_4 = \frac{35x^4 - 30x^2 + 3}{8}. \]

不截斷的那條鏈成為第二類解 \(Q_n(x)\) ——在 \(x = \pm 1\) 對數發散（正則奇點，3.3 節），物理上（兩極）必須拋棄。物理條件剔除一半的數學解，多項式倖存並統治。

正交性：Legendre 多項式是球面的「座標軸」#

正交關係： \(\displaystyle\int_{-1}^{1} P_m(x)\, P_n(x)\, dx = 0 \quad (m \neq n)\) ； \(\displaystyle\int_{-1}^{1} P_n^2\, dx = \frac{2}{2n+1}\) .

證明思想：方程可寫為自伴形式（4.1 節場景二）

\[ \frac{d}{dx}\Big[(1 - x^2) P_n'\Big] + n(n+1) P_n = 0. \]

以 \(P_m\) 乘 \(P_n\) 的方程、 \(P_n\) 乘 \(P_m\) 的方程，相減積分：

\[ \big[n(n+1) - m(m+1)\big] \int P_m P_n\, dx = \Big[(1 - x^2)(P_m P_n' - P_n P_m')\Big]_{-1}^{1}. \]

邊界項為零（ \(1 - x^2 = 0\) ），故 \(m \neq n\) 時積分為零。∎

這就是 Sturm–Liouville 理論（4.1 節）的教科書演出：自伴形式 + 邊界條件 ⇒ 正交特徵函數族。 \(P_n\) 之於球面函數，正如 \(e^{inx}\) 之於週期函數——任何（足夠好的）函數可展開為 Legendre 級數：

\[ f(x) = \sum_{n \ge 0} c_n P_n(x), \qquad c_n = \frac{2n+1}{2} \int_{-1}^{1} f(x) P_n(x)\, dx. \]

（Fourier–Legendre 展開；理論地位與傅立葉級數平行，9.1 節。）

遞推關係與生成函數#

計算工具箱（由母函數 \(G(x, t) = (1 - 2xt + t^2)^{-1/2} = \sum P_n(x) t^n\) 對 \(t\) 求導比對而得）：

\[ (n+1) P_{n+1} = (2n+1)\, x\, P_n - n\, P_{n-1}. \]

三項遞推讓計算如流水線——Legendre 1785 年論文已在使用。「多項式族 + 遞推 + 正交」的三件套，是所有經典正交多項式（Hermite、Laguerre、Chebyshev、Jacobi）的家族特徵，統攝於超幾何方程（5.1 節）之下。

為何母函數長這樣： \((1 - 2xt + t^2)^{-1/2}\) 正是兩點距離倒數的展開——點 \(\mathbf{r}\) （ \(r, \theta\) ）與單位向量上點 \(\mathbf{r}'\) （ \(t\) ）的重力/靜電位勢：

\[ \frac{1}{|\mathbf{r} - \mathbf{r}'|} = \frac{1}{r'} \sum_{n \ge 0} \Big(\frac{r}{r'}\Big)^n P_n(\cos\theta) \qquad (r < r'). \]

** \(P_n\) 是「把 3D 位勢按角向分解」的基**——這一句話就是它在物理中的全部身份。

一個具體例題#

例題 5.2.1：把函數 \(f(x) = |x|\) （ \(-1 \le x \le 1\) ）展開為 Legendre 級數至前三項（ \(P_0, P_2\) ），並與直接多項式逼近比較。

解： \(|x|\) 是偶函數，故只含偶多項式 \(P_{2k}\) 。

係數 \(c_0 = \frac12 \int_{-1}^{1} |x| \cdot 1\, dx = \frac12 \cdot 1 = \frac12\) 。

\(c_2 = \frac52 \int_{-1}^{1} |x|\, \frac{3x^2 - 1}{2}\, dx\) 。被積函數偶， \(\int_{-1}^{1} = 2 \int_0^1\) ：

\[ \int_0^1 x \cdot \frac{3x^2 - 1}{2}\, dx = \frac12 \Big[\frac{3x^4}{4} - \frac{x^2}{2}\Big]_0^1 = \frac12 \Big(\frac34 - \frac12\Big) = \frac18. \]

故 \(\int_{-1}^{1} |x| P_2\, dx = 2 \cdot \frac18 = \frac14\) ， \(c_2 = \frac52 \cdot \frac14 = \frac58\) 。

二項近似：

\[ |x| \approx \frac12 P_0 + \frac58 P_2 = \frac12 + \frac58 \cdot \frac{3x^2 - 1}{2} = \frac18 + \frac{15}{16} x^2. \]

評論： - 這是最小平方意義下的最佳二次逼近（正交展開的通性——Fourier 係數使誤差平方積分最小，9.1 節同證）。 - 比較 Chebyshev 或 Taylor 逼近：Taylor 在 \(x = 0\) 精確但兩端失真大（ \(|x|\) 在 0 不可微，Taylor 根本不存在！）；Legendre 均勻「攤平」誤差。正交多項式逼近 = 全域最優，這是數值積分（Gauss–Legendre 求積，1814）的原理：以 \(n\) 個 Legendre 零點為節點的求積公式達到 \(2n\) 次代數精確度——正交性轉化為計算效率。

歷史啟示#

物理正則性挑選數學解： \(Q_n\) 在兩極發散，被物理剔除；留下的 \(P_n\) 恰好構成完備正交族。應用問題的邊界條件，是數學結構的建築師。
參數離散化的先聲： \(\lambda = n(n+1)\) 由「解在兩極正則」強制，是量子化（角動量 \(l(l+1)\hbar^2\) ）的十九世紀彩排。
逼近論的誕生：Legendre 級數的均勻攤平誤差性質，開啟了最佳逼近理論與 Gauss 求積——解方程的副產品成為計算數學的支柱。

與其他章節及領域的關聯#

5.1 節：級數解法與遞推的直接應用；超幾何家族成員。
5.3 節：Bessel 是柱面（圓柱對稱）的對應物；Legendre 是球面的對應物。
4.1 節：Sturm–Liouville 理論與自伴形式；正交性證明的模板。
9.1、9.3 節：分離變數解熱方程/Laplace 方程時，角向部分就是 \(P_n\) ；球諧函數 \(Y_{lm} = P_l^m(\cos\theta) e^{im\phi}\) 。
量子力學：角動量算子的本徵函數即 \(P_l\) ；氫原子的 \(p, d\) 軌道形狀由 Legendre 多項式決定。
大地測量與衛星：地球重力場的球諧展開（EGM 模型）以 \(P_n\) 為基——GPS 高層模型裡都是它。
數值分析：Gauss–Legendre 求積；譜方法（spectral methods）用 \(P_n\) 展開解 PDE，指數收斂（10.1 節的「高階兄弟」）。

練習#

由遞推關係生成 \(P_2, P_3\) （從 \(P_0 = 1, P_1 = x\) 出發），並用 Rodrigues 公式核對 \(P_3\) 。
驗證 \(P_2 = \frac{3x^2 - 1}{2}\) 滿足 \(\lambda = 6\) 的 Legendre 方程（直接代入）。
計算 \(\int_{-1}^{1} x\, P_3(x)\, dx\) （提示： \(x = P_1\) ，用正交性一步得到；再用直接積分核對）。
證明 \(P_n\) 在 \((-1, 1)\) 內恰有 \(n\) 個實根（提示：Rodrigues 公式 + \((x^2 - 1)^n\) 在 \(\pm 1\) 有 \(n\) 重零點 + Rolle 定理）。
思考題：地球重力場按 \(P_n\) 分解時， \(n = 0\) 項是球對稱（平均質量）， \(n = 2\) 主導項反映扁率（赤道隆起）。請從「偶函數、低階主導」的角度解釋為什麼旋轉扁球的引力場以 \(P_0, P_2\) 為主、奇數階幾乎消失。

進階研讀：Legendre 理論的五個深化場景#

場景一：連帶 Legendre 函數。 \(x = \cos\theta\) 且帶方位角依賴時，方程升級為

\[ (1 - x^2) y'' - 2x y' + \Big[l(l+1) - \frac{m^2}{1 - x^2}\Big] y = 0, \]

解為 \(P_l^m(x) = (1 - x^2)^{m/2} \frac{d^m}{dx^m} P_l(x)\) 。** \(m\) 是方位角節數， \(l\) 是總節數——球諧函數 \(Y_l^m\) 的兩個量子數在此成形。角動量的代數（升降算子）對應這些函數間的遞推——微分方程的遞推關係，在量子力學中化身為角動量代數**。

場景二：零點與數值求積。 \(P_n\) 的零點（練習 4）是 Gauss–Legendre 求積節點： \(\int_{-1}^{1} f \approx \sum w_i f(x_i)\) 以 \(n\) 點達 \(2n - 1\) 次精確。證明核心：正交性使誤差泛函對所有 \(2n - 1\) 次多項式消失。Gauss 1814 年的發明至今是所有數值積分函式庫（QUADPACK、scipy）的核心。多項式的根，是積分的最優座標。

場景三：譜方法。把 PDE 的解展開為 \(P_n\) 級數、在譜空間做微分（利用遞推把導數變為係數矩陣），得到譜方法——光滑問題上誤差指數衰減，遠勝有限差分（10.1 節的多項式衰減）。氣象數值模式（球面上的 Navier–Stokes）傳統上以球諧函數為基底——你每天看的天氣預報，內部全是 Legendre 變換。

場景四：Legendre 函數的非整數階。 \(\lambda\) 非整數時解 \(P_\nu(x), Q_\nu(x)\) 仍是良定義函數（超幾何參數 \(\nu, -\nu - 1\) ），在勢論（旋轉橢球的引力）、流體力學（球外流場）中出現。參數連續化——從「多項式家族」到「函數家族」——是特殊函數論的標準動作，也是本徵值離散與連續譜的分界示範。

場景五：球面上的調和分析。 \(P_l^m(\cos\theta) e^{im\phi}\) （球諧函數）構成球面平方可積函數的完備正交基；球面上的「傅立葉分析」由此建立（Herglotz、Funaro 等 20 世紀標準化）。地球物理（重力、地磁、宇宙微波背景的多極展開 \(a_{lm}\) ——2019 諾貝爾獎級的宇宙學數據分析基礎）都以它為語言。Legendre 多項式是球面文明的通用文字。

補充：Legendre 的人生速寫#

Adrien-Marie Legendre（1752–1833）出身巴黎富裕家庭，大革命中失去財產，靠教職與科學院職位度日。貢獻：數論（二次互反律的猜想者，Gauss 證明）、橢圓積分（《橢圓函數論》1825–1832，後被 Abel 與 Jacobi 的反演革命超越——Legendre 晚年信中對新理論的優雅表示折服，是科學史上優雅的謝幕）、最小平方法（1805，與 Gauss 的優先權之爭）、本節的多項式。他站在多場革命的隔壁：互反律由 Gauss 證明、橢圓函數由 Abel–Jacobi 反演、最小平方法優先權與 Gauss 爭執——Legendre 是「打開門的人」的典型。

補充：Rodrigues 公式的證明思想#

\(P_n = \frac{1}{2^n n!} \frac{d^n}{dx^n}(x^2 - 1)^n\) 。驗證它滿足方程：令 \(u = (x^2 - 1)^n\) ，則 \((1 - x^2) u' + 2nx u = 0\) （直接驗）。對此式求導 \(n + 1\) 次（Leibniz 公式），出現 \(v = u^{(n)} = 2^n n! P_n\) ：

\[ (1 - x^2) v'' - 2x v' + n(n+1) v = 0. \]

「對一個容易的方程求導 \(n\) 次，得到難的方程的解」——微分恆等式的高階槓桿，這是 Rodrigues（1816）與後來 Rodrigues 型公式（Hermite、Laguerre）的共同套路。

名詞速查表（本節出現的核心術語）#

Legendre 方程： \((1 - x^2) y'' - 2x y' + n(n+1) y = 0\) 。
** \(P_n\) **： \(n\) 次 Legendre 多項式；級數截斷的產物。
** \(Q_n\) **：第二類解；在 \(\pm 1\) 對數發散，物理剔除。
Rodrigues 公式： \(P_n = \frac{1}{2^n n!} \frac{d^n}{dx^n}(x^2 - 1)^n\) 。
正交關係： \(\int_{-1}^1 P_m P_n = \frac{2}{2n+1} \delta_{mn}\) 。
母（生成）函數： \((1 - 2xt + t^2)^{-1/2}\) ；兩點位勢的展開。
三項遞推： \((n+1) P_{n+1} = (2n+1) x P_n - n P_{n-1}\) 。
連帶 Legendre 函數： \(P_l^m\) ；球諧函數的成分。
Gauss–Legendre 求積：以 \(P_n\) 零點為節點的高精度數值積分。
球諧函數： \(Y_l^m = P_l^m(\cos\theta) e^{im\phi}\) ；球面上的正交基。

年表速查（本節涉及的關鍵年份）#

1784–1787：Legendre 行星形狀論文， \(P_n\) 問世。
1793：Legendre 引力吸引研究。
1799–1825：Laplace《天體力學》，球面位勢理論。
1805：Legendre 最小平方法。
1814：Gauss 求積。
1816：Rodrigues 公式。
1825–1832：Legendre 橢圓函數論。
1926：量子力學角動量理論（ \(P_l\) 成為本徵函數）。
1970s：譜方法與球諧模式在天氣預報普及。

深化問答三則#

問：為什麼 \(\lambda\) 必須是 \(n(n+1)\) 形式才能得到多項式？ 答：遞推 \(a_{k+2} = \frac{k(k+1) - \lambda}{(k+2)(k+1)} a_k\) 的截斷需要分子在偶數或奇數鏈的某處歸零： \(\lambda = n(n+1)\) 使 \(k = n\) 處係數為零。參數的本徵值由「級數截斷」條件決定——這與 Sturm–Liouville 理論（特徵值離散化）完全一致；量子力學的能量量子化是同一機制的物理顯影。

問： \(Q_n\) 對數發散「壞」在哪裡？ 答：在球面問題中， \(x = \pm 1\) 對應兩極——單點而非邊界。位勢在單點發散意味著那裡有（無限集中的）源，違反「極軸上無特殊物質」的假設。數學上， \(Q_n\) 的對數項也破壞球諧函數的完備性（非平方可積）。「剔除發散解」不是數學潔癖，而是「空間沒有無限細的針」的物理承諾。

問：Legendre 級數與傅立葉級數是什麼關係？ 答：同一思想的兩個座標系。 \(x = \cos\theta\) 換元後， \(P_n(\cos\theta)\) 近似 \(\cos\) 函數的變形（加權正交：權重 \(dx = -\sin\theta\, d\theta\) ），區間 \([-1, 1]\) 變為 \([0, \pi]\) 。Legendre 展開是「球面經線方向」的傅立葉分析；加上 \(e^{im\phi}\) 就是完整球諧展開。一切正交展開都是一個思想：用特徵函數當座標（4.1 節場景二、9.1 節）。

（本節完）