next up previous
Next: 6 練習 Up: 計算天文学 II 第4回 安定性と陰解法 Previous: 4 陰解法

5 高精度の方法

さて、前節の簡単な陰解法は無条件安定で大変結構なのであるが、時間刻みを 大きくとれるとなると、今度は計算精度のほうが気になってくる。空間方向の 微分は2次精度で、テイラー展開の2次の項までが正しく入っているのに対し、 時間微分は1次精度にしかなっていないからである。

なお、上の1次の陽解法のことを1次の前進差分、陰解法のことを1次の後退差 分ということがある。

時間方向の精度をあげる一つの方法は、これまでに見てきた2つの方法を混ぜ て使うことである。つまり、

この方法をクランク・ニコルソン法という。これは時間方向に対称な形になっ ている。この方法では、時間方向も2次精度にすることができる。



Jun Makino
Sun Oct 28 18:39:25 JST 2001