Skip to content

方法Method

谱配置法

Spectral collocation method

在全局插值空间中求导,并用内部配点和边界方程组成微分方程的离散系统。

形式陈述 ​

谱配点法怎样把微分方程变成一组有限维方程?先选 N+1 个不同节点 x0,…,xN,以节点值 uj 为未知量,用插值构造全局多项式

pN(x)=∑j=0Nujℓj(x),ℓj(xi)=δij.

定义 Dij=ℓj′(xi),则 (Du)i=pN′(xi)。二阶导数可用 D2u 表示。对边值问题 Lu=f,在指定内部节点要求 LpN(xi)=f(xi),剩余方程用于施加边界条件。

以 −u″=f、u(−1)=a,u(1)=b 为例,若节点含两个端点,就在 N−1 个内部节点使用 −D2u=f,再加两个端点方程。总计仍是 N+1 个未知数和方程。高维PDE可在张量网格或分域空间中作类似构造。

直觉

有限差分常用附近几个点描述局部斜率;全局谱配点把所有节点连成一张光滑曲面,再精确微分这张有限维曲面。微分矩阵通常稠密,因为一个节点值变化会影响整张插值多项式。

“配点”约束的是残差在选定节点处为零;Galerkin 方法约束的是残差对检验空间正交。两种方法都先选有限维空间,但让方程成立的方式不同。

例子与边界

三个点解一个边值问题 ​

解 −u″=2、u(−1)=u(1)=0。取节点 −1,0,1,唯一内部未知值记为 v。满足端点条件的二次插值多项式只能是 p2(x)=v(1−x2)。在内部节点 0 配点给出 −p2″(0)=2v=2,所以 v=1,得到 p2=1−x2。

这次解恰好落在试探空间中,因此方法精确恢复它。若右端换成 ex,同样三个点只能产生一个二次近似,需增加节点并检查误差。节点残差为零表示离散方程被解开,并不直接保证节点之间或真实解的误差小。

边界条件怎样进入矩阵 ​

令 D(2)=D2,把端点未知值固定后,内部系统为

−DII(2)uI=fI+DIB(2)uB.

这里 I、B 分别是内部与边界索引。先删掉微分矩阵端点行列,再对缩小后的 D 平方,通常得不到同一个二阶算子;边界值参与微分的作用会丢失。

推论与应用

算法依次选择适合边界的节点、构造导数算子、施加边界条件、求解、在独立检验点评估残差与误差。直接存储稠密矩阵需 O(N2) 空间,通用线性求解约 O(N3);快速变换和结构化求解可减少成本。

光滑度决定谱空间的逼近能力,离散稳定性决定小残差能否转成小误差。间断、边界层或复杂几何可采用分段、映射或局部谱元处理。Chebyshev 微分矩阵是区间问题的一种具体实现。

参考资料
  • Lloyd N. Trefethen, Spectral Methods in MATLAB, SIAM, 2000, Chapters 6–7,Chebyshev 微分与边值问题。
  • John P. Boyd, Chebyshev and Fourier Spectral Methods, 2nd ed., Dover, 2001,collocation 与 Galerkin 离散的比较。
关系图谱8 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组