Principal axis theorem

Principal axis theorem In the mathematical fields of geometry and linear algebra, a principal axis is a certain line in a Euclidean space associated with an ellipsoid or hyperboloid, generalizing the major and minor axes of an ellipse or hyperbola. The principal axis theorem states that the principal axes are perpendicular, and gives a constructive procedure for finding them.

Mathematically, the principal axis theorem is a generalization of the method of completing the square from elementary algebra. In linear algebra and functional analysis, the principal axis theorem is a geometrical counterpart of the spectral theorem. It has applications to the statistics of principal components analysis and the singular value decomposition. In physics, the theorem is fundamental to the studies of angular momentum and birefringence.

Contenu 1 Motivation 2 Déclaration formelle 3 Voir également 4 References Motivation The equations in the Cartesian plane R2: {style d'affichage {commencer{aligné}{frac {x^{2}}{9}}+{frac {y ^{2}}{25}}&=1\[3pt]{frac {x^{2}}{9}}-{frac {y ^{2}}{25}}&=1end{aligné}}} définir, respectivement, an ellipse and a hyperbola. In each case, the x and y axes are the principal axes. This is easily seen, given that there are no cross-terms involving products xy in either expression. Cependant, the situation is more complicated for equations like {style d'affichage 5x^{2}+8xy+5y^{2}=1.} Here some method is required to determine whether this is an ellipse or a hyperbola. The basic observation is that if, by completing the square, the quadratic expression can be reduced to a sum of two squares then the equation defines an ellipse, whereas if it reduces to a difference of two squares then the equation represents a hyperbola: {style d'affichage {commencer{aligné}tu(X,y)^{2}+v(X,y)^{2}&=1qquad {texte{(ellipse)}}\tu(X,y)^{2}-v(X,y)^{2}&=1qquad {texte{(hyperbola)}}.fin{aligné}}} Ainsi, in our example expression, the problem is how to absorb the coefficient of the cross-term 8xy into the functions u and v. Officiellement, this problem is similar to the problem of matrix diagonalization, where one tries to find a suitable coordinate system in which the matrix of a linear transformation is diagonal. The first step is to find a matrix in which the technique of diagonalization can be applied.

The trick is to write the quadratic form as {style d'affichage 5x^{2}+8xy+5y^{2}={commencer{bmatrice}x¥d{bmatrice}}{commencer{bmatrice}5&4\4&5end{bmatrice}}{commencer{bmatrice}x\yend{bmatrice}}= mathbf {X} ^{textsf {J}}Amathbf {X} } where the cross-term has been split into two equal parts. The matrix A in the above decomposition is a symmetric matrix. En particulier, by the spectral theorem, it has real eigenvalues and is diagonalizable by an orthogonal matrix (orthogonally diagonalizable).

To orthogonally diagonalize A, one must first find its eigenvalues, and then find an orthonormal eigenbasis. Calculation reveals that the eigenvalues of A are {style d'affichage lambda _{1}=1,quad lambda _{2}=9} with corresponding eigenvectors {style d'affichage mathbf {v} _{1}={commencer{bmatrice}1\-1fin{bmatrice}},quad mathbf {v} _{2}={commencer{bmatrice}1\1fin{bmatrice}}.} Dividing these by their respective lengths yields an orthonormal eigenbasis: {style d'affichage mathbf {tu} _{1}={commencer{bmatrice}1/{sqrt {2}}\-1/{sqrt {2}}fin{bmatrice}},quad mathbf {tu} _{2}={commencer{bmatrice}1/{sqrt {2}}\1/{sqrt {2}}fin{bmatrice}}.} Now the matrix S = [u1 u2] is an orthogonal matrix, since it has orthonormal columns, and A is diagonalized by: {displaystyle A=SDS^{-1}=SDS^{textsf {J}}={commencer{bmatrice}1/{sqrt {2}}&1/{sqrt {2}}\-1/{sqrt {2}}&1/{sqrt {2}}fin{bmatrice}}{commencer{bmatrice}1&0\0&9end{bmatrice}}{commencer{bmatrice}1/{sqrt {2}}&-1/{sqrt {2}}\1/{sqrt {2}}&1/{sqrt {2}}fin{bmatrice}}.} This applies to the present problem of "diagonalizing" the quadratic form through the observation that {style d'affichage 5x^{2}+8xy+5y^{2}= mathbf {X} ^{textsf {J}}Amathbf {X} = mathbf {X} ^{textsf {J}}la gauche(SDS^{textsf {J}}droit)mathbf {X} =gauche(S^{textsf {J}}mathbf {X} droit)^{textsf {J}}Dleft(S^{textsf {J}}mathbf {X} droit)=1left({frac {x-y}{sqrt {2}}}droit)^{2}+9la gauche({frac {x+y}{sqrt {2}}}droit)^{2}.} Ainsi, the equation {style d'affichage 5x^{2}+8xy+5y^{2}=1} is that of an ellipse, since the left side can be written as the sum of two squares.

It is tempting to simplify this expression by pulling out factors of 2. Cependant, it is important not to do this. The quantities {displaystyle c_{1}={frac {x-y}{sqrt {2}}},quad c_{2}={frac {x+y}{sqrt {2}}}} have a geometrical meaning. They determine an orthonormal coordinate system on R2. Autrement dit, they are obtained from the original coordinates by the application of a rotation (and possibly a reflection). Par conséquent, one may use the c1 and c2 coordinates to make statements about length and angles (particularly length), which would otherwise be more difficult in a different choice of coordinates (by rescaling them, par exemple). Par exemple, the maximum distance from the origin on the ellipse c12 + 9c22 = 1 occurs when c2 = 0, so at the points c1 = ±1. De la même manière, the minimum distance is where c2 = ±1/3.

It is possible now to read off the major and minor axes of this ellipse. These are precisely the individual eigenspaces of the matrix A, since these are where c2 = 0 or c1 = 0. Symboliquement, the principal axes are {style d'affichage E_{1}={texte{span}}la gauche({commencer{bmatrice}1/{sqrt {2}}\-1/{sqrt {2}}fin{bmatrice}}droit),quad E_{2}={texte{span}}la gauche({commencer{bmatrice}1/{sqrt {2}}\1/{sqrt {2}}fin{bmatrice}}droit).} To summarize: The equation is for an ellipse, since both eigenvalues are positive. (Autrement, if one were positive and the other negative, it would be a hyperbola.) The principal axes are the lines spanned by the eigenvectors. The minimum and maximum distances to the origin can be read off the equation in diagonal form.

Using this information, it is possible to attain a clear geometrical picture of the ellipse: to graph it, par exemple.

Formal statement The principal axis theorem concerns quadratic forms in Rn, which are homogeneous polynomials of degree 2. Any quadratic form may be represented as {style d'affichage Q(mathbf {X} )= mathbf {X} ^{textsf {J}}Amathbf {X} } where A is a symmetric matrix.

The first part of the theorem is contained in the following statements guaranteed by the spectral theorem: The eigenvalues of A are real. A is diagonalizable, and the eigenspaces of A are mutually orthogonal.

En particulier, A is orthogonally diagonalizable, since one may take a basis of each eigenspace and apply the Gram-Schmidt process separately within the eigenspace to obtain an orthonormal eigenbasis.

For the second part, suppose that the eigenvalues of A are λ1, ..., λn (possibly repeated according to their algebraic multiplicities) and the corresponding orthonormal eigenbasis is u1, ..., un. Alors {style d'affichage Q(mathbf {X} )=lambda _{1}c_{1}^{2}+lambda _{2}c_{2}^{2}+dots +lambda _{n}c_{n}^{2},} where the ci are the coordinates with respect to the given eigenbasis. Par ailleurs, The i-th principal axis is the line determined by the n − 1 equations cj = 0, j ≠ i. This axis is the span of the vector ui. See also Sylvester's law of inertia References Strang, Gilbert (1994). Introduction to Linear Algebra. Wellesley-Cambridge Press. ISBN 0-9614088-5-5. Catégories: Theorems in geometryTheorems in linear algebra

Si vous voulez connaître d'autres articles similaires à Principal axis theorem vous pouvez visiter la catégorie Théorèmes de géométrie.

Laisser un commentaire

Votre adresse email ne sera pas publiée.

Monter

Nous utilisons nos propres cookies et ceux de tiers pour améliorer l'expérience utilisateur Plus d'informations