La déclinaison du soleil et la longueur de la nuit

[curieuxdenature]

Bonjour VEGA302

je te le disais dans mon post précédent, avec excel il faut inverser les termes.
Chez microsoft, le gars qui a programmé la fonction ATAN2 devais être accro au beaujolais nouveau parce que pour tout le monde c'est ArcTan2(y,x) alors que chez eux c'est ATAN2(x,y)

Prend un exemple dans le bouquin de Jean Meeus où L = atn(y/x) ce qui revient à poser L = ArcTan2(y,x)

Avec ce classique je retrouve bien le bon résultat

alpha0 = ArcTan2(Cos(0.40907) * Sin(5.478403), Cos(5.478403)) ' = ArcTan2(y,x)

List1.AddItem "RA = " & alpha0 & " rd" ' RA = 5.52144549136948 rd

List1.AddItem "RA = " & alpha0 * 180# / Pi & " °" ' RA = 316.355523467008 °

et 'RA = +5.52144549136948 rd' se conjuge de la même façon que 'RA = -0.761739815810101 rd' si tu lui enlèves 2 Pi
-----

[VEGAS302]

bonjour c'est beau mais je n ai rien su du faite que la virgule qui se trouve après Sin(5.478403) ,

alpha0 = ArcTan2(Cos(0.40907) * Sin(5.478403), Cos(5.478403)) '

que voilà

=Atan^2(Cos(0.40907) * Sin(5.478403) (ici + OU - ou bien / ) , Cos(5.478403)) pour trouver -0.761739 RAD
MERCI

[ansset]

grrrr,
j'aimerai comprendre
atan2 ou atan^2 ???

[VEGAS302]

BONJOUR

atan-2 Merci

[curieuxdenature]

grrrr,
j'aimerai comprendre
atan2 ou atan^2 ???

Bonjour ansset

ArcTan2 est une fonction informatisée de l'Atn des calculettes.
Elle permet de ne pas avoir à manipuler d'autres calculs pour mettre les réponses dans les bons quadrants du cercle.

[curieuxdenature]

bonjour c'est beau mais je n ai rien su du faite que la virgule qui se trouve après Sin(5.478403) ,

alpha0 = ArcTan2(Cos(0.40907) * Sin(5.478403), Cos(5.478403)) '

que voilà

=Atan^2(Cos(0.40907) * Sin(5.478403) (ici + OU - ou bien / ) , Cos(5.478403)) pour trouver -0.761739 RAD
MERCI

Bonjour VEAG302

Si tu n'es pas un peu attentif on va pas en sortir là.
Utilise donc :

Atn(Cos(0.40907) * Sin(5.478403) / Cos(5.478403))

dans ce cas précis ArcTan2 ne sert à rien (mais si tu veux faire d'autres calculs ce ne sera pas toujours le cas).

[ansset]

Bonjour ansset

ArcTan2 est une fonction informatisée de l'Atn des calculettes.
Elle permet de ne pas avoir à manipuler d'autres calculs pour mettre les réponses dans les bons quadrants du cercle.

oui, j'avais compris,
je réagissais à "l'écriture" précédente de VEGA qui mélangeait les deux sans ses équations.
( l'écriture , pas lui forcement )

[curieuxdenature]

OK il me semblais bien que tu avais déjà posé la question.

[VEGAS302]

bonjour

dans ce site http://forums.futura-sciences.com/at...soleil-001.jpg
en ligne 21 22 23 Le resultat que j'ai trouve n est pas le meme comment on t il fait merci

[curieuxdenature]

Bonjour VEAGS302

laisse tomber cette référence, elle est truffée d'erreurs, déjà quand on voit que le gars n'a pas été fichu de placer la longitude et la latitude aux bons endroits, ça ne laisse rien présager de bon pour la suite.
Même la correction de la réfraction de l'air est 60 fois trop grande...

[VEGAS302]

BONJOUR

donne moi un exemple avec cette formule merci
GMSThrs = (6.697374558hrs + 0,06570982441908 × (JJ - J2000) + 1,00273790935 × heures UTC + 0,000026 × T^2)
par exemple pour le 02/02/2013 à 11h 30 TU
MERCI

[curieuxdenature]

Bonsoir

J = 2
M = 2
A = 2013
he = 11
mn = 30

aaa = -1
bbb = 8239.5
ccc = -3522
ddd = 63
JJ-J2000 = 4779.5
T = 0.130868697239334

GMST = 20.354 676 596 238

[VEGAS302]

BONSOIR
Tous mes remerciements et a la prochaine question bonne nuit

[VEGAS302]

bonjour
pourquoi -1 d'ou vient il , j'ai trouvé 9.237959...aprs simplication avec mod 24
Si cela ne vous gene pas donne moi la formule exact avec ses mod24 c'est une formule en heure merci

[curieuxdenature]

Bonjour

si tu nous donnais les tiens, vu que c'est toi qui fais l'erreur, on y verrait plus clair.

[VEGAS302]

Bonjour
(((6,697374558+0,0657098244190 8*(4779,5)+1,00273790935*11.50 +0,000026*0,130855578370979^2)/24)-13)*24 = 20,28896677

je n'ai pas trouvé comme vous T MERCI

[curieuxdenature]

Salut VEGAS302

je ne sais pas où tu as trouvé cette formule.

Essaye plutôt ça :
GMST = 6.697374558 + 0.06570982441908 * D0 + Tu * 1.0027379093508 + 0.000026 * T ^ 2
GMST = GMST - Int(GMST / 24) * 24 (pour ramener l'heure entre 0 et 24)

tu dois retrouver
GMST = 20.354 676 596 425

(un petit conseil, soit tu utilises le point soit la virgule comme séparateur décimal, mais pas les deux dans la même équation.
Et il n'y a que les français qui se servent de la virgule, ce n'est pas une riche idée...)

Sinon tu avais raison, le -1 est une erreur, suivant le cas on doit trouver +1, 0, -1 mais là c'était 0
ce qui donne
D0 = 4780.5 (nombre de jours julien écoulés depuis le 1.5 janvier 2000)
T = 0.130896075747205 (jour julien séculaire)
ce qui ne change que la 9eme décimale de l'heure sidérale de Greenwich.

[VEGAS302]

bonjour

pour le 2/2/2013 à 11.50 h ils ont trouvé pour le GSMT = 18hr 19min 02sec en page 7 dans ce site => http://www.precisedirections.co.uk/S...0Astronomy.pdf
meme là je n'ai pas trouve comme vous pour le T qui est 4780.50/36525 = 0,13088295687885 et pour GMST => 20,3546765963260 qui donne 20H21MN16.84SEC merci

[curieuxdenature]

Re

je ne me suis pas amusé à refaire tous ses calculs, voilà les miens

 Cliquez pour afficher

avec sa méthode :
J = 2
M = 2
A = 2015 ----> attention dans le document il est noté 2013...
he = 11
mn = 30
Zone = 2
été/hiver = 0
longitude = -1.58111133333333 h; = -23.71667° ---> attention au signe
latitude = 37.96667°

Tu = 9.5h du lieu (TU Greenwich - Zone - heure d'été)
aaa = 0
bbb = 8973.5
ccc = 3526
ddd = 63
JJ-J2000 = 5510.5
T = 0.150 882 386 493 269

GMST = 18.317 372 750 076 h heure sidérale de Greenwich
GMST = 18h 19 mn 02.542

LMST = 19.898 484 083 409 h heure sidérale locale
----------------------------------------------

avec la méthode VSOP87 pour le reste :

RA = 21.044 629 3082 h
RA = 21h 02m 40.666s
RA = 315.669 439 6231 °
RA = 315° 40' 09.983"

Dec = -16.851 796 3432 °
Dec = -16° 51' 06.467"

Sidereal time at Greenwich :
Mean: 18.317 372 7469h heure sidérale moyenne
True: 18.317 468 5130h heure sidérale vraie
Diff: +0.0000957661h

Soleil à Athènes :

Hauteur : 32.851 384 772 1064 °
: 32°51'04.985"

Azimut : 160.323 554 205 258 °
: 160°19'24.795"

Longitude Soleil
313.204 370 1127°

A priori, il n'y a pas d'erreur grossière à relever, c'est correct.

[VEGAS302]

bonjour merci de ton aide

voilà mon nouveau probleme j'ai fait tous les calculs mais je n'est pas trouve le même resultat je demande ton aide une nouvelle fois

T = 0.08537 est juste

Coordonnées équatoriales (J2000) J2000 :
RA : 19 :25 :27.913
DEC : 42 :47 : 3.700

Coordonnées équatoriales vraies (J2008.53) :
Précession variation des paramètres de précession (” par siècle) :
m = 4612.4362 + 2.79312* T − 0.000278*T^2
n = 2004.3109 − 0.85330 *T − 0.000217 T^2
∆α = m/100 + n/100 sin α tan δ
∆δ = n/100 cos α coordonnées a une époque t seront :
α(t) = α(J2000.0) + ∆α (t − 2000.0)
δ(t) = δ(J2000.0) + ∆δ (t − 2000.0)
Précession annuelle,
m : 1.92349 secondes = 0.00053 heure
n : 7.30202” = 0.00203 degré

résultat de la précession :
RA : 19 :25 :44.330
DEC : 42 :48 : 6.025

le site c'est celui ci

http://userpages.irap.omp.eu/~jlebor...s/corriges.pdf

merci pour tout pas seulement le resultat mais la façon de calculer

[curieuxdenature]

Bonjour VEGAS302

je retrouve bien les mêmes résultats.

pour le calcul des precessions, pas de problème particulier, yaka suivre les formules.
La feinte à éviter, c'est de ne pas oublier de convertir m (qui est donné en secondes d'arc) en secondes de temps (à diviser par 15)
je retrouve m=28.85234 " qui devient bien 1.9234958 secondes pour Delta AD PAR AN

enfin, en multipliant m et n par le nombre d'années (8.53) tu as bien ce qui es donné en solution.

soit m=1.92349 / 3600 * 8.53 = 0.0045576 heure; ajouté à 19.42442028 H = 19.42897788 -> 19h25mn44.32 s

idem pour n = 7.302024" / 3600 * 8.53 = 0.01730173 °; +42.7843611 ° = 42°48'06""

[VEGAS302]

BONJOUR

mais vous n avez pas pose l operations comment vous aviez trouvé pour
m = 4612.4362 + 2.79312 T − 0.000278 T 2 :
et pour
n = 2004.3109 − 0.85330 T − 0.000217 T 2 pour une operation integrale merci

[curieuxdenature]

T est donnée dans l'énoncé : 8.53 ans après 2000, puisque les équations sont en " par siècle, tu poses T=0.0853
et attention T 2 c'est T au carré = T²

[VEGAS302]

BONJOUR

en toute franchise moi je cale ici (4612,4362+2,79312*0,085366188-0,000278*0,085366188^2) =4612,674636
je ne sais si je dois /3600 où quoi faire merci
bonne journée

[curieuxdenature]

Bonjour VEGAS302

Primo, c'est T=0.85300000 dans l'énoncé.
Toujours dans l'énoncé, il est question d'une dérive donnée en secondes d'angle par siècle.
Automatiquement le résultat devra être converti en degrés à la fin des calculs (et en heures pour AD).

 Cliquez pour afficher

date : 13/07/2008 01:58:48.000
Jour Julien saisi : 2454660.5825
T2000 saisi : 8.52999999999878E-02
An : 2008
N° du jour 195

Données du problème :
291.366304167; 291°21'58.6950"; 19:25:27.9130
42.78436111; 42°47'03.7000"; 02:51:08.2467

DeltaT : 8.53000 ans

m : 4612.67455111325
n : 2004.23811193109

DeltaAD : 1.923 4896 s = 0.00053 H par an ( à * par 8.53 )
DeltaDec : 7.302 0240 " = 0.00203 ° par an ( à * par 8.53 )

291.434668193905; 291°26'04.8055"; 19:25:44.3204
42.801662850116; 42°48'05.9863"

on a en ligne de visée l'étoile V* RR Lyrae invisible à l'oeil nu.
encore appelée HIP 95497 située entre Vega de la Lyre et Deneb

http://en.wikipedia.org/wiki/RR_Lyrae

[VEGAS302]

REBONJOUR

m : 4612.67455111325
n : 2004.23811193109
J'ai trouve ça mais comment convertir pour trouver pour m= 1.923 4896 s et pour n= 7.302 0240 " merci

[curieuxdenature]

pour n tu le trouves directement en appliquant l'équation donné dans ton lien :
DeltaDEC = n/100 * cos(AD) = 7.302024 "

Pour m tu divises le résultat de l'équation par 15 :
DeltaAD = m/100 + n/100 * Sin(AD) * Tan(DEC) = 28.85234 " d'angle
28.85234 " / 15 = 1.92349 seconde de temps

si tu avais essayé de diviser ton résultat par 1.92349 donné en solution, tu aurais vu que le coeff était 15, faut tout tester quand c'est ainsi, c'est comme ça qu'on avance...

[curieuxdenature]

je me suis un peu attardé sur ton lien, et en effet la question fait partie d'un ensemble, la donnée T=0.0853 est un résumé, je pensais que c'était une question à part alors que référence est bien le 15/07/2008 à 22h TU (ce qui donne T= 0.0853775952543877)


Mais bon, ça ne change pas grand chose aux résultats, des fractions de secondes qui ne risquent pas de faire rater le pointage au télescope...
comme l'indique la correction suivante :

 Cliquez pour afficher

date : 15/07/2008 22:00:00.000
Jour Julien saisi : 2454663.41666667
T2000 saisi : 8.53775952543877E-02
An : 2008
N° du jour 198

Données du problème :
291.366304167; 291°21'58.6950"; 19:25:27.9130
42.78436111; 42°47'03.7000"

DeltaT : 8.53776 ans
m : 4612.67476784242
n : 2004.23804571618
DeltaAD : 1.923 4898 s
DeltaDec : 7.302 0237 "

291.434730389426; 291°26'05.0294"; 19:25:44.3353
42.8016785884996; 42°48'06.0429"

[VEGAS302]

bonsoir

moi j'ai trouve pour m = 4612,674668 et pour n = 2004,238046
comment aviez vous fait les conversion m et n merci

[curieuxdenature]

Bonjour VEGAS302

je ne fais pas les calculs à la main, c'est trop périlleux, je laisse ça à l'ordi qui ne se trompe pas avant la 15ième décimale.
Pour la conversion, je ne fais que suivre les formules données dans ton lien, je t'ai donné la marche au message #57.
Tu remplaces AD par 291.366304167°
et DEC par 42.78436111°

Aucune astuce à part la conversion d'angles vers les heures. (360° = 24H donc 15° = 1H), c'est basique.