中点法
応用数学、特に数値計算における中点法(ちゅうてんほう、英: midpoint method)は、一段階法(one-step method)による常微分方程式
- {\displaystyle y'(t)=f(t,y(t)),\quad y(t_{0})=y_{0}}
の数値解析的な解法の一つである。
陽的中点法(explicit midpoint method)は以下の公式に基づく。
- {\displaystyle y_{n+1}=y_{n}+hf\left(t_{n}+{\frac {h}{2}},y_{n}+{\frac {h}{2}}f(t_{n},y_{n})\right)\qquad \qquad (1e)}
陰的中点法(implicit midpoint method)は以下の公式に基づく。
- {\displaystyle y_{n+1}=y_{n}+hf\left(t_{n}+{\frac {h}{2}},{\frac {1}{2}}(y_{n}+y_{n+1})\right)\qquad \qquad (1i)}
ここで {\displaystyle n=0,1,2,\dots }、{\displaystyle h} はステップ幅(step size)と呼ばれる微小量、{\displaystyle t_{n}=t_{0}+nh}、また {\displaystyle y_{n}} は算出された {\displaystyle y(t_{n})} の近似値。陽的中点法は修正オイラー法(modified Euler method)としても知られる[1] 。陰的中点法は最も単純な選点法であり、ハミルトン力学におけるシンプレクティック積分法に応用される。
この手法の名称は、上掲の公式において傾きを与える関数 {\displaystyle f} の値が、{\displaystyle y(t)} が既知の点 {\displaystyle t_{n}} と未知の(これから算出したい)点 {\displaystyle t_{n+1}} との中点 {\displaystyle t=t_{n}+h/2={\frac {t_{n}+t_{n+1}}{2}}} において計算されることに由来する。
幾何学的に解釈することで直感的な理解がしやすくなるかもしれない(上掲の図を参照)。基本的なオイラー法では、曲線の {\displaystyle (t_{n},y_{n})} での接線を {\displaystyle f(t_{n},y_{n})} から計算し、次の値 {\displaystyle y_{n+1}} を接線と垂線 {\displaystyle t=t_{n+1}} との交点から求める。しかし、2階微分が {\displaystyle t_{n}} と {\displaystyle t_{n+1}} の間で常に正、または(上掲の図のように)常に負であれば、曲線は接線から逸れ続けて {\displaystyle h} の増大とともに誤差は拡大する。この図からは、この区間においては中点における接線(上側の緑色の線分)のほうが、ほとんどいつもより正確な曲線の近似となることが見てとれる。この接線は、まず元のオイラー法を使って中点での {\displaystyle y(t)} を推定し、次に {\displaystyle f} から傾きを計算することで求まる。この改良された接線を使って {\displaystyle y_{n}} から {\displaystyle y_{n+1}} の値が計算される(図中の赤色の弦)。中点での {\displaystyle y(t)} の推定に誤差が伴うため、赤色の弦は正確には緑色の線分(真の接線)と平行ではないことに注意する。
中点法における各ステップでの局所的な誤差は {\displaystyle O\left(h^{3}\right)} であり、大域的な誤差は {\displaystyle O\left(h^{2}\right)} となる。このように、中点法はオイラー法と比べて数値計算的に強力でありながら、{\displaystyle h\to 0} のときの誤差の収束は一般的により速い。
中点法は、ルンゲ=クッタ法として知られるより高次の手法の一クラスである。
導出
[編集 ]中点法はオイラー法
- {\displaystyle y_{n+1}=y_{n}+hf(t_{n},y_{n})}
を改良したものであり、手順も類似している。 オイラー法を導く鍵は近似式
- {\displaystyle y(t+h)\approx y(t)+hf(t,y(t))\qquad \qquad (2)}
であり、これは傾きの公式
- {\displaystyle y'(t)\approx {\frac {y(t+h)-y(t)}{h}}\qquad \qquad (3)}
と {\displaystyle y'=f(t,y)} に注意することで得られる。
中点法の場合、(3) をより正確な
- {\displaystyle y'\left(t+{\frac {h}{2}}\right)\approx {\frac {y(t+h)-y(t)}{h}}}
で置き換える。すると (2) は次の式に変わる。
- {\displaystyle y(t+h)\approx y(t)+hf\left(t+{\frac {h}{2}},y\left(t+{\frac {h}{2}}\right)\right)\qquad \qquad (4)}
この式を使って {\displaystyle y(t+h)} を求めることはできない。なぜなら {\displaystyle y} の {\displaystyle t+h/2} での値が分からないからである。そのため、ちょうどオイラー法で {\displaystyle y(t+h/2)} を求めるのと同じように、テイラー展開を使って
- {\displaystyle y\left(t+{\frac {h}{2}}\right)\approx y(t)+{\frac {h}{2}}y'(t)=y(t)+{\frac {h}{2}}f(t,y(t))}
と近似する。これを (4) に代入すると
- {\displaystyle y(t+h)\approx y(t)+hf\left(t+{\frac {h}{2}},y(t)+{\frac {h}{2}}f(t,y(t))\right)}
のように陽的中点法 (1e) の公式が得られる。
陰的中点法 (1i) は、ステップの中間点 {\displaystyle t+h/2} での値を {\displaystyle y(t)} と {\displaystyle y(t+h)} とを結ぶ線分の中点
- {\displaystyle y\left(t+{\frac {h}{2}}\right)\approx {\frac {1}{2}}{\bigl (}y(t)+y(t+h){\bigr )}}
で近似し、
- {\displaystyle {\frac {y(t+h)-y(t)}{h}}\approx y'\left(t+{\frac {h}{2}}\right)\approx k=f\left(t+{\frac {h}{2}},{\frac {1}{2}}{\bigl (}y(t)+y(t+h){\bigr )}\right)}
とすることで得られる。{\displaystyle y_{n}+h,円k} で {\displaystyle y(t_{n}+h)} を近似すると陰的ルンゲ=クッタ法
- {\displaystyle {\begin{aligned}k&=f\left(t_{n}+{\frac {h}{2}},y_{n}+{\frac {h}{2}}k\right)\\y_{n+1}&=y_{n}+h,円k\end{aligned}}}
が得られるが、前半部分はステップ幅を {\displaystyle h/2} とした後退オイラー法を含んでいる。
陰的中点法では時間についての対称性から {\displaystyle h} の偶数乗の局所誤差項は全てキャンセルし、誤差は自動的に {\displaystyle O\left(h^{3}\right)} となる。{\displaystyle k} の決定を後退オイラー法ではなく(前進)オイラー法で行うと、再び陽的中点法が得られる。
関連項目
[編集 ]脚注
[編集 ]- ^ Süli & Mayers 2003, p. 328
参考文献
[編集 ]- Griffiths,D. V.; Smith, I. M. (1991). Numerical methods for engineers: a programming approach. Boca Raton: CRC Press. p. 218. ISBN 0-8493-8610-1
- Süli, Endre; Mayers, David (2003), An Introduction to Numerical Analysis, Cambridge University Press, ISBN 0-521-00794-1 .