ITPASS 実習レポート 2 ～数値計算実習課題～

     名前 : 藤田 哲也

担当情報実験機 : joho13

その１．

< tex ファイル >

修正後

tex ファイルのダウンロード

txt 形式での表示

修正前

tex ファイルのダウンロード

txt 形式での表示

< pdf ファイル >

修正後

修正前

< 工夫したこと, 感想 >

< 参考資料 >

 数式の書き方について参考にした。

・ 土屋勝 著 ｢LATEX はじめの一歩｣ カットシステム 2007

 数式の書き方や文書の構成の方法について参考にした。

・ C. D. Murray and S. F. Dermott 著 ｢Solar System Dynamics｣ Cambridge University Press 2000

 運動方程式の導出とその意味合いについて参考にした。

・ 福島登志夫 著 ｢シリーズ現代の天文学13 天体の位置と運動｣ 日本評論者 2009

 運動方程式の導出とその意味合いについて参考にした。

< 共同作業した人の名前 >

 坂本 大樹、赤松 紀美

その２．

< 実習の目的 >

• 三体問題における惑星の運動について、fortran を用いて数値計算を行う。
• 求めた値を gnuplot でグラフとして表示させ、惑星がどのような軌道を描くか見る。
• 実習およびレポート作成を行うなかで、fortran や gnuplot の扱いに慣れる。

< まとめ・考察 >

　表の解説でも書いたように、中心星に対して惑星の質量が小さく距離が遠いほど軌道は安定になり、質量が大きく距離が近いほど軌道は不安定になる。二体問題では、惑星をプロットした点において、その位置での速度が初期速度と等しい安定な楕円軌道を描くが、三体問題においてはもう一つの惑星の影響があるために複雑な、もしくは不安定な軌道となる。今回は三体だけだったが、それでも太陽系の質量の大部分を占める三つについて計算していることを思えば、最初に xp(3)=1.873 , factor=1.0 で計算を行うことで描かれた図は、実際の三体の運動に近いものになっていると考えられる。他の惑星の質量や距離を考えれば、それらによる影響は小さな摂動として現れるだろう。

< 工夫したこと, 感想 >

　実習中で用いた値とは違う値を用いて、安定と不安定の境界が出るようにしてみた。惑星の様々な軌道が見られるのは楽しく、より複雑な状況や N 体にしたときにどうなるのかということにも興味が出てくる。

< 共同作業した人の名前 >

 坂本 大樹　赤松 紀美

< 参考資料 >

 gnuplot でグラフを保存する方法について参考にした。

 html で表を作成する方法について参考にした。