← 返回 资料库 › 分析学 › 级数 分析学
傅里叶级数 几乎任何周期波,无论多么参差不齐,都可以通过叠加频率递增的纯正弦与余弦音调来构造。微妙之处不在于构造本身,而在于理解这个无穷和在何种意义下收敛回原来的波。
直观 用纯音构造一个波 音叉发出的乐音是纯正弦波。同样音高的小提琴音听起来却完全不同——因为它其实是基频加上更弱的谐波 之和,谐波频率是基频的 2 , 3 , 4 , … 2, 3, 4, \dots 2 , 3 , 4 , … 倍,各自的比例造就了小提琴特有的音色。约瑟夫·傅里叶(1807年)大胆的想法是:这并非乐器独有的性质——几乎任何 周期信号,包括像方波这样带尖角或跳变的信号,都可以写成频率为某个基频整数倍的正弦与余弦之和。
仅用 3 3 3 个正弦谐波逼近的方波。每增加一个谐波,起伏的曲线就更靠近方波的尖角一些。 大学 傅里叶系数:求出振幅 定义: 傅里叶级数
设 f f f 是周期为 2 π 2\pi 2 π 的函数(即对一切 x x x 有 f ( x + 2 π ) = f ( x ) f(x+2\pi)=f(x) f ( x + 2 π ) = f ( x ) ),在 [ − π , π ] [-\pi,\pi] [ − π , π ] 上可积。它的傅里叶级数 是 f ( x ) ∼ a 0 2 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty}\left(a_n\cos nx + b_n\sin nx\right) f ( x ) ∼ 2 a 0 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) ,而部分和 S N f ( x ) = a 0 2 + ∑ n = 1 N ( a n cos n x + b n sin n x ) S_N f(x) = \frac{a_0}{2} + \sum_{n=1}^{N}\left(a_n\cos nx + b_n\sin nx\right) S N f ( x ) = 2 a 0 + ∑ n = 1 N ( a n cos n x + b n sin n x ) 是仅用前 N N N 个谐波构造出的对 f f f 的最佳逼近。符号 ∼ \sim ∼ 是刻意使用的:这个级数是否真的收敛回 f ( x ) f(x) f ( x ) 是另一个问题,将在后文讨论。
a n = 1 π ∫ − π π f ( x ) cos n x d x ( n ≥ 0 ) , b n = 1 π ∫ − π π f ( x ) sin n x d x ( 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) a n = π 1 ∫ − π π f ( x ) cos n x d x ( n ≥ 0 ) , b n = π 1 ∫ − π π f ( x ) sin n x d x ( n ≥ 1 ) 这些公式并非凭空得来。函数 1 , cos x , sin x , cos 2 x , sin 2 x , … 1,\cos x,\sin x,\cos 2x,\sin 2x,\dots 1 , cos x , sin x , cos 2 x , sin 2 x , … 在 [ − π , π ] [-\pi,\pi] [ − π , π ] 上是正交 的:两个不同函数之积的积分为 0 0 0 ,而当 n ≥ 1 n\ge1 n ≥ 1 时 ∫ − π π cos 2 n x d x = ∫ − π π sin 2 n x d x = π \int_{-\pi}^{\pi}\cos^2 nx\,dx = \int_{-\pi}^{\pi}\sin^2 nx\,dx = \pi ∫ − π π cos 2 n x d x = ∫ − π π sin 2 n x d x = π 。把 f f f 的级数乘以 cos m x \cos mx cos m x 后逐项积分,除了与 m m m 匹配的那一项外其余全部消失,从而分离出 a m a_m a m ——这正是在向量空间中利用标准正交基求坐标的手法,只不过这里的"向量"是函数。
f ( x ) ∼ ∑ n = − ∞ ∞ f ^ ( n ) e i n x , f ^ ( n ) = 1 2 π ∫ − π π f ( x ) e − i n x d x f(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 f ( x ) ∼ n = − ∞ ∑ ∞ f ^ ( n ) e in x , f ^ ( n ) = 2 π 1 ∫ − π π f ( x ) e − in x d x 例题: 方波的傅里叶级数
求周期为 2 π 2\pi 2 π 的方波 f ( x ) = 1 f(x) = 1 f ( x ) = 1 (0 < x < π 0 < x < \pi 0 < x < π )、f ( x ) = − 1 f(x) = -1 f ( x ) = − 1 (− π < x < 0 -\pi < x < 0 − π < x < 0 )的傅里叶级数。
解答 因为 f f f 是奇函数,所以对一切 n n n 都有 a n = 0 a_n=0 a n = 0 。对于 b n b_n b n :b n = 1 π ∫ − π π f ( x ) sin n x d x = 2 π ∫ 0 π sin n x d x = 2 n π ( 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) b n = π 1 ∫ − π π f ( x ) sin n x d x = π 2 ∫ 0 π sin n x d x = nπ 2 ( 1 − cos nπ ) 。当 n n n 为偶数时此值为 0 0 0 ,当 n n n 为奇数时为 4 n π \frac{4}{n\pi} nπ 4 。于是 f ( x ) ∼ 4 π ∑ k = 0 ∞ sin ( ( 2 k + 1 ) x ) 2 k + 1 = 4 π ( sin x + sin 3 x 3 + sin 5 x 5 + ⋯ ) 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) f ( x ) ∼ π 4 ∑ k = 0 ∞ 2 k + 1 s i n ( ( 2 k + 1 ) x ) = π 4 ( sin x + 3 s i n 3 x + 5 s i n 5 x + ⋯ ) 。
同一个方波用 25 25 25 个谐波逼近。逼近效果精细得多,但在跳变附近,无论加多少谐波,都会超出目标值约跳变高度的 9 % 9\% 9% 。这种持续的过冲称为吉布斯现象 。 大学 帕塞瓦尔恒等式:能量守恒 若 f f f 的傅里叶系数为 a n , b n a_n, b_n a n , b n ,则 1 π ∫ − π π ∣ f ( x ) ∣ 2 d x = a 0 2 2 + ∑ n = 1 ∞ ( a n 2 + b n 2 ) \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 ∫ − π π ∣ f ( x ) ∣ 2 d x = 2 a 0 2 + ∑ n = 1 ∞ ( a n 2 + b n 2 ) 。
为什么成立? 这恰恰是周期函数无穷维空间中的勾股定理:由于 1 , cos n x , sin n x 1,\cos nx,\sin nx 1 , cos n x , sin n x 两两正交,f f f 的"长度平方" ∫ ∣ f ∣ 2 \int |f|^2 ∫ ∣ f ∣ 2 分解为它在各个方向上分量长度平方之和,积分之后交叉项全部消失。
证明 用 f f f 的级数展开 ∣ f ( x ) ∣ 2 |f(x)|^2 ∣ f ( x ) ∣ 2 ,再在 [ − π , π ] [-\pi,\pi] [ − π , π ] 上逐项积分。由正交性,当 n ≠ m n\ne m n = m 时 ∫ cos n x cos m x = ∫ sin n x sin m x = 0 \int \cos nx\cos mx = \int \sin nx\sin mx = 0 ∫ cos n x cos m x = ∫ sin n x sin m x = 0 ,且总有 ∫ cos n x sin m x = 0 \int \cos nx\sin mx = 0 ∫ cos n x sin m x = 0 ,因此只剩下"对角"项 a n 2 ∫ cos 2 n x = a n 2 π a_n^2\int\cos^2 nx = a_n^2\pi a n 2 ∫ cos 2 n x = a n 2 π 、b n 2 ∫ sin 2 n x = b n 2 π b_n^2\int\sin^2 nx = b_n^2\pi b n 2 ∫ sin 2 n x = b n 2 π 以及 ( a 0 / 2 ) 2 ∫ 1 = a 0 2 π / 2 (a_0/2)^2\int 1 = a_0^2\pi/2 ( a 0 /2 ) 2 ∫ 1 = a 0 2 π /2 留存;两边除以 π \pi π 即得所述公式。
例题: 由帕塞瓦尔恒等式导出巴塞尔问题
把帕塞瓦尔恒等式应用于 ( − π , π ] (-\pi, \pi] ( − π , π ] 上周期延拓的锯齿波 f ( x ) = x f(x) = x f ( x ) = x ,求 ∑ n = 1 ∞ 1 n 2 \sum_{n=1}^{\infty} \frac{1}{n^2} ∑ n = 1 ∞ n 2 1 ——即[巴塞尔问题](/problems/basel-problem)。
解答 因为 f f f 是奇函数,所以 a n = 0 a_n=0 a n = 0 ;简单计算得 b n = 2 π ∫ 0 π x sin n x d x = 2 ( − 1 ) n + 1 n b_n = \frac{2}{\pi}\int_0^\pi x\sin nx\,dx = \frac{2(-1)^{n+1}}{n} b n = π 2 ∫ 0 π x sin n x d x = n 2 ( − 1 ) n + 1 。帕塞瓦尔恒等式变为 1 π ∫ − π π x 2 d x = ∑ n = 1 ∞ b n 2 \frac{1}{\pi}\int_{-\pi}^{\pi} x^2\,dx = \sum_{n=1}^{\infty} b_n^2 π 1 ∫ − π π x 2 d x = ∑ n = 1 ∞ b n 2 。左边等于 1 π ⋅ 2 π 3 3 = 2 π 2 3 \frac{1}{\pi}\cdot\frac{2\pi^3}{3} = \frac{2\pi^2}{3} π 1 ⋅ 3 2 π 3 = 3 2 π 2 ,右边等于 ∑ n = 1 ∞ 4 n 2 = 4 ∑ n = 1 ∞ 1 n 2 \sum_{n=1}^{\infty} \frac{4}{n^2} = 4\sum_{n=1}^{\infty}\frac{1}{n^2} ∑ n = 1 ∞ n 2 4 = 4 ∑ n = 1 ∞ n 2 1 。令两边相等并求解,得 ∑ n = 1 ∞ 1 n 2 = π 2 6 \sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6} ∑ n = 1 ∞ n 2 1 = 6 π 2 ——这正是莱昂哈德·欧拉在1735年用 sin x \sin x sin x 的无穷乘积、以一种完全不同且更精细的论证首次公布的值;这种傅里叶分析式的推导是后来理论被严格化之后才出现的。
进阶 这个级数真的收敛回 f f f 吗? 在帕塞瓦尔恒等式的意义下——即按均方误差收敛,称为 L 2 L^2 L 2 收敛——对任意满足 ∫ ∣ f ∣ 2 < ∞ \int |f|^2 < \infty ∫ ∣ f ∣ 2 < ∞ 的 f f f ,答案总是肯定的。而在某个具体的 x x x 处逐点收敛则微妙得多。狄利克雷在1829年证明,只要 f f f 分段光滑,就有 S N f ( x ) → f ( x ) S_N f(x) \to f(x) S N f ( x ) → f ( x ) (在跳跃点处,级数收敛到左右极限的平均值)。人们曾普遍相信这可以推广到每一个连续函数,直到杜布瓦-雷蒙构造出一个傅里叶级数在某一点发散的连续函数(1873年)。柯尔莫哥洛夫在1923年更进一步,构造出一个傅里叶级数处处 发散的可积函数。
若 f ∈ L 2 ( − π , π ) f \in L^2(-\pi,\pi) f ∈ L 2 ( − π , π ) ,则对几乎处处的 x x x 都有 S N f ( x ) → f ( x ) S_N f(x) \to f(x) S N f ( x ) → f ( x ) (亨特于1968年将此推广到一切满足 p > 1 p > 1 p > 1 的 f ∈ L p f \in L^p f ∈ L p )。
为什么成立? 柯尔莫哥洛夫的反例存在于 L 1 L^1 L 1 中,这是使傅里叶系数仍有意义的最大空间;卡尔松定理(1966年)表明,只要 f f f 具有有限能量(L 2 L^2 L 2 ),这种麻烦就真的消失了。证明并不直接展示收敛性,而是通过一种精细的时间-频率分解,把函数拆成局部化的波包,从而控制最大部分和 sup N ∣ S N f ( x ) ∣ \sup_N |S_N f(x)| sup N ∣ S N f ( x ) ∣ 。
证明 完整的论证被认为是二十世纪分析学中最难的证明之一,远超本页范围;查尔斯·费弗曼在1973年给出了一个更清晰(但仍然技术性很强)的版本,利用了关于时间-频率"瓦片"的组合论证。主题[调和分析](/giai-tich-dieu-hoa)勾勒了其背后的图景。
历史注记
三角级数在傅里叶之前几十年就已出现,源于18世纪关于振动弦的一场争论。莱昂哈德·欧拉与让·勒朗·达朗贝尔找到了弦方程的行波解,但与丹尼尔·伯努利意见不合——伯努利在1753年断言,被拨动的弦任何 可能的形状都可以写成正弦模态 sin ( n x ) \sin(nx) sin ( n x ) 之和,欧拉认为这对于带尖角的形状而言是荒谬的。事实证明伯努利在本质上是对的。约瑟夫·傅里叶在1807年一篇关于热传导的论文(最初因缺乏严谨性被巴黎科学院拒绝)以及1822年的完整专著中重新提出这一想法,断言这种展开对几乎任意函数都成立;把这一断言变得精确,占据了此后一个半世纪分析学家的工作。
莱昂哈德·欧拉 约瑟夫·傅里叶
常见错误. f f f 的连续性不足以 保证 S N f ( x ) → f ( x ) S_N f(x) \to f(x) S N f ( x ) → f ( x ) 在每一点都成立——杜布瓦-雷蒙 1873 年的例子处处连续,但其傅里叶级数在某一点发散。而在跳跃间断点处,级数从不收敛到左右两个单侧值中的任何一个:它收敛到它们的平均值 f ( x − ) + f ( x + ) 2 \frac{f(x^-)+f(x^+)}{2} 2 f ( x − ) + f ( x + ) (对上面的方波而言,级数在 x = 0 x=0 x = 0 处恰好等于 0 0 0 ,而不是 ± 1 \pm1 ± 1 )。对于周期为 2 π 2\pi 2 π 的函数,傅里叶系数 a 0 a_0 a 0 定义为
1 π ∫ − π π f ( x ) d x \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\,dx π 1 ∫ − π π f ( x ) d x 1 2 π ∫ − π π f ( x ) d x \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\,dx 2 π 1 ∫ − π π f ( x ) d x 1 π ∫ − π π f ( x ) cos x d x \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos x\,dx π 1 ∫ − π π f ( x ) cos x d x 2 ∫ − π π f ( x ) d x 2\int_{-\pi}^{\pi} f(x)\,dx 2 ∫ − π π f ( x ) d x 若 f f f 是奇函数(即 f ( − x ) = − f ( x ) f(-x)=-f(x) f ( − x ) = − f ( x ) ),哪些傅里叶系数会自动为零?
余弦系数 a n a_n a n 正弦系数 b n b_n b n a n a_n a n 与 b n b_n b n 都为零两者都不为零 吉布斯现象是指,在跳跃间断点附近,部分和 S N f S_N f S N f
当 N N N 足够大时一致收敛到 f f f 无论 N N N 多大,都会超出跳变约 9 % 9\% 9% 的高度 在所有点都总是收敛到左右极限的平均值 对任意 N N N 都发散 对 ( − π , π ] (-\pi,\pi] ( − π , π ] 上 f ( x ) = x f(x)=x f ( x ) = x 的傅里叶级数使用帕塞瓦尔恒等式,∑ n = 1 ∞ 1 n 2 \sum_{n=1}^{\infty} \frac{1}{n^2} ∑ n = 1 ∞ n 2 1 等于
π 2 6 \frac{\pi^2}{6} 6 π 2 π 2 8 \frac{\pi^2}{8} 8 π 2 π 2 \frac{\pi}{2} 2 π π 2 4 \frac{\pi^2}{4} 4 π 2