谱方法
概述
谱方法（Spectral Methods）是求解偏微分方程（PDE）的一类高精度数值方法，其核心思想是将未知函数展开为全局基函数（如切比雪夫多项式或傅里叶基函数）的有限级数，利用展开系数代替点值来离散化方程。对于光滑解，谱方法收敛速度是指数级（谱精度），远优于有限差分或有限元方法的代数收敛速度。代价是仅适用于规则计算域和光滑问题。
关键内容
核心思想
有限差分 / 有限元：函数用局部分段低次多项式表示，误差 $O(h^p)$（$h$ 为网格尺寸）。
谱方法：函数用 $N$ 个全局基函数展开 $u(x) \approx \sum_{k=0}^N \hat{u}_k\phi_k(x)$，误差随 $N$ 指数衰减：
$$|u - u_N|_\infty \leq C e^{-\sigma N}$$
（对解析函数，$\sigma > 0$ 取决于解的光滑半径）。
| 方法 | 误差阶 | 适用条件 |
|---|---|---|
| 有限差分 | $O(h^2)$ 或 $O(h^4)$ | 通用 |
| 有限元 | $O(h^p)$ | 通用，复杂几何 |
| 谱方法 | $e^{-\sigma N}$（指数） | 光滑解，简单域 |
两类基函数
傅里叶谱方法：适用于周期边界（流体周期通道流、信号处理），基函数 $e^{ikx}$；系数由 FFT 计算。
切比雪夫谱方法：适用于非周期边界，基函数 $T_k(x)$；节点取切比雪夫-高斯-洛巴托点（Chebyshev-Gauss-Lobatto nodes）$x_j = \cos(j\pi/N)$；系数由 DCT/FFT 计算。
Clenshaw-Curtis 求积
以切比雪夫极值点 $x_k = \cos(k\pi/n)$ 为节点，权重通过对切比雪夫展开积分得到：
$$\int_{-1}^1 f(x)dx \approx \sum_{k=0}^n w_k f(x_k)$$
优势：节点嵌套（加倍节点数时旧节点保留），权重可用 FFT 以 $O(n\log n)$ 计算，实践精度可与高斯求积公式媲美（Trefethen, 2008）。
配点法与 Galerkin 法
配点法（Pseudospectral）：要求方程在 $N+1$ 个配点处精确满足，等价于在物理空间操作，可直接利用 FFT/DCT。
Galerkin 谱方法：要求残差与每个基函数正交，在频率空间操作，便于误差分析。
代表性应用与软件
| 阶段 | 代表成果 |
|---|---|
| 理论奠基 | Canuto, Hussaini, Quarteroni, Zang《Spectral Methods in Fluid Dynamics》1988 |
| 教学实现 | Trefethen《Spectral Methods in MATLAB》2000 |
| Chebfun | 2004年牛津大学，用切比雪夫多项式表示和操作连续函数，"函数的 MATLAB" |
| 应用 | 流体力学直接数值模拟（DNS）、气候模拟、天体物理 |
局限性
- 几何限制：全局基函数只适合规则域（矩形、球壳）；复杂几何需谱元法（Spectral Element Method）+ 有限元思想结合
- 不适合非光滑问题：间断解出现吉布斯（Gibbs）现象，收敛速度退化为代数阶
- 稠密离散矩阵：全局基函数产生稠密（非稀疏）矩阵，时间步积分的稳定性条件苛刻
来源
- raw/books/数值分析/07_chebyshev_approximation.md