Re: Gibt es ein "Sonnensystemflash" für die Homepa
Hallo Hendrik,
Wo finde ich denn Formeln zur korekten Berechnung der Planeten Positionen?
Ich würde sowas gerne mal als Java Applet schreiben.
Kommt drauf an, was für eine Genauigkeit du gerne rechnen möchtest. Ich finde für die meisten "kleinen" Sachen (wie z.B. meine Flashs) oder für Applets fürs Netz reicht eine "Truncated Version" (wie z.B. von Meeus) völlig aus. So genau darstellbar, wie die komplette VSOP87 Planetenpositionen rechnet, ist es in einem Applet/Flash ohnehin nicht. Wenn du natürlich irgendwie über Anklicken die Daten textmäßig ausgeben möchtest (z.B. die heliozentrischen Werte oder auch die daraus errechneten RA/DEC), dann kann man freilich auch die ganze VSOP nehmen.
Im Grunde genommen berechnt man - wie du sicherlich weißt - die JDE des gesuchten Zeitpunktes und summiert dann einfach die tabellenartig gegebenen Cosinus-Terme auf.
Quellen z.B.:
http://www.freevbcode.com/ShowCode.asp?ID=464
Download als .zip, die Terme stehen jeweils in den Dateien "
Planet_LBR.bas"
Hier kann man relativ gut die Termreihen "abbrechen", weil die Terme mit absteigenden Koeffizienten sortiert sind.
Oder hier der Quellcode, den Celestia verwendet: (Achtung, gekürzte Version, aber nicht ident mit der von Meeus!)
http://celestia.teyssier.org/source-documentation/vsop87_8cpp-source.html
Nicht so schön ist hier die ungeordnete Anordnung der Terme. Aber man kann ja nicht alles haben....
Die gesamte VSOP87 (Versionen A, B, C, ....) kannst du via ftp z.B. hier beziehen:
ftp://ftp.bdl.fr/pub/ephem/planets/vsop87/
Eine Fehlerabschätzung, wenn man die Reihen abbricht, ergibt sich wie folgt: (Empirisch nach Bretagnon&Francou)
Wenn n die Anzahl der verbliebenen Terme ist und A die Amplitude (Koeffizient) des letzten Terms, dann ist der Fehler ca.
q*Sqrt(n)*A, wobei q<=2 ist.
Viel Spaß beim Programmieren!
greets,
harald
--
