Unter Verwendung von Gleichung 28 in Abschnitt 3.8.1 von Solar Position Algorithm (Reda, I. und A. Andreas, Jan. 2008) habe ich die mittlere Sternzeit in Greenwich für den 17. Oktober 2003 für vier verschiedene Zeiten (UT) mit einem Casio berechnet fx-991MS. Die berechneten Werte wurden online mit der URL- Berechnung von GAST, GMST und ERA gegengeprüft
Gleichung 28 zur Berechnung von GMST lautet:
v0=280.46061837 + 360.98564736629*(JD - 2451545) + 0.000387933*JC^2 - JC^3/38710000
JC wird berechnet als:
JC = (JD - 2451545)/36525
Die berechneten Werte von JD und JC lauten wie folgt:
1. 00:00:00 UT: 2452929.5 (JD) & 0.037905544 (JC)
2. 05:30:00 UT: 2452929.729 (JD) & 0.037911818 (JC)
3. 11:00:00 UT: 2452929.958 (JD) & 0.037918083 (JC)
4. 19:21:00 UT: 2452930.306 (JD) & 0.037927611 (JC)
Es gibt jedoch einen deutlichen Versatz zwischen den berechneten Werten und den Werten, die von der oben genannten URL bereitgestellt werden. Zum Beispiel:
1. 00:00:00 UT: 22.32738 hour (Computed) & 1.6726265 hour (URL)
2. 05:30:00 UT: 16.7231 hour (Computed) & 7.1960077 hour (URL)
3. 11:00:00 UT: 11.30528667 hour (Computed) & 12.702743 hour (URL)
4. 19:21:00 UT: 2.93042 hour (Computed) & 21.075605 hour (URL)
Wenn ich außerdem 22,32738 Stunde von 24,0 subtrahiere, erhalte ich 1,67262, was der von der URL bereitgestellte Wert ist.
Wenn ich zur gleichen Zeit mit demselben Gleichungssystem die GMST für den 10. April 1987 um 00:00:00 berechne, sind die berechneten und die URL-Werte dieselben.
10. April um 00:00:00:
JD=2446895.5
JC=-0.127296372
GMST = 13.18 hour (Computed) & 13.17954 hour (URL)
Bitte helfen Sie mir zu verstehen, wo genau das Problem liegt. Wenn das Problem bei der Interpretation liegt, helfen Sie bitte zum besseren Verständnis.
Hinweis: Knew gmst ist jetzt veraltet, aber da ich ein Anfänger bin, dachte ich in erster Linie daran, gmst zu verstehen, und werde mich dann mit dem aktuellen Szenario befassen.
Vielen Dank.
Die Formel 28 von SPA, die Sie verlinkt haben, ist mit Einheiten und in einer Form angegeben, die nicht mit offiziellen Versionen übereinstimmt. Es erfordert ein wenig Arbeit, es zu „übersetzen“: Außerdem ist es unnötig anfällig für Rundungsfehler (weil alle signifikanten Zahlen für Datum und Uhrzeit in einer sehr langen Zahl aneinandergereiht sind). Die Autoren von SPA scheinen nicht nach vollständiger offizieller Genauigkeit gestrebt zu haben, wahrscheinlich erforderte ihre Anwendung dies nicht.
Ich kann aus den Daten in der Frage nicht ganz ableiten, was mit den Berechnungen passiert sein könnte, online oder auf andere Weise, aber vielleicht könnten die folgenden Referenzdaten hilfreich sein.
GMST für 0h UT ist immer noch im Astronomical Almanac tabelliert, und es gibt auch zwei ziemlich einfache Formeln, die von den offiziellen Definitionen abgeleitet sind, die Sie verwenden können (siehe unten), sie ermöglichen eine etwas bessere Erhaltung der Endzahlen als eine Berechnung im SPA-Formular .
In der Zwischenzeit, nur als Referenz, gab der Astronomical Almanac für 2003 die GMST für 0 Uhr UT am 17. Oktober 2003 als 1 Uhr 40 Minuten 21,4554 Sekunden (1,6726265 Stunden) an. Das ist für 0h UT, und während der folgenden 24h ist die GMST für x h UT leicht zu finden, sie liegt sehr nahe bei .
2003 war ein Umstellungsjahr für die offiziellen (IAU) Berechnungen, die sich auf GMST und UT1 beziehen. Vor 2003 wurde die Formel von Aoki et al. (1982), „The new definition of universal time“ ( http://adsabs.harvard.edu/abs/1982A&A...105..359A ), angenommen. Ab 1. Januar 2003 gilt die Formel von N. Capitaine et al. (2003), "Expressions to implement the IAU 2000 definition of UT1", Astronomy & Astrophysics 406, 1135-1149 ( http://adsabs.harvard.edu/abs/2003A%26A...406.1135C ) war (und ist immer noch ist) angenommen. Der Unterschied zwischen den Ergebnissen im Jahr 2003 war gering (+0,0003 Zeitsekunden, um das Ergebnis der neueren Berechnung zu erhalten), aber es sieht so aus, als ob AA für 2003 immer noch die ältere Formel von 1982 verwendete.
Die Formel von Aoki et al (1982) ergab:
GMST von 0h UT (in Zeit-Sek.) =
plus eine sehr kleine Korrektur für DeltaT: Diese kann aus der Version der Aoki-Formel entnommen werden, die in der Capitaine-Referenz bei B.3 auf Seite 1149 angegeben ist (die anscheinend diejenige ist, die in offiziellen Berechnungen verwendet wurde), aber sie kann es praktisch bei der Genauigkeit der Daten in der Frage ausgelassen werden. Tcy bedeutet das Zeitintervall von J2000 (julianisches Datum 2451545.0) reduziert auf Jahrhunderte von 36525 Tagen. (Der JD für 0 Stunden am 17. Oktober 2003, 0 Uhr UT war JD(UT) 2452929,5 Tage und DeltaT war nahe bei +64,55 Sekunden.)
Die aktuelle Formel aus der Referenz von Capitaine (2003) (am bequemsten aus Formel B.2 auf Seite 1149) ergibt
plus einen sehr kleinen Term vierter Ordnung (in der Referenz angegeben, aber hier weggelassen) und wiederum eine sehr kleine Korrektur für den leichten Effekt von DeltaT, die weggelassen oder auf Wunsch eingeschlossen werden kann
(Zeit-Sek.), wobei DeltaT immer noch in Sekunden ist.
Reda und Andreas 2008 Gleichung 28 stammt von Jean Meeus, Astronomical Algorithms , Gleichung 12.4. Er leitete es von Aoki et al. 1982 Gleichung 13, die die IAU von 1984 bis 2002 verwendete.
Lassen Sie es uns aufschlüsseln:
Der Satz von Referenzrahmen, den die IAU im Jahr 2000 annahm und 2003 implementierte, wurde durch die extremen Präzisionsanforderungen der Funkinterferometrie mit sehr langer Basislinie (VLBI) motiviert. Der Unterschied zwischen den alten und neuen Ausdrücken für GMST ist am 1.1.2003 null und innerhalb eines Jahrhunderts nach diesem Datum vernachlässigbar, es sei denn, Sie machen selbst VLBI.
Die Ergebnisse, die Sie zuvor für den 17.10.2003 erhalten haben, könnten darauf zurückzuführen sein, dass irgendwo eine Addition oder Subtraktion invertiert wurde, vielleicht bei der Normalisierung auf [0, 360]. Wenn Sie es noch einmal versuchen, wird es möglicherweise so aussehen wie am 10.04.1987.
Schlau
Mike G