Bonjour à tous
je désire calculer la position du soleil, et parr hazard j'ai trouvé un algorithme sur le net.
quelqu'un peu m'expliquer cet algorithme .
que signifie alpha et delta? et est ce qui'ils sont suffisant pour déterminer la position du soleil??
merci infiniment d'avance
Code:
Position of the Sun at 11:00 UT on 1997 August 7th
1. Find the days before J2000.0 (d) from the table
d = 11/24 + 212 + 7 - 1096.5 = -877.04167
2. Find the Mean Longitude (L) of the Sun
L = 280.461 + 0.9856474 * d
= -583.99284 + 720
(add multiples of 360 to bring in range 0 to 360)
= 136.00716
3. Find the Mean anomaly (g) of the Sun
g = 357.528 + 0.9856003 * d
= -506.88453 + 720
= 213.11547
4. Find the ecliptic longitude (lambda) of the sun
lambda = L + 1.915 * sin(g) + 0.020 * sin(2*g)
= 134.97925
(note that the sin(g) and sin(2*g) terms constitute an
approximation to the 'equation of centre' for the orbit
of the Sun)
beta = 0 (by definition as the Sun's orbit defines the
ecliptic plane. This results in a simplification
of the formulas below)
5. Find the obliquity of the ecliptic plane (epsilon)
epsilon = 23.439 - 0.0000004 * d
= 23.439351
6. Find the Right Ascension (alpha) and Declination (delta) of
the Sun
Y = cos(epsilon) * sin(lambda)
X = cos(lambda)
a = arctan(Y/X)
If X < 0 then alpha = a + 180
If Y < 0 and X > 0 then alpha = a + 360
else alpha = a
Y = 0.6489924
X = -0.7068507
a = -42.556485
alpha = -42.556485 + 180 = 137.44352 (degrees)
delta = arcsin(sin(epsilon)*sin(lambda))
= 16.342193 degrees
Final result
Right ascension is usually given in hours of time, and both
figures need to be rounded to a sensible number of decimal
places.
alpha = 9.163 hrs or 9h 09m 46s
delta = +16.34 degrees or +16d 20' 32"
The Interactive Computer Ephemeris gives
alpha = 9h 09m 45.347s and
delta = +16d 20' 30.89"
Pour t'expliquer tu dois avoir des notions d'astronomie (anomalie vraie, anomalie moyenne, équation du centre, longitude du périgée, obliquité, excentricité, nutation, aberration, temps de propagation, réfraction, équinoxe moyen...) , je te conseille vivement le livre Astronomie Générale (André Danjan).
sinon alpha : c'est l’ascension droite.
et delta: la déclinaison.
Oui c'est suffisant c'est les coordonnées équatoriales du soleil.
Ciao
10/04/2012 - 23h11
raikko21
Date d'inscription
juin 2011
Âge
32
Messages
193
Re : Position du soleil
L'algorithme que tu as donné fait les opérations suivantes:
1- calcul du nombre de jour julien écoulé depuis j2000 (1/1/2000). Dans la 1ere equation 11/24 correspond a l'heure (11:00) convertie en fraction de jour, sinon les autres parametres (année, mois, jour) ont été remplacés dans l'équation;
2- Calcul de longitude moyenne du soleil;
3- Calcul de l'anomalie moyenne du soleil;
4- Calcul de la longitude apparente du soleil;
5- Calcul de l'obliquité;
6- Calcul de l'ascension droite et de la déclinaison (coordonnées équatoriales desquelles tu peux déduire l'azimut et la hauteur).
Par contre, je n'aime pas trop cet algorithme et les equations choisies je te conseille Astronomical ALgorithms de Jean Meeus.
10/04/2012 - 23h25
narakphysics
Date d'inscription
février 2011
Messages
813
Re : Position du soleil
Merci beaucoup pour votre réponse
en ce qui concerne le livre Astronomical ALgorithms de Jean Meeus, je doit lire le chapitre Equation of Kepler?
11/04/2012 - 08h35
raikko21
Date d'inscription
juin 2011
Âge
32
Messages
193
Re : Position du soleil
Bonjour,
Ben d'abord tu jettes un coup d’œil aux chapitres Julian Day, Nutation and the Obliquity of th Ecliptic, et Solar Coordinates, et tu dois savoir par exemple comment lire les table des éléments de l'orbite.
Par exemple, pour le chapitre Solar Coordinates la première équation donne la longitude moyenne du soleil rapportée a équinoxe moyen de la date.
Si tu prends la table 30.A (page 200) "Orbital Elements for the mean equinox of the date", dans les paramètres de la Terre tu retrouve les coefficients utilisés.
A la première ligne il y a marquer "L" tu retrouve les coefficients de l’équation de "" (au premier paramètre tu dois rajouter 180 degrés).
On peut voir par exemple que le coefficient de est si tu le rajoute tu aura l’équation: (Tout dépend de la précision que tu souhaites avoir).
Pour obtenir les coefficients de l'anomalie moyenne, toujours dans la même table tu soustrais la dernière ligne a la première ligne et tu auras tes coefficients...
Voila, juste une petite précision T est en siècle juliens.