Fourier 的 Gibbs 现象 | Codeboy's Blog

形式陈述

增加 Fourier 模态，为什么间断旁的振荡越来越窄，最高超调却不消失？对具有跳跃的分段光滑周期函数，普通 Fourier 部分和在跳点处收敛到左右极限的平均值；在跳点附近按模态数缩放距离后，出现一个非零的极限超调形状。

以周期方波 $f(x)=\operatorname{sign}(x)$ （ $-\pi<x<\pi$ ）为例，它的Fourier 部分和可写为

S_m(x)=\frac4\pi\sum_{j=0}^{m}\frac{\sin((2j+1)x)}{2j+1}.

令 $x=z/(2m+2)$ 。对固定 $z$ ， $m\to\infty$ 时

S_m\left(\frac{z}{2m+2}\right)\longrightarrow \frac2\pi\operatorname{Si}(z),\qquad \operatorname{Si}(z)=\int_0^z\frac{\sin t}{t}\,dt.

第一处极大值在缩放位置 $z=\pi$ ，其高度趋于 $2\operatorname{Si}(\pi)/\pi\approx1.17898$ 。方波跳跃大小为二，超出右侧平台一的高度约为 $0.17898$ ，即跳跃的 $8.949\%$ 。

直觉

截断 Fourier 展开等于用 Dirichlet 核平滑原函数。这个核有振荡的正负旁瓣，跳跃被它涂开时就会出现超调。增加模态把核压窄，却不使缩放后的旁瓣形状消失。

上述极限还可以直接看成Riemann 和的极限：把 $\sin((2j+1)x)/(2j+1)$ 改写成缩放频率上的 $\sin(zs)/s$ ，奇数频率的网格逐渐填满 $[0,1]$ ，积分便给出 sine integral。

例子与边界

第一处峰的位置是 $x_m=\pi/(2m+2)$ ，因此 $m$ 从 $9$ 增加到 $99$ ，峰离跳点的距离由 $\pi/20$ 缩为 $\pi/200$ ，约缩小十倍。峰高却趋近同一个 $1.17898$ 。这就是“更窄”与“更低”两种误差概念的差别。

在固定 $0<x<\pi$ 处，部分和最终仍可收敛到一；但在整个包含跳点邻域上，连续三角多项式无法一致逼近这个间断函数。逐点收敛与一致收敛因此给出不同结论。对二值图像或守恒律，超调还可能生成负浓度或超过物理上界的值，即使整体积分误差已经很小。

把跳点的单个函数值改成零或一，不会去掉 Gibbs 现象。Fourier 系数对单点修改不敏感，振荡来自左右极限的跳跃。这个反例也帮助区分“跳点取值约定”与“邻域近似形状”。

推论与应用

可以平均相邻部分和、使用适当滤波，或把计算域在间断处分开。Fejér 平均使用非负核，减轻此类超调，代价是改变逼近算子与局部分辨率。对真正不连续的解，局部网格、限制器或其他非振荡方法往往更符合物理约束。

判断数值结果时应注明误差范数。一个振荡带宽趋零的近似，可能在积分意义上收敛，却仍保留不消失的最大误差。增加分辨率后只看整体缩略图，容易漏掉这个局部现象。

参考资料

NIST DLMF, §6.16(i) The Gibbs Phenomenon，方波部分和、sine integral 及峰值极限。
Trefethen, Approximation Theory and Approximation Practice, Chapter 9，Gibbs 现象与不同收敛意义。