nathan frad

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

Recommended Posts

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

Share this post


Link to post
Share on other sites
Publicité
En vous inscrivant sur Astrosurf,
ce type d'annonce ne sera plus affiché.
Photographier la Lune
Guide complet pour la photographier de la Lune.
Information et commande sur www.photographierlalune.com

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 ?

 

Share this post


Link to post
Share on other sites

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 )

Edited by nathan frad
Apres de longue recherche

Share this post


Link to post
Share on other sites

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 : 20171130_105554.thumb.jpg.e27787d12f63495523930527df0e2e46.jpg

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...

Share this post


Link to post
Share on other sites

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.

Edited by Bruno-
  • Like 1

Share this post


Link to post
Share on other sites

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 ! 

Share this post


Link to post
Share on other sites

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.

 

 

Share this post


Link to post
Share on other sites

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

 

 

Share this post


Link to post
Share on other sites

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...

Edited by Bruno-

Share this post


Link to post
Share on other sites

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 .

Share this post


Link to post
Share on other sites

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

Share this post


Link to post
Share on other sites

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°

Share this post


Link to post
Share on other sites

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 ?

 

Edited by Bruno-
  • Confused 1

Share this post


Link to post
Share on other sites

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.

Edited by nathan frad

Share this post


Link to post
Share on other sites

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...

 

 

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now