|
|
| 1行目: |
1行目: |
| − | [[数学]]の[[数値解析]]の分野における'''ガウス＝クロンロッド求積法'''（ガウス＝クロンロッドきゅうせきほう、{{Lang-en-short|Gauss–Kronrod quadrature formula}}）とは、（[[積分]]の近似値を計算するための）[[数値積分]]法の一種である。[[ガウス求積法]]の変形版であり、精度の低い近似での計算結果から得られる情報を再利用することで、より精度の高い近似を行うことが出来るように評価点を選ぶ求積法である。入れ子型求積則（nested quadrature rule）の一例で、函数の評価点の集合の中に高位と低位の二種類の求積則が存在する（後者は「埋め込み則」（embedded rule）と呼ばれる）。それら二つの近似の差は、積分の計算誤差を推定するために用いられる。
| + | {{テンプレート:20180815sk}} __NOINDEX__ |
| − | | |
| − | ガウス＝クロンロッド求積法は、1960年代にこの求積法を発見した{{仮リンク|アレクサンダー・クロンロッド|en|Alexander Kronrod}}と、[[カール・フリードリヒ・ガウス]]の名にちなむ。{{仮リンク|QUADPACK|en|QUADPACK}}ライブラリや、[[GNU Scientific Library]]、[[NAG数値計算ライブラリ|NAG Numerical Libraries]]、[[R言語]]で用いられている<ref>[http://stat.ethz.ch/R-manual/R-patched/library/stats/html/integrate.html http://stat.ethz.ch/R-manual/R-patched/library/stats/html/integrate.html]</ref>。
| |
| − | | |
| − | == 解説 ==
| |
| − | | |
| − | 数値積分の問題では、次の形式の定積分の近似値を求める。
| |
| − | | |
| − | :<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>
| |
| − | | |
| − | によって求めることが出来る。ここで ''w''<sub>''i''</sub> は重みであり、''x''<sub>''i''</sub> は函数 ''f''(''x'') の評価点である。
| |
| − | | |
| − | 区間 [''a'', ''b''] が細分されるとき、新しい区間のガウスの評価点は決して以前の評価点とは一致しない（奇数個の評価点の中央の点を除く）。したがって積分はそのような全ての点において評価される。ガウス＝クロンロッド求積法は、上述のガウス求積法にさらに <math>n+1</math> 個の評価点を加えることで、位数 <math>2n+1</math> となるように拡張された求積法である。そのような新たな点は、[[スティルチェス多項式]]の零点で与えられる。このような方法によって、函数の低位の推定値を再利用することにより、高位の推定を行うことが可能となる。ガウス求積法とガウス＝クロンロッド求積法の差は、しばしば近似誤差の推定に用いられる。
| |
| − | | |
| − | == 例 ==
| |
| − | | |
| − | ある有名な例では、7-点ガウス則と 15-点クロンロッド則が組み合わされる{{Harv|Kahaner|Moler|Nash|1989|loc=§5.5}}。ガウスの点はクロンロッドの点に組み込まれるため、求積および誤差推定に必要な函数の評価の総数は 15 となる。
| |
| − | | |
| − | :{| class="wikitable" style="background-color:white"
| |
| − | |+ [−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
| |
| − | |-
| |
| − | | style="text-align:right" | 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
| |
| − | |-
| |
| − | | style="text-align:right" | 0.00000 00000 00000 || ∗ || 0.20948 21410 84728
| |
| − | |}
| |
| − | | |
| − | 推奨される誤差推定は <math>(200 |G7 - K15|)^{1.5}</math> である。
| |
| − | | |
| − | {{harvtxt|Patterson|1968}} では、このタイプのさらなる拡張を見つける方法が示されている。
| |
| − | | |
| − | == 関連項目 ==
| |
| − | * {{仮リンク|クレンショウ＝カーティス求積法|en|Clenshaw–Curtis quadrature}} 同程度の精度を備える別の入れ子型求積則
| |
| − | | |
| − | == 注釈 ==
| |
| − | {{Reflist}}
| |
| − | | |
| − | == 参考文献 ==
| |
| − | * {{springer|title=Gauss–Kronrod quadrature formula|id=p/g120030}}
| |
| − | * {{citation | last1=Kahaner | first1=David | last2=Moler | first2=Cleve | author2-link=:en:Cleve Moler | last3=Nash | first3=Stephen | title=Numerical Methods and Software | year=1989 | publisher=[[:en:Prentice–Hall]] | isbn=978-0-13-627258-8 }}
| |
| − | * {{citation | last=Kronrod | first=Aleksandr Semenovish | author-link=:en:Alexander Kronrod | title=Nodes and weights of quadrature formulas. Sixteen-place tables | year=1965 | place=New York | publisher=Consultants Bureau | isbn= }} (Authorized translation from the Russian)
| |
| − | * {{citation | last1=Piessens | first1=Robert |last2=de Doncker-Kapenga | first2=Elise | last3=Überhuber | first3=C. W. | author3-link=:de:Christoph Überhuber | last4=Kahaner | first4=D. K. | title=QUADPACK, A subroutine package for automatic integration | year=1983 | publisher=[[Springer-Verlag]] | isbn=978-3-540-12553-2 }} (Reference guide for QUADPACK)
| |
| − | * {{citation | last=Patterson | first=T. N. L. | title=The Optimum Addition of Points to Quadrature Formulae | journal= Math. Comput. | volume=22 | pages=847–856 and C1–C11 | year=1968 | doi=10.2307/2004583 | issue=104 | publisher=American Mathematical Society | jstor=2004583}}. Erratum in ''Math. Comput.'' '''23''': 892.
| |
| − | | |
| − | == 外部リンク ==
| |
| − | * [http://www.netlib.org/slatec/src/qpdoc.f QUADPACK (part of SLATEC)], source code [http://www.netlib.org/slatec/src]. QUADPACK はガウス＝クロンロッド求積法に基づく数値積分のための [[Fortran]] でのアルゴリズムを集めたものである。（[[Netlib]] における）[[SLATEC]] は、数値計算のための広範なパブリック・ドメイン・ライブラリである。
| |
| − | * [http://www.alglib.net/integral/autogk61.php ALGLIB source code in C#, C++, Delphi & Visual Basic]
| |
| − | | |
| − | {{DEFAULTSORT:かうすくろんろつときゆうせきほう}}
| |
| − | [[Category:数値積分]]
| |
| − | [[Category:数学に関する記事]]
| |