[论文解读] A unified framework for the computation of polynomial quadrature weights and errors

Mário M. Graça, M. Esmeralda Sousa-Dias|arXiv (Cornell University)|Mar 21, 2012

Iterative Methods for Nonlinear Equations参考文献 15被引用 5

一句话总结

本文提出了一种统一框架，通过选择特定的多项式基底，利用待定系数法同时计算多项式求积规则中的求积权重与误差表达式。通过结合正交多项式与上三角线性系统，该方法实现了权重与误差项的递归、稳定计算——在牛顿-柯特斯、亚当斯-巴什福斯、亚当斯-莫尔顿及高斯求积规则中均成功验证，包括精确恢复著名的高斯-勒让德误差公式。

ABSTRACT

For the class of polynomial quadrature rules we show that conveniently chosen bases allow to compute both the weights and the theoretical error expression of a $n$-point rule via the undetermined coefficients method. As an illustration, the framework is applied to some classical rules such as Newton-Cotes, Adams-Bashforth, Adams-Moulton and Gaussian rules.

研究动机与目标

• 开发一个统一的计算框架，用于同时确定多项式求积规则中的求积权重与理论误差表达式。
• 克服传统方法中权重与误差计算分离的问题，后者通常涉及冗长且规则特异的推导。
• 表明选择合适的多项式基底可使待定系数法系统性地、递归地同时获得权重与误差表达式。
• 展示该框架在各类经典求积规则中的广泛适用性，包括高斯求积。
• 提供一种计算稳定且兼容符号计算与数值计算的算法，用于高精度计算权重与误差系数。

提出的方法

• 该方法采用一种特别构造的多项式基底，使得求积权重作为待定系数法求解上三角线性系统的结果被计算得出。
• 通过扩展该基底以包含更高次多项式，推导出精度阶次的可计算判据，该判据可直接导出理论误差表达式。
• 使用正交多项式（如勒让德多项式）以闭式表达误差项，尤其适用于高斯求积规则。
• 该框架依赖内积性质与正交性，将误差表达式导出为包含函数导数与类似佩亚诺核结构的积分形式。
• 该算法在 Mathematica 中实现，支持符号计算与高精度数值计算，并与已发表的表格数据进行验证。
• 通过监控权重之和（等于区间长度）与权重的正性，确保数值稳定性。

实验结果

研究问题

• RQ1 能否开发出一个单一、系统性的框架，用于同时计算多项式求积规则中的求积权重与误差表达式？
• RQ2 待定系数法如何扩展以同时获得理论误差表达式，而不仅限于权重？
• RQ3 何种多项式基底选择可实现权重与误差项的递归且稳定计算？
• RQ4 该框架在牛顿-柯特斯与高斯求积等多样化经典求积规则中，能在多大程度上实现统一应用？
• RQ5 该方法能否以高精度重现已知的闭式误差表达式，如高斯-勒让德规则的误差表达式？

主要发现

• 所提出的框架成功采用单一统一方法，计算了牛顿-柯特斯、亚当斯-巴什福斯、亚当斯-莫尔顿及高斯求积规则的权重与误差表达式。
• 对于高斯求积，该方法精确恢复了著名的误差表达式 $ E_{G_n}(f) = c_n f^{(2n)}(\bar{\xi}) $，其中 $ c_n = \frac{2^{2n+1}(n!)^4}{(2n+1)[(2n)!]^3} $，证实了其正确性。
• 该算法可计算高达 256 个节点的高斯-勒让德求积规则下 100 位精度的权重与误差系数，并与发表的表格数据（最高至 512 个节点）一致。
• 权重之和始终在高精度容差范围内等于 2（区间长度），确认了数值可靠性。
• 对于较小的 $ n $（最大至 5），在牛顿-柯特斯、亚当斯-巴什福斯与亚当斯-莫尔顿规则中，该方法获得精确的符号结果；对于较大的 $ n $，则获得高精度近似值。
• 该框架展示了计算上的可行性与稳定性，已在 Mathematica 中成功实现，并与多个来源的既有文献结果一致。