Bonjour,
Voilà je créé un programme informatique sur les phases lunaires (Calendrier lunaire)
Je créé mon programme en langage c sous GNU/Linux debian
Il éxiste déjà sous GNU/Linux des programmes sur la lune "glunarclock", "wmoonclock" et
"xplanet"
Vous allez peut-être me prendre pour un fou de créer un programme comme celui-ci.
Et vous allez me dire pourquoi créer un programme vu qu'il y a tout sur internet.
Des centaines de paramètres à prendre en compte... Oups
Je suis tétu et continuerai à vouloir espérer arriver un jour à un programme "fiable"
Dans un premier temps, j'ai créé un programme en utilisant la formule de copernic (Le mois synodique)
Sauf que ms est une moyenne.
Concrètement la lunaison peux varier entre 29j 6h et 29j 20h (+ ou - 7 heures).
Je ne peux donc pas utiliser Copernic.
Ce premier programme me donnait une date décalée d'un jour pour certaines phases lunaires.
J'ai trouvé un algorithme sur le calcul des phases lunaires qui semble bizarre aussi
http://zpag.net/Calendrier/calculer_phases_lune.htm
Vous pouvez vous amuser à vérifier
Pour Mai 2009, il est noté Plaine Lune = 10 Mai (Plaine ???)
Faites le test sur ce puissant site
http://portail.imcce.fr/page.php?nav=fr ... /index.php
Il vous donnera
PLEINE LUNE 9 mai 2009 à 04h01m UTC
Hallucinant non ?
Cet algorithme donne donc une fausse date pour cette plEine lune.
J'ai donc le même genre de bug (un écart d'une journée).
Je me suis décidé à créer un nouveau programme qui intègre le mouvement de la lune.
Merci donc à Newton, Kepler et Euler.
Grâce à l'équation du temps je trouve donc la position de la terre (xT,yT,zT) par rapport au soleil à un instant tn+1.
La terre accélére et ralenti de + ou - 15 minutes pour une rotation autour du soleil donc sur 365.257 jours.
Il n'éxiste pas d'équation du temps pour celle de la lune !
Le mouvement de la lune peux être décrit avec cette équation
C'est beau, non ?
Donc cette équation me donne la vitesse (ou bien l'accélération) en fonction de :
1- La position en coordonnées cartésiennes de la lune par rapport au repère de la terre (xL,yL,zL)
2- La position en coordonnées cartésiennes du soleil par rapport au repère de la terre (xS,yS,zS)
3- La distance Lune-Terre et Terre Soleil (rL, rT)
4- La distance Lune-Soleil (delta)
Est-il possible d'écrire une sorte d'équation du temps de la lune autour de la terre ?
J'utilise cette équation pour calculer la vitesse de la lune et Euler me permet de la calculer à un instant tn+1. C'est déjà un bon début
L'équation du temps me donnait une valeur d'angle tétaT en fonction du temps.
Plutôt intéressant pour calculer la position de la terre autour du soleil.
xT=rT*cos(tétaT)
yT=rT*sin(tétaT)
Pour connaître LA position (xL,yL,zL) il me faut intégrer la rotation de la lune autour de la terre
Donc il me faudrait un angle (tétaL) en fonction du temps.
TétaL=wL*t avec wL vitesse de rotation autour de la terre.
En plus l'orbite de la lune est vraiment complexe.
J'ai découvert en fait que sa trajectoire elliptique est orientée à 5°9' en moyenne.
En plus il éxiste 2 rotations suplémentaires de pertubations.
J'ai comme l'impression de tourner en rond ???? hem heu,... plutôt accéléré sur une trajectoire elliptique
A mon avis l'équation vue plus haut est la clé.
La troisième loi de Kepler peux sûrement m'aider.
J'ai besoin d'y voir plus clair (Comme en phase de pleine lune la nuit ).
Je suis un peu désorienté (Dans la lune quoi )
Pouvez-vous me dire s'il est possible de connaître la position "éxacte" à l'instant tn+1 de la lune J'ai besoin de conaître sa position (xL,yL,zL) à l'instant tn+1
Ou bien est-il possible de visualiser la courbe d'accélération de la lune ?
Je ne sais pas si j'arrive vraiment à me faire comprendre
Merci
damien
Voilà je créé un programme informatique sur les phases lunaires (Calendrier lunaire)
Je créé mon programme en langage c sous GNU/Linux debian
Il éxiste déjà sous GNU/Linux des programmes sur la lune "glunarclock", "wmoonclock" et
"xplanet"
Vous allez peut-être me prendre pour un fou de créer un programme comme celui-ci.
Et vous allez me dire pourquoi créer un programme vu qu'il y a tout sur internet.
Des centaines de paramètres à prendre en compte... Oups
Je suis tétu et continuerai à vouloir espérer arriver un jour à un programme "fiable"
Dans un premier temps, j'ai créé un programme en utilisant la formule de copernic (Le mois synodique)
Code : Tout sélectionner
ms=(1/((1/27,322)-(1/365,257)))
ms=29.530589 jours
Sauf que ms est une moyenne.
Concrètement la lunaison peux varier entre 29j 6h et 29j 20h (+ ou - 7 heures).
Je ne peux donc pas utiliser Copernic.
Ce premier programme me donnait une date décalée d'un jour pour certaines phases lunaires.
J'ai trouvé un algorithme sur le calcul des phases lunaires qui semble bizarre aussi
http://zpag.net/Calendrier/calculer_phases_lune.htm
Vous pouvez vous amuser à vérifier
Pour Mai 2009, il est noté Plaine Lune = 10 Mai (Plaine ???)
Faites le test sur ce puissant site
http://portail.imcce.fr/page.php?nav=fr ... /index.php
Il vous donnera
PLEINE LUNE 9 mai 2009 à 04h01m UTC
Hallucinant non ?
Cet algorithme donne donc une fausse date pour cette plEine lune.
J'ai donc le même genre de bug (un écart d'une journée).
Je me suis décidé à créer un nouveau programme qui intègre le mouvement de la lune.
Merci donc à Newton, Kepler et Euler.
Grâce à l'équation du temps je trouve donc la position de la terre (xT,yT,zT) par rapport au soleil à un instant tn+1.
La terre accélére et ralenti de + ou - 15 minutes pour une rotation autour du soleil donc sur 365.257 jours.
Il n'éxiste pas d'équation du temps pour celle de la lune !
Le mouvement de la lune peux être décrit avec cette équation
Code : Tout sélectionner
/* Le repère de la terre est le repère de réfèrence */
/* en x */
/* d²x/dt²=-G*(MT+ML)*xL/rL³+G*MS*Ax */
/* Avec Ax=((xS-xL)/delta³)-(xS/rT³) */
/* Avec d²x/dt² accélération de la lune */
/* G Constante gravitationnelle en m³.Kg-¹.s-² (Newton) */
/* MT Masse de la terre */
/* ML Masse de la lune */
/* MS Masse du soleil */
/* rL Distance Lune-Terre */
/* rT Distance Terre-Soleil */
/* xL position de la lune par rapport à la terre */
/* xS position du soleil par rapport à la terre */
/* delta Distance Lune-Soleil */
/* Cette équation s'écrit de façon analogue en y et z */
C'est beau, non ?
Donc cette équation me donne la vitesse (ou bien l'accélération) en fonction de :
1- La position en coordonnées cartésiennes de la lune par rapport au repère de la terre (xL,yL,zL)
2- La position en coordonnées cartésiennes du soleil par rapport au repère de la terre (xS,yS,zS)
3- La distance Lune-Terre et Terre Soleil (rL, rT)
4- La distance Lune-Soleil (delta)
Est-il possible d'écrire une sorte d'équation du temps de la lune autour de la terre ?
J'utilise cette équation pour calculer la vitesse de la lune et Euler me permet de la calculer à un instant tn+1. C'est déjà un bon début
L'équation du temps me donnait une valeur d'angle tétaT en fonction du temps.
Plutôt intéressant pour calculer la position de la terre autour du soleil.
xT=rT*cos(tétaT)
yT=rT*sin(tétaT)
Pour connaître LA position (xL,yL,zL) il me faut intégrer la rotation de la lune autour de la terre
Donc il me faudrait un angle (tétaL) en fonction du temps.
TétaL=wL*t avec wL vitesse de rotation autour de la terre.
En plus l'orbite de la lune est vraiment complexe.
J'ai découvert en fait que sa trajectoire elliptique est orientée à 5°9' en moyenne.
En plus il éxiste 2 rotations suplémentaires de pertubations.
J'ai comme l'impression de tourner en rond ???? hem heu,... plutôt accéléré sur une trajectoire elliptique
A mon avis l'équation vue plus haut est la clé.
La troisième loi de Kepler peux sûrement m'aider.
J'ai besoin d'y voir plus clair (Comme en phase de pleine lune la nuit ).
Je suis un peu désorienté (Dans la lune quoi )
Pouvez-vous me dire s'il est possible de connaître la position "éxacte" à l'instant tn+1 de la lune J'ai besoin de conaître sa position (xL,yL,zL) à l'instant tn+1
Ou bien est-il possible de visualiser la courbe d'accélération de la lune ?
Je ne sais pas si j'arrive vraiment à me faire comprendre
Merci
damien
Dernière modification par damien le mardi 5 janvier 2010 à 18:24, modifié 1 fois.