Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen der Seite angezeigt.

--- loesung_der_keplergleichung [2024/05/14 13:50] – hcgreier
+++ loesung_der_keplergleichung [2025/08/24 02:20] (aktuell) – hcgreier
@@ Zeile 7: / Zeile 7: @@
   - Berechnen der heliozentrischen Positionen aus den mittleren Bahnelementen des Körpers (Keplerbahnen). Dies ist zwar die ungenaueste Methode, liefert aber sehr schnelle Ergebnisse und ist für manche Anwendungen genau genug. Man könnte z.B. eine Grafik mit den aktuellen Planetenpositonen erzeugen, wobei natürlich nicht auf Bogensekunden genau gerechnet werden muss, da dies in der Grafik gar nicht darstellbar ist.
-In diesem Kapitel geht es um den dritten Fall. Es soll die wahre Anomalie $\nu$ des Objekts berechnet werden. Dies kann entweder durch Lösung der Keplergleichung erreicht werden oder – wenn die Bahnexzentrizität $\epsilon$ nicht zu groß ist – durch eine Reihenentwicklung (Mittelpunktsgleichung) oder eine einfache Näherungsformel.
+In diesem Kapitel geht es um den dritten Fall. Es soll die **wahre Anomalie** $\nu$ des Objekts berechnet werden. Dies kann entweder durch Lösung der Keplergleichung erreicht werden oder – wenn die Bahnexzentrizität $\epsilon$ nicht zu groß ist – durch eine Reihenentwicklung (Mittelpunktsgleichung) oder eine einfache Näherungsformel.
 ===== Begriffe und Zusammenhänge =====
@@ Zeile 40: / Zeile 40: @@
 Diese Größe ist dimensionslos und liegt bei einer Ellipse immer zwischen 0 und 1. Die Exzentrizität $\epsilon = 0$ wäre bei einem Kreis gegeben. Der Perihelabstand und Aphelabstand werden mit $ q $ bzw. $ Q $ bezeichnet. Befindet sich der Körper im Perihel, ist $ \nu = 0^{\circ} $ und $ r = q $, während im Aphel $ \nu = 180^{\circ} $ ist und $ r = Q $.
-Nun betrachtet man in **Abb. 2** einen fiktiven Planeten oder Kometen $K'$, der eine //Kreisbahn// um die Sonne beschreibt, und sich mit **konstanter** Geschwindigkeit bewegt und die gleiche Umlaufzeit hat, wie der reale Planet oder Komet $K$ auf der elliptischen Bahn. Darüber hinaus wird angesetzt, dass sich das fiktive Objekt im Punkt $P'$ befindet, wenn der reale Körper im Perihel $P$ steht. Als Hilfskreis wird der Scheitelkreis der Hautpachse $a$  genommen, auf dem sich der fiktive Planet bewege. Einige Zeit später, wenn sich der wahre Körper in $K$ befindet, ist das fiktive Objekt bis nach $K'$ gelaufen. Wie bereits angemerkt ist der Winkel $\nu = \angle PSK$ die **wahre Anomalie** $\nu$ des Objekts zum gegebenen Zeitpunkt. Zur gleichen Zeit durchläuft $K'$ den Winkel $\angle P'ZK'$, und dies ist die **mittlere Anomalie**, die im allgemeinen mit $M$ bezeichnet wird.
+Nun betrachtet man in **Abb. 2** einen fiktiven Planeten oder Kometen $K'$, der eine //Kreisbahn// um die Sonne beschreibt, und sich mit **konstanter** Geschwindigkeit bewegt und die **gleiche Umlaufzeit** hat, wie der reale Planet oder Komet $K$ auf der elliptischen Bahn. Darüber hinaus wird angesetzt, dass sich das fiktive Objekt im Punkt $P'$ befindet, wenn der reale Körper im Perihel $P$ steht. Als Hilfskreis wird der Scheitelkreis der Hauptachse $a$ genommen, auf dem sich der fiktive Planet bewege. Einige Zeit später, wenn sich der wahre Körper in $K$ befindet, ist das fiktive Objekt bis nach $K'$ gelaufen. Wie bereits angemerkt ist der Winkel $\nu = \angle PSK$ die **wahre Anomalie** $\nu$ des Objekts zum gegebenen Zeitpunkt. Zur gleichen Zeit durchläuft $K'$ den Winkel $\angle P'ZK'$, und dies ist die **mittlere Anomalie**, die im allgemeinen mit $M$ bezeichnet wird.
 <imgcaption image2|Maßgebliche Winkel für die Keplergleichung>{{ :keplergleichung_winkel.png?800 |}}</imgcaption>
@@ Zeile 48: / Zeile 48: @@
 Definitionsgemäß wächst der Winkel $M$ linear (gleichförmig) mit der Zeit. Der Wert von $M$ ist leicht zu finden, da im Perihel $M = 0^{\circ}$ gilt und der Winkel während eines vollständigen Umlaufes um exakt $360^{\circ}$ zunimmt.
-Das eigentliche Problem besteht aber darin, die wahre Anomalie $\nu$ zu finden, wenn die mittlere Anomalie $M$ und die Bahnexzentrizität $\epsilon$ bekannt sind. Falls man nicht auf eine Reihenentwicklung zurückgreifen will, ist man gezwungen, die Keplergleichung zu lösen. Hierzu ist es nötig, einen Hilfswinkel $E$, die **exzentrische Anomalie**, einzuführen. Dessen geometrische Definition ist in **Abb. 2** angegeben. Der äußere gestrichelte Hilfskreis hat den Durchmesser $2a = \overline{AP}$ der großen Achse. Man zeichnet $\overline{HK}$ rechtwinklig zu $\overline{AP}$. Den Winkel $\angle PZH$ nennt man die **exzentrische Anomalie** $E$.
+Das eigentliche Problem besteht aber darin, die wahre Anomalie $\nu$ zu finden, wenn die mittlere Anomalie $M$ und die Bahnexzentrizität $\epsilon$ bekannt sind. Falls man nicht auf eine Reihenentwicklung zurückgreifen will, ist man gezwungen, die Keplergleichung zu lösen. Hierzu ist es nötig, einen Hilfswinkel $E$, die **exzentrische Anomalie**, einzuführen. Deren geometrische Definition ist in **Abb. 2** angegeben. Der äußere gestrichelte Hilfskreis hat den Durchmesser $2a = \overline{AP}$ der großen Achse. Man zeichnet $\overline{HK}$ rechtwinklig zu $\overline{AP}$. Den Winkel $\angle PZH$ nennt man die **exzentrische Anomalie** $E$.
-Steht der Planet im Perihel, dann sind die Winkel $\nu = E = M = 0^{\circ}$. In der Nähe des Perihels hat der wahre Planet eine größere Geschwindigkeit als der mittlere, fiktive Planet. Diese ist dem 2. Keplerschen Gesetz geschuldet, wonach ein Körper in gleichen Zeiten mit dem Radiusvektor $r$ dieselben Flächen überfährt, und sich der Körper deshalb in Perihelnähe schneller bewegt als in Aphelnähe. Deshalb gilt zwischen Perihel und Aphel (wenn sich der Planet von der Sonne entfernt)
+Steht der Planet im Perihel, dann sind die Winkel $\nu = E = M = 0^{\circ}$. In der Nähe des Perihels hat der wahre Planet eine größere Geschwindigkeit als der mittlere, fiktive Planet. Dies ist dem 2. Keplerschen Gesetz geschuldet, wonach ein Körper in gleichen Zeiten mit dem Radiusvektor $r$ dieselben Flächen überfährt, und sich der Körper deshalb in Perihelnähe schneller bewegt als in Aphelnähe. Deshalb gilt zwischen Perihel und Aphel (wenn sich der Planet von der Sonne entfernt)
 $$\nu \gt M$$
@@ Zeile 170: / Zeile 170: @@
 Es zeigt sich, dass derart viele Iterationen nur bei sehr großen Werten für $\epsilon$ auftreten, die sehr nahe an $1$ liegen. Die höchste Exzentrizität der Planeten hat Merkur* mit etwa $\epsilon = 0.205631^{\circ}$ (für Epoche $J2000$). Die Berechnung für alle Planeten mit dieser Methode sollte daher keine Probleme bereiten. Für Kometen, die eine sehr hohe Exzentrizität haben können, empfiehlt sich diese Methode nicht.
+{{anchor:kepler1_iter}}
 <WRAP center round info 100%>
 * Pluto hat eine größere Exzentrizität als Merkur mit $\epsilon_{Pluto} = 0.2488^{\circ}$, gilt aber seit 2006 nicht mehr als Planet (Zwergplanet).
@@ Zeile 293: / Zeile 293: @@
 <WRAP center round tip 100%>
-Es sei darauf hingewiesen, dass eine Genauigkeit von $10^{-9}$ absurd ist, das entspricht etwa 0.000004 Bogensekunden! Der Grenzwert von $10^{-9}$ wurde hier zu Anschauungszwecken gewählt. Er lässt sich aber in der Funktion nach Belieben setzen.
+Es sei darauf hingewiesen, dass eine Genauigkeit von $10^{-9}$ Grad absurd ist, das entspricht etwa 0.000004 Bogensekunden! Der Grenzwert von $10^{-9}$ wurde hier zu Anschauungszwecken gewählt. Er lässt sich aber in der Funktion nach Belieben setzen.
 </WRAP>
@@ Zeile 451: / Zeile 451: @@
 <WRAP center round tip 100%>
-Gegebenenfalls kann diese Näherungsformel einen **Startwert** für die Methoden 1/2 liefern.
+Gegebenenfalls kann diese Näherungsformel einen **Startwert** für die Methoden 1 bzw. 2 liefern.
 </WRAP>
@@ Zeile 480: / Zeile 480: @@
 <imgcaption image9|Abweichung der Methode 4 von der Lösung mittels Keplergleichung für die Planeten>{{ :keplergleichung_diff_kepler_approx_planets_big.png?800 |}}</imgcaption>
-Größere Abweichungen gibt es für die Planeten Merkur, Mars und Saturn. Man sieht, das vor allem Merkur ausreißt, er hat auch die größte Exzentrizität. Die größte Abweichung liegt bei $306\overset{''}{.}91$ beim Winkel $70\overset{\circ}{.}39$. Zur Erstellung einer Grafik wäre das aber immer noch ausreichend, da $300'' = 5' = 0\overset{\circ}{.}08333$ immer noch klein ist und in einer Grafik praktisch nicht sichtbar wäre.
+Größere Abweichungen gibt es für die Planeten Merkur, Mars und Saturn. Man sieht, dass vor allem Merkur ausreißt, er hat auch die größte Exzentrizität. Die größte Abweichung liegt bei $306\overset{''}{.}91$ beim Winkel $70\overset{\circ}{.}39$. Zur Erstellung einer Grafik wäre das aber immer noch ausreichend, da $300'' = 5' = 0\overset{\circ}{.}08333$ immer noch ein recht kleiner Gradwert ist und in einer Grafik praktisch nicht sichtbar wäre.