開普勒方程

開普勒方程（）把軌道力學中受到連心力影響的軌道中天體多種幾何性質聯繫起來。

這條方程最早由約翰內斯·開普勒得出，推導過程見於他在1609年出版的著作《新天文學》的第60章，而他在著作《》（1621年出版）的第五卷中還提出了此方程的一個迭代解。這條方程在物理學史和數學史上都扮演着重要的角色，特別是在古典天體力學史上。

方程

開普勒方程如下：

其中為平近點角，為偏近點角，而則為離心率.

偏近點角對於計算開普勒軌道上移動點的位置相當有用。比方說，天體在座標、和時間時經過近星點，那麼要找出天體在任意時間的位置的話，首先可由時間和平均運動用方程計算出平近點角，然後解上述的開普勒方程得，便能由下列方程求座標：

其中為半長軸，而則為半短軸。

由於正弦是超越函數，所以開普勒方程是超越方程，即是說不能用代數方法解出。一般求需要用到數值分析和級數。

 其他形式 

開普勒方程共有幾種形式。每一種形式都同一種特定軌道類型有關。標準的開普勒方程用於橢圓軌道（0 ≤<1）。雙曲開普勒方程用於雙曲軌跡（≫1）。徑向開普勒方程用於線性（徑向）軌跡 (=1)。抛物線軌跡（=1）則用。

當e = 0時，軌道是圓形的。增加e會導致圓形變成橢圓。當e = 1時共有三種可能性：
 抛物線軌跡；
 沿着從吸引中心起始無限長射線向內或向外的軌跡；
 或沿着由吸引中心和距其一定距離的點所形成的線段來回移動的軌跡。

把從1輕微增加會形成轉角剛好低於180度的雙曲軌道。繼續增加數值會導致轉角變小，而當e趨向無限時，軌道則變成長度無限的直線。

 雙曲開普勒方程 

雙曲開普勒方程如下：

其中為雙曲偏近點角。
這條方程是通過把重新定義為−1的平方根乘上橢圓方程的右邊而得：

（其中現為虛數），然後以取代。

徑向開普勒方程

徑向開普勒方程如下：

其中與時間成正比，而則與射線上與吸引中心的距離成正比。下式是通過把開普勒方程成並把定為1而成：

。

代入得

。

 逆問題 

可以利用已知的值直接求出。但用已知值來求卻要複雜得多。因為不存在解析解。

可以使用寫出開普勒方程解的級數表達式，但所有和的組合都不能使級數收斂（見下文）。

文獻中對開普勒方程可解性的混淆已存在了四個世紀開普勒方程常被聲稱“沒有解釋解”；例子見這裏 。至於這是否真確取決於是否認為無限級數（或不一定收歛的級數）算解釋解。也有其他作者無理地聲稱此方程無解；例子見M. V. K. Chari, Sheppard Joel Salon 2000 Technology & Engineering.開普勒本人曾對找得到通解的可能性表示懷疑。

 逆開普勒方程 

逆開普勒方程是所有實數值的開普勒方程解：



展開得：



以上的級數可用Wolfram Mathematica的InverseSeries運算得出。
 InverseSeries[Series[M - Sin[M], {M, 0, 10}]]
 InverseSeries[Series[M - e Sin[M], {M, 0, 10}]]

這些函數只是簡單的麥克勞林級數。這樣的超越函數泰勒級數寫法可被視為那些函數的定義。因此，這個解是逆開普勒方程的正式定義。然而，在非零的時候，並不是的整函數。其導數

在<1時於無限複數集合中趨向零。在有解，而在這些值時

（其中逆cosh取正值），而在這些點則趨向無限。這意味着此麥克勞林級數的收歛半徑為，並且當比這大時級數不會收歛。此級數還可用於雙曲情況，此時其收歛半徑為。當=1時此級數只在時收歛。

還可以寫出用冪表示的麥克勞林級數。這級數在大於拉普拉斯極限（約為0.66）時不會收歛，與值無關（除非為的倍數），但若小於拉普拉斯極限時則級數在任何值時都會收歛。級數的係數除了第一個之外（只是而已），都與成週期式關係，週期為。

 逆徑向開普勒方程 

逆徑向開普勒方程（e = 1）可被寫成：


展開得：

上述結果可用Wolfram Mathematica求得:
InverseSeries[Series[ArcSin[Sqrt[t]] - Sqrt[(1 - t) t], {t, 0, 15}]]

 逆問題的數值近似 

逆問題在大部份的應用中都能以數值方法求得函數的根：



可以經牛頓法來進行迭代：



注意在這個計算和的單位是弧度角。重覆迭代直到已經達到想要的準確度（例如，當 < 所需的準確度）。對大部份的橢圓軌道而言，起始值取E0 = M(t) 已經足夠。而對e > 0.8的軌道而言，則應取起始值。相近的手法還可以用於開普勒方程的雙曲形式。而在面對抛物線軌跡的個案時則使用。

 定點迭代 

另一種有關的解法由考慮開始。重覆地將的值代入右方式子就能得到簡單的求的算法。此方法與開普勒1621年的解相同。

function E(e,M,n)
    E = M
    for k = 1 to n
        E = M + e*sin E
    next k
    return E

迭代次數取決於的數值。雙曲形式則用相近的。

這個方法與上文牛頓法解的關係如下：



若及的數值小的話，取一次項得近似：

。

 另見 
開普勒定律
開普勒問題
廣義相對論中的開普勒問題

 參考資料 

 外部連結 










 Wolfram Mathworld的開普勒方程頁面 

Category:軌道
Category:约翰内斯·开普勒