数値積分

出典: フリー百科事典『ウィキペディア(Wikipedia)』
ナビゲーションに移動 検索に移動
数学 > 数値解析 > 数値積分

数値積分(すうちせきぶん、: numerical integration)とは、狭義には与えられる関数の定積分の値を、解析的にではなく数値的に求める求積法のことであり、広義には与えられる導関数から原関数を求める手法、また微分方程式を数値的に解く手法 (常微分方程式の数値解法偏微分方程式の数値解法) を含む[1]数値解析の一分野である。

狭義の数値積分(関数の定積分の値を求める方法)は有限要素法などで応用されている[1]。 以下では、狭義の数値積分について述べる。

乱数を使わない積分[編集]

1変数の定積分の数値積分としては乱数を使わない方法として、

などがある。

ニュートン・コーツの公式の場合、誤差項は中点則と台形公式は同じ2階導関数、シンプソンの公式とシンプソンの3/8公式は同じ4階導関数なので、同じ誤差のグループ同士は滑らかな関数の場合は大きな差はなく、基本的にはシンプソンの公式の方が誤差が小さいが、場合によってはそうならない場合もある。

二重指数関数型数値積分公式、IMT積分[2]などの変数変換を用いた公式を適用すれば、被積分関数の端点に特異性がある場合でも、積分値を計算することが可能な場合もある。

多重積分[編集]

2変数以上の多重積分の場合は、外側から積分し、外の変数を定数として内側の積分を数値積分すれば良い。ただし一般に、変数が増えると、モンテカルロ法や準モンテカルロ法の方が計算効率が良くなる。 なお、1983年当時における多重積分の理論、アルゴリズムの状況は次元に応じて次のように分類されている[3][4]

  • Range O(2次元):満足できる状況.
  • Range I(3−7次元):プロダクトルールがその変形でなんとかなる.
  • Range II(7−15次元):Ranges I と IIIの 境界領域.
  • Range III(15次元以上):モンテカルロ法, 準モンテカルロ法が必要になる.

多次元の数値積分は金融工学などで必要とされているため、活発に研究されている[4][5]

精度保証付き数値積分[編集]

定積分の精度保証付き数値計算は、被積分関数の多項式近似と多項式の値を精度保証付き数値計算する技術の組み合わせで実現される。多項式の値を計算する際はホーナー法が使われ、多項式近似にはテイラー展開多項式補間が使われる[6]

特殊関数の零点を活用する近似公式[編集]

ガウス求積直交多項式の零点を活用する積分公式である[1]ベッセル関数の零点を活用する積分公式も開発されており、被積分関数がベッセル関数を含む場合に有効とされている[7][8]

超関数を使った積分[編集]

特異点を持つ積分の場合、佐藤超函数を経由して複素数値積分に持ち込むことで特異点を避け、誤差を抑えられる手法が提案されている[9][10][11][12]

部分積分を使った数値積分[編集]

振動積分に対しては部分積分を組み合わせた数値積分が提案されている[13][14]

乱数を使った積分[編集]

被積分関数の定義域が高次元であったりして近似公式がうまく計算できないような場合(次元の呪いが発生する場合)、モンテカルロ法が上手く行く場合もある。

  • 単純なモンテカルロ積分
  • 加重サンプリングを行う VEGAS 法
  • 層化抽出法を行うよう改良を加えた MISER 法

などがあり、このような数値積分法はモンテカルロ積分とよばれる。

一様乱数の代わりに超一様分布列英語版を使うと準モンテカルロ法英語版になり、より速く収束する場合がある。

確率論ランダム行列理論などを用いて上記の近似公式とモンテカルロ積分を融合させたアルゴリズムも提案されている[15]

積分範囲が無限区間の場合[編集]

積分範囲が無限区間の場合は、下記の方法で置換積分で変換して数値積分する方法がある。ただし関数によっては違う方法を利用した方が良い場合もある。無限大に近づくときに急激に0に収束することが解析的に分かっている場合は、積分範囲を有界で区切ってしまえば良い場合もある。

数値積分が困難な場合[編集]

数値積分のアルゴリズムは、滑らかである、連続である、特異点の数や場所が限られている、など、それぞれアルゴリズムが何らかの前提条件を関数にかけていて、それに沿った物が正常に積分できる。例えば、無数に不連続点がある関数は、基本的に数値積分が困難である。更に、広義積分可能であることも通常は前提条件に入っていて、コーシーの主値積分としては積分可能であっても、広義積分不可能な関数は、前提が崩れるためおかしな結果になる場合がある。

例えば

はコーシーの主値積分では 0 であるが、広義積分不可能であり[16]、これを数値積分すると(そもそも正しい結果が無いが)おかしな結果が返ってくることがある。定積分が広義積分可能であるかどうかは簡単には判断が付かない場合もあるため注意が必要である。

実装[編集]

近似計算[編集]

精度保証付き[編集]

出典[編集]

  1. ^ a b c 山本哲朗『数値解析入門』サイエンス社〈サイエンスライブラリ 現代数学への入門 14〉、2003年6月、増訂版。ISBN 4-7819-1038-6
  2. ^ 伊理正夫, 森口繁一, & 高澤嘉光. (1970). ある数値積分公式について (科学計算基本ライブラリーのアルゴリズムの研究会報告集).
  3. ^ Davis,P. J.,and Rabinowitz,P.,Methods of Nunlerical Integration,2nd Edition, Academic Press,San Diego,1984.
  4. ^ a b 手塚集、「数値多重積分に関する話題(<特集>数値計算)」 『応用数理』 1998年 8巻 4号 p.267-276, doi:10.11540/bjsiam.8.4_267, 日本応用数理学会
  5. ^ 量子コンピュータを用いた数値積分計算について, UNISYS TECHNOLOGY REVIEW 第90号, AUG. 2006.
  6. ^ 大石進一:「精度保証付き数値計算」、コロナ社、(1999年)
  7. ^ 緒方秀教, 杉原正顯、「Bessel 関数を含む振動積分に対する数値積分公式」 『数理解析研究所講究録』 1995年 915巻 p.171-181, 京都大学数理解析研究所
  8. ^ 緒方秀教, 杉原正顯、「Bessel関数の零点を標本点に持つ補間および数値積分公式」 『日本応用数理学会論文誌』 1996年 6巻 1号 p.39-66, doi:10.11540/jsiamt.6.1_39, 日本応用数理学会
  9. ^ 緒方秀教. (2018). 佐藤超函数論に基づく数値解析. 応用数理, 27(4), 8-15.
  10. ^ 緒方秀教. (2017). 佐藤超函数論に基づく数値積分 (現象解明に向けた数値解析学の新展開 (2)). 京都大学数理解析研究所講究録, (2037), 57-60.
  11. ^ 緒方秀教, & 平山弘. (2016). 数値積分に対する超函数法. 日本応用数理学会論文誌, 26(1), 33-43.
  12. ^ Ogata, H. (2019). A numerical method of computing oscillatory integral related to hyperfunction theory. arXiv preprint arXiv:1909.04911.
  13. ^ 平山弘. (1997). 部分積分法による半無限区間振動型積分の数値計算法. 日本応用数理学会論文誌, 7(2), 131-138.
  14. ^ 平山弘, 館野裕文, & 平野照比古. (2004). 部分積分法による数値積分法. 京都大学数理解析研究所講究録 1395: 190-195.
  15. ^ Nakatsukasa, Y. (2018). Approximate and integrate: Variance reduction in Monte Carlo integration via function approximation. arXiv preprint arXiv:1806.05492.
  16. ^ integrate from -1 to 1 1/x dx - Wolfram|Alpha
  17. ^ QUADPACK (nines)
  18. ^ quadpack
  19. ^ int_fcn_sing - IMSL C Math Library
  20. ^ Chapter 11. Quadrature and Differentiation - 1.69.0
  21. ^ S.M. Rump: INTLAB - INTerval LABoratory. In Tibor Csendes, editor, Developments in Reliable Computing, pages 77-104. Kluwer Academic Publishers, Dordrecht, 1999.
  22. ^ Johansson, F. (2017). Arb: efficient arbitrary-precision midpoint-radius interval arithmetic. IEEE Transactions on Computers, 66(8), 1281-1292.
  23. ^ Johansson, F. (2018, July). Numerical integration in arbitrary-precision ball arithmetic. In International Congress on Mathematical Software (pp. 255-263). Springer, Cham.
  24. ^ Johansson, F., & Mezzarobba, M. (2018). Fast and Rigorous Arbitrary-Precision Computation of Gauss--Legendre Quadrature Nodes and Weights. en:SIAM Journal on Scientific Computing, 40(6), C726-C747.

関連項目[編集]

関連文献[編集]

外部リンク[編集]