Position du soleil

**narakphysics** · 03/04/2012, 22h17

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"

invitefe6f47fa · 10/04/2012, 13h21

Bonjour,

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

invitefe6f47fa · 10/04/2012, 23h11

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.

**narakphysics** · 10/04/2012, 23h25

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?

A voir en vidéo sur Futura · Aujourd'hui

invitefe6f47fa · 11/04/2012, 08h35

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.
$\text{[math]}$

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 " $\text{[math]}$ " (au premier paramètre tu dois rajouter 180 degrés).

On peut voir par exemple que le coefficient de $\text{[math]}$ est $\text{[math]}$ si tu le rajoute tu aura l’équation:
$\text{[math]}$ (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.

Ciao

Position du soleil