两个家族,两个物理问题
在本阶梯前面,你已经学会了那台引擎:当一个微分方程拒绝给出闭式解时,你猜一个幂级数,代进去,让方程碾出一条递推关系,一个接一个地定下系数——这就是幂级数法与弗罗贝尼乌斯法。多数时候这台机器产出一个无名的无穷级数。但每隔一阵就会发生一件美妙的事:对参数的某些特殊取值,递推会突然停下,无穷级数在有限多项后啪地合上,掉出一个多项式。整个物理学里最重要的两个多项式家族,正是这样诞生的。
第一个是厄米家族。把一个粒子丢进一个尽可能温柔的陷阱——一个把它往中心拉回、拉力正比于它偏离多远的力,即量子谐振子——支配它允许态的方程便约化为 y'' - 2x y' + 2n y = 0,即厄米方程。级数解恰好在 2n 为偶整数时合成一个多项式,而那个整数就是能级:能量量子化的最初登场,无非就是「级数必须终止,否则波函数会炸掉」。最低的几个是 H_0 = 1,H_1 = 2x,H_2 = 4x^2 - 2,H_3 = 8x^3 - 12x。每个 H_n 都是次数恰为 n 的多项式,并且恰有 n 个实零点——也就是那个振子态绝不可能被找到的 n 个位置。
第二个是拉盖尔家族。现在把一个电子束缚到一个质子上——氢原子——在你用勒让德多项式剥掉角向部分之后,径向部分服从 x y'' + (1 - x) y' + n y = 0,即拉盖尔方程。同样地,级数只在 n 为整数时终止,而那个整数同样是一个量子数;拉盖尔多项式 L_n(x) 描述电子的概率如何沿半径铺开,其中 L_0 = 1,L_1 = 1 - x,L_2 = (x^2 - 4x + 2)/2。于是每门量子力学入门课的两块基石——振子与氢原子——在数学上是同一个故事里的两位成员:一个被逼着终止的弗罗贝尼乌斯级数,留下一个干净的多项式。
罗德里格斯公式:一道配方造出整个家族
靠递推一项一项地碾出每个多项式当然行得通,但既繁琐,又对整个家族什么都说不出。有一把远为优雅的钥匙,能一次打开每一位成员:罗德里格斯公式。这些家族每一个都有一条——一个单一的表达式,你取某个简单的基底函数,把它求导 n 次,第 n 个多项式连同归一化就一起出来了。对厄米,它写作 H_n(x) = (-1)^n e^{x^2} (d^n/dx^n) e^{-x^2}:从高斯鼓包 e^{-x^2} 出发,对它施加第 n 阶导数,再乘以 e^{x^2} 把多出来的指数抵消掉,剩下的就是那个次数为 n 的纯多项式。
RODRIGUES FORMULAS (one differentiation recipe per family)
Hermite H_n(x) = (-1)^n e^{x^2} d^n/dx^n [ e^{-x^2} ]
Laguerre L_n(x) = (1/n!) e^{x} d^n/dx^n [ x^n e^{-x} ]
Legendre P_n(x) = (1/(2^n n!)) d^n/dx^n [ (x^2 - 1)^n ]
WORKED: Hermite, n = 2
e^{-x^2} -> base
d/dx e^{-x^2} = -2x e^{-x^2}
d^2/dx^2 e^{-x^2} = (4x^2 - 2) e^{-x^2}
H_2 = (+1) e^{x^2} (4x^2 - 2) e^{-x^2} = 4x^2 - 2 [matches]
The pattern: base function = (weight w(x)) times (a simple factor),
and differentiating n times pumps out the degree-n polynomial.看一看这三条罗德里格斯公式共有的东西,你就看见了那个深层的模式。每一条里你所求导的基底函数,都是该家族的权函数乘以一个简单的代数因子:厄米是 e^{-x^2},拉盖尔是 x^n e^{-x},勒让德是 (x^2 - 1)^n。那个权不是装饰——它是整个家族的指纹,而在下一节里,它正是让这些多项式正交的东西。罗德里格斯公式就是那座桥:它取来权函数——那个编码了问题所居之处的物理的东西——并从中制造出多项式,一次一个干净的高阶导数。
生成函数:以多项式为系数的幂级数
这是第二把万能钥匙,可以说还是更魔幻的那一把:生成函数。它的想法是一个记账员的把戏。与其拎着一整张无穷长的多项式清单 H_0, H_1, H_2, ...,不如把它们全部塞进一个普通的双变量函数的系数里——一个变量是 x,外加一个记账变量 t——再让一个老老实实的、对 t 的泰勒展开在你需要时把它们交还给你。对厄米,整个家族就住在一个利落的指数里:e^{2xt - t^2} = sum over n of H_n(x) t^n / n!。把左边展成对 t 的幂级数,读出 t^n 的系数,乘以 n!,那就是 H_n(x)——家族里的每一个多项式,都藏在这一个闭式函数里。
把整个家族挤进一个函数,妙处在于你现在能整批地对家族做运算。把厄米生成函数对 t 求导,你得到相邻系数之间的一条关系;对 x 求导,又得到另一条。把这两条关系翻译回系数的语言,你就白白拿到了家族的递推关系:H_{n+1} = 2x H_n - 2n H_{n-1},它让你从任意两个相邻多项式爬到下一个;以及 H_n'(x) = 2n H_{n-1},它说求导不过是让你顺梯子往下走一级。那些一个多项式一个多项式地证要写好几页的性质,这里一两行就掉了出来,因为你是对所有 n 同时证好了的。
带权正交:为什么这些才是对的基
现在轮到那条让这些多项式不可或缺、而不只是可爱的性质。它们是正交的——但是是在一种微妙的、带权的意义下正交,这一点你必须分毫不差地弄对。两个不同的厄米多项式,在它们乘积的普通积分下并不正交;你必须先插进家族的权函数 e^{-x^2}。具体说,H_m(x) H_n(x) e^{-x^2} 从负无穷到正无穷的积分,只要 m 不等于 n 就为零。这个权不是事后补上的——它正是当初坐在罗德里格斯基底里的那个 e^{-x^2},也正是首先让这个积分能在整条无穷直线上收敛的那个因子。
每个家族都带着自己的权和自己的定居区间,而且这种配对是刚性的。厄米住在整条实轴上,权为 e^{-x^2};拉盖尔住在从 0 到无穷的半轴上,权为 e^{-x};勒让德住在从 -1 到 1 的有限区间上,权干脆就是 1;切比雪夫住在同一区间上,但权为 1/sqrt(1 - x^2)。这些都不是随意的联姻。每个权恰恰就是出现在对应量子波函数里的那个 e^{-x^2} 或 e^{-x}——振子的高斯包络、束缚电子的指数衰减——所以正交性正是由多项式所从来的那个物理本身搭起来的。
正交性为什么如此要紧?因为它把一个无穷的多项式家族变成了一套坐标系。正如空间里的一个向量沿着相互垂直的轴分解成独立的分量,直线上一个像样的函数可以展开成厄米多项式之和,而每个系数都由一个单独的带权积分求得——把你的函数乘以 H_n,配上权 e^{-x^2},积分,再除以已知的范数。不用解联立方程,模与模之间互不干扰:正交意味着每个分量都能独立地读出来。这与傅里叶级数是同一套逻辑,只不过由 H_n 扮演那里正弦余弦所扮演的角色,而它正是你在量子力学与信号分析里将不断用到的那些展开背后的引擎。
超几何模板:藏在它们全体背后的一个主方程
退后一步,一种惊人的统一便映入眼帘。厄米、拉盖尔、勒让德、切比雪夫、贝塞尔——每一个原是一个独立的方程,各有自己的多项式、自己的权、自己那一隅物理。可它们其实并不独立。它们几乎全都是同一个主方程——超几何方程——的特例:x(1 - x) y'' + [c - (a + b + 1)x] y' - a b y = 0。它的解是超几何函数,记作 2F1(a, b; c; x)——一个幂级数,其系数遵循一条由仅仅三个旋钮(参数 a、b、c)支配的通用递推。把这三个旋钮拨到对的位置,主级数就一个接一个地坍缩到勒让德、坍缩到切比雪夫、坍缩到其余各位身上。一个方程,三个旋钮,一整座特殊函数的动物园。
厄米与拉盖尔在这套体系里坐在哪儿?它们由一位近亲支配,即合流超几何方程:x y'' + (c - x) y' - a y = 0,其解 1F1(a; c; x) 就是合流超几何函数。「合流」是一个精确的术语:它意味着完整超几何方程的两个奇点被推到一起,直到合并为一个——就像你让一个参数趋向某个极限时,一片地貌上两个不同的特征会汇流到一处。正是这种合并,把有界区间上的问题带进了半轴与全轴上的问题。拉盖尔是一个 a 取整数的 1F1;厄米也是由 1F1 搭起来的。于是振子与氢原子都住在这个合流模板里,是同一组旋钮的两种设定。
你真正该带走的,是一种视角的转变。从前你看到的是一摞互不相干、需要死记的命名函数,如今你能看见一棵家谱树:根上是超几何方程,朝厄米与拉盖尔生长的枝杈是合流方程,而那些命名的多项式则是叶子,每一片都由让级数终止的特殊参数取值拣选出来。同样那三套机器——终止成多项式、罗德里格斯公式、生成函数,全都由带权正交性编织在一起——在每一根枝条上反复出现。学会这个模板,你学到的就不是又多了一个特殊函数;你学到的,是它们全体所讲的那门语法。