› ケプラーの方程式 シミュレーター 戻る
天体力学シミュレーター

ケプラーの方程式 シミュレーター — 楕円軌道の数値解

超越方程式 M = E − e·sin(E) を Newton 法でリアルタイムに解き、離心近点角・真近点角・焦点距離・公転周期を表示。軌道形状と二次収束の様子を直感的に学べます。

パラメータ設定
半長径 a
AU
離心率 e
平均近点角 M
deg
中心質量 M_central
M_sun

許容誤差 1e-10、最大 30 反復で収束判定。Newton 法は二次収束のため、通常 3〜5 回で 10 桁の精度に達します。

計算結果
—
離心近点角 E
—
真近点角 ν
—
焦点からの距離 r
—
公転周期 T
楕円軌道と衛星位置

青楕円＝軌道／黄丸＝中心天体（焦点）／白丸＝衛星位置／緑線＝焦点-衛星間距離 r／灰破線＝離心近点角 E に対応する補助円上の点

Newton 反復の収束履歴

横軸＝反復回数 n／縦軸＝残差 log10|f(E_n)|／二次収束のため反復ごとに桁数がほぼ倍増するのが特徴です（許容誤差ライン 1e-10）。

理論・主要公式

ケプラーの方程式は、楕円軌道上の天体位置を時刻から求めるための超越方程式です。平均近点角 M は時間に比例しますが、実際の位置を与える離心近点角 E は M と非線形に結びついており、数値解が必要です。

ケプラーの方程式と Newton 反復：

$$M = E - e\,\sin E,\qquad E_{n+1} = E_n - \frac{E_n - e\,\sin E_n - M}{1 - e\,\cos E_n}$$

真近点角 ν と焦点からの距離 r：

$$\tan\!\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}}\,\tan\!\frac{E}{2},\qquad r = a\,(1 - e\,\cos E)$$

ケプラーの第三法則による公転周期（中心質量は太陽質量単位）：

$$T = 2\pi\sqrt{\frac{a^{3}}{G\,M_{c}}}\;\;\Longrightarrow\;\; T_{\mathrm{yr}} = \sqrt{\frac{a_{\mathrm{AU}}^{3}}{M_{c,\odot}}}$$

$a$ は半長径 [AU]、$e$ は離心率（0〜1）、$M$ は平均近点角 [rad]、$E$ は離心近点角 [rad]、$\nu$ は真近点角 [rad]、$M_c$ は中心質量 [太陽質量]、$T$ は公転周期 [年] です。

ケプラーの方程式 シミュレーターとは

🙋
惑星って楕円軌道上を動くって習ったんですけど、時刻を入れると位置が出てくる、みたいな公式ってあるんですか？
🎓
惜しい——「公式」と言いたいところだけど、実は解析的に解けないのがミソだ。時刻に比例する量を平均近点角 M とすると、楕円中心から測った離心近点角 E と $M = E - e\sin E$ で結ばれる。これがケプラーの方程式で、初等関数では E について解けないので Newton 法で数値的に解く。本シミュレーターはまさにそれをやっている。
🙋
えっ、惑星の位置すら解析的に出せないんですか？
🎓
そうなんだ。だから初学者は驚くんだけど、これは「超越方程式」の典型例として今でも数値計算の教科書冒頭に出てくる。デフォルト値 e=0.20, M=60° なら、E_0 = 60° から始めて Newton 反復が E_1 ≈ 71.03°, E_2 ≈ 70.823° と二次収束して、3 反復で 10 桁の精度に達する。右下の収束グラフで残差が一気に下がっていくのが見えるはずだ。
🙋
離心率を 0.9 とかに上げてみたらどうなりますか？
🎓
いい着眼だ。e を 0.9 に上げると軌道が極端に細長くなって、特に M が π 付近のとき初期値 E_0 = M がよくない近似になって反復が増える。e=0.967 のハレー彗星クラスだと専用の初期値選択（Battin の式など）が必要になることもある。本ツールでは安全のため最大 30 反復、許容誤差 1e-10 にしてあって、e=0.95 でもちゃんと収束する。スライダーで e を変えながら収束履歴グラフを見ると、二次収束が見える条件の違いがよくわかる。
🙋
右上の真近点角 ν は、E とどう違うんですか？
🎓
いい質問だ。E は楕円中心から測った「幾何学的な角度」、ν は焦点（太陽）から測った「実際に観測される角度」だ。デフォルト e=0.20, M=60° では E=70.82°, ν=82.10° と少しずれる。離心率 0 の円軌道なら M=E=ν で全部一致するけど、楕円ではこの 3 つの角度の違いがいわば「軌道のアイデンティティ」と言える。位置を観測量に変換するときは ν の方を使う。
🙋
公転周期はケプラーの第三法則ですよね？
🎓
そう、$T^2 \propto a^3$ だ。太陽質量・AU 単位で書くと $T_{\mathrm{yr}} = \sqrt{a^3/M_c}$ とシンプルになる。地球（a=1, M_c=1）なら T=1 年、ハレー彗星（a=17.8）なら T ≈ 75 年。本シミュレーターは中心質量も変えられるから、白色矮星（0.6 M_sun）の周りの惑星や、巨星（10 M_sun）の周りの軌道周期もすぐ確認できる。

よくある質問

離心率 e が小さいとき、ケプラーの方程式 M = E − e·sin(E) は M ≈ E に近く、E_0 = M は良い近似となります。Newton 法は二次収束なので、初期値が真値から ε だけずれていると、1 反復後に O(ε²) のずれに、2 反復後に O(ε⁴) のずれになります。e=0.2 の場合 ε ≈ 0.18 rad 程度ですが、2 反復で 10⁻⁴ rad、3 反復で 10⁻¹⁰ rad と急速に縮まります。一方 e > 0.6 や M が π に近い場合は初期値の質が落ちるため、本ツールでも反復回数が増える様子が観察できます。
これがケプラーの第一法則の核心で、惑星は楕円軌道を描き、太陽はその焦点の一つに位置するという経験則です。Newton 力学から導かれる中心力（万有引力）下の二体問題を解くと、軌道は円錐曲線（楕円・放物線・双曲線）になり、引力源は必ずその焦点に来ます。直感的には、焦点では太陽からの距離 r が点ごとに変わる（楕円の幾何学性質）ため、近日点では速く、遠日点では遅く動く（面積速度一定 = ケプラーの第二法則）。シミュレーターでも、M を 0〜360° スイープすると衛星が近日点付近で速く、遠日点付近で遅く動くのが確認できます。
そのままでは使えません。e ≥ 1 では軌道は閉じず周期が無限大になり、本ツールが扱う「楕円ケプラー方程式」とは別の双曲線版・放物線版（Barker の方程式など）が必要になります。彗星のように放物線・双曲線軌道で太陽系を横切る天体は、$M_h = e\sinh F - F$（双曲線）や $D + D^3/3 = 2M_p$（放物線、Barker）といった専用方程式を使います。本シミュレーターでは安全のため e の上限を 0.95 にしており、楕円軌道の領域を超えない設計です。
はい、GPS 受信機内部でリアルタイムにケプラー方程式を解いています。GPS 衛星は地球を中心とする半径約 26,600 km・離心率約 0.01 のほぼ円軌道を描いており、各衛星は航法メッセージで自分の軌道要素（半長径・離心率・近点引数・昇交点経度・軌道傾斜・近点通過時刻など）を放送します。受信機はそれと現在時刻から M を計算し、Newton 法で E を解いて衛星位置を決定し、複数衛星までの距離から自分の位置を三角測量で求めます。e が小さいので通常 2〜3 反復で十分です。

実世界での応用

衛星軌道計画と運用：あらゆる人工衛星・宇宙ステーション・宇宙望遠鏡の軌道予測でケプラー方程式が中核となります。NASA の SGP4/SDP4 や ESA の Orekit といった軌道伝播ライブラリは、初期軌道要素から将来の位置を計算するためにケプラー方程式を反復的に解きます。火星探査機の運用、月や小惑星への着陸シーケンス、デブリ衝突予測、これらすべてが基本原理として本シミュレーターと同じ計算を行っています。

系外惑星検出と特性化：トランジット法・視線速度法のデータ解析では、観測された光度・速度変化を「ケプラー軌道モデル」にフィッティングして惑星パラメータ（質量・離心率・公転周期など）を推定します。Kepler 宇宙望遠鏡（皮肉にも同名）が発見した 2,000 個以上の系外惑星も、本シミュレーターと同じ方程式を MCMC や最小二乗で繰り返し解くことで特性化されました。

GPS・GNSS 測位：GPS・Galileo・BeiDou・GLONASS すべての GNSS（衛星測位システム）は、各衛星が放送するエフェメリスデータ（軌道要素）から受信機がリアルタイムでケプラー方程式を解いて衛星位置を計算します。スマートフォンの GPS 機能は毎秒これを実行しており、計算量は小さいながら正確さがメートル級測位の前提となります。

小惑星・彗星の軌道決定：小惑星探査ミッション（はやぶさ2、OSIRIS-REx 等）の標的天体軌道は、地上望遠鏡の観測データから軌道要素を最小二乗フィットで推定し、ケプラー方程式で将来位置を予測します。地球接近天体（NEO）監視ネットワークも同じ枠組みで、潜在的に危険な小惑星の軌道を数十年先まで計算し続けています。

よくある誤解と注意点

最も多い誤解は、「M（平均近点角）と E（離心近点角）は同じ角度」と考えてしまうことです。両者は離心率 0 の円軌道でのみ一致し、楕円では一般に異なります。M は時間に比例する仮想的な角度で、円軌道なら衛星位置の角度に対応しますが、楕円ではあくまで「時刻の表現」にすぎません。実際の幾何学的位置は E が、太陽から見た角度は ν が与えます。シミュレーターで M=60° と入力したとき出てくる E=70.82° は、衛星が近日点を出発してから 60°/360° = 1/6 周期分の時間経過後に「中心から見て 70.82°、太陽から見て 82.10°」の位置にいることを意味します。

次に多いのが、「Newton 法は常に高速収束する」と思い込むことです。確かに楕円ケプラー方程式は f'(E) = 1 − e·cos(E) > 0 なので Newton 法が単調収束しますが、e が 1 に近く M が π 付近の組み合わせでは初期値 E_0 = M の質が悪くなり、二次収束に入るまで数反復必要となります。本ツールでは最大反復回数を 30 に設定していますが、本格的な軌道計算では Halley 法（三次収束）や Laguerre 法のような高次法を使うこともあります。スライダーで e=0.95、M=170° を試すと反復回数の増加が観察できます。

最後に、「ケプラー方程式は楕円軌道の本質を表している」と過大評価することです。実際にはケプラー方程式は二体問題（太陽 + 惑星 1 個）の純粋なケースに限定されます。実際の太陽系では他の惑星の重力摂動が小さいながら存在し、相対論効果（水星の近日点移動が典型例）も加わります。GPS や深宇宙探査では、ケプラー軌道を「ベース」として、摂動項を逐次加える Cowell 法・Encke 法といった軌道計算法に発展させます。本シミュレーターはあくまで「純粋ケプラー軌道」の出発点として位置づけてください。

関連ツール

ケプラーの法則シミュレーター ケプラーの第1・2・3法則を軌道アニメーションと面積速度で視覚的に確認するシミュレーター。離心率・軌道周期・惑星間距離の関係をインタラクティブに探索できる天文学ツール。 ケプラー軌道シミュレーター ケプラー軌道シミュレーターで惑星の運動を可視化。離心率や長半径を操作し、楕円軌道や速度変化をリアルタイムアニメーションで体感。万有引力によるハレー彗星や地球・火星の軌道をシミュレーシ... シェラーの式 シミュレーター — X 線回折からの結晶子サイズ推定 シェラーの式 D = Kλ/(β cosθ) のリアルタイムシミュレーター。X 線波長 λ、ブラッグ角 θ、実測 FWHM β_total、装置 FWHM β_inst から結晶子サ... 軌道力学シミュレーター ケプラーの法則に基づく楕円軌道をリアルタイムアニメーション。離心率・軌道長半径・周期の相互関係をインタラクティブに探索できる天文学ツール。 軌道力学シミュレーター 軌道力学シミュレーターで、半長軸・離心率・傾斜角を操作し、ケプラー楕円軌道をリアルタイムに可視化。ケプラー方程式をニュートン・ラフソン法で数値的に解き、衛星の正確な位置と動きをアニメ... 軌道・投射物シミュレーター 軌道・投射物シミュレーター「NovaSolver」は、初速度の調整だけで弾道飛行から円軌道、楕円軌道、さらには双曲線脱出軌道までをリアルタイムに可視化。ケプラーの法則や軌道周期、離心...