[论文解读] A corrected quadrature formula and applications

Nenad Ujević, A. J. Roberts|ArXiv.org|Mar 26, 2003

Aerospace Engineering and Control Systems被引用 16

一句话总结

本文提出了一种校正的三点闭合求积公式，通过引入端点导数信息改进了Simpson法则，实现了与网格间距相关的六阶精度。该方法同时修改了权重并增加了导数校正项，通过更紧致的误差界和优于标准复合Simpson法则的收敛速率，显著降低了误差——尤其在光滑被积函数下表现更优。

ABSTRACT

A straightforward 3-point quadrature formula of closed type is derived that improves on Simpson's rule. Just using the additional information of the integrand's derivative at the two endpoints we show the error is sixth order in grid spacing. Various error bounds for the quadrature formula are obtained to quantify more precisely the errors. Applications in numerical integration are given. With these error bounds, which are generally better than the usual Peano bounds, the composite formulas can be applied to integrands with lower order derivatives.

研究动机与目标

• 开发一种更高阶的求积规则，通过引入端点处的导数信息，改进Simpson法则。
• 推导一种校正求积公式，其误差界更优，尤其适用于低阶可微的被积函数。
• 提供具有更高精度和更紧误差估计的复合公式，优于标准复合Simpson法则。
• 通过数值例子和收敛性分析，展示其实际优越性。

提出的方法

• 通过修改Simpson法则的权重并增加涉及 $ f'(b) - f'(a) $ 的导数校正项，推导出校正的三点闭合求积规则。
• 利用有限差分法和算子微积分（涉及 $ \bar{P}_{ki}, \bar{Q}_{ki} $ 和 $ \beta(h\bar{\partial}) $）推导校正项。
• 通过Peano核表示和余项的积分形式建立误差界，其中 $ \bar{R}_k(f) = \sum \int \bar{T}_{ki}(x) f^{(k)}(x) dx $。
• 在 $ 2n $ 个子区间上以复合形式应用该校正规则，保持与标准Simpson法则相当的计算成本。
• 利用高阶导数的上确界范数和常数 $ D_k, E_k, \Gamma_k, \gamma_k $ 推导多种误差界。
• 通过数值例子验证该方法，包括 $ \int_0^1 e^{-x^2} dx $，并在双对数坐标系下绘制收敛速率图。

实验结果

研究问题

• RQ1 能否构造一种校正求积规则，通过引入端点导数值来改进Simpson法则？
• RQ2 此类校正规则的收敛阶是多少？与标准Simpson法则相比有何差异？
• RQ3 能否为校正规则推导出比标准Peano误差界更紧致的误差界？
• RQ4 与复合Simpson法则相比，校正规则的复合版本在实际应用中表现如何？
• RQ5 该校正规则是否在保持更高精度的同时维持计算效率？

主要发现

• 校正求积规则在网格间距 $ h $ 上达到六阶收敛，误差为 $ \bar{R}(f) \approx \frac{(b-a)^6}{302400} [f^{(5)}(b) - f^{(5)}(a)] $。
• 该规则对五次及以下的多项式精确成立，而标准Simpson法则仅对三次及以下多项式精确。
• 对于复合规则，误差界 $ |\bar{R}_4(f)| \leq 2\frac{\Gamma_4 - \gamma_4}{3645} h^4 (b-a) $ 比标准Simpson法则的误差界 $ \frac{\|f^{(4)}\|_\infty}{180} h^4 (b-a) $ 更为紧致。
• 数值结果表明，对于 $ \int_0^1 e^{-x^2} dx $，该校正规则仅用64个子区间即可达到六位有效数字精度，展现出快速收敛特性。
• 与复合Simpson法则相比，该方法仅增加一个额外项——$ -\frac{h^2}{15}[f'(b) - f'(a)] $——从而保持了计算效率。
• 校正规则在理论误差界和实际数值例子中均优于标准Simpson法则，尤其在光滑函数上表现更优。