一个思想,许多戏服
到现在为止,你已经在三四套戏服里见过同一个把戏了。普通的傅里叶级数把函数写成正弦余弦之和;勒让德级数把它写成勒让德多项式之和;傅里叶-贝塞尔级数则用贝塞尔函数。每一次,配方都感觉一模一样:把函数投影到一族构件上,每个构件一个积分。广义傅里叶级数不过是把这个配方只陈述一遍,而把构件留作一个待填的空白。无论你的几何递给你哪一组正交本征基,记账方式都相同。
为什么这件事会一再发生?因为所有这些基都源自同一处。前几篇指南已说明:斯图姆-刘维尔问题——带合适边界条件的本征值方程 (d/dx)(p y') + q y + lambda w y = 0——是自伴的,而自伴性逼着它的本征函数彼此正交。正弦、勒让德多项式、贝塞尔函数并不是三个互不相干的奇迹;它们是同一类方程在三个不同区间、带三种不同权函数下的三个解。正交性使展开变得廉价,而正交性是有保证的。所以「在一组正交本征基里展开任意函数」并非侥幸的特例——它正是整套理论的要旨所在。
通用的系数公式
下面就是统御一切的那一个公式。设 {y_n} 是你的斯图姆-刘维尔问题在 [a, b] 上的正交本征函数,带权函数 w(x)。你想把 f(x) 写成 sum over n of c_n y_n(x)。要钉住某个特定系数 c_m,把两边同乘 y_m(x) w(x) 再从 a 到 b 积分。在右边,正交性让除 n = m 那一项外的每一项都消失——恰如把向量与 x 轴作点积会抹掉 y、z 分量。存活下来的是 c_m 乘以 y_m^2 w dx 的积分。解出便得通用法则:c_m =(integral from a to b of f y_m w dx)除以(integral from a to b of y_m^2 w dx)。
看清两个部件。分子 integral of f y_m w dx 是投影——「f 中含有多少模式 m」。分母 integral of y_m^2 w dx 是 y_m 长度的平方,只起归一化作用;若你预先把每个 y_m 标度成单位带权范数,分母就成 1,c_m 直接是投影。决定性的角色是权 w,而它恰是初学者最爱漏掉的那一个。它不是可以无视的装饰:它正是自伴形式中坐在 lambda y 前面的那个系数,而本征函数是关于 w 正交的,从来不是关于裸的 dx。一旦漏掉 w,你的系数就会出错,级数也重建不出 f。
Generalized Fourier series: f(x) = sum_n c_n y_n(x)
Project onto y_m, integrate with weight w over [a,b]:
integral_a^b f y_m w dx = sum_n c_n * integral_a^b y_n y_m w dx
|
orthogonality kills every n != m
v
integral_a^b f y_m w dx = c_m * integral_a^b y_m^2 w dx
==> c_m = ( integral_a^b f y_m w dx ) / ( integral_a^b y_m^2 w dx )
Special cases (just change a,b,w):
sines on [0,L], w=1 -> ordinary Fourier coefficient
Legendre on [-1,1], w=1 -> c_n = (2n+1)/2 * integral f P_n dx
Bessel on [0,a], w=r -> Fourier-Bessel coefficient一幅做出来的图:球面上的一个台阶
想象一个球,上半球被涂成热的、下半球被涂成冷的——赤道以上温度为 +1,以下为 -1。在变量 x = cos(theta) 下,这就是 [-1, 1] 上的阶跃函数 f(x) = sign(x),而天然的基是勒让德多项式(权 w = 1)。由对称性,f 是奇函数,所以只有奇数次的 P_1, P_3, P_5, ... 存活;每个偶数系数在你动手算之前就恰为零,这正是「奇向量在偶坐标轴上没有分量」在函数空间里的回声。对第一个非零项套用那个通用公式,你得到 f 约等于 (3/2) P_1(x) 加上更小的修正,而其后每个系数都由各自独立的积分算出。
- 辨认几何的基与权:球面给出 [-1, 1] 上的勒让德多项式 P_n,权 w = 1。
- 先利用对称性:f(x) = sign(x) 是奇函数,所以所有偶数次系数无需任何积分就为零。
- 套用通用公式 c_n = ((2n+1)/2) 乘以 从 -1 到 1 对 f(x) P_n(x) dx 的积分,每个存活的 n 算一个积分。
- 把投影求和:f(x) = c_1 P_1 + c_3 P_3 + c_5 P_5 + ...,即冷热球的勒让德级数。
现在留意它与幂级数有何不同、又有何相同。泰勒级数也把 f 写成无穷和,但它的系数通过某一点处的各阶导数彼此耦合,而像 sign(x) 这样的单个跳跃会瞬间把它摧毁——你在 x = 0 处连导数都求不出。广义傅里叶级数没有这种脆弱:每个系数都是 f 对 y_n 的一个全局平均,由一个货真价实的定积分算出,它根本不在乎 f 有折角还是有跳跃。这种稳健——欣然吞下不连续的数据——正是这些展开(而非泰勒级数)成为处理带粗糙初值的边值问题之正确工具的原因。
级数在何种意义下等于 f?
下面这个问题,把仔细的读者和马虎的读者分开。当我们写 f(x) = sum of c_n y_n(x) 时,那个等号许诺了什么?诚实的回答是均方意义下收敛(也叫均方收敛或 L^2 收敛):随着你保留更多项,误差平方的带权积分 integral of (f - 部分和)^2 w dx 一路奔向零。这是一个关于累计总误差的陈述,而不是关于任何单独一点的陈述。级数可以在零散的几个点上有一丁点偏差却依然在均方意义下完美收敛,因为孤立的几处误差对积分毫无贡献。
级数在均方意义下收敛到 f、不留任何残余,恰恰就是本征基的完备性——对正则斯图姆-刘维尔问题而言,那条深刻的保证:本征函数触及函数空间中的每一个方向、一个都不漏。完备性是比正交性更强、更难的定理,但对正则问题它成立。要把分工看清楚:正交性使系数容易算出;完备性使所得级数真正等于 f。两者你都需要,而它们是不同的事实。
帕塞瓦尔:能量的账本
每个正交展开里都藏着一条守恒律,而它不过是被拉伸到无穷多维的勾股定理。在普通空间里,向量长度的平方等于其各分量平方之和。帕塞瓦尔恒等式对函数说的是同一句话:f 带权范数的平方等于其系数平方之和(每个再按其本征函数的范数平方加权)。用符号写,integral from a to b of f^2 w dx 等于 sum over n of c_n^2 乘以 integral of y_n^2 w dx——而若你已把 y_n 归一化为单位范数,它就干净地写作「f 的能量等于 c_n^2 之和」。函数的整个大小被逐模式精确核算,无所遗失,亦无所杜撰。
这不只是优雅——它是你的误差表。只保留前 N 项,帕塞瓦尔告诉你残余的均方误差恰是从 N+1 起的 c_n^2 尾和。于是你能预先决定为达到某容差需要多少模式:盯着尾和缩小,足够小时就停。帕塞瓦尔还是完备性的精确凭证。其单边版本贝塞尔不等式(c_n^2 之和至多等于范数)对任何正交集都成立;当 f 没有任何部分逃出基底时,它恰好升级为帕塞瓦尔等式。所以若帕塞瓦尔成为等式,你的基底就是完备的;若它有所亏欠,你便漏了一些方向。
它还附赠一份令人愉快的红利。因为帕塞瓦尔是一个精确的等式,你可以把它倒过来用,去求那些表面上与本征函数毫不相干的和。把 f(x) = x 在 [0, pi] 上展为正弦级数,把系数喂进帕塞瓦尔的傅里叶版本,那条著名的恒等式 sum of 1/n^2 = pi^2/6 便掉了出来——一个关于整数的事实,竟由一个关于振动弦的事实变出。这与本征函数情形下的同一恒等式,即本征函数的帕塞瓦尔恒等式,应用于所有斯图姆-刘维尔基中最简单的那一组,是同一回事。
为什么这是整卷书的引擎
退一步,看清下游的每个方法为何总要去抓正交的碎片。当你用分离变量法求解偏微分方程时,空间部分坍缩成一个斯图姆-刘维尔问题,它的本征基便白白送到你手上;于是本征函数展开法说:把你的初始数据展为广义傅里叶级数,又因为每个本征函数都是一个按自己简单时间表演化的纯模式,难解的偏微分方程便裂成一摞各自易解、再加回去的单模式问题。热量在棒中扩散、带电球周围的势、鼓面振响——全是同一招。
两条诚实的提醒,好让你怀着正确的期待。其一,这一切都立在完备性之上,而完备性对有限区间上的正则斯图姆-刘维尔问题有保证,在别处却可能失效;在无界域上,离散谱可能铺展成连续谱,求和变为积分,你就跨进了变换方法(下一阶梯)。其二,均方收敛确实弱于逐点收敛:级数在能量意义下是正确答案,跳跃处的吉布斯过冲永不消失,而对级数逐项求导后的收敛可能远比级数本身慢得多——所以对广义傅里叶级数求导务必小心。这些都不是缺陷;它是这套方法所给出的、精确而诚实的契约。