Ellipse (mathématiques)
Modèle:Infobox Polytope{\sqrt{1-u^2}}\mathrm du</math>
| propriétés =
}}
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<ref>voir problème à deux corps et problème à N corps.</ref>, 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.
Définitions géométriques
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
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 Modèle:Math. 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 :
- <math>\frac{\operatorname{d}(\mathrm{M}, \mathrm{F})}{\operatorname{d}(\mathrm{M}, (d))} = e </math>
où Modèle:Math mesure la distance du point M au point F et Modèle:Math 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
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 :
- <math>\operatorname{d}(\mathrm{M}, \mathrm{F}) + \operatorname{d}(\mathrm{M}, \mathrm{F}') = 2a = 2\sqrt{c^2 + b^2},\quad a \in\R^*_+,\quad b \in\R^*_+</math>
où Modèle:Math est la longueur du grand axe, Modèle:Math et Modèle:Math 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 Modèle:Math.
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é
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 Modèle:Sfrac.
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,
- <math>\frac{m\mathrm{M}}{m\mathrm{M}_2} = \frac{\mathrm{OM}_1}{ \mathrm{OM}_2 } = \frac{b}{a}</math>
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 <math>\widehat{\mathrm{AOM}}_1 = \widehat{\mathrm{AOM}}_2</math> 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
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<ref>Modèle:Harvsp.</ref>, 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
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 Modèle:Math et |Modèle:Math|. 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
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'). 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<ref>Les différents projecteurs sur le site educauto.org.</ref>.
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 Modèle:Lien siècleModèle:Vérification siècle<ref>Voir par exemple le livre d'Athanasius Kircher, [[Musurgia universalis|Modèle:Lang]], 1650 Modèle:Nobr,Modèle:P..</ref>. Elle serait à l'origine de la construction de chambres à écho<ref>La chambre des échos de La Chaise Dieu, sur le site de l'Union des professeurs de physique et de chimie.</ref>, de confessionnaux pour lépreux<ref name="Auriol">Modèle:Dr, Le confessionnal des lépreux.</ref> et expliquerait les qualités sonores des stations de métro parisiennes<ref name="Auriol" />. La rotonde du Capital Building à Washington et le Mormon Tabernacle à Salt Lake City sont des exemples de cette sorte de galeries<ref>Swokowski (trad. Micheline Citta), Analyse, Modèle:5e, Modèle:P..</ref>.
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<ref>Association Française d'Urologie, Qu'est-ce que la lithotritie extra corporelle (LEC) ?.</ref>.
Cercle orthoptique
Modèle:Article détaillé 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 <math> \sqrt{a^2+b^2}</math>, appelé cercle orthoptique, ou cercle de Monge de l'ellipse<ref>Voir par exemple Modèle:Harvsp, ou Modèle:Note autre projet.</ref>,<ref>Modèle:Lien web</ref>.
Relations entre les grandeurs
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 Modèle:Mvar ;
- la longueur du petit rayon (ou demi-petit axe), généralement notée Modèle:Mvar ;
- la distance séparant le centre de l'ellipse et un des foyers, généralement notée Modèle:Mvar ;
- la distance séparant un foyer F de sa directrice (d) associée, généralement notée Modèle:Mvar ;
- la distance séparant le centre de l’ellipse et une de ses deux directrices, généralement notée Modèle:Mvar ;
- l'excentricité de l’ellipse (strictement comprise entre 0 et 1), généralement notée Modèle:Mvar ;
- le « paramètre » de l’ellipse, généralement noté Modèle:Mvar, 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é Modèle:Mvar et la distance Modèle:Mvar entre le foyer F et la directrice (d), alors :
- <math>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 \,;</math>
- si l’ellipse est donnée par ses rayons Modèle:Mvar et Modèle:Mvar où Modèle:Math :
- <math>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} \,;</math>
- lorsque l’on connait le grand rayon Modèle:Mvar et l’excentricité Modèle:Mvar :
- <math>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) \,.</math>
- enfin, dans la définition bifocale de l'ellipse où sont connues la longueur Modèle:Math et la distance Modèle:Math entre les foyers :
- <math>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} \,.</math>
É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) :
- <math>\left ( \frac{x}{a} \right ) ^2 + \left ( \frac{y}{b} \right ) ^2 = 1 </math> ;
avec Modèle:Math. La distance du centre de l'ellipse à l'un des foyers vaut :
- <math>c = \sqrt{a^2 - b^2}</math>
et donc l'excentricité vaut :
- <math>e = \frac{c}{a} = \sqrt{1 - \left ( \frac{b}{a} \right ) ^2}</math>
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 :
- <math>\frac{(x - u)^2}{a^2} + \frac{(y - v)^2}{b^2} = 1</math>
où les paramètres Modèle:Mvar et Modèle:Mvar 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 :
- <math>\frac{[(x - u) \cos \theta + (y - v) \sin \theta]^2}{a^2} + \frac{[(x - u) \sin \theta - (y - v) \cos \theta]^2}{b^2} = 1</math>
Comme toute conique, une ellipse possède une équation de la forme :
- <math>\mathrm{A} x^2 + \mathrm{B} xy + \mathrm{C} y^2 + \mathrm{D} x + \mathrm{E} y + \mathrm{F} = 0</math>
avec la contrainte Modèle:Math. Modèle:Math étant réel, on en déduit immédiatement que Modèle:Math et Modèle:Math sont non nuls et de même signe (Modèle:Math).
On a alors la relation :
- <math>u = \frac{2CD-BE}{B^{2}-4AC}</math>
- <math>v = \frac{2AE-BD}{B^2-4AC}</math>
Le coefficient Modèle:Math vaut 0 si les axes de l'ellipse sont parallèles à ceux de coordonnées. On a alors :
- <math>\mathrm{A} x^2 + \mathrm{C} y^2 + \mathrm{D} x + \mathrm{E} y + \mathrm{F} = 0</math>
avec, à un coefficient de proportionnalité près :
Réciproquement, une équation de la forme
- <math>\mathrm{A} x^2 + \mathrm{C} y^2 + \mathrm{D} x + \mathrm{E} y + \mathrm{F} = 0</math>
dans laquelle Modèle:Math, Modèle:Math et <math>\mathrm{K} =\frac{\mathrm{D}^2}{4\mathrm{A}} + \frac{\mathrm{E}^2}{4\mathrm{C}} - \mathrm{F}</math> sont de même signe, est l'équation d'une ellipse pour laquelle
- <math>a = \sqrt{\mathrm {K}/\mathrm{A}}</math> ;
- <math>b = \sqrt{\mathrm{K}/\mathrm{C}}</math> ;
- <math>u = -\mathrm{D}/{2\mathrm{A}}</math> ;
- <math>v = -\mathrm{E}/{2\mathrm{C}}</math>.
Forme matricielle
Nous pouvons exprimer l'équation cartésienne sous forme matricielle<ref name="gander">Modèle:Article.</ref> :
- <math>^{\mathrm{t}}\mathbf{x}\mathbf{A}\mathbf{x} + ^{\mathrm{t}}\mathbf{b}\mathbf{x} + c = 0</math>
où
- <math>\mathbf{x} = \begin{pmatrix} x \\ y \end{pmatrix}</math> ; <math>^{\mathrm{t}}\mathbf{x}</math> est la transposée de Modèle:Math ;
- <math>\mathbf{A}</math> est une matrice symétrique réelle définie positive 2×2, <math>\mathbf{A} = \begin{pmatrix} \mathrm{A} & \mathrm{B}/2 \\ \mathrm{B}/2 & \mathrm{C} \end{pmatrix}</math> ;
- <math>\mathbf{b} = \begin{pmatrix} \mathrm{D} \\ \mathrm{E} \end{pmatrix}</math> ; <math>^{\mathrm{t}}\mathbf{b}</math> est la transposée de <math>\mathbf{b}</math> ;
- Modèle:Math.
On peut effectuer une rotation de matrice Modèle:Math et une translation de matrice Modèle:Math, en effectuant un changement de coordonnées :
- <math>\mathbf{x} = \mathbf{Q} \bar{\mathbf{x}} + \mathbf{t}</math>
l'équation devient :
- <math>^{\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</math>.
En posant :
- <math>\bar{\mathbf{A}} = ^{\mathrm{t}}\mathbf{Q}\mathbf{A}\mathbf{Q}</math>
- <math>\bar{\mathbf{b}} = (2^{\mathrm{t}}\mathbf{t}\mathbf{A} + ^{\mathrm{t}}\mathbf{b})\mathbf{Q}</math>
- <math>\bar{c} = ^{\mathrm{t}}\mathbf{t}\mathbf{A}\mathbf{t} + ^{\mathrm{t}}\mathbf{b}\mathbf{t} + c</math>
on obtient l'équation
- <math>^{\mathrm{t}}\bar{\mathbf{x}}\bar{\mathbf{A}}\bar{\mathbf{x}} + ^{\mathrm{t}}\bar{\mathbf{b}}\bar{\mathbf{x}} + \bar{c} = 0</math>.
On peut alors choisir Modèle:Math de sorte que <math>\bar{\mathbf{A}}</math> soit une matrice diagonale
- <math>\bar{\mathbf{A}} = \begin{pmatrix} \lambda_1 & 0 \\ 0 & \lambda_2\end{pmatrix}</math>
et <math>\mathbf{t}</math> de sorte que <math>\bar{\mathbf{b}}</math> soit nulle. L'équation cartésienne s'écrit alors
- <math>\lambda_1 \bar{x}^2 + \lambda_2 \bar{y}^2 + \bar{c} = 0</math>.
Le centre de l'ellipse a alors pour coordonnées Modèle:Math, et les grand et petit axes valent :
- <math>a = \sqrt{-\bar{c}/\lambda_1}</math> ;
- <math>b = \sqrt{-\bar{c}/\lambda_2}</math>.
Modèle:Démonstration(1,2) = \bar{\mathbf{A}}(2,1) = \frac{\mathrm{B}}{2}(\cos^2 \varphi - \sin^2 \varphi) + (\mathrm{C} - \mathrm{A}) \cos \varphi \sin \varphi</math> On désire déterminer φ de sorte que Modèle:Surligner(1, 2) = 0. En utilisant les formules de trigonométrie classiques on peut réécrire
- <math> \frac{\mathrm{B}}{2}\cos \left (2\varphi \right ) + \frac{\mathrm{C}-\mathrm{A}}{2}\sin \left (2\varphi \right ) = 0</math>
On en tire
- <math>\varphi=\frac{1}{2}\arctan\left ( \frac{\mathrm{B}}{\mathrm{A}-\mathrm{C}}\right)</math>
Détermination de t
La matrice Q n'est pas nulle. Donc, pour que Modèle:Surligner soit nulle, il faut et il suffit que
- <math>2^{\mathrm{t}}\mathbf{t}\mathbf{A} + ^{\mathrm{t}}\mathbf{b} = 0</math>.
Posons
- <math>\mathbf{t} = \begin{pmatrix} x_\mathrm{c} \\ y_\mathrm{c} \end{pmatrix}</math>
cette condition devient
- <math>\left \{ \begin{align}
2 \mathrm{A} x_\mathrm{c} + \mathrm{B} y_\mathrm{c} + \mathrm{D} = 0 \\ \mathrm{B} x_\mathrm{c} + 2 \mathrm{C} y_\mathrm{c} + \mathrm{E} = 0 \end{align} \right .</math> soit
- <math>\left \{ \begin{align}
x_\mathrm{c} = \frac{2 \mathrm{C} \mathrm{D} - \mathrm{B} \mathrm{E}}{\delta} \\ y_\mathrm{c} = \frac{2 \mathrm{A} \mathrm{E} - \mathrm{B} \mathrm{D}}{\delta} \end{align} \right .</math> avec δ = B2 - 4AC (non nul pour une ellipse). }}
É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 :
- <math>\begin{cases}x = x_\mathrm{C} + a\cos t \\
y = y_\mathrm{C} + b\sin t \end{cases} \quad t \in \R</math>
Ce paramétrage est périodique, de période Modèle:Math, ce qui signifie qu'on peut aussi le restreindre à un intervalle semi-ouvert de longueur minimale égale à cette période, par exemple à <math>t \in [0, 2\pi[</math>. 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
- <math>r(t) = \frac{a^2}{b} \left(1-e^2 \cos^2 t\right)^{3/2} \qquad t \in\R</math>
D'où, en particulier, aux sommets :
- <math> r(0) = \frac{b^2}{a} </math>
- <math> r(\pi/2) = \frac{a^2}{b} </math>
Ces deux rayons de courbures permettent une construction approchée de l'ellipse à l'aide de Modèle:Nombre de cercle<ref>Modèle:Ouvrage.</ref>.
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
- <math> z=z_\mathrm{C} + \frac{a + b}2\mathrm e^{\mathrm i t}+ \frac{a - b}2\mathrm e^{-\mathrm i t}\quad t \in \R</math>
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 :
- <math>\begin{pmatrix} \cos \theta & -\sin \theta \\
\sin \theta & \cos \theta \end{pmatrix} \times \begin{pmatrix} a\cos t \\ b\sin t \end{pmatrix} = \begin{pmatrix} a\cos \theta \cos t - b\sin \theta \sin t \\ a\sin \theta \cos t + b\cos \theta \sin t \end{pmatrix}</math> et l'équation paramétrique devient
- <math>\begin{cases}x = x_\mathrm{C} + a \cos \theta \cos t - b\sin \theta \sin t \\
y = y_\mathrm{C} + a \sin \theta \cos t + b\cos \theta \sin t \end{cases} \quad t \in \R</math> que l'on peut aussi noter<ref>Modèle:Ouvrage.</ref> :
- <math>\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 t \end{cases} \quad t \in \R</math> avec
- ax = a cos θ ;
- ay = a sin θ ;
- bx = –b sin θ ;
- by = b cos θ.
Équation polaire
Dans le repère défini par le foyer F et l'axe focal, avec l'origine des angles du côté opposé de la directrice, l'équation polaire de l'ellipse est :
- <math>\qquad r_1(\theta _1) = \frac{p}{1-e \cos \theta _1} = \frac{(a-c)(a+c)}{a-c \cos \theta _1} \qquad \theta \in\R.</math>
La Modèle:1re forme permet de retrouver la définition de p (le demi Modèle:Lang) avec <math>r _1\left({\pi \over 2}\right) = p</math>
La Modèle:2e forme permet de voir immédiatement
- <math>r_1(0) = a+c \quad \text{ et } \quad r_1(\pi) = a-c</math>
Remarque : <math>\qquad r_1(\theta _1) = \frac{p}{1+e \cos \theta _1}</math> quand l'origine des angles est prise du côté de la directrice.
Dans le repère défini par le centre 0 et l'axe focal, l'équation polaire est :
- <math>\qquad r_2^2(\theta _2) = \frac{b^2}{1-e^2 \cos ^2 \theta _2} = \frac{a^2b^2}{a^2\sin ^2\theta _2+b^2 \cos ^2 \theta _2} \qquad \theta \in\R.</math>
La Modèle:2e forme permet de voir immédiatement <math>r_2(0) =r_2(\pi) = a</math> et <math>r_2\left({\pi \over 2}\right) = r_2\left({3\pi \over 2}\right) = b</math>
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.
| 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 Modèle:Mvar, de demi petit axe Modèle:Mvar et d'excentricité Modèle:Mvar, compris dans le premier quadrant conduit au calcul de
- <math>L(B,M)=a \int_{0}^{x/a}\frac{\sqrt{1-e^2u^2}}{\sqrt{1-u^2}}\,\mathrm du</math>
ou à celui de
- <math>L(B,M)=a \int_{0}^{\theta}\sqrt{1-e^2\sin^2t}\,\mathrm dt</math>
pour la longueur de l'arc entre le point B (0, b) et le point <math>M(x=a\sin(\theta) , y=b\cos(\theta))</math>. 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 Modèle:Lien siècleModèle:Vérification 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
- <math> L=4a \int_{0}^{1}\frac{\sqrt{1-e^2u^2}}{\sqrt{1-u^2}}\mathrm du</math>
ou encore
- <math> L=4a \int_{0}^{\pi/2}\sqrt{1-e^2\sin^2t}\mathrm dt = 4a\rm E(e)</math>
où Modèle:Math est la fonction intégrale elliptique complète de deuxième espèce<ref>Attention aux confusions sur la notation de la fonction Modèle:Math ; par exemple dans le programme Mathematica et sur le site Wolfram, on trouve définie une fonction Modèle:Math similaire, mais différente, <math>{\rm EllipticE}(x) = {\rm E}(\sqrt x)</math>, de sorte que la formule à utiliser est alors <math>L = 4a \, {\rm EllipticE}(e^2)</math>.</ref>.
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<ref name="AB14">Modèle:Harvsp.</ref> ou Euler en 1749<ref name="Euler1749">{{#invoke:Langue|indicationDeLangue}} Euler, Animadversiones in rectificationem ellipsis, sur le site [Eulerarchive].</ref> et 1773<ref name="NumericanaExact">{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref> et des approximations dont celles de Johannes Kepler de 1609<ref name="AB15">Modèle:Harvsp.</ref>, Giulio Fagnano de 1750<ref name="Fagnano">Serge Mehl, FAGNANO Giulio Carlo sur le site Chronomath.</ref> et de Ramanujan du début du Modèle:Lien siècleModèle:Vérification 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<ref name="Euler1749" />. Il s'obtient en utilisant le développement en série de <math>(1-x)^{1/2}</math> pour Modèle:Mvar plus petit que 1 puis en intégrant par récurrence des expressions comme <math>\int_0^1\frac{u^{2n}}{\sqrt{1-e^2u^2}}\mathrm du</math> pour Euler<ref name="Euler1749" /> ou bien en utilisant les intégrales de Wallis, si c'est la forme trigonométrique qui est privilégiée<ref>Serge Mehl, Longueur de l'ellipse sur le site Chronomath.</ref>. Il s'exprime en fonction de Modèle:Mvar, demi grand axe et Modèle:Mvar excentricité par :
- <math>\begin{align}
L & = 2 \pi a \sum_{n=0}^{+\infty} \left[ \dfrac{ (2n)! }{ (2^n \cdot n!)^2 } \right]^2 \frac{-e^{2n}}{2n-1} \\
& = 2 \pi a \left[ 1 - \left(\dfrac{1}{2}\right)^2 e^2 - \left(\dfrac{1 \cdot 3}{2 \cdot 4}\right)^2 \dfrac{e^4}{3} - \left(\dfrac{1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6}\right)^2 \dfrac{e^6}{5} - \cdots \right] \\
& = 2 \pi a \left[ 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{align}</math>
On peut aussi l'exprimer par :
- <math>L = 2 \pi a \cdot {}_2{\rm F}_1\left(- \frac12, \frac12 ; 1 ; e^2 \right)</math>,
où <math>{}_2{\rm F}_1</math> 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<ref name="NumericaCircum">{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref>. 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<ref>{{#invoke:Langue|indicationDeLangue}} James Ivory, Modèle:Lang Modèle:Harv.</ref>. 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<ref name="NumericanaExact" />. C'est un développement symétrique en Modèle:Mvar et Modèle:Mvar s'exprimant à l'aide de la quantité <math>h = \left( {a - b \over a + b} \right)^2</math> :
- <math>\begin{align}
L & = \pi (a + b) \cdot{}_2{\rm F}_1\left(- \frac12, -\frac12 ; 1 ; h \right)\\
& = \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 \\
& = \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{align}</math>
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<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref>.
On peut également citer la série d'Euler<ref name="AB14" /> de 1773, développement symétrique en Modèle:Mvar et Modèle:Mvar s'exprimant à l'aide de la quantité <math>d = \left( {a^2 - b^2 \over a^2 + b^2} \right)^2</math>.
- <math> L=\pi \sqrt{2(a^2+b^2)}\cdot{}_2{\rm F}_1\left(- \frac14, -\frac14 ; 1 ; d \right)</math>
Cette série converge légèrement plus vite que la série de Maclaurin mais moins vite que celle de Gauss-Kummer<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref>.
Pour les ellipses à forte excentricité, la série d’Arthur Cayley (établie en 1876)<ref name="NumericaCayley">{{#invoke:Langue|indicationDeLangue}}Numericana, Modèle:Lang.</ref> offre une meilleure convergence et s'exprime à l'aide de la quantité <math>x = 1 - e^2 = \left({b \over a}\right)^2</math> par :
- <math>L = 4a \left( 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)
</math> qui se développe en :
- <math>\begin{align}
{L \over 4a} - 1 & = \frac{x}{4} \left[ \ln\left(\frac{16}{x}\right) - 1 \right]
+ \frac{3x^2}{32} \left[ \ln\left(\frac{ 16}{x }\right) - \frac{13}{ 6} \right]
+ { 15 x^3 \over 256} \left[ \ln\left({ 16 \over x }\right) - \frac{ 12}{ 5} \right]
+ {175 x^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 1\times 2}
- {4 \over 3\times 4}
- \cdots
- {4 \over (2n - 3)(2n - 2)}
- {2 \over (2n - 1)(2n)}
\right]
+ \cdots
\end{align}</math> 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 %)<ref name="NumericaCayley" />.
Valeurs approchées
Une première approximation très simple est présentée par Kepler en 1609<ref name="AB15" /> 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 Modèle:Mvar et Modèle:Mvar :
- <math> 2\pi \sqrt{ab} \le L \quad \mathrm{et}\quad L \approx \pi (a+b)</math>
Une seconde approximation date de 1750 et provient de Giulio Carlo Fagnano<ref name="Fagnano" />:
- <math> L \approx \pi\left(\frac32 (a+b) - \sqrt{ab}\right)</math>
Euler, quant à lui, propose<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref> une approximation par excès d'une précision équivalente à celle par défaut obtenue en prenant la moyenne des deux demi-axes :
- <math> L \approx \pi\sqrt{2(a^2+b^2)}</math>
On trouve également, dans des formulaires, cette approximation quadratique<ref name="NumericaCircum" />, moyenne quadratique des approximations de Kepler et d'Euler :
- <math>L \approx \pi \sqrt{ 2 (a^2 + b^2) - \frac12 (a - b)^2 }</math>
En exprimant Modèle:Mvar en fonction de Modèle:Mvar et en développant en série
- <math>\begin{matrix}
L \approx 2 \pi a \sqrt{ \frac 34 + \frac 14\sqrt{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] \end{matrix}</math> on obtient un développement dont les 4 premiers termes correspondent à la série de Maclaurin<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Modèle:Lang.</ref> ; l'erreur commise par les coefficients incorrects des termes de degré supérieur est alors très faible.
À partir du Modèle:Lien siècleModèle:Vérification siècle les formules d'approximations se multiplient <ref>Modèle:Harvsp.</ref> , 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<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Ramanujan (Modèle:Rom-maj) & Lindner.</ref>:
- <math>L \approx \pi \left[ 3(a + b) - \sqrt{(3a + b)(a + 3b)} \right] = \pi (a + b) \left[ 3 - \sqrt{4 - h} \right]</math>,
- où <math>h = \left( \frac{a - b}{a + b} \right)^2</math>.
- Deuxième formule de Srinivasa Ramanujan<ref>{{#invoke:Langue|indicationDeLangue}} Numericana, Ramanujan (Modèle:II).</ref> :
- <math>L \approx \pi \left( a + b \right) \left[ 1 + {3h \over 10 + \sqrt{4 - 3h}} \right]</math>
Ces dernières formules sont très précises pour des excentricités modérées<ref>À titre d'exemple, Gerard Michon indique dans Modèle:Lang, les erreurs pour le calcul d'un méridien terrestre (e proche de 0,082) de l'ordre de <math>1,75\times 10^{-12}</math> mètre pour la première formule et <math>1,63\times 10^{-25}</math> mètre pour la seconde.</ref>. En revanche pour b = 0, la deuxième formule de Ramanujan conduit à une approximation de π de 22/7<ref>Modèle:Harvsp.</ref>.
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 :
- <math>\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1</math>
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 :
- <math>y= b \sqrt{1 - \left(\frac xa\right)^2} </math>
pour x dans [0,a]. D'où l'aire du quart supérieur droit d'ellipse :
- <math>I = \int_0^a b \sqrt{1-\left(\frac xa\right)^2}\,\mathrm dx = ab \int_0^1 \sqrt{1- t^2}\,\mathrm dt = ab \int_0^{\frac\pi2} \cos^2 u\,\mathrm du</math>
la dernière réécriture obtenue avec le changement de variable <math>u \mapsto \sin u = t</math> de <math>[0,\pi/2]</math> sur <math>[0,1]</math>. Reste à linéariser <math>\cos^2 u</math> pour trouver le quart de l'aire d'une ellipse :
- <math>I= ab \int_0^{\frac\pi2} \frac{1+ \cos 2u}2\,\mathrm du = \frac{\pi ab}4</math>
et pour l'aire de toute l'ellipse :
- <math> S= \pi a b </math>
On peut remarquer que pour le cas Modèle:Math, on revient à un disque et la formule permet bien de retrouver l'aire du disque.
Tracer une ellipse
Méthode des deux points et de la corde : selon la définition bifocale, la somme <math>AF + AF'</math> des distances entre un point <math>A</math> 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 ».
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.
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é à 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).
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<ref>{{#invoke:Langue|indicationDeLangue}} Modèle:Lang, Modèle:Lang - Modèle:P..</ref> : 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ù Modèle:Formule 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
Voir aussi
Bibliographie et sources
- Modèle:Ouvrage
- Jean-Denis Eiden, Géométrie analytique classique, Calvage & Mounet, 2009, Modèle:ISBN
- Méthodes modernes en géométrie de Jean Fresnel
- Bruno Ingrao, Coniques affines, euclidiennes et projectives, C&M, Modèle:ISBN
- Modèle:Article
- {{#invoke:Langue|indicationDeLangue}} Gerard P. Michon, Modèle:Lang, sur Numericana.com
Articles connexes
- Coniques
- Régression elliptique
- Transformée de Hough > Détection d'ellipses
- Ovale
- Théorème de Marden
- Ellipsographe
- Ellipsoïde
