Ellipse

PD Image
Fig. 1. Ellipse (black closed curve). The sum of the lengths of the two red line segments (ending in P₁) is equal to the sum of the lengths of the two blue line segments (ending in P₂).

In mathematics, an ellipse is a planar locus of points characterized by having a constant sum of distances to two given fixed points in the plane. In figure 1 two fixed points F₁ and F₂ are shown, these are the foci of the ellipse. An arbitrary point P₁ on the ellipse has distance F₁P₁ to F₁ and distance F₂P₁ to F₂. Let d be the sum of distances of P₁ to the foci,

$d\equiv \mathrm {F} _{1}\mathrm {P} _{1}+\mathrm {F} _{2}\mathrm {P} _{1},$

then all points of the ellipse have the constant sum of distances d. Thus, another arbitrary point P₂ on the ellipse has distance F₁P₂ to F₁ and distance F₂P₂ to F₂. By definition the sum of distances of P₂ to the foci is equal to d,

$\mathrm {F} _{1}\mathrm {P} _{2}+\mathrm {F} _{2}\mathrm {P} _{2}=d.\,$

The horizontal line segment S₁–S₂ in figure 1, going through the foci, is known as the major axis of the ellipse. Traditionally, the length of the major axis is indicated by 2a. The vertical dashed line segment, drawn halfway the foci perpendicular to the major axis, is referred to as the minor axis of the ellipse; its length is usually indicated by 2b. The major and minor axes are distinguished by a ≥ b. When a = b the ellipse is a circle—a special case of an ellipse. The foci of a circle coincide with the center of the circle. Clearly both ellipse axes are symmetry axes, reflection in either of them transforms the ellipse into itself. Basically, this a consequence of the fact that reflection conserves (sums of) distances. The intersection of the axes is the center of the ellipse.

The two foci and the points S₁ and S₂ are connected by reflection in the minor axis. Hence the distance S₂F₂ ≡ p is by symmetry equal to the distance S₁F₁. The distance of S₂ to F₁ is equal to 2a − p. By the definition of the ellipse the sum is equal to d, hence

$d=2a-p+p=2a.\;$

The sum d of distances from any point on the ellipse to the foci is equal to the length of the major axis.

Conic section

PD Image
Fig. 2. Upper shaded (green) section: ellipse; lower shaded (red) section: circle.

In the work of the Greek mathematician Apollonius (c. 262–190 BC) the ellipse arose as the intersection of a plane with a cone. Apollonius gave the ellipse its name, though the term ???????? (elleipsis, meaning "falling short") was used earlier by Euclid (c. 300 BC) in the construction of parallelograms with areas that "fell short". Apollonius applied the word to the conic section that at present we call ellipse. See Ref. for the—in modern eyes—complicated reasoning by which Apollonius tied the shape of certain conic sections to Euclid's concept of deficient areas.

In figure 2 a cone with a circular base is shown. It has a vertical symmetry axis, an axis of revolution. A cone can be generated by revolving around the axis a line that intersects the symmetry axis under an angle. A horizontal intersecting plane (plane perpendicular to the symmetry axis of the cone) gives a circle (a special ellipse), that is, the intersection of a horizontal plane with the cone is a circle. Planes that make an angle less than, or equal to, 90^° (but more than half the top angle of the cone) with the axis have an ellipse as intersection.

Eccentricity

The eccentricity e of an ellipse (usually denoted by e or ε) is the ratio of the distance OF₂ (cf. figure 3) to the length a (the semi-major axis), that is, e ≡ OF₂ / a. Let ${\vec {a}}$ be a vector of length a along the x-axis, then

$e{\vec {a}}\equiv {\overrightarrow {\mathrm {OF} }}_{2}.$

The following two vectors have common endpoint at P, see figure 3,

${\vec {r}}\equiv {\overrightarrow {\mathrm {OP} }}={\begin{pmatrix}x\\y\end{pmatrix}}\quad {\hbox{and}}\quad {\vec {g}}\equiv {\overrightarrow {\mathrm {F} _{2}\mathrm {P} }}={\overrightarrow {\mathrm {F} _{2}\mathrm {O} }}+{\overrightarrow {\mathrm {O} \mathrm {P} }}=-e{\vec {a}}+{\vec {r}}.$

Move P now to the positive y-axis; its new position vector is:

${\vec {r}}={\begin{pmatrix}0\\b\end{pmatrix}}.$

By symmetry, the distance of the moved P to either focus is equal to the semi-major axis a and equal to the length of the new vector ${\vec {g}}$ (with endpoint on the y-axis). For the following two inner products (indicated by a centered dot) we find,

${\vec {r}}\cdot (e{\vec {a}})=0\quad {\hbox{and}}\quad {\vec {r}}\cdot {\vec {r}}=b^{2}.$

PD Image
Fig. 3. An ellipse situated such that the major and minor axis are along Cartesian axes. The center of the ellipse coincides with the origin O.

Hence, (in fact by the Pythagoras theorem applicable for P on the y-axis),

$|{\vec {g}}\;|^{2}=a^{2}=({\vec {r}}-e{\vec {a}}\;)\cdot ({\vec {r}}-e{\vec {a}}\;)=b^{2}+e^{2}a^{2},$

so that the eccentricity is given by

$e={\sqrt {\frac {a^{2}-b^{2}}{a^{2}}}}\quad {\hbox{with}}\quad 0\leq e\leq 1.$

Algebraic form

Consider an ellipse that is located with respect to a Cartesian frame as in figure 3 (major axis on x-axis, minor axis on y-axis). For a point P=(x,y) of the ellipse it holds that

${\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1.$

Note that this equation is reminiscent of the equation for a unit circle. An ellipse may be seen as a unit circle in which the x and the y coordinates are scaled independently, by 1/a and 1/b, respectively.

Proof

Introduce the vectors

${\begin{aligned}{\overrightarrow {\mathrm {F} _{1}\mathrm {P} }}={\overrightarrow {\mathrm {F} _{1}\mathrm {O} }}+{\overrightarrow {\mathrm {O} \mathrm {P} }}=&\;e{\vec {a}}+{\vec {r}}\\{\overrightarrow {\mathrm {F} _{2}\mathrm {P} }}={\overrightarrow {\mathrm {F} _{2}\mathrm {O} }}+{\overrightarrow {\mathrm {O} \mathrm {P} }}=&-e{\vec {a}}+{\vec {r}}={\vec {g}}\\\end{aligned}}$

By definition of ellipse, the sum of the lengths is 2a

$|{\vec {r}}+e{\vec {a}}|+|{\vec {r}}-e{\vec {a}}|=2a\qquad \qquad \qquad \qquad (1)$

Multiply Eq. (1) by

$|{\vec {r}}+e{\vec {a}}|-|{\vec {r}}-e{\vec {a}}|$

and work out the left-hand side:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\vec{r} + e \vec{a}|^2 - |\vec{r} - e \vec{a}|^2 = 4e\vec{r}\cdot\vec{a} }

Hence

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 4e\vec{r}\cdot\vec{a} = 2a(|\vec{r} + e \vec{a}| - |\vec{r} - e \vec{a}|) }

Use

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\vec{r}\cdot\vec{a}}{a} = x }

and one obtains

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\vec{r} + e \vec{a}| - |\vec{r} - e \vec{a}| = 2ex \qquad\qquad\qquad\qquad(2) }

Add and subtract Eqs (1) and (2) and we find expressions for the distance of P to the foci,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} |\overrightarrow{\mathrm{F}_1\mathrm{P}}| &= |\vec{r} + e \vec{a}| &= a + ex &\\ |\overrightarrow{\mathrm{F}_2\mathrm{P}}| &= |\vec{r} - e \vec{a}| &= a - ex & \qquad\qquad\qquad\quad(3)\\ \end{align} }

Square both equations

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} r^2 +e^2a^2 +2 e\vec{r}\cdot\vec{a} &= a^2 + e^2x^2 + 2eax\\ r^2 +e^2a^2 -2 e\vec{r}\cdot\vec{a} &= a^2 + e^2x^2 - 2eax\\ \end{align} }

Adding, using the earlier derived value for e², and reworking gives

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r^2 +e^2a^2 = a^2 + e^2x^2 \Longrightarrow r^2 + a^2-b^2 = a^2 + \frac{x^2}{a^2} (a^2-b^2) \Longrightarrow x^2+y^2 = b^2 + x^2 - x^2\frac{b^2}{a^2}\Longrightarrow y^2 = b^2 - x^2\frac{b^2}{a^2} }

Division by b² gives finally

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{x^2}{a^2} + \frac{y^2}{b^2} = 1. }

Second degree equation

Under certain conditions the following general equation of second degree in x and y represents an ellipse:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(x,y) \equiv Ax^2 +2Bxy + Cy^2 + 2Dx + 2Ey + F = 0, }

(all variables are real).

The first condition is AC − B² > 0.

For the second condition one needs to solve a set of two linear equations yielding t₁ and t₂

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{pmatrix} A & B \\B & C\\ \end{pmatrix} \begin{pmatrix} t_1 \\t_2 \end{pmatrix} = -\begin{pmatrix} D\\E\\ \end{pmatrix}. }

Define σ ≡ f(t₁, t₂), then the second condition is σ ≠ 0.

The third condition is

If σ > 0 then A+C < 0

If σ < 0 then A+C > 0

Proof

In order to find the conditions we switch to matrix-vector notation and write f(x,y) as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}) = \mathbf{r}^\mathrm{T} \mathbb{A} \mathbf{r} + \mathbf{a}^\mathrm{T}\mathbf{r} +\mathbf{r}^\mathrm{T}\mathbf{a}+ F=0 }

with

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r} \equiv \begin{pmatrix} x\\y\\ \end{pmatrix}, \quad \mathbb{A} \equiv \begin{pmatrix} A & B \\B & C\\ \end{pmatrix},\quad \mathbf{a} \equiv \begin{pmatrix} D\\E\\ \end{pmatrix}. }

The superscript T stands for the transpose (row vector becomes column vector and vice versa). The expression can be rewritten by introducing the inverse of the matrix 𝔸. This gives the first condition: in order that the matrix be invertible its determinant det(𝔸) ≡ AC−B² ≠ 0.

PD Image
Fig. 4. r′ = r − t

Then

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}) =\left(\mathbf{r} + \mathbb{A}^{-1} \mathbf{a} \right)^\mathrm{T} \mathbb{A}\left(\mathbf{r} + \mathbb{A}^{-1} \mathbf{a} \right) - \mathbf{a}^\mathrm{T}\mathbb{A}^{-1}\mathbf{a} +F. }

Define

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{t} = - \mathbb{A}^{-1} \mathbf{a} \quad\Longrightarrow\quad \mathbb{A} \mathbf{t} = -\mathbf{a}. }

The minus sign in the definition of t is introduced to get the translation of the origin as depicted in figure 4.

Substitute t in the expression for f:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}') = \left(\mathbf{r}'\right)^\mathrm{T} \mathbb{A} \mathbf{r}' + f(\mathbf{t}) \quad\hbox{with}\quad f(\mathbf{t}) = -\mathbf{a}^\mathrm{T}\mathbb{A}^{-1} \mathbf{a}+ F=0, }

where we used

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{A}^\mathrm{T} = \mathbb{A}\quad\Longrightarrow \quad \left(\mathbb{A}^{-1}\right) ^\mathrm{T} = \mathbb{A}^{-1}, }

that is, the matrix 𝔸 and its inverse are symmetric.

By translation of the origin over t the linear terms in f have been eliminated, only two quadratic terms, one bilinear term, and one constant term appear in

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}') = \left(\mathbf{r}'\right)^\mathrm{T} \mathbb{A} \mathbf{r}' + f(\mathbf{t}). }

The price paid for it is the requirement det(𝔸) ≠ 0.

The next step is rotation of the x′ and y′ axis (with origin in O'), this will eliminate the bilinear term and decouple x′ and y′, the components of r′. Let us recall that a real symmetric matrix may be diagonalized by an orthogonal matrix. For the 2×2 case:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{R}^\mathrm{T} \mathbb{A} \mathbb{R} = \begin{pmatrix} \alpha_1 & 0 \\ 0 &\alpha_2\end{pmatrix}\quad \hbox{with}\quad \mathbb{R}^\mathrm{T}\mathbb{R} = \mathbb{R}\mathbb{R}^\mathrm{T}= \mathbb{E}, }

where the last matrix on the right is the identity matrix. Now

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}') = \left(\mathbf{r}'\right)^\mathrm{T}\; \mathbb{R}\mathbb{R}^\mathrm{T}\;\mathbb{A}\;\mathbb{R}\mathbb{R}^\mathrm{T}\; \mathbf{r}' + f(\mathbf{t}) = \left(\mathbf{r}''\right)^\mathrm{T} \boldsymbol{\alpha} \mathbf{r}'' + f(\mathbf{t}) }

with

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \boldsymbol{\alpha} \equiv \begin{pmatrix} \alpha_1 & 0 \\ 0 &\alpha_2\end{pmatrix},\quad\hbox{and}\quad \mathbf{r}'' \equiv \mathbb{R}^\mathrm{T}\mathbf{r}'. }

A (non-degenerate) ellipse is obtained if in the expression

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(\mathbf{r}'') =\left(\mathbf{r}''\right)^\mathrm{T} \boldsymbol{\alpha} \mathbf{r}'' + f(\mathbf{t}) = \alpha_1 (x'')^2 + \alpha_2 (y'')^2 + f(\mathbf{t}) = 0 }

the parameters α₁, α₂, and f(t) are non-zero and if the signs of α₁ and α₂ are equal and opposite to the sign of f(t). It is known that the determinant of a matrix is invariant under a matrix similarity transformation, hence equal signs of the α's gives

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \det(\mathbb{A}) = AC-B^2 = \alpha_1\alpha_2 > 0. }

So, the condition det(𝔸) ≠ 0 is narrowed to det(𝔸) > 0.

Also the trace A+C of the matrix is invariant under a similarity transformation. There are two possibilities in order that the equation f(x,y) = 0 represents an ellipse

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \alpha_1, \alpha_2 < 0& \Longrightarrow A+C=\alpha_1+\alpha_2 < 0 \quad\hbox{then}\quad f(\mathbf{t}) > 0 \\ \alpha_1, \alpha_2 > 0& \Longrightarrow A+C=\alpha_1+\alpha_2 > 0 \quad\hbox{then}\quad f(\mathbf{t}) < 0 .\\ \end{align} }

Clearly, it is necessary to solve t and determine the sign of f(t) in order to determine a priori whether the quadratic equation represents an ellipse, a diagonalization of 𝔸 is unnecessary.

Polar representation relative to focus

PD Image
Fig. 5. Polar representation

The length g of the vector (cf. figure 5) with endpoint on the ellipse

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overrightarrow{\mathrm{F}_2\mathrm{P}} \equiv \vec{g} = (ea + g\cos\theta,\; g\sin\theta) }

is given by the polar equation of an ellipse

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g = \frac{\ell}{1 + e\cos\theta}\quad \hbox{with}\quad \ell \equiv \frac{b^2}{a}, }

where 2l is known as the latus rectum (lit. right side) of the ellipse; it is equal to 2g for θ = 90^° (twice the length of the vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{g}} when it makes a right angle with the major axis).

Proof

Earlier [Eq. (3)] it was derived for the distance from the right focus F₂ to P that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\overrightarrow{\mathrm{F}_2\mathrm{P}}| \equiv |\vec{g}| \equiv g = a -ex . }

Elimination of x from

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = ea + g\cos\theta = ea + (a-ex)\cos\theta \, }

gives

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = \frac{ea + a\cos\theta}{1+e\cos\theta}, }

so that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g = a-ex = a - \frac{e^2 a + ea \cos\theta}{1+e\cos\theta} = \frac{a(1-e^2)}{1+e\cos\theta}. }

Substitute

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a(1-e^2) = a(1 - \frac{a^2-b^2}{a^2}) = \frac{b^2}{a} \equiv \ell }

and the polar equation for the ellipse follows.

Trammel construction

PD Image
Fig. 6. A trammel in theory

Before drafting was done almost exclusively by the aid of computers, draftsmen used a simple device for drawing ellipses, a trammel. Basically, a trammel is a rigid bar of length a (semi-major axis). In the top drawing of figure 6 the bar is shown as a blue line segment bounded by a black and a blue bead. On this bar a segment of length b (semi-minor axis) is marked; this is the red segment on the bar. Two beads fixed to the rigid bar move back and forth along the x-axis and y-axis, respectively. The blue bead fixed at one end of the bar moves along the y-axis, the red bead, which marks the beginning of the red segment of length b, moves along the x-axis. The endpoint of the bar (the black bead in figure 6) moves along an ellipse with semi-major axis a and semi-minor axis b and typically has a pen fixed to it.

(PD)  : http://chestofbooks.com/
Fig. 7. A trammel in practice

The fact that the trammel construction works is proved very easily, cf. the bottom drawing in figure 6,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = a\; \cos\theta \quad \hbox{and}\quad y = b\;\sin\theta . }

Hence

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{x^2}{a^2} + \frac{y^2}{b^2} = \cos^2\theta+\sin^2\theta = 1, }

which indeed is the equation for an ellipse.

A device called a trammel point is used to guide a woodworking router in making elliptical cuts.

Gardener's construction

(PD)  : http://chestofbooks.com/
Fig. 8. Gardener's construction

It is possible to construct an ellipse of given major and minor axis by the aid of a compass, a ruler, three thumbtacks, and a piece of string, see figure 8.

First draw the major axis AB, and then obtain with the compass its perpendicular bisector intersecting AB in the midpoint E. Along the bisector one measures off the length of the minor axis CD. Given that the distances CF and CG are the semi-major axis (AB/2), one can determine the foci by drawing an arc with the compass using C as center and AB/2 as radius. One now pins the thumbtacks in the foci and the point C and fixes a piece of string of length 2AB as shown in figure 8. Keeping the string taut one draws the upper part of the ellipse by moving the pencil from A to B. Then one removes the pin in point C and places it in D and repeats the procedure for the lower part of the ellipse.

Clearly this procedure can be used in the garden to create an elliptic lawn or flowerbed, which is why the procedure is sometimes referred to gardener's construction.

Notes

Figures 6 and 7 are from George Watson Kittredge, The New Metal Worker Pattern Book, David Williams Company, New York, (1901) Online