Ellipsoid: Difference between revisions

From testwiki
Jump to navigation Jump to search
imported>Fgnievinski
ย 
(No difference)

Latest revision as of 04:02, 14 February 2025

Template:Short description

Examples of ellipsoids with equation Template:Math: Template:Ubl

An ellipsoid is a surface that can be obtained from a sphere by deforming it by means of directional scalings, or more generally, of an affine transformation.

An ellipsoid is a quadric surface;โ€Š that is, a surface that may be defined as the zero set of a polynomial of degree two in three variables. Among quadric surfaces, an ellipsoid is characterized by either of the two following properties. Every planar cross section is either an ellipse, or is empty, or is reduced to a single point (this explains the name, meaning "ellipse-like"). It is bounded, which means that it may be enclosed in a sufficiently large sphere.

An ellipsoid has three pairwise perpendicular axes of symmetry which intersect at a center of symmetry, called the center of the ellipsoid. The line segments that are delimited on the axes of symmetry by the ellipsoid are called the principal axes, or simply axes of the ellipsoid. If the three axes have different lengths, the figure is a triaxial ellipsoid (rarely scalene ellipsoid), and the axes are uniquely defined.

If two of the axes have the same length, then the ellipsoid is an ellipsoid of revolution, also called a spheroid. In this case, the ellipsoid is invariant under a rotation around the third axis, and there are thus infinitely many ways of choosing the two perpendicular axes of the same length. In the case of two axes being the same length:

  • If the third axis is shorter, the ellipsoid is a sphere that has been flattened (called an oblate spheroid).
  • If the third axis is longer, it is a sphere that has been lengthened (called a prolate spheroid).

If the three axes have the same length, the ellipsoid is a sphere.

Standard equation

The general ellipsoid, also known as triaxial ellipsoid, is a quadratic surface which is defined in Cartesian coordinates as:

x2a2+y2b2+z2c2=1,

where a, b and c are the length of the semi-axes.

The points (a,0,0), (0,b,0) and (0,0,c) lie on the surface. The line segments from the origin to these points are called the principal semi-axes of the ellipsoid, because Template:Math are half the length of the principal axes. They correspond to the semi-major axis and semi-minor axis of an ellipse.

In spherical coordinate system for which (x,y,z)=(rsinθcosφ,rsinθsinφ,rcosθ), the general ellipsoid is defined as:

r2sin2θcos2φa2+r2sin2θsin2φb2+r2cos2θc2=1,

where θ is the polar angle and φ is the azimuthal angle.

When a=b=c, the ellipsoid is a sphere.

When a=bc, the ellipsoid is a spheroid or ellipsoid of revolution. In particular, if a=b>c, it is an oblate spheroid; if a=b<c, it is a prolate spheroid.

Parameterization

The ellipsoid may be parameterized in several ways, which are simpler to express when the ellipsoid axes coincide with coordinate axes. A common choice is

x=asinθcosφ,y=bsinθsinφ,z=ccosθ,

where

0θπ,0φ<2π.

These parameters may be interpreted as spherical coordinates, where Template:Mvar is the polar angle and Template:Mvar is the azimuth angle of the point Template:Math of the ellipsoid.[1]

Measuring from the equator rather than a pole,

x=acosθcosλ,y=bcosθsinλ,z=csinθ,

where

π2θπ2,0λ<2π,

Template:Mvar is the reduced latitude, parametric latitude, or eccentric anomaly and Template:Mvar is azimuth or longitude.

Measuring angles directly to the surface of the ellipsoid, not to the circumscribed sphere,

[xyz]=R[cosγcosλcosγsinλsinγ]

where

R=abcc2(b2cos2λ+a2sin2λ)cos2γ+a2b2sin2γ,π2γπ2,0λ<2π.

Template:Mvar would be geocentric latitude on the Earth, and Template:Mvar is longitude. These are true spherical coordinates with the origin at the center of the ellipsoid.Template:Citation needed

In geodesy, the geodetic latitude is most commonly used, as the angle between the vertical and the equatorial plane, defined for a biaxial ellipsoid. For a more general triaxial ellipsoid, see ellipsoidal latitude.

Volume

The volume bounded by the ellipsoid is

V=43πabc.

In terms of the principal diameters Template:Math (where Template:Math, Template:Math, Template:Math), the volume is

V=16πABC.

This equation reduces to that of the volume of a sphere when all three elliptic radii are equal, and to that of an oblate or prolate spheroid when two of them are equal.

The volume of an ellipsoid is Template:Sfrac the volume of a circumscribed elliptic cylinder, and Template:Sfrac the volume of the circumscribed box. The volumes of the inscribed and circumscribed boxes are respectively:

Vinscribed=833abc,Vcircumscribed=8abc.

Surface area

Template:See also

The surface area of a general (triaxial) ellipsoid is[2]

S=2πc2+2πabsin(φ)(E(φ,k)sin2(φ)+F(φ,k)cos2(φ)),

where

cos(φ)=ca,k2=a2(b2c2)b2(a2c2),abc,

and where Template:Math and Template:Math are incomplete elliptic integrals of the first and second kind respectively.[3]

The surface area of this general ellipsoid can also be expressed in terms of Template:Tmath, one of the Carlson symmetric forms of elliptic integrals:[4]

S=4πbcRG(a2b2,a2c2,1).

Simplifying the above formula using properties of Template:Math,[5] this can also be expressed in terms of the volume of the ellipsoid Template:Math:

S=3VRG(a2,b2,c2).

Unlike the expression with Template:Math and Template:Math, the equations in terms of Template:Math do not depend on the choice of an order on Template:Math, Template:Math, and Template:Math.

The surface area of an ellipsoid of revolution (or spheroid) may be expressed in terms of elementary functions:

Soblate=2πa2(1+c2ea2artanhe),where e2=1c2a2 and (c<a),

or

Soblate=2πa2(1+1e2eartanhe)

or

Soblate=2πa2 +πc2eln1+e1e

and

Sprolate=2πa2(1+caearcsine)where e2=1a2c2 and (c>a),

which, as follows from basic trigonometric identities, are equivalent expressions (i.e. the formula for Template:Math can be used to calculate the surface area of a prolate ellipsoid and vice versa). In both cases Template:Mvar may again be identified as the eccentricity of the ellipse formed by the cross section through the symmetry axis. (See ellipse). Derivations of these results may be found in standard sources, for example Mathworld.[6]

Approximate formula

S4πapbp+apcp+bpcp3p.

Here Template:Math yields a relative error of at most 1.061%;[7] a value of Template:Math is optimal for nearly spherical ellipsoids, with a relative error of at most 1.178%.

In the "flat" limit of Template:Mvar much smaller than Template:Mvar and Template:Mvar, the area is approximately Template:Math, equivalent to Template:Math.

Plane sections

Template:See also

Plane section of an ellipsoid

The intersection of a plane and a sphere is a circle (or is reduced to a single point, or is empty). Any ellipsoid is the image of the unit sphere under some affine transformation, and any plane is the image of some other plane under the same transformation. So, because affine transformations map circles to ellipses, the intersection of a plane with an ellipsoid is an ellipse or a single point, or is empty.[8] Obviously, spheroids contain circles. This is also true, but less obvious, for triaxial ellipsoids (see Circular section).

Determining the ellipse of a plane section

Plane section of an ellipsoid (see example)

Given: Ellipsoid Template:Math and the plane with equation Template:Math, which have an ellipse in common.

Wanted: Three vectors Template:Math (center) and Template:Math, Template:Math (conjugate vectors), such that the ellipse can be represented by the parametric equation

๐ฑ=๐Ÿ0+๐Ÿ1cost+๐Ÿ2sint

(see ellipse).

Plane section of the unit sphere (see example)

Solution: The scaling Template:Math transforms the ellipsoid onto the unit sphere Template:Math and the given plane onto the plane with equation

 nxau+nybv+nzcw=d.

Let Template:Math be the Hesse normal form of the new plane and

๐ฆ=[mumvmw]

its unit normal vector. Hence

๐ž0=δ๐ฆ

is the center of the intersection circle and

ρ=1δ2

its radius (see diagram).

Where Template:Math (i.e. the plane is horizontal), let

 ๐ž1=[ρ00],๐ž2=[0ρ0].

Where Template:Math, let

๐ž1=ρmu2+mv2[mvmu0],๐ž2=๐ฆ×๐ž1 .

In any case, the vectors Template:Math are orthogonal, parallel to the intersection plane and have length Template:Mvar (radius of the circle). Hence the intersection circle can be described by the parametric equation

๐ฎ=๐ž0+๐ž1cost+๐ž2sint.

The reverse scaling (see above) transforms the unit sphere back to the ellipsoid and the vectors Template:Math are mapped onto vectors Template:Math, which were wanted for the parametric representation of the intersection ellipse.

How to find the vertices and semi-axes of the ellipse is described in ellipse.

Example: The diagrams show an ellipsoid with the semi-axes Template:Math which is cut by the plane Template:Math. Template:Clear

Pins-and-string construction

Pins-and-string construction of an ellipse:
Template:Math, length of the string (red)
Pins-and-string construction of an ellipsoid, blue: focal conics
Determination of the semi axis of the ellipsoid

The pins-and-string construction of an ellipsoid is a transfer of the idea constructing an ellipse using two pins and a string (see diagram).

A pins-and-string construction of an ellipsoid of revolution is given by the pins-and-string construction of the rotated ellipse.

The construction of points of a triaxial ellipsoid is more complicated. First ideas are due to the Scottish physicist J. C. Maxwell (1868).[9] Main investigations and the extension to quadrics was done by the German mathematician O. Staude in 1882, 1886 and 1898.[10][11][12] The description of the pins-and-string construction of ellipsoids and hyperboloids is contained in the book Geometry and the imagination written by D. Hilbert & S. Vossen,[13] too.

Steps of the construction

  1. Choose an ellipse Template:Mvar and a hyperbola Template:Mvar, which are a pair of focal conics: E(φ)=(acosφ,bsinφ,0)H(ψ)=(ccoshψ,0,bsinhψ),c2=a2b2 with the vertices and foci of the ellipse S1=(a,0,0),F1=(c,0,0),F2=(c,0,0),S2=(a,0,0) and a string (in diagram red) of length Template:Mvar.
  2. Pin one end of the string to vertex Template:Math and the other to focus Template:Math. The string is kept tight at a point Template:Mvar with positive Template:Mvar- and Template:Mvar-coordinates, such that the string runs from Template:Math to Template:Mvar behind the upper part of the hyperbola (see diagram) and is free to slide on the hyperbola. The part of the string from Template:Mvar to Template:Math runs and slides in front of the ellipse. The string runs through that point of the hyperbola, for which the distance Template:Math over any hyperbola point is at a minimum. The analogous statement on the second part of the string and the ellipse has to be true, too.
  3. Then: Template:Mvar is a point of the ellipsoid with equation x2rx2+y2ry2+z2rz2=1rx=12(la+c),ry=rx2c2,rz=rx2a2.
  4. The remaining points of the ellipsoid can be constructed by suitable changes of the string at the focal conics.

Semi-axes

Equations for the semi-axes of the generated ellipsoid can be derived by special choices for point Template:Mvar:

Y=(0,ry,0),Z=(0,0,rz).

The lower part of the diagram shows that Template:Math and Template:Math are the foci of the ellipse in the Template:Mvar-plane, too. Hence, it is confocal to the given ellipse and the length of the string is Template:Math. Solving for Template:Mvar yields Template:Math; furthermore Template:Math.

From the upper diagram we see that Template:Math and Template:Math are the foci of the ellipse section of the ellipsoid in the Template:Mvar-plane and that Template:Math.

Converse

If, conversely, a triaxial ellipsoid is given by its equation, then from the equations in step 3 one can derive the parameters Template:Mvar, Template:Mvar, Template:Mvar for a pins-and-string construction.

Confocal ellipsoids

If Template:Overline is an ellipsoid confocal to Template:Mathcal with the squares of its semi-axes

rx2=rx2λ,ry2=ry2λ,rz2=rz2λ

then from the equations of Template:Mathcal

rx2ry2=c2,rx2rz2=a2,ry2rz2=a2c2=b2

one finds, that the corresponding focal conics used for the pins-and-string construction have the same semi-axes Template:Math as ellipsoid Template:Mathcal. Therefore (analogously to the foci of an ellipse) one considers the focal conics of a triaxial ellipsoid as the (infinite many) foci and calls them the focal curves of the ellipsoid.[14]

The converse statement is true, too: if one chooses a second string of length Template:Math and defines

λ=rx2rx2

then the equations

ry2=ry2λ,rz2=rz2λ

are valid, which means the two ellipsoids are confocal.

Limit case, ellipsoid of revolution

In case of Template:Math (a spheroid) one gets Template:Math and Template:Math, which means that the focal ellipse degenerates to a line segment and the focal hyperbola collapses to two infinite line segments on the Template:Mvar-axis. The ellipsoid is rotationally symmetric around the Template:Mvar-axis and

rx=12l,ry=rz=rx2c2.

Properties of the focal hyperbola

Top: 3-axial Ellipsoid with its focal hyperbola.
Bottom: parallel and central projection of the ellipsoid such that it looks like a sphere, i.e. its apparent shape is a circle
True curve
If one views an ellipsoid from an external point Template:Mvar of its focal hyperbola, then it seems to be a sphere, that is its apparent shape is a circle. Equivalently, the tangents of the ellipsoid containing point Template:Mvar are the lines of a circular cone, whose axis of rotation is the tangent line of the hyperbola at Template:Mvar.[15][16] If one allows the center Template:Mvar to disappear into infinity, one gets an orthogonal parallel projection with the corresponding asymptote of the focal hyperbola as its direction. The true curve of shape (tangent points) on the ellipsoid is not a circle.Template:Paragraph The lower part of the diagram shows on the left a parallel projection of an ellipsoid (with semi-axes 60, 40, 30) along an asymptote and on the right a central projection with center Template:Mvar and main point Template:Mvar on the tangent of the hyperbola at point Template:Mvar. (Template:Mvar is the foot of the perpendicular from Template:Mvar onto the image plane.) For both projections the apparent shape is a circle. In the parallel case the image of the origin Template:Mvar is the circle's center; in the central case main point Template:Mvar is the center.
Umbilical points
The focal hyperbola intersects the ellipsoid at its four umbilical points.[17]

Property of the focal ellipse

The focal ellipse together with its inner part can be considered as the limit surface (an infinitely thin ellipsoid) of the pencil of confocal ellipsoids determined by Template:Math for Template:Math. For the limit case one gets

rx=a,ry=b,l=3ac.

In higher dimensions and general position Template:Anchor

A hyperellipsoid, or ellipsoid of dimension n1 in a Euclidean space of dimension n, is a quadric hypersurface defined by a polynomial of degree two that has a homogeneous part of degree two which is a positive definite quadratic form.

One can also define a hyperellipsoid as the image of a sphere under an invertible affine transformation. The spectral theorem can again be used to obtain a standard equation of the form

x12a12+x22a22++xn2an2=1.

The volume of an Template:Mvar-dimensional hyperellipsoid can be obtained by replacing Template:Mvar by the product of the semi-axes Template:Math in the formula for the volume of a hypersphere:

V=πn2Γ(n2+1)a1a2an1πn(2eπn)n/2a1a2an

(where Template:Math is the gamma function).

As a quadric

If Template:Mvar is a real, symmetric, Template:Mvar-by-Template:Mvar positive-definite matrix, and Template:Mvar is a vector in โ„n, then the set of points Template:Math that satisfy the equation

(๐ฑ๐ฏ)๐–ณ๐‘จ(๐ฑ๐ฏ)=1

is an n-dimensional ellipsoid centered at Template:Mvar. The expression (๐ฑ๐ฏ)๐–ณ๐‘จ(๐ฑ๐ฏ) is also called the ellipsoidal norm of Template:Math. For every ellipsoid, there are unique Template:Mvar and Template:Math that satisfy the above equation.[18]Template:Rp

The eigenvectors of Template:Mvar are the principal axes of the ellipsoid, and the eigenvalues of Template:Mvar are the reciprocals of the squares of the semi-axes (in three dimensions these are Template:Math, Template:Math and Template:Math).[19] In particular:

  • The diameter of the ellipsoid is twice the longest semi-axis, which is twice the square-root of the reciprocal of the largest eigenvalue of Template:Mvar.
  • The width of the ellipsoid is twice the shortest semi-axis, which is twice the square-root of the reciprocal of the smallest eigenvalue of Template:Mvar.

An invertible linear transformation applied to a sphere produces an ellipsoid, which can be brought into the above standard form by a suitable rotation, a consequence of the polar decomposition (also, see spectral theorem). If the linear transformation is represented by a symmetric 3 ร— 3 matrix, then the eigenvectors of the matrix are orthogonal (due to the spectral theorem) and represent the directions of the axes of the ellipsoid; the lengths of the semi-axes are computed from the eigenvalues. The singular value decomposition and polar decomposition are matrix decompositions closely related to these geometric observations.

For every positive definite matrix

๐‘จ

, there exists a unique positive definite matrix denoted Template:Math, such that

๐‘จ=๐‘จ1/2๐‘จ1/2;

this notation is motivated by the fact that this matrix can be seen as the "positive square root" of

๐‘จ.

The ellipsoid defined by

(๐ฑ๐ฏ)๐–ณ๐‘จ(๐ฑ๐ฏ)=1

can also be presented as[18]Template:Rp

A1/2S(๐ŸŽ,1)+๐ฏ

where S(0,1) is the unit sphere around the origin.

Parametric representation

ellipsoid as an affine image of the unit sphere

The key to a parametric representation of an ellipsoid in general position is the alternative definition:

An ellipsoid is an affine image of the unit sphere.

An affine transformation can be represented by a translation with a vector Template:Math and a regular 3 ร— 3 matrix Template:Math:

๐ฑ๐Ÿ0+๐‘จ๐ฑ=๐Ÿ0+x๐Ÿ1+y๐Ÿ2+z๐Ÿ3

where Template:Math are the column vectors of matrix Template:Math.

A parametric representation of an ellipsoid in general position can be obtained by the parametric representation of a unit sphere (see above) and an affine transformation:

๐ฑ(θ,φ)=๐Ÿ0+๐Ÿ1cosθcosφ+๐Ÿ2cosθsinφ+๐Ÿ3sinθ,π2<θ<π2,0φ<2π.

If the vectors Template:Math form an orthogonal system, the six points with vectors Template:Math are the vertices of the ellipsoid and Template:Math are the semi-principal axes.

A surface normal vector at point Template:Math is

๐ง(θ,φ)=๐Ÿ2×๐Ÿ3cosθcosφ+๐Ÿ3×๐Ÿ1cosθsinφ+๐Ÿ1×๐Ÿ2sinθ.

For any ellipsoid there exists an implicit representation Template:Math. If for simplicity the center of the ellipsoid is the origin, Template:Math, the following equation describes the ellipsoid above:[20]

F(x,y,z)=det(๐ฑ,๐Ÿ2,๐Ÿ3)2+det(๐Ÿ1,๐ฑ,๐Ÿ3)2+det(๐Ÿ1,๐Ÿ2,๐ฑ)2det(๐Ÿ1,๐Ÿ2,๐Ÿ3)2=0

Applications

The ellipsoidal shape finds many practical applications:

Geodesy
Mechanics
Crystallography

Computer science

Lighting
Medicine

Dynamical properties

The mass of an ellipsoid of uniform density Template:Mvar is

m=Vρ=43πabcρ.

The moments of inertia of an ellipsoid of uniform density are

Ixx=15m(b2+c2),Iyy=15m(c2+a2),Izz=15m(a2+b2),Ixy=Iyz=Izx=0.

For Template:Math these moments of inertia reduce to those for a sphere of uniform density.

Artist's conception of Template:Dp, a Jacobi-ellipsoid dwarf planet, with its two moons

Ellipsoids and cuboids rotate stably along their major or minor axes, but not along their median axis. This can be seen experimentally by throwing an eraser with some spin. In addition, moment of inertia considerations mean that rotation along the major axis is more easily perturbed than rotation along the minor axis.[22]

One practical effect of this is that scalene astronomical bodies such as Template:Dp generally rotate along their minor axes (as does Earth, which is merely oblate); in addition, because of tidal locking, moons in synchronous orbit such as Mimas orbit with their major axis aligned radially to their planet.

A spinning body of homogeneous self-gravitating fluid will assume the form of either a Maclaurin spheroid (oblate spheroid) or Jacobi ellipsoid (scalene ellipsoid) when in hydrostatic equilibrium, and for moderate rates of rotation. At faster rotations, non-ellipsoidal piriform or oviform shapes can be expected, but these are not stable.

Fluid dynamics

The ellipsoid is the most general shape for which it has been possible to calculate the creeping flow of fluid around the solid shape. The calculations include the force required to translate through a fluid and to rotate within it. Applications include determining the size and shape of large molecules, the sinking rate of small particles, and the swimming abilities of microorganisms.[23]

In probability and statistics

The elliptical distributions, which generalize the multivariate normal distribution and are used in finance, can be defined in terms of their density functions. When they exist, the density functions Template:Mvar have the structure:

f(x)=kg((๐ฑμ)Σ1(๐ฑμ)๐–ณ)

where Template:Mvar is a scale factor, Template:Math is an Template:Mvar-dimensional random row vector with median vector Template:Math (which is also the mean vector if the latter exists), Template:Math is a positive definite matrix which is proportional to the covariance matrix if the latter exists, and Template:Mvar is a function mapping from the non-negative reals to the non-negative reals giving a finite area under the curve.[24] The multivariate normal distribution is the special case in which Template:Math for quadratic form Template:Mvar.

Thus the density function is a scalar-to-scalar transformation of a quadric expression. Moreover, the equation for any iso-density surface states that the quadric expression equals some constant specific to that value of the density, and the iso-density surface is an ellipsoid.

See also

Notes

  1. โ†‘ Template:Harvtxt
  2. โ†‘ F.W.J. Olver, D.W. Lozier, R.F. Boisvert, and C.W. Clark, editors, 2010, NIST Handbook of Mathematical Functions (Cambridge University Press), Section 19.33 Template:Cite web
  3. โ†‘ Template:Cite web
  4. โ†‘ Template:Cite web
  5. โ†‘ Template:Cite web
  6. โ†‘ Template:Cite web
  7. โ†‘ Final answers Template:Webarchive by Gerard P. Michon (2004-05-13). See Thomsen's formulas and Cantrell's comments.
  8. โ†‘ Template:Citation
  9. โ†‘ W. Bรถhm: Die FadenKonstruktion der Flรคchen zweiter Ordnung, Mathemat. Nachrichten 13, 1955, S. 151
  10. โ†‘ Staude, O.: Ueber Fadenconstructionen des Ellipsoides. Math. Ann. 20, 147โ€“184 (1882)
  11. โ†‘ Staude, O.: Ueber neue Focaleigenschaften der Flรคchen 2. Grades. Math. Ann. 27, 253โ€“271 (1886).
  12. โ†‘ Staude, O.: Die algebraischen Grundlagen der Focaleigenschaften der Flรคchen 2. Ordnung Math. Ann. 50, 398 - 428 (1898).
  13. โ†‘ D. Hilbert & S Cohn-Vossen: Geometry and the imagination, Chelsea New York, 1952, Template:ISBN, p. 20 .
  14. โ†‘ O. Hesse: Analytische Geometrie des Raumes, Teubner, Leipzig 1861, p. 287
  15. โ†‘ D. Hilbert & S Cohn-Vossen: Geometry and the Imagination, p. 24
  16. โ†‘ O. Hesse: Analytische Geometrie des Raumes, p. 301
  17. โ†‘ W. Blaschke: Analytische Geometrie, p. 125
  18. โ†‘ 18.0 18.1 Template:Cite Geometric Algorithms and Combinatorial Optimization
  19. โ†‘ Template:Cite web pp. 17โ€“18.
  20. โ†‘ Computerunterstรผtzte Darstellende und Konstruktive Geometrie. Template:Webarchive Uni Darmstadt (PDF; 3,4 MB), S. 88.
  21. โ†‘ Template:Cite journal
  22. โ†‘ Goldstein, H G (1980). Classical Mechanics, (2nd edition) Chapter 5.
  23. โ†‘ Dusenbery, David B. (2009).Living at Micro Scale, Harvard University Press, Cambridge, Massachusetts Template:Isbn.
  24. โ†‘ Frahm, G., Junker, M., & Szimayer, A. (2003). Elliptical copulas: applicability and limitations. Statistics & Probability Letters, 63(3), 275โ€“286.

References

Template:Commons category

ta:เฎจเฏ€เฎณเฏเฎฐเฏเฎฃเฏเฎŸเฏˆ