Kepler-Gleichung

Die Kepler-Gleichung i​st eine transzendente Gleichung z​ur Berechnung elliptischer Bahnen v​on Himmelskörpern. Sie ergibt s​ich aus d​en ersten beiden keplerschen Gesetzen, d​ie Johannes Kepler 1609 publizierte, u​nd lautet

${\displaystyle M=E-e\cdot \sin E}$.

Sie verknüpft d​en Ort e​ines Himmelsobjekts m​it der Zeit.

Die Kepler-Gleichung k​ann z. B. b​ei der Ermittlung d​er Zeitgleichung angewendet werden. Eine Teilaufgabe d​abei ist, d​ie wahre Anomalie d​er Erde a​uf ihrer elliptischen Bahn u​m die Sonne anzugeben.

Herleitung

Zur Kepler-Gleichung auf elliptischer Keplerbahn

Längen:	Punkte:
${\displaystyle a\!:}$ große Halbachse	${\displaystyle \mathrm {C} \!:}$ Mittelpunkt
${\displaystyle b\!:}$ kleine Halbachse	${\displaystyle \mathrm {S} \!:}$ Brennpunkt
${\displaystyle e\!\cdot \!a\!:}$ lineare Exzentrizität	${\displaystyle \mathrm {Z} \!:}$ Periapsis
Winkel:
${\displaystyle T\!:}$ wahre Anomalie	${\displaystyle \mathrm {P} \!:}$ Objekt
${\displaystyle E\!:}$ exzentrische       Anomalie	${\displaystyle \mathrm {X} \!:}$ Hilfspunkt zum       Objekt
${\displaystyle M\!:}$ mittlere Anomalie	${\displaystyle \mathrm {Y} \!:}$ fiktives Objekt

Das zweite keplersche Gesetz, der Flächensatz, folgt aus der Drehimpulserhaltung im Zweikörperproblem, welch letzteres in der Astronomie auch Kepler-Problem genannt wird. Nur eine radiale Kraft wirkt hier im Fahrstrahl vom Schwerezentrum ${\displaystyle \mathrm {S} }$ zum Himmelsobjekt ${\displaystyle \mathrm {P} }$. Gehorcht diese Kraft überdies einem ${\displaystyle 1/r^{2}}$-Gesetz (wie die Newtonsche Gravitationskraft), ist also der Gesamt-Kraftfluss durch alle Kugeloberflächen gleich (d. h. unabhängig vom Kugelradius ${\displaystyle r}$), dann ist die Planetenbahn ein Kegelschnitt, wie es das erste keplersche Gesetz besagt. Im Fall der periodischen Wiederkehr des Himmelsobjekts, dem Fall der Ellipse, ist die Kepler-Gleichung

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

die in eine Rechenformel gegossene Aussage des Flächensatzes. Sie bringt die Zeit ${\displaystyle t}$ in Form der (von Kepler so genannten) mittleren Anomalie ${\displaystyle M}$ mit der Position des Himmelsobjekts ${\displaystyle \mathrm {P} }$ auf seiner Umlaufbahn (Kepler-Ellipse „Orbit“) in Form der (von Kepler so genannten) wahren Anomalie ${\displaystyle T}$, d. i. sein Winkelabstand von der Periapsis ${\displaystyle \mathrm {Z} }$, (über die Hilfsgröße der exzentrischen Anomalie ${\displaystyle E}$) in einen eindeutigen formelmäßigen Zusammenhang.

Dabei ist ${\displaystyle e}$ die numerische Exzentrizität der Ellipse.

Mittlere Anomalie

Die gleichmäßig vergehende Zeit lässt sich mit der Bewegung eines fiktiven Körpers (${\displaystyle \mathrm {Y} }$ in den Abbildungen) auf einer Kreisbahn mit konstanter Winkelgeschwindigkeit gleichsetzen. Hierfür wird ein „Umkreis“ als Hilfskreis, auf dem ${\displaystyle \mathrm {Y} }$ umläuft, um die Kepler-Ellipse gelegt. Zum Zeitpunkt ${\displaystyle t_{P}}$ seien sowohl ${\displaystyle \mathrm {Y} }$ wie das wahre Objekt ${\displaystyle \mathrm {P} }$ als in der Periapsis ${\displaystyle \mathrm {Z} }$ stehend angenommen. Beide Punkte haben dieselbe Umlaufzeit und stehen bei jedem ganzzahligen Vielfachen in der Periapsis und bei jedem halbzahligen in der Apoapsis (zusammen).

zu Gleichung ${\displaystyle \mathrm {(2)} }$

Die momentane Lage des Punktes ${\displaystyle \mathrm {Y} }$ wird als Winkel (alle folgenden Winkel werden mit Bogenmaß dargestellt) im Hilfskreis- (und Ellipsen-)Mittelpunkt ${\displaystyle \mathrm {C} }$ im Bezug zur Periapsis ${\displaystyle \mathrm {Z} }$ angegeben und als mittlere Anomalie ${\displaystyle M}$ bezeichnet:

${\displaystyle \mathrm {(1)} \quad M=2\pi {\frac {t-t_{P}}{U}}}$ .

Dabei ist ${\displaystyle U}$ die Bahnperiode und ${\displaystyle {2\pi }/U}$ die mittlere Winkelgeschwindigkeit. Im Zeitpunkt ${\displaystyle t_{P}}$ befindet sich das Himmelsobjekt in der Periapsis, wo es den geringsten Abstand zu seinem Schwerezentrum ${\displaystyle \mathrm {S} }$ hat.

Gemäß dem zweiten keplerschen Gesetz überstreicht der Fahrstrahl ${\displaystyle {\overline {\mathrm {SP} }}}$ des Körpers ${\displaystyle \mathrm {P} }$ im gleichen Zeitabschnitt die gleiche Fläche. Da der Zeitanteil (am Umlauf) proportional ist zum Anteil des Kreissektors am Umkreis, ist der Anteil der elliptischen Teilfläche ${\displaystyle \mathrm {SPZ} }$ an der Ellipse gleich groß wie der des Kreissektors ${\displaystyle \mathrm {CYZ} }$ am Umkreis:

${\displaystyle \mathrm {(2)} \quad {\frac {\operatorname {area} \,\mathrm {CYZ} }{\operatorname {area} \,\mathrm {SPZ} }}={\frac {\pi a^{2}}{\pi ab}}={\frac {a}{b}}}$ .

Dabei ist ${\displaystyle a}$ die große Halbachse der Ellipse und gleichzeitig der Radius des Umkreises, ${\displaystyle b}$ die kleine Halbachse der Ellipse. Ellipse und Umkreis sind im Verhältnis ${\displaystyle b/a}$ affin zueinander, d. h., die Ellipse ist in jeder Parallele zur kleinen Halbachse der mit diesem Verhältnis „gestauchte“ Umkreis.

zu Gleichung ${\displaystyle \mathrm {(4)} }$

zu Gleichung ${\displaystyle \mathrm {(3)} }$

Exzentrische Anomalie

Durch eine zur kleinen Halbachse parallele Projektion des Punktes ${\displaystyle \mathrm {P} }$ auf den Umkreis entsteht der Hilfspunkt ${\displaystyle \mathrm {X} }$, dessen Winkel im Mittelpunkt ${\displaystyle \mathrm {C} }$ zur Periapsis ${\displaystyle \mathrm {Z} }$ von Kepler exzentrische Anomalie ${\displaystyle E}$ genannt wurde. Die Affinität begründet folgenden Zusammenhang:

${\displaystyle \mathrm {(3)} \quad \operatorname {area} \,\mathrm {SXZ} ={\frac {a}{b}}\operatorname {area} \,\mathrm {SPZ} }$ .

Nach Einsetzen von Gleichung ${\displaystyle \mathrm {(2)} }$ in Gleichung ${\displaystyle \mathrm {(3)} }$ folgt:

${\displaystyle \mathrm {(4)} \quad \operatorname {area} \,\mathrm {SXZ} =\operatorname {area} \,\mathrm {CYZ} }$ .

Keplergleichung

zu den Gleichungen ${\displaystyle \mathrm {(5)} }$ und ${\displaystyle \mathrm {(6)} }$

zu Gleichung ${\displaystyle \mathrm {(7)} }$

Mit der Gleichung ${\displaystyle \mathrm {(4)} }$ ist die gesuchte, das zweite keplersche Gesetz erfüllende Beziehung zwischen der exzentrischen Anomalie (Punkt ${\displaystyle \mathrm {X} }$) und der mittleren Anomalie (Punkt ${\displaystyle \mathrm {Y} }$) implizit gefunden. Eine explizite Beziehung ergibt sich durch folgende Schritte:

Wenn der Fahrstrahl ${\displaystyle {\overline {\mathrm {CY} }}}$ in einer Periode ${\displaystyle U}$ den Winkel ${\displaystyle 2\pi }$ zurücklegt und die Fläche ${\displaystyle \pi a^{2}}$ überstreicht, so überstreicht er bis zum Zeitpunkt ${\displaystyle t}$ den Winkel ${\displaystyle M}$ und eine um den Faktor ${\displaystyle M/2\pi }$ kleinere Fläche:

${\displaystyle \mathrm {(5)} \quad \displaystyle \operatorname {area} \,\mathrm {CYZ} ={\frac {a^{2}}{2}}M}$ .

Die analoge Betrachtung für den Fahrstrahl ${\displaystyle {\overline {\mathrm {CX} }}}$ über den Winkel ${\displaystyle E}$ ergibt:

${\displaystyle \mathrm {(6)} \quad \displaystyle \operatorname {area} \,\mathrm {CXZ} ={\frac {a^{2}}{2}}E}$ .

Die Fläche ${\displaystyle \mathrm {CXZ} }$ besteht aus den Teilflächen ${\displaystyle \mathrm {CXS} }$ und ${\displaystyle \mathrm {SXZ} }$:

${\displaystyle \mathrm {(7)} \quad \displaystyle \operatorname {area} \,\mathrm {CXZ} =\operatorname {area} \,\mathrm {CXS} +\operatorname {area} \,\mathrm {SXZ} }$ .

Die Teilfläche ${\displaystyle \mathrm {CXS} }$ (hellblau umrandet in der Abbildung) ist ein geradlinig begrenztes Dreieck mit der Basis ${\displaystyle e\cdot a}$ und der Höhe ${\displaystyle a\cdot \sin E}$ :

${\displaystyle \mathrm {(8)} \quad \displaystyle \operatorname {area} \,\mathrm {CXS} ={\frac {ea\cdot a\sin E}{2}}={\frac {a^{2}}{2}}\,e\sin E}$ .

${\displaystyle e}$ ist die numerische Exzentrizität der Ellipse und ${\displaystyle ea={\sqrt {a^{2}-b^{2}}}}$ die lineare, die den Abstand zwischen Mittelpunkt und Brennpunkt angibt.

Die Teilfläche ${\displaystyle \mathrm {SXZ} }$ ist nach Gleichung ${\displaystyle \mathrm {(4)} }$ gleich groß wie die Fläche ${\displaystyle \mathrm {CYZ} }$, deren Wert in Gleichung ${\displaystyle \mathrm {(5)} }$ angegeben ist.

Durch Einsetzen der Gleichungen ${\displaystyle \mathrm {(6)} }$, ${\displaystyle \mathrm {(8)} }$ und ${\displaystyle \mathrm {(5)} }$ wird aus Gleichung ${\displaystyle \mathrm {(7)} }$ die Gleichung

${\displaystyle \mathrm {(9)} \quad \displaystyle {\frac {a^{2}}{2}}E={\frac {a^{2}}{2}}e\sin E+{\frac {a^{2}}{2}}M}$ .

Daraus ergibt s​ich schließlich d​ie Kepler-Gleichung:

${\displaystyle E-e\cdot \sin E=M}$ .

Lösung der Kepler-Gleichung

Die Kepler-Gleichung ist nicht in geschlossener Form nach der exzentrischen Anomalie ${\displaystyle E(t)}$ auflösbar. Beispiele dafür, wie ${\displaystyle E(t)}$ mit ihr aus der mittleren Anomalie ${\displaystyle M(t)=2\pi {\frac {t-t_{P}}{U}}}$ ermittelt werden kann:

1. ${\displaystyle E-M}$ ist eine ungerade, mit ${\displaystyle 2\pi }$ periodische Funktion in ${\displaystyle M}$. Als solche lässt sie sich in eine Fourierreihe entwickeln, die für alle ${\displaystyle M\in \mathbb {R} }$ und ${\displaystyle e\in \mathbb {R} }$ konvergiert, und zwar ist

   ${\displaystyle F(M):=E-M=e\cdot \sin E=2\cdot \sum _{n=1}^{\infty }{\frac {J_{n}(ne)}{n}}\sin(nM)}$

   mit ${\displaystyle J_{n}}$ als Bessel-Funktion erster Gattung ${\displaystyle n}$-ter Ordnung.[1][2]
   

   Aus den Werten ${\displaystyle F(M^{\prime })}$ für ${\displaystyle M^{\prime }\in [0,\pi ]}$ lassen sich alle anderen Werte ${\displaystyle F(M)}$ leicht berechnen:

   ${\displaystyle F(M)=s\cdot F{\bigl (}s\cdot M^{\prime }{\bigr )}}$

   mit ${\displaystyle k:={\bigl \lfloor }{\tfrac {M}{\pi }}{\bigr \rfloor }\in \mathbb {Z} }$ (Gaußklammer), ${\displaystyle s:=(-1)^{k}\in \{1,-1\}}$ und ${\displaystyle M^{\prime }:=M-{\bigl (}k+{\tfrac {1-s}{2}}{\bigr )}\pi \in [-\pi ,+\pi ]}$, sodass ${\displaystyle s\cdot M^{\prime }\in [0,+\pi ]}$.

2. Eine Nullstelle ${\displaystyle E}$ der Funktion

   ${\displaystyle f(E)=E-e\cdot \sin E-M}$

   ist eine Lösung der Keplergleichung. Die Nullstelle kann etwa mit dem Newton-Verfahren wie folgt numerisch berechnet werden:

   ${\displaystyle E_{n+1}=E_{n}-{\frac {f(E_{n})}{f'(E_{n})}}=E_{n}-{\frac {E_{n}-e\sin(E_{n})-M}{1-e\cos(E_{n})}}}$ .

   Für die meisten elliptischen Bahnen ist der Anfangswert ${\displaystyle E_{0}=M}$ geeignet. Für Exzentrizitäten ${\displaystyle 0{,}8<e<1}$ kann ${\displaystyle E_{0}=\pi }$ genommen werden.

3. Ein stabileres, aber langsamer konvergierendes Verfahren beruht auf dem banachschen Fixpunktsatz:[3]

   ${\displaystyle E_{0}=M,\qquad E_{n+1}=M+e\cdot \sin E_{n}}$ .

4. Für kleine Exzentrizität ${\displaystyle e}$ kann ${\displaystyle E}$ auch folgendermaßen approximiert werden:[4]

   ${\displaystyle E=M+e\cdot \sin M+{\frac {1}{2}}e^{2}\cdot \sin 2M}$

   Der Fehler ist hierbei von der Größenordnung ${\displaystyle {\mathcal {O}}(e^{3})}$. Bei der Erde und ihrer Exzentrizität ${\displaystyle e=0{,}0167}$ liegt der Fehler für begrenzte Zeiträume hinter der 5. Kommastelle.

5. Eine Auflösung für ${\displaystyle e<1}$ nach Art der Lagrangeschen Inversionsformel ist die Maclaurin-Reihe in ${\displaystyle M}$   ${\displaystyle {\begin{array}{ll}\textstyle E={\frac {1}{1-e}}M\!\!\!\!&-{\frac {e}{(1-e)^{4}}}{\frac {M^{3}}{3!}}+{\frac {(e+9e^{2})}{(1-e)^{7}}}{\frac {M^{5}}{5!}}-{\frac {(e+54e^{2}+225e^{3})}{(1-e)^{10}}}{\frac {M^{7}}{7!}}\\&+{\frac {(e+243e^{2}+4131e^{3}+11025e^{4})}{(1-e)^{13}}}{\frac {M^{9}}{9!}}-{\frac {(e+1008e^{2}+50166e^{3}+457200e^{4}+893025e^{5})}{(1-e)^{16}}}{\frac {M^{11}}{{11}!}}\\&+{\frac {(e+4077e^{2}+520218e^{3}[phone]e^{4}[phone]e^{5}[phone]e^{6})}{(1-e)^{19}}}{\frac {M^{13}}{{13}!}}\mp \cdots ,\end{array}}}$
   die für ${\displaystyle |M|<\cosh ^{-1}(e^{-1})-{\sqrt {1-e^{2}}}}$ linear konvergiert. Ist also ${\displaystyle 0\leq e\lessapprox 0{,}031803066}$, dann konvergiert sie für ${\displaystyle |M|\leq \pi }$ linear.
   Die Koeffizienten der Zähler-Polynome in ${\displaystyle e}$ sind in der Folge A306557 in OEIS festgehalten.

Lösung einiger Teilaufgaben im Kepler-Problem

Wahre Anomalie

Für einen Himmelskörper auf einer Keplerbahn ist für den Zeitpunkt ${\displaystyle t}$ beziehungsweise für die zugehörige mittlere Anomalie ${\displaystyle M(t)}$ der Ort beziehungsweise die wahre Anomalie ${\displaystyle T(t)}$ anzugeben. Mit Hilfe der Kepler-Gleichung wird zuerst die exzentrische Anomalie ${\displaystyle E(t)}$ ermittelt (siehe oben). Aus Letzterer folgt die wahre Anomalie ${\displaystyle T(t)}$ nach einer der folgenden Beziehungen:[5]

${\displaystyle \tan {\frac {T}{2}}={\sqrt {\frac {1+e}{1-e}}}\cdot \tan {\frac {E}{2}}}$

oder

${\displaystyle \cos T={\frac {a\cos E-ae}{a-ae\,\cos E}}={\frac {\cos E-e}{1-e\cos E}}}$

Hier ist ${\displaystyle ae={\sqrt {a^{2}-b^{2}}}}$ die lineare Exzentrizität der Bahnellipse. Zum Auflösen nach ${\displaystyle T}$ ist jeweils eine Unterscheidung der Fälle ${\displaystyle 0\leq E\leq \pi }$ und ${\displaystyle \pi \leq E\leq 2\pi }$ nötig.

Bemerkungen

• Der Nenner der zweiten Formel gibt gerade den Abstand ${\displaystyle r}$ des Himmelsobjekts zum Brennpunkt ${\displaystyle s}$ an:

${\displaystyle r=a-ae\,\cos E}$

• Die Formeln können leicht nach ${\displaystyle \tan {\tfrac {E}{2}}}$ oder ${\displaystyle \cos E}$ aufgelöst werden, es ergibt sich:[6]

${\displaystyle \tan {\frac {E}{2}}={\sqrt {\frac {1-e}{1+e}}}\cdot \tan {\frac {T}{2}}}$

und

${\displaystyle \cos E={\frac {a\cos T+ae}{a+ae\cos T}}={\frac {\cos T+e}{1+e\cos T}}}$

Zwischen der wahren Anomalie ${\displaystyle T,}$ der exzentrischen Anomalie ${\displaystyle E}$ und der mittleren Anomalie ${\displaystyle M}$ bestehen noch zahlreiche weitere Zusammenhänge,[7] die in der langen Geschichte der Himmelsmechanik entwickelt wurden. Insbesondere lässt sich die wahre Anomalie – ohne Umweg über die Keplergleichung – direkt aus einer speziellen Differenzialgleichung in ${\displaystyle M}$ errechnen,[8] was für numerische Näherungsverfahren von Interesse ist.

Insbesondere kann auch hier die wahre Anomalie ${\displaystyle T}$ durch die mittlere Anomalie ${\displaystyle M}$ für kleine Exzentrizitäten genähert werden, es ergibt sich die nützliche Näherung

${\displaystyle T=M+2e\sin(M)+{\frac {5}{4}}e^{2}\sin(2M)+{\mathcal {O}}(e^{3}).}$

Die Differenz ${\displaystyle T}$ − ${\displaystyle M}$ heißt Mittelpunktsgleichung.[8]

Mittlere Anomalie

Für einen Himmelskörper auf einer Keplerbahn mit der wahren Anomalie ${\displaystyle T}$ ist die zugehörende mittlere Anomalie ${\displaystyle M(T)}$ beziehungsweise der zugehörende Zeitpunkt ${\displaystyle t(T)}$ anzugeben. Es handelt sich um die zur obigen umgekehrte Aufgabenstellung.

Ausgehend von ${\displaystyle T}$ ergibt sich die exzentrische Anomalie zu

${\displaystyle E=2\arctan _{\tfrac {T}{2}}\left({\sqrt {\frac {1-e}{1+e}}}\cdot \tan {\frac {T}{2}}\right)}$ .

Der Lageparameter-Index ${\displaystyle {\tfrac {T}{2}}}$ bei ${\displaystyle \arctan }$ gibt denjenigen Wert des Arkustangens zurück, der diesem (${\displaystyle {\tfrac {T}{2}}}$) am nächsten liegt (siehe Arkustangens mit Lageparameter). Die Kepler-Gleichung liefert die ${\displaystyle E}$ zugehörige mittlere Anomalie

${\displaystyle M(t)=E(t)-{\frac {180^{\circ }}{\pi }}\cdot e\cdot \sin E(t)}$ .

Aus d​er linearen Gleichung für d​as Bahnelement f​olgt schließlich:

${\displaystyle t={\frac {M-M_{0}}{\dot {M}}}}$

Beispiel

Passagezeiten der vier Erdbahnellipsen-Scheitel:
Die für die Erde gültigen Bahnelemente sind unter mittlere Kepler-Elemente angegeben. Die im Verweisartikel verwendete Zeit ${\displaystyle T}$ ist in Julianischen Jahrhunderten gerechnet. Hier wird ${\displaystyle t}$ in Tagen gemessen, sodass die linearen Koeffizienten der Zeit ${\displaystyle T}$ durch 36525 zu teilen sind, um ${\displaystyle {\dot {M}}}$ und ${\displaystyle {\dot {e}}}$ zu erhalten. Die sehr langsame Änderung der numerischen Exzentrizität wird allerdings vernachlässigt (${\displaystyle {\dot {e}}=0}$). Der Nullpunkt der Zeit ${\displaystyle T}$ – und damit auch von ${\displaystyle t}$ – ist der 1. Januar 2000, 12:00 UT. Die wahre Anomalie bei Perihelpassage der Erde im Jahr 2000 ist gleich 360° (nicht null!), im Jahr 2001 gleich 720° usw.

Passagezeiten der vier Erdbahnellipsen-Scheitel

	Perihel 2000	Frühlings- Nebenscheitel	Aphel	Herbst- Nebenscheitel	Perihel 2001
Wahre Anomalie ${\displaystyle T/^{\circ }}$	360	450	540	630	720
Zeit ${\displaystyle t/{\text{d}}}$	2,511	91,883	185,140	278,398	367,770
Zeitabstand ${\displaystyle \Delta t/{\text{d}}}$		89,372	93,258	93,258	89,372

Der Abstand zwischen den mittleren Perihelpassagen (anomalistisches Jahr) beträgt ${\displaystyle J_{an}=360^{\circ }/{\dot {M}}=365,260{\text{ d}}.}$ Die so berechneten mittleren Perihelzeiten können sich um mehrere Tage vom realen (vor allem mondgestörten) Wert unterscheiden.

Bahnradius

→ Hauptartikel: Keplerbahn

Mit der wahren Anomalie wird die Richtung eines Himmelskörpers auf seiner Keplerbahn für eine Zeit ${\displaystyle t}$ angegeben. Die zugehörende Entfernung – der Bahnradius – ist wie folgt berechenbar:

${\displaystyle r=r(T(t))=r(t)=a\cdot {\frac {1-e^{2}}{1+e\cdot \cos T}}}$

${\displaystyle r\!:}$ Entfernung (Bahnradius)

${\displaystyle a\!:}$ große Halbachse der Ellipse

${\displaystyle e\!:}$ numerische Exzentrizität

${\displaystyle T\!:}$ wahre Anomalie

Bahngeschwindigkeit

Die zeitliche Änderung der wahren Anomalie entspricht der Winkelgeschwindigkeit ${\displaystyle \omega }$ in Bezug auf das Gravizentrum. Die Normalkomponente der Geschwindigkeit folgt also direkt aus

${\displaystyle v_{\perp }={\dot {T}}\cdot r.}$

Die Radialgeschwindigkeit i​st die Änderung d​es Bahnradius m​it der Zeit:

${\displaystyle v_{r}={\dot {r}}}$

Für die Bahngeschwindigkeit oder Orbitalgeschwindigkeit ${\displaystyle v}$ folgt dann ${\displaystyle v^{2}=v_{\perp }^{2}+v_{r}^{2}.}$

${\displaystyle v=v(T(t),r(t))=v(t)={\sqrt {({\dot {T}}\cdot r)^{2}+{\dot {r}}^{2}}}}$

${\displaystyle v\!:}$ Bahngeschwindigkeit

${\displaystyle T\!:}$ wahre Anomalie

${\displaystyle r\!:}$ Bahnradius

Einfacher lässt sich die Bahngeschwindigkeit über den Hodograph ${\displaystyle {\vec {\dot {r}}}}$ aus dem Flächensatz ableiten:[9]

${\displaystyle v^{2}={\frac {C^{2}}{p}}\left({\frac {2}{r}}-{\frac {1}{a}}\right)}$

${\displaystyle C\!:}$ spezifischer Drehimpuls als zentrale Kenngröße der Bewegung

${\displaystyle C=v_{\mathrm {max} }\cdot r_{\mathrm {min} }=v_{\mathrm {min} }\cdot r_{\mathrm {max} }}$

${\displaystyle p\!:}$ Halbparameter als kennzeichnendes Bahnelement

${\displaystyle p=2\cdot {\frac {r_{\mathrm {min} }\cdot r_{\mathrm {max} }}{r_{\mathrm {min} }+r_{\mathrm {max} }}}={\frac {b^{2}}{a}}}$

${\displaystyle a\!:}$ große Halbachse

${\displaystyle b\!:}$ kleine Halbachse

${\displaystyle {\frac {C^{2}}{p}}=G\cdot M}$ mit Gravitationskonstante ${\displaystyle G}$ und Masse ${\displaystyle M}$ des Zentralkörpers

Daraus folgen d​ie Minimal- u​nd Maximalgeschwindigkeit i​m Apozentrum u​nd Perizentrum e​iner Ellipsenbahn:[9]

${\displaystyle v_{\mathrm {max} }^{2}={\frac {C^{2}}{p\cdot a}}\cdot {\frac {1+e}{1-e}}\qquad v_{\mathrm {min} }^{2}={\frac {C^{2}}{p\cdot a}}\cdot {\frac {1-e}{1+e}}}$

${\displaystyle e\!:}$ numerische Exzentrizität

Anwendung der Kepler-Gleichung bei der Zeitgleichung

Die quantitative, a​lso rechnerische Behandlung[10] d​er Zeitgleichung i​st im Wesentlichen – nämlich b​eim aus d​er elliptischen Bahnbewegung d​er Erde resultierenden Zeitgleichungsanteil – e​ine Anwendung d​er Kepler-Gleichung. Insbesondere w​ird damit d​er Ort d​er Erde a​uf ihrer elliptischen Bahn (auch Keplerbahn) z​u einem vorgegebenen Zeitpunkt bestimmt.

Definitionen der Zeitgleichung

Erste Definition:

${\displaystyle \mathrm {(10)} \quad \mathrm {ZG} =\mathrm {WOZ} -\mathrm {MOZ} }$

Dem Wert der wahren Ortszeit (WOZ) bzw. mittleren Ortszeit (MOZ) entspricht der jeweilige Stand der wahren bzw. einer fiktiven mittleren Sonne am Himmel. Da die Tageszeit im Zusammenhang mit der Drehung der Erde um ihre Achse steht, interessiert nur die jeweilige Rektaszension (nicht die Deklination) der Sonne(n). Anders gesagt: Von den in zwei zueinander rechtwinklig erfolgten scheinbaren Jahresbewegungen der wahren Sonne interessiert nur die auf dem Himmelsäquator stattfindende, jedoch nicht das periodische An- und Absteigen. Die die gleichmäßig vergehende Zeit repräsentierende mittlere Sonne läuft auf dem Himmelsäquator um. Die Zeitgleichung ist proportional zur Differenz zwischen den Rektaszensionen ${\displaystyle \alpha _{M}}$ der fiktiven mittleren und ${\displaystyle \alpha }$ der realen wahren Sonne.

Zweite Definition:

${\displaystyle \mathrm {(11)} \quad {\text{ZG}}=4(\alpha _{M}-\alpha )\quad [{\text{min}}]}$

Der Faktor 4 ergibt sich daraus, dass zwei Himmelskörper mit 1° Rektaszensionsdifferenz den Meridian im zeitlichen Abstand von 4 Minuten passieren. Die Reihenfolge der beiden Subtraktionsterme hat sich umgekehrt, weil die Richtungen für den Stundenwinkel ${\displaystyle \tau }$ (ihm entsprechen WOZ und MOZ) und die Rektaszension ${\displaystyle \alpha }$ zueinander entgegengesetzt definiert sind.

Vorgehensweise

Der zu einem bestimmten Zeitpunkt ${\displaystyle t}$ zu ermittelnden Rektaszension ${\displaystyle \alpha }$ (Gleichung ${\displaystyle \mathrm {(11)} }$) der Sonne entspricht in heliozentrischer Betrachtung die äquatoriale Länge der Erde, die auf einfache Weise aus ihrer ekliptikalen Länge ${\displaystyle \lambda }$ (zweite der nebenstehenden Abbildungen) errechenbar ist. Mit Hilfe der Kepler-Gleichung wird die wahre Anomalie ${\displaystyle V}$ (erste der nebenstehenden Abbildungen) ermittelt, aus der dann durch Bezugspunktänderung ${\displaystyle \lambda }$ bestimmt wird.

Anwendung der Kepler-Gleichung

Momentane Anomalien der Erde (zum Zeitpunkt t) auf ihrer elliptischen Bahn um die Sonne:
V – wahre, M – mittlere, E – exzentrische Anomalie
B – Sonne, X – Erde, Y – fiktive Erde, P – Perihel, A – Aphel, K – 1.-Jan.-Punkt
unten links: V und M als Funktionen der Zeit

Übergang von der heliozentrischen elliptischen Erdbahn (links, mit wahrer Erde X und fiktiver Erde Y) zur geozentrischen Sonnenbahn (Ekliptikkreis, rechts, mit wahrer Sonne S und fiktiver Sonne S')

„Herunterholen“ der wahren Sonne auf den Äquator: Ermittlung ihres Rektaszensionswinkels α aus ihrem ekliptikalen Längenwinkel λ
S″: mittere Sonne auf dem Äquator

Mittlere Anomalie:

Die in Gleichung ${\displaystyle \mathrm {(1)} }$ allgemein formulierte mittlere Anomalie lautet im Zusammenhang mit der Zeitgleichung:

${\displaystyle \mathrm {(12)} \quad M(t)={\frac {360^{\circ }}{J_{\text{an}}}}\cdot (t-t_{P})}$

${\displaystyle J_{\text{an}}}$: anomalistisches Jahr zwischen zwei Passagen des Perihels

${\displaystyle t_{P}}$: Zeitpunkt der Perihel-Passage

Bei Periheldurchgang h​at die mittlere Anomalie folgenden Wert:

${\displaystyle \mathrm {(13)} \quad M_{0}=-{\frac {360^{\circ }}{J_{\text{an}}}}\cdot t_{P}}$

Bei der Zeitgleichung ist es üblich, die Werte eines Kalenderjahres im entsprechenden Astronomischen Jahrbuch zu veröffentlichen. Der 1. Januar 12:00 (UT) des entsprechenden Jahres wird als Nullpunkt für ${\displaystyle t}$ verwendet, sodass gegenwärtig für ${\displaystyle t_{P}}$ etwa 2 bis 3 Tage und daraus für ${\displaystyle M_{0}}$ etwa 2° bis 3° gelten.[11] Es hat sich bequemerweise eingebürgert, den jeweils neuen Wert für ${\displaystyle M_{0}}$ als eine sogenannte Jahreskonstante im Voraus zu veröffentlichen.

Mit ${\displaystyle M_{0}}$ und ${\displaystyle t}$ ab 1. Januar 12:00 (UT) wird aus Gleichung (12):

${\displaystyle M(t)=M_{0}+{\frac {360^{\circ }}{J_{\text{an}}}}\cdot t}$

Kepler-Gleichung:

${\displaystyle M(t)=E(t)-{\frac {180^{\circ }}{\pi }}\cdot e\cdot \sin E(t)}$

Mit der dem vorgegebenen Zeitpunkt entsprechenden mittleren Anomalie ${\displaystyle M}$ und der Erdbahn-Exzentrizität ${\displaystyle e}$ wird mit Hilfe der Kepler-Gleichung die exzentrische Anomalie ${\displaystyle E}$ ermittelt.

Wahre Anomalie:

Bei der Behandlung der Zeitgleichung wird für die wahre Anomalie meistens das Formelzeichen ${\displaystyle V}$ (anstatt ${\displaystyle T}$ wie oben) verwendet.

Die exzentrische Anomalie ${\displaystyle E}$ führt in einer rein geometrischen Betrachtung in der Ellipse und in ihrem Umkreis (erste der nebenstehenden Abbildungen) wie folgt zur wahren Anomalie ${\displaystyle V}$:[12]

${\displaystyle \tan \left({\frac {V(t)}{2}}\right)=\kappa \cdot \tan \left({\frac {E(t)}{2}}\right)}$

${\displaystyle \mathrm {(14)} \quad \kappa ={\sqrt {\frac {1+e}{1-e}}}}$ … eine Ellipsenkonstante

Das Kepler-Problem i​st mit d​er Ermittlung d​er wahren Anomalie d​er Erde gelöst. Im Folgenden w​ird die Ermittlung d​er Zeitgleichung abgeschlossen.

Wahre Anomalie der Erde → Rektaszension der Sonne

Wahre Anomalie d​er Erde → ekliptikale Länge d​er Erde → ekliptikale Länge d​er Sonne:

Von der Erde aus gesehen spiegelt sich die Bewegung der Erde um die Sonne wider in der scheinbaren Bewegung der Sonne in der Ekliptik, dem Schnitt der Erdbahnebene mit der um die Erde als Mittelpunkt geschlagenen Richtungskugel (siehe zweite der nebenstehenden Abbildungen).[13][14] Die ekliptikale Länge der Erde und die ekliptikale Länge der Sonne sind somit Synonyme mit dem Formelzeichen ${\displaystyle \lambda .}$

Bezugspunkt für die ekliptikale Länge (und auch der Rektaszension) ist gemäß allgemeinem Brauch der Frühlingspunkt. Die ekliptikale Länge ${\displaystyle \lambda (t)}$ der Sonne wird erhalten, indem dem auf das Perihel der Erdbahn bezogenen Winkel ${\displaystyle V(t)}$ der Winkel ${\displaystyle L}$ zwischen Perihel P und dem dem Frühlingspunkt entsprechenden Ort (F) addiert wird:[15]

${\displaystyle \mathrm {(15)} \quad \lambda (t)=V(t)+L}$

Der Wert von ${\displaystyle L}$ ist negativ. Unter den nahezu konstanten Grundgrößen ist ${\displaystyle L}$ diejenige, die sich mit der Zeit wegen der langsamen Annäherung zwischen Frühlingspunkt bzw. Punkt (F) und Perihel am stärksten verändert. Sie wird deshalb nicht nur jährlich als sogenannte Jahreskonstante ${\displaystyle L_{0}}$ neu gesetzt, sondern mit folgender Gleichung permanent verändert:

${\displaystyle \mathrm {(16)} \quad L(t)=L_{0}+{\tfrac {0{,}0172^{\circ }}{J_{\text{tr}}}}\cdot t}$

Frühlingspunkt und Perihel nähern sich mit ${\displaystyle \approx {\tfrac {0{,}0172^{\circ }}{J_{\text{tr}}}}.}$ ${\displaystyle J_{\text{tr}}}$ ist das tropische Jahr (Zeit für zwei aufeinanderfolgende Passagen des Frühlingspunkts bzw. des Punktes (F)). Unter Beachtung der Gleichung ${\displaystyle \mathrm {(16)} }$ ist statt Gleichung ${\displaystyle \mathrm {(15)} }$ zu schreiben:

${\displaystyle \mathrm {(17)} \quad \lambda (t)=V(t)+L_{0}+{\tfrac {0{,}0172^{\circ }}{J_{\text{tr}}}}\cdot t}$

Der Wert von ${\displaystyle L_{0}}$ ist negativ.

Ekliptikale Länge d​er Sonne → Rektaszension d​er Sonne:

Neben d​er Elliptizität d​er Erdbahn verursacht d​ie zur Erdbahnebene n​icht rechtwinklige Lage d​er Erdachse u​nd ihre Richtungsänderung relativ z​ur Sonne d​ie Zeitgleichung.

Die Rektaszension ${\displaystyle \alpha }$ der Sonne lässt sich z. B. mit allgemein bekannten Transformationsgleichungen oder mit folgender einfachen Beziehung im entsprechenden rechtwinkligen sphärischen Dreieck (siehe dritte der nebenstehenden Abbildungen) aus der ekliptikalen Länge ${\displaystyle \lambda }$ ermitteln:

${\displaystyle \mathrm {(18)} \quad \alpha (t)=\arctan(\tan \lambda (t)\cdot \cos \varepsilon )}$

${\displaystyle \varepsilon }$ ist die Schiefe der Erdachse: ${\displaystyle \varepsilon =23{,}44^{\circ }}$ .

Rektaszension der mittleren Sonne

Die Bewegung d​er mittleren Sonne S″ (dritte d​er rechts stehenden Abbildungen) a​uf dem Äquator m​acht die gleichmäßig vergehende Zeit gleich w​ie die d​er auf d​er Erdbahn umlaufenden fiktiven Erde (Punkt Y) anschaulich. Ihr Lauf i​st möglichst e​ng an d​en der wahren Sonne z​u koppeln, d​amit sie d​eren Lauf e​twa „mittelt“. Das w​urde mit folgender Definition erreicht:[16]

${\displaystyle \mathrm {(19)} \quad \alpha _{M}(t)=L(t)+M(t)}$

Wenn man die zeitliche Änderung von ${\displaystyle L}$ vernachlässigt, gilt auch:

${\displaystyle \left(20\right)\quad \alpha _{M}(t)=L_{0}+M_{0}+{\tfrac {360^{\circ }}{J_{\text{tr}}}}\cdot t}$

Zeitgleichung

Die beiden zur Anwendung der Zeitgleichung ${\displaystyle \mathrm {(11)} }$ erforderlichen Rektaszensionen ${\displaystyle \alpha _{M}}$ und ${\displaystyle \alpha }$ sind gefunden.

${\displaystyle \mathrm {(11)} \quad {\text{ZG}}=4(\alpha _{M}-\alpha )\quad [{\text{min}}]}$

Rechenbeispiel

Die Zeitgleichung für d​en 2. April 2015, 12:00 UT (t = 91 Tage) s​ei zu berechnen.

Die Jahreskonstanten 2015 sind:[16][17][18]

${\displaystyle M_{0}=-2{,}3705^{\circ }}$

${\displaystyle J_{an}=365{,}259991{\text{ Tage}}}$

${\displaystyle J_{tr}=365{,}242907{\text{ Tage}}}$

${\displaystyle e=0{,}016703}$

${\displaystyle \varepsilon =23{,}43734^{\circ }}$

${\displaystyle L_{0}=-76{,}8021^{\circ }}$

Die Rechnungen sind:

${\displaystyle M(t)=M_{0}+{\frac {360^{\circ }}{J_{\text{an}}}}\cdot t=87{,}3190^{\circ }}$

${\displaystyle \mathrm {(16)} \quad L(t)=L_{0}+{\tfrac {0{,}0172^{\circ }}{J_{\text{tr}}}}\cdot t=-76{,}7978^{\circ }}$

${\displaystyle M(t)=E(t)-{\frac {180^{\circ }}{\pi }}\cdot e\cdot \sin E(t)\quad \rightarrow \quad E(t)=88{,}2756^{\circ }}$

${\displaystyle \textstyle \kappa ={\sqrt {\frac {1+e}{1-e}}}=1{,}0168445}$

${\displaystyle \tan \left({\frac {V(t)}{2}}\right)=\kappa \cdot \tan \left({\frac {E(t)}{2}}\right)\quad \rightarrow \quad V(t)=89{,}2325^{\circ }}$

${\displaystyle \mathrm {(17)} \quad \lambda (t)=V(t)+L(t)=12,4347^{\circ }}$

${\displaystyle \mathrm {(18)} \quad \alpha (t)=\arctan(\tan \lambda (t)\cdot \cos \varepsilon )=11{,}4369^{\circ }}$

${\displaystyle \mathrm {(19)} \quad \alpha _{M}(t)=L(t)+M(t)=10{,}5212^{\circ }}$

${\displaystyle \mathrm {(11)} \quad {\text{ZG}}(t)=4{\frac {\text{min}}{^{\circ }}}\cdot (\alpha _{M}(t)-\alpha (t))=-3{,}6629{\text{ min}}=-3{\text{ min}}{\text{ 40}}{\text{ sec}}}$

Die Zeitgleichung h​at am 2. April 2015, 12:00 UT d​en Wert:

${\displaystyle {\text{ZG}}(t=91{\text{ Tage}})=-3{\text{ min}}{\text{ 40}}{\text{ sec}}}$

Zeitgleichungswerte für die Passage ausgezeichneter Bahnpunkte

Vom Kalender und damit von der Jahreskonstanten ${\displaystyle M_{0}}$ unabhängig sind Zeitgleichungswerte für die Passage ausgezeichneter Punkte durch die Erde auf ihrer Bahn (beziehungsweise durch die Sonne auf der Ekliptik): Frühlings-, Sommer-, Herbst- und Winteranfangspunkt, Perihel und Aphel.

Zeitgleichungswerte und Zeitpunkte für die Passage ausgezeichneter Bahnpunkte *)

	F-Anfang	S-Anfang	H-Anfang	W-Anfang	Perihel	Aphel
λ/°	0	90	180	270	L0	L0 + 180
ZG/min	−7,44	−1,74	+7,48	+1,70	−4,50	−4,50
tP/d **)	76,234	168,990	262,641	352,485	0	182,621

 *) Die Werte gelten für das Jahr 2004 mit L₀ = −76,99° und J_tr =365,2428 Tage.[16]
**) Die angegebenen Zeiten beziehen sich auf den Periheldurchgang, nicht wie in obigem Beispiel auf den 1. Januar 12:00 UT.

Ihre Berechnung ist einfacher, als die für beliebige Zeitpunkte, weil die Kepler-Gleichung ${\displaystyle E=f(M)}$ nicht gelöst werden muss. Von der vorgegebenen ekliptikalen Länge ${\displaystyle \lambda }$ eines der ausgezeichneten Punkte ist leicht zur wahren (Gl. ${\displaystyle \mathrm {(15)} }$)[19] und weiter zur exzentrischen Anomalie zu finden. Aus Letzterer folgt mit der umgestellten Kepler-Gleichung ${\displaystyle M=f(E)}$ die mittlere Anomalie, also der Bahnpunkt der fiktiven mittleren Erde. Die ekliptikale Länge des Perihels[19] zu Letzterer addiert (Gl. ${\displaystyle \mathrm {(19)} }$) ist die gesuchte mittlere Rektaszension ${\displaystyle \alpha _{M}}$ (Minuend in der Zeitgleichung ${\displaystyle \mathrm {(11)} }$). Die wahre Rektaszension ${\displaystyle \alpha }$ (Subtrahend) ist bei den Punkten Frühling bis Winter mit deren ekliptikaler Länge ${\displaystyle \lambda }$ identisch. Nur bei den Punkten Perihel und Aphel ergibt die Koordinatentransformation (Gl. ${\displaystyle \mathrm {(18)} }$) kleine Werteunterschiede.

Bei der Vorgehensweise, die Berechnung mit einer vorgegebenen ekliptikalen Länge bzw. einer vorgegebenen wahren Anomalie zu beginnen, erhält man neben der Zeitgleichung auch die seit der Perihelpassage der Erde vergangene Zeit. Das ist die Zeit, die die mittlere Anomalie repräsentiert und sie wird aus dem Zwischenergebnis für die mittlere Anomalie ${\displaystyle M}$ mit Hilfe der entsprechend umzustellenden Gleichung ${\displaystyle \mathrm {(12)} }$ errechnet.

Diese Vorgehensweise w​ird gelegentlich a​uch für d​ie allgemeine Arbeit empfohlen, Zeitgleichungstabellen z​u ermitteln.[20] Man erspart s​ich dabei d​as aufwändige Lösen d​er Kepler-Gleichung, findet z​u Werten für gewünschte Zeitpunkte a​ber nur d​urch Probieren o​der bei genügender Ergebnisdichte d​urch Interpolieren.

Literatur

• Andreas Guthmann: Einführung in die Himmelsmechanik und Ephemeridenrechnung. BI-Wiss.-Verl., Mannheim 1994, ISBN 3-411-17051-4.
• Peter Colwell: Solving Kepler's equation over three centuries. Hrsg.: Willmann-Bell. Richmond, VA 1993, ISBN 0-943396-40-9, S. 202.

Weblinks

• J. M. Danby, T. M. Burkardt: The solution of Kepler's equation. I (= Cel. Mech. Band 31). 1983, S. 95–107, doi:10.1007/BF01686811, bibcode:1983CeMec..31...95D.
• B. A. Conway: An improved algorithm due to Laguerre for the solution of Kepler's equation. 1986, doi:10.2514/6.1986-84.
• Seppo Mikkola: A cubic approximation for Kepler's equation (= Cel. Mech. Band 40, Nr. 3). 1987, doi:10.1007/BF01235850, bibcode:1987CeMec..40..329M.
• Albert Nijenhuis: Solving Kepler's equation with high efficiency and accuracy (= Cel. Mech. Dyn. Astr. Band 51, Nr. 4). 1991, S. 319–330, doi:10.1007/BF00052925, bibcode:1991CeMDA..51..319N.
• Toshio Fukushima: A method solving kepler's equation without transcendental function evaluations (= Cel. Mech. Dyn. Astron. Band 66, Nr. 3). 1996, S. 309–319, doi:10.1007/BF00049384, bibcode:1996CeMDA..66..309F.
• E. D. Charles, J. B. Tatum: The convergence of Newton-Raphson iteration with Kepler's equation (= Cel. Mech. Dyn. Astr. Band 69, Nr. 4). 1997, S. 357–372, doi:10.1023/A:1008200607490, bibcode:1997CeMDA..69..357C.
• Laura Stumpf: Chaotic behaviour in the newton iterative function associated with kepler's equation (= Cel. Mech. Dyn. Astr. Band 74, Nr. 2). 1999, S. 95–109, doi:10.1023/A:1008339416143.
• M. Palacios: Kepler equation and accelerated Newton method (= J. Comp. Appl. Math. Band 138). 2002, S. 335–346, doi:10.1016/S0377-0427(01)00369-7, bibcode:2002JCoAM.138..335P.
• John P. Boyd: Rootfinding for a transcendental equation without a first guess: Polynomialization of Kepler's equation through Chebyshev polynomial equation of the sine (= Appl. Num. Math. Band 57, Nr. 1). 2007, S. 12–18, doi:10.1016/j.apnum.2005.11.010.
• Eric W. Weisstein: Kepler's Equation. In: MathWorld (englisch).

Einzelnachweise

1. J.-L. Lagrange, Sur le problème de Kepler, in Mémoires de l'Académie Royale des Sciences de Berlin, vol. 25, 1771, Seiten 204–233
2. Peter Colwell: Bessel functions and Kepler's equation (= Amer. Math. Monthly. Band 99, Nr. 1). Januar 1992, S. 45–48 (englisch).
3. § II.6.67 Numerische Verfahren. Guthmann, S. 128 f.
4. § II.6.66 Reihenentwicklung der exzentrischen Anomalie. Guthmann, S. 125 ff.
5. Siegfried Wetzel: Die Zeitgleichung für Nicht-Astronomen. Deutsche Gesellschaft für Chronometrie, Mitteilungen Nr. 111, Herbst 2007, Anhang 3.
6. R. Strebel: Die Keplersche Gleichung. (Memento des Originals vom 13. August 2011 im Internet Archive)  Info: Der Archivlink wurde automatisch eingesetzt und noch nicht geprüft. Bitte prüfe Original- und Archivlink gemäß Anleitung und entferne dann diesen Hinweis. Oktober 2001, Kap. 1.3 und 5.1.
7. Aufgaben zu § II.5. Guthmann, S. 122 f.
8. 10. und 11. Aufgabe zu § II.5. Guthmann, S. 123.
9. § II.5.58 Der Hodograph. Guthmann, S. 114 f.
10. Die hier verwendeten Formelzeichen sind die gleichen wie in:
    Sonnenuhren-Handbuch, Berechnung der Zeitgleichung. Deutsche Gesellschaft für Chronometrie e.V., Fachkreis Sonnenuhren, 2006, S. 43–49.
11. Wegen der Schalttagregelung im Kalender schwanken beide Werte innerhalb der Vierjahresperiode schwach: Δt_P ≈ ¾Tag, ΔM₀ ≈ ¾°.
12. Siegfried Wetzel: Die Zeitgleichung für Nicht-Astronomen. (Memento vom 7. April 2014 im Internet Archive) Deutsche Gesellschaft für Chronometrie, Mitteilungen Nr. 111, Herbst 2007, Anhang 3.
13. Manfred Schneider: Himmelsmechanik, Band II: Systemmodelle. BI-Wissenschaftsverlag, 1993, ISBN 3-411-15981-2, S. 507.
14. Dieser Zusammenhang erlaubt umgekehrt, die ekliptikale Länge ${\displaystyle \lambda }$ und den Frühlingspunkt F als Bezugspunkt (sowohl für ${\displaystyle L}$ als auch für ${\displaystyle \alpha }$) auf die Erdbahn zurückzuspiegeln (siehe nebenstehende Abbildung, rechts → links).
15. Zeichen für Winkeldifferenz und Ort in nebenstehender Abbildung in Klammern gesetzt, da Winkel und Ort für den Gebrauch auf der Erdbahn nicht definiert sind.
16. Sonnenuhren-Handbuch, 3.3 Berechnung der Zeitgleichung. Deutsche Gesellschaft für Chronometrie e.V., Fachkreis Sonnenuhren, 1900.
17. Diese „Basiswerte“ gelten für den 1. Januar 2015 12:00 UT. Ihre langsame Veränderung wird im Folgenden während des gesamten Jahres 2015 nicht beachtet. Die in dieser Zeit kumulierte Veränderung schlägt sich erst in den Jahreskonstanten 2016 nieder. Ausnahme ist ${\displaystyle L_{0}}$. Gleichung ${\displaystyle \mathrm {(8)} }$ enthält die permanente Veränderung ${\displaystyle L(t)}$.
    Die Hochrechnung der Jahreskonstanten erfolgt mit den Basiswerten der Jahre 2000 bzw. 1900 wie folgt (DGC-Handbuch, S. 47):

    ${\displaystyle M_{0}=357{,}5256^{\circ }+35999{,}0498^{\circ }\cdot T/36525}$

    ${\displaystyle J_{tr}=(365{,}24219878+6{,}16\cdot 10^{-8}\cdot J){\text{ Tage}}}$

    ${\displaystyle J_{an}=(365{,}25964124+3{,}04\cdot 10^{-8}\cdot J){\text{ Tage}}}$

    ${\displaystyle e_{0}=0{,}016709-4{,}2\cdot 10^{-7}\cdot T/36525}$

    ${\displaystyle \varepsilon _{0}=23{,}439291^{\circ }-0{,}013004^{\circ }\cdot T/36525}$

    ${\displaystyle L_{0}=282{,}9400^{\circ }+1{,}7192^{\circ }\cdot T/36525}$

    ${\displaystyle T}$ ist die Zahl der Tage seit 1. Januar 2000 12:00 UT; ${\displaystyle J}$ ist die Zahl der Jahre seit 1900. Bei den Winkeln ${\displaystyle M_{0}}$ und ${\displaystyle L_{0}}$ ist modulo 360° zu rechnen, und sie müssen zwischen −180° und +180° liegen.
18. Die Jahreskonstanten (z. B. für 2015) werden hier so bezeichnet, weil sie nur für das eine Jahr benutzt werden, auf das sie sich beziehen. Darüber hinaus gelten sie ohne bedeutsame Einbuße an Genauigkeit der Zeitgleichung auch für Termine in fernliegenden Jahren (z. B. für 2050 oder 1950). Die Zeit ${\displaystyle t}$ nimmt dann entsprechend hohe positive bzw. negative Werte an; das gegebene Rechenschema bleibt aber unverändert anwendbar. Bei der Bestimmung von ${\displaystyle V}$ und ${\displaystyle \alpha }$ sind die Nebenwerte des Arkustangens zu verwenden, die ${\displaystyle E/2}$ bzw. ${\displaystyle \lambda }$ am nächsten liegen.
19. Dabei wird mit der ekliptikalen Länge L = L₀ des Perihels gerechnet, was ausreichend genau und wegen der nicht bekannten Zeit t auch nicht anders möglich ist.
20. Heinz Schilt: Zur Berechnung der mittleren Zeit für Sonnenuhren. Schriften der Freunde alter Uhren, 1990.