http://www.mumumu.org/~pascal/kaguya/kaguya_energy.png
データ見てみようかなと思ってちょっとだけ遊んだよ。
グラフにしてみたのは、「1kgあたりの地球をベースにした重力エネルギー」と「1kgあたりの運動エネルギー」の時間変化。
かぐやの重さはわからないので1kgあたり。
目的は、もしかしたら軌道修正なんかのマヌーバを行った時間と軌道変化が見られるかも…というものだったけれど、もともとのエネルギーがすごいこともあってわかるとは言えない。

あ、参考までにかぐやの軌道修正に関係しそうな日時も。
http://www.selene.jaxa.jp/ja/event/timeline_j.htm

グラフに関しては以下のとおり。
・時間はJSTにしてある…つもりだったけどUTCだ。あうあう。
・最初をふくめて3回突き抜けてるのは、近地点で運動エネルギーがものすごいことになってるとき。
・さすがにマヌーバなんかが読み取れるほどの変化はなかった。というか運動エネルギーがもともと大きすぎるのでどうやっても埋もれる。
・9月19日の午前0時56分(UTC)、エネルギー総量が急激に変化している(-1.62e+06[J/kg]から-1.02+06[J/kg])。原因は？
　・19日には周期調整マヌーバが予定されていた
　　・そんなにすごい軌道変更があったわけじゃないし、そもそも自発的にこんなに大きなエネルギーは出せない…んじゃないかな？
　・0時56分は、ちょうどログファイルが別のものになった時間
　　・2つのログファイルの両方に0時56分のデータがある
　　　・位置情報はほぼ同じ値。
　　　・速度情報がかなり異なる
　　・値がおかしくて別の計測器に切り替えたとか？
　・教えてえらい人
・10月3日の夜からエネルギー総量が変化
　・月に捕まったため
　　・このグラフでは月の重力エネルギーを考慮していないもんね
　　　・だって自分で月の軌道をこの座標系に入れるの面倒だし…
　・月周回の円軌道に入っても運動エネルギーが振動
　　・地球を基準とした運動エネルギーなので当たり前
・月に投入してからは位置エネルギーは安定
　・確かに、遠地点で軌道修正して月に投入したしたことがわかる
　　・エネルギーが一番小さくてすむらしいね。何軌道っていうんだっけ。
　　　・ホーマン軌道、だ。

他に「位置と速度ベクトルを用いてdtずつ時間を進めてみるっていうシミュレーション」と実際の位置を比較して軌道修正の有無を見てみようーってのもやってみたんだけど、なんとなくそれっぽい日時は出たけど定量的じゃないのでまだ載せない。

月の位置の関数を知りたいんだけどどうするのが一番いいのかね。かぐやが月に捕まって以降の軌跡から月の存在する平面はわかって、…っていう方法かねぇ。あとは、エネルギー総量が変わらないようにかぐやと動きを合わせる…いいのかそれ。