JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

帕塞瓦尔、能量与正交性

为什么三角函数系彼此垂直是让系数公式得以成立的那一个事实——以及帕塞瓦尔定理如何把这份垂直化作一笔毕达哥拉斯式的能量账,把信号的总能量逐个谐波地清算出来。

把函数当向量,把积分当点积

在本阶梯的第一篇里,我们用一个看起来像魔法的动作求出了傅里叶系数:乘上某一个波、在一个周期上积分,然后眼看着其余每一项都消失。那个动作并不是魔法——它是乔装打扮的几何,而本篇就把这层几何挑明,好让你看清为什么这一切都立得住。要抓住的那一个想法是:定义在一段区间上的函数,行为恰恰像一个向量,只不过它有无穷多个分量。要把这句口号变成能拿来计算的东西,你需要一种给两个函数取点积的办法。

对空间中的普通向量 u 和 v,点积是 u.v = u_1 v_1 + u_2 v_2 + u_3 v_3——你把对应的分量两两配对再求和。对一段区间上的两个函数 f 和 g,把「对应分量两两配对再求和」换成「把对齐的取值相乘、再在整段区间上加起来」,这正好就是一个定积分:定义内积为 <f, g> = f(x) g(x) 在一个周期上的积分。离散的求和变成连续的积分,是因为函数在每一点都有取值,而不只是在第 1、2、3 号槽位上有。有了这一个定义,几何的全部词汇——长度、夹角、垂直、投影——就整套地搬到了函数身上。

正交性:三角波是一组互相垂直的坐标轴

现在就为我们目录里的两个波算一算这个内积。那个核心事实——三角函数系的正交性——是说,任意两个不同的波,其内积恰好为零。在从 -L 到 L 的区间上:只要 m 不等于 n,sin(m pi x / L) sin(n pi x / L) dx 的积分就是 0;一个余弦配上另一个不同的余弦也是 0;而每一个正弦配上每一个余弦,全是 0。原因是具体的,并不神秘:两个不同频率的波之乘积本身就是一段振荡,它待在零线之上的时间恰与待在零线之下的时间一样多,于是在一整个周期里,正面积与负面积相消为乌有。

内积唯一为零的时候,是一个波遇上它自己的孪生兄弟。cos^2(n pi x / L) 在 -L 到 L 上的积分是 L,sin^2 也一样——因为一个被平方的波从不取负值,它的面积无从相消。用向量的话来说:每个波的平方长度都是 L(它不是单位向量,只是有一个固定长度),而任意两个不同的波彼此垂直。于是函数族 { 1, cos(pi x / L), sin(pi x / L), cos(2 pi x / L), sin(2 pi x / L), …… } 就是函数空间里一组互相垂直的坐标轴,正像标准基向量分别指向 x、y、z 那样——只不过这里有无穷多根。

无穷维里的毕达哥拉斯:帕塞瓦尔定理

正交性解锁的回报就在这里。对一个以垂直分量写出的普通向量,它的平方长度等于各分量平方之和:|v|^2 = v_1^2 + v_2^2 + v_3^2。这无非就是毕达哥拉斯定理。因为三角波是一组互相垂直的坐标轴,完全相同的定理对函数也成立——而且它有个名字,叫帕塞瓦尔定理。f 的「平方长度」是 <f, f> = [f(x)]^2 在一个周期上的积分,帕塞瓦尔说它等于每个谐波各自平方贡献之和,没有任何交叉项幸存。正是这份垂直杀死了那些交叉项,正如 |v|^2 里没有 x.y 这一项才让毕达哥拉斯定理这么干净。

PARSEVAL'S THEOREM  (trig Fourier series, period 2L)

  (1/L) * integral_{-L}^{L} [f(x)]^2 dx
          =  a_0^2 / 2  +  sum_{n>=1} ( a_n^2 + b_n^2 )

WHY it falls out of orthogonality:
  square the series  f = a_0/2 + sum [ a_n cos_n + b_n sin_n ]
  integrate term by term over one period.
    - cross terms  cos_m cos_n (m != n), sin_m sin_n, cos sin  -> integrate to 0
    - only the SELF terms survive:  integral cos_n^2 = L,  integral sin_n^2 = L
  divide by L  ->  each a_n, b_n contributes its own square, nothing else.

ENERGY reading:  left side = average power of the signal over a period;
                 right side = power summed harmonic by harmonic.
一框之内的帕塞瓦尔。把级数平方会刮起一阵交叉项的风暴,但正交性把它们一个不剩地抹去,只留下一个干净的系数平方和——把 f 的能量逐个频率清点出来。

为什么把左边叫做「能量」?在物理里,极多的量都按振幅的平方变化:振动弦的能量、电流耗散的功率、声音的强度。[f(x)]^2 在一个周期上的时间平均,恰恰就是那个平均功率。所以帕塞瓦尔说:一个信号的总功率等于它各个谐波所携带功率之和。从弯弯曲曲的时间图景过渡到柱状图般的频率图景——也就是离散频谱——其间没有任何损耗,也没有任何重复计数。正是这一句话,让工程师可以谈论「60 赫兹处坐着多少能量」,仿佛各个频率是彼此分隔、互不相干的格子。是正交性让这些格子彼此独立。

一笔做出来的账,以及白送给级数的一份礼物

我们把这笔账落实到第一篇里的那个方波上——那个奇波,在 (0, L) 上等于 +1、在 (-L, 0) 上等于 -1,它的系数是:n 为奇数时 b_n = 4/(n pi),其余为零。帕塞瓦尔的左边平凡得很:[f(x)]^2 处处就是 1(因为 f 不是 +1 就是 -1),于是它的平均就是 1。右边是 b_n^2 对奇数 n 求和,也就是 (4/pi)^2 乘以 [ 1 + 1/9 + 1/25 + 1/49 + …… ]。令两边相等就逼出 1 + 1/9 + 1/25 + …… = pi^2 / 8。从一个方波的能量平衡里,竟掉出一个著名无穷和的精确值——一个本来没道理这么容易求到的数。

这是帕塞瓦尔货真价实的一项超能力,值得停下来体会:它把一个你也许永远凑不出初等技巧去求和的难无穷级数,转化成一个你能算的常规积分。挑一个函数、求出它的傅里叶系数、写下能量平衡,你就把一个和式以闭形式求出来了。三角波给你 1/n^4 之和;别的选择能够到 1/n^6 乃至更高。这一整窝偶数次幂的和(也就是偶数 p 处那个收敛的 p 级数之值,历史上的巴塞尔问题及其表亲)都是这样掉出来的。你并没有发明什么新的求和技巧;你只是换个角度重新用了那条能量恒等式。

  1. 挑一个其平方容易积分的周期函数 f(常常是一个分段常数的波,或像 x、|x| 这样的简单多项式)。
  2. 照普通办法把它的傅里叶系数 a_n、b_n 算出来一次。
  3. 写下帕塞瓦尔:(1/L) 乘以 f^2 的积分 = a_0^2/2 + sum (a_n^2 + b_n^2)。左边是一个容易的单积分。
  4. 解出右边那个未知的级数——它现在等于一个你通过积分算出的、以闭形式给出的数。

完备性:诚实的细则

现在轮到课本常常一笔带过的那条警示。单凭正交性只能保证一个不等式——贝塞尔不等式——即各谐波能量之和至多等于总能量:a_0^2/2 + sum(a_n^2 + b_n^2) 小于或等于 (1/L) 乘以 f^2 的积分。从几何上看这是显然的:把一个向量投影到若干根垂直轴上,所捕获的长度绝不可能多过这个向量本身的长度。帕塞瓦尔把这个「小于或等于」升级成精确的等式,而这次升级是一个独立的、更强的论断。它说三角函数系不只是正交的,而且是完备的:没有任何一块剩余的 f 藏在你的坐标轴指不到的某个方向里。你在所有谐波上收集到的能量,一点都没漏。

能量和里还埋着一条实用的教训。因为各谐波能量加起来是一个有限的总数,所以项 a_n^2 + b_n^2 必须随 n 增大趋于零——它们不可能永远持续贡献下去。f 越光滑,它们衰减得越快(连续函数的系数比带跳跃的函数衰减得更快),于是绝大部分能量都住在头几个谐波里。这就是有损压缩背后的全部原理:留下响亮的低次谐波、丢掉微弱的高次谐波,而帕塞瓦尔预先精确地告诉你,你正在丢弃多少能量——丢弃多少保真度。同一套账目会在非周期信号那里以帕塞瓦尔-普朗歇尔定理的面貌重现,那里谐波之和变成了对一段连续频谱的积分。

为什么这是本阶梯的承重思想

退一步,看看你手里握着什么。整套傅里叶方法立在一个几何事实——正交性——之上,它让你独立地读出每一个系数;以及一个账目事实——帕塞瓦尔——它让你在从不离开频域的情况下追踪能量。两者合起来说:一个信号和它的频谱,是同一个对象的两份忠实描述;换一组基,什么也不丢。这与普通几何里挑选一组方便的垂直坐标轴完全是同一招,只不过「坐标轴」如今是一些波,而「空间」是无穷维的。

而且这个想法并不止步于正弦和余弦。正交性加完备性这一格局,是本卷往后一切内容的模板:当你求解热方程或一张振动的膜时,天然的搭建砌块未必是三角波,而是别的正交族——勒让德多项式、贝塞尔函数、某个边值问题的本征函数。它们每一个都带着自己的帕塞瓦尔式能量恒等式,而其抽象版本——本征函数的完备性——正是保证那些展开能够捕获整个函数的东西。在这里把三角情形当作纯几何来吃透,往后每一个「广义傅里叶级数」都不过是同一个故事,换了一套字母表来讲。