# Eigenvalues and eigenvectors

# Eigenvalues and eigenvectors

In linear algebra, an **eigenvector** (/ˈaɪɡənˌvɛktər/) or **characteristic vector** of a linear transformation is a nonzero vector that changes at most by a scalar factor when that linear transformation is applied to it.

There is a direct correspondence between *n*-by-*n*n*-dimensional vector space basis of the vector space. For this reason, in a finite-dimensional vector space, it is equivalent to define eigenvalues and eigenvectors using either the language of matrices or the language of linear transformations.^{[2]}^{[3]}

Geometrically, an eigenvector, corresponding to a real nonzero eigenvalue, points in a direction in which it is stretched by the transformation and the eigenvalue is the factor by which it is stretched. If the eigenvalue is negative, the direction is reversed.^{[4]} Loosely speaking, in a multidimensional vector space, the eigenvector is not rotated. However, in a one-dimensional vector space, the concept of rotation is meaningless.

Formal definition

If T is a linear transformation from a vector space V over a field F into itself and **v** is a vector in V that is not the zero vector, then **v** is an eigenvector of T if *T* (**v**) is a scalar multiple of **v**. This condition can be written as the equation

where λ is a scalar in the field F, known as the **eigenvalue**, **characteristic value**, or **characteristic root** associated with the eigenvector **v**.

If the vector space V is finite-dimensional, then the linear transformation T can be represented as a square matrix *A*, and the vector **v** by a column vector, rendering the above mapping as a matrix multiplication on the left-hand side and a scaling of the column vector on the right-hand side in the equation

Overview

Eigenvalues and eigenvectors feature prominently in the analysis of linear transformations.

The prefix eigen- is adopted from the German word eigen for "proper", "characteristic".^{[5]} Originally utilized to study principal axes of the rotational motion of rigid bodies, eigenvalues and eigenvectors have a wide range of applications, for example in stability analysis, vibration analysis, atomic orbitals, facial recognition, and matrix diagonalization.

In essence, an eigenvector **v** of a linear transformation *T* is a nonzero vector that, when *T* is applied to it, does not change direction. Applying *T* to the eigenvector only scales the eigenvector by the scalar value *λ*, called an eigenvalue. This condition can be written as the equation

referred to as the **eigenvalue equation** or **eigenequation**. In general, *λ* may be any scalar. For example, *λ* may be negative, in which case the eigenvector reverses direction as part of the scaling, or it may be zero or complex.

The Mona Lisa example pictured here provides a simple illustration. Each point on the painting can be represented as a vector pointing from the center of the painting to that point. The linear transformation in this example is called a shear mapping. Points in the top half are moved to the right and points in the bottom half are moved to the left proportional to how far they are from the horizontal axis that goes through the middle of the painting. The vectors pointing to each point in the original image are therefore tilted right or left and made longer or shorter by the transformation. Points *along* the horizontal axis do not move at all when this transformation is applied. Therefore, any vector that points directly to the right or left with no vertical component is an eigenvector of this transformation because the mapping does not change its direction. Moreover, these eigenvectors all have an eigenvalue equal to one because the mapping does not change their length, either.

Alternatively, the linear transformation could take the form of an *n* by *n* matrix, in which case the eigenvectors are *n* by 1 matrices that are also referred to as eigenvectors. If the linear transformation is expressed in the form of an *n* by *n* matrix *A*, then the eigenvalue equation above for a linear transformation can be rewritten as the matrix multiplication

where the eigenvector *v* is an *n* by 1 matrix. For a matrix, eigenvalues and eigenvectors can be used to decompose the matrix, for example by diagonalizing it.

Eigenvalues and eigenvectors give rise to many closely related mathematical concepts, and the prefix *eigen-* is applied liberally when naming them:

The set of all eigenvectors of a linear transformation, each paired with its corresponding eigenvalue, is called the

**eigensystem**of that transformation.^{[6]}^{[7]}The set of all eigenvectors of

*T*corresponding to the same eigenvalue, together with the zero vector, is called an**eigenspace**or**characteristic space**of*T*.^{[8]}If a set of eigenvectors of

*T*forms a basis of the domain of*T*, then this basis is called an**eigenbasis**.

History

Eigenvalues are often introduced in the context of linear algebra or matrix theory. Historically, however, they arose in the study of quadratic forms and differential equations.

In the 18th century Euler studied the rotational motion of a rigid body and discovered the importance of the principal axes.^{[9]} Lagrange realized that the principal axes are the eigenvectors of the inertia matrix.^{[10]} In the early 19th century, Cauchy saw how their work could be used to classify the quadric surfaces, and generalized it to arbitrary dimensions.^{[11]} Cauchy also coined the term *racine caractéristique* (characteristic root) for what is now called *eigenvalue*; his term survives in characteristic equation.^{[12]}

Fourier used the work of Laplace and Lagrange to solve the heat equation by separation of variables in his famous 1822 book Théorie analytique de la chaleur.^{[13]} Sturm developed Fourier's ideas further and brought them to the attention of Cauchy, who combined them with his own ideas and arrived at the fact that real symmetric matrices have real eigenvalues.^{[14]} This was extended by Hermite in 1855 to what are now called Hermitian matrices.^{[15]} Around the same time, Brioschi proved that the eigenvalues of orthogonal matrices lie on the unit circle,^{[14]} and Clebsch found the corresponding result for skew-symmetric matrices.^{[15]} Finally, Weierstrass clarified an important aspect in the stability theory started by Laplace by realizing that defective matrices can cause instability.^{[14]}

In the meantime, Liouville studied eigenvalue problems similar to those of Sturm; the discipline that grew out of their work is now called Sturm–Liouville theory.^{[16]} Schwarz studied the first eigenvalue of Laplace's equation on general domains towards the end of the 19th century, while Poincaré studied Poisson's equation a few years later.^{[17]}

At the start of the 20th century, Hilbert studied the eigenvalues of integral operators by viewing the operators as infinite matrices.^{[18]} He was the first to use the German word *eigen*, which means "own", to denote eigenvalues and eigenvectors in 1904,^{[19]} though he may have been following a related usage by Helmholtz. For some time, the standard term in English was "proper value", but the more distinctive term "eigenvalue" is standard today.^{[20]}

The first numerical algorithm for computing eigenvalues and eigenvectors appeared in 1929, when Von Mises published the power method. One of the most popular methods today, the QR algorithm, was proposed independently by John G.F. Francis^{[21]} and Vera Kublanovskaya^{[22]} in 1961.^{[23]}^{[24]}

Eigenvalues and eigenvectors of matrices

Eigenvalues and eigenvectors are often introduced to students in the context of linear algebra courses focused on matrices.^{[25]}^{[26]} Furthermore, linear transformations over a finite-dimensional vector space can be represented using matrices,^{[27]}^{[3]} which is especially common in numerical and computational applications.^{[28]}

Consider n-dimensional vectors that are formed as a list of n scalars, such as the three-dimensional vectors

These vectors are said to be scalar multiples of each other, or parallel or collinear, if there is a scalar λ such that

In this case λ = −1/20.

Now consider the linear transformation of n-dimensional vectors defined by an n by n matrix A,

or

where, for each row,

- .

If it occurs that v and w are scalar multiples, that is if

then v is an **eigenvector** of the linear transformation A and the scale factor λ is the **eigenvalue** corresponding to that eigenvector. Equation (1) is the **eigenvalue equation** for the matrix A.

Equation (1

where I is the n by n identity matrix and 0 is the zero vector.

Eigenvalues and the characteristic polynomial

Equation (2) has a nonzero solution *v* if and only if the determinant of the matrix (*A* − *λI*) is zero. Therefore, the eigenvalues of *A* are values of *λ* that satisfy the equation

Using Leibniz' rule for the determinant, the left-hand side of Equation (3) is a polynomial function of the variable *λ* and the degree of this polynomial is *n*, the order of the matrix *A*. Its coefficients depend on the entries of *A*, except that its term of degree *n* is always (−1)*n λ n*. This polynomial is called the* characteristic polynomial* of

*A*. Equation (3) is called the

*characteristic equation*or the

*secular equation*of

*A*.

The fundamental theorem of algebra implies that the characteristic polynomial of an *n*-by-*n*matrix* A*, being a polynomial of degree* nfactored into the product of

where each *λ**i*may be real but in general is a complex number. The numbers λ*, *λ*2,... *λ**n*, which may not all have distinct values, are roots of the polynomial and are the eigenvalues of* A*.

As a brief example, which is described in more detail in the examples section later, consider the matrix

Taking the determinant of (*M* − *λI*), the characteristic polynomial of *M* is

Setting the characteristic polynomial equal to zero, it has roots at λ = 1 and λ = 3, which are the two eigenvalues of *M*. The eigenvectors corresponding to each eigenvalue can be found by solving for the components of *v* in the equation *Mv* = *λv*. In this example, the eigenvectors are any nonzero scalar multiples of

If the entries of the matrix *A* are all real numbers, then the coefficients of the characteristic polynomial will also be real numbers, but the eigenvalues may still have nonzero imaginary parts. The entries of the corresponding eigenvectors therefore may also have nonzero imaginary parts. Similarly, the eigenvalues may be irrational numbers even if all the entries of *A* are rational numbers or even if they are all integers. However, if the entries of *A* are all algebraic numbers, which include the rationals, the eigenvalues are complex algebraic numbers.

The non-real roots of a real polynomial with real coefficients can be grouped into pairs of complex conjugates, namely with the two members of each pair having imaginary parts that differ only in sign and the same real part. If the degree is odd, then by the intermediate value theorem at least one of the roots is real. Therefore, any real matrix with odd order has at least one real eigenvalue, whereas a real matrix with even order may not have any real eigenvalues. The eigenvectors associated with these complex eigenvalues are also complex and also appear in complex conjugate pairs.

Algebraic multiplicity

Let *λ**i*be an eigenvalue of an* n*by* n*matrix* A*. The* **algebraic multiplicity***μ**A* (λ*imultiplicity as a root of the characteristic polynomial, that is, the largest integer k*such that (*λ* − *λ**i*)*k* divides evenly that polynomial.^{[8]}^{[29]}^{[30]}

Suppose a matrix *A* has dimension *n* and *d* ≤ *n* distinct eigenvalues. Whereas Equation (4) factors the characteristic polynomial of *A* into the product of *n* linear terms with some terms potentially repeating, the characteristic polynomial can instead be written as the product of *d* terms each corresponding to a distinct eigenvalue and raised to the power of the algebraic multiplicity,

If *d* = *n* then the right-hand side is the product of *n* linear terms and this is the same as Equation (4). The size of each eigenvalue's algebraic multiplicity is related to the dimension *n* as

If *μ**A*(*λ**i*) = 1, then* λ**i is said to be a simple eigenvalue*.^{[30]} If* μ**A*(

*λ**i*) equals the geometric multiplicity of* λ

**i*,* γ**A* (

*λ**i*), defined in the next section, then* λ**i

*is said to be a*semisimple eigenvalue*.

Eigenspaces, geometric multiplicity, and the eigenbasis for matrices

Given a particular eigenvalue *λ* of the *n* by *n* matrix *A*, define the set *E* to be all vectors * v* that satisfy Equation (2

On one hand, this set is precisely the kernel or nullspace of the matrix (*A* − *λI*). On the other hand, by definition, any nonzero vector that satisfies this condition is an eigenvector of *A* associated with *λ*. So, the set *E* is the union of the zero vector with the set of all eigenvectors of *A* associated with *λ*, and *E* equals the nullspace of (*A* − *λI*). *E* is called the **eigenspace** or **characteristic space** of *A* associated with *λ*.^{[31]}^{[8]} In general *λ* is a complex number and the eigenvectors are complex *n* by 1 matrices. A property of the nullspace is that it is a linear subspace, so *E* is a linear subspace of ℂ*n*.

Because the eigenspace *E* is a linear subspace, it is closed under addition. That is, if two vectors *u* and *v* belong to the set *E*, written (*u*,*v) ∈ E, then (u + v) ∈ E or equivalently A(u + v) = λ*s can be checked using the distributive property of matrix multiplication. Similarly, because *E* is a linear subspace, it is closed under scalar multiplication. That is, if *v* ∈ *E* and α is a complex number, (*αv*) ∈ *E* or equivalently *A*(α*v*) =* λ* (*αv*). This can be checked by noting that multiplication of complex matrices by complex numbers is commutative. As long as *u* + *v* and *αv* are not zero, they are also eigenvectors of *A* associated with *λ*.

The dimension of the eigenspace *E* associated with *λ*, or equivalently the maximum number of linearly independent eigenvectors associated with *λ*, is referred to as the eigenvalue's **geometric multiplicity** *γA(λ* Because *E* is also the nullspace of (*A* − *λI*), the geometric multiplicity of *λ* is the dimension of the nullspace of (*A* − *λI*), also called the *nullity* of (*A* − *λI*), which relates to the dimension and rank of (*A* − *λI*) as

Because of the definition of eigenvalues and eigenvectors, an eigenvalue's geometric multiplicity must be at least one, that is, each eigenvalue has at least one associated eigenvector.

Furthermore, an eigenvalue's geometric multiplicity cannot exceed its algebraic multiplicity.

Additionally, recall that an eigenvalue's algebraic multiplicity cannot exceed *n*.

The direct sum of the eigenspaces of all of 's eigenvalues is the entire vector space.

A basis of can be formed from linearly independent eigenvectors of; such a basis is called an

**eigenbasis**Any vector in can be written as a linear combination of eigenvectors of.

Additional properties of eigenvalues

The trace of, defined as the sum of its diagonal elements, is also the sum of all eigenvalues,

^{[32]}^{[33]}^{[34]}The determinant of is the product of all its eigenvalues,

^{[32]}^{[35]}^{[36]}The eigenvalues of the th power of; i.e., the eigenvalues of, for any positive integer, are.

The matrix is invertible if and only if every eigenvalue is nonzero.

If is invertible, then the eigenvalues of are and each eigenvalue's geometric multiplicity coincides.

Moreover, since the characteristic polynomial of the inverse is the reciprocal polynomial of the original, the eigenvalues share the same algebraic multiplicity.

If is equal to its conjugate transpose, or equivalently if is Hermitian, then every eigenvalue is real. The same is true of any symmetric real matrix.

If is not only Hermitian but also positive-definite, positive-semidefinite, negative-definite, or negative-semidefinite, then every eigenvalue is positive, non-negative, negative, or non-positive, respectively.

If is unitary, every eigenvalue has absolute value.

if is a matrix and are its eigenvalues, then the eigenvalues of matrix (where is the identity matrix) are.

Moreover, if, the eigenvalues of are.

More generally, for a polynomial the eigenvalues of matrix are.

Left and right eigenvectors

Diagonalization and the eigendecomposition

Suppose the eigenvectors of *A* form a basis, or equivalently *A* has *n* linearly independent eigenvectors *v*1, *v*2,..., *v**n*with associated eigenvalues λ*, *λ*2,..., *λ**n*. The eigenvalues need not be distinct. Define a square matrix* Q*whose columns are the* n*linearly independent eigenvectors of* A*,

Since each column of *Q* is an eigenvector of *A*, right multiplying *A* by *Q* scales each column of *Q* by its associated eigenvalue,

With this in mind, define a diagonal matrix Λ where each diagonal element Λ*ii*is the eigenvalue associated with the* i*th column of* Q*. Then

Because the columns of *Q* are linearly independent, Q is invertible. Right multiplying both sides of the equation by *Q*−1,

or by instead left multiplying both sides by *Q*−1,

*A* can therefore be decomposed into a matrix composed of its eigenvectors, a diagonal matrix with its eigenvalues along the diagonal, and the inverse of the matrix of eigenvectors. This is called the eigendecomposition and it is a similarity transformation. Such a matrix *A* is said to be *similar* to the diagonal matrix Λ or diagonalizable. The matrix *Q* is the change of basis matrix of the similarity transformation. Essentially, the matrices *A* and Λ represent the same linear transformation expressed in two different bases. The eigenvectors are used as the basis when representing the linear transformation as Λ.

Conversely, suppose a matrix *A* is diagonalizable. Let *P* be a non-singular square matrix such that *P*−1*AP* is some diagonal matrix *D*. Left multiplying both by *P*, *AP* = *PD*. Each column of *P* must therefore be an eigenvector of *A* whose eigenvalue is the corresponding diagonal element of *D*. Since the columns of *P* must be linearly independent for *P* to be invertible, there exist *n* linearly independent eigenvectors of *A*. It then follows that the eigenvectors of *A* form a basis if and only if *A* is diagonalizable.

A matrix that is not diagonalizable is said to be defective. For defective matrices, the notion of eigenvectors generalizes to generalized eigenvectors and the diagonal matrix of eigenvalues generalizes to the Jordan normal form. Over an algebraically closed field, any matrix *A* has a Jordan normal form and therefore admits a basis of generalized eigenvectors and a decomposition into generalized eigenspaces.

Variational characterization

Matrix examples

Two-dimensional matrix example

Consider the matrix

The figure on the right shows the effect of this transformation on point coordinates in the plane.

The eigenvectors *v* of this transformation satisfy Equation (1), and the values of *λ* for which the determinant of the matrix (*A* − *λI*) equals zero are the eigenvalues.

Taking the determinant to find characteristic polynomial of *A*,

Setting the characteristic polynomial equal to zero, it has roots at *λ* = 1 and *λ* = 3, which are the two eigenvalues of *A*.

For *λ* = 1, Equation (2

Any nonzero vector with *v*1 = −*v*2 solves this equation. Therefore,

is an eigenvector of *A* corresponding to *λ* = 1, as is any scalar multiple of this vector.

For *λ* = 3, Equation (2

Any nonzero vector with *v*1 = *v*2 solves this equation. Therefore,

is an eigenvector of *A* corresponding to *λ* = 3, as is any scalar multiple of this vector.

Thus, the vectors *vλ*=1 and *vλ*=3 are eigenvectors of *A* associated with the eigenvalues *λ* = 1 and *λ* = 3, respectively.

Three-dimensional matrix example

Consider the matrix

The characteristic polynomial of *A* is

Three-dimensional matrix example with complex eigenvalues

Consider the cyclic permutation matrix

This matrix shifts the coordinates of the vector up by one position and moves the first coordinate to the bottom.

Its characteristic polynomial is 1 − *λ*3, whose roots are

For the real eigenvalue *λ*1 = 1, any vector with three equal nonzero entries is an eigenvector. For example,

For the complex conjugate pair of imaginary eigenvalues,

Then

and

Diagonal matrix example

Matrices with entries only along the main diagonal are called diagonal matrices

The characteristic polynomial of *A* is

which has the roots *λ*1 = 1, *λ*2 = 2, and *λ*3 = 3. These roots are the diagonal elements as well as the eigenvalues of *A*.

Each diagonal element corresponds to an eigenvector whose only nonzero component is in the same row as that diagonal element.

In the example, the eigenvalues correspond to the eigenvectors,

respectively, as well as scalar multiples of these vectors.

Triangular matrix example

A matrix whose elements above the main diagonal are all zero is called a triangular matrix, while a matrix whose elements below the main diagonal are all zero is called an *upper triangular matrix*. As with diagonal matrices, the eigenvalues of triangular matrices are the elements of the main diagonal.

Consider the lower triangular matrix,

The characteristic polynomial of *A* is

which has the roots *λ*1 = 1, *λ*2 = 2, and *λ*3 = 3. These roots are the diagonal elements as well as the eigenvalues of *A*.

These eigenvalues correspond to the eigenvectors,

respectively, as well as scalar multiples of these vectors.

Matrix with repeated eigenvalues example

As in the previous example, the lower triangular matrix

has a characteristic polynomial that is the product of its diagonal elements,

The roots of this polynomial, and hence the eigenvalues, are 2 and 3.

The *algebraic multiplicity* of each eigenvalue is 2; in other words they are both double roots. The sum of the algebraic multiplicities of each distinct eigenvalue is *μ**A*= 4 =* n*, the order of the characteristic polynomial and the dimension of* A*.

Eigenvalues and eigenfunctions of differential operators

The definitions of eigenvalue and eigenvectors of a linear transformation *T* remains valid even if the underlying vector space is an infinite-dimensional Hilbert or Banach space. A widely used class of linear transformations acting on infinite-dimensional spaces are the differential operators on function spaces. Let *D* be a linear differential operator on the space **C**∞ of infinitely differentiable real functions of a real argument *t*. The eigenvalue equation for *D* is the differential equation

The functions that satisfy this equation are eigenvectors of *D* and are commonly called **eigenfunctions**.

Derivative operator example

This differential equation can be solved by multiplying both sides by *dt*/*f* (*t*) and integrating. Its solution, the exponential function

is the eigenfunction of the derivative operator.

In this case the eigenfunction is itself a function of its associated eigenvalue.

In particular, for *λ* = 0 the eigenfunction *f* (*t*) is a constant.

The main eigenfunction article gives other examples.

General definition

The concept of eigenvalues and eigenvectors extends naturally to arbitrary linear transformations on arbitrary vector spaces. Let *V* be any vector space over some field *K* of scalars, and let *T* be a linear transformation mapping *V* into *V*,

We say that a nonzero vector **v** ∈ *V* is an **eigenvector** of *T* if and only if there exists a scalar *λ* ∈ *K* such that

This equation is called the eigenvalue equation for *T*, and the scalar *λ* is the **eigenvalue** of *T* corresponding to the eigenvector **v**. *T* (**v**) is the result of applying the transformation *T* to the vector **v**, while *λ***v** is the product of the scalar *λ* with **v**.^{[37]}^{[38]}

Eigenspaces, geometric multiplicity, and the eigenbasis

Given an eigenvalue *λ*, consider the set

which is the union of the zero vector with the set of all eigenvectors associated with *λ*. *E* is called the **eigenspace** or **characteristic space** of *T* associated with *λ*.

By definition of a linear transformation,

for (**x**,**y**) ∈ *V* and α ∈ *K*. Therefore, if **u** and **v** are eigenvectors of *T* associated with eigenvalue *λ*, namely **u**,**v** ∈ *E*, then

So, both **u** + **v** and α**v** are either zero or eigenvectors of *T* associated with *λ*, namely **u** + **v**, α**v** ∈ *E*, and *E* is closed under addition and scalar multiplication. The eigenspace *E* associated with *λ* is therefore a linear subspace of *V*.^{[39]} If that subspace has dimension 1, it is sometimes called an **eigenline**.^{[40]}

The **geometric multiplicity** *γ**T*(λ* of an eigenvalue *λ* is the dimension of the eigenspace associated with *λ*, i.e., the maximum number of linearly independent eigenvectors associated with that eigenvalue.^{[8]}^{[30]} By the definition of eigenvalues and eigenvectors, *γ**T* (*λ*) ≥ 1 because every eigenvalue has at least one eigenvector.

The eigenspaces of *T* always form a direct sum. As a consequence, eigenvectors of *different* eigenvalues are always linearly independent. Therefore, the sum of the dimensions of the eigenspaces cannot exceed the dimension *n* of the vector space on which *T* operates, and there cannot be more than *n* distinct eigenvalues.^{[41]}

Any subspace spanned by eigenvectors of *T* is an invariant subspace of *T*, and the restriction of *T* to such a subspace is diagonalizable. Moreover, if the entire vector space *V* can be spanned by the eigenvectors of *T*, or equivalently if the direct sum of the eigenspaces associated with all the eigenvalues of *T* is the entire vector space *V*, then a basis of *V* called an **eigenbasis** can be formed from linearly independent eigenvectors of *T*. When *T* admits an eigenbasis, *T* is diagonalizable.

Zero vector as an eigenvector

While the definition of an eigenvector used in this article excludes the zero vector, it is possible to define eigenvalues and eigenvectors such that the zero vector is an eigenvector.^{[42]}

Consider again the eigenvalue equation, Equation (5). Define an **eigenvalue** to be any scalar *λ* ∈ *K* such that there exists a nonzero vector **v** ∈ *V* satisfying Equation (5). It is important that this version of the definition of an eigenvalue specify that the vector be nonzero, otherwise by this definition the zero vector would allow any scalar in *K* to be an eigenvalue. Define an **eigenvector** **v** associated with the eigenvalue *λ* to be any vector that, given *λ*, satisfies Equation (5). Given the eigenvalue, the zero vector is among the vectors that satisfy Equation (5

Spectral theory

If *λ* is an eigenvalue of *T*, then the operator (*T* − *λI*) is not one-to-one, and therefore its inverse (*T* − *λI*)−1 does not exist. The converse is true for finite-dimensional vector spaces, but not for infinite-dimensional vector spaces. In general, the operator (*T* − *λI*) may not have an inverse even if *λ* is not an eigenvalue.

For this reason, in functional analysis eigenvalues can be generalized to the spectrum of a linear operator *T* as the set of all scalars *λ* for which the operator (*T* − *λI*) has no bounded inverse. The spectrum of an operator always contains all its eigenvalues but is not limited to them.

Associative algebras and representation theory

One can generalize the algebraic object that is acting on the vector space, replacing a single operator acting on a vector space with an algebra representation – an associative algebra acting on a module. The study of such actions is the field of representation theory.

The representation-theoretical concept of weight is an analog of eigenvalues, while *weight vectors* and *weight spaces* are the analogs of eigenvectors and eigenspaces, respectively.

Dynamic equations

The simplest difference equations have the form

The solution of this equation for *x* in terms of *t* is found by using its characteristic equation

A similar procedure is used for solving a differential equation of the form

Calculation

The calculation of eigenvalues and eigenvectors is a topic where theory, as presented in elementary linear algebra textbooks, is often very far from practice.

Classical method

The classical method is to first find the eigenvalues, and then calculate the eigenvectors for each eigenvalue.

It is in several ways poorly suited for non-exact arithmetics such as floating-point.

Eigenvalues

Eigenvectors

Once the (exact) value of an eigenvalue is known, the corresponding eigenvectors can be found by finding nonzero solutions of the eigenvalue equation, that becomes a system of linear equations with known coefficients. For example, once it is known that 6 is an eigenvalue of the matrix

This matrix equation is equivalent to two linear equations

- that is

Simple iterative methods

Modern methods

Efficient, accurate methods to compute eigenvalues and eigenvectors of arbitrary matrices were not known until the QR algorithm was designed in 1961.^{[43]} Combining the Householder transformation with the LU decomposition results in an algorithm with better convergence than the QR algorithm. For large Hermitian sparse matrices, the Lanczos algorithm is one example of an efficient iterative method to compute eigenvalues and eigenvectors, among several other possibilities.^{[43]}

Most numeric methods that compute the eigenvalues of a matrix also determine a set of corresponding eigenvectors as a by-product of the computation, although sometimes implementors choose to discard the eigenvector information as soon as it is no longer needed.

Applications

Eigenvalues of geometric transformations

The following table presents some example transformations in the plane along with their 2×2 matrices, eigenvalues, and eigenvectors.

A linear transformation that takes a square to a rectangle of the same area (a squeeze mapping) has reciprocal eigenvalues.

Schrödinger equation

Molecular orbitals

In quantum mechanics, and in particular in atomic and molecular physics, within the Hartree–Fock theory, the atomic and molecular orbitals can be defined by the eigenvectors of the Fock operator. The corresponding eigenvalues are interpreted as ionization potentials via Koopmans' theorem. In this case, the term eigenvector is used in a somewhat more general meaning, since the Fock operator is explicitly dependent on the orbitals and their eigenvalues. Thus, if one wants to underline this aspect, one speaks of nonlinear eigenvalue problems. Such equations are usually solved by an iteration procedure, called in this case self-consistent field method. In quantum chemistry, one often represents the Hartree–Fock equation in a non-orthogonal basis set. This particular representation is a generalized eigenvalue problem called Roothaan equations.

Geology and glaciology

In geology, especially in the study of glacial till, eigenvectors and eigenvalues are used as a method by which a mass of information of a clast fabric's constituents' orientation and dip can be summarized in a 3-D space by six numbers. In the field, a geologist may collect such data for hundreds or thousands of clasts in a soil sample, which can only be compared graphically such as in a Tri-Plot (Sneed and Folk) diagram,^{[44]}^{[45]} or as a Stereonet on a Wulff Net.^{[46]}

Principal component analysis

The eigendecomposition of a symmetric positive semidefinite (PSD) matrix yields an orthogonal basis of eigenvectors, each of which has a nonnegative eigenvalue. The orthogonal decomposition of a PSD matrix is used in multivariate analysis, where the sample covariance matrices are PSD. This orthogonal decomposition is called principal components analysis (PCA) in statistics. PCA studies linear relations among variables. PCA is performed on the covariance matrix or the correlation matrix (in which each variable is scaled to have its sample variance equal to one). For the covariance or correlation matrix, the eigenvectors correspond to principal components and the eigenvalues to the variance explained by the principal components. Principal component analysis of the correlation matrix provides an orthonormal eigen-basis for the space of the observed data: In this basis, the largest eigenvalues correspond to the principal components that are associated with most of the covariability among a number of observed data.

Principal component analysis is used to study large data sets, such as those encountered in bioinformatics, data mining, chemical research, psychology, and in marketing. PCA is also popular in psychology, especially within the field of psychometrics. In Q methodology, the eigenvalues of the correlation matrix determine the Q-methodologist's judgment of *practical* significance (which differs from the statistical significance of hypothesis testing; cf. criteria for determining the number of factors). More generally, principal component analysis can be used as a method of factor analysis in structural equation modeling.

Vibration analysis

Eigenvalue problems occur naturally in the vibration analysis of mechanical structures with many degrees of freedom. The eigenvalues are the natural frequencies (or **eigenfrequencies**) of vibration, and the eigenvectors are the shapes of these vibrational modes. In particular, undamped vibration is governed by

or

leads to a so-called quadratic eigenvalue problem,

This can be reduced to a generalized eigenvalue problem by algebraic manipulation at the cost of solving a larger system.

The orthogonality properties of the eigenvectors allows decoupling of the differential equations so that the system can be represented as linear summation of the eigenvectors.

The eigenvalue problem of complex structures is often solved using finite element analysis, but neatly generalize the solution to scalar-valued vibration problems.

Eigenfaces

In image processing, processed images of faces can be seen as vectors whose components are the brightnesses of each pixel.^{[49]} The dimension of this vector space is the number of pixels. The eigenvectors of the covariance matrix associated with a large set of normalized pictures of faces are calledeigenfaces; this is an example of principal component analysis. They are very useful for expressing any face image as a linear combination of some of them. In the facial recognition branch of biometrics, eigenfaces provide a means of applying data compression to faces for identification purposes. Research related to eigen vision systems determining hand gestures has also been made.

Similar to this concept, **eigenvoices** represent the general direction of variability in human pronunciations of a particular utterance, such as a word in a language. Based on a linear combination of such eigenvoices, a new voice pronunciation of the word can be constructed. These concepts have been found useful in automatic speech recognition systems for speaker adaptation.

Tensor of moment of inertia

In mechanics, the eigenvectors of the moment of inertia tensor define the principal axes of a rigid body. The tensor of moment of inertia is a key quantity required to determine the rotation of a rigid body around its center of mass.

Stress tensor

In solid mechanics, the stress tensor is symmetric and so can be decomposed into a diagonal tensor with the eigenvalues on the diagonal and eigenvectors as a basis. Because it is diagonal, in this orientation, the stress tensor has no shear components; the components it does have are the principal components.

Graphs

The principal eigenvector is used to measure the centrality of its vertices. An example is Google's PageRank algorithm. The principal eigenvector of a modified adjacency matrix of the World Wide Web graph gives the page ranks as its components. This vector corresponds to the stationary distribution of the Markov chain represented by the row-normalized adjacency matrix; however, the adjacency matrix must first be modified to ensure a stationary distribution exists. The second smallest eigenvector can be used to partition the graph into clusters, via spectral clustering. Other methods are also available for clustering.

Basic reproduction number

See also

Antieigenvalue theory

Eigenoperator

Eigenplane

Eigenvalue algorithm

Introduction to eigenstates

Jordan normal form

List of numerical analysis software

Nonlinear eigenproblem

Quadratic eigenvalue problem

Singular value

## References

*Topics In Algebra*, Waltham: Blaisdell Publishing Company, ISBN 978-1114541016, pp. 228, 229.

*Linear Algebra and Matrix Theory*(2nd ed.), New York: Wiley, LCCN 76091646, p. 38.

*Numerical Analysis*(5th ed.), Boston: Prindle, Weber and Schmidt, ISBN 0-534-93219-3, p. 401.

*The New Cassell's German Dictionary*, New York: Funk & Wagnall, LCCN 58-7924.

*mathworld.wolfram.com*. Retrieved 2019-08-04.

*Histoire de l'Académie royale des sciences et des belles lettres de Berlin*, pp. 176–227. On p. 212, Euler proves that any body contains a principal axis of rotation:

*"Théorem. 44. De quelque figure que soit le corps, on y peut toujours assigner un tel axe, qui passe par son centre de gravité, autour duquel le corps peut tourner librement & d'un mouvement uniforme."*(Theorem. 44. Whatever be the shape of the body, one can always assign to it such an axis, which passes through its center of gravity, around which it can rotate freely and with a uniform motion.) In 1755, Johann Andreas Segner proved that any body has three principal axes of rotation: Johann Andreas Segner,

*Specimen theoriae turbinum*[Essay on the theory of tops (i.e., rotating bodies)] ( Halle ("Halae"), (Germany) : Gebauer, 1755). p. xxviiii [29], Segner derives a third-degree equation in

*t*, which proves that a body has three principal axes of rotation. He then states (on the same page):

*"Non autem repugnat tres esse eiusmodi positiones plani HM, quia in aequatione cubica radices tres esse possunt, et tres tangentis t valores."*(However, it is not inconsistent [that there] be three such positions of the plane HM, because in cubic equations, [there] can be three roots, and three values of the tangent t.) The relevant passage of Segner's work was discussed briefly by Arthur Cayley. See: A. Cayley (1862) "Report on the progress of the solution of certain special problems of dynamics,"

*Report of the Thirty-second meeting of the British Association for the Advancement of Science; held at Cambridge in October 1862*,

**32**: 184–252 ; see especially 225–226.

*Historia Mathematica*,

**2**: 1–29, doi:10.1016/0315-0860(75)90032-4, §2.

*Mathematical thought from ancient to modern times*, Oxford University Press, ISBN 0-19-501496-0, pp. 807–808 Augustin Cauchy (1839) "Mémoire sur l'intégration des équations linéaires" (Memoir on the integration of linear equations),

*Comptes rendus*,

**8**: 827–830, 845–865, 889–907, 931–937. From p. 827:

*"On sait d'ailleurs qu'en suivant la méthode de Lagrange, on obtient pour valeur générale de la variable prinicipale une fonction dans laquelle entrent avec la variable principale les racines d'une certaine équation que j'appellerai l'*(One knows, moreover, that by following Lagrange's method, one obtains for the general value of the principal variable a function in which there appear, together with the principal variable, the roots of a certain equation that I will call the "characteristic equation", the degree of this equation being precisely the order of the differential equation that must be integrated.)

*équation caractéristique*, le degré de cette équation étant précisément l'order de l'équation différentielle qu'il s'agit d'intégrer."*Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse*(News of the Philosophical Society at Göttingen, mathematical-physical section), pp. 49–91. From p. 51:

*"Insbesondere in dieser ersten Mitteilung gelange ich zu Formeln, die die Entwickelung einer willkürlichen Funktion nach gewissen ausgezeichneten Funktionen, die ich*Eigenfunktionen

*nenne, liefern: …*(In particular, in this first report I arrive at formulas that provide the [series] development of an arbitrary function in terms of some distinctive functions, which I call

*eigenfunctions*: … ) Later on the same page:

*"Dieser Erfolg ist wesentlich durch den Umstand bedingt, daß ich nicht, wie es bisher geschah, in erster Linie auf den Beweis für die Existenz der Eigenwerte ausgehe, … "*(This success is mainly attributable to the fact that I do not, as it has happened until now, first of all aim at a proof of the existence of eigenvalues, … ) For the origin and evolution of the terms eigenvalue, characteristic value, etc., see: Earliest Known Uses of Some of the Words of Mathematics (E)

*Earliest Known Uses of Some of the Words of Mathematics*, retrieved 2006-08-22.