一架命名曲线的参考书架
你早就懂得如何在理解一个函数之前先认识它。在任何人给你看证明之前很久,你就认识正弦了:一道升到 1、落到 -1、在 pi 的整数倍处穿过零点、永远重复的波。你能想象它、勾画它、用它推理。直到后来,泰勒级数和微分方程 y'' = -y 才解释了它为何如此。本篇要为一架著名的曲线书架做同样的事——正交多项式与贝塞尔函数——它们出现在高等物理与工程所到之处。我们现在把它们当作有已知图像和几条鲜明事实的命名对象来认识;以后的台阶才会从它们的微分方程把它们推导出来。
下面这个统一的想法让整架书架连成一体,在任何名字出现之前就值得直白地说清楚。你在线性代数里学过用互相垂直的坐标轴方向来搭一个向量,因为互相垂直的分量彼此不干扰:每个坐标都能单独读出。同样的把戏对函数也奏效,只要我们决定「垂直」对两条曲线意味着什么。答案是:把它们相乘、积分,当这个积分为零时就称它们垂直——正交。这架书架上的每个家族,不过是一组多项式(对贝塞尔来说是一组振荡曲线),它们被精心挑选得彼此正交。它们就是函数空间的垂直坐标轴。
四个多项式家族,按它们的栖息地
区分这四个家族最干净的方式,不是看它们的公式,而是看每一个住在哪里——它的区间和它的权。把它们想成适应了四种栖息地的四个物种。勒让德多项式住在有限区间 [-1, 1] 上,配最朴素的权 w(x) = 1;它们是球面的天然函数,所以才在球坐标下主宰物理。厄米多项式住在整条直线 (-infinity, infinity) 上,配钟形权 e^{-x^2}——正是前几篇里那个高斯形状;它们是量子谐振子的函数。参见勒让德多项式与厄米多项式。
拉盖尔多项式住在半直线 [0, infinity) 上,配衰减权 e^{-x};它们是氢原子的径向函数,那里电子的概率向外指数式地消退。切比雪夫多项式回到 [-1, 1] 上,但配权 1 over square root of (1 - x^2),这个权在两个端点附近堆起极大的强调。这种端点强调使它们成为数值逼近的冠军:它们把误差均匀地铺开,驯服了朴素高次拟合中那种狂野的振荡。参见拉盖尔多项式与切比雪夫多项式。四种栖息地、四种权、四个家族——而栖息地预示了用途。
Family Interval Weight w(x) Famous home
---------- --------------- ------------------ -----------------------
Legendre [-1, 1] 1 sphere / spherical coords
Chebyshev [-1, 1] 1/sqrt(1 - x^2) numerical approximation
Hermite (-inf, +inf) e^{-x^2} quantum harmonic oscillator
Laguerre [0, +inf) e^{-x} hydrogen atom (radial part)
First few Legendre: P0 = 1, P1 = x, P2 = (3x^2 - 1)/2, P3 = (5x^3 - 3x)/2
First few Chebyshev: T0 = 1, T1 = x, T2 = 2x^2 - 1, T3 = 4x^3 - 3x
Beautiful Chebyshev fact: T_n(cos theta) = cos(n theta) (a disguised cosine!)把图像想出来,它们就不再抽象。在 [-1, 1] 上,P_n 是一条起伏的曲线,归一化得使 P_n(1) = 1,它在区间内部恰好穿过零点 n 次——P_2 在轴下方下沉一次,P_3 摆过零点三次,依此类推;下标越大,起伏越多,正像弦上更高的谐波。切比雪夫的 T_n 更加生动:因为 T_n(cos theta) = cos(n theta),它的图像简直就是一道被弯到 [-1, 1] 上的余弦波,在 +1 与 -1 之间以等高的涟漪振荡。一旦你能勾画出它们,这些多项式就和它们所推广的三角函数一样具体了。
两台能生成任何家族的机器
你不必背这些多项式,也不必靠试错去找它们。有两台标准机器能把它们摇出来,现在把两台都见识一遍,神秘感就荡然无存了。第一台是一行求导公式。对勒让德它写作 P_n(x) = (1 over (2^n n!)) times (x^2 - 1)^n 的 n 阶导数。把简单的多项式 (x^2 - 1)^n 恰好求导 n 次,再缩放,P_n 就掉了出来——完美正交,无需运气。这就是一条罗德里格斯公式,书架上的每个家族都有一条(各有自己的内函数与权)。参见罗德里格斯公式。
第二台机器更漂亮:一个单独的封闭形式函数,它关于一个哑变量 t 的幂级数,把整个家族都装在系数里。对勒让德,把 1 over square root of (1 - 2 x t + t^2) 展开成关于 t 的幂级数,t^n 的系数恰好就是 P_n(x)。一个整洁的公式,悄悄存下了无穷多个多项式。这就是一个生成函数,它是一举证明家族性质的主力——把生成函数求导或相乘,整条整条的恒等式立刻掉出来。
为何正交性才是全部的要点
下面就是为整架书架正名的回报。因为家族是正交的,你可以把区间上几乎任何函数展开成家族成员之和——一个广义傅里叶级数——并通过单独一个积分读出每个系数,完全不必解代数方程组。这正是前面某个台阶上普通傅里叶级数的精确回响:那里正弦余弦彼此正交,每个系数都来自一个积分。它之所以奏效,靠的是那幅垂直坐标轴图景:当你向 P_k 投影时,其余每个成员都贡献为零,因为它与 P_k 的重叠积分等于零。
- 选出区间与权和你的问题相匹配的家族(球面 -> 勒让德,振子 -> 厄米,半直线衰减 -> 拉盖尔,逼近 -> 切比雪夫)。
- 要找到展开式中 P_k 的系数 c_k,就把你的函数与 P_k 带权相乘并积分,再除以 P_k 平方的已知积分(该家族的归一化)。一个积分,一个系数。
- 把 c_0 P_0 + c_1 P_1 + c_2 P_2 + ... 加起来,在项变小处截断。因为各部分正交,一个截断和就是该长度下最好的逼近——当你添加更多项时,没有哪个系数需要重新调整。
在你太过信任它之前,有两条诚实的提醒。第一,展开只在对函数加上条件时才好好收敛——粗略地说,函数应相对该权平方可积,并且和普通傅里叶级数一样,一处跳跃间断会在跳跃附近产生一个顽固的过冲,它不会随着你添项而消失(吉布斯现象)。第二,单凭正交性并不能保证家族丰富到足以表示每一个函数;那个额外的性质叫完备性,是一条你应当知道正在被假定的独立定理,而不是正交性免费奉送给你的东西。
贝塞尔函数:一道衰减的鼓面正弦
现在来认识这架书架上最著名的、却不是多项式的居民。往圆形池塘里丢一颗石子,或敲一面圆鼓面,涟漪会以一圈圈的环向外扩散。这种圆几何的天然形状,就是第一类贝塞尔函数,记作 J_n(x)。这幅心象精确而易记:J_n(x) 看上去像一道正弦波,永远振荡下去,但振幅缓慢缩小,大致按 1 over square root of x 衰减。它是一道在向外行进时耗尽能量的正弦——正是水面上一圈涟漪在它那一环的高度铺到越来越长的圆周上时所做的事。参见第一类贝塞尔函数。
几条鲜明的事实让 J_n 变得具体。J_0 在 x = 0 处从值 1 出发(像余弦),而每个下标 n 至少为 1 的 J_n 都从 0 出发(像正弦);之后它们全都振荡并衰减。它们的零点——使 J_n(x) = 0 的那些值——并不像正弦的 pi 整数倍那样等间距,但在远处会趋于近乎等间距,而这些零点恰好就是被允许的鼓的频率,正是圆鼓听起来不和谐、不像弦那样有确定音高的原因。同一个问题还有第二个解,即第二类贝塞尔函数 Y_n(x),它是那位在原点处爆向 -infinity 的搭档;参见第二类贝塞尔函数。只有当你区域的中心被排除时你才保留 Y_n——对一整面鼓面,你把它丢掉,好让答案在正中保持有限。
别把这层家族相似读过头。正交多项式是货真价实的多项式——有限的、精确的、你能写全的初等表达式。贝塞尔函数不是多项式,也没有初等封闭形式;J_n(x) 由一个无穷幂级数(或等价地由一个积分)定义,是一个非初等函数,正如前几篇里的伽马函数与误差函数。统一整架书架的,不是它们的代数类型,而是它们的角色:每一个都是为某种特定几何量身定制的一组正交积木。多项式用于平面与球面问题,贝塞尔函数用于圆形与柱形问题。
把它们全都系在一起的东西
退后一步,这架书架便显出一条统一的组织律,以后的台阶会把它说精确。这些家族中的每一个——勒让德、厄米、拉盖尔、切比雪夫、贝塞尔——都是某个特定二阶微分方程的解集,而这些方程并不是五个互不相干的偶然。它们全是同一个母模式(斯图姆-刘维尔形式)的实例,这个母模式自动保证:解会带着恰当的权正交地呈现出来。正交性不是我们手工强加的;它是从微分方程的结构里掉出来的定理,正如振动模态的相互垂直从一个对称矩阵里掉出来一样。
而且这些特定家族不断出现、而非某个别的随机清单,是有一个干净的理由的。当你在一个对称区域里取一个波动、热传导或量子问题并分离变量时,每个坐标里剩下的常微分方程,恰恰就是这些命名方程之一:球面角给你勒让德方程,柱面半径给你贝塞尔方程,振子给你厄米方程。问题的几何替你选定了家族。这就是为何同样这五个名字会在静电学、声学、热流与量子力学中反复出现——它们是那些真实物理可解的对称区域的天然词汇。
所以,带着这幅可用的图景而非那些公式离开本篇。你现在能一眼认出五条命名曲线:四个在各自区间上像弯折谐波般起伏的多项式家族,以及像渐弱正弦般向外荡漾的贝塞尔函数。你知道每一个都住在一处栖息地——一个区间和一个权——它使家族成员彼此正交,成为垂直坐标轴,你只需为每个系数算一个积分,就能沿着它们展开任何合理的函数。这恰好足够你跟上后面的台阶,在那里这些同样的对象将从它们的微分方程被推导出来,而那个母级的斯图姆-刘维尔模式终将被彻底揭开。