MathLabs

分析学

傅里叶级数

几乎任何周期波,无论多么参差不齐,都可以通过叠加频率递增的纯正弦与余弦音调来构造。微妙之处不在于构造本身,而在于理解这个无穷和在何种意义下收敛回原来的波。

直观用纯音构造一个波

音叉发出的乐音是纯正弦波。同样音高的小提琴音听起来却完全不同——因为它其实是基频加上更弱的谐波之和,谐波频率是基频的 2,3,4,…2, 3, 4, \dots 倍,各自的比例造就了小提琴特有的音色。约瑟夫·傅里叶(1807年)大胆的想法是:这并非乐器独有的性质——几乎任何周期信号,包括像方波这样带尖角或跳变的信号,都可以写成频率为某个基频整数倍的正弦与余弦之和。

方波及其前三个奇次正弦谐波之和的图像,已粗略勾勒出方波的形状。
仅用 33 个正弦谐波逼近的方波。每增加一个谐波,起伏的曲线就更靠近方波的尖角一些。

大学傅里叶系数:求出振幅

定义: 傅里叶级数

设 ff 是周期为 2π2\pi 的函数(即对一切 xx 有 f(x+2π)=f(x)f(x+2\pi)=f(x)),在 [−π,π][-\pi,\pi] 上可积。它的傅里叶级数是 f(x)∼a02+∑n=1∞(ancos⁡nx+bnsin⁡nx)f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty}\left(a_n\cos nx + b_n\sin nx\right),而部分和 SNf(x)=a02+∑n=1N(ancos⁡nx+bnsin⁡nx)S_N f(x) = \frac{a_0}{2} + \sum_{n=1}^{N}\left(a_n\cos nx + b_n\sin nx\right) 是仅用前 NN 个谐波构造出的对 ff 的最佳逼近。符号 ∼\sim 是刻意使用的:这个级数是否真的收敛回 f(x)f(x) 是另一个问题,将在后文讨论。

an=1π∫−ππf(x)cos⁡nx dx (n≥0),bn=1π∫−ππf(x)sin⁡nx dx (n≥1)a_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos nx\,dx\ (n\ge 0), \qquad b_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\sin nx\,dx\ (n\ge 1)

这些公式并非凭空得来。函数 1,cos⁡x,sin⁡x,cos⁡2x,sin⁡2x,…1,\cos x,\sin x,\cos 2x,\sin 2x,\dots 在 [−π,π][-\pi,\pi] 上是正交的:两个不同函数之积的积分为 00,而当 n≥1n\ge1 时 ∫−ππcos⁡2nx dx=∫−ππsin⁡2nx dx=π\int_{-\pi}^{\pi}\cos^2 nx\,dx = \int_{-\pi}^{\pi}\sin^2 nx\,dx = \pi。把 ff 的级数乘以 cos⁡mx\cos mx 后逐项积分,除了与 mm 匹配的那一项外其余全部消失,从而分离出 ama_m——这正是在向量空间中利用标准正交基求坐标的手法,只不过这里的"向量"是函数。

f(x)∼∑n=−∞∞f^(n) einx,f^(n)=12π∫−ππf(x) e−inx dxf(x) \sim \sum_{n=-\infty}^{\infty} \hat f(n)\, e^{inx}, \qquad \hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\, e^{-inx}\, dx

例题: 方波的傅里叶级数

求周期为 2π2\pi 的方波 f(x)=1f(x) = 1(0<x<π0 < x < \pi)、f(x)=−1f(x) = -1(−π<x<0-\pi < x < 0)的傅里叶级数。

解答

因为 ff 是奇函数,所以对一切 nn 都有 an=0a_n=0。对于 bnb_n:bn=1π∫−ππf(x)sin⁡nx dx=2π∫0πsin⁡nx dx=2nπ(1−cos⁡nπ)b_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\sin nx\,dx = \frac{2}{\pi}\int_0^\pi \sin nx\,dx = \frac{2}{n\pi}\left(1-\cos n\pi\right)。当 nn 为偶数时此值为 00,当 nn 为奇数时为 4nπ\frac{4}{n\pi}。于是 f(x)∼4π∑k=0∞sin⁡((2k+1)x)2k+1=4π(sin⁡x+sin⁡3x3+sin⁡5x5+⋯ )f(x) \sim \frac{4}{\pi}\sum_{k=0}^{\infty} \frac{\sin\left((2k+1)x\right)}{2k+1} = \frac{4}{\pi}\left(\sin x + \frac{\sin 3x}{3} + \frac{\sin 5x}{5} + \cdots\right)。

方波及其25次谐波傅里叶逼近的图像,整体呈现细密波纹,每个跳跃间断点之后都有一个持续的过冲尖峰。
同一个方波用 2525 个谐波逼近。逼近效果精细得多,但在跳变附近,无论加多少谐波,都会超出目标值约跳变高度的 9%9\%。这种持续的过冲称为吉布斯现象。

大学帕塞瓦尔恒等式:能量守恒

若 ff 的傅里叶系数为 an,bna_n, b_n,则 1π∫−ππ∣f(x)∣2 dx=a022+∑n=1∞(an2+bn2)\frac{1}{\pi}\int_{-\pi}^{\pi} |f(x)|^2\,dx = \frac{a_0^2}{2} + \sum_{n=1}^{\infty}\left(a_n^2 + b_n^2\right)。

为什么成立?

这恰恰是周期函数无穷维空间中的勾股定理:由于 1,cos⁡nx,sin⁡nx1,\cos nx,\sin nx 两两正交,ff 的"长度平方" ∫∣f∣2\int |f|^2 分解为它在各个方向上分量长度平方之和,积分之后交叉项全部消失。

证明

用 ff 的级数展开 ∣f(x)∣2|f(x)|^2,再在 [−π,π][-\pi,\pi] 上逐项积分。由正交性,当 n≠mn\ne m 时 ∫cos⁡nxcos⁡mx=∫sin⁡nxsin⁡mx=0\int \cos nx\cos mx = \int \sin nx\sin mx = 0,且总有 ∫cos⁡nxsin⁡mx=0\int \cos nx\sin mx = 0,因此只剩下"对角"项 an2∫cos⁡2nx=an2πa_n^2\int\cos^2 nx = a_n^2\pi、bn2∫sin⁡2nx=bn2πb_n^2\int\sin^2 nx = b_n^2\pi 以及 (a0/2)2∫1=a02π/2(a_0/2)^2\int 1 = a_0^2\pi/2 留存;两边除以 π\pi 即得所述公式。

例题: 由帕塞瓦尔恒等式导出巴塞尔问题

把帕塞瓦尔恒等式应用于 (−π,π](-\pi, \pi] 上周期延拓的锯齿波 f(x)=xf(x) = x,求 ∑n=1∞1n2\sum_{n=1}^{\infty} \frac{1}{n^2}——即[巴塞尔问题](/problems/basel-problem)。

解答

因为 ff 是奇函数,所以 an=0a_n=0;简单计算得 bn=2π∫0πxsin⁡nx dx=2(−1)n+1nb_n = \frac{2}{\pi}\int_0^\pi x\sin nx\,dx = \frac{2(-1)^{n+1}}{n}。帕塞瓦尔恒等式变为 1π∫−ππx2 dx=∑n=1∞bn2\frac{1}{\pi}\int_{-\pi}^{\pi} x^2\,dx = \sum_{n=1}^{\infty} b_n^2。左边等于 1π⋅2π33=2π23\frac{1}{\pi}\cdot\frac{2\pi^3}{3} = \frac{2\pi^2}{3},右边等于 ∑n=1∞4n2=4∑n=1∞1n2\sum_{n=1}^{\infty} \frac{4}{n^2} = 4\sum_{n=1}^{\infty}\frac{1}{n^2}。令两边相等并求解,得 ∑n=1∞1n2=π26\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6}——这正是莱昂哈德·欧拉在1735年用 sin⁡x\sin x 的无穷乘积、以一种完全不同且更精细的论证首次公布的值;这种傅里叶分析式的推导是后来理论被严格化之后才出现的。

进阶这个级数真的收敛回 ff 吗?

在帕塞瓦尔恒等式的意义下——即按均方误差收敛,称为 L2L^2 收敛——对任意满足 ∫∣f∣2<∞\int |f|^2 < \infty 的 ff,答案总是肯定的。而在某个具体的 xx 处逐点收敛则微妙得多。狄利克雷在1829年证明,只要 ff 分段光滑,就有 SNf(x)→f(x)S_N f(x) \to f(x)(在跳跃点处,级数收敛到左右极限的平均值)。人们曾普遍相信这可以推广到每一个连续函数,直到杜布瓦-雷蒙构造出一个傅里叶级数在某一点发散的连续函数(1873年)。柯尔莫哥洛夫在1923年更进一步,构造出一个傅里叶级数处处发散的可积函数。

若 f∈L2(−π,π)f \in L^2(-\pi,\pi),则对几乎处处的 xx 都有 SNf(x)→f(x)S_N f(x) \to f(x)(亨特于1968年将此推广到一切满足 p>1p > 1 的 f∈Lpf \in L^p)。

为什么成立?

柯尔莫哥洛夫的反例存在于 L1L^1 中,这是使傅里叶系数仍有意义的最大空间;卡尔松定理(1966年)表明,只要 ff 具有有限能量(L2L^2),这种麻烦就真的消失了。证明并不直接展示收敛性,而是通过一种精细的时间-频率分解,把函数拆成局部化的波包,从而控制最大部分和 sup⁡N∣SNf(x)∣\sup_N |S_N f(x)|。

证明

完整的论证被认为是二十世纪分析学中最难的证明之一,远超本页范围;查尔斯·费弗曼在1973年给出了一个更清晰(但仍然技术性很强)的版本,利用了关于时间-频率"瓦片"的组合论证。主题[调和分析](/giai-tich-dieu-hoa)勾勒了其背后的图景。

对于周期为 2π2\pi 的函数,傅里叶系数 a0a_0 定义为

若 ff 是奇函数(即 f(−x)=−f(x)f(-x)=-f(x)),哪些傅里叶系数会自动为零?

吉布斯现象是指,在跳跃间断点附近,部分和 SNfS_N f

对 (−π,π](-\pi,\pi] 上 f(x)=xf(x)=x 的傅里叶级数使用帕塞瓦尔恒等式,∑n=1∞1n2\sum_{n=1}^{\infty} \frac{1}{n^2} 等于

参考文献

  1. Joseph Fourier (1822). Théorie analytique de la chaleur · DOI:10.1017/cbo9780511693229
  2. Elias M. Stein, Rami Shakarchi (2003). Fourier Analysis: An Introduction
  3. T. W. Körner (1988). Fourier Analysis
  4. Lennart Carleson (1966). On convergence and growth of partial sums of Fourier series