Keplerin yhtälö

Keplerin yhtälö on eräs taivaanmekaniikan perustavimmista yhtälöistä. Se ilmaisee yhteyden kahden taivaankappaleiden (planeetta, asteroidi, komeetta, tms.) liikkeiden laskemisessa tarvittavan suureen, keskianomalian ${\displaystyle M}$ ja eksentrisen anomalian ${\displaystyle E}$, välillä. Keplerin yhtälö on näennäisesti yksinkertainen:

${\displaystyle M=E-e\sin E}$,

missä ${\displaystyle e}$ on taivaankappaleen radan eksentrisyys. Yhtälössä sekä ${\displaystyle E}$ että ${\displaystyle M}$ on lausuttu radiaaneina. Varsinainen ongelma on, että tyypillisesti ${\displaystyle M}$ ja ${\displaystyle e}$ tunnetaan, ja eksentrinen anomalia ${\displaystyle E}$ täytyisi saada selville. Tämän ratkaiseminen ei ole mitenkään mahdollista analyyttisesti vaan ratkaisussa joudutaan aina turvautumaan numeerisiin menetelmiin.

Yllä oleva muoto Keplerin yhtälöstä on voimassa ellipsiradoille. Hyperbeliradoille se pätee myös, mutta tällöin kulma ${\displaystyle E}$ on imaginäärinen. Paraabeliradalle yhtälöstä on olemassa oma erikoistapauksensa, jota kutsutaan Barkerin yhtälöksi.

Jos taivaankappaleen rata halutaan esittää ajan funktiona, lasku vaatii käytännössä eksentrisen anomalian ratkaisemista jokaista aika-askelta varten. Tästä syystä Keplerin yhtälön mahdollisimman tehokas ratkaisu on tärkeää. Yhtälön muoto tuottaa kuitenkin hankaluutta, sillä jos rata on hyvin soikea eli ${\displaystyle e}$ on lähellä ykköstä ja etenkin jos tämän lisäksi ${\displaystyle M}$ on pieni (lähellä nollaa), yhtälön iterointi suppenee varsin hitaasti ja voi olla numeerisesti epätarkka. Tästä syystä erilaisia menetelmiä yhtälön ratkaisemiseksi on esitetty satoja.

Eräs kätevimpiä menetelmiä on seuraava. Merkitään ${\displaystyle f=E-e\sin E-M}$, ja lasketaan tämän ensimmäinen ja toinen derivaatta ${\displaystyle E}$:n suhteen eli ${\displaystyle f\,'=1-e\cos E}$ ja ${\displaystyle f\,''=e\sin E}$. Nyt yhtälö voidaan ratkaista valitsemalla alkuarvoksi ${\displaystyle E_{0}=M}$ ja iteroimalla kaavasta

${\displaystyle E_{n+1}=E_{n}-{\frac {f_{n}}{\sqrt {(f'_{n})^{2}-f_{n}f''_{n}}}}}$,

kunnes haluttu tarkkuus on saavutettu (eli kunnes ${\displaystyle |E_{n+1}-E_{n}|<\epsilon }$, missä ${\displaystyle \epsilon }$ on yleensä luokkaa 10⁻⁶ tai pienempi). Tämä menetelmä toimii aina ja suppenee hyvin nopeasti.