Résolution de l'équation de Kepler
779971
28137466
2008-04-02T15:22:20Z
Badmood
4507
nouveau style de bandeau de portail, voir [[Wikipédia:Prise de décision/Bandeaux de portail]]
{{ébauche|mathématiques}}
Dans le [[mouvement keplerien]], l''''[[équation de Kepler]]''' relie l'[[anomalie moyenne]] M = nt à l'[[anomalie excentrique]] E par l'équation
<center>
:<math> M = E - e \cdot \sin E</math>
</center>
où ''e'' est l'[[excentricité orbitale|excentricité]] de la planète.
'''Résoudre cette équation, c'est trouver E(e,M) :'''
* comme série de Fourier puisque c'est une fonction périodique impaire de M
* comme série de puissance de e, si e < eo := 0.6627..., rayon de convergence de la série.
* comme une valeur numérique avec un nombre de chiffres (d), pour un temps de calcul tc(d) optimisé.
== Série de Fourier ==
C'est [[Joseph-Louis Lagrange|Lagrange]] qui trouve l'expression, bien que le nom J<sub>n</sub>(x) soit associé au nom de [[Friedrich Wilhelm Bessel|Bessel]].
* E-M = fonction impaire périodique de M :
<center>
<math> E-M = e\cdot sinE = 2 \cdot \Sigma_{n=1} \frac{J_n(ne)}{n} \sin(nM)</math>
</center>
*il est facile d'en déduire (a/r)-1 = 2<math> \Sigma_{n=1}J_n(ne)\cdot \cos (nM)</math>
== Série entière de l'excentricité ==
C'est encore Lagrange qui trouve la solution ; et [[Pierre-Simon Laplace|Laplace]] donnera le rayon de convergence : mais [[Augustin Louis Cauchy|Cauchy]], pas content du tout, fonde la théorie des séries analytiques pour résoudre ce problème épineux, qui verra son aboutissement avec les travaux de [[Victor Puiseux|Puiseux]].
<center>
<math> E-M = \Sigma_{n=1} {e^n \over n!}\cdot a_n(M)</math>
</center>
avec <math>\ a_n(M) = D^{n-1} (\sin^n M)</math> et D := opérateur dérivée.
C'est l'application du théorème d'inversion de série de Lagrange.
*Le rayon de convergence de la série est : eo = 0.6627434193
tel qu'indiqué par Laplace ([[1823]]) et démontré par Cauchy et Puiseux :
*<math>e_o = 2 {\sqrt{x(1-x)}\over {1-2x}} </math> et x tel que <math>{1-x \over x} = exp ({2 \over 1-2x})</math>.
Note: cela peut s'écrire différemment : cf discussion.
=== Cas des comètes : <math>e > \ e_o</math> ===
Le premier à se confronter au problème est [[Jeremiah Horrocks|Horrocks]], puis surtout [[Edmond Halley|Halley]] ([[1705]]), pour les calculs sur sa comète d'excentricité e = 0,9673.
Il faut modifier légèrement la solution de Barker (e = 1). Et Bessel([[1805]]) résout ce cas, mais pour e > 0.997
[[Carl Friedrich Gauss|Gauss]] ([[1809]]) s'illustra en donnant une belle solution pour 0,2 < e < 0,95
Autant dire que le voisinage de (0,95 ; 0,98) est fertile à problèmes, en cas d'itération !
== Calcul numérique ==
Les calculs via les [[intégrateur symplectique]]s exigent de rester toujours en butée du nombre de digits, dans le moindre coût de calcul.<br/>
Depuis 300 ans, on cherche la « meilleure » méthode. Elle reste à trouver !
Bien sûr, cela dépend beaucoup du doublet (M,e), M compris entre 0 et Pi et de e, surtout quand e est voisin de 1.
Nijenhuis (1991) adopte la methode de Mikkola (1987) qui est la méthode de Newton d'ordre 4, en choisissant « adéquatement » le germe Eo en fonction du doublet (M,e).
Il est clair que dans les calculs numériques, le volume de calculs est essentiel, autant que le nombre de décimales, vu l'instabilité du système solaire évaluée à un [[exposant de Lyapunov|coefficient de Liapunov]] de 10^(t/5Myr). On se heurte à une muraille exponentielle : difficile d'aller plus loin que 25 Myr, même avec un traitement 128 bits.
Ce sont ces calculs (astronomiques... mais informatisés) qui tournent sur les machines de l'IMCCE-Paris. Le calcul de l'ensoleillement terrestre à la latitude 65°Nord, I(65,t) est calculé et on essaie d'en déduire la corrélation avec le climat passé : l'échelle géologique jusqu'au Néogène (25M ans) en est déduite(échelle géologique Gradstein 2004). Prochaine étape prévue : les 65 M ans.
== Histoire des sciences ==
Avant Kepler, l'équation est déjà étudiée ! bien sûr, pas pour le même problème, mais pour la même équation :
c'est le problème de la réduction des coordonnées locales aux cordonnées géocentriques : il faut réduire la correction de parallaxe. Habash al Hasib s'y est déjà attaqué.
Avant 1700, il y a déjà beaucoup de tentatives : Kepler naturellement, Curtz (1626), Niele, [[Ismaël Bouillau|Bouillau]] (1645, 1657), [[Seth Ward]] (1653), Paganus (1657), Horrebow (1717), [[Jean-Dominique Cassini (Cassini I)|Cassini]] (1669), Newton (1665?), [[Christopher Wren|Wren]] (1658), [[John Wallis|Wallis]] (1659),... De toutes, celle de [[Jeremiah Horrocks]] (1638) est de plus grande beauté. Cf le Colwell, déjà cité.
== Voir aussi ==
*[[mouvement keplerien]]
*[[intégrateur symplectique]]
*[[Jeremiah Horrocks]], cf discussion.
== Bibliographie ==
*Colwell (1993) : ''Solving Kepler's equation over three centuries'', ed Willmann-Bell, ISBN 0-943396-40-9
*Brinkley (1803) : ''Trans roy irish ac'', 7,321-356.
{{Portail|astronomie|mathématiques}}
[[Catégorie:Mécanique céleste]]
[[Catégorie:Loi en physique]]
[[Catégorie:Équation|Kepler2]]