Helmholtz方程一维Dirichlet问题有限差分格式的傅里叶分析：尖锐估计与相对误差
数值分析
2026-05-01 v3 数值分析
摘要
我们考察一维不可定向Helmholtz方程的Dirichlet问题，即在区间(0,1)内求解u''+k^2u=f，边界条件为u(0)=g_0，u(1)=g_1，其中波数k为常数，满足k∈(0,∞)\πℕ，源项f属于H^p_0(0,1)，且p≥4。我们提出一种基于傅里叶分析的方法，推导有限差分方法在绝对误差和相对误差上的具有波数显式性质的尖锐估计。此类结果在有限元方法(FEM)中已为人所知。我们采用该方法分析经典的中心化有限差分格式。对于在齐次(或非齐次)Dirichlet边界条件下的离散解的傅里叶插值函数，我们在假设k>20且k(kh)^2/σ_k≤4/(π-2)（其中σ_k:=dist(k,πℕ)）的条件下，严格地证明了当∑_{p=0}^4k^{-p}||f^{(p)}||_{L^2}=O(1)（或|g_i|≍k^{-1}）时，绝对误差在L^2-范数下的最坊情况可达收敛阶为(kh)^2/σ_k^2，在H^1-半范数下为k(kh)^2/σ_k^2；相对误差在L^2-范数和H^1-半范数下的收敛阶均为k(kh)^2/σ_k。特别地，我们在同一阶数上严格建立了这些误差估计的下界，这是本文的主要新颖之处。我们还指出，傅里叶分析方法可作为评估具有源项情况下有限差分格式的便捷可视化工具，这超出了色散分析的范畴。理论与可视化分析的结果通过数值实验得到验证。
引用
@article{arxiv.2501.16696,
title = {Fourier Analysis of Finite Difference Schemes for the Helmholtz Equation in 1D with Dirichlet Conditions: Sharp Estimates and Relative Errors},
author = {Martin J. Gander and Hui Zhang and Haiyang Zhou},
journal= {arXiv preprint arXiv:2501.16696},
year = {2026}
}
备注
39 pages, 35 pictures grouped to 10 figures, significantly revised to clarify the dependence on $\sigma_k$