2017年11月9日木曜日

ホイン(Heun)法

ホイン(Heun)法


§1 オイラー法

微分方程式
の解は
である。
(2)式の右辺の積分を
で近似すると、(2)式は
と近似することにより、x₀を計算の起点とし、x₁x₂、・・・と築時的に計算をすることができる。
とあらわすことにすると、
になる。
とくに、点が等間隔hに並んでいるとき、
である。
この方法をオイラー法という。

問題 オイラー法を用いて、次の微分方程式のx=0.10.20.3におけるyの値を求めよ。
【解】
x₀=0y₀=y(x₀)=1h=0.1とおくと、
x₁=0.1y₁=1.05だから
x₂=0.2y₂=1.163だから
したがって、
x=0.1のとき、y=1.05
x=0.2のとき、y=1.1106
x=0.3のとき、y=1.1846
(解答終)

表計算ソフトを用いた計算結果を以下に示す。





オイラー法だと、計算を進めると、誤差が蓄積して、その結果、厳密解
との食い違いが大きくなっていくことが理解できると思う。


§2 ホイン法

まず、オイラー法
を用いてを計算し、次に台形公式
と修正する。
以下同様に、所定の精度が得られるまで
と反復計算をする方法をホイン法という。

ホイン法を用いて、問題の微分方程式をh=0.1と同一の条件で解いた結果は次の通り。


0 件のコメント:

コメントを投稿