数値積分は、積分を解析的に解くことが難しい場合に役立つ方法です。特にSimpson則や台形則は、数値積分でよく使用される手法ですが、これらの方法が解析解とどのように比較されるのかについて、具体例を使って解説します。この記事では、関数f(x)=1/xを対象に、解析解、Simpson則、台形則を比較し、得られる誤差の原因を考察します。
1. 数値積分とは?
数値積分は、積分を解析的に解く代わりに、数値的に解く手法です。特に、複雑な関数や解析解が得られない場合に、近似的な解を求めるために使用されます。一般的な数値積分の手法には、台形則、Simpson則、リーマン和などがあり、それぞれが異なる精度を持っています。
2. 解析解の計算
関数f(x) = 1/xの積分を1≦x≦2の範囲で求めると、解析解はF(x) = log|x|となります。範囲1から2までの積分値は、log2 – log1であり、計算すると約0.693147180559945となります。これはEXCELのLN(2)関数で求めることができる値です。
3. Simpson則と台形則による数値積分
Simpson則と台形則は、どちらも数値積分の代表的な方法です。Simpson則はより高精度な近似を提供することが知られています。例として、n=100で計算した場合、Simpson則の結果は0.693147180872367であり、解析解との差は非常に小さい(誤差がほぼゼロ)です。
一方、台形則で同様に計算した場合、結果は0.693153430481824となり、解析解との差は0.000006249921878809240000000000となり、Simpson則の誤差よりも大きいことがわかります。
4. なぜSimpson則や台形則で誤差が発生するのか
数値積分における誤差は、主に近似方法の精度に依存します。Simpson則は、関数を2次関数で近似するため、台形則よりも高精度な結果を得ることができます。しかし、数値積分は常に近似であり、厳密な解析解との差は誤差として現れます。
5. まとめ:Simpson則と台形則の比較と精度向上の方法
Simpson則と台形則は、数値積分において重要な手法ですが、それぞれの精度には差があります。特に、Simpson則は台形則よりも精度が高いですが、nの値を増やすことでさらに精度が向上します。実際にnを1000、10000、100000に増やすことで、誤差はさらに小さくなり、より正確な結果を得ることができます。


コメント