Astronomische Berechnungen für Amateure/ Himmelsmechanik/ Sonne

Aus Wikibooks

Wechseln zu: Navigation, Suche

Die VSOP87-Theorie gibt – mit Ausnahme der Serie E – keine Entwicklung für die Koordinaten der Sonne bezogen auf ein Koordinatensystem, das fest mit der Erde verbunden ist. Dies ist auch nicht nötig: ist für einen bestimmten Zeitpunkt t der Positionsvektor r SE von der Sonne zur Erde bekannt, so ist die Position der Sonne im erdbezogenen Koordinatensystem durch r ES = –r SE gegeben. Allerdings ist sorgfältig darauf zu achten, welches Koordinatensystem gemeint ist. r SE wird in der Regel entweder auf den Schwerpunkt des Erd-Mond-Systems (Baryzentrum) oder aber auf den Erdmittelpunkt (Geozentrum) bezogen, während wir für Beobachtungen in der Regel den Beobachtungsort auf der Erdoberfläche benötigen (Topozentrum). Im weiteren ist darauf zu achten, für welches Äquinoktium die Koordinatensysteme festgelegt sind – sie müssen darin übereinstimmen. In einem späteren Kapitel werden wir zeigen, wie auf diese Weise Sonnenephemeriden berechnet werden können.

Die Erdbahn – und damit nach dem eben Gesagten auch die scheinbare Sonnenbahn – ist durch eine kleine Exzentrizität gekennzeichnet: e beträgt nur 0.0167. In diesem Fall ist es nicht nötig, die Keplergleichung zu lösen, um die wahre Anomalie v bestimmen zu können. Durch eine Näherungsformel kann sie direkt berechnet werden (auf der Basis einer Reihenentwicklung):


 \begin{align} v \; - \; M \quad  \equiv \quad C \quad & = \quad \left( 2 e \; - \; \frac{e^{3}}{4} \; + \; \frac{5}{96} e^{5} \right) \cdot \sin M \quad \\ &+ \quad \left( \frac{5}{4} \; e^{2} \; - \; \frac{11}{24} \; e^{4}  \right) \cdot \sin (2 M) 
 \qquad \\ &+ \; \left( \frac{13}{12} \; e^3 \; - \; \frac{43}{64} \; e^5 \right) \cdot \sin (3 M) \quad \\ &+ \quad \frac{103}{96} \; e^4 \cdot \sin (4 M) \quad + \quad \frac{1097}{960} \; e^5 \cdot \sin (5 M) \end{align}



Diese Gleichung für die Grösse C ist bekannt unter dem Namen Mittelpunktsgleichung, das Ergebnis ist im Bogenmass. Die Grösse C gibt an, wieviel ein Planet auf einer Keplerbahn einem fiktiven Himmelskörper vorauseilt oder hinterher rennt, der die Sonne auf einer Kreisbahn gleichmässig umläuft.

Für den Bahnradius r existiert ebenfalls eine Reihenentwicklung:


\begin{align} \frac{r}{a} \quad &= \quad 1 \quad + \quad \frac{e^2}{2} \quad - \quad \left( e \; - \; \frac{3}{8} \; e^3 \; + \; \frac{5}{192} \; e^5 \right) \; \cos M  \quad \\ & - \quad \left( \frac{e^2}{2} \; - \; \frac{e^4}{3} \right) \; \cos (2 M) 
\; \\ & - \; \left( \frac{3}{8} \; e^3  \; - \; \frac{45}{128} \; e^5 \right) \; \cos (3 M) \quad \\ & -  \quad \frac{e^4}{3} \; \cos (4 M) \quad - \quad \frac{125}{384} \; e^5 \; \cos (5 M) \end{align}


Ist die Länge des Perigäums der Sonnenbahn in der Ekliptik bekannt, lassen sich damit zumindest genäherte Sonnenpositionen berechnen.

Persönliche Werkzeuge