next up previous
Next: 5 レポート課題 Up: 計算天文学 II 第2回 安定性と陰解法 Previous: 3 陰解法

4 高精度の方法

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

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

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

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



Jun Makino
Thu Oct 19 14:00:54 JST 2000