勒让德方程与勒让德多项式

这个方程从何而来

在本阶梯前面的指南里，你学会了用一招对付没有闭形式解的常微分方程：把解写成一个幂级数，再一项一项地把系数磨出来。勒让德方程正是这套方法的招牌问题，而它的盛名当之无愧：每当一个问题落在球面上，你就躲不开这个方程。把它写成 (1 - x^2) y'' - 2x y' + l(l+1) y = 0，其中 l 是一个可由你拨动的常数。那个看起来怪怪的组合 l(l+1) 并非偶然——正如我们将看到的，它恰恰是从几何中掉出来的形状。

下面是它如何出现的简短版本。取拉普拉斯方程 nabla^2 V = 0，这是静电学与引力的主方程，把它写到球坐标里。然后用分离变量法——把未知函数拆成径向部分、极角部分和方位角部分的乘积。径向与方位角部分都能轻松解出；而极角部分，在代换 x = cos(theta) 之后，恰好就是勒让德方程。所以这里的 x 不是长度——它是从北极量起的夹角的余弦，取值落在从 -1 到 +1 的区间上。

用级数求解——以及一桩幸运的意外

x = 0 是这个方程的常点（那里没有任何东西发散），所以朴素的幂级数法就适用——不必动用你为奇点所学的那套更笨重的弗罗贝尼乌斯机器。代入 y = sum of a_k x^k，塞进 (1 - x^2) y'' - 2x y' + l(l+1) y = 0，再按 x 的幂次归并。对齐下标之后，你会得到一个把 a_{k+2} 与 a_k 相联系的两项递推关系：a_{k+2} = a_k 乘 [k(k+1) - l(l+1)] / [(k+1)(k+2)]。自由地选取 a_0 与 a_1，就给出任何二阶方程都该有的两个独立解。

现在盯紧那个递推关系，看好戏上演。分子是 k(k+1) - l(l+1)。当 k 恰好等于 l 的那一刻，这个因子正好为零——于是 a_{l+2} = 0，接着 a_{l+4} = 0，整条级数的尾巴全部坍缩为零。如果 l 是非负整数，两条级数中就有一条干脆终止：它不再是无穷级数，而成了一个货真价实的、次数为 l 的多项式。这正是本主题核心处那桩幸运的意外。把每个这样的多项式归一化，使它在 x = 1 处取值为 1，你就得到了勒让德多项式 P_l(x)。

Recurrence:  a_{k+2} = a_k * [ k(k+1) - l(l+1) ] / [ (k+1)(k+2) ]

Numerator vanishes when k = l  =>  series stops at degree l.

First few Legendre polynomials (normalized P_l(1) = 1):
  P_0(x) = 1
  P_1(x) = x
  P_2(x) = (1/2)(3 x^2 - 1)
  P_3(x) = (1/2)(5 x^3 - 3 x)
  P_4(x) = (1/8)(35 x^4 - 30 x^2 + 3)

递推关系在 l 次处把级数截断；由此得到、并归一化为 P_l(1) = 1 的多项式，就是勒让德多项式。

不过对另一个解也要诚实。对每个整数 l，两条级数里只有一条会终止；与之配对的解（记作 Q_l，即第二类勒让德函数）仍是无穷级数，并在 x = +1 与 x = -1 处发散——恰恰是北极和南极。物理把它丢弃，因为真实的势在那里必须保持有限。所以多项式并非全部故事；它们只是经受住有限性要求而存活下来的那一半。

同一对象的三副面孔

一个特殊函数一旦安顿下来，数学家们就会收集召唤它的不同方式，而每一副面孔都擅长一类不同的活儿。我们刚搭好的「级数加递推」这副面孔，适合用来逐项地证明命题。第二副面孔是罗德里格斯公式，一个简洁的单一表达式：P_l(x) = (1 / (2^l l!)) 乘以 (x^2 - 1)^l 的第 l 阶导数。把 (x^2 - 1)^l 恰好求导 l 次，再乘上那个常数，整个多项式就完整地蹦出来了——不需要任何递推记账。

第三副面孔是一个生成函数，也是物理学家私下里最钟爱的那一副：1 / sqrt(1 - 2 x t + t^2) = sum over l of P_l(x) t^l。把所有勒让德多项式打包进一个关于 t 的幂级数的系数里。这不是一个乔装成公式的巧合——那个左端正是出现在每个静电学问题里的倒数距离 1/|r - r'|，按距离比的幂次展开的结果。所以这个生成函数实际上就是在说「点电荷的势，用球谐展开写出来」。三副面孔——级数、罗德里格斯、生成函数——三者描述的都是同一族 P_l(x)。

正交性：让它们好用的那条性质

关于这些多项式，最重要的一个事实是：它们在从 -1 到 +1 的区间上彼此正交。具体地说：从 -1 到 1 对 P_m(x) P_n(x) dx 的定积分，在 m 与 n 不同时为零，在 m = n 时等于 2/(2n+1)。把每个 P_l 想成函数无穷维空间里的一个「方向」，把这个积分想成点积；正交性是在说这些方向彼此垂直，就像 x、y、z 三个坐标轴。这种垂直不是运气——它之所以成立，是因为勒让德方程是一个斯图姆–刘维尔问题，而这类问题总会交给你一族正交的解。

正交性为什么这么要紧？因为它让你能把 -1 到 1 上任意一个像样的函数 f(x)，写成一个勒让德级数：f(x) = sum over l of c_l P_l(x)。要求出某个特定的系数 c_n，你不必去解一个庞大的线性方程组——你只要做投影，正如你靠把向量与 x 轴作点积来读出它的 x 分量一样。把两边同乘 P_n(x)，从 -1 到 1 积分，右端每一项都因正交性而消失，只剩第 n 项。存活下来的便给出 c_n = [(2n+1)/2] 乘以从 -1 到 1 对 f(x) P_n(x) dx 的积分。每个系数一个积分，各自独立算出。这与傅里叶级数是同一个思想——正弦余弦是一族正交函数，勒让德多项式则是适配于球面的另一族。

它们究竟出现在哪里

从静电学说起，那是它的主场。假设你知道一个球面上的电压随纬度的变化，并想求出球外各处的势。一旦手里有了勒让德多项式，配方就是机械的。把边界数据按 cos(theta) 展开成勒让德级数；给每个 P_l(cos theta) 配上分离后的径向方程所提供的径向解 r^l（球内）或 1/r^{l+1}（球外）；再把它们加起来。l = 0 项是单极子（总电荷），l = 1 是偶极子，l = 2 是四极子——著名的多极展开，无非就是一个关于角度的勒让德级数，再附上那些径向幂次。

在球坐标下分离拉普拉斯方程；其中的极角因子，令 x = cos(theta)，满足勒让德方程。
要求解在两极 x = +1 与 x = -1 处保持有限；这迫使 l 取非负整数，并选出 P_l，丢弃发散的搭档 Q_l。
把边界数据展开成勒让德级数来匹配；用正交积分作投影，逐个读出系数。
给每个角向项 P_l(cos theta) 配上它的径向搭档（球内 r^l，球外 1/r^{l+1}）再求和——这就是完整的势。

再看更深一层的出场：氢原子。描述束缚在质子上的电子的薛定谔方程，又是一个球对称问题，于是它以同样的方式分离，极角因子又一次由勒让德方程支配（确切地说，一旦方位角指标进来，是它的连带表亲）。那个让级数终止的整数 l，恰恰就是轨道角动量的量子数，而组合 l(l+1)，至多差个常数，正是角动量平方算符的本征值。电子角动量之所以以整数步长量子化，原因与勒让德级数之所以必须终止是同一个：只有整数 l 才能让波函数在两极保持有限。你见过的那些标着 s、p、d、f 轨道的形状，就是 P_l 及其连带亲属的图像。

可以带走的模式

退一步，看清刚刚发生之事的形状，因为接下来的指南会几乎一字不差地用贝塞尔和埃尔米特重演它。一个带对称性的物理问题逼出某个特定的常微分方程。用级数去解它，会产生一个递推关系；一个量子化条件（此处是 l 取非负整数）让级数终止或保持良好行为；存活的解组成一族正交函数；而那种正交性把「在这组基里展开任意函数」变成了每个系数一个投影积分的例行公事。勒让德是这台机器最干净的头一个实例，因为它的解是货真价实的多项式。

留一个诚实的提醒让你脚踏实地：终止成多项式是整数 l 的特权。对非整数的 l，勒让德方程仍然有完全合理的解，但它们是真正的无穷级数（勒让德函数，而非多项式），通常在两极处行为不佳。你一路游历的那个干净的多项式世界之所以存在，恰恰是因为物理在一个闭合球面上要求有限性，从而选出了整数值——几何悄然地选择了算术。真正的教训不是多项式本身，而是这条选择原理。