Ellipse (mathématiques)





Page d'aide sur l'homonymie Pour les articles homonymes, voir Ellipse.



L'ellipse est le lieu des points dont la somme des distances à deux points fixes, dits foyers, est constante.




Section du cône ou projection du cercle.


En géométrie, une ellipse est une courbe plane fermée obtenue par l’intersection d’un cône de révolution avec un plan, à condition que celui-ci coupe l'axe de rotation du cône ou du cylindre : c'est une conique d'excentricité strictement comprise entre 0 et 1. On peut également la définir comme le lieu des points dont la somme des distances à deux points fixes, dits foyers, est constante (sa construction par la méthode du jardinier est très simple).


Dans la vie courante, l’ellipse est la forme qu'on perçoit en regardant un cercle en perspective, ou la figure formée par l’ombre d'un disque sur une surface plane.


On retrouve aussi, en première approximation[1], des ellipses dans les trajectoires des corps célestes (planètes, comètes ou satellites artificiels) en orbite autour d'une étoile ou d’une autre planète. La Terre parcourt approximativement une ellipse dont le Soleil est un foyer.


Les différentes définitions de l'ellipse peuvent conduire, dans certains cas extrêmes à la construction d'un point, d'un segment ou d'un cercle, qui sont alors considérés comme des ellipses dégénérées n'en possédant pas toutes les propriétés géométriques.




Sommaire






  • 1 Définitions géométriques


    • 1.1 Section d’un cône


    • 1.2 Directrice et foyer


    • 1.3 Définition bifocale de l'ellipse


    • 1.4 Image d'un cercle par une affinité


    • 1.5 Construction par un cercle directeur


    • 1.6 Hypotrochoïde




  • 2 Propriétés géométriques


    • 2.1 Éléments de symétrie


    • 2.2 Tangente et bissectrice


    • 2.3 Cercle orthoptique




  • 3 Relations entre les grandeurs


  • 4 Équations caractéristiques


    • 4.1 Équation cartésienne


    • 4.2 Forme matricielle


    • 4.3 Équation paramétrique


    • 4.4 Équation polaire


    • 4.5 Paramètres d'une ellipse




  • 5 Circonférence


    • 5.1 Développements en séries


    • 5.2 Valeurs approchées




  • 6 Aire du domaine intérieur à une ellipse


  • 7 Tracer une ellipse


  • 8 Notes et références


  • 9 Bibliographie et sources


  • 10 Voir aussi


    • 10.1 Articles connexes


    • 10.2 Liens externes







Définitions géométriques |



Section d’un cône |




Construction d'une ellipse par section d'un cône.


L’ellipse est une courbe plane qui fait partie de la famille des coniques. Elle est obtenue par l’intersection d'un plan avec un cône de révolution (non dégénéré à une droite ou un plan) lorsque ce plan traverse de part en part le cône.


Le cercle est alors un cas particulier de l'ellipse (quand le plan de coupe est perpendiculaire à l'axe du cône, sans passer toutefois par son sommet).




Directrice et foyer |




Construction d'une ellipse par foyer et directrice. Ici avec une excentricité de 1/2.


Le cadre est l'espace affine euclidien de dimension 2. Soient (d) une droite, F un point n'appartenant pas à (d), e un réel dans ]0,1[. Soit P le plan affine déterminé par (d) et F. On appelle ellipse de droite directrice (d), de foyer F et d'excentricité e l'ensemble des points M du plan P vérifiant :


d⁡(M,F)d⁡(M,(d))=e{displaystyle {frac {operatorname {d} (mathrm {M} ,mathrm {F} )}{operatorname {d} (mathrm {M} ,(d))}}=e}{displaystyle {frac {operatorname {d} (mathrm {M} ,mathrm {F} )}{operatorname {d} (mathrm {M} ,(d))}}=e}

d(M, F) mesure la distance du point M au point F et d(M, (d)) = d(M, H) celle de M à la droite (d).


La constante e de proportionnalité des deux distances, appelée excentricité de l’ellipse, est sans dimension ; elle est caractéristique de la forme de l'ellipse, indépendamment de ses isométries (translations et/ou rotations) ou homothéties (de rapport non nul) dans le plan affine, et donc du choix arbitraire de repère orthonormé pour ce plan ; elle détermine tous les autres rapports de distances (et toutes les différences angulaires) mesurés sur l'ellipse.


Notons K le projeté orthogonal de F sur (d). La droite (KF) est alors clairement un axe de symétrie de l'ellipse appelé axe focal. Les points d'intersection de l'ellipse avec son axe focal sont appelés des sommets de l'ellipse. La médiatrice du segment dont les extrémités sont ces deux sommets est également un axe de symétrie de l'ellipse, appelé petit axe de l'ellipse.


Une telle ellipse a donc deux droites directrices, symétriques l'une de l'autre par rapport au petit axe de l’ellipse et deux foyers associés également symétriques par rapport au petit axe.


Une ellipse est également entièrement déterminée par la position de ses foyers et son excentricité, ou bien encore par la position de ses droites directrices (parallèles entre elles), et son excentricité.


En géométrie euclidienne, cette définition de l'ellipse exclut le cercle, le segment de droite et le point.



Définition bifocale de l'ellipse |




Construction de l'ellipse par deux foyers et une ficelle de longueur constante.


Soient F et F’ deux points distincts du plan. On appelle ellipse de foyers F et F’, l'ensemble des points M du plan vérifiant la propriété suivante :


d⁡(M,F)+d⁡(M,F′)=2a=2c2+b2,a∈R+∗,b∈R+∗{displaystyle operatorname {d} (mathrm {M} ,mathrm {F} )+operatorname {d} (mathrm {M} ,mathrm {F} ')=2a=2{sqrt {c^{2}+b^{2}}},quad ain mathbb {R} _{+}^{*},quad bin mathbb {R} _{+}^{*}}{displaystyle operatorname {d} (mathrm {M} ,mathrm {F} )+operatorname {d} (mathrm {M} ,mathrm {F} ')=2a=2{sqrt {c^{2}+b^{2}}},quad ain mathbb {R} _{+}^{*},quad bin mathbb {R} _{+}^{*}}

où 2a est la longueur du grand axe, 2c = d(F, F’) et 2b est la longueur du petit axe (perpendiculaire au grand axe). Cette relation exprime que la somme des distances d'un point M aux foyers est constante et vaut la longueur du grand axe.


Cette définition de l'ellipse permet de tracer un cercle de rayon a lorsque les deux foyers sont confondus mais on ne peut plus alors parler d'axe focal ni de petit axe. Elle permet aussi de tracer un segment de droite lorsque 2a = 2c.


Le théorème de Dandelin permet de placer géométriquement les deux foyers de l'ellipse obtenue par section d'un cône de révolution par un plan, à l'aide des deux sphères inscrites dans le cône et tangentes au plan de section. Un raisonnement analogue permet de prouver que la section d'un cylindre ou d'un hyperboloïde à une nappe par un plan coupant son axe est également une ellipse dont les foyers sont les points de tangence des sphères inscrites dans la surface de révolution et tangentes au plan.




Image d'un cercle par une affinité |




L'ellipse et les deux cercles de l'affinité.


Soient (C1) un cercle de centre O et de rayon b, (C2) un cercle de centre O et de rayon a (b < a), appelé cercle principal, et (xx’) une droite passant par O. On appelle ellipse de centre O, de demi-grand axe a et de demi-petit axe b l'image du cercle (C2) par l'affinité d'axe (xx’), de direction perpendiculaire à (xx’) et de rapport b/a.


Pour construire le point M de l'ellipse, image du point M2 du grand cercle, on construit le point M1 du cercle (C1) situé sur [OM2]. On mène par M2 une perpendiculaire à (xx’) et par M1 une parallèle à (xx’). Les droites se coupent en M. En effet, si m est le projeté orthogonal de M2 sur (xx’), on a, d'après le théorème de Thalès,


mMmM2=OM1OM2=ba{displaystyle {frac {mmathrm {M} }{mmathrm {M} _{2}}}={frac {mathrm {OM} _{1}}{mathrm {OM} _{2}}}={frac {b}{a}}}{displaystyle {frac {mmathrm {M} }{mmathrm {M} _{2}}}={frac {mathrm {OM} _{1}}{mathrm {OM} _{2}}}={frac {b}{a}}}

Un telle définition de l'ellipse permet de tracer un cercle lorsque le rapport de l'affinité est 1 ou un segment quand le rapport de l'affinité est 0.


Si A est un point de tangence entre l'ellipse et le grand cercle, l'angle AOM^1=AOM^2{displaystyle {widehat {mathrm {AOM} }}_{1}={widehat {mathrm {AOM} }}_{2}}{displaystyle {widehat {mathrm {AOM} }}_{1}={widehat {mathrm {AOM} }}_{2}} est appelé l'anomalie excentrique et peut servir de paramètre simple dans une équation paramétrée de l'ellipse.




Construction par un cercle directeur |




Ellipse et une portion de son cercle directeur.


Soient F et F’ deux points distincts, (C) un cercle de centre F’ et de rayon 2a (2a > FF’).


On appelle ellipse de cercle directeur (C) et de foyer F, l'ensemble des centres des cercles tangents intérieurement à (C) et passant par F.


Pour construire le point M sur l'ellipse, centre du cercle tangent à (C) en m, on trace la médiatrice du segment [Fm] : cette médiatrice rencontre le rayon [F’m] en M (le centre du cercle tangent) ; cette médiatrice est aussi une tangente en M de l’ellipse. On peut également obtenir l'ellipse par pliages : sur une feuille de papier, on dessine le cercle directeur et le foyer F et on plie la feuille de papier de manière à superposer un point du cercle avec le foyer F, l'ensemble des pliures obtenues dessine le faisceau des tangentes à l'ellipse.


La figure obtenue est symétrique par rapport au segment [FF’]. L'ellipse possède donc un autre cercle directeur de centre F et de rayon 2a.


Lorsque les deux points F et F’ sont confondus, une telle définition permet de tracer un cercle de centre F et de rayon a. Lorsque le rayon du cercle directeur est exactement égal à la distance FF’, la construction conduit au tracé du segment [FF’] : lorsque le point m est en F, les centres des cercles tangents intérieurement (C) décrivent le segment [FF’], lorsque m est distinct de F, le seul cercle tangent à (C) et passant par F est le cercle (C) lui-même et son centre est F’.


On peut également remarquer que le foyer F se projette orthogonalement sur la tangente en M en un point m’ appartenant au cercle principal de l'ellipse (le cercle de centre O et de rayon a). Plus précisément, l'ensemble des projetés de F sur les tangentes à l'ellipse dessine le cercle principal de l'ellipse[2], ce qui fait du cercle principal la podaire de l'ellipse par rapport à son foyer, et de l'ellipse l'anti-podaire du cercle principal par rapport au foyer.




Hypotrochoïde |


Articles détaillés : Hypotrochoïde et Théorème de La Hire.



Construction d'une ellipse comme hypotrochoïde : le cercle noir tourne à l'intérieur du cercle bleu, le point rouge décrit alors une ellipse.


Si on fait rouler sans glissement un cercle de centre O’ de rayon R à l'intérieur d'un cercle de centre O et de rayon 2R, le lieu parcouru par un point M, solidaire du petit cercle, situé à une distance d du centre O’ parcourt une ellipse. Son centre est O, ses demi-axes sont R + d et |R – d|. L'ellipse est donc un cas particulier d'hypotrochoïde.




Propriétés géométriques |



Éléments de symétrie |


« L’axe focal » (FF'), aussi appelé « axe principal », passant par les foyers et perpendiculaire aux directrices, est un axe de symétrie de l'ellipse ; il en est de même pour « l'axe secondaire », perpendiculaire à l'axe focal et passant par le « centre de l'ellipse » au milieu du « segment focal » [FF'].


L'intersection de l'axe focal et de l'axe secondaire, au centre de l'ellipse, est aussi un centre de symétrie.


Les points d'intersection de l'ellipse avec son axe focal sont appelés sommets principaux, ceux de l'ellipse avec son axe secondaire sont les sommets secondaires.


Le segment diamétral qui joint les sommets principaux (respectivement secondaires) en passant par le centre de l'ellipse est appelée « grand axe » (respectivement « petit axe »), comme aussi toute mesure de sa longueur une fois fixée l'unité arbitraire de longueur.



Tangente et bissectrice |




La bissectrice du secteur angulaire formé par les droites reliant un point de l'ellipse aux foyers est perpendiculaire à la tangente en ce point.


Soit une ellipse dont les foyers sont F et F'. En un point M de cette ellipse, considérons la bissectrice du secteur angulaire (FMF′){displaystyle (FMF')}(FMF'). Alors, cette bissectrice est perpendiculaire à la tangente en M.


Cette propriété, aussi appelée « propriété de réflexivité », est utilisée en optique géométrique dans les miroirs elliptiques : un rayon lumineux qui passe par un des foyers, lorsqu'il est réfléchi, passe par l'autre foyer. Ainsi, si l'on met une ampoule à un foyer d'un miroir elliptique, le faisceau lumineux se concentre sur l'autre foyer. Elle est également utilisée dans la construction de réflecteur de phare associé à une lentille convergente[3].




Chambre à écho, gravure de Friedrich Schultes dans Neue Hall- und Thon-Kunst d'Athanasius Kircher et Tobias Nislen (1684).


On la retrouve en acoustique : dans une pièce de forme elliptique, une personne placée à l'un des foyers peut entendre aisément ce que murmure une personne placée à l'autre. En effet, si la source d'un son se trouve à l’un des foyers, cette propriété assure que tous les sons réfléchis vont converger vers l'autre foyer, compensant la perte d'énergie sonore due à la distance parcourue, par la concentration vers le même point des multiples faisceaux réfléchis. Cependant, cette propriété de convergence spatiale ne suffit pas à assurer que le son perçu soit identique au son émis au foyer source car les signaux suivant différents trajets doivent pour cela arriver en phase à l'autre foyer. Si la nature du milieu est homogène (air ambiant par exemple), la vitesse de propagation du son y est constante, il faut alors pour que les signaux convergent en phase que les distances parcourues selon les différents trajets soient identiques. Or la somme des distances d'un point de l'ellipse aux foyers est constante pour tout point ce qui assure donc la convergence temporelle des signaux. C'est le cumul de ces deux propriétés de l'ellipse qui permet à la fois de comprendre à un foyer des paroles prononcées à l’autre (égale longueur des trajets maintenant les signaux en phase), et de les entendre même si elles ne sont que murmurées (convergence des trajets concentrant l’énergie des signaux). Cette propriété acoustique des voûtes elliptiques est connue au moins depuis le XVIIe siècle[4]. Elle serait à l'origine de la construction de chambres à écho[5], de confessionnaux pour lépreux[6] et expliquerait les qualités sonores des stations de métro parisiennes[6]. La rotonde du Capital Building à Washington et le Mormon Tabernacle à Salt Lake City sont des exemples de cette sorte de galeries[7].


Cette propriété est également utilisée en urologie dans les lithotriteurs extra-corporels : les ondes de choc émises à l'un des foyers de l'ellipse se concentrent, après réflexion sur un réflecteur ellipsoïde, sur l'autre foyer placé judicieusement sur le calcul à détruire[8].



.



Cercle orthoptique |




Ellipse et son cercle orthoptique.


L'ensemble des points M qui « regardent » l'ellipse suivant un angle droit, c'est-à-dire dont partent deux tangentes à l'ellipse perpendiculaires, dessine un cercle de même centre que l'ellipse et de rayon a2+b2{displaystyle {sqrt {a^{2}+b^{2}}}} sqrt{a^2+b^2}, appelé cercle orthoptique de l'ellipse[9].




Relations entre les grandeurs |




Une ellipse avec ses axes, son centre, un foyer et la directrice associée.


Les grandeurs (géométriques ou numériques) d’une ellipse sont :



  • la longueur du grand rayon (ou demi-grand axe), généralement notée a ;

  • la longueur du petit rayon (ou demi-petit axe), généralement notée b ;

  • la distance séparant le centre de l'ellipse et un des foyers, généralement notée c ;

  • la distance séparant un foyer F de sa directrice (d) associée, généralement notée h ;

  • la distance séparant le centre de l’ellipse et une de ses deux directrices, généralement notée f ;

  • l'excentricité de l’ellipse (strictement comprise entre 0 et 1), généralement notée e ;

  • le « paramètre » de l’ellipse, généralement noté p, représentant le demi latus rectum (corde parallèle à la directrice et passant par le foyer).


Des relations existent entre ces grandeurs :



  • si l’ellipse est définie par son excentricité e et la distance h entre le foyer F et la directrice (d), alors :
    a=eh1−e2,b=eh1−e2,c=e2h1−e2,f=h1−e2,p=eh;{displaystyle a={eh over 1-e^{2}},,b={eh over {sqrt {1-e^{2}}}},,c={e^{2}h over 1-e^{2}},,f={h over 1-e^{2}},,p=eh,;}a = {e h over 1 - e^2} ,, b = {e h over sqrt{1 - e^2}} ,, c = {e^2 h over 1 - e^2}  ,, f = {h over 1 - e^2} ,, p = e h ,;


  • si l’ellipse est donnée par ses rayons a{displaystyle a}a et b{displaystyle b}ba>b{displaystyle a>b}a > b :
    c=a2−b2,h=b2a2−b2,f=a2a2−b2,e=a2−b2a,p=b2a;{displaystyle c={sqrt {a^{2}-b^{2}}},,h={b^{2} over {sqrt {a^{2}-b^{2}}}},,f={a^{2} over {sqrt {a^{2}-b^{2}}}},,e={{sqrt {a^{2}-b^{2}}} over a},,p={b^{2} over a},;}c = sqrt{a^2 - b^2} ,, h = {b^2 over sqrt{a^2 - b^2}} ,, f = {a^2 over sqrt{a^2 - b^2}} ,, e = {sqrt{a^2 - b^2} over a} ,, p = {b^2 over a} ,;


  • lorsque l’on connait le grand rayon a{displaystyle a}a et l’excentricité e{displaystyle e}e :
    b=a1−e2,c=ae,h=a(1−e2)e,f=ae,p=a(1−e2).{displaystyle b=a{sqrt {1-e^{2}}},,c=ae,,h={a(1-e^{2}) over e},,f={a over e},,p=a(1-e^{2}),.}b = a sqrt{1 - e^2} ,, c = a e ,, h = {a(1 - e^2) over e} ,, f = {a over e} ,, p = a(1 - e^2) ,.


  • enfin, dans la définition bifocale de l'ellipse où sont connues la longueur 2a et la distance 2c entre les foyers :
    b=a2−c2,e=ca,h=a2−c2c,f=a2c,p=a2−c2a.{displaystyle b={sqrt {a^{2}-c^{2}}},,e={frac {c}{a}},,h={a^{2}-c^{2} over c},,f={a^{2} over c},,p={a^{2}-c^{2} over a},.}b =  sqrt{a^2 - c^2} ,, e = frac ca ,, h = {a^2-c^2 over c} ,, f = {a^2 over c} ,, p = {a^2-c^2 over a} ,.




Équations caractéristiques |



Équation cartésienne |


Dans le repère défini par le demi grand axe et le demi petit axe de l'ellipse, son équation est (si l'axe focal est x) :



(xa)2+(yb)2=1{displaystyle left({frac {x}{a}}right)^{2}+left({frac {y}{b}}right)^{2}=1}left ( frac{x}{a} right ) ^2 + left ( frac{y}{b} right ) ^2 = 1  ;

avec a > b > 0. La distance du centre de l'ellipse à l'un des foyers vaut


c=a2−b2{displaystyle c={sqrt {a^{2}-b^{2}}}}c = sqrt{a^2 - b^2}

et donc l'excentricité vaut


e=ca=1−(ba)2{displaystyle e={frac {c}{a}}={sqrt {1-left({frac {b}{a}}right)^{2}}}}{displaystyle e={frac {c}{a}}={sqrt {1-left({frac {b}{a}}right)^{2}}}}

Si l'axe focal est y alors a et b sont inversés.


Si une ellipse n'est pas centrée à l'origine d'un système de coordonnées, mais que son grand axe et son petit axe restent parallèles aux axes des coordonnées, celle-ci peut être spécifiée par l'équation suivante :


(x−u)2a2+(y−v)2b2=1{displaystyle {frac {(x-u)^{2}}{a^{2}}}+{frac {(y-v)^{2}}{b^{2}}}=1}frac{(x - u)^2}{a^2} + frac{(y - v)^2}{b^2} = 1

où les paramètres u et v sont les coordonnées du centre de l'ellipse.


Dans le cas où les axes de l'ellipse ne sont plus alignés avec les axes du plan mais tournés d'un angle θ par rapport au repère initial, on obtient l'équation cartésienne suivante :


((x−u)cos⁡θ+(y−v)sin⁡θ)2a2+((x−u)sin⁡θ(y−v)cos⁡θ)2b2=1{displaystyle {frac {left((x-u)cos theta +(y-v)sin theta right)^{2}}{a^{2}}}+{frac {left((x-u)sin theta -(y-v)cos theta right)^{2}}{b^{2}}}=1}{displaystyle {frac {left((x-u)cos theta +(y-v)sin theta right)^{2}}{a^{2}}}+{frac {left((x-u)sin theta -(y-v)cos theta right)^{2}}{b^{2}}}=1}

Comme toute conique, une ellipse possède une équation de la forme :


Ax2+Bxy+Cy2+Dx+Ey+F=0{displaystyle mathrm {A} x^{2}+mathrm {B} xy+mathrm {C} y^{2}+mathrm {D} x+mathrm {E} y+mathrm {F} =0}mathrm{A} x^2 + mathrm{B} xy + mathrm{C} y^2 + mathrm{D} x + mathrm{E} y + mathrm{F} = 0

avec la contrainte B² – 4AC < 0.
B étant réel, on en déduit immédiatement que A et C sont non nuls et de même signe (AC > (B/2)² ≥ 0).


Le coefficient B vaut 0 si les axes de l'ellipse sont parallèles à ceux de coordonnées. On a alors :


Ax2+Cy2+Dx+Ey+F=0{displaystyle mathrm {A} x^{2}+mathrm {C} y^{2}+mathrm {D} x+mathrm {E} y+mathrm {F} =0}mathrm{A} x^2 + mathrm{C} y^2 + mathrm{D} x + mathrm{E} y + mathrm{F} = 0

avec, à un coefficient de proportionnalité près :




  • A = 1/a2 ;


  • C = 1/b2 ;


  • D = –2u/a2 ;


  • E = –2v/b2 ;


  • F = u2/a2 + v2/b2 – 1.


Réciproquement, une équation de la forme


Ax2+Cy2+Dx+Ey+F=0{displaystyle mathrm {A} x^{2}+mathrm {C} y^{2}+mathrm {D} x+mathrm {E} y+mathrm {F} =0}mathrm{A} x^2 + mathrm{C} y^2 + mathrm{D} x + mathrm{E} y + mathrm{F} = 0

dans laquelle A, C et K=D24A+E24C−F{displaystyle mathrm {K} ={frac {mathrm {D} ^{2}}{4mathrm {A} }}+{frac {mathrm {E} ^{2}}{4mathrm {C} }}-mathrm {F} }mathrm{K} =frac{mathrm{D}^2}{4mathrm{A}} + frac{mathrm{E}^2}{4mathrm{C}} - mathrm{F} sont de même signe, est l'équation d'une ellipse pour laquelle




  • a=K/A{displaystyle a={sqrt {mathrm {K} /mathrm {A} }}}a = sqrt{mathrm {K}/mathrm{A}} ;


  • b=K/C{displaystyle b={sqrt {mathrm {K} /mathrm {C} }}}b = sqrt{mathrm{K}/mathrm{C}} ;


  • u=−D/2A{displaystyle u=-mathrm {D} /{2mathrm {A} }}u = -mathrm{D}/{2mathrm{A}} ;


  • v=−E/2C{displaystyle v=-mathrm {E} /{2mathrm {C} }}v = -mathrm{E}/{2mathrm{C}}.



Forme matricielle |


Nous pouvons exprimer l'équation cartésienne sous forme matricielle[10] :


txAx+tbx+c=0{displaystyle ^{mathrm {t} }mathbf {x} mathbf {A} mathbf {x} +^{mathrm {t} }mathbf {b} mathbf {x} +c=0}^{mathrm{t}}mathbf{x}mathbf{A}mathbf{x} + ^{mathrm{t}}mathbf{b}mathbf{x} + c = 0





  • x=(xy){displaystyle mathbf {x} ={begin{pmatrix}x\yend{pmatrix}}}mathbf{x} = begin{pmatrix} x \ y end{pmatrix} ; tx{displaystyle ^{mathrm {t} }mathbf {x} }^{mathrm{t}}mathbf{x} est la transposée de x ;


  • A{displaystyle mathbf {A} }mathbf{A} est une matrice symétrique réelle définie positive 2×2, A=(AB/2B/2C){displaystyle mathbf {A} ={begin{pmatrix}mathrm {A} &mathrm {B} /2\mathrm {B} /2&mathrm {C} end{pmatrix}}}mathbf{A} = begin{pmatrix} mathrm{A} & mathrm{B}/2 \ mathrm{B}/2 & mathrm{C} end{pmatrix} ;


  • b=(DE){displaystyle mathbf {b} ={begin{pmatrix}mathrm {D} \mathrm {E} end{pmatrix}}}mathbf{b} = begin{pmatrix} mathrm{D} \ mathrm{E} end{pmatrix} ; tb{displaystyle ^{mathrm {t} }mathbf {b} }^{mathrm{t}}mathbf{b} est la transposée de b{displaystyle mathbf {b} }mathbf {b}  ;


  • c= F.


On peut effectuer une rotation de matrice Q et une translation de matrice t, en effectuant un changement de coordonnées :


x=Qx¯+t{displaystyle mathbf {x} =mathbf {Q} {bar {mathbf {x} }}+mathbf {t} }mathbf{x} = mathbf{Q} bar{mathbf{x}} + mathbf{t}

l'équation devient :



tx¯(tQAQ)tx¯+(2ttA+tb)Qx¯+ttAt+tbt+c=0{displaystyle ^{mathrm {t} }{bar {mathbf {x} }}(^{mathrm {t} }mathbf {Q} mathbf {A} mathbf {Q} )^{mathrm {t} }{bar {mathbf {x} }}+(2^{mathrm {t} }mathbf {t} mathbf {A} +^{mathrm {t} }mathbf {b} )mathbf {Q} {bar {mathbf {x} }}+^{mathrm {t} }mathrm {t} mathbf {A} mathbf {t} +^{mathrm {t} }mathbf {b} mathbf {t} +c=0}^{mathrm{t}}bar{mathbf{x}}(^{mathrm{t}}mathbf{Q}mathbf{A}mathbf{Q})^{mathrm{t}}bar{mathbf{x}} + (2^{mathrm{t}}mathbf{t}mathbf{A} + ^{mathrm{t}}mathbf{b})mathbf{Q}bar{mathbf{x}} + ^{mathrm{t}}mathrm{t}mathbf{A}mathbf{t} + ^{mathrm{t}}mathbf{b}mathbf{t} + c = 0.

En posant :



=tQAQ{displaystyle {bar {mathbf {A} }}=^{mathrm {t} }mathbf {Q} mathbf {A} mathbf {Q} }bar{mathbf{A}} = ^{mathrm{t}}mathbf{Q}mathbf{A}mathbf{Q}

=(2ttA+tb)Q{displaystyle {bar {mathbf {b} }}=(2^{mathrm {t} }mathbf {t} mathbf {A} +^{mathrm {t} }mathbf {b} )mathbf {Q} }bar{mathbf{b}} = (2^{mathrm{t}}mathbf{t}mathbf{A} + ^{mathrm{t}}mathbf{b})mathbf{Q}

=ttAt+tbt+c{displaystyle {bar {c}}=^{mathrm {t} }mathbf {t} mathbf {A} mathbf {t} +^{mathrm {t} }mathbf {b} mathbf {t} +c}{displaystyle {bar {c}}=^{mathrm {t} }mathbf {t} mathbf {A} mathbf {t} +^{mathrm {t} }mathbf {b} mathbf {t} +c}


on obtient l'équation



tx¯+tb¯+c¯=0{displaystyle ^{mathrm {t} }{bar {mathbf {x} }}{bar {mathbf {A} }}{bar {mathbf {x} }}+^{mathrm {t} }{bar {mathbf {b} }}{bar {mathbf {x} }}+{bar {c}}=0}^{mathrm{t}}bar{mathbf{x}}bar{mathbf{A}}bar{mathbf{x}} + ^{mathrm{t}}bar{mathbf{b}}bar{mathbf{x}} + bar{c} = 0.

On peut alors choisir Q de sorte que {displaystyle {bar {mathbf {A} }}}bar{mathbf{A}} soit une matrice diagonale


=(λ100λ2){displaystyle {bar {mathbf {A} }}={begin{pmatrix}lambda _{1}&0\0&lambda _{2}end{pmatrix}}}bar{mathbf{A}} = begin{pmatrix} lambda_1 & 0 \ 0 & lambda_2end{pmatrix}

et t{displaystyle mathbf {t} }mathbf {t} de sorte que {displaystyle {bar {mathbf {b} }}}bar{mathbf{b}} soit nulle. L'équation cartésienne s'écrit alors



λ1x¯2+λ2y¯2+c¯=0{displaystyle lambda _{1}{bar {x}}^{2}+lambda _{2}{bar {y}}^{2}+{bar {c}}=0}lambda_1 bar{x}^2 + lambda_2 bar{y}^2 + bar{c} = 0.

Le centre de l'ellipse a alors pour coordonnées t, et les grand et petit axes valent :




a=−1{displaystyle a={sqrt {-{bar {c}}/lambda _{1}}}}a = sqrt{-bar{c}/lambda_1} ;


b=−2{displaystyle b={sqrt {-{bar {c}}/lambda _{2}}}}b = sqrt{-bar{c}/lambda_2}.




Équation paramétrique |


Dans un repère orthonormé du plan affine, dont les vecteurs directeurs sont parallèles aux axes de l'ellipse et où (xC, yC) sont les coordonnées du centre de l'ellipse, un paramétrage possible de l'ellipse est :


{x=xC+acos⁡ty=yC+bsin⁡tt∈R{displaystyle {begin{cases}x=x_{mathrm {C} }+acos t\y=y_{mathrm {C} }+bsin tend{cases}}quad tin mathbb {R} }{displaystyle {begin{cases}x=x_{mathrm {C} }+acos t\y=y_{mathrm {C} }+bsin tend{cases}}quad tin mathbb {R} }

Ce paramétrage est périodique, de période , ce qui signifie qu'on peut aussi le restreindre à un intervalle semi-ouvert de longueur minimale égale à cette période, par exemple à t∈[0,2π[{displaystyle tin [0,2pi [}t in [0, 2pi[. Toute restriction du paramètre à un intervalle fermé de longueur inférieure donne un arc fermé ne parcourant qu'une partie de l'ellipse.


Pour ce paramétrage, le rayon de courbure r au point M de paramètre t est donné par


r(t)=a2b(1−e2cos2⁡t)3/2t∈R{displaystyle r(t)={frac {a^{2}}{b}}left(1-e^{2}cos ^{2}tright)^{3/2}qquad tin mathbb {R} }r(t) = frac{a^2}{b} left(1-e^2 cos^2 tright)^{3/2} qquad t inR

D'où, en particulier, aux sommets :



  • r(0)=b2a{displaystyle r(0)={frac {b^{2}}{a}}} r(0) = frac{b^2}{a}

  • r(π/2)=a2b{displaystyle r(pi /2)={frac {a^{2}}{b}}} r(pi/2) = frac{a^2}{b}


Ces deux rayons de courbures permettent une construction approchée de l'ellipse à l'aide de 4 arcs de cercle[11].


Si l'on travaille dans le plan complexe, la représentation paramétrique de l'ellipse de centre C et de demi-axes a et b est


z=zC+a+b2eit+a−b2e−itt∈R{displaystyle z=z_{mathrm {C} }+{frac {a+b}{2}}mathrm {e} ^{mathrm {i} t}+{frac {a-b}{2}}mathrm {e} ^{-mathrm {i} t}quad tin mathbb {R} }{displaystyle z=z_{mathrm {C} }+{frac {a+b}{2}}mathrm {e} ^{mathrm {i} t}+{frac {a-b}{2}}mathrm {e} ^{-mathrm {i} t}quad tin mathbb {R} }

Si maintenant les axes de l'ellipse ne sont plus alignés avec les axes du plan mais tourné d'un angle θ, on obtient l'équation paramétrique en multipliant les coordonnées relatives au centre par la matrice de rotation :


(cos⁡θsin⁡θsin⁡θcos⁡θ(acos⁡tbsin⁡t)=(acos⁡θcos⁡t−bsin⁡θsin⁡tasin⁡θcos⁡t+bcos⁡θsin⁡t){displaystyle {begin{pmatrix}cos theta &-sin theta \sin theta &cos theta end{pmatrix}}times {begin{pmatrix}acos t\bsin tend{pmatrix}}={begin{pmatrix}acos theta cos t-bsin theta sin t\asin theta cos t+bcos theta sin tend{pmatrix}}}{displaystyle {begin{pmatrix}cos theta &-sin theta \sin theta &cos theta end{pmatrix}}times {begin{pmatrix}acos t\bsin tend{pmatrix}}={begin{pmatrix}acos theta cos t-bsin theta sin t\asin theta cos t+bcos theta sin tend{pmatrix}}}

et l'équation paramétrique devient


{x=xC+acos⁡θcos⁡t−bsin⁡θsin⁡ty=yC+asin⁡θcos⁡t+bcos⁡θsin⁡tt∈R{displaystyle {begin{cases}x=x_{mathrm {C} }+acos theta cos t-bsin theta sin t\y=y_{mathrm {C} }+asin theta cos t+bcos theta sin tend{cases}}quad tin mathbb {R} }{displaystyle {begin{cases}x=x_{mathrm {C} }+acos theta cos t-bsin theta sin t\y=y_{mathrm {C} }+asin theta cos t+bcos theta sin tend{cases}}quad tin mathbb {R} }

que l'on peut aussi noter[12] :


{x=xC+axcos⁡t+bxsin⁡ty=yC+aycos⁡t+bysin⁡tt∈R{displaystyle {begin{cases}x=x_{mathrm {C} }+a_{x}cos t+b_{x}sin t\y=y_{mathrm {C} }+a_{y}cos t+b_{y}sin tend{cases}}quad tin mathbb {R} }{displaystyle {begin{cases}x=x_{mathrm {C} }+a_{x}cos t+b_{x}sin t\y=y_{mathrm {C} }+a_{y}cos t+b_{y}sin tend{cases}}quad tin mathbb {R} }

avec




  • ax = a cos θ ;


  • ay = a sin θ ;


  • bx = –b sin θ ;


  • by = b cos θ.



Équation polaire |


Dans le repère défini par le foyer et l'axe focal, en prenant l'origine des angles du côté de la directrice, l'équation polaire de l'ellipse de demi-axes a et b est :


[3a]r1(θ)=p1+ecos⁡θθR.{displaystyle [3a]qquad r_{1}(theta )={frac {p}{1+ecos theta }}qquad theta in mathbb {R} .}[3a] qquad r_1(theta) = frac{p}{1+e cos theta} qquad theta inR.

Remarque : r1(θ)=p1−ecos⁡θ {displaystyle qquad r_{1}(theta )={frac {p}{1-ecos theta }}qquad ~}qquad r_1(theta) = frac{p}{1-e cos theta} qquad~ quand l'origine des angles est prise du côté opposé à la directrice (r1(0)=a+c{displaystyle r_{1}(0)=a+c}r_1(0) = a+c , c étant la demi-distance entre les foyers).


Rappel : le « paramètre » de l’ellipse, généralement noté p, représentant le demi latus rectum (corde parallèle à la directrice et passant par le foyer).


Dans le repère défini par le centre et l'axe focal, cette équation est :


[3b]r22(θ)=b21−e2cos2⁡θθR.{displaystyle [3b]qquad r_{2}^{2}(theta )={frac {b^{2}}{1-e^{2}cos ^{2}theta }}qquad theta in mathbb {R} .}[3b] qquad r_2^2(theta) = frac{b^2}{1-e^2 cos ^2 theta} qquad theta inR.


Paramètres d'une ellipse |


Le paramétrage d'une ellipse est un point important pour les algorithmes de régression ou de reconnaissance de forme. Une ellipse est décrite par cinq paramètres.































Paramétrisation de l'ellipse
Définition
Paramètres
Détails
Foyer, directrice
et excentricité

x, y, a, b, e
Coordonnées (x, y) de F, paramètres de la droite provenant de l'équation cartésienne — pente et ordonnées à l'origine (a, b) — ou bien de l'équation normale — (ρ, θ) —, et excentricité e.
Foyers et
grand axe

x1, y1, x2, y2, a
Coordonnées (x1, y1) de F1, (x2, y2) de F2 et longueur du demi grand axe a.
Équation cartésienne
B, C, D, E, F
ou bien
A, B, D, E, F
l'équation quadratique Ax2 + Bxy + Cy2 + Dx + Ey + F = 0 définit l'ellipse ; A et C sont non-nuls, on peut donc diviser l'équation par l'un des deux pour la normaliser : x2 + Bxy + Cy2 + Dx + Ey + F = 0 ou bien Ax2 + Bxy + y2 + Dx + Ey + F = 0.
Équation paramétrique

x, y, a, b, θ
ou bien
x, y, ax, ay, bx, by
Coordonnées (x, y) du centre C, longueurs des demi axes a et b, inclinaison θ du grand axe par rapport à l'axe des x. L'équation paramétrique peut aussi s'écrire avec ax = a cos θ, ay = bsin θ, bx = –bsin θ et by = bcos θ.


Circonférence |


Le calcul de la longueur d'un arc de l'ellipse de demi grand axe a, de demi petit axe b et d'excentricité e, compris dans le premier quadrant conduit au calcul de


L(B,M)=a∫0x/a1−e2u21−u2du{displaystyle L(B,M)=aint _{0}^{x/a}{frac {sqrt {1-e^{2}u^{2}}}{sqrt {1-u^{2}}}}mathrm {d} u}L(B,M)=a int_{0}^{x/a}frac{sqrt{1-e^2u^2}}{sqrt{1-u^2}}mathrm du

ou à celui de


L(B,M)=a∫1−e2sin2⁡tdt{displaystyle L(B,M)=aint _{0}^{theta }{sqrt {1-e^{2}sin ^{2}t}}mathrm {d} t}L(B,M)=a int_{0}^{theta}sqrt{1-e^2sin^2t}mathrm dt

pour la longueur de l'arc entre le point B (0, b) et le point M(x=asin⁡),y=bcos⁡)){displaystyle M(x=asin(theta ),y=bcos(theta ))}M(x=asin(theta) , y=bcos(theta))
Cependant, ces intégrales ne s'expriment pas à l'aide des fonctions classiques (fonctions algébriques, trigonométriques, logarithmes ou exponentielles).
Cette rectification, entreprise dès le début du XVIIIe siècle après la mise en place du calcul infinitésimal, conjointement à celle de la lemniscate de Bernoulli, a conduit à une nouvelle classe d'intégrales. Une telle intégrale est appelée intégrale elliptique incomplète de seconde espèce.


La circonférence de l'ellipse, c'est-à-dire son périmètre ou son orbite, s'exprime alors sous forme intégrale par


L=4a∫011−e2u21−u2du{displaystyle L=4aint _{0}^{1}{frac {sqrt {1-e^{2}u^{2}}}{sqrt {1-u^{2}}}}mathrm {d} u} L=4a int_{0}^{1}frac{sqrt{1-e^2u^2}}{sqrt{1-u^2}}mathrm du

ou encore


L=4a∫/21−e2sin2⁡tdt=4aE(e){displaystyle L=4aint _{0}^{pi /2}{sqrt {1-e^{2}sin ^{2}t}}mathrm {d} t=4a{rm {E(e)}}} L=4a int_{0}^{pi/2}sqrt{1-e^2sin^2t}mathrm dt = 4arm E(e)

E{displaystyle {rm {E}}}{rm E} est la fonction intégrale elliptique complète de deuxième espèce[13].


Il n'existe pas d'expression simple du résultat final mais des développements en série à convergence plus ou moins rapide dont les premiers sont l'œuvre de Colin Maclaurin de 1742[14]ou Euler en 1749[15] et 1773[16] et des approximations dont celles de Johannes Kepler de 1609[17], Giulio Fagnano de 1750[18] et de Ramanujan du début du XXe siècle.



Développements en séries |


Le premier développement en série est attribué à Colin Maclaurin mais on en trouve également une démonstration chez Euler en 1749[15]. Il s'obtient en utilisant le développement en série de (1−x)1/2{displaystyle (1-x)^{1/2}}(1-x)^{1/2} pour x plus petit que 1 puis en intégrant par récurrence des expressions comme 01u2n1−e2u2du{displaystyle int _{0}^{1}{frac {u^{2n}}{sqrt {1-e^{2}u^{2}}}}mathrm {d} u}int_0^1frac{u^{2n}}{sqrt{1-e^2u^2}}mathrm du pour Euler[15] ou bien en utilisant les intégrales de Wallis, si c'est la forme trigonométrique qui est privilégiée[19]. Il s'exprime en fonction de a, demi grand axe et e excentricité par :


L=2πa∑n=0+∞[(2n)!(2n⋅n!)2]2−e2n2n−1=2πa[1−(12)2e2−(1⋅32⋅4)2e43−(1⋅3⋅52⋅4⋅6)2e65−]=2πa[1−14e2−364e4−5256e6−17516384e8−44165536e10⋯]{displaystyle {begin{aligned}L&=2pi asum _{n=0}^{+infty }left[{dfrac {(2n)!}{(2^{n}cdot n!)^{2}}}right]^{2}{frac {-e^{2n}}{2n-1}}\&=2pi aleft[1-left({dfrac {1}{2}}right)^{2}e^{2}-left({dfrac {1cdot 3}{2cdot 4}}right)^{2}{dfrac {e^{4}}{3}}-left({dfrac {1cdot 3cdot 5}{2cdot 4cdot 6}}right)^{2}{dfrac {e^{6}}{5}}-cdots right]\&=2pi aleft[1-{dfrac {1}{4}}e^{2}-{dfrac {3}{64}}e^{4}-{dfrac {5}{256}}e^{6}-{dfrac {175}{16384}}e^{8}-{dfrac {441}{65536}}e^{10}cdots right]end{aligned}}}{begin{aligned}L&=2pi asum _{{n=0}}^{{+infty }}left[{dfrac  {(2n)!}{(2^{n}cdot n!)^{2}}}right]^{2}{frac  {-e^{{2n}}}{2n-1}}\&=2pi aleft[1-left({dfrac  {1}{2}}right)^{2}e^{2}-left({dfrac  {1cdot 3}{2cdot 4}}right)^{2}{dfrac  {e^{4}}{3}}-left({dfrac  {1cdot 3cdot 5}{2cdot 4cdot 6}}right)^{2}{dfrac  {e^{6}}{5}}-cdots right]\&=2pi aleft[1-{dfrac  {1}{4}}e^{2}-{dfrac  {3}{64}}e^{4}-{dfrac  {5}{256}}e^{6}-{dfrac  {175}{16384}}e^{8}-{dfrac  {441}{65536}}e^{{10}}cdots right]end{aligned}}

On peut aussi l'exprimer par :



L=2πa⋅2F1(−12,12;1;e2){displaystyle L=2pi acdot {}_{2}{rm {F}}_{1}left(-{frac {1}{2}},{frac {1}{2}};1;e^{2}right)}L = 2 pi a cdot {}_2{rm F}_1left(- frac12, frac12 ; 1 ; e^2 right),

2F1{displaystyle {}_{2}{rm {F}}_{1}}{}_2{rm F}_1 est la série hypergéométrique de Gauss.


Le principal défaut de la série ci-dessus est que sa convergence est très lente pour des excentricités proches de 1, c'est-à-dire pour des ellipses aplaties[20]. Elle reste toutefois utile sur le plan théorique car elle permet de comparer les autres formules donnant des approximations de la circonférence, puisqu'on peut toujours les développer sous forme de série de Taylor, et par différence déterminer la série de Taylor donnant les erreurs relatives commises par ces formules alternatives, en fonction de l'excentricité.


Le second développement intéressant est démontré par Sir James Ivory en 1796[21]. Il porte également le nom de série de Gauss-Kummer pour leur travail consistant à prouver l'équivalence entre plusieurs séries hypergéométriques[16]. C'est un développement symétrique en a et b s'exprimant à l'aide de la quantité h=(a−ba+b)2{displaystyle h=left({a-b over a+b}right)^{2}}h = left( {a - b over a + b} right)^2 :


L=π(a+b)⋅2F1(−12,−12;1;h)=π(a+b)∑n=0+∞(12n)2hn=π(a+b)∑n=0+∞[(2n)!(n!)21(1−2n)(−4)n]2hn=π(a+b)[1+h4+h264+h3256+25h416384+49h565536+⋯]{displaystyle {begin{aligned}L&=pi (a+b)cdot {}_{2}{rm {F}}_{1}left(-{frac {1}{2}},-{frac {1}{2}};1;hright)\&=pi (a+b)sum _{n=0}^{+infty }{{frac {1}{2}} choose n}^{2}h^{n}=pi (a+b)sum _{n=0}^{+infty }left[{{dfrac {(2n)!}{(n!)^{2}}}{frac {1}{(1-2n)(-4)^{n}}}}right]^{2}h^{n}\&=pi (a+b)left[1+{dfrac {h}{4}}+{dfrac {h^{2}}{64}}+{dfrac {h^{3}}{256}}+{dfrac {25h^{4}}{16384}}+{dfrac {49h^{5}}{65536}}+cdots right]end{aligned}}}begin{align}<br />
L & = pi (a + b) cdot{}_2{rm F}_1left(- frac12, -frac12 ; 1 ; h right)\<br />
 & = pi (a + b) sum_{n=0}^{+infty} {frac12 choose n}^2 h^n = pi (a + b) sum_{n=0}^{+infty} left[{ dfrac{(2n)!}{(n!)^2} frac{1}{(1 - 2n) (-4)^n} }right]^2 h^n \<br />
  & = pi (a + b) left[ 1 + dfrac{h}{4} + dfrac{h^2}{64} + dfrac{h^3}{256} + dfrac{25h^4}{16384} + dfrac{49h^5}{65536} + cdots right]<br />
end{align}

Cette série a meilleure convergence que celle de Maclaurin sur la totalité de la plage des excentricités mais l'erreur relative pour les ellipses aplaties laisse encore à désirer[22].


On peut également citer la série d'Euler[14] de 1773, développement symétrique en a et b s'exprimant à l'aide de la quantité d=(a2−b2a2+b2)2{displaystyle d=left({a^{2}-b^{2} over a^{2}+b^{2}}right)^{2}}d = left( {a^2 - b^2 over a^2 + b^2} right)^2.


L=π2(a2+b2)⋅2F1(−14,−14;1;d){displaystyle L=pi {sqrt {2(a^{2}+b^{2})}}cdot {}_{2}{rm {F}}_{1}left(-{frac {1}{4}},-{frac {1}{4}};1;dright)} L=pi sqrt{2(a^2+b^2)}cdot{}_2{rm F}_1left(- frac14, -frac14 ; 1 ; d right)

Cette série converge légèrement plus vite que la série de Maclaurin mais moins vite que celle de Gauss-Kummer[23].


Pour les ellipses à forte excentricité, la série d’Arthur Cayley (établie en 1876)[24] offre une meilleure convergence et s'exprime à l'aide de la quantité x=1−e2=(ba)2{displaystyle x=1-e^{2}=left({b over a}right)^{2}}x = 1 - e^2 = left({b over a}right)^2 par :


L=4a(1+∑n=1+∞xn2[(2n)!(2nn!)2]22n2n−1[ln⁡(16x)−(∑k=1n4(2k−1)2k)+2(2n−1)2n]){displaystyle L=4aleft(1+sum _{n=1}^{+infty }{x^{n} over 2}left[{(2n)! over (2^{n}n!)^{2}}right]^{2}{2n over 2n-1}left[ln left({16 over x}right)-left(sum _{k=1}^{n}{4 over (2k-1)2k}right)+{2 over (2n-1)2n}right]right)}L = 4a left( 1 + sum_{n = 1}^{+infty}<br />
                      {x^n over 2} left[{ (2n)! over (2^n n!)^2 }right]^2 { 2n over 2n - 1 }<br />
                      left[ lnleft({ 16 over x }right) - left(sum_{k = 1}^n{ 4 over (2k - 1)2k }right) + {2 over (2n - 1)2n} right]<br />
               right)<br />

qui se développe en :


L4a−1=x4[ln⁡(16x)−1]+3x232[ln⁡(16x)−136]+15x3256[ln⁡(16x)−125]+175x44096[ln⁡(16x)−1051420]+⋯ +xn212223242⋯(2n−3)2(2n−2)2(2n−1)2n[ln⁡(16x)−41×2−43×4−4(2n−3)(2n−2)−2(2n−1)(2n)]+⋯{displaystyle {begin{aligned}{L over 4a}-1&={frac {x}{4}}left[ln left({frac {16}{x}}right)-1right]+{frac {3x^{2}}{32}}left[ln left({frac {16}{x}}right)-{frac {13}{6}}right]+{15x^{3} over 256}left[ln left({16 over x}right)-{frac {12}{5}}right]+{175x^{4} over 4096}left[ln left({16 over x}right)-{frac {1051}{420}}right]+cdots \& +{x^{n} over 2}{1^{2} over 2^{2}}{3^{2} over 4^{2}}cdots {(2n-3)^{2} over (2n-2)^{2}}{(2n-1) over 2n}left[ln left({16 over x}right)-{4 over 1times 2}-{4 over 3times 4}-cdots -{4 over (2n-3)(2n-2)}-{2 over (2n-1)(2n)}right]+cdots end{aligned}}}{begin{aligned}{L over 4a}-1&={frac  {x}{4}}left[ln left({frac  {16}{x}}right)-1right]+{frac  {3x^{2}}{32}}left[ln left({frac  {16}{x}}right)-{frac  {13}{6}}right]+{15x^{3} over 256}left[ln left({16 over x}right)-{frac  {12}{5}}right]+{175x^{4} over 4096}left[ln left({16 over x}right)-{frac  {1051}{420}}right]+cdots \& +{x^{n} over 2}{1^{2} over 2^{2}}{3^{2} over 4^{2}}cdots {(2n-3)^{2} over (2n-2)^{2}}{(2n-1) over 2n}left[ln left({16 over x}right)-{4 over 1times 2}-{4 over 3times 4}-cdots -{4 over (2n-3)(2n-2)}-{2 over (2n-1)(2n)}right]+cdots end{aligned}}

Cette série exacte converge pour toute la plage des excentricités, mais n’a une convergence excellente et bien plus rapide que les séries de Maclaurin et de Gauss-Kummer que pour les excentricités élevées, où elle est également stable numériquement. Dans les évaluations numériques rapides, on peut donc associer cette série de Cayley à la série de Gauss-Kummer plus simple à calculer et convergeant plus rapidement pour les excentricités faibles et modérées, en fixant arbitrairement la coupure entre les deux méthodes à une valeur d’excentricité voisine de 96 % (qui correspond aussi à un rapport d’aspect b/a voisin de 28 % ou une valeur de x voisine de 7,84 %)[24].



Valeurs approchées |


Une première approximation très simple est présentée par Kepler en 1609[17] et consiste à approcher la circonférence de l'ellipse par la circonférence de deux cercles dont les rayons sont respectivement la moyenne géométrique et la moyenne arithmétique des quantités a et b :



ab≤L{displaystyle 2pi {sqrt {ab}}leq L} 2pi sqrt{ab} le L et L≈π(a+b){displaystyle Lapprox pi (a+b)} L approx  pi (a+b)

Une seconde approximation date de 1750 et provient de Giulio Carlo Fagnano[18]:


L≈π(32(a+b)−ab){displaystyle Lapprox pi left({frac {3}{2}}(a+b)-{sqrt {ab}}right)} L approx pileft(frac32 (a+b) - sqrt{ab}right)

Euler, quant à lui, propose[25] une approximation par excès d'une précision équivalente à celle par défaut obtenue en prenant la moyenne des deux demi-axes :


L≈π2(a2+b2){displaystyle Lapprox pi {sqrt {2(a^{2}+b^{2})}}} L approx pisqrt{2(a^2+b^2)}

On trouve également, dans des formulaires, cette approximation quadratique[20], moyenne quadratique des approximations de Kepler et d'Euler :


L≈π2(a2+b2)−12(a−b)2{displaystyle Lapprox pi {sqrt {2(a^{2}+b^{2})-{frac {1}{2}}(a-b)^{2}}}}L approx pi sqrt{ 2 (a^2 + b^2) - frac12 (a - b)^2 }

En exprimant b en fonction de e et en développant en série


L≈a34+141−e2−38e2=2πa[1−14e2−364e4−5256e6−898192e8−23132768e10⋯]{displaystyle {begin{matrix}Lapprox 2pi a{sqrt {{frac {3}{4}}+{frac {1}{4}}{sqrt {1-e^{2}}}-{frac {3}{8}}e^{2}}}=2pi aleft[1-{frac {1}{4}}e^{2}-{frac {3}{64}}e^{4}-{frac {5}{256}}e^{6}-{frac {89}{8192}}e^{8}-{frac {231}{32768}}e^{10}cdots right]end{matrix}}}begin{matrix}<br />
L approx 2 pi a sqrt{ frac 34 + frac 14sqrt{1 - e^2} - frac{3}{8}e^2 } = 2 pi a left[ 1 - frac{1}{4}e^2 - frac{3}{64}e^4 - frac{5}{256}e^6 - frac{89}{8192}e^8 - frac{231}{32768}e^{10} cdots right]<br />
end{matrix}

on obtient un développement dont les 4 premiers termes correspondent à la série de Maclaurin[26] ; l'erreur commise par les coefficients incorrects des termes de degré supérieur est alors très faible.


À partir du XIXe siècle les formules d'approximations se multiplient [27]
, chacun cherchant à trouver des formules simples ayant le maximum de précision. Une mention doit être faite pour les deux formules de Srinivasa Ramanujan de 1914




Première formule de Srinivasa Ramanujan[28]:


L≈π[3(a+b)−(3a+b)(a+3b)]=π(a+b)[3−4−h]{displaystyle Lapprox pi left[3(a+b)-{sqrt {(3a+b)(a+3b)}}right]=pi (a+b)left[3-{sqrt {4-h}}right]}L approx pi left[ 3(a + b) - sqrt{(3a + b)(a + 3b)} right] = pi (a + b) left[ 3 - sqrt{4 - h} right],

h=(a−ba+b)2{displaystyle h=left({frac {a-b}{a+b}}right)^{2}}h = left( frac{a - b}{a + b} right)^2.


Deuxième formule de Srinivasa Ramanujan[29] :

L≈π(a+b)[1+3h10+4−3h]{displaystyle Lapprox pi left(a+bright)left[1+{3h over 10+{sqrt {4-3h}}}right]}L approx pi left( a + b right) left[ 1 + {3h over 10 + sqrt{4 - 3h}} right]


Ces dernières formules sont très précises pour des excentricités modérées[30]. En revanche pour b = 0, la deuxième formule de Ramanujan conduit à une approximation de π de 22/7[31].



Aire du domaine intérieur à une ellipse |


Il existe différentes manières de calculer l'aire d'une ellipse. On peut se placer dans le repère porté par les axes où l'équation de l'ellipse s'écrit :


x2a2+y2b2=1{displaystyle {frac {x^{2}}{a^{2}}}+{frac {y^{2}}{b^{2}}}=1}frac{x^2}{a^2} + frac{y^2}{b^2} = 1

Avec les symétries établies plus haut, il suffit de calculer par exemple l'aire de la portion d'ellipse dans le quart supérieur droit du plan rapporté à ce repère. L'équation de la portion d'ellipse correspondante est :


y=b1−(xa)2{displaystyle y=b{sqrt {1-left({frac {x}{a}}right)^{2}}}}y= b sqrt{1 - left(frac xaright)^2}

pour x dans [0,a].
D'où l'aire du quart supérieur droit d'ellipse :


I=∫0ab1−(xa)2dx=ab∫011−t2dt=ab∫2cos2⁡udu{displaystyle I=int _{0}^{a}b{sqrt {1-left({frac {x}{a}}right)^{2}}},mathrm {d} x=abint _{0}^{1}{sqrt {1-t^{2}}},mathrm {d} t=abint _{0}^{frac {pi }{2}}cos ^{2}u,mathrm {d} u}I = int_0^a b sqrt{1-left(frac xaright)^2},mathrm dx = ab int_0^1 sqrt{1- t^2},mathrm dt = ab int_0^{fracpi2} cos^2 u,mathrm du

la dernière réécriture obtenue avec le changement de variable u↦sin⁡u=t{displaystyle umapsto sin u=t}u mapsto sin u = t de [0,π/2]{displaystyle [0,pi /2]}[0,pi/2] sur [0,1]{displaystyle [0,1]}[0,1]. Reste à linéariser cos2⁡u{displaystyle cos ^{2}u}cos^2 u pour trouver le quart de l'aire d'une ellipse :


I=ab∫21+cos⁡2u2du=πab4{displaystyle I=abint _{0}^{frac {pi }{2}}{frac {1+cos 2u}{2}},mathrm {d} u={frac {pi ab}{4}}}I= ab int_0^{fracpi2} frac{1+ cos 2u}2,mathrm du = frac{pi ab}4

et pour l'aire de toute l'ellipse :


S=πab{displaystyle S=pi ab} S= pi a b

Remarquer que pour a=b, on retrouve l'aire du disque.



Tracer une ellipse |


Article détaillé : Ellipsographe.

Méthode des deux points et de la corde : selon la définition bifocale, la somme AF+AF′{displaystyle AF+AF'}AF + AF' des distances entre un point A{displaystyle A}A de l'ellipse et ses deux foyers F et F' est constante. Ainsi, on plante deux piquets dans le sol (les deux foyers), on prend une corde non élastique de longueur donnée (la somme constante) que l'on attache aux piquets ; le trajet que l'on parcourt en maintenant la corde tendue est une ellipse. On nomme cette technique celle de « l'ellipse du jardinier ».




Ellipse servant à représenter un perçage droit vu en perspective (figure de droite) ; le trait d'axe vertical figure l'axe du perçage.


En dessin industriel, une ellipse est en général un cercle vu en perspective (une pièce est rarement elliptique même si ce n'est pas exclu), ou bien un perçage en biais par rapport à la surface de la pièce.




Ellipse servant à représenter un perçage oblique vu de face (figure de droite).


L'ellipse se représente donc avec les mêmes traits d'axe que pour le cercle. Dans le cas d'un cercle vu en perspective, ces traits d'axe sont inclinés et suivent les directions de référence. Dans le cas d'une forme réellement elliptique, les traits d'axes sont perpendiculaires.





Tracé d'une ellipse à main levée à l'aide d'un parallélogramme.


Tracé à main levée, méthode du parallélogramme exinscrit : on a vu ci-dessus qu'une ellipse pouvait être considérée comme un cercle vu en perspective. De même qu'un cercle est inscrit dans un carré, une ellipse est inscrite dans un parallélogramme qui n'est autre que ce carré vu en perspective cavalière (notez qu'il existe une infinité de parallélogrammes circonscrits, il suffit d'en choisir un). On trace d'abord un parallélogramme, on le divise en quatre quartiers selon les parallèles aux côtés passant par les milieux des autres côtés ; dans chaque quartier, on trace un arc passant par les milieux des côtés et tangent aux côtés en ces milieux (certaines caractéristiques de sécantes dans le cercle permettent de trouver d'autres points de passage intermédiaires de ces arcs).




Construction de l'ellipse par la méthode du parallélogramme.


On peut également tracer une ellipse en continu à l'aide d'un instrument mécanique appelé ellipsographe. Il en existe de plusieurs types exploitant chacun une propriété spécifique de l'ellipse.



Tracé à la règle, méthode du parallélogramme[32] : dans un rectangle de dimensions 2a et 4b, on découpe les côtés de longueurs 2a en 2n parties égales et les côtés de longueurs 4b, en 4n parties égales. On trace alors les segments comme sur la figure ci-contre. Leurs points d'intersection sont sur l'ellipse de demi-axes a et b. Le principe en est le suivant : une telle construction appliquée au cas où b = a conduit à tracer des segments perpendiculaires. Leur intersection est donc sur le cercle de diamètre 2a. La figure ci-contre est la déformation par une affinité de rapport b/a, de la figure qui conduit à la construction du cercle. Cette figure conduit donc à la construction de l'ellipse.




Notes et références |




  1. voir problème à deux corps et problème à N corps.


  2. Tauvel 2005, p. 384.


  3. Les différents projecteurs sur le site educauto.org.


  4. Voir par exemple le livre d'Athanasius Kircher, Musurgia universalis sive ars magna consoni et dissoni, 1650 tome II,p. 300.


  5. La chambre des échos de La Chaise Dieu, sur le site de l'Union des professeurs de physique et de chimie.


  6. a et bDr Bernard Auriol, Le confessionnal des lépreux.


  7. Swokowski (trad. Micheline Citta), Analyse, 5e édition, p. 621.


  8. Association Française d'Urologie, Qu'est-ce que la lithotritie extra corporelle (LEC) ?.


  9. Voir par exemple Tauvel 2005, p. 387, ou cet exercice corrigé sur la Wikiversité..


  10. (en) Walter Gander, Gene H. Golub et Rolf Strebel, « Least-Squares Fitting of Circles and Ellipses », BIT Numerical Mathematics, Springer, vol. 34, no 4,‎ décembre 1994, p. 558-578 (ISSN 0006-3835 et 1572-9125, lire en ligne).


  11. W Gellert, H. Küstner, M. Hellwich et H. Kästner, Petite encyclopédie des mathématiques, Didier, 1980, p. 197.


  12. (en) A. S. Aguado et M. S. Nixon, A new Hough transform mapping for ellipse detection, 1995(lire en ligne).


  13. Attention aux confusions sur la notation de la fonction E{displaystyle {rm {E}}}{rm E} ; par exemple dans le programme Mathematica et sur le site Wolfram, on trouve définie une fonction EllipticE{displaystyle {rm {EllipticE}}}{rm EllipticE} similaire, mais différente, EllipticE(x)=E(x){displaystyle {rm {EllipticE}}(x)={rm {E}}({sqrt {x}})}{rm EllipticE}(x) = {rm E}(sqrt x), de sorte que la formule à utiliser est alors L=4aEllipticE(e2){displaystyle L=4a,{rm {EllipticE}}(e^{2})}L = 4a , {rm EllipticE}(e^2).


  14. a et bAlmkvist et Berndt 1988, p. 14-598.


  15. a b et c(la) Euler, Animadversiones in rectificationem ellipsis, sur le site [Eulerarchive].


  16. a et b(en) Numericana, Exact Expansions for the Perimeter of an Ellipse.


  17. a et bAlmkvist et Berndt 1988, p. 15-599.


  18. a et bSerge Mehl, FAGNANO Giulio Carlo sur le site Chronomath.


  19. Serge Mehl, Longueur de l'ellipse sur le site Chronomath.


  20. a et b(en) Numericana, Circumference of an Ellipse.


  21. (en) James Ivory, A new séries for the rectification of the ellipsis (Almkvist et Berndt 1988, p. 14-598 note 25).


  22. (en) Numericana, Gauss-Kummer Series.


  23. (en) Numericana, Euler's Expansion.


  24. a et b(en)Numericana, Cayley's series.


  25. (en) Numericana, Euler's formula and the naive formula.


  26. (en) Numericana, Best quadratic formula.


  27. Almkvist et Berndt 1988, p. 16-600; 17-601.


  28. (en) Numericana, Ramanujan (I) & Lindner.


  29. (en) Numericana, Ramanujan (II).


  30. À titre d'exemple, Gerard Michon indique dans From Kepler's lower bound to Muir's lower bound, les erreurs pour le calcul d'un méridien terrestre (e proche de 0,082) de l'ordre de 1,75×10−12{displaystyle 1,75times 10^{-12}}1,75times 10^{-12} mètre pour la première formule et 1,63×10−25{displaystyle 1,63times 10^{-25}}1,63times 10^{-25} mètre pour la seconde.


  31. Almkvist et Berndt 1988, p. 18-602.


  32. (en) Computer Drafting, Rex Bookstore, Inc., Drawing ellipse using parallelogram method - p. 46.



Bibliographie et sources |



  • Patrice Tauvel, Géométrie : Agrégation-2e cycle/master, Dunod, 2005

  • Jean-Denis Eiden, Géométrie analytique classique, Calvage & Mounet, 2009, (ISBN 978-2-91-635208-4)

  • Méthodes modernes en géométrie de Jean Fresnel

  • Bruno Ingrao, Coniques affines, euclidiennes et projectives, C&M, (ISBN 978-2-916352-12-1)

  • (en) Gert Almkvist et Bruce Berndt, « Gauss, Landen, Ramanujan, the Arithmetic-Geometric Mean, Ellipses, Pi and the Ladies Diary », Amer. Math. Monthly, vol. 95, no 1,‎ août-septembre 1988(lire en ligne)


  • (en) Gerard P. Michon, Perimeter of an Ellipse - Final Answers, sur Numericana.com



Voir aussi |


.mw-parser-output .autres-projets ul{margin:0;padding:0}.mw-parser-output .autres-projets li{list-style-type:none;list-style-image:none;margin:0.2em 0;text-indent:0;padding-left:24px;min-height:20px;text-align:left}.mw-parser-output .autres-projets .titre{text-align:center;margin:0.2em 0}.mw-parser-output .autres-projets li a{font-style:italic}

Sur les autres projets Wikimedia :





Articles connexes |




  • Coniques

    • Hyperbole

    • Parabole



  • Régression elliptique

  • Transformée de Hough > Détection d'ellipses

  • Ovale

  • Théorème de Marden

  • Ellipsographe



Liens externes |



  • Bernard Gisin, « Propriétés de l'ellipse avec démonstrations » [PDF], août 2006

  • Sections d'un cône et théorème de Dandelin

  • Robert Ferréol, Jacques Mandonnet et Alain Esculier, « Ellipse », sur Encyclopédie des formes mathématiques remarquables, 2012


  • (en) Eric W. Weisstein, « Ellipse », MathWorld




  • Portail de la géométrie Portail de la géométrie



Popular posts from this blog

Quarter-circle Tiles

build a pushdown automaton that recognizes the reverse language of a given pushdown automaton?

Mont Emei