Voir la source de Mouvement képlérien

En [[astronomie]], plus précisément en [[mécanique céleste]], le '''mouvement képlérien''' correspond à une description du mouvement d'un astre par rapport à un autre respectant les trois [[lois de Kepler]]. Pour cela il faut que l'interaction entre les deux astres puisse être considérée comme purement newtonienne, c'est-à-dire qu'elle varie en raison inverse du carré de leur distance, et que l'influence de tous les autres astres soit négligée. Il est fréquent que l'un des deux corps (qualifié alors de « central ») soit de masse (beaucoup) plus importante que l'autre, et soit donc dominant gravitationnellement, tel que celui d'une [[étoile]] par rapport à l'une de ses [[planète|planètes]], d'une planète par rapport à l'un de ses satellites [[satellite naturel|naturels]] ou [[satellite artificiel|artificiels]], etc. Le mouvement képlérien décrit alors le mouvement du corps relativement moins massif.

Le mouvement képlérien est en fait un cas particulier du [[problème à deux corps]], dans le cas d'une force variant comme l'inverse du carré de la distance, ce qui correspond à la force d'[[Loi universelle de la gravitation|interaction gravitationnelle]] entre les deux corps, s'ils sont supposés [[Point matériel|ponctuels]] ou encore à symétrie sphérique <ref group="N" name="APM">L'approximation qui consiste à considérer les corps comme ponctuels est valable si leurs dimensions respectives peuvent être considérées comme petites devant la distance qui les sépare. Si l'on considère qu'ils possèdent la symétrie sphérique, ce qui est vérifié en première approximation pour de nombreux corps célestes, notamment les planètes, le champ de gravitation créé à l'extérieur de chacun des astres est le même que celui créé par un point matériel placé au centre de corps, et affecté de sa masse totale (cf. par exemple Perez, ''op. cit.'', chapitre 6, pages 77-78).</ref>, et constituant un système isolé<ref name="goldstein">Voir par exemple: {{Goldstein}}, {{2e|édition}}, chapitre 3.</ref>. Or le problème à deux corps peut toujours être ramené à celui du mouvement d'une particule (dite ''fictive'') dans le référentiel du centre de masse des deux corps, les trajectoires des corps « réels » étant homothétiques de celles de cette particule fictive.

Dans le cas du [[Champ gravitationnel|champ de gravitation]] « newtonien », l'étude du mouvement de ce corps fictif dans le champ de gravitation « créé » par un centre fixe à l'origine est appelé '''''problème de Kepler'''''<ref name="duriez">Luc Duriez, ''Cours de mécanique céleste classique'', université Lille-I / IMCCE, 2002, ([http://www.imcce.fr/fr/formations/cours/CoursMC_Duriez/mc/CoursMCecr_Duriez.pdf lire en ligne]), chapitre 3.</ref>{{,}}<ref name="perez">Cf. Perez, ''Cours de physique : mécanique'', {{4e|édition}}, Masson, Paris, 2001, chapitre 12.</ref> et le mouvement résultant est précisément le mouvement képlérien<ref name="duriez" />. Il est intéressant de noter qu'en dépit de ces appellations, [[Johannes Kepler|Kepler]] ne formula pas ainsi la situation, mais déduisit ses [[Lois de Kepler|lois du mouvement des planètes]] en utilisant les observations réalisées par [[Tycho Brahe]], notamment sur le mouvement de la planète [[Mars (planète)|Mars]], alors que ni la loi de la gravitation universelle ni le principe fondamental de la dynamique n'avaient encore été découverts. De fait, ce problème dit de Kepler n'a été formulé correctement qu'avec [[Isaac Newton|Newton]], sous une forme géométrique<ref>Cf. Isaac Newton, [http://www.e-rara.ch/doi/10.3931/e-rara-1743 |''Principes mathématiques de la Philosophie Naturelle''] , livre 1, édition de 1756, trad. [[Émilie du Châtelet]], livre 1, notamment les sections II à III (numérisation [[e-rara]]).</ref>.

Il est alors possible de montrer que les trajectoires des corps célestes sont des [[conique|coniques]], en particulier des ''[[ellipse (mathématiques)|ellipses]]''<ref group="N">Du moins bien sûr pour la situation où les deux corps demeurent liés entre eux, qui correspond sur le plan énergétique à une [[énergie mécanique]] (notée ''H'') négative du système, autrement dit que l'[[énergie potentielle]] d'interaction gravitationnelle entre les deux corps, notée ''V'' (négative par convention) « l'emporte » sur l'énergie cinétique totale du système, notée ''T'', toujours positive, évaluée dans le référentiel barycentrique. Dans ce cas il est clair qu'aucun des deux corps ne pourra aller à l'infini l'un de l'autre, puisque la condition <math>T=H-V\ge 0 </math> ne sera plus vérifiée dans le cas où ''H<0'' au-delà d'une certaine distance ''r'' telle que <math>V(r)=H</math>, ''V'' tendant vers ''0'' quand ''r'' croît. En revanche si ''H>0'' ou si ''H=0'', cette condition restera vérifiée même à l'infini.</ref>, ce dernier cas très important étant souvent désigné sous le nom d{{'}}''orbite de Kepler''. Dans le cas fréquent où l'un des deux corps est de masse beaucoup plus importante que l'autre, il peut être considéré comme immobile à l'un des foyers de la conique, qui représente alors la trajectoire du corps « périphérique ». 

Comme tout modèle, le problème de Kepler correspond bien sûr à une situation idéalisée, puisque l'influence de tous les autres corps célestes est négligée, tout comme les déviations au caractère non-newtonien de la force de gravitation<ref group="N">Ces déviations sont à la fois celles dues au caractères non-ponctuel et non-sphérique (aplatissement lié à la rotation, irrégularité, etc.) des corps célestes en jeu, mais aussi à la présence de forces non-gravitationnelles, telles que la [[Traînée|traînée atmosphérique]], ou la [[pression de radiation solaire]] pour un [[satellite artificiel]], où les deux forces peuvent avoir une grande importance, ou encore pour une [[comète]], en ce qui concerne la pression de radiation.</ref>. Toutefois cette approximation est souvent vérifiée avec une très bonne approximation dans de nombreux cas. Ainsi dans le [[Système solaire]] le [[Soleil]], qui concentre 99,86 % de la masse totale du système, domine gravitationnellement tous les autres corps. Il est donc possible de considérer avec une bonne approximation le mouvement de chaque planète ou de tout autre corps du système comme képlérien. L'influence des autres corps peut ensuite être prise en compte comme des ''perturbations'' de l'orbite képlérienne<ref group="N">Au besoin en tenant compte des corrections liées à la [[relativité générale]], comme pour [[Mercure (planète)|Mercure]], cf. [[Tests expérimentaux de la relativité générale#Avance du périhélie de Mercure]]</ref>, celle-ci constituant ainsi « l'orbite de base » à partir de laquelle il est possible de construire par approximations successives l'orbite « réelle » du corps céleste considéré. Par suite, l'étude du mouvement képlérien est d'une importance capitale en mécanique céleste.
<gallery widths="220" heights="200&quot;">
Image:Terrestrial planet size comparisons.jpg|<center>Les planètes telluriques, de gauche à droite :<br/> Mercure, Vénus, Terre, Mars</center>
</gallery>

== Aspects historiques ==

{{Article détaillé|Histoire du mouvement képlérien}}

De nombreux modèles furent proposés durant l'[[Antiquité]] pour représenter les mouvements des [[planète]]s<ref group="N">À cette époque, cette notion inclut également le [[Soleil]] et la [[Lune]]. Au sens moderne, cinq planètes étaient alors connues : [[Mercure (planète)|Mercure]], [[Vénus (planète)|Vénus]], [[Mars (planète)|Mars]], [[Jupiter (planète)|Jupiter]] et [[Saturne (planète)|Saturne]]. Le mot ''planète'' dérive de l'expression πλανήτης αστήρης ''planêtês astêrês'' « astre (ou étoile) errant », en effet celles-ci se distinguent des étoiles « fixes » par leur mouvement apparent sur la sphère céleste au cours du temps.</ref>. Ces modèles sont à la fois des « ''[[Histoire des cosmologies scientifiques pré-relativistes|cosmologies]]'' », c'est-à-dire des « [[représentations du monde]] » (ou de l'Univers, ou du κόσμος, ''kósmos'', « monde ordonné » des Grecs<ref>Ainsi selon [[Socrate]] : {{Citation bloc|À ce qu’assurent les doctes [[École pythagoricienne|pythagoricien]]s, [[Calliclès]], le ciel et la terre, les Dieux et les hommes sont liés entre eux par une communauté, faite d’amitié et de bon arrangement, de sagesse et d’esprit de justice, et c’est la raison pour laquelle, à cet univers, ils donnent, mon camarade, le nom de cosmos, d’arrangement, et non celui de dérangement non plus que de dérèglement.|[[Platon]], ''[[Gorgias (Platon)|Gorgias]]'', 507e - 508a}}.</ref>) et des tentatives de comprendre et de prévoir les phénomènes astronomiques, notamment pour déterminer les dates des [[éclipse|éclipses]] ou la construction des [[calendrier|calendriers]]. Tous ces systèmes sont '''''[[géocentrisme|géocentriques]]''''', c'est-à-dire qu'ils placent la Terre au centre de l'Univers. La seule exception, qui n'aura aucune postérité immédiate, est le modèle [[héliocentrisme|héliocentrique]] d'[[Aristarque de Samos]], proposé vers -280. 
[[Image:Cellarius ptolemaic system.jpg|thumb|left|Représentation du système de Ptolémée, selon l'ouvrage ''Andreas Cellarius Harmonia Macrocosmica'', 1660/61.]]
Selon [[Simplicius (philosophe)|Simplicius]] (fin {{s-|V}} - début {{sap-|VI}})<ref>Ainsi dans son ''Commentaire à la Physique d'Aristote'', celui-ci écrit : {{citation bloc|Platon […] pose alors ce problème aux mathématiciens : quels sont les mouvements circulaires et uniformes et parfaitement réguliers qu'il convient de prendre pour hypothèses, afin que l'on puisse sauver les apparences que les astres errants présentent ?|[[Simplicius (philosophe)|Simplicius]]|''Commentaires sur la physique d'Aristote''}}.</ref> c'est [[Platon]] (427-327 {{av JC}}) qui aurait proposé à son élève [[Eudoxe de Cnide]] (408-355 {{av JC}}) d'étudier le mouvement des planètes en n'utilisant que des ''mouvements circulaires et uniformes'', considérés comme parfaits <ref name="EU">Voir par exemple ''Encyclopédia Universalis'', édition 2002, volume 3, article « Astronomie et astrophysique » {{ISBN|2-85229-550-4}}, ou encore ''Notionnaires Universalis - Idées'', {{ISBN|2-85229-562-8}}, Encyclopedia Universalis France SA, Paris, 2005.</ref>. 
[[Image:Claudius Ptolemaeus.jpg|thumb|Portrait allégorique de Ptolémée]]
[[Image:Epicycle et deferent.png|left|thumb|Mouvement d'une planète dans le système géocentrique de Ptolémée, décrit à l'aide d'un épicycle et d'un déférent.]]

La difficulté de décrire précisément les mouvements des planètes, notamment les phénomènes de ''[[rétrogradation (astronomie)|rétrogradation]]'', conduit à des représentations complexes. Les connaissances astronomiques du monde gréco-romain sont résumées au {{s-|II}} de notre ère par [[Ptolémée]] (vers 90 - 168 ap. JC), dans un ouvrage en Grec transmis par  les Arabes sous le nom de l'[[Almageste]]<ref group="N">Le titre original est Μαθηματική σύνταξις (Mathématikế sýntaxis), ''Composition mathématique''. Le nom courant de cet ouvrage via d'une latinisation de l'arabe al-majisṭī (المجسطي). Cet ouvrage est disponible en trois volumes sur Internet :[https://gallica.bnf.fr/ark:/12148/bpt6k3974x première partie], [https://gallica.bnf.fr/ark:/12148/bpt6k29016x seconde partie] et [https://gallica.bnf.fr/ark:/12148/bpt6k290178 troisième partie].</ref> Connu sous le nom de '''modèle de Ptolémée''', la représentation du système solaire et du mouvement des [[planète|planètes]] (ainsi que de la [[Lune]] et du [[Soleil]]) utilise comme ses prédécesseurs un modèle [[géocentrisme|géocentrique]], et un système élaboré de sphères en rotation circulaire et uniforme, les [[épicycle|épicycles et déférents]], introduit par [[Hipparque (astronome)|Hipparque]] ({{-s-|II}}), qu'il perfectionne en introduisant la notion d'[[équant]], qui est un point distinct du centre du cercle déférent par rapport auquel une planète, ou le centre d'un épicycle, se déplace à vitesse uniforme. 

Le système de Ptolémée va dominer l'astronomie pendant quatorze siècles. Il donne des résultats satisfaisants malgré sa complexité, au besoin en modifiant et raffinant le modèle des épicycles, déférents, et points équants. Considéré comme compatible avec la philosophie d'Aristote, le géocentrisme devient doctrine officielle de l’Église en Europe au cours du [[Moyen Âge]]. 

[[Image:Nikolaus Kopernikus.jpg|thumb|Nicolas Copernic]][[Image:CopernicSystem.png|thumb| Version simplifiée du système héliocentrique de Copernic, image extraite de son ouvrage ''De revolutionibus Orbium Coelestium''.]]
On doit à [[Nicolas Copernic]] (1473-1543) la remise en cause du dogme géocentrique. Il publie en 1543, l'année de sa mort,  son ouvrage majeur ''[[De revolutionibus|De revolutionibus Orbium Coelestium]]'' où il propose un système héliocentrique, dans lequel les planètes et la Terre se déplacent selon des orbites circulaires, parcourues à vitesses constantes, la Lune étant le seul astre tournant autour de la Terre. Bien qu'imparfaite, cette vision s'avère très féconde : les mouvements des planètes sont plus simples à décrire dans un référentiel héliocentrique<ref name="EU" />. L'ensemble des irrégularités de mouvements telles que les rétrogradations s'explique uniquement par le mouvement de la Terre sur son orbite, plus précisément en termes modernes par l'effet du passage du référentiel héliocentrique au [[référentiel géocentrique]]. Le système de Copernic permet également de supposer que les étoiles de la « sphère des fixes » sont à une distance bien plus grande de la Terre (et du Soleil) que l'on le supposait jusqu'alors, pour expliquer l'absence d'effet observé alors (''[[parallaxe]]'') du mouvement de la Terre sur la position des étoiles. Il est à noter qu'initialement le système de Copernic, qui sur la pratique astronomique consistait à échanger les positions de la Terre et du Soleil, ne suscita pas une opposition de principe de l’Église, jusqu'à ce que celle-ci s'aperçut que ce modèle remettait en cause la philosophie d'Aristote<ref name="EU" />.

[[Image:Portrait Confused With Johannes Kepler 1610.jpg|thumb|left|Johannes Kepler vers 1610.]]
[[Johannes Kepler]] (1571-1630) perfectionnera ce modèle, grâce à l'analyse soigneuse des observations précises de son maître [[Tycho Brahe]] (1541-1601), notamment concernant le mouvement de la planète [[Mars (planète)|Mars]] publie ses trois célèbres lois (Cf. [[Lois de Kepler]]) en 1609, 1611, 1618<ref name="EU" /> :
* '''première loi''' : {{citation|Les planètes décrivent des ellipses dont le Soleil occupe l'un des foyers.}}
* '''deuxième loi''' : {{citation|Le rayon vecteur reliant le centre de la planète au foyer décrit des aires égales en des temps égaux.}}
* '''troisième loi''' : {{citation|Les cubes des demi-grands axes des orbites sont proportionnels au carré des périodes de révolution.}}

Le modèle de Kepler est d'une importance capitale en astronomie : il rompt avec le dogme géocentrique, mais également avec celui de l'utilisation du mouvement circulaire et uniforme pour décrire le mouvement des planètes, car non seulement les orbites ne sont plus des cercles, mais la deuxième loi implique que les vitesses de révolution sur celles-ci ne sont plus constantes, les planètes devant se mouvoir plus rapidement au voisinage de leur passage au plus près du foyer (le ''périastre'') qu'à celui du point le plus éloigné (l{{'}}''apoastre'').

[[Isaac Newton]] (1642-1727), démontre ces lois empiriques en 1687 (voir ''{{lang|la|[[Philosophiae Naturalis Principia Mathematica]]}}'' ; [[Démonstration des lois de Kepler]]), en proposant une loi de la gravitation dite ''universelle'', car s'appliquant aussi bien pour les mouvements des planètes que pour l'explication de la chute d'un corps sur la Terre. Ceci marque le début d'une nouvelle ère : la mécanique céleste et la mécanique classique fondées sur le [[calcul infinitésimal]].

Dans sa théorie, Newton utilise la notion controversée d'action ''instantanée à distance'' : la gravitation agit instantanément en 1/''r''². C'est le célèbre {{Citation|[[hypotheses non fingo]]}}. La « résistance cartésienne », plus la très grande difficulté mathématique des ''Principia'' provoquera un temps de réception assez long de son œuvre, et il faudra attendre les ouvrages d'[[Leonhard Euler|Euler]], de [[Colin Maclaurin|McLaurin]] et de [[Alexis Claude Clairaut|Clairaut]], pour éclaircir la situation. Néanmoins les plus grands esprits de l'époque ([[Christian Huygens|Huygens]], [[Gottfried Wilhelm Leibniz|Leibniz]], etc.) reconnaissent immédiatement la portée des ''Principia''.

La [[Constante gravitationnelle|constante de gravitation]] sera évaluée par [[Henry Cavendish|Cavendish]] en 1785 : G = {{unité|6.674 e-11 N m2 kg-2}}. Il convient de noter que la valeur de cette constante tout comme celle de la masse M de nombreux objet du système solaire, n'est connue qu'avec une assez faible précision, alors que les valeurs de μ = GM, appelé [[paramètre gravitationnel standard]] sont souvent déterminées avec une grande précision. Ainsi pour la Terre GM = {{unité|398600.441,8 ± 0.000,8|km3 s−2}}, soit une précision de 1 pour {{formatnum:500000000}} alors que M = {{unité|5.973,6 e24 kg}}, soit une précision de l'ordre de 1 pour {{formatnum:7000}}.

== Orbites et lois de Kepler ==

=== Résultats généraux du problème à deux corps ===
{{Article détaillé|Problème à deux corps}}
==== Notion de particule fictive ====
Le mouvement képlérien est un cas particulier important d'une situation générale dans lequel deux corps, considérés comme ponctuels<ref group="N" name="APM" /> et formant un système isolé<ref group="N">Typiquement, cette approximation consiste à négliger l'influence des autres corps pouvant être au voisinage du système des deux corps considérés. Ceci est valable si un des deux corps est de masse prépondérante par rapport à tous les autres, comme le Soleil dans le Système solaire, ou encore si les autres corps sont à des distances telles que leur influence peut être négligée : par exemple, la Terre n'exerce qu'une faible influence sur le système constitué par Neptune et ses satellites. Ceci étant, même en présence d'un corps extérieur à l'influence gravitationnelle importante il est possible de se placer dans la situation d'un problème à deux corps. Ainsi pour le système Terre-Lune, il est facile de vérifier que le Soleil exerce sur cette dernière une influence gravitationnelle près de 3,7 fois supérieure à celle de la Terre. Il pourrait alors sembler illusoire de traiter le système Terre-Lune comme un problème à deux corps. Toutefois, ceci est possible en se plaçant dans le référentiel du barycentre du système (voire de celui lié au centre de la Terre), en translation par rapport au référentiel de Copernic (héliocentrique). Si ce dernier référentiel est galiléen, ce n'est pas le cas du premier, qui est accéléré par l'attraction du Soleil, qui ne peut être négligée. Toutefois dans ce référentiel il est possible de montrer (cf. l'article [[force de marée]]) que l'influence du Soleil sur la Lune se réduit à un terme différentiel, dit de marée, bien plus faible que celui lié à l'attraction gravitationnelle Terre-Lune, laquelle varie en première approximation comme l'inverse du carré de la distance. Il sera ainsi possible de traiter le mouvement de la Lune autour de la Terre comme étant képlérien, puis de tenir compte des corrections liées aux différents termes (de marée, de non-sphéricité du potentiel…) dans un deuxième temps.</ref>, sont en interaction gravitationnelle mutuelle<ref group="N">Dans le cas le plus général, en interaction mutuelle via une force dérivant d'un potentiel ''V(r)''.</ref>. Il est alors possible (cf. [[Problème à deux corps|article détaillé]]) de montrer que ce problème général peut toujours être ramené à celui du mouvement d'un seul corps, dénommé ''particule fictive'', affecté de la masse <math>m=\frac{m_1m_2}{m_1+m_2}</math>, appelée [[masse réduite]]<ref group="N">La notation ''m'' est utilisée ici pour la masse réduite de préférence à la notation habituelle ''μ'' généralement utilisée, qui sera réservée à celle du [[paramètre gravitationnel standard]] du système, défini par <math>\mu=GM</math>, où <math>M=m_1+m_2</math> est la masse totale du système, et ''G'' la [[constante gravitationnelle]].</ref>, se déplaçant dans le champ de gravitation créé à l'origine par un corps ponctuel immobile affecté de la masse totale du système <ref name="duriez" />, à condition de se placer dans le [[référentiel (physique)|référentiel]] du centre de masse des deux corps (référentiel barycentrique). Les [[orbite]]s réelles de chacun des corps célestes dans ce référentiel sont alors [[homothétie|homothétiques]] de celle de cette particule fictive.  

Le problème de Kepler correspond simplement au cas où ce champ de gravitation est de type « newtonien », c'est-à-dire variant comme l'inverse du carré de la distance<ref name="duriez" />. Les résultats généraux de l'étude du problème à deux corps permettent de préciser les caractéristiques des trajectoires possibles, le champ newtonien possédant une constante du mouvement spécifique qui permet d'obtenir aisément l'équation polaire des orbites correspondantes.

==== Propriétés générales du mouvement ====
La résolution du problème de Kepler se ramène donc à résoudre l'équation du mouvement de la particule fictive ''P'' dans le champ de gravitation « causé » par le centre de force, (affecté de la masse totale ''M'' du système), donnée par l'application du [[principe fondamental de la dynamique]], soit après réarrangement des termes :
:<math>\frac{d^2\vec{r}}{dt^2}=-\mu\frac{\vec{r}}{r^3}</math>
avec <math>\vec{r}=\overrightarrow{OP}</math> et <math>\mu=GM</math>, appelé ''[[paramètre gravitationnel standard]]'' (''M'' étant la masse totale du système).

Le caractère central de la force implique la conservation du [[moment cinétique]] de la particule fictive (voir aussi [[Moment cinétique spécifique]]), ce qui a deux conséquences :
* '''planéité des orbites''', du fait que le rayon-vecteur <math>\vec{r}</math> et le vecteur vitesse <math>\vec{v}</math> de la particule sont à tout instant orthogonaux à un vecteur de direction constante ;
* que le mouvement vérifie la '''loi des aires''' <ref name="goldstein" />{{,}}<ref name="duriez"/>{{,}}<ref name="perez" />, aussi connue sous le nom de '''''deuxième loi de Kepler''''' (1611): {{Citation|Le rayon vecteur reliant le centre de la planète au foyer décrit des aires égales en des temps égaux.}}.

De façon plus précise, si <math>d\mathcal{A}</math> désigne l'aire élémentaire balayée pendant ''dt'' par le rayon vecteur, cette loi implique que la vitesse aréolaire du corps <math>\dot{\mathcal{A}}=\frac{d\mathcal{A}}{dt}</math> est constante, avec <math>\dot{\mathcal{A}}=C/2</math>, où <math>C=\|\vec{r}\wedge\vec{v}\|</math>.

Ces propriétés sont générales pour le problème à deux corps, quelle que soit la forme du potentiel d'interaction entre les deux corps. Elles sont valables que le corps possède ou non une trajectoire fermée, ou qu'il puisse ou non aller à l'infini.

=== Orbites de Kepler ===
[[Image:OrbitalEccentricityDemo.svg|thumb|Représentation des allures d'une orbite de Kepler pour diverses valeurs de l'excentricité ''e''. En bleu, trajectoire hyperbolique (''e > 1''). En vert, ''e = 1'', trajectoire parabolique. Le cas de l'orbite elliptique (''0< e <1'') est en rouge, l'orbite circulaire d'excentricité nulle étant en gris.]]
Dans le cas du problème de Kepler, il existe une constante du mouvement spécifique, appelée [[vecteur de Runge-Lenz]], ou encore (sous forme adimensionnelle) '''''[[vecteur excentricité]]'''''<ref name="duriez" /> :
:<math>\vec{e}=\frac{\dot{\vec{r}}\wedge\vec{L}}{GM\mu}-\vec{e}_r</math> où <math>\vec{e}_r=\vec{r}/r</math>.

La projection de cette grandeur vectorielle, contenue dans le plan de l'orbite, sur le rayon vecteur <math>\vec{r}</math> permet d'obtenir aussitôt l'équation de la trajectoire en [[coordonnées polaires]] :
:<math>r(\theta)=\frac{p}{1+e\cos{\theta}}</math>,
qui est celle d'une ''[[conique]]'', avec :
* <math>\theta=\theta(t)</math> angle polaire entre <math>\vec{r}</math> et <math>\vec{e}</math>, appelé '''''anomalie vraie''''' ;
* <math>e=\|\vec{e}\|</math>, ''[[excentricité orbitale|excentricité]]'' de l'orbite ;
* <math>p=\frac{C^2}{\mu}</math>, ''paramètre'' de l'orbite.

La forme de l'orbite correspondante dépend de la valeur de l'excentricité ''e'', qui elle-même dépend de l'[[énergie mécanique]] ''H'' du système (par nature conservatif), puisqu'il est possible d'établir la relation<ref name="duriez" /> :
:<math>e=\sqrt{1+\frac{2HC^2}{\mu^2 m}}</math> (''C'' étant la [[constante des aires]], ''m'' la masse réduite du système, toutes deux définies plus haut).

Il est possible de distinguer les cas suivants (cf. figure ci-contre) :
* ''e > 1'' : orbite (ou plutôt trajectoire) '''''hyperbolique'''''. Physiquement ce cas correspond à ''H > 0'', et implique que le corps peut aller en l'infini avec une vitesse (donc une [[énergie cinétique]]) non nulle, l'énergie potentielle de gravitation tendant vers zéro dans cette limite ;

* ''e = 1'' : trajectoire '''''parabolique'''''. Il s'agit du cas limite ou ''H = 0'', le corps peut donc aller en l'infini, mais avec une vitesse nulle ;

* ''0 < e < 1'' : orbite '''''elliptique''''' (dite de Kepler). Dans ce cas ''H <0'', telle que <math>H > - \frac{\mu^2 m}{2C^2}</math> pour que la relation précédente entre ''e'' et ''H'' ait un sens physique ;

* ''e = 0'' : orbite '''''circulaire''''', cas limite du précédent où <math>H=- \frac{\mu^2 m}{2C^2}</math>.

Le cas de l'orbite ''elliptique'' (et à la limite, circulaire) correspond à la '''''première loi de Kepler''''' (1609), qui s'énonce originellement : {{Citation|Les planètes décrivent des ellipses dont le Soleil occupe l'un des foyers.}}. De façon plus rigoureuse, l'un des foyers de l'orbite est occupé par le centre de masse du système, qui pour les planètes du Système solaire est pratiquement confondu avec le centre du Soleil<ref group="N">En revanche, la différence peut être importante si la masse d'un deux corps n'est négligeable par rapport à celle de l'autre, par exemple pour une [[étoile]] binaire.</ref>.

Dans tous les cas, le vecteur excentricité est dirigé vers le point d'approche minimale entre l'un des foyers (en pratique, l'astre « central », de masse dominante) et la particule fictive (en pratique, le corps « en orbite » autour du corps « central »), appelé ''périapse'' (ou périastre)<ref group="N">Il s'agit là du terme général, qui peut prendre des noms particuliers suivant le corps céleste considéré : ainsi pour le mouvement de révolution de la Terre autour du Soleil, le périastre s'appellera ''périhélie'', pour une satellite de la Terre, le ''périgée''.</ref>. La distance correspondante est <math>r_{min}=\frac{p}{1+e}</math>.

Du fait de la loi des aires, la vitesse du corps céleste sur son orbite n'est ''pas constante'' (sauf si le mouvement est circulaire). Il est facile de voir par exemple dans le cas d'une orbite elliptique que la vitesse du corps devra être plus rapide au périastre qu'à l'apoastre, pour que le rayon vecteur puisse balayer des aires égales en des temps égaux.

{{souligner|'''Remarques :'''}}
* L'existence du vecteur excentricité comme [[intégrale première]], en plus du moment cinétique ou de l'énergie totale, est ''spécifique'' au problème de Kepler, plus précisément à la dépendance en <math>1/r</math> du potentiel d'interaction gravitationnelle. Comme il est souligné dans l'article relatif au [[problème à deux corps]] ou au [[vecteur excentricité]], ceci implique de par le [[Théorème de Noether (physique)|théorème de Noether]] que le problème de Kepler possède une symétrie additionnelle. En effet, la [[conservation du moment cinétique]] est liée à l'invariance par rotation du problème à deux corps<ref name="landau">Cf. notamment à ce sujet {{Landau|tome 1}}, chapitre II.</ref>, tandis que celle de l'énergie totale est liée à l'invariance par translation dans le temps<ref name = "landau" />. La symétrie associée au vecteur excentricité est plus complexe, car elle ne s'interprète correctement qu'en dimension 4<ref>Cf. à ce sujet {{Goldstein}}, {{2e|édition}}, chapitre 9, paragraphe 8, pour une explication détaillée montrant l'identité entre les commutateurs des matrices de rotation générant le groupe ''SO(4)'' des rotations propres dans un espace quadridimensionnel et des [[crochet de Poisson|crochets de Poisson]] entre les composantes du vecteur de Runge-Lenz et du moment cinétique.</ref>.
* Le fait que les orbites de Kepler dans le cas lié soient fermées est également spécifique du mouvement képlérien, n'étant absolument pas valable pour un problème à deux corps où le potentiel d'interaction ''U(r)'' est quelconque, mais uniquement pour le potentiel newtonien en ''1/r'' ou le potentiel harmonique isotrope en ''r<sup>2</sup>'' ([[théorème de Bertrand]]).
* L'importance du vecteur excentricité fut soulignée par [[Jakob Hermann|Hermann]] (1710 & 1713), redécouverte par [[Pierre-Simon de Laplace|Laplace]], puis réutilisée par [[Carl Runge|Runge]] et [[Wilhelm Lenz|Lenz]] (Cf. [[vecteur de Runge-Lenz]]).

=== Cas de l'orbite elliptique ===
[[Image:Ellipse orbit parameters.png|thumb|left|Caractéristique principale de l'orbite de Kepler. Le vecteur excentricité est en vert sur la figure, l'anomalie vraie étant notée ν plutôt que θ.]]
Dans le cas de l'orbite de Kepler elliptique, <math>\vec{e}</math> est dirigé selon le grand axe de l'orbite, et il existe un point d'éloignement maximal entre le foyer « centre » et le corps céleste (en toute rigueur la particule fictive), appelé ''apoastre'', qui correspond à ''θ=π'', soit <math>r_{max}=\frac{p}{1-e}</math>. Le paramètre ''p'' de l'ellipse correspond à la distance entre le foyer et le corps en orbite quand celui-ci est en quadrature (''θ=±π/2''). L'ellipse possède deux foyers symétriques, par rapport à son centre ''C'', son [[petit axe]] étant perpendiculaire au grand axe en ce point.

Il est possible de définir les grandeurs caractéristiques suivantes (cf. l'article [[Ellipse (mathématiques)|Ellipse]]) :
* '''demi-grand axe''' ''a'' : il correspond au « rayon » de l'orbite, et est égal à la moitié de la longueur du grand axe. Il est facile de vérifier les relations :
:<math>p=a(1-e^2)</math>, <math>r_{min}=a(1-e)</math>, <math>r_{max}=a(1+e)</math> ;
* '''demi-petit axe''' ''b'' : moins important que le précédent, il lui est lié par la relation <math>b=a\sqrt{1-e^2}</math>, par ailleurs <math>b=p/(\sqrt{1-e^2})</math>, ce qui donne également la relation entre ''a'', ''p'' et ''b'' : <math>p=b^2/a</math> ;
* la distance du foyer au centre, notée ''c'', donnée par la relation <math>c=ea=\sqrt{a^2-b^2}</math>, il est facile de vérifier également que <math>r_{max}-r_{min}=2c</math>.

Dans le cas limite de l'orbite circulaire de rayon ''R'' et de centre ''O'', ''e = 0'', et donc ''a=b=p=R'', ''c = 0'', le point ''C'' et les foyers ''F'' et ''F''' sont confondus avec ''O''.

Pour les huit planètes du Système solaire, il est possible de considérer avec une bonne approximation que le Soleil en raison de sa masse est immobile au foyer de chacune des orbites, lesquelles sont pour la plupart quasi circulaires vu les faibles excentricités, sauf pour Mercure (''e=0,206'').

=== Troisième loi de Kepler ===
L'orbite elliptique est fermée, par suite le mouvement de l'astre sur celle-ci est ''périodique'', de période ''T'' appelée ''période de révolution''. Durant une période, le rayon-vecteur balaie la totalité de l'aire de l'ellipse, soit <math>\pi ab</math>. Or du fait de la constance de la vitesse aréolaire (deuxième loi de Kepler), avec <math>\dot{\mathcal{A}}=C/2</math>, cette aire est également égale à <math>CT/2</math>, ce qui donne en utilisant la relation <math>b^2=ap</math> avec <math>p=C^2/\mu</math> :
:<math>\pi^2a^3p=C^2T^2/4</math>, soit encore <math>\pi^2a^3C^2/\mu=C^2T^2/4</math>.
Finalement il vient une relation entre le carré de la période de révolution et le cube du demi-grand axe de l'orbite de Kepler :
:<math>\frac{T^2}{a^3}=\frac{4\pi^2}{\mu}</math>,
qui correspond à la '''''troisième loi de Kepler''''' (1618) :{{citation|Les cubes des demi-grands axes des orbites sont proportionnels au carré des périodes de révolution.}}.

Cette loi, spécifique au mouvement képlérien, est d'une importance considérable : il est en effet facile de déterminer par des observations la période de révolution d'un astre tel qu'une planète ou un satellite. Il est alors possible d'en déduire la valeur du demi-grand axe de l'orbite, ainsi que celle de <math>\mu=GM</math>.

== Description des orbites de Kepler - Équation de Kepler ==
Si les résultats de l'étude du problème de Kepler permettent de connaître la forme des orbites possibles, ils ne sont pas directement applicables en astronomie. En premier lieu, il faut pouvoir repérer dans l'espace l'orbite du corps céleste considéré, en relation avec un plan de référence, en général l'[[écliptique]], qui est le plan de l'orbite terrestre. Les paramètres qui permettent de le faire sont appelés les ''éléments d'orbite'' ou ''[[éléments orbitaux]]'', dont la détermination résulte des observations. En second lieu, une fois ces éléments d'orbite connus, il faut disposer d'une relation donnant la position de l'astre au cours du temps sur cette orbite, ce qui permettra notamment de déterminer les [[Éphéméride (astronomie)|éphémérides]] du corps céleste considéré. C'est l''''équation de Kepler''' qui permet d'obtenir une telle relation entre position de l'astre sur son orbite et temps.

Dans la suite, seul le cas de l'orbite de Kepler elliptique est considéré.

=== Éléments d'orbite ===
Il a été indiqué que dans le cas du mouvement képlérien le problème à deux corps possède trois intégrales premières, dont deux vectorielles (<math>\vec{C}</math>, lié au moment cinétique, et <math>\vec{e}</math>), et une scalaire (''H'', l'énergie totale), ce qui fait au total 6 constantes scalaires<ref group="N">Et non pas 7, car <math>\vec{e}</math>, contraint à être dans le plan de l'orbite, n'a que deux composantes indépendantes.</ref>.

Toutefois de par la relation <math>e=\sqrt{1+\frac{2HC^2}{\mu m}}</math> entre les valeurs de ces intégrales premières, une orbite de Kepler donnée (quelle que soit sa forme) est caractérisée par seulement 5 grandeurs scalaires indépendantes, qui peuvent être choisies comme étant<ref name="duriez" />{{,}}<ref group="N">Bien entendu, on suppose là que ''C#0'', le cas ''C=0'', dit dégénéré, correspondant à un mouvement rectiligne.</ref>:
* les trois composantes de <math>\vec{C}</math> (qui définissent l'orientation dans l'espace du plan de l'orbite) et les deux de <math>\vec{e}</math> (qui définissent l'orientation de l'orbite dans ce plan);
* les trois composantes de <math>\vec{C}</math>, la valeur de ''H'' (liée à celle de ''e'' et de ''a'', donc définissant la forme de l'orbite), et un angle donnant la direction du périapse dans le plan orbital;
* ou encore deux angles repérant la direction de <math>\vec{C}</math>, la valeur du demi-grand axe ''a'', celle de l'excentricité ''e'', un angle déterminant la direction du périapse dans le plan orbital…

Ces différentes possibilités sont à la base de la détermination des '''[[orbite|éléments d'orbite]]''', ensemble de 5 données géométriques permettant de caractériser une orbite donnée. En effet, la détermination de l'orbite d'un corps céleste donné résulte des observations et mesures réalisées.

En pratique les éléments d'orbite sont donnés par les deux valeurs du demi-grand axe ''a'' et de l'excentricité ''e'', ainsi que de trois angles définis sur la figure ci-après<ref name="dumoulin">Cf. Dumoulin et Parisot, ''Astronomie pratique et informatique'', Masson, Paris, 1987, {{ISBN|2-225-81142-3}} .</ref>:
[[Image:Orbite.png|Définition des différents éléments d'une orbite de Kepler elliptique.]]
* l''''[[Inclinaison orbitale|inclinaison]]''' ''i'' (entre 0 et 180 degrés) est l'angle que fait le plan orbital avec un ''plan de référence'', qui pour les orbites planétaires dans le système solaire est en général le plan de l'[[écliptique]], c'est-à-dire le plan de l'orbite terrestre;
* la '''[[longitude du nœud ascendant]]''' ''Ω'': il s'agit de l'angle entre la direction du point vernal (position du soleil lors de l'[[équinoxe]] de printemps<ref group="N">Du fait de la [[précession des équinoxes]] la position de ce point, qui sert également d'origine pour l'[[ascension droite]], varie légèrement au cours du temps, par suite il faut préciser pour quel date on donne la valeur de ''Ω''.</ref>) et la ligne des nœuds, c'est-à-dire de l'intersection entre le plan de l'orbite considéré et du plan de référence, le nœud ascendant étant le point où le corps céleste passe au nord du plan de référence.
* l''''[[argument du périastre]]''' ''ω'': il s'agit de l'angle formé par la ligne des nœuds et la direction du périastre, dans le plan orbital.

Les deux premiers de ces angles permettent de fixer la position du plan de l'orbite dans l'espace, par rapport au plan de référence, le dernier angle permettant d'en préciser l'orientation, et les données de ''a'' et de ''e'' la forme. Bien entendu ces différents paramètres peuvent être obtenus par les observations.

Un autre angle indiqué sur la figure, l'[[anomalie moyenne]], verra sa définition et son utilité précisés plus bas.

=== Mouvement d'un corps sur une orbite de Kepler - Équation de Kepler ===
{{article détaillé|Équation de Kepler}}
En astronomie, on cherche souvent à obtenir des '''[[éphéméride (astronomie)|éphémérides]]''' du mouvement des différents corps célestes, tels que les planètes du système solaire. L'obtention de celles-ci nécessite non seulement de caractériser l'orbite de ce corps à partir de ses éléments, déterminés notamment par les observations astronomiques, mais également de connaître son mouvement sur cette orbite, ce qui permettra de connaître sa position à un moment donné.

Or l'équation polaire <math>r(\theta)=\frac{a(1-e^2)}{1+e\cos{\theta}}</math>, liant la distance ''r'' du corps par rapport au foyer de l'orbite à l'anomalie vraie ''θ'', ne contient aucune information dynamique. C'est la deuxième loi de Kepler qui contient cette information: la constance de la vitesse aréolaire au cours du déplacement de la planète sur l'orbite implique que si l'instant de passage au périastre est pris pour origine des temps, l'aire balayée par le rayon vecteur au foyer du corps à un instant ''t'' quelconque, soit <math>\mathcal{A}(t)=\dot{\mathcal{A}}t=\frac{Ct}{2}</math> est égale à celle du « secteur d'ellipse » correspondant à l'anomalie vraie au même instant <math>\theta (t)</math>. Ceci permettrait en théorie d'obtenir une relation entre ''θ'' et ''t'' et donc la position de l'astre sur son orbite, supposée bien déterminée géométriquement, à tout instant.

Cependant, il n'est pas aisé en pratique d'exprimer l'aire d'un secteur d'ellipse d'angle ''θ'' donné, contrairement au cas d'un cercle, et le problème doit être abordée en introduisant d'autres paramètres:
[[Image:Diagram Anomalies Kepler orbit.svg|lang=fr|thumb|Représentation des différentes « anomalies » introduites pour l'orbite de Kepler elliptique. <math>\nu</math> est ici l'anomalie vraie, ''E'' et ''M'' étant respectivement l'anomalie excentrique et l'anomalie moyenne.]]
* le premier est l''''[[anomalie excentrique]]''', notée ''E''. Cet angle est défini en introduisant le cercle dit auxiliaire (ou exinscrit) de rayon ''a'', le demi-grand axe, et de centre celui de l'ellipse. Si ''P'' est la position de l'astre d'anomalie vraie ''θ'' à un instant ''t'' donné, le projeté orthogonal de ''P'' sur le cercle auxiliaire dans une direction perpendiculaire au demi-grand axe est localisé par son angle au centre ''E'' (cf. figure ci-contre).
:La raison de l'introduction de cet angle est simple: comme il repère un point sur un cercle, il est nettement plus aisé de calculer l'aire d'une portion de cercle qu'il définit que de la portion d'ellipse associée à la valeur correspondante de l'anomalie vraie, qui seule possède une réalité physique. L'anomalie excentrique est un paramètre géométrique, donnant la position d'un astre sur son orbite de Kepler, plus aisé à utiliser que l'anomalie vraie ''θ'', à laquelle elle est liée par la relation:
:<math>\tan{\frac{E}{2}}=\sqrt{\frac{1-e}{1+e}}\tan{\frac{\theta}{2}}</math>.
:Il est également possible d'obtenir facilement les relations suivantes entre ''r'', ''θ'', et ''E''<ref name="dumoulin" />:
:<math>\begin{cases} r\cos{\theta}=a\left(\cos{E}-e\right) \\ r\sin{\theta}=a\sqrt{1-e^2}\sin{E} \\ r=a\left(1-e\cos{E}\right) \end{cases}</math>,
:ces relations montrent que la connaissance de ''E'' pour une orbite donnée (donc avec ''a'' et ''e'' connus) permet de connaître ''r'' et ''θ'', et donc la position de l'astre sur son orbite.

* le second paramètre est l''''[[anomalie moyenne]]''', notée ''M''. Il s'agit d'un paramètre dynamique, directement lié à l'aire balayée par le rayon vecteur reliant l'astre au foyer de l'orbite. La vitesse angulaire moyenne du corps sur l'orbite est appelée ''moyen mouvement'', et est définie par <math>n=\frac{2\pi}{T}</math>, ''T'' étant la période de révolution, or d'après la deuxième loi de Kepler il vient aussitôt <math>n=\sqrt{\frac{\mu}{a^3}}=\sqrt{\frac{GM}{a^3}}</math>. Si l'origine des temps est prise au dernier passage au périastre un corps hypothétique parcourant l'orbite à vitesse angulaire constante aurait à l'instant ''t'' parcouru l'angle <math>M=nt</math>, qui par définition est l''''anomalie moyenne'''. Or, d'après la loi des aires l'aire balayée par le rayon vecteur reliant le corps « réel » au foyer, entre l'instant de son dernier passage au périastre et l'instant ''t'', est également donnée par <math>\mathcal{A}(t)=\frac{\pi ab}{T}t=\frac{ab}{2}M</math>. Par suite, l'anomalie moyenne est donc directement liée, à une constante multiplicative près, à l'aire balayée par le rayon vecteur reliant le foyer et l'astre, d'où l'intérêt de cette grandeur.

L'équation de Kepler s'obtient alors en exprimant géométriquement l'aire <math>\mathcal{A}</math> du secteur d'ellipse correspondant à ''θ'' donné en fonction de ''E'': il vient (cf. [[Équation de Kepler|article détaillé]]) <math>\mathcal{A}=\frac{ab}{2}\left(E-e\sin{E}\right)</math>.
En éliminant ''ab'' entre cette expression de <math>\mathcal{A}</math> et celle obtenue avec l'anomalie moyenne, il vient aussitôt l'équation reliant le paramètre géométrique ''E'' au paramètre dynamique ''M'', appelée '''[[équation de Kepler]]''':
:<math>M=E-e\sin{E}</math>.
Il convient de préciser qu'en dépit de sa forme, qui ne fait apparaître que des angles, l'équation de Kepler est bien une relation entre ''E'', qui fixe la position de l'astre sur son orbite, et le temps ''t'', l'origine des dates étant le dernier passage au périastre. Il est d'ailleurs possible de substituer dans cette équation l'expression de l'anomalie moyenne ''M'' en fonction de ''t'' et de la période de révolution ''T'', ce qui permet d'obtenir: <math>\frac{2\pi}{T}t=E-e\sin{E}</math>.

Connaissant les différents éléments de l'orbite d'un corps donné ainsi que la période de révolution de celui-ci (et donc son moyen mouvement), la détermination de l'instant de passage au périastre par les observations permet de connaître ensuite par l'équation de Kepler la position de l'astre sur son orbite à tout instant ultérieur, ou antérieur. Bien entendu, cela suppose que l'orbite réelle du corps est bien képlérienne, ce qui ne sera pas rigoureusement exact, du fait des perturbations induites par les autres corps dont l'influence a été négligée: toutefois celles-ci peuvent être prises en compte par des techniques appropriées.

== Résolution de l'équation de Kepler ==
Depuis trois siècles, les astronomes et les mathématiciens se sont préoccupés de trouver ''θ(t)'', puisque l'observable la plus facile à repérer est la position du Soleil dans le ciel. Passer de ''θ'' à ''E'' et réciproquement est aisé. En dépit de sa simplicité apparente, l'équation de Kepler n'a pas de solution analytique connue, et sa résolution requiert donc d'utiliser soit des développements en série, soit des méthodes numériques. De nombreuses méthodes ont été proposées au cours des siècles<ref>Cf. Peter Colwell, ''Solving Kepler's equation over three centuries'', Ed. Willmann-Bell, Richmond, 1993, {{ISBN|978-0943396408}} .</ref>, et les principales sont rappelées ci-après.

=== Développements usuels du mouvement képlérien ===
Les considérations des paragraphes précédents montrent qu'étant donnée une date ''t''<ref group="N" name="origine">Il est rappelé que cette date a pour origine le dernier passage au périastre.</ref>, et donc une anomalie moyenne ''M'', l'inversion de l'équation de Kepler <math>M=E-e\sin{E}</math> permet d'obtenir l'anomalie excentrique ''E'', qui ensuite permet d'obtenir ''r'' et ''θ'' en utilisant les relations fondamentales:
:<math>\begin{cases} r\cos{\theta}=a\left(\cos{E}-e\right) \\ r\sin{\theta}=a\sqrt{1-e^2}\sin{E} \\ r=a\left(1-e\cos{E}\right) \end{cases}</math>,
<math>X=r\cos{\theta}</math> et <math>Y=r\sin{\theta}</math> étant d'ailleurs les coordonnées de l'astre sur son orbite, par rapport au foyer de celle-ci.

Il est possible d'obtenir des développements en séries de trigonométriques en fonction de l'anomalie moyenne ''M'' (donc du temps<ref group="N" name="origine" />) de ces différentes grandeurs, dont les principaux sont donnés ci-après.

==== Développement de l'anomalie excentrique ''E'' ====
L'équation de Kepler permet d'exprimer la différence entre ''E'' et ''M'': <math>E-M=e\sin{E}</math>. Or le second membre est une fonction de ''M'', qui est donc de par les propriétés du sinus périodique (période ''2π'') et ''impaire''. Il est donc possible d'exprimer <math>f(M)=E-M=e\sin{E}</math> sous la forme de [[série de Fourier|séries de Fourier]], soit:
:<math>f(M)=\sum_{k=1}^{+\infty}{b_k \sin{kM}}</math>, en tenant compte du caractère impair de <math>f(M)</math>,
les coefficients <math>b_k</math> étant donnés par:
:<math>b_k=\frac{2}{\pi}\int_{0}^{2\pi} f(M) \sin{kM}\,dM=\frac{2}{\pi}\int_{0}^{2\pi} e\sin{E}\sin{kM}\,dM</math>.

À partir de là, différentes manipulations permettent d'obtenir le développement en séries de ''E'', dans lequel interviennent les [[fonction de Bessel|fonctions de Bessel]] <math>J_k(ke)</math>:
:<math>E=M+2\sum_{k=1}^{+\infty} \frac{J_k(ke)}{k}\sin{kM}</math>
{{démonstration|contenu=
L'intégration par parties de l'expression de <math>b_k</math>, avec <math>\sin{kM}\,dM=-\frac{d\left(\cos{kM}\right)}{k}</math> permet d'obtenir:
:<math>b_k=\frac{e}{k\pi}\int_{0}^{2\pi} \cos{E}\cos{kM}\,dE</math>,
or <math>\cos{E}\cos{kM}=\frac{1}{2}\left(\cos{(E+kM)}+\cos{(E-kM)}\right)=\frac{1}{2}\left(\cos{((k+1)E-(ke)\sin{E})}+\cos{((1-k)E+(ke)\sin{E})}\right)</math>, ce qui donne par substitution dans l'expression précédente, et compte tenu de l'expression intégrale de la fonction de Bessel de première espèce <math>J_k(x)=\tfrac{1}{2\pi}\int_{0}^{2\pi} \cos{\left(kt-x\sin{t}\right)}\,dt</math> ainsi que de la relation de récurrence entre ces fonctions, <math>\frac{2k}{x}J_k(x)=J_{k-1}(x)+J_{k+1}(x)</math>:
:<math>b_k=\frac{e}{k}\left(J_{k-1}(ke)+J_{k+1}(ke)\right)=\frac{2}{k}J_k(ke)</math>, d'où le résultat.}}

Comme <math>r=a\frac{dM}{dE}</math>, on en déduit <math>\frac{dE}{dM}=\frac{a}{r}</math>, il est possible d'obtenir par dérivation du développement précédent celui de <math>\frac{a}{r}</math>:
:<math>\frac{a}{r}=\frac{a}{r}=1+2\sum_{k=1}^{+\infty} J_k(ke)\cos{kM}</math>.

==== Développements de cos(E) et sin(E) ====
Le développement de <math>\sin{E}=\frac{E-M}{e}</math> s'obtient facilement à partir de celui de ''E'' :
:<math>\sin{E}=2\sum_{k=1}^{+\infty} \frac{J_k(ke)}{ke}\sin{kM}</math>.

Le développement en série de <math>\cos{E}</math> s'obtient lui en tenant compte du fait qu'il s'agit d'une fonction paire de ''M'', de période ''2π''. Comme précédemment il est possible de la décomposer en séries de Fourier, faisant intervenir les termes en <math>\cos{kM}</math> du fait de la parité, il vient en procédant de même que pour le développement de ''E'':
:<math>\cos{E}=-\frac{e}{2}+\sum_{k=1}^{+\infty} {\frac{1}{k}\left(J_{k-1}(ke)+J_{k+1}(ke)\right)\cos{kM}}</math>.

De ce développement, il est facile de déduire celui du rayon vecteur ''r'' puisque <math>r=a\left(1-e\cos{E}\right)</math>, il vient :
:<math>r=a\left(1+\frac{e^2}{2}-\sum_{k=1}^{+\infty}{\frac{e}{k}\left(J_{k-1}(ke)+J_{k+1}(ke)\right)\cos{kM}}\right)</math>,
il est possible de voir que la valeur moyenne de ''r'' est <math>a(1+e^2/2)</math>. De la même façon il sera possible d'exprimer les coordonnées au foyer de l'axe sur son orbite, <math>X=a\left(\cos{E}-1\right)</math> et <math>Y=a\sqrt{1-e^2}\sin{E}</math>.

Les développements en séries de <math>\sin{kE}</math> et <math>\cos{kE}</math> peuvent également être obtenus, en procédant de la même façon. Celui de <math>\sin{kE}</math> présente un intérêt particulier pour exprimer l'anomalie vraie en fonction de ''M'':
:<math>\sin{kE}=k\sum_{n=1}^{+\infty}{\frac{\left(J_{k+n}(ne)+J_{n-k}(ne)\right)}{n}\sin{nM}}</math>.
{{démonstration|contenu=
Comme il s'agit d'une fonction impaire de ''E'', donc de ''M'', qui s'annule en même temps et a même période ''2π'' que cette dernière, <math>\sin{nE}</math> avec ''n'' entier peut être décomposé en séries de Fourier, faisant intervenir uniquement les coefficients de sinus des multiples de ''M'':
:<math>\sin{kE}=\sum_{n=1}^{+\infty}{b_n \sin{nM}}</math>, avec <math>b_n=\frac{1}{\pi}\int_{0}^{2\pi} {\sin{kE}\sin{nM}\,dM}</math>.
Comme <math>\sin{nM}=\frac{-d(\cos{nM}}{n}</math>, il est possible d'intégrer par partie, le terme tout intégré étant nul, il vient :
:<math>b_n=\frac{k}{n\pi}\int_{0}^{2\pi} {\cos{kE}\cos{nM}\,dE}</math>.
Or comme <math>\cos{kE}\cos{nM}=\left(\cos{(kE+nM)}+\cos{(kE-nM)}\right)/2</math>, soit compte tenu de l'équation de Kepler <math>\cos{kE}\cos{nM}=\left(\cos{\left((k+n)E-ne\sin{E}\right)}+\cos{\left((n-k)E-ne\sin{E}\right)}\right)/2</math>,
par substitution dans l'intégrale et compte tenu de l'expression <math>J_n(x)=\frac{1}{2\pi}\int_{0}^{2\pi}{\cos{\left(nt-x\sin{t}\right)}\,dt}</math> il vient finalement le résultat cherché.}}

==== Développements de l'anomalie vraie - équation du centre ====
La relation entre anomalie vraie ''θ'' et l'anomalie excentrique ''E'' est une relation entre les tangentes des angles moitié :
:<math>\tan{\frac{\theta}{2}}=\sqrt{\frac{1+e}{1-e}}\tan{\frac{E}{2}}</math>.

Il est possible de montrer que cette relation peut s'écrire sous la forme d'un développement en série de la forme :
:<math>\theta=E+\sum_{k=1}^{+\infty}{\frac{q^k}{k}\sin{(2kE)}}</math>
Il est possible à partir de cette relation d'exprimer ''θ'' en séries entières, il vient:
:<math>\tan{\theta}=E+2\sum_{k=1}^{+\infty} {\frac{q^k}{k}\sin{kE}}</math>, avec <math>q=\frac{1-\sqrt{1-e^2}}{e}</math>.
{{démonstration|contenu=
La relation <math>\tan{\frac{\theta}{2}}=\sqrt{\frac{1+e}{1-e}}\tan{\frac{E}{2}}</math> est de la forme <math>\tan{x/2}=p\tan{y/2}</math> avec <math>p=\sqrt{\frac{1+e}{1-e}}</math>. Comme <math>\tan{x/2}=\frac{e^{ix/2}-e^{-ix/2}}{i(e^{ix/2}+e^{-ix/2})}</math>, il vient après réarrangement:
:<math>\frac{e^{ix}-1}{e^{ix}+1}=p\frac{e^{iy}-1}{e^{iy}+1}</math>, soit encore <math>e^{ix}=\frac{(p+1)e^{iy}-(p-1)}{(p+1)-(p-1)e^{iy}}=e^{iy}\frac{(p+1)-(p-1)e^{-iy}}{(p+1)-(p-1)e^{iy}} =e^{iy}\frac{1-qe^{-iy}}{1-qe^{iy}}</math>.
En prenant le logarithme de cette expression, il vient <math>ix=iy+\ln{(1-qe^{-iy})}-\ln{(1-qe^{iy})}</math>, or <math>\ln{(1-X)}=\sum_{k=1}^{+\infty}{\frac{X^n}{n}}</math>, par suite et compte tenu de <math>\sin{ny}=\frac{e^{iny}-e^{-iny}}{2i}</math>, il vient finalement le résultat cherché.
}}
La substitution dans cette expression des développements de ''E'' et ''sin kE'' permet d'obtenir l'expression de <math>\theta-M</math>, appelée '''équation du centre'''<ref name="dumoulin" />{{,}}<ref name="duriez" />:
:<math>\theta-M=2\sum_{n=1}^{+\infty}{\frac{1}{n}\left(J_n(ne)+\sum_{k=1}^{+\infty}{\left[q^k\left(J_{n-k}(ne)+J_{n+k}(ne)\right)\right]}\right)\sin{nM}}</math>


==== Développement en série de l'excentricité e ====
Les développements en séries précédents, obtenus à partir des séries de Fourier, sont uniformément convergents. Les coefficients de ces séries, qui font intervenir des fonctions de Bessel de type <math>J_n(ne)</math> peuvent également être exprimés sous la forme de séries entières de l'excentricité ''e'', et les termes réarrangés de façon à obtenir des développements en séries de chacune des diverses grandeurs en fonction de ''e''. Toutefois le réarrangement des termes peut conduire à des séries qui ne seront plus convergentes pour toutes les valeurs de l'excentricité.

Par exemple, pour l'anomalie excentrique ''E'' le développement correspondant se met sous la forme:
<math>E = M + \sum_{k=1}^{+\infty} A_k(e) \sin{kM}</math>,
avec:
:<math>A_1(e)=e-\frac{e^3}{8}+\frac{e^5}{192}+...</math>, <math>A_2(e)=\frac{e^2}{2}-\frac{e^4}{6}+\frac{e^6}{48}+...</math>, etc.
Ces séries convergent si ''e < 0.662''<ref group="N">
''Remarque d'histoire des sciences'' : ce résultat fut trouvé par Cauchy qui inventa exprès la méthode des séries entières de la variable complexe z pour trouver ce résultat, obtenant E, tel que E = <math>2 \sqrt{x(1-x) \frac{1}{1-2x}}</math> et x tel que (1-x)/x = exp (2/(1-2x)), soit E = 0.662 743 419 349 182 </ref>{{,}}<ref group="N"> Remarque numérique : à l'heure actuelle, nombreux sont les astronomes qui voudraient avoir des [[Éphéméride (astronomie)|éphémérides]] très précises, ne serait-ce que pour la Terre. Un exemple parmi d'autres : en [[climatologie]], il est intéressant de regarder le cycle des glaciations, ce qui implique une bonne connaissance des [[paramètres de Milanković]]. Ceci a été calculé par l'IMCCE de Paris pour le Néogène (les 25 derniers millions d'années), ce qui a servi pour calibrer l’échelle géologique mondiale 2004. Toutefois remonter plus loin dans le passé (par exemple jusqu'à 100 Myr) requiert des méthodes rapides de résolution de l'équation de Kepler.</ref>. L'intérêt de telles séries est que dans le cas où ''e'' est petit il est possible de tronquer le développement tout en obtenant une bonne précision<ref group="N">Par exemple pour ''e < 0,1'' un développement à l'ordre 6 permettra d'obtenir une précision de 10<sup>-6</sup> sur ''E'', cf. Duriez, ''op. cit.'', pages 171-172.</ref>.

== Sources ==
* Danjon : cours de cosmographie de [[mathématiques]] élémentaires, 1950
* Maillard et Millet : cours de cosmographie de mathématiques élémentaires, 1956
* Wintner : analytical foundations of celestial mechanics, 1941, Princeton
* {{Ouvrage|langue=en|prénom1=Isaac|nom1=Newton|auteur2=Dana Densmore|traducteur=William H. Donahue|titre=Newton's Principia : the central argument : translation, notes, and expanded proofs|lieu=Santa Fe, N.M|éditeur=Green Lion Press|année=1995|pages totales=425|isbn=978-1-888009-01-9|isbn2=978-1-888-00900-2|oclc=33415300}}
* Brackenridge, Key to Newton's dynamics, 1995, u California p, {{ISBN|0-520-20065-9}}
* {{Ouvrage|langue=en|prénom1=Niccolò|nom1=Guicciardini|lien auteur1=Niccolò Guicciardini|titre=Reading the Principia : the debate on Newton's mathematical methods for natural philosophy from 1687 to 1736|lieu=Cambridge New York|éditeur=[[Cambridge University Press]]|année=1999|pages totales=285|isbn=978-0-521-64066-4|oclc=954623586}}
* Cordani, the Kepler problem, 2003, ed Birkhauser, {{ISBN|3-7643-6902-7}}
* J. Keill, philosophical transactions, 26/1 (1708), 174.
* Hermann, Histoires de l'Académie Royale des Sciences, paris (1710), 519-544.
* Leibniz, Tentamen de Motuum, [[Acta Eruditorum]] (1689), 82-96
* Varignon, mémoires de l'Ac Roy de Paris, (1700), 280.

== Notes et références ==
=== Notes ===
{{Références|group=N|colonnes=2}}

=== Références ===
{{Références|colonnes=2}}

== Voir aussi ==

=== Articles connexes ===
{{Colonnes|nombre = 2|
* [[Accélération de Siacci]]
* [[Balistique extérieure]]
* [[Diagramme horaire#Mouvement de Kepler selon Leibniz|Mouvement de Kepler selon Leibniz]]
* [[Échelle de temps en mécanique classique]]
* [[Histoire du mouvement képlérien]]
* [[Lois de Kepler]]
* [[Lois de Kepler, démonstration]]
* [[Mécanique spatiale]]
* [[Perturbation du mouvement képlérien]]
* [[Philosophiae Naturalis Principia Mathematica]]
* [[Principia et Calculus]]
* [[Problème à deux corps]]
* [[Résolution de l'équation de Kepler]]
* [[Satellite artificiel]]
* [[Théorème de Newton-Hamilton]]
* [[Vecteur de Runge-Lenz]]
}}

=== Liens externes ===

* [http://florenaud.free.fr/Kepler.php Définition et démonstration des lois de Kepler, sur le site Astrophy (en français)]
* [http://www.perso.ch/bernard.gisin/physique/KeplerLois/Cours_mecanique_3e_4e_Annexe_II.pdf Démonstration des trois lois de Kepler et propriétés d'une ellipse]

{{autres projets
  | wikiversity = L'univers en mouvement et le temps
  | wikiversity titre = L'univers en mouvement et le temps
}}
{{Palette|Johannes Kepler|Orbites}}
{{Portail|astronomie|physique}}

[[Catégorie:Mécanique céleste]]
[[Catégorie:Loi en physique|Mouvement Kepler]]
[[Catégorie:Johannes Kepler]]