HELP : éphémérides des planètes

[nathan frad 1]

Je viens de me rendre contre que dans la théorie vsop87, il n'y a pas la Lune...

Du coup je vais calculer la position de la lune avec la méthode décrite dans le livre de Jean MEEUS. Dans le livre, il utilise la théorie ELP2000 de M.Chapront et H.Chapront. Ce pendant, pour l'instant, je ne retrouve pas sur internet, les termes de la théorie ELP2000 présente dans le livre ...

Exemple de source trouvé : ftp://cdsarc.u-strasbg.fr/pub/cats/VI/79/

ftp://ftp.imcce.fr/pub/ephem/moon/elp82b/

http://astronomie.baillet.org/logiciels/heclium/telechargement.php

.... Si vous avez des informations sur le sujet n’hésiter pas.

Merci d'avance

[Bruno- 3 965]

Grâce à toi, j'ai ré-écrit un vieux programme d'éphémérides planétaires qui datait de la décennie précédente (je n'avais jamais utilisé les éléments précis de la théorie VSOP87) mais... zut, et la Lune ? Je me suis posé la même question, et je n'ai rien trouvé de mieux que l'algorithme de Meeus.

 

Hé, mais le lien FTP du C.D.S  semble contenir tout ce qu'il faut, il y a même un programme Fortran. J'ai juste regardé vite fait, il y a tout ce qu'il faut, non ? Du coup quel est le problème ?

 

[nathan frad 1]

Ah super ! J'ai trouvé une autre source : http://cdsarc.u-strasbg.fr/viz-bin/Cat?cat=VI%2F79&target=http&

Dans celui ci il y a un fichier exemple tout en bas . Et il y aussi un fichier nommé elp82b.ps.gz ( vers le bas ) ( je te le met aussi en fichier joint),  c’est un pdf explicatif mais j'avoue que je ne comprend pas bien tout ^^

elp82b.ps.gz

 

Modification :

J'ai trouvé le calcul du livre de Meeus sur internet. Les terme sont presque identique. ( http://www.stargazing.net/kepler/ch30.txt )

Modifié 30 novembre 2017 par nathan frad
Apres de longue recherche

[nathan frad 1]

Enfaîte mon problème , c'est que dans le livre pour le calcul de la position de la lune il se limite aux termes periodique les plus important de la théorie ELP2000 pour calculer sa longitude et latitude . Pour la longitude, il fait la somme des termes cela donne ça :

Et moi du coup je voudrais faire la somme de la longitude ( par example )  mais avec tout les termes et non que les plus important. Et dans la theorie elp200 il y a plusieurs fichier contenant des termes longitude :

ELP1 :  Main problem. Longitude periodic terms (sine)

ELP4 :  Earth figure perturbations. Longitude

ELP7 : Earth figure perturbations. Longitude/t

ELP10 : Planetary perturbations. Table 1 Longitude

ELP13 : Planetary perturbations. Table 1 Longitude/t

ELP16 : Planetary perturbations. Table 2 Longitude

ELP19 : Planetary perturbations. Table 2 Longitude/t

ELP22 : Tidal effects. Longitude

ELP25 : Tidal effects. Longitude/t

ELP28 :  Moon figure perturbations. Longitude

ELP31 :  Relativistic perturbations. Longitude

ELP34 :  Planetary perturbations - solar eccentricity. Longitude/t2

 

Il y a trop de fichier  comportant des " longitude" , du coup je ne sais pas le qu’elle prendre pour faire la somme de tout les termes qui sont dans ce projet...

[Bruno- 3 965]

Je regarderai ça quand j'aurais le temps (ça m'intéresse aussi). Ce n'est pas urgent ?

 

---

Re. J'ai regardé un peu le fichier “elp82b.ps” mais je n'ai pas compris comment utiliser les données des fichiers EPL*. Par exemple dès les premiers fichiers EPL1, EPL2, EPL3, je comprends bien qu'il faut calculer une série, mais je ne sais pas à quoi servent B,1, B2, B3, B4, B5 et B6, ni où calculer D, l, l' et F. Peut-être qu'il faut se coltiner la théorie complète... Mais je n'abandonne pas : j'essaierai de décortiquer le programme Fortran.

Modifié 30 novembre 2017 par Bruno-

[nathan frad 1]

Salut, Salut ! Bon bah me voila de retour . 

Après avoir testé pour chaque planètes leurs positions, je me rend compte que  mon calcul me sort les bonnes positions ( RA , DEC ) que  pour Saturne et Venus ^^... Bizarre le calcul étant pourtant universel " normalement" ... Je n'ai plus qu'a chercher l'erreur dans tout les calculs... Let's go ! 

[Bruno- 3 965]

Cherche encore les problèmes d'atan2, peut-être qu'il y a encore des endroits où tu as utilisé un arc tangente ou un arc sinus au lieu de la fonction atan2. Ça se trouve Vénus et Saturne sont dans des portions d'orbite où le problème de tangente n'intervient pas (genre : leur longitude est inférieure à 180°) contrairement aux autres planètes.

 

Sinon, tu peux indiquer les résultats des calculs intermédiaires, je pourrai comparer avec mes calculs et ainsi on pourra identifier l'endroit où il y a une erreur.

 

 

[nathan frad 1]

J'ai commencé a chercher les problème d'atan2, mais je continuerais demain. C'est bizarre car les longitudes de venus et saturne sont supérieur a 180°. ( venus longitude en deg : 233.40442810862078 ,  saturn longitude en deg : 269.7095372598329 ) ...

 

C'est très gentil de proposer de comparer nos résultat : voici mes résultat pour mars : Un grand merci d'avance !!

 

Debut constructeur greenwichPosition --   
Julian day : 2458092.7674421296
t : 0.01792681024539234

 heliocentrique L (deg) : 182.02328635700997
 heliocentrique B  (deg): 1.3677899490808634
 heliocentrique R : 1.6493478011978169

 heliocentrique  terre Lo (deg) : 73.31189105773137
 heliocentrique terre Bo  (deg): 2.2184043058413104E-5
 heliocentrique terre Ro : 0.9854272119530312
x : -1.9308268537879374
y : -1.0021378363422846
z : 0.03936983200723439
Vrai distance : 2.175757929092687
temps lumiere : 0.012566129735844917

Coordonnées heliocentrique de venus pour l'instant t-T en degres

 heliocentrique L1 (deg  pour t-T): 182.01769026187776
 heliocentrique B1 (deg  pour t-T): 1.3679115509958741
 heliocentrique R1  : 0.7242975917764743

x1 : -1.0066192270923175
y2 : -0.969416945334475
z3 : 0.01729026238490708
Vrai distance1 : 1.3976231378404187
temps lumiere1 : 0.00807199800910076

longitude en deg: 223.92143576175135
latitude en deg: 0.7088350970353545
l'exentricité de le l'orbite de la terre (e) : 0.016701094035012586
longitude du perihelie ( pi deg ) : 103.24561090717062
piRad : 1.8019758485630968

long geo moyenne du soleil ( LO) : 6734.256164048508
anomalie moyenne du soleil ( M) : 6811.010540692197
centre du soleil : 0.07360506468399858
LongVraiSoleil : 254.32976911319201
LongVraiSoleilRad : 4.438891856862179
K en radian: 9.936508497454117E-5
Aberration_Longitude_deg : -0.026625793469303344
Aberration_Latitude_deg : -0.014763885925450427

nutation_longitude en deg  : -0.0034513864105610603
nutation_latitude : -7.477089645817657
Longitude_nutation :  223.91798437534078
rad en deg : 221.67866015827406
rad en heure : 14.785966632556878
dec en deg : -15.339496520557677
Fin constructeur greenwichPosition -- *******************************************   
 pi : 2.5849317577312025E-5
TSM :170.52953186631203

long geo moyenne du soleil ( L0) : 1.5764727194520578E9
anomalie moyenne du soleil ( M) : 118.87455932311805

nutation_longitude en deg  : -0.0030281514152984018
nutation_latitude : -6.507131460476364
ObliquitM : 23.003626543129922
ObliquitVrai : 16.49649508265356
TSM : 170.51844129819844
TSM heure : 11.367896086546562
longitude : -0.5805
longitude heure : -0.0387
H_rad : -0.9030459014570468
u degres : 44.744230177359114
 
ra_position : 221.67951667644542
ra_position en heure : 14.786023762318909
dec_position : -15.33864000238631

 

 

[Bruno- 3 965]

Voilà ce que j'obtiens avec cette valeur du jour julien (ce sont des coordonnées moyennes, je n'ai pas tenu compte de la nutation et de l'aberration) :

Coordonnées cartésiennes héliocentriques :
    planète 3 : X0 = 0.282964   Y0 = 0.943927   Z0 = 0.000000
    planète 4 : Xh = -1.647849  Yh = -0.058225  Zh = 0.039370
Coordonnées écliptiques héliocentriques :
    planète 4 : λ = 182.023638° β = 1.367782°   ρ = 1.649347
Coordonnées cartésiennes relatives 4 / 3 :
    planète 4 : X = -1.930813   Y = -1.002152   Z = 0.039370
Coordonnées écliptiques relatives 4 / 3 :
    planète 4 : l = 207.430714° b = 1.036807°   r = 2.175752
Coordonnées équatoriales moyennes pour l'équinoxe :
    α = 205.842131°     δ = -9.589192°

(planète 3, c'est la Terre, et planète 4, c'est Mars)

Reste à comparer...

 

Ah, j'ai relevé un truc bizarre dans ton message précédent :

Citation

Coordonnées heliocentrique de venus pour l'instant t-T en degres

 heliocentrique L1 (deg  pour t-T): 182.01769026187776
 heliocentrique B1 (deg  pour t-T): 1.3679115509958741
 heliocentrique R1  : 0.7242975917764743

 

 

C'est Vénus ou Mars ? Du coup R1 devient celui de Vénus, alors qu'au début c'était celui de Mars (et j'avais le même : 1,649 et quelques). Par contre les longitude et latitude sont toujours celles de Mars (182 et quelques et 1,36 et quelques). Tu devrais revoir ça...

 

Je n'arrive pas à comparer la suite car tu ne dis pas explicitement ce que sont X1, Y2 et Z3, puis « longitude en deg », « latitude en deg », c'est un peu le fouillis...

Modifié 5 décembre 2017 par Bruno-

[nathan frad 1]

Merci , je remet de l'ordre dans mes résultat afin de pouvoir comparer mes résultat avec les tiens et je reviens vers toi en modifiant ce message.

 

C'est bien Mars et non venus, juste une erreur d'affichage .

[Bruno- 3 965]

Et ça aussi c'est une erreur d'affichage : « heliocentrique R1  : 0.7242975917764743 » ? On dirait la distance de Vénus, justement.

[nathan frad 1]

Je viens de tout afficher correctement et cela est correcte pour toute les planètes, c’était un problème de variable non mis a jours. Merci beaucoup !

Exemple pour mars aujourd’hui :

Julian day : 2458093.5606712964
t : 0.01792898198849127
Coordonnées écliptiques héliocentriques :
    Planète Mars  : λ = 182.376633438691°  β = 1.3600854103658517°  ρ = 1.6488959513779677°
    Planète Terre :  λo = 74.11666839850386°  βo = 359.9999974415216°  ρo = 0.9853179018538474°
Coordonnées cartésiennes héliocentriques :
    Planète Mars :  X = -1.9166746951200193  Y = -1.016056876472827  Z = 0.03913780831125706
Coordonnées écliptiques héliocentriques pour l'instant t-T (aberration) :
    Planète Mars  : λ1 = 182.3710499309832°  β1 = 1.360207558857363°  ρ1 = 1.6489031342337546°
Coordonnées cartésiennes héliocentriques pour l'instant t-T (aberration)  :
    Planète Mars :  X1 = -1.9166884400884336  Y1 = -1.0158966673340497  Z1 = 0.03914149309682227
Coordonnées écliptiques héliocentriques pour l'instant t-T  :
    Planète Mars  : λ1 = 207.9247984504529°  β1 = 1.033710442737094°  ρ1 = 1.6489031342337546°
Coordonnées écliptiques héliocentriques pour l'instant t-T + abberation  :
    Planète Mars  : λ1 = 1.009780488577874°  β1 = 1.0085960442814097°  ρ1 = 1.6489031342337546°
Coordonnées écliptiques héliocentriques pour l'instant t-T + abberation + nutation  :
    Planète Mars  : λ1 = 207.92138108076188°  β1 = 1.0085960442814097°  ρ1 = 1.6489031342337546°
Coordonnées équatoriales moyennes pour l'équinoxe : :
    Planète Mars  : α = 206.29608276909988°  α heure = 13.759948720698961 h  δ = -9.792715014994402°

[Bruno- 3 965]

OK. Donc il ne reste plus qu'à s'occuper de la Lune, je regarderai ça un de ces jours parce que ça m'intéresse.

 

Une question : quand tu parles d'instant t-T, T est la durée qu'il faut à la lumière pour faire le trajet entre la planète et la Terre, c'est bien ça ?

 

Modifié 6 décembre 2017 par Bruno-

[nathan frad 1]

Oui tout a fait, c'est l'effet du Temps-lumière.

Et pour la Lune j'ai tout les calculs, il me manque juste la somme des termes de la théorie ELP2000 . Merci d'avance .

 

New:

Je suis en entreprise jusqu'au 6 janvier (avant de retourner en formation)  donc je vais continuer a rechercher des nouvelles informations. Et je vais peut être mettre un autre post sur le forum pour la théorie ELP2000.

Modifié 26 décembre 2017 par nathan frad

[nathan frad 1]

Holà tout le monde ! Voici le résultat des éphémérides des planètes ! : https://github.com/nathanfrad/PlanetaryEphemeris

 

[Bruno- 3 965]

Ah, merci d'avoir posté la version finale ! Je regarderai ça avec intérêt vu que je n'ai pas encore réussi à calculer les positions de la Lune...