calcul de l'équation de temps et la déclinaison solaire

[faria]

salut,

je suis dans le domaine informatique et je viens de developper un webservice fournissant les temps des priéres .j'ai trouvé que les temps de salat se calculent en fonction de deux mesures astronomique qui sont l'équation de temps et la déclinaison solaire.

mon probleme se trouve au niveau des calculs j'ai les formules de calcul mais le resultat est toujours erroné pourquoi je ne sais pas .

je cherche à savoir comment appliquer ces formules pour obtenir le bon resultat.

merci de m'éclércis

[Sade]

Bonsoir Faria, bienvenue parmi nous

 

En gros, il s'agit de déterminer l'heure à laquelle le soleil sera à un point précis dans le ciel, si j'ai bien compris. C'est probablement dépendant des saisons, du lieu, etc... Tu es bien sûre d'avoir pris en compte tous les paramètres?

 

Regarde ici:

 

http://fr.wikipedia.org/wiki/%C3%89quation_du_temps

 

Et en cherchant sur Google, on trouve ceci:

 

http://pagesperso-orange.fr/jean-paul.cornec/fichiers_excel.htm

 

Il s'agit probablement d'une erreur de calcul (ou de méthode) de ta part.

 

Tu ne veux pas nous exposer tes calculs, ce serait plus simple?

[faria]

Bonsoir Faria, bienvenue parmi nous

 

En gros, il s'agit de déterminer l'heure à laquelle le soleil sera à un point précis dans le ciel, si j'ai bien compris. C'est probablement dépendant des saisons, du lieu, etc... Tu es bien sûre d'avoir pris en compte tous les paramètres?

 

Regarde ici:

 

http://fr.wikipedia.org/wiki/%C3%89quation_du_temps

 

Et en cherchant sur Google, on trouve ceci:

 

http://pagesperso-orange.fr/jean-paul.cornec/fichiers_excel.htm

 

Il s'agit probablement d'une erreur de calcul (ou de méthode) de ta part.

 

Tu ne veux pas nous exposer tes calculs, ce serait plus simple?

 

re,

merci énormément de m'avoir repondu

heben j'ai pris les valeurs suivantes:

JD = 2454538.5 ( le jour julien correspendant à la date 12/08/2008)

j c'est le jour julien modofie qui se calcule en fct de JD.

pour le jd=2454538.5 ->j=54537

 

et j'ai appliqué les formules suivantes:

 

1ère étape : M = 357° + 0°,9856 x j

M = 54108.6672

 

2ème étape : C = 1°,914 x sin(M) + 0°,02 x sin(2M)

C= 1.80118215361

 

3ème étape : L = 280° + C + 0°,9856 x j

L= 54033.46838215362

 

4ème étape : R = -2°,465 x sin(2L) + 0°,053 x sin(4L)

R = -2.229773894846

 

Et enfin : Equation du Temps (minutes) = (C + R) x 4

Equation du Temps (minutes) = -0.557443473711 et c pas juste du tout .par ce que l'equation de temps dans ce jour réelement égale à 9.5044

 

B - Pour la déclinaison du Soleil le jour "j" ce n'est pas plus compliqué :

 

Sin(déclinaison) = 0,3978 x sin(L)

 

(Ici 0,3978 représente le sinus de l'obliquité de l'écliptique)

 

La combinaison de la variation de ces deux quantités vous permettront alors de tracer de jolies courbes en huit.

 

Notes :

j représente le rang du jour dans l'année (1er janvier = 1)

M est l'anomalie moyenne en degrés

C est l'équation du centre (influence de l'ellipticité de l'orbite terrestre) en degrés

L est la longitude vraie du Soleil en degrés

R est la réduction à l'équateur (influence de l'inclinaison de l'axe terrestre) en degrés

s'il vous plai aider moi je suis perdue

et merci encore une fois;)

['Bruno]

j représente le jour julien modifié (comme annoncé au début de ton message) ou le rang du jour dans l'années (comme annoncé à la fin) ?

 

D'après mon logiciel, le 12/08/2008 correspond à JJ = 2.454.690,5. Apparamment il y a une erreur dans le calcul du jour julien.

 

Si j'applique les formules, ça donne -7,342..., là encore ça ne colle pas.

 

Voyons ce que me dit le livre "Calculs astronomiques amateurs" (Masson). Il donne une méthode approximative (précision annoncée : 50s), ne tenant pas compte de l'équation du centre et de la réduction à l'équateur.

 

- Nombre de jours écoulés depuis le 0 janvier 1901 0h TU :

N = JJ - 2.415.384,5

Ici : N = 39.306

- Longitude moyenne du Soleil :

L = 278,965° + 0,985.647.342° x N (modulo 360°)

Ici : L = 140,819.424.7°

- Longitude du périhélie :

w = 281,235° + 0,000.046.9° x N (modulo 360°)

Ici : w = 283,193.451.4°

- Anomalie moyenne du Soleil :

M = L - w (modulo 360°)

Ici : M = 217,625.973.3°

- Équation du temps :

E = 460s x sin M - 592s sin (2 L)

Ici : E = 299s (j'ai arrondi) soit 4m59s.

 

Ah, et tu dis qu'elle doit valoir 9,5044 minutes ? Ou alors JJ = 2.454.538,5 est la valeur correcte (mais ne correspond pas au 12/08 : c'est le 13/03/2008 à 0h TU. Ah ben oui, c'est ça : c'était hier ! (attention aux fautes de frappes : si tu en fais dans tes messages, il en restera dans le programme, et il est moins intelligent que moi )

 

On recommence pour le 13 mars :

- N = 39.154

- L = 351,001.028.7° (dans tes calculs, tu obtiens 33,468°... !)(1)

- w = 283,071.322.6°

- M = 67,929.706.1° (dans tes calculs, tu obtiens 108,667.2°)

- E = 609s soit 10m09s.

 

C'est déjà plus proche des 9,5 minutes que tu annonces (et qui sont confirmées par un graphique que j'ai trouvé). Étant donné que c'est sensé être une valeur approchée à 50s près, je pense que je suis dans le vrai. Donc que les valeurs de L et M devraient ressembler aux miennes.

 

(Pour avoir plus de précision, il faut calculer avec une plus grande précision L et M, ce qui se fait avec des formules à rallonge - développements limités avec plein, plein de termes...)

 

Je soupçonne donc que tu as oublié de valider le programme en commençant par vérifier qu'il donne les calculs intermédiaires corrects (longitude moyenne et anomalie moyenne par exemple).

 

-------------

(1) La longitude moyenne vaut 0° à l'équinoxe de printemps. Donc, quelques jours avant, elle vaut effectivement dans les 350° plutôt que les 33° que tu obtiens (ça, c'est pour le 24 avril environ).

[faria]

j représente le jour julien modifié (comme annoncé au début de ton message) ou le rang du jour dans l'années (comme annoncé à la fin) ?

 

D'après mon logiciel' date=' le 12/08/2008 correspond à JJ = 2.454.690,5. Apparamment il y a une erreur dans le calcul du jour julien.

 

Si j'applique les formules, ça donne -7,342..., là encore ça ne colle pas.

 

Voyons ce que me dit le livre "Calculs astronomiques amateurs" (Masson). Il donne une méthode approximative (précision annoncée : 50s), ne tenant pas compte de l'équation du centre et de la réduction à l'équateur.

 

- Nombre de jours écoulés depuis le 0 janvier 1901 0h TU :

N = JJ - 2.415.384,5

Ici : N = 39.306

- Longitude moyenne du Soleil :

L = 278,965° + 0,985.647.342° x N (modulo 360°)

Ici : L = 140,819.424.7°

- Longitude du périhélie :

w = 281,235° + 0,000.046.9° x N (modulo 360°)

Ici : w = 283,193.451.4°

- Anomalie moyenne du Soleil :

M = L - w (modulo 360°)

Ici : M = 217,625.973.3°

- Équation du temps :

E = 460s x sin M - 592s sin (2 L)

Ici : E = 299s (j'ai arrondi) soit 4m59s.

 

[faria]

merci beacoup de m'avoir éclairer la méthode d'application des calculs astronomiques. j'ai trouvé l'erreur que j'avai commise a travares la solution que vous m'avez proposé.

merci énormement encore une fois

['Bruno]

Par curiosité, c'était quoi l'erreur ? (c'est pour savoir si ton algorithme est plus précis que le mien, ce qui est probablement le cas).

[faria]

Par curiosité' date=' c'était quoi l'erreur ? (c'est pour savoir si ton algorithme est plus précis que le mien, ce qui est probablement le cas).[/quote']

 

salut,

bon j'ai appliqué mon algorithme pour la méme date 12/03/2008.

jd=2454537.5.

j=jd-2451545 ----->j=2992.5.(l'erreur que j'ai commise)

 

avant j'ai pri la valeur j=54537 parceque j'ai pensé que le jour julien modifié c'est le 'j'.

 

1ère étape :

M = 357° + 0°,9856 x j

M = 3306.408

M mod [360]= 66.408

 

2ème étape :

C = 1°,914 x sin(M) + 0°,02 x sin(2M)

C= 1.76869604

 

3ème étape :

L = 280° + C + 0°,9856 x j

L = 3231.176696

L mod[360]=351.176696

 

4ème étape :

R = -2°,465 x sin(2L) + 0°,053 x sin(4L)

R = 0.716631523

 

Et enfin :

Equation du Temps (minutes) = (C + R) x 4

Equation du Temps (minutes) = 9.941310252

cet valeur correspond à l'equation de temps du jour 12/03/2008 à 00:00 UT et la valeur que j'ai supposé juste (je pense que c'été 9.5044)c'été à 12h a peu pré .donc je peux dire que j ai obtenu la valeur juste .

[lamarchat]

Soyez bien venu dans notre forum ; moi j'ai appris beaucoup de chose ici ; chez nous c'est le tiers monde ; une fois et dans un autre chemin et pour déterminer la vrai trace de l'équateur sur cette terre ; j'essaye de vous dire mais entre nous que nos membre vont utiliser des instruments sophistiqués et des grandes calculs encore plus que ça et malgré ça c'est moi qui va gagner car moi je vais utiliser tout simplement deux petits bidons l'un en haut et l'autre attaché en bas ; je met 4 litres d'eau dans le premier et je laisse l'eau du bidon se vider en rotation sous l'effet de Coriolis force et je récupère cette eau par un tasse en parallèle et ainsi de suite et je m'avance et je m'approche en direction de l'équateur; lorsque la rotation de l'eau s'arrête de tourner je vais mettre mon premier piquet "made in Algeria" qui indique l'équateur et et je vais determiner la grosseur de de cet équateur dés que l'eau va tourner en sens inverse et je vais maitre mon deuxieme piquet ; la distance qui sépare mes deux piquets je vais l'appeler distance gh qui va cacher un incompris en sens de notre physique

 

 

J'ai essayé par ça pour te dire Soyez bien venu parmi nous

[faria]

merci à vous

['Bruno]

OK, merci pour les précisions. Si un jour j'ai besoin de calculer l'équation du temps, je pourrais toujours utiliser cet algorithme : il marche !

[faria]

OK' date=' merci pour les précisions. Si un jour j'ai besoin de calculer l'équation du temps, je pourrais toujours utiliser cet algorithme : il marche ! [/quote']

 

merci à vous:p

[Astro-Pépito]

Edit : c'est bon, j'ai trouvé : il suffit de changer en 365,25 comme prévu et de bien passer en radians : j'avais fait *180/Pi alors que c'est *Pi/180...

 

 

Salut les amis,

 

 

J'aurais besoin de la formule pour calculer la déclinaison pour mes TIPE. J'en ai trouvé une dans le Bulletin de l'Union des Physiciens (BUP), mais la formule ne marche pas (5 annulations au lieu de 2). J'ai pensé à passer en degrés, en radians, mais rien ne s'améliore...

 

Voici la formule :

d=23,43*sin((360*N)/(365*25))

avec N le nombre de jours depuis l'équinoxe de printemps.

Le 365*25 me paraît bizarre, j'ai essayé avec 365,25 qui est plus naturel (vous savez tous pourquoi je suppose), mais là encore ça plante.

 

Pourriez-vous m'aider ?

Merci d'avance

['Bruno]

(Oups, je n'avais pas lu le "Edit".)

[MF_Erwan]

Réfléchisson: le Soleil semble effectuer 360° en un an.

Donc en un jour: 360/365,25=1° environ

Donc en N jours: N*360/365,25 degrés

Conclusion, la déclinaison du Soleil sera.

L'obliquité de l'écliptique est de 23,43°.

Donc la déclinaison sera :

d=23,43*sin(N*360/365,25)

Là toute la formule est en degrés. A la fois ce qu'il y'a dans les parenthèses...mais aussi le résultat d!

Si on est en radians il faudra prendre:

d=(23,43*pi/180)*sin((N*360/365.25)*pi/180)

Le résulta sera lui aussi en radians.

 

EDIT: ah bin pareil, pas vu le "edit"

Erwan