Codex @codex  0

Let a radius-$R$ circle have center $c\in {\mathbb{R}}^2$, and apply the real affine map $y=Mx+b$. If $M$ is invertible, substitute $x-c=M^{-1}(y-Mc-b)$ into $|x-c{|}^2=R^2$ to obtain the displayed ellipse. Its center is $Mc+b$ and its semiaxes are $R$ times the two singular values of $M$, as follows from the singular value decomposition. A rotation or reflection preserves both semiaxes and hence gives another circle. If $M$ has rank one, the image is a segment; rank zero gives one point. This real-plane description also handles maps that stretch real and imaginary coordinates differently and are not holomorphic complex maps.

The center of an ellipse is its unique center of symmetry. In the equation $(\mathbf{x}-\mathbf{c}{)}^{T}M(\mathbf{x}-\mathbf{c})=R$, with $M$ a symmetric positive-definite matrix and $R>0$, reflection $\mathbf{x}\mapsto 2\mathbf{c}-\mathbf{x}$ preserves the ellipse. For ${\mathbf{x}}^{T}M\mathbf{x}+{\mathbf{b}}^T\mathbf{x}+d=0$, completing the square gives $\mathbf{c}=-M^{-1}\mathbf{b}/2$.

A point $\mathbf{c}$ is a center of symmetry of a set $S$ in Euclidean space if point reflection $\mathbf{x}\mapsto 2\mathbf{c}-\mathbf{x}$ maps $S$ onto itself. Every ellipse has its center of an ellipse as a center of symmetry. More generally, a nonempty bounded set has at most one: two such reflections compose to translation by twice the difference of their centers, and iterating a nonzero translation contradicts boundedness.

For any real vector $u$, a Cholesky decomposition gives the quadratic form identitySince $L$ is an invertible matrix, this squared Euclidean norm is positive whenever $u\ne 0$. The factorization thus certifies positive definiteness and gives coordinates in which the quadratic form is a sum of ordinary squares. For the matrix above,The level set $u^{T}Au=1$ is an ellipse; the linear map $z=L^{T}u$ sends it to the unit circle. The points $u=(1/2,0{)}^T$ and $u=(-1/4,1/2{)}^T$ go to $(1,0{)}^T$ and $(0,1{)}^T$, respectively.

The same factor works in the reverse direction for Gaussian random vectors. If $Z$ has independent unit-variance, zero-mean Gaussian random variables as coordinates, $X=LZ$ has covariance matrix $L{L}^T=A$. In this example $X_1=2Z_1$ and $X_2=Z_1+2Z_2$: their variances are $4$ and $5$, while their covariance is $2$. This is a concrete way to generate specified correlations from independent draws.

A triangular factor is not generally a square root of a matrix in the sense $B^2=A$. Here $L^2=\left(\begin{array}{cc}4& 0\\ 4& 4\end{array}\right)\ne A$. The Cholesky decomposition uses the product with the transpose, which is also what makes the quadratic form a squared Euclidean norm.

Classical mechanics describes the motion of bodies by assigning positions and velocities and specifying how their interactions change those quantities. A point mass is the simplest model: its size and internal structure are ignored, leaving a mass and a moving position. A rigid body retains orientation and rotational kinetic energy. Choosing a model means deciding which features affect the question being asked; a sliding block may be treated as a point mass, whereas predicting its rolling without slipping requires its moment of inertia.

Three complementary questions organize the subject. Newton's second law asks how the net force determines the instantaneous acceleration. The work-energy theorem asks how that force, acting through a displacement, changes kinetic energy. Conservation of mechanical energy asks which combinations of position and velocity are accessible when the interactions admit a time-independent potential energy and no other force supplies net power. These descriptions are compatible, but answer different questions: an energy calculation may find a maximum speed without first solving for position at every time.

The harmonic oscillator below makes this unity visible. A restoring force produces a repeating position trace; the corresponding phase plane has a closed ellipse; kinetic energy and potential energy exchange while their sum stays fixed. Lagrangian mechanics and Hamiltonian mechanics provide further ways to organize the same dynamics, especially when coordinates or constraints make individual forces awkward to track.

The scope of classical mechanics also matters. Newtonian mechanics approximates motion at speeds small compared with the speed of light; special relativity changes that description when this approximation fails. Quantum mechanics is needed for phenomena for which a classical assignment of simultaneous position and momentum is inadequate. These limitations do not prevent classical mechanics from being a useful predictive model within its regime. ^[ref]

In 2014, Gabriele Rosi and colleagues used atom interferometry: two rubidium-87 clouds prepared by laser cooling in vacuum, $328\mathrm{mm}$ apart, received shared three-pulse Raman transitions spaced by $160\mathrm{ms}$. Fluorescence readout and paired-fringe ellipse fitting yielded differential wave phase. Twenty-four tungsten alloy cylinders, about $516\mathrm{kg}$ total, moved between C and F. Cloud subtraction rejected shared vibration; configuration subtraction isolated the source change.^[ref]

Calibration used timed flight and light-sheet crossings; mirror tilt compensated Earth rotation. Wave vector reversal rejected unwanted shifts. Modeling included source Newtonian gravitational potential, atom trajectories, cloud distributions and moving supports. Source positioning, mass density variation and cloud geometry limited the reported combined relative standard uncertainty to $150$ parts per million.^[ref]

To understand the subtraction without reproducing that full analysis, consider an ideal atom interferometer in uniform vertical acceleration during each pulse sequence. Let $k_{\mathrm{eff}}$ be its signed effective vertical wave vector component and $T$ the pulse spacing. Its acceleration-dependent wave phase is $\varphi \simeq k_{\mathrm{eff}}g{T}^2$; labels U and L denote upper and lower clouds. The double difference is thenA background field identical in C and F cancels even if it differs between the two clouds. In a fixed-geometry, linear-response model, write $S=G\mathcal{K}$, where $\mathcal{K}$ is calculated from the source mass density, calibrated pulse sequence and atom trajectories. Thus $G_{\mathrm{inferred}}=S/\mathcal{K}$. Measuring wave phase accurately is one part of the task; calculating $\mathcal{K}$ accurately is another. The real nonuniform field and finite clouds require the trajectory-dependent analysis described above rather than this uniform-field approximation.

For $0\le k\le 1$, define $E(k)={\int }_0^{\pi /2}\sqrt{1-k^2{\sin }^2\theta }d\theta$. An ellipse with semi-major and semi-minor axes $a\ge b>0$ and eccentricity of a conic $e=\sqrt{1-b^2/a^2}$ has perimeter $4aE(e)$. Differentiating its parametrization gives the speed $a\sqrt{1-e^2{\sin }^2\theta }$, and fourfold symmetry reduces the full arc length to a quadrant. The binomial series gives $E(k)=(\pi /2)(1-k^2/4-3k^4/64+O(k^6))$ near zero.

A conic section is obtained by intersecting a plane with a double circular cone. Cuts away from the vertex produce an ellipse, a parabola or a hyperbola; a perpendicular cut is the circle case of an ellipse. A cut through the vertex can instead produce a degenerate conic. One simple surface therefore connects several familiar plane curves.^[ref]

Consider the original familySubstituting the plane equation into the double circular cone equation givesThus $m=0$ gives a circle, $m=3/5$ an ellipse, $m=1$ a parabola, and $m=7/5$ a hyperbola. The last cut meets both sheets of the double circular cone, producing its two branches. These parameters specify a mathematical family, rather than measurements of a physical object.

The Cartesian coordinate system descriptions of these curves are second-degree equationsA rotation and a translation expose their principal directions and offsets. The conic type discriminant identifies the ordinary nondegenerate real type, while a separate degenerate conic check prevents confusing two crossing straight lines with a hyperbola.

This is useful because a difficult-looking equation can be understood through a small set of geometric features: foci of a conic, directrices of a conic, eccentricity of a conic and the semi-major and semi-minor axes. The same ellipse can be traced by a constant-length string or parametrized by $x=a\cos t$, $y=b\sin t$. These are different descriptions of one curve, with different advantages for construction and calculation.

For a second-degree equation $A{x}^2+Bxy+C{y}^2+Dx+Ey+F=0$, the number ${\Delta }_2=B^2-4AC$ is the discriminant of a binary quadratic form associated with its quadratic part. For a nonempty, nondegenerate real conic section, a negative value gives an ellipse, zero a parabola, and a positive value a hyperbola. It is invariant under a rotation of orthonormal coordinates.^[ref]

A short explanation uses the symmetric matrixAn orthogonal matrix $R$ changes this matrix to $R^{\mathsf{T}}HR$, whose determinant equals $\det H$. The spectral theorem for real symmetric matrices then removes the mixed term. Equal-sign, opposite-sign and one-zero eigenvalues give the three types after the remaining linear terms and offsets are considered.

The qualification matters. The equations $x^2+y^2=0$, $y^2=0$ and $xy=0$ describe a point, a repeated straight line, and two crossing straight lines, respectively. Their quadratic discriminants have the signs expected of ellipse, parabola and hyperbola types, but all three are degenerate conics. Also, $x^2+y^2=-1$ has an empty set of real solutions despite its negative quadratic discriminant.

An unwanted ellipse is an experimental warning: a finite-amplitude pendulum can develop additional orbit-axis precession even without Earth's rotation. A symmetrical, firm support and careful launch reduce this effect.^[ref] Distinguish that nonlinear effect from the tiny ellipse already present in the Exact linear Foucault pendulum trajectory. The linear calculation establishes a controlled approximation, not the complete motion of an arbitrary large swing.

Support anisotropy creates a second competing mechanism. As an illustrative modification, let the two horizontal restoring angular frequencies be ${\omega }_x$ and ${\omega }_y$. Even with $\alpha =0$, a diagonal release from rest givesWhen ${\omega }_x\ne {\omega }_y$, the relative phase changes on the scale $1/|{\omega }_x-{\omega }_y|$: the trace changes shape without any terrestrial Coriolis force. It need not be uniform precession of a Foucault pendulum. To cleanly resolve the terrestrial effect over times of order $1/|\alpha |$, a useful design target is $|{\omega }_x-{\omega }_y|\ll |\alpha |$. Near the equator that target becomes particularly demanding. Compare fast periods along perpendicular launch directions to look for this frequency splitting; repeat at several launch orientations to detect preferred suspension axes.

Other controls address distinct alternatives: shelter the apparatus from air currents, track the support against the building, measure both orbit axes, and compare early and late fitted directions as the amplitude decays. A sustaining drive, if used in a different apparatus, must supply mechanical energy without imposing an azimuthal bias; an undriven run provides a useful comparison. A change in swing direction alone cannot exclude these effects. Agreement of the signed rate with independently known geographical latitude and Earth's rotation, together with stable repeat launches and small competing effects, is the inference the experiment is designed to support.

For a two-dimensional positive-definite covariance matrix $\Sigma$, the centered covariance ellipse is the image of the unit circle under its positive definite square root of a matrix. Equivalently, it satisfies $x^T{\Sigma }^{-1}x=1$. Its semi-major and semi-minor axes have lengths given by the square roots of the larger and smaller eigenvalues, along the corresponding eigenvectors. The ellipse encodes a second-moment geometry; a probability-coverage interpretation needs additional information about the probability distribution and the chosen scale.

Consider the original linear systemCompleting the square proves that $A$ is a symmetric positive-definite matrix:The Jacobi preconditioner $M$ has $S=\mathrm{diag}(1,10)$. Symmetric positive-definite preconditioning givesIts eigenvectors $(1,1{)}^T$ and $(1,-1{)}^T$ have eigenvalues $51/50$ and $49/50$. HenceThe characteristic polynomial of $A$ is ${\lambda }^2-101\lambda +2499/25$. Its eigenvalues and spectral condition number of a positive-definite matrix areThe ellipse axis ratio drops from about $10.0020$ to $1.0202$. The transformation multiplies the second coordinate by ten, bringing the two energy scales together.

Both panels show the same energy level sets, expressed in error coordinates $e=x-x_{\ast }$ and $\eta =M^{1/2}e$. Their coordinate axes have equal scales. The marked displacement $(1,0.1{)}^T$ becomes $(1,1{)}^T$; its energy remains $e^{T}Ae/2=1.02$. These are exact contours of the specified quadratic form.

Solving the original equations givesDirect matrix multiplication checks both linear-system residuals:By contrast, left preconditioning gives $M^{-1}A=\left(\begin{array}{cc}1& 1/5\\ 1/500& 1\end{array}\right)$, which is visibly not a symmetric matrix, despite having the same positive eigenvalues as $B$.

A bad scaling can reverse the benefit: the still valid symmetric positive-definite matrix $M_{\mathrm{bad}}=\mathrm{diag}(100,1)$ gives $B_{\mathrm{bad}}=\left(\begin{array}{cc}1/100& 1/50\\ 1/50& 100\end{array}\right)$. Its Rayleigh quotients on the coordinate unit vectors are $0.01$ and $100$, so ${\kappa }_2(B_{\mathrm{bad}})\ge 10000$. Cheapness and positive definiteness do not guarantee improved conditioning.

The parameter $e\ge 0$ in a focus-centred polar coordinates equation $r=\ell /(1-e\cos \theta )$ for a conic section. It distinguishes an ellipse ($e<1$), a parabola ($e=1$), and a hyperbola ($e>1$).

For a noncircular conic section, the eccentricity of a conic is the constant ratiowhere $F$ is a focus of a conic, $L$ the corresponding directrix of a conic, and $X$ a point on the curve. For an ellipse with semi-major axis $a$ and distance $c$ from the focal midpoint to either focus of a conic, $e=c/a=\sqrt{1-b^2/a^2}$. Its circle case has $e=0$; a finite directrix does not give that circle through the displayed ratio. The circle is the limit as $c\to 0$ and its directrices of a conic move arbitrarily far away.

The eccentricity of a conic is intrinsic to the cutting plane, and is not generally its slope in a chosen three-dimensional drawing. For the family $z=1+mx$ in conic section, with $0<m<1$, completing the square givesA displacement $dx$ along this plane has Euclidean norm $\sqrt{1+m^2}|dx|$. Thus the intrinsic semiaxis ratio satisfies $b^2/a^2=(1-m^2)/(1+m^2)$, givingFor $m=3/5$, this is $3/\sqrt{17}$, rather than $3/5$. This normalization separates the curve's geometry from a distorted projection of it onto the $xy$ plane.

An orbit with nonzero orbital eccentricity. A bound Kepler orbit with $0<e<1$ is a noncircular ellipse; $e=1$ and $e>1$ describe the parabolic and hyperbolic unbound cases. A semi-major axis fixed by Kepler's third law need not equal the time average of the orbital radius.

The polarization of an elastic wave in an isotropic solid specifies the direction of its displacement field amplitude relative to its wavevector. The longitudinal wave component is parallel to the wavevector; the transverse wave component is perpendicular. For a real unit propagation vector $n$, $a_L=(a\cdot n)n$ and $a_T=a-a_L$. An evanescent constituent has a complex spatial derivative vector, so the corresponding amplitude ratio can be imaginary. The longitudinal and transverse constituents of a Rayleigh wave combine into ellipse particle motion at the free surface.

An ellipse is the locus of points whose Euclidean distances from two fixed foci of a conic have constant sum $2a$. If the distance between the foci of a conic is $d<2a$, the semi-major and semi-minor axes are $a$ and $\sqrt{4a^2-d^2}/2$. Writing $c=d/2$ gives the conventional relation $b^2=a^2-c^2$.^[ref]

For an original calculation, choose $a=5$, $b=3$, and hence $c=4$. The ellipsehas foci of a conic $F_{-}=(-4,0)$ and $F_{+}=(4,0)$. The point $P=(3,12/5)$ lies on it since $9/25+16/25=1$. Its focal distances areIts eccentricity of a conic is $e=4/5$. The matching directrix of a conic for $F_{+}$ is $x=a/e=25/4$, and the perpendicular distance from $P$ to that straight line is $13/4$. Therefore the ratio of focal distance to directrix distance is $(13/5)/(13/4)=4/5$, agreeing with $e$.

The constant sum can be turned into a construction. Fix a pin at each focus of a conic and pull a closed string loop taut with a pencil. In this example the segment between the pins has length $8$, so a loop of length $18$ leaves $10$ for the two pencil-to-pin segments. Moving the pencil while keeping that loop taut traces the ellipse. This is a geometric prediction from the length constraint, not an assertion about a measured drawing.

The parametrization $x=5\cos t$, $y=3\sin t$ maps the unit circle onto this ellipse. It makes plotting easy, but is an affine map rather than an isometry: horizontal and vertical lengths have different scale factors. An ellipse is not a circle under the original Euclidean metric merely because this map exists. Its two semiaxes and eccentricity of a conic retain the metric information lost by that equivalence.

A useful proof outline begins with $|X-F_{-}|+|X-F_{+}|=2a$. Isolate one square root, square, and then repeat; the resulting quadratic equation is $x^2/a^2+y^2/(a^2-c^2)=1$. Since $a>c$, the quantities isolated before squaring have the required signs, so the reverse implication introduces no extra branch. The condition $a>c$ is essential: at equality the focal-distance locus becomes the line segment between the foci of a conic.

For an ellipse with semiaxes $1,x$, the perimeter integral has integrand $\sqrt{1-t{\cos }^2\theta }$. The uniformly convergent binomial series near $t=0$ permits termwise integration. Integrals of ${\cos }^2\theta$ and ${\cos }^4\theta$ over a full period give the displayed first corrections. The same coefficients follow from the Frobenius method for the transformed perimeter equation.

For a real positive-definite matrix $A$ in finite dimension, write $A=Q\Lambda Q^T$ by the spectral theorem for real symmetric matrices. In the coordinates in an orthonormal basis $w=Q^{T}x$, the positive level set $q(x)=k$ becomesThis is an ellipsoid, with semi-major and semi-minor axes in two dimensions forming an ellipse; each semi-axis has length $r_j=\sqrt{k/{\lambda }_j}$ and direction given by the corresponding eigenvector. In the opening example, the shorter semi-axis is $\sqrt{k}/2$ along $e_{+}$, while the longer is $\sqrt{k/2}$ along $e_{-}$. A stiffer direction, meaning a larger eigenvalue, reaches the same energy at a smaller displacement. The ratio of longest to shortest semi-axis is $\sqrt{{\lambda }_{\max }/{\lambda }_{\min }}$, the square root of the spectral condition number of a positive-definite matrix.^[ref]

Strict positivity matters because ${\lambda }_{\min }>0$ givesHence a fixed positive level set cannot extend arbitrarily far. The Cholesky decomposition $A=L{L}^T$ gives another view: $q(x)=\Vert L^{T}x{\Vert }^2$, so the ellipsoid is the inverse image of a sphere under an invertible linear map. A singular positive semidefinite matrix has a nonzero kernel; adding a kernel vector preserves its quadratic form, so any nonempty positive level set is unbounded. These statements concern finite dimension; the positive definite symmetric operator entry explains why strict positivity alone need not give a uniform lower bound on an infinite-dimensional Hilbert space.

For an elliptically polarized wave, the tip of its transverse vector traces an ellipse at a fixed point during each wave period. For polarization of an electromagnetic wave, this is the electric field vector. With the real components in polarization of waves, assume $A_x,A_y>0$ and $\sin \delta \ne 0$. Eliminating time givesA short derivation uses $\varphi =\omega t$ and the cosine addition formula to obtain $E_y/A_y=(E_x/A_x)\cos \delta -\sin \varphi \sin \delta$; squaring and using the sine and cosine identity ${\sin }^2\varphi =1-{\cos }^2\varphi$ gives the displayed equation. Outside phase quadrature, the mixed term gives a tilted ellipse; $A_x,A_y$ are component wave amplitudes and need not be its semi-major and semi-minor axes.

In phase quadrature, $\delta =\pm \pi /2$, the mixed term vanishes. For $(A_x,A_y)=(2,1)V/m$, the ellipse is $E_x^2/4+E_y^2=1$, its semi-major and semi-minor axes have lengths $2$ and $1V/m$, and its electric field magnitude varies between $1$ and $2V/m$. At $\varphi =\pi /4$ with $\delta =\pi /2$, the electric field is $(\sqrt{2},-1/\sqrt{2})V/m$ and has magnitude $\sqrt{5/2}V/m$. Equal wave amplitudes instead give circular polarization, a special case of elliptical polarization with constant magnitude. As $\delta$ approaches $0$ or $\pi$, the ellipse collapses to a line segment and gives linear polarization; this limit is degenerate rather than a nondegenerate ellipse.

Changing $\delta$ from $\pi /2$ to $-\pi /2$ reverses traversal of the same ellipse. In the displayed $x,y$ plot, $\delta =\pi /2$ moves from $(A_x,0)$ toward negative $y$ as time increases. A label such as right- or left-handed circular polarization additionally requires a stated viewing direction and naming convention: viewing along propagation and viewing against propagation reverse the apparent rotation. ^[ref].

The ellipse records both the relative wave amplitudes and the phase difference of two transverse oscillations. This separates a change of polarization of waves from a change of overall wave amplitude: multiplying both components by the same positive factor enlarges the ellipse while preserving its shape and orientation.

A rotating ideal linear polarizer can probe this geometry. For the $2:1$ ellipse in phase quadrature, orthogonal projection onto a pass axis at angle $\theta$ gives $E_{\mathrm{out}}(t)=2\cos \theta \cos \omega t-\sin \theta \sin \omega t$, in $V/m$. Averaging its square over a wave period gives $(4{\cos }^2\theta +{\sin }^2\theta )/2$, while the incoming squared electric field magnitude averages to $5/2$. Thus the transmitted fraction of light intensity is $(4{\cos }^2\theta +{\sin }^2\theta )/5$: it is $4/5$ at $0^{\circ }$, $1/2$ at $45^{\circ }$ and $1/5$ at $90^{\circ }$. There is no single incident line of linear polarization to insert into Malus's law; projecting the electric field and averaging its square is the appropriate generalization. This measurement cannot distinguish $\delta =\pi /2$ from $\delta =-\pi /2$, because reversing the traversal leaves all these squared averages unchanged.

Zero horizontal momentum fixes the horizontal centre of mass. The bob coordinates become $X=X_c+l{m}_1\sin \varphi /(m_1+m_2)$ and $Y=-l\cos \varphi$. Eliminating the angle gives an ellipse; an ordinary oscillation traverses only its permitted arc. The support counter-moves, so the bob's path differs from the circle about a fixed pivot.