[论文解读] A Note on Regularized Shannon's Sampling Formulae

Liwen Qian, Guangyue Wei|ArXiv.org|Apr 30, 2000

Numerical methods in inverse problems被引用 14

一句话总结

本文为一种正则化Shannon采样公式提供了严格的误差估计，该公式通过结合高斯滤波器与Shannon采样核，提升了求解偏微分方程时的数值精度与鲁棒性。该方法减少了截断误差与正则化误差，在有限网格点数下实现了更高的精度。

ABSTRACT

Error estimation is given for a regularized Shannon's sampling formulae, which was found to be accurate and robust for numerically solving partial differential equations.

研究动机与目标

• 理论证明先前数值研究中观察到的正则化Shannon采样公式的高精度与鲁棒性。
• 推导使用正则化公式进行函数插值及其导数计算的精确误差界。
• 建立正则化显著降低实际计算中截断误差与近似误差的条件。
• 指导在有限计算域中选择参数(r, M)以实现期望的精度。
• 将正则化公式与标准Shannon采样进行比较，展示其更优的收敛性与更低的网格需求。

提出的方法

• 通过将sinc函数与高斯包络相乘引入正则化采样核：$ \text{exp}(-\frac{(t-n\Delta)^2}{2\sigma^2}) $。
• 利用Hermite多项式$ H_k $建模正则化核导数的行为，并推导误差界。
• 应用傅里叶分析表达正则化核的变换，并推导基于卷积的误差表示。
• 将总误差分解为三部分：正则化误差$ E_1 $，以及两个截断误差$ E_2, E_3 $，分别对应求和尾部。
• 采用$ L^2 $-范数估计与涉及Hermite多项式衰减和高斯尾部的不等式，对每个误差分量进行有界。
• 推导出复合误差界（定理），其依赖于$ \|f\|_{L^2} $，$ \|f^{(s)}\|_{L^2} $，网格间距$ \Delta $，以及参数$ r, M $，其中$ r $控制频带限制，$ M $控制截断范围。

实验结果

研究问题

• RQ1 如何从误差界的角度理论证明Shannon采样公式的正则化在函数插值与导数计算中的有效性？
• RQ2 函数与核参数需满足何种精确条件，才能确保正则化公式实现高精度与稳定性？
• RQ3 在截断误差与收敛率方面，正则化公式相较于标准Shannon采样有何差异？
• RQ4 为在实际PDE计算中实现期望精度$ \eta $，参数$ r $与$ M $的最优选择是什么？
• RQ5 $ L^2 $-范数误差估计如何随网格间距$ \Delta $，带宽$ B $，与正则化宽度$ \sigma $变化？

主要发现

• 正则化公式的误差界以$ \mathcal{O}(1/\sqrt{M\Delta}) $的速率衰减，显著小于标准Shannon采样在有限域中的$ \mathcal{O}(1/\sqrt{M\Delta}) $量级。
• 对于机器精度$ \eta \sim 15 $，所需正则化参数满足$ r > 2.8 $且$ r < 4.0 $，且$ M \sim 30 $，与先前的数值测试结果一致。
• 误差界主要由三项主导：导数范数$ \|f^{(s)}\|_{L^2} $，以及来自$ E_2 $与$ E_3 $的两个尾部贡献，两者均被$ \exp(-(M\Delta)^2/(2\sigma^2)) $指数抑制。
• 该方法可准确逼近几乎带限于$ B < \pi/\Delta $的函数，其中$ B $由$ r $确定，并在$ f $非严格带限时仍能实现高精度。
• Hermite多项式约束确保了正则化核导数的单调性，从而可通过微分不等式实现严格的误差估计。
• 计算效率得到提升：由于正则化核的更快衰减，实现相同精度所需的网格点数少于标准采样。