[论文解读] A Gaussian quadrature rule for oscillatory integrals on a bounded interval

Andreas Asheim, Alfredo Deaño+2 more|arXiv (Cornell University)|Dec 6, 2012

Mathematical functions and polynomials参考文献 16被引用 3

一句话总结

本文提出了一种针对区间 $[-1,1]$ 上带权 $e^{i\theta x}$ 的振荡积分的高斯求积规则，利用针对振荡核量身定制的正交多项式。该方法实现了最优渐近阶 $\mathcal{O}(\omega^{-2n-1})$，并在 $\omega \to 0$ 时平滑地退化为高斯-勒让德求积，数值证据证实了偶数阶规则的存在性，并表明在高频时收敛至超插值点。

ABSTRACT

We investigate a Gaussian quadrature rule and the corresponding orthogonal polynomials for the oscillatory weight function $e^{iωx}$ on the interval $[-1,1]$. We show that such a rule attains high asymptotic order, in the sense that the quadrature error quickly decreases as a function of the frequency $ω$. However, accuracy is maintained for all values of $ω$ and in particular the rule elegantly reduces to the classical Gauss-Legendre rule as $ω o 0$. The construction of such rules is briefly discussed, and though not all orthogonal polynomials exist, it is demonstrated numerically that rules with an even number of points are always well defined. We show that these rules are optimal both in terms of asymptotic order as well as in terms of polynomial order.

研究动机与目标

• 开发一种适用于有界区间上带权 $e^{i\omega x}$ 的振荡积分的高斯求积规则，确保在所有频率下保持高精度。
• 研究与振荡权 $e^{i\omega x}$ 相关的正交多项式的存在性与结构，特别是奇数次与偶数次多项式。
• 建立在高度振荡求积背景下，多项式精确性与渐近阶之间的联系。
• 证明该求积规则实现了最优渐近阶 $\mathcal{O}(\omega^{-2n-1})$，且当 $\omega \to \infty$ 时收敛至超插值点。
• 证明该规则在 $\omega \to 0$ 的极限下退化为经典高斯-勒让德规则，从而确保在整个频率范围内具有鲁棒性。

提出的方法

• 该方法基于关于内积 $\langle f, g \rangle = \int_{-1}^1 f(x)g(x)e^{i\omega x} dx$ 的正交多项式 $p_n^\omega$ 构建高斯求积规则。
• 求积节点被定义为 $2n$ 次正交多项式 $p_{2n}^\omega$ 的零点，从而确保对次数 $\leq 2n-1$ 的多项式精确积分。
• 分析了此类多项式的存在性，数值证据表明，对于所有 $\omega$，仅偶数次多项式存在，而奇数次多项式在某些临界 $\omega$ 值处因双线性形式非正定而失效。
• 渐近分析表明，当 $\omega \to \infty$ 时，求积节点收敛至数值最速下降法的超插值点，节点位置的误差为 $\mathcal{O}(\omega^{-2})$。
• 该方法利用克莱姆法则与范德蒙德矩阵分析，推导出求积权重的渐近行为，表明其衰减为 $\mathcal{O}(\omega^{-n-2})$。
• 建立了多项式精确性与渐近阶之间的正式联系，证明高次多项式精确性意味着高渐近阶 $\mathcal{O}(\omega^{-2n-1})$。

实验结果

研究问题

• RQ1 能否为有界区间上带权 $e^{i\omega x}$ 的振荡积分构造一个对所有 $\omega \geq 0$ 都有效的高斯求积规则？
• RQ2 当 $\omega \to \infty$ 时，求积节点与权重的渐近行为如何？它们与已知的超插值点有何关联？
• RQ3 为何某些 $\omega$ 值下奇数次正交多项式无法存在？在何种条件下偶数次多项式仍保持良好定义？
• RQ4 该求积规则如何在低频与高频下均保持最优精度，特别是在 $\omega \to 0$ 与 $\omega \to \infty$ 的极限下？
• RQ5 在此非标准正交多项式设置下，多项式精确性与渐近阶之间存在何种关系？

主要发现

• 所提出的高斯求积规则实现了渐近阶 $\mathcal{O}(\omega^{-2n-1})$，与高度振荡积分的最佳已知方法一致。
• 当 $\omega \to \infty$ 时，求积节点收敛至数值最速下降法的超插值点，节点位置的误差为 $\mathcal{O}(\omega^{-2})$。
• 当 $\omega \to 0$ 时，该规则退化为经典高斯-勒让德求积，确保了低频下的一致性与高精度。
• 数值实验表明，所有具有偶数个节点的高斯求积规则在所有 $\omega$ 下均存在，尽管奇数次多项式存在理论障碍。
• 对于 $2n$ 点规则，求积权重衰减为 $\mathcal{O}(\omega^{-n-2})$，与方法的高渐近阶一致。
• 该方法建立了多项式精确性（次数 $2n-1$）与渐近阶 $\mathcal{O}(\omega^{-2n-1})$ 之间的直接联系，证明了在两项度量下的最优性。