welche Genauigkeit benötigst du. Hier eine einfache Methode. Falls du aber eine höhere Genauigkeit benötigst muss ich noch mal zu Hause scheuen.
Gruß
Lots
Ich benötige eine Formel zur möglichst einfachen Berechnung der Rektaszension und Deklination der Sonne.
Es sei frac(x) die Funktion, die die Nachkommastellen einer Zahl zurueck gibt: also frac(7.21) = 0.21 aber frac(-0.1) = 0.9, frac(-1.2) = 0.8; Es sei T die Anzahl Jahrhunderte seit dem 1.1.2000, ist also im Moment noch negativ.
Naeherung: T = 0.01*(Jahreszahl - 2000 + (Tagnummer-1)/365.25)
Tagnummer: 1.1 = #1, 2.1 = #2 , u.s.w
Wenn es genauer sein muss, muss man das Julianische Datum (JD) gebrauchen: T = (JD - 2451545.0)/36525.0
P2 = 6.28318.. ( = 2*pi)
M = P2 * frac(0.993133+99.97361*T)
DL = 6893.0 * sin(M) + 72.0*sin(M+M)
L = P2 * frac(0.7859453 + M/P2 + (6191.2*T +DL)/1.296e+6)
X = cos(L)
Y = 0.91748*sin(L)
Z = 0.39778*sin(L)
r = sqrt(1.0 - Z*Z)
Deklination = (360/P2)*arctan(Z/r), in Gradmass
Rektaszension = (48.0/P2)*arctan(Y/(X+r)),
in Stunden, falls < 0 muss man 24 Stunden hinzu addieren.
Wenn Sie die Naeherung fuer T gebrauchen, ist die Position der Sonne auf ein Grad genau, mit dem Julianischen Datum erreichen Sie eine Bogenminute, also 60 x genauer. Die errechneten Koordinaten beziehen sich auf das Aequinoktikum des Datums, also das was man braucht, um z.B. die Hoehe ueber dem Horizont zu berechnen.
Literatur:
Oliver Montenbruck, Thomas Pfleger
Astronomie mit dem Pesonal Computer
schau auch mal
hier
