next up previous
Next: 5 次週予告 Up: 計算天文学 II 第5回 常微分方程式の初期値問題(2) Previous: 3 ハミルトン系向けの方法

4 レポート課題

1は必須。後はオプション。 プログラムが必要なものはプログラムを提出 すること。

  1. 2次元ケプラー問題

    について、古典ルンゲクッタと leap frog の両方で、適当な初期条件(解が 楕円になるものをとること)から出発して 10, 1000, 周期後の解の厳 密解からの誤差が刻み幅のどのような関数になるか調べよ。ここでは刻み幅は 固定とする。

  2. 4次のシンプレクティック公式は、 が刻み幅 h の leap frog を表す作用素であるとして、以下のように書ける。

    つまり、リープフロッグでまずちょっと行き過 ぎて、次に逆に戻って、最後に最初と同じだけ進んで次の時刻に行くというも のである。これをプログラムし、 1 と同様な解析を行なえ。

  3. 4次のエルミート公式についても同様な解析を行なえ。説明されている ものは陰的公式なので、予測子に何を使えるかは各自検討すること。

  4. 2.1 で述べた古典的 RK 法で誤差の調整をする方法をプログラムし、離 心率の大きなケプラー問題について刻み一定の場合にくらべてどの程度得にな るか調べてみよ。



Jun Makino
Wed Dec 6 20:43:19 JST 2000