[論文レビュー] A corrected quadrature formula and applications

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

Aerospace Engineering and Control Systems被引用数 16

ひとこと要約

本稿では、端点微分情報を取り入れることでシンプソン則を改善した修正3点閉型求積公式を提案する。グリッド間隔に関して6次精度を達成する。この手法は重みを修正するとともに、微分補正項を追加し、滑らかな被積分関数に対して特に顕著な誤差低減を実現する。標準の合成シンプソン則と比較して、誤差の上限が厳しく、収束速度も優れている。

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.

研究の動機と目的

• 端点における微分情報を取り入れることで、シンプソン則を上回る高次数の求積則を開発すること。
• 滑らかさが低い関数に対しても優れた誤差上限を得られる修正求積公式を導出すること。
• 標準の合成シンプソン則と比較して、より高い精度とタイトな誤差推定値を持つ合成公式を提供すること。
• 数値例と収束解析を通じて、実用的優位性を示すこと。

提案手法

• シンプソン則の重みを修正し、$ f'(b) - f'(a) $ を含む微分補正項を追加することで、修正3点閉型求積則を導出する。
• 有限差分および演算子計算（$ \bar{P}_{ki}, \bar{Q}_{ki} $、および $ \beta(h\bar{\partial}) $ を含む）を用いて補正項を導出する。
• ペアノ核表現および剰余項の積分形を用いて誤差上限を確立し、$ \bar{R}_k(f) = \sum \int \bar{T}_{ki}(x) f^{(k)}(x) dx $ とする。
• 2n個の部分区間における合成形として修正則を適用し、標準のシンプソン則と同等の計算コストを維持する。
• 高階微分のノルムと定数 $ D_k, E_k, \Gamma_k, \gamma_k $ を用いて複数の誤差上限を導出する。
• 数値例（$ \int_0^1 e^{-x^2} dx $ を含む）を用いて手法を検証し、対数-対数スケールで収束速度をプロットする。

実験結果

リサーチクエスチョン

• RQ1 端点微分値を組み込むことで、シンプソン則を改善した修正求積則を構築できるか？
• RQ2 このような修正則の収束次数は何か？ また、標準のシンプソン則と比較してどうなるか？
• RQ3 標準のペアノ上限よりも優れた、修正則に適したタイトな誤差上限を導出できるか？
• RQ4 修正則の合成版は、合成シンプソン則と比較して実際の性能で優れているか？
• RQ5 より高い精度を達成しつつ、計算効率を維持できるか？

主な発見

• 修正求積則はグリッド間隔 $ h $ に対して6次収束を達成し、誤差 $ \bar{R}(f) \approx \frac{(b-a)^6}{302400} [f^{(5)}(b) - f^{(5)}(a)] $ となる。
• この則は5次多項式に対して正確であるが、標準のシンプソン則は3次までしか正確でない。
• 合成則の場合、誤差上限 $ |\bar{R}_4(f)| \leq 2\frac{\Gamma_4 - \gamma_4}{3645} h^4 (b-a) $ は、標準のシンプソン則の上限 $ \frac{\|f^{(4)}\|_\infty}{180} h^4 (b-a) $ よりもタイトである。
• 数値結果では、$ \int_0^1 e^{-x^2} dx $ に対して64個の部分区間で6桁の精度が達成され、収束が著しく速いことが示された。
• 合成シンプソン則と比較して、追加の項 $ -\frac{h^2}{15}[f'(b) - f'(a)] $ のみを追加すればよく、計算効率が維持される。
• 滑らかな関数に対して、理論的誤差上限および実際の数値例の両面で、修正則は標準のシンプソン則を上回る性能を示した。