bonjour merci pour tout
j'attends les formules pour le calcul de la position du soleil
a bientot et merci pour les explications
Bonjour,
Dans mon premier post sur ce sujet tu trouveras comment calculer les coordonnées du soleil. Il y a toutes les étapes a suivre.
Ciao.
Bonjour VEGA302
Voilà donc la méthode que j'utilise, toutes les étapes sont nécessaires et serviront aussi pour tous les autres calculs.
Inutile de penser faire les opérations à la calculette...
Pour la position du soleil selon le lieu on a besoin de ces infos.Pi = 3.14159265358979
' D'abord calculer JD en référence à J1900
T1900 = (JD - 2415020.0) / 36525.0
LongSol = 279.696 677 7 + (36 000.768 92 * T1900) + (0.000 302 5 * T1900^2
' ne pas oublier :
LongSol = LongSol - Int(LongSol / 360) * 360 ' débarrassée des révolutions exédentaires
' corrections séculaire
Aa = 153.23 + 22 518.754 1 * T1900
Ba = 216.57 + 45 037.508 1 * T1900
Ca = 312.69 + 32 964.357 7 * T1900
Da = 350.74 + 445 267.114 2 * T1900 - 0.000 144 * T1900^2
Ea = 231.19 + 20.2 * T1900
Ha = 353.4 + 65 928.715 5 * T1900
LongSol = LongSol + 0.001 34 + Cos(Aa)
+ 0.00154 * Cos(Ba)
+ 0.002 * Cos(Ca)
+ 0.001 79 * Sin(Da)
+ 0.001 78 * Sin(Ea)
' à débarrasser des révolutions exédentaires : LongSol = LongSol - Int(LongSol / 360) * 360
AnomSol = 358.478 664 659 167
+ 35 999.049 8 * T1900
- 0.000 15 * T1900^2
+ 0.000 003 3 * T1900^3
' à débarrasser des révolutions exédentaires : AnomSol = AnomSol - Int(AnomSol / 360) * 360
C1 = SIn(AnomSol) * [1.919 460 2 - 0.004 817 * T1900 - 0.000 014 * T1900^2]
+ Sin(2 * AnomSol ) * [0.019 993 - 0.000 101 * T1900]
+ SIn(3 * AnomSol) * 0.000 289
' à débarrasser des révolutions exédentaires
Lg = LongSol + C1
' à débarrasser des révolutions exédentaires
Omega = 259.18 - 1934.142 * T1900
' à débarrasser des révolutions exédentaires
L = Lg - 0.005 69 - 0.004 79 * Sin(Omega)
' à débarrasser des révolutions exédentaires
Rz = -2.465 * Sin(2 * L) + 0.052 * SIn(4 * L)
EquaTemps = 4 * (Rz + C1) * 60 ' donné en secondes
' ****************************** *********
' Application :
' Heure du passage au méridien:
H = 12 + (Longitude * 12 / Pi) + EquaTemps / 3600 + Décalage
' 1) recalculer EquaTemps au moins une fois puisque JD devient plus précis
' ex : Jour Julien du 11 juin 2012 : 2 456 090.000 000 000 000 (à 12H00)
' JJ à H 13.74 h (13h44mn24s) : 2 456 090.072 500 000 000
' 2) recalculer H avec la nouvelle EquaTemps
' on gagne en précision, au 1er passage on trouve Et = -18.0432 s
' au deuxieme, Et = -17.1496 s
' avec Longitude en Radians = */- Longitude en ° / 57.295 779 513 082 320 876 798 154 814 105
' Attention à la saisie des nombres : 445 267.114 2 = 445267.1142
Evidemment il y a toujours moyen de simplifier, mais ces formules n'ont pas cours en vraie astronomie
Voilà une de ces horreurs par exemple, tout juste bonne à établir une courbe à regarder de (très) loin :
J = NumJour ' 1 à 366 le cas échéant
For J = 1 To NumJour ' là on affiche tout
M = 357 + (0.9856 * J) ' Anomalie moyenne en °
' : 1° * Pi / 180 = 1 Radian pour VB6, PHP, etc..
C = (1.914 * Sin(M * PI / 180) + 0.02 * Sin(2 * M * PI / 180))
L = 280 + C + (0.9856 * J) ' en °
R = (-2.465 * Sin(2 * L * PI / 180) + 0.053 * Sin(4 * L * PI / 180))
ET = 4 * (R + C) ' x4 = en minutes, x240 = en sec.;
SinusDeclinaison = Sin(23.440723055 * PI / 180) * Sin(L * PI / 180)
Declinaison = Atn(SinusDeclinaison / Sqr(-SinusDeclinaison * SinusDeclinaison + 1)) * 180 / PI
Print Round(ET, 2); " mn; "; Round(Declinaison, 3); " °"
Next
Expliquer ça pas à pas revient à faire un cours complet, ici ce n'est pas possible.
L'electronique, c'est fantastique.
bonjour j'attends toujours les calcul de la position du soleil dans le ciel de toi Mr. Raikko 21 merci
rebonjour Mr curieuxdenature merci pour tout mais il seulement les résultat mais pas les formules et le chemin à suivre pour s'y retrouver merci pour tout
Bonjour
exprime toi plus clairement parce que là, ou bien le "il" est de trop ou bien il manque quelque chose après ...
Tu as les formules, yapukafoc...
Raikko21 t'as donné un lien pour la marche à suivre, je t'ai donné les équations nécessaires et suffisantes pour mener à bien les calculs de l'équation du temps, le reste coule de source à condition de faire un petit effort.
Donne un exemple étape par étape de ce que tu as trouvé, je te dirais où ça merdouille.
Mais même ça tu es capable de le faire, tu compares mes formules avec les tiennes et tu trouveras tout de suite.
Tu as bien remarqué que mes équations comportent une dizaine de correctifs pourLongSol, ce que tu n'as surement pas...
L'electronique, c'est fantastique.
bonjour
j'ai calculé pour le 13/10/1992 à 0h j'ai trouvé comme Mr raikko21 -13mn42.65 sec mais par logiciel c'est -13mn51sec j'ai aussi cnangé de date pour le 11/06/2012 à 0h cela a donnée -0.366925566 mn à0h et non comme sur le logiciel -0.254933
se que je cherche c'est avoir des formules conforme à la sec pour calculer la même chose que le logiciel merci
bonjour Mr raikko21 merci pour tout mais sachez que je ne saurais jamais si mes calculs sont juste du faite que je n'ai pas quelqu'un à me corriger
si vous voulez on fera un exemple pour la date du 12/06/2012 à 12h TU avec les coordonnées de bourges en cherchant l'équation du temps et la position du soleil pour ce même lieu et même date merci
Bonjour VEGAS302bonjour
j'ai calculé pour le 13/10/1992 à 0h j'ai trouvé comme Mr raikko21 -13mn42.65 sec mais par logiciel c'est -13mn51sec j'ai aussi cnangé de date pour le 11/06/2012 à 0h cela a donnée -0.366925566 mn à0h et non comme sur le logiciel -0.254933
se que je cherche c'est avoir des formules conforme à la sec pour calculer la même chose que le logiciel merci
la précision dépend essentiellement de celle qu'on a sur la longitude apparente du soleil
- raikko21 donne Et pour les équations fournies
- moi je trouve Et = 13 mn 49.67 s avec les algorithmes de Jean Meeus récentes (cette année)
- et 13 mn 41.726 avec celles que je t'ai données au post #95
Pour la longitude apparente du soleil je trouve 199.897 086 196 304 (199° 53' 49.51")
alors que la théorie VSOP42 donne 199.901 716 350 988 (199° 54' 6.179 ") (elle fournie la date des éclipses entre -3000 et +8000...)
Pour ce que tu cherches cela joue sur des fractions négligeables, même avec cet écart tu prévoies une éclipse à 1 mn près dans un siècle.
Dans toutes ces théories, ce qui est important n'est pas la précision absolue mais d'avoir un écart négligeable avec les positions relatives de tous les autres corps célestes.
Autrement dit, qu'importe l'imprécision constatée aujourd'hui entre la longitude apparente du soleil et le calcul, si les autres corps comportent la même erreur alors la théorie fera des prévisions exactes pendant encore de nombreuses années.
Et en 2050 tu verras sortir une autre série de références astronomiques centrée sur J2050.
L'electronique, c'est fantastique.
A propos, tu parles de quel logiciel ?
L'electronique, c'est fantastique.
bonjour le logiciel c'est celui de l'équation du temps
http://www.toocharger.com/fiches/win...temps/3399.htm
merci
bonjour Mr raikko21
je vous ai posé une question sur ce site http://forums.futura-sciences.com/as...-soleil-2.html et j'attends toujours ta réponse merci.
bonjour Mr curieudenature c'est quelle date "T1900 = (JD - 2415020.0) / 36525.0" pour JD moi j'ai pris le 13/10/1992 à 0h et cela m'a donné ET = -813.775532sec pas compris
oui quelle est la date de l'instant du calcul
et aussi voici les formules pour le calcul de la position du soleil d'un instant donné
vérifie si elles sont juste merci
- Calcul de Longitude moyenne du soleil;
L = 279.69668 + 36000.76892*T + 0.0003025 * T^2 en deg
L = 280.46645 36000.76983*T + 0.00030327*T^2 + 0.000000022* T^3
- Calcul de l'anomalie moyenne du soleil;
M = 358.47583 35999.04975 *T – 0.000150* T^2 – 0.0000033* T^3 en deg
- Calcul de l’excentricité de l'orbite terrestre, et l’obliquité de l’écliptique;
e = 0.01675104 – 0.0000418*T – 0.000000126*T^2
- Calcul de l’équation du centre;
C = + (1.919460-0.004789*T -0.000014*T^2 * SIN(M) + (0.020094 -0.000100*T) * SIN(2*M) +0.000293*sin(3*M)
- Calcul de la longitude vraie du soleil;
thétav =L + C
- Calcul de la longitude apparente du soleil (c'est la longitude vraie corrigée de la nutation et de l'aberration);
Omega = 259.18 – 1934.142 * T
thetapp = thétav -0.00569 – 0.00479 * sin(Omega)
- Calcul de l’obliquité de l’écliptique
epsilon = 23.452294- 0.0130125 *T – 0.00000164*T^2 + 0.000000503*T^3 + 0.00256 * cos(Omega)
- Calcul de l’obliquité de l’écliptique
= 23.452294- 0.0130125 *T – 0.00000164*T^2 + 0.000000503*T^3 + 0.00256 * cos()
Bonsoir
à cette date JD par rapport à T1900 est 2 448 908.5
et les routines écrites en C de Jean Meeus donne Et = -829.771 s
par contre je n'ai pas le détail des datas qu'il utilise (je peux décortiquer ses sources mais je ne me souviens pas s'il les a mise en ligne, je dois vérifier...)
J'ai finalement préférer utiliser ses données parce qu'elles sont plus précises que tous les résultats que j'ai eu avec les autres formules, même de celles trouvées dans ses livres.
Cela vient probablement de la précision supérieure du C par rapport aux autres langages.
Voilà un exemple qu'il fournit pour tester:
tu peux constater que ça simplifie sérieusement le problème.' What is the Solar Equation of Time at 6am, 21st January, 2001 Amsterdam time?:
' Amsterdam 's timezone is -1
EOT = -sdxEOT(sdMJDx(Annee, mois, jov, pmmDegre, pmmMinutes, pmmSecondes, -fh))
'equals -2.8175 degrees.
' And converted to minutes:
MsSS = -sdD2Unit(sdxEOT(sdMJDx(Annee, mois, jov, pmmDegre, pmmMinutes, pmmSecondes, -fh)), sdUnitTMin)
'equals -11.2700 minutes.
L'electronique, c'est fantastique.
Si je compare ce que fournit J.Meeus et le programme, je trouve
11 mn 20.89 s (J.M)
11mn 18 s (prog.)
(21/1/2001 à 6h00 Amsterdam)
il n'y a pas de quoi se frapper on ne risque pas de louper le rendez-vous.
Il ne faut pas oublier que Et varie au cours de la journée, celle que fournit mon programme est à la date du 'midi'.
L'electronique, c'est fantastique.
Bonjour,
D’où ce que tu as tiré ces formules? Tu travailles avec le J1900? ou J2000?
Maintenant, tu peux calculer la déclinaison du soleil ainsi que son ascension droite comme suit:
: déclinaison du soleil;
: longitude apparente du soleil;
: obliquité de l’écliptique.
Tu peux aisément passer du système de coordonnées équatorial (déclinaison, ascension droite) vers le système de coordonnées horaire (déclinaison, angle horaire) ou le système de coordonnées horizontal (azimut, hauteur).
Ciao.
bonjour
j'ai travaillé avec J1900 et les formules se trouvent dans la page 55 chapitre 15 du livre de Jean Meeus (calculs astronomiques à l'usage des amateurs) merci
re
donne moi la formule de la longitude apparente du soleil merci
rebonjour merci je l'avais devant les yeux merci encore
et les formules de nutation de la longitude et de l'obliquité merci donne moi un site ou bien autre chose
Je t'ai donné les formules ici.
Ben je travaille avec "Astronomical Algorithms" de Jean Meeus.
Bonjour VEGAS302
les sources y sont bien, il y a environ une centaine de termes pour corriger la longitude solaire dans ses calculs...
Dans les calculs de l'astronome Mr Bretagnon ce n'est pas plus facile, une cinquantaine.
Le problème est qu'ils utilisent une référence différente, J2000 pour J.Meeus et J1900 pour Bretagnon.
Cela donne des formules longues comme un jour sans pain.
L'electronique, c'est fantastique.
bonjour
c'est pour MR. Rakkoi21 moi aussi je travaille avec "Astronomical Algorithms" de Jean Meeus. maios en anglais mais je ne sais pas l'anglais alors je traduit sur google traduction et c'est pas tres forts .
je vous dema,de s'ilm y a un sitye ou je peut trouver le livre de "Astronomical Algorithms" de Jean Meeus. en PDF en français merci.
Bonjour,
Aucune idée, je ne sais pas si ça existe.bonjour
c'est pour MR. Rakkoi21 moi aussi je travaille avec "Astronomical Algorithms" de Jean Meeus. maios en anglais mais je ne sais pas l'anglais alors je traduit sur google traduction et c'est pas tres forts .
je vous dema,de s'ilm y a un sitye ou je peut trouver le livre de "Astronomical Algorithms" de Jean Meeus. en PDF en français merci.
Ciao.
merci et bon apetit mais traduit tout un livre sur google traduction c'est pas donnée