ロケット観望シミュレーション

というわけで、ロケット打上の見え方を計算するExcel書類を作ってみました。
「Tracking.xls」をダウンロード
観測地の緯度と経度を入力すれば、任意の地点からの見え方をプロットできます。計算にあたって、機体の緯度・経度・高度の３次元情報が欲しいのですが、あいにく今回の７号機で十分な情報は見つけられませんでした。代わりに、似た経路で打ち上げた３号機のデータを入れてみました。参考程度にはなるでしょう。
鹿児島、大阪、東京の計算結果をお見せします。途中の曲線がいびつなのはデータ点が少ないためです。ほんとうはもっと滑らかになります。また、グラフの縦横の縮尺は自動にしていて、横の方が縮められています。実際の見た目の形状はもっと横長になります。
オレンジの線がブースター燃焼中。炎は明るく、白煙が残ります。
太い白線がその後の第１段の燃焼中。排気は透明に近く、彗星のように見えます。
細い白線が２段目の飛行。自分で見ていないのですが、たぶん暗いでしょう。

図中の用語の説明
SRB：固体ロケットブースター（大きい方）
MECO：第１段主エンジン停止（１段目分離はこの直後）
SECO1：第２段エンジン１回目停止
方位角：東が90度、南が180度、西が270度
仰角：水平線が0度、真上が90度、マイナスは水平線下で見えません。

３号機と７号機の主な違いについて、以下に留意してご覧下さい。
（１）打ち上げた方位が、３号機は方位角90度（真東）、７号機は102度（やや南）なので、３号機の方が若干本土の近くを飛行し、見やすくなっている。
（２）３号機は固体補助ロケット(SSB)が４本で、７号機は２本。
（３）SRBの分離が３号機の方がだいぶ早い。
（４）３号機は衛星２基搭載で、上部フェアリングと下部フェアリングがある。図中「フェアリング」と書いたのは上部フェアリングのほうで、このタイミングは７号機のフェアリング分離とだいたい同じ。

【技術的解説】
３号機の「追跡管制計画書」には、各イベント時の飛行距離と高度の表があります。このデータをもとに以下の仮定をして、球面三角法で緯度・経度に換算します。
　仮定１：打ち上げた方位から変針しない→軌道面内のみの運動とする
　仮定２：打上直後の短時間に限る→地球の自転は無視
　仮定３：地球は完全な球体モデルで偏平性は無視
　仮定４：データが少ない区間は直線で補間
それから、観測地の局所水平面を基準にした［東，北，上］のデカルト座標での相対位置に変換。最後に、相対位置座標値から、方位角・仰角を決定してプロット。方位角の0度と360度をまたぐ場合の処理は省略しているので、飛行経路の南側（小笠原など）から見ると変になります。

飛行データのシートを更新すれば、他の打上にも使えるので、今後のロケット見物の支援用にもいかが？
距離・高度を３号機のデータのままでも、打上方位を102度にすれば、７号機に近くなるでしょう。南に向けて打てばALOSやIGSに近くなるでしょう。