「数値積分」の版間の差分

削除された内容追加された内容

インライン

2019年9月1日 (日) 12:54時点における版

数値積分（すうちせきぶん、英: 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]。

乱数を使った積分

→詳細は「モンテカルロ法 § 数値積分」を参照

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

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

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

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

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

積分範囲が無限区間の場合

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

\int _{-\infty }^{\infty }f(x)\,dx=\int _{-1}^{1}f\left({\frac {t}{1-t^{2}}}\right){\frac {1+t^{2}}{(1-t^{2})^{2}}}\,dt

\int _{-\infty }^{\infty }f(x)\,dx=\int _{0}^{1}\left(f\left({\frac {1-t}{t}}\right)+f\left(-{\frac {1-t}{t}}\right)\right){\frac {1}{t^{2}}}\,dt

\int _{a}^{\infty }f(x)\,dx=\int _{0}^{1}f\left(a+{\frac {1-t}{t}}\right){\frac {1}{t^{2}}}\,dt

\int _{-\infty }^{b}f(x)\,dx=\int _{0}^{1}f\left(b-{\frac {1-t}{t}}\right){\frac {1}{t^{2}}}\,dt

数値積分が困難な場合

数値積分のアルゴリズムは、滑らかである、連続である、特異点の数や場所が限られている、など、それぞれアルゴリズムが何らかの前提条件を関数にかけていて、それに沿った物が正常に積分できる。例えば、無数に不連続点がある関数は、基本的に数値積分が困難である。更に、広義積分可能であることも通常は前提条件に入っていて、コーシーの主値積分としては積分可能であっても、広義積分不可能な関数は、前提が崩れるためおかしな結果になる場合がある。例えば $\int _{-1}^{1}{\frac {1}{x}}\,dx$ はコーシーの主値積分では 0 であるが、広義積分不可能であり^[11]、これを数値積分すると（そもそも正しい結果が無いが）おかしな結果が返ってくることがある。定積分が広義積分可能であるかどうかは簡単には判断が付かない場合もあるため注意が必要である。

実装

QUADPACK - 乱数を使わない1変数の定積分が FORTRAN 77 で実装されている^[12]^[13]。
- GNU Scientific Library ではこれをC言語に移植している。モンテカルロ積分や超一様分布列も実装されている。
- SciPy、R言語、GNU Octave、NAG数値計算ライブラリ、IMSL^[14] などが QUADPACK を使用している。
Boost Math Toolkit^[15]
Apache Commons Math では Java で実装されている。
Cuba - 主にモンテカルロ積分が実装されている
Mathematica・Maple・MATLAB などの数式処理システムおよび数値解析ソフトウェア

出典

^ ^a ^b ^c 山本哲朗『数値解析入門』（増訂版）サイエンス社〈サイエンスライブラリ現代数学への入門 14〉、2003年6月。ISBN 4-7819-1038-6。
^ 伊理正夫, 森口繁一, & 高澤嘉光. (1970). ある数値積分公式について (科学計算基本ライブラリーのアルゴリズムの研究会報告集).
^ Davis，P. J.，and Rabinowitz，P.，Methods of Nunlerical Integration，2nd Edition， Academic Press，San Diego，1984．
^ ^a ^b 手塚集. (1998). 数値多重積分に関する話題 (<特集> 数値計算). 応用数理, 8(4), 267-276. 日本応用数理学会より発行
^ 量子コンピュータを用いた数値積分計算について, UNISYS TECHNOLOGY REVIEW 第90号, AUG. 2006.
^ 大石進一：「精度保証付き数値計算」、コロナ社、（1999年）
^ 緒方秀教, 杉原正顯, Bessel 関数を含む振動積分に対する数値積分公式, 数理解析研究所講究録915 巻 1995 年 171-181.
^ 緒方秀教, 杉原正顯, Bessel関数の零点を標本点に持つ補間および数値積分公式, 日本応用数理学会論文誌, 6 巻 (1996) 1 号.
^ 緒方秀教. (2017). 佐藤超函数論に基づく数値積分 (現象解明に向けた数値解析学の新展開 (2)). 数理解析研究所講究録, (2037), 57-60.
^ Nakatsukasa, Y. (2018). Approximate and integrate: Variance reduction in Monte Carlo integration via function approximation. arXiv preprint arXiv:1806.05492.
^ integrate from -1 to 1 1/x dx - Wolfram|Alpha
^ QUADPACK (nines)
^ quadpack
^ int_fcn_sing - IMSL C Math Library
^ Chapter 11. Quadrature and Differentiation - 1.69.0

外部リンク

Lobatto Quadrature Wolfram Mathworld の記事
Lobatto quadrature formula Encyclopedia of Mathematics の記事

この項目は、解析学に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています（プロジェクト:数学／Portal:数学）。

[Yamamoto1-1] 山本哲朗『数値解析入門』（増訂版）サイエンス社〈サイエンスライブラリ現代数学への入門 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．

[JSIAM-4] 手塚集. (1998). 数値多重積分に関する話題 (<特集> 数値計算). 応用数理, 8(4), 267-276. 日本応用数理学会より発行

[5] 量子コンピュータを用いた数値積分計算について, UNISYS TECHNOLOGY REVIEW 第90号, AUG. 2006.

[6] 大石進一：「精度保証付き数値計算」、コロナ社、（1999年）

[7] 緒方秀教, 杉原正顯, Bessel 関数を含む振動積分に対する数値積分公式, 数理解析研究所講究録915 巻 1995 年 171-181.

[8] 緒方秀教, 杉原正顯, Bessel関数の零点を標本点に持つ補間および数値積分公式, 日本応用数理学会論文誌, 6 巻 (1996) 1 号.

[9] 緒方秀教. (2017). 佐藤超函数論に基づく数値積分 (現象解明に向けた数値解析学の新展開 (2)). 数理解析研究所講究録, (2037), 57-60.

[10] Nakatsukasa, Y. (2018). Approximate and integrate: Variance reduction in Monte Carlo integration via function approximation. arXiv preprint arXiv:1806.05492.

[11] tegrate from -1 to 1 1/x dx - Wolfram|Alpha

[12] QUADPACK (nines)

[13] quadpack

[14] t_fcn_sing - IMSL C Math Library

[15] Chapter 11. Quadrature and Differentiation - 1.69.0

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

2019年8月28日 (水) 00:47時点における版編集 MathXYZ (会話 \| 投稿記録) 2,902 回編集 m →乱数を使わない積分 ← 古い編集		2019年9月1日 (日) 12:54時点における版編集取り消し MathXYZ (会話 \| 投稿記録) 2,902 回編集リンクを追加タグ: モバイル編集モバイルアプリ編集 Androidアプリ編集新しい編集 →
1行目:		1行目:
	'''数値積分'''（すうちせきぶん、{{lang-en-short\|numerical integration}}）とは、狭義には与えられる関数の[[定積分]]の値を、解析的にではなく数値的に求める[[求積法]]のことであり、広義には与えられる[[導関数]]から原関数を求めること、また[[微分方程式]]を数値的に解くことを含む<ref name="Yamamoto1">{{Cite book \|和書 \|author=山本哲朗 \|title=数値解析入門 \|edition=増訂版 \|date=2003-06 \|publisher=[[サイエンス社]] \|series=サイエンスライブラリ現代数学への入門 14 \|ISBN=4-7819-1038-6}}</ref>。[[数値解析]]の一分野である。		'''数値積分'''（すうちせきぶん、{{lang-en-short\|numerical integration}}）とは、狭義には与えられる関数の[[定積分]]の値を、解析的にではなく数値的に求める[[求積法]]のことであり、広義には与えられる[[導関数]]から原関数を求める手法、また[[微分方程式]]を数値的に解く手法 ([[常微分方程式の数値解法]]、[[偏微分方程式の数値解法]]) を含む<ref name="Yamamoto1">{{Cite book \|和書 \|author=山本哲朗 \|title=数値解析入門 \|edition=増訂版 \|date=2003-06 \|publisher=[[サイエンス社]] \|series=サイエンスライブラリ現代数学への入門 14 \|ISBN=4-7819-1038-6}}</ref>。[[数値解析]]の一分野である。

	狭義の数値積分（関数の定積分の値を求める方法）は[[有限要素法]]などで応用されている<ref name="Yamamoto1"/>。		狭義の数値積分（関数の定積分の値を求める方法）は[[有限要素法]]などで応用されている<ref name="Yamamoto1"/>。