[论文解读] Error Estimates with Explicit Constants for Sinc Quadrature and Sinc Indefinite Integration over Infinite Intervals

Tomoaki Okayama|arXiv (Cornell University)|Feb 6, 2013

Numerical Methods and Algorithms参考文献 3被引用 4

一句话总结

本文针对半无穷和无穷区间上的Sinc积分与Sinc不定积分，利用双指数（DE）变换，提供了显式且可计算的误差界。文中提出一种新型DE变换 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$，其收敛速度更优且误差界中包含显式常数，从而实现了此前未被解决情形下的验证性自动积分。

ABSTRACT

The Sinc quadrature and the Sinc indefinite integration are approximation formulas for definite integration and indefinite integration, respectively, which can be applied on any interval by using an appropriate variable transformation. Their convergence rates have been analyzed for typical cases including finite, semi-infinite, and infinite intervals. In addition, for verified automatic integration, more explicit error bounds that are computable have been recently given on a finite interval. In this paper, such explicit error bounds are given in the remaining cases on semi-infinite and infinite intervals.

研究动机与目标

• 为此前未提供误差界的半无穷和无穷区间上的Sinc积分与Sinc不定积分，提供显式且可计算的误差界。
• 通过引入一种新型DE变换，提升第3类情形（半无穷区间上指数衰减被积函数）的收敛速度与误差估计精度。
• 建立新型DE变换 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$，其误差衰减速率优于现有变换。
• 通过在误差估计中提供显式常数，实现验证性自动积分，这对严格数值计算至关重要。

提出的方法

• 基于无穷与半无穷区间上的双指数变换，推导Sinc积分的显式误差界。
• 针对第3类情形，提出一种新型DE变换 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$，其保证反函数为初等函数，并提升收敛性。
• 建立形式为 $|\text{Error}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$ 的误差界，其中常数 $C$ 显式给出。
• 应用复分析技术，包括解析延拓与最大模原理，以在条状区域中界定变换后被积函数的上界。
• 使用函数类 $\mathbf{L}_{L,R,2\alpha,2\beta}^{\text{DE}}(\mathscr{D}_d)$ 来刻画变换后函数的衰减性与解析性特征。
• 通过在复平面上对对数与双曲函数的估计，特别是关键有理指数表达式的 $\sup$ 范数，证明误差界。

实验结果

研究问题

• RQ1 能否为半无穷与无穷区间上的Sinc积分推导出显式且可计算的误差界，而此前仅存在渐近界？
• RQ2 在第3类情形中，新型DE变换 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$ 的收敛速度与现有DE变换相比如何？
• RQ3 在新型变换的误差估计 $|\text{Error}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$ 中，显式常数 $C$ 的值是多少？
• RQ4 能否证明该新型变换的收敛速度优于 $\psi_{\text{DE3}}$ 与 $\psi_{\text{DE3$\dagger$}}$？
• RQ5 该新型变换在Sugihara不存在性定理的意义下是否为最优？即在第3类情形中，是否为该类变换中最佳可能？

主要发现

• 新型DE变换 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$ 的收敛速度优于 $\psi_{\text{DE3}}$ 与 $\psi_{\text{DE3$\dagger$}}$，其误差界为 $|\text{Error}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$。
• 误差界中的显式常数 $C$ 已推导并可计算，从而实现了半无穷与无穷区间上验证性自动积分。
• 该变换的误差界优于现有界 $\mathrm{e}^{-2\pi d n / \log(8dn/\mu)}$ 与 $\mathrm{e}^{-2\pi d n / \log(2\pi d n /\mu)}$。
• 函数 $F(\zeta) = f(\psi_{\text{DE3$\ddag$}}(\zeta)) \psi_{\text{DE3$\ddag$}}'(\zeta)$ 属于类 $\mathbf{L}_{L,R,2\alpha,2\beta}^{\text{DE}}(\mathscr{D}_d)$，其中 $L = 2(\tilde{c}_d)^{1-\alpha}K$ 且 $R = 2(\mathrm{e}^{\pi/12})^{1-\alpha}K$，确保其在复条状区域内的解析性与衰减性。
• 证明依赖于对复条状区域内表达式 $\left| \frac{1+\log(1+\mathrm{e}^{\pi \sinh \zeta})}{\log(1+\mathrm{e}^{\pi \sinh \zeta})} \cdot \frac{1}{1+\mathrm{e}^{-\pi \sinh \zeta}} \right| \leq \tilde{c}_d$ 的精确估计。
• 该新型变换在Sugihara不存在性定理的意义下被证明为最优，表明其为第3类情形中该类变换的最优可能。