ロンバーグ積分

ロンバーグ積分（ロンバーグせきぶん、Romberg integration）またはロンベルク積分は、関数の数値積分アルゴリズムのひとつである。この方法では台形公式とリチャードソンの補外を組み合わせ離散化幅 $h$ をゼロとする極限として数値積分を評価する。他の数値積分法に比べ、少ない回数の被積分関数の評価によって高精度の結果が得られる。

1955年にヴェルナー・ロンベルクによって考案された。

概要

ロンバーグ積分は区間 $[a,b]$ で定義された関数 $f$ の定積分

$\int _{a}^{b}f(x)dx$

を数値的に求めるアルゴリズムである。積分区間 $[a,b]$ を幅 $h=(b-a)/N$ の $N$ 個の小区間 $[x_{i},x_{i+1}]$ ($x_{i}=a+ih$, $i=0,1,\cdots ,N-1$) に分割するとき, 数値積分アルゴリズムのひとつである台形公式は求める定積分を

$T(h)=h\left[{\frac {1}{2}}f(x_{0})+f(x_{1})+\cdots +f(x_{N-1})+{\frac {1}{2}}f(x_{N})\right]$

により近似する。台形公式による値 $T(h)$ と真の定積分の値の差は ${\mathcal {O}}(h^{2})$ であるが、オイラー・マクローリンの和公式（オイラーの和公式）によりそれをより詳しく評価することができる。すなわち、$f$ が $C^{2m+2}$-級関数であるならば

$T(h)=\tau _{0}+\tau _{1}h^{2}+\tau _{2}h^{4}+\cdots +\tau _{m}h^{2m}+\alpha _{m+1}(h)h^{2m+2}$

$\tau _{0}=\int _{a}^{b}f(x)dx$

$\tau _{k}={\frac {B_{2k}}{(2k)!}}\left[f^{(2k-1)}(a)-f^{(2k-1)}(b)\right]\ \ (k=1,2,\cdots ,m)$

$\alpha _{m+1}(h)={\frac {B_{2m+2}}{(2m+2)!}}(b-a)f^{(2m-2)}(\xi (h))$

が成立する。ここに $B_{2k}$ はベルヌーイ数 であり、$\xi (h)$ は $a<\xi (h)<b$ を満足するある数である。それ故に台形公式を $h$ に関する漸近級数と見るとき、誤差項としては $h$ の偶数次の項だけが現れる。そこで複数の $h$ について台形公式を適用し、それらの結果について誤差項を打ち消すようにリチャードソンの補外を行うことで、より精度の高い定積分の評価が得られる。これがロンバーグ積分である。

アルゴリズム

ロンバーグ積分のアルゴリズムは以下のように記述される。

1. いくつかの幅 $h_{0}$, $h_{1}=h_{0}/N_{1}$, ..., $h_{m}=h_{0}/N_{m}$ に対して台形公式を適用することにより定積分の近似値 $T_{0,0}=T(h_{0})$, ..., $T_{m,0}=T(h_{m})$ を求める.
2. 次の漸化式（ネヴィルのアルゴリズム）
   $\quad T_{i,k}=T_{i,k-1}+{\frac {T_{i,k-1}-T_{i-1,k-1}}{\left[{\frac {h_{i-k}}{h_{i}}}\right]^{2}-1}}$
   により帰納的に $T_{i,k}$ を計算し、得られた $T_{m,m}$ が求める補外値を与える。

幅 $h_{0},\cdots ,h_{m}$ としては、ロンバーグは

$h_{0}=b-a,\ \ h_{1}={\frac {h_{0}}{2}},\ \ \cdots ,h_{i}={\frac {h_{i-1}}{2}},\cdots$

を用いていたが、Roland Bulirschは次数が急激に増加し計算コストが増大することを抑えるために

$h_{0}=b-a,\ \ h_{1}={\frac {h_{0}}{2}}h_{2}={\frac {h_{0}}{3}},\ \ \cdots ,\ \ h_{i}={\frac {h_{i-2}}{2}},\cdots$

を提案している。

補外ステップにおいてリチャードソン補外テーブルの各項の差分に基づいて積分の収束性を判定することで、計算をどこまで進めればよいか判定することも可能である。

特徴

ロンバーグ積分 $T_{i,k}$ の誤差は、評価に用いたステップ $h_{i-k},\cdots ,h_{i}$ から ${\mathcal {O}}(h_{i-k}^{2}h_{i-k+1}^{2}\cdots h_{i}^{2})$ と見積もられる。特に $h_{i}$ をロンバーグの提案に従って選ぶときには誤差の厳密な評価

$T_{i,k}-\int _{a}^{b}f(x)dx=(b-a)h_{i-k}^{2}h_{i-k+1}^{2}\cdots h_{i}^{2}{\frac {(-1)^{k}B_{2k+2}}{(2k+2)!}}f^{(2k+2)}(\xi )$

が得られている。

$h_{2}=h_{1}/2$ に選ぶとき、リチャードソン補外において現れる $T_{2,1}$ や $T_{3,1}$ はシンプソンの公式によって問題の定積分を評価した値に一致する。その意味で、より多くの積分幅 $h$ での評価を行うことはニュートン・コーツの公式において積分公式を拡張することに、漸化式を解き補外を行うことは積分の次数を上げることに対応していると解釈できる。

Press らは著書「Numerical Recipes in C」においてロンバーグ積分が異なる次数の数値積分法を組み合わせた一般的な方法であることから数値積分においてこの方法を用いることを強く推奨している。特に、被積分関数が積分区間において十分に滑らかであれば、ロンバーグ積分はニュートン・コーツ型の積分公式に比べて極めて少ない回数の被積分関数の評価だけで十分な精度の結果を与える。

なお、ネヴィルのアルゴリズムはラグランジュ多項式による補外であるが、有理関数補外を行うこともでき、この場合ロンバーグ積分における漸化式は

$\quad T_{i,k}=T_{i,k-1}+{\frac {T_{i,k-1}-T_{i-1,k-1}}{\left[{\frac {h_{i-k}}{h_{i}}}\right]^{2}\left[1-{\frac {T_{i,k-1}-T_{i-1,k-1}}{T_{i,k-1}-T_{i-1,k-2}}}\right]-1}}$

へと変更される。

具体例

次の定積分をロンバーグ積分により求めることを考える。

$\int _{0}^{1}x^{5}dx={\frac {1}{6}}$

積分ステップとして $h_{0}=1$, $h_{1}=1/2$, $h_{2}=1/4$ を取り補外を行う。その結果、次の表に示すようにロンバーグ積分によって正確な積分値を得ることができる。

ロンバーグ積分におけるリチャードソン補外

$h_{0}=1$	$T_{0,0}=0.500000={\frac {1}{2}}$
$h_{1}=1/2$	$T_{1,0}=0.265625={\frac {17}{64}}$	$T_{1,1}=0.187500={\frac {3}{16}}$
$h_{2}=1/4$	$T_{2,0}=0.192383={\frac {197}{1024}}$	$T_{2,1}=0.167969={\frac {43}{256}}$	$T_{2,2}=0.166667={\frac {1}{6}}$

ライブラリ

C言語により実装されているオープンソースの科学技術計算ライブラリGNU Scientific Library (GSL) にはロンバーグ積分を行う関数が含まれている。また、Pythonの数値解析ライブラリSciPyにはロンバーグ積分により数値積分を行う scipy.integrate.romberg および scipy.integrate.romb というふたつの関数がある.

参考文献

• Press, William H.; Teukolsky, Saul A.; Vetterling, William T.; Flannery, Brian P. (1992). Numerical Recipes in C: The Art of Scientific Computing (2nd ed.). Cambridge University Press. doi:10.2277/0521431085. ISBN 978-0-521-43108-8 
• Stoer, Josef; Bulirsch, R. (2002). Introduction to Numerical Analysis. Springer. doi:10.1007/978-0-387-21738-3. ISBN 978-0-387-21738-3 

関連項目

• 台形公式
• シンプソンの公式
• ニュートン・コーツの公式
• ブリルシュ・ストア法 - 常微分方程式の数値解法のひとつで、修正中点法の結果に対してロンバーグ積分と同様にリチャードソンの補外を行う。