ニュートン・コーツの公式

ニュートン・コーツの公式(ニュートン・コーツのこうしき、: Newton–Cotes formulae, Newton-Cotes rules)とは、等間隔の点における被積分関数の値に基づく数値積分法の総称である。名前はアイザック・ニュートンロジャー・コーツに由来する。

ニュートン・コーツの公式は、等間隔の点での被積分関数の値が与えられた場合に有用である。もし他の点での値も求められるならば、ガウス求積クレンショー・カーチス求積英語版などの他の方法の方が適している場合もある。

概要

編集

ニュートン・コーツの公式は、端点を使う「閉じた」ものと、端点を使わない「開いた」ものの 2 種類に大別できる。

n 次の閉じたニュートン・コーツの公式は次のようになる。

 

ここで     である。

wi は重みと呼ばれる。重みは以下のようにラグランジュ補間による補間多項式から導かれる。

 

また、以上の導出から重みは関数 f によらず、xi のみによって決まることがわかる。

n 次の開いたニュートン・コーツの公式の場合は、    とし、重みは閉じたものと同様である。

ニュートン・コーツの公式の一覧

編集
閉じたニュートン・コーツの公式
次数 名前 誤差項
1 台形公式    
2 シンプソンの公式    
3 シンプソンの3/8公式    
4 ブールの公式英語版    
開いたニュートン・コーツの公式
次数 名前 誤差項
0 中点則    
1 台形法    
2 ミルンの公式    
3    

ここで、fif (xi) の略記である。

誤差項 E  となる ξ ∈ (a, b) が存在することを意味する。 また、f の導関数の次数は、それ未満の次数の多項式が正確に積分できる(即ち、誤差が 0 になる)ことを示している。なお、(ba) の次数と f の導関数の階数は、1 つおきに 2 ずつ増加することに注意。

重みの計算

編集

ニュートン・コーツの公式の重みは線形方程式系の解として求めることもできる。 これは補間多項式の一意性より f (x)n 次以下の多項式の場合 L(x) = f (x) となることに基づく。係数行列ファンデルモンド行列である。

 

高次における不安定性

編集
 
ルンゲ関数 f(x) = 1/(1 + 25x2) (赤)と等分点に基づく5次(青)と9次(緑)の補間多項式。高次多項式の方が端点付近での誤差は大きくなっている。

ニュートン・コーツの公式は、任意の次数で構築できる。しかし大きな次数 n においてはルンゲ現象により誤差n の増加するにつれて指数関数的に大きくなる。そのため、通常は大きな次数ではガウス求積クレンショー・カーチス求積英語版などの非等分点法の方が、安定してより正確な値を求められる。もしもそれらの方法を使えないならば、合成積分公式を使うことでルンゲ現象を避けることができる。高次の公式には積分の重みの中に負のものが含まれるなどの不自然さが伴う(Gaussの積分公式の重みは常に正である)。

合成積分公式

編集

ニュートン・コーツの公式の精度を良くするには、ステップ長 ba/n は小さくする必要がある。つまり、積分区間 [a, b] 自体が小さくなければならない。このため、積分区間 [a, b] を小さな部分区間に分割し、各部分区間ごとにニュートン・コーツの公式を使い、その結果を足し合わせるという方法が使われる。これは合成積分公式と呼ばれる。

関連項目

編集

参考文献

編集
  • Abramowitz, M.; Stegun, I. A. (1972). “Section 25.4”. Handbook of Mathematical Functions with Formulae, Graphs, and Mathematical Tables. New York: Dover 
  • Forsythe, George E.; Malcolm, Michael A.; Moler, Cleve B. (1977). “Section 5.1”. Computer Methods for Mathematical Computations. Englewood Cliffs, NJ: Prentice–Hall 
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007), “Section 4.1. Classical Formulas for Equally Spaced Abscissas”, Numerical Recipes: The Art of Scientific Computing (3rd ed.), New York: Cambridge University Press, ISBN 978-0-521-88068-8, http://apps.nrbook.com/empanel/index.html?pg=156 
  • Stoer, Josef; Bulirsch, Roland (1980). “Section 3.1”. Introduction to Numerical Analysis. New York: Springer-Verlag 

外部リンク

編集