ガウス=クロンロッド求積法
数学の数値解析の分野におけるガウス=クロンロッド求積法(ガウス=クロンロッドきゅうせきほう、英: Gauss–Kronrod quadrature formula)とは、(積分の近似値を計算するための)数値積分法の一種である。ガウス求積法の変形版であり、精度の低い近似での計算結果から得られる情報を再利用することで、より精度の高い近似を行うことが出来るように評価点を選ぶ求積法である。入れ子型求積則(nested quadrature rule)の一例で、函数の評価点の集合の中に高位と低位の二種類の求積則が存在する(後者は「埋め込み則」(embedded rule)と呼ばれる)。それら二つの近似の差は、積分の計算誤差を推定するために用いられる。
ガウス=クロンロッド求積法は、1960年代にこの求積法を発見したアレクサンダー・クロンロッドと、カール・フリードリヒ・ガウスの名にちなむ。QUADPACKライブラリや、GNU Scientific Library、NAG Numerical Libraries、R言語で用いられている[1]。
解説
数値積分の問題では、次の形式の定積分の近似値を求める。
- [math]\int_a^b f(x)\,dx.[/math]
このような積分の近似値は、例えば n-点ガウス求積法
- [math]\int_a^b f(x)\,dx \approx \sum_{i=1}^n w_i f(x_i) [/math]
によって求めることが出来る。ここで wi は重みであり、xi は函数 f(x) の評価点である。
区間 [a, b] が細分されるとき、新しい区間のガウスの評価点は決して以前の評価点とは一致しない(奇数個の評価点の中央の点を除く)。したがって積分はそのような全ての点において評価される。ガウス=クロンロッド求積法は、上述のガウス求積法にさらに [math]n+1[/math] 個の評価点を加えることで、位数 [math]2n+1[/math] となるように拡張された求積法である。そのような新たな点は、スティルチェス多項式の零点で与えられる。このような方法によって、函数の低位の推定値を再利用することにより、高位の推定を行うことが可能となる。ガウス求積法とガウス=クロンロッド求積法の差は、しばしば近似誤差の推定に用いられる。
例
ある有名な例では、7-点ガウス則と 15-点クロンロッド則が組み合わされる{{#invoke:Footnotes | harvard_citation }}。ガウスの点はクロンロッドの点に組み込まれるため、求積および誤差推定に必要な函数の評価の総数は 15 となる。
[−1,1] 上の (G7,K15) ガウス点 重み ±0.94910 79123 42759 ∗ 0.12948 49661 68870 ±0.74153 11855 99394 ∗ 0.27970 53914 89277 ±0.40584 51513 77397 ∗ 0.38183 00505 05119 0.00000 00000 00000 ∗ 0.41795 91836 73469 クロンロッド点 重み ±0.99145 53711 20813 0.02293 53220 10529 ±0.94910 79123 42759 ∗ 0.06309 20926 29979 ±0.86486 44233 59769 0.10479 00103 22250 ±0.74153 11855 99394 ∗ 0.14065 32597 15525 ±0.58608 72354 67691 0.16900 47266 39267 ±0.40584 51513 77397 ∗ 0.19035 05780 64785 ±0.20778 49550 07898 0.20443 29400 75298 0.00000 00000 00000 ∗ 0.20948 21410 84728
推奨される誤差推定は [math](200 |G7 - K15|)^{1.5}[/math] である。
Patterson (1968) では、このタイプのさらなる拡張を見つける方法が示されている。
関連項目
- クレンショウ=カーティス求積法 同程度の精度を備える別の入れ子型求積則
注釈
参考文献
- テンプレート:Springer
- Kahaner, David; Moler, Cleve; Nash, Stephen (1989), Numerical Methods and Software, en:Prentice–Hall, ISBN 978-0-13-627258-8
- Kronrod, Aleksandr Semenovish (1965), Nodes and weights of quadrature formulas. Sixteen-place tables, New York: Consultants Bureau (Authorized translation from the Russian)
- Piessens, Robert; de Doncker-Kapenga, Elise; Überhuber, C. W.; Kahaner, D. K. (1983), QUADPACK, A subroutine package for automatic integration, Springer-Verlag, ISBN 978-3-540-12553-2 (Reference guide for QUADPACK)
- Patterson, T. N. L. (1968), “The Optimum Addition of Points to Quadrature Formulae”, Math. Comput. (American Mathematical Society) 22 (104): 847–856 and C1–C11, doi:10.2307/2004583, JSTOR 2004583. Erratum in Math. Comput. 23: 892.
外部リンク
- QUADPACK (part of SLATEC), source code [1]. QUADPACK はガウス=クロンロッド求積法に基づく数値積分のための Fortran でのアルゴリズムを集めたものである。(Netlib における)SLATEC は、数値計算のための広範なパブリック・ドメイン・ライブラリである。
- ALGLIB source code in C#, C++, Delphi & Visual Basic