In linear algebra, a square matrix is called diagonalizable or nondefective if it is similar to a diagonal matrix, i.e., if there exists an invertible matrix such that is a diagonal matrix. If is a finite-dimensional vector space, then a linear map is called diagonalizable if there exists an ordered basis of with respect to which is represented by a diagonal matrix. Diagonalization is the process of finding a corresponding diagonal matrix for a diagonalizable matrix or linear map. A square matrix that is not diagonalizable is called defective.
Diagonalizable matrices and maps are of interest because diagonal matrices are especially easy to handle; once their eigenvalues and eigenvectors are known, one can raise a diagonal matrix to a power by simply raising the diagonal entries to that same power, and the determinant of a diagonal matrix is simply the product of all diagonal entries. Geometrically, a diagonalizable matrix is an inhomogeneous dilation (or anisotropic scaling) — it scales the space, as does a homogeneous dilation, but by a different factor in each direction, determined by the scale factors on each axis (diagonal entries).
A square matrix over a field is called diagonalizable or nondefective if there exists an invertible matrix such that is a diagonal matrix. Formally,
The fundamental fact about diagonalizable maps and matrices is expressed by the following:
- An matrix over a field is diagonalizable if and only if the sum of the dimensions of its eigenspaces is equal to , which is the case if and only if there exists a basis of consisting of eigenvectors of . If such a basis has been found, one can form the matrix having these basis vectors as columns, and will be a diagonal matrix. The diagonal entries of this matrix are the eigenvalues of .
- A linear map is diagonalizable if and only if the sum of the dimensions of its eigenspaces is equal to , which is the case if and only if there exists a basis of consisting of eigenvectors of . With respect to such a basis, will be represented by a diagonal matrix. The diagonal entries of this matrix are the eigenvalues of .
Another characterization: A matrix or linear map is diagonalizable over the field if and only if its minimal polynomial is a product of distinct linear factors over . (Put in another way, a matrix is diagonalizable if and only if all of its elementary divisors are linear.)
The following sufficient (but not necessary) condition is often useful.
- An matrix is diagonalizable over the field if it has distinct eigenvalues in , i.e. if its characteristic polynomial has distinct roots in ; however, the converse may be false. Consider
which has eigenvalues 1, 2, 2 (not all distinct) and is diagonalizable with diagonal form (similar to )
- A linear map with is diagonalizable if it has distinct eigenvalues, i.e. if its characteristic polynomial has distinct roots in .
Let be a matrix over . If is diagonalizable, then so is any power of it. Conversely, if is invertible, is algebraically closed, and is diagonalizable for some that is not an integer multiple of the characteristic of , then is diagonalizable. Proof: If is diagonalizable, then is annihilated by some polynomial , which has no multiple root (since ) and is divided by the minimal polynomial of .
As a rule of thumb, over almost every matrix is diagonalizable. More precisely: the set of complex matrices that are not diagonalizable over , considered as a subset of , has Lebesgue measure zero. One can also say that the diagonalizable matrices form a dense subset with respect to the Zariski topology: the complement lies inside the set where the discriminant of the characteristic polynomial vanishes, which is a hypersurface. From that follows also density in the usual (strong) topology given by a norm. The same is not true over .
The Jordan–Chevalley decomposition expresses an operator as the sum of its semisimple (i.e., diagonalizable) part and its nilpotent part. Hence, a matrix is diagonalizable if and only if its nilpotent part is zero. Put in another way, a matrix is diagonalizable if each block in its Jordan form has no nilpotent part; i.e., each "block" is a one-by-one matrix.
If a matrix can be diagonalized, that is,
Writing as a block matrix of its column vectors
the above equation can be rewritten as
So the column vectors of are right eigenvectors of , and the corresponding diagonal entry is the corresponding eigenvalue. The invertibility of also suggests that the eigenvectors are linearly independent and form a basis of . This is the necessary and sufficient condition for diagonalizability and the canonical approach of diagonalization. The row vectors of are the left eigenvectors of .
When a complex matrix [a] is a Hermitian matrix (or a real matrix[b], a symmetric matrix), eigenvectors of can be chosen to form an orthonormal basis of (or for a real matrix). Under such circumstance will be a unitary matrix (resp. orthogonal matrix) and equals the conjugate transpose of (if real, then the transpose of ).
For most practical work matrices are diagonalized numerically using computer software. Many algorithms exist to accomplish this.
A set of matrices is said to be simultaneously diagonalizable if there exists a single invertible matrix such that is a diagonal matrix for every in the set. The following theorem characterises simultaneously diagonalisable matrices: A set of diagonalizable matrices commutes if and only if the set is simultaneously diagonalizable.:pp. 51-53
The set of all diagonalizable matrices (over ) with is not simultaneously diagonalizable. For instance, the matrices
are diagonalizable but not simultaneously diagonalizable because they do not commute.
- Involutions are diagonalizable over the reals (and indeed any field of characteristic not 2), with ±1 on the diagonal.
- Finite order endomorphisms are diagonalisable over (or any algebraically closed field where the characteristic of the field does not divide the order of the endomorphism) with roots of unity on the diagonal. This follows since the minimal polynomial is separable, because the roots of unity are distinct.
- Projections are diagonalizable, with 0s and 1s on the diagonal.
- Real symmetric matrices are diagonalizable by orthogonal matrices; i.e., given a real symmetric matrix , is diagonal for some orthogonal matrix . More generally, matrices are diagonalizable by unitary matrices if and only if they are normal. In the case of the real symmetric matrix, we see that , so clearly holds. Examples of normal matrices are real symmetric (or skew-symmetric) matrices (e.g. covariance matrices) and Hermitian matrices (or skew-Hermitian matrices). See spectral theorems for generalizations to infinite-dimensional vector spaces.
Matrices that are not diagonalizableEdit
In general, a rotation matrix is not diagonalizable over the reals, but all rotation matrices are diagonalizable over the complex field. Even if a matrix is not diagonalizable, it is always possible to "do the best one can", and find a matrix with the same properties consisting of eigenvalues on the leading diagonal, and either ones or zeroes on the superdiagonal – known as Jordan normal form.
Some matrices are not diagonalizable over any field, most notably nonzero nilpotent matrices. This happens more generally if the algebraic and geometric multiplicities of an eigenvalue do not coincide. For instance, consider
This matrix is not diagonalizable: there is no matrix such that is a diagonal matrix. Indeed, has one eigenvalue (namely zero) and this eigenvalue has algebraic multiplicity 2 and geometric multiplicity 1.
Some real matrices are not diagonalizable over the reals. Consider for instance the matrix
The matrix does not have any real eigenvalues, so there is no real matrix such that is a diagonal matrix. However, we can diagonalize if we allow complex numbers. Indeed, if we take
then is diagonal. It is easy to find that B is the rotation matrix which rotates counterclockwise by angle
Note that the above examples show that the sum of diagonalizable matrices need not be diagonalizable.
How to diagonalize a matrixEdit
Consider a matrix
This matrix has eigenvalues
is a matrix with 3 different eigenvalues; therefore, it is diagonalizable. Note that if there are exactly n distinct eigenvalues in an matrix then this matrix is diagonalizable.
These eigenvalues are the values that will appear in the diagonalized form of matrix A, so by finding the eigenvalues of A we have diagonalized it. We could stop here, but it is a good check to use the eigenvectors to diagonalize A.
The eigenvectors of are
One can easily check that
Now, let be the matrix with these eigenvectors as its columns:
Then diagonalizes , as a simple computation confirms, having calculated using any suitable method:
In fact, this follows abstractly from the fact that for the standard basis we have
where we make use of the fact that is the k-th column of and hence . Note that the eigenvalues appear in the diagonal matrix.
Diagonalization can be used to compute the powers of a matrix efficiently, provided the matrix is diagonalizable. Suppose we have found that
where is a diagonal matrix. Then, as the matrix product is associative,
and the latter is easy to calculate since it only involves the powers of a diagonal matrix. This approach can be generalized to matrix exponential and other matrix functions that can be defined as power series.
For example, consider the following matrix:
Calculating the various powers of reveals a surprising pattern:
The above phenomenon can be explained by diagonalizing . To accomplish this, we need a basis of consisting of eigenvectors of . One such eigenvector basis is given by
where ei denotes the standard basis of Rn. The reverse change of basis is given by
Straightforward calculations show that
Thus, a and b are the eigenvalues corresponding to u and v, respectively. By linearity of matrix multiplication, we have that
Switching back to the standard basis, we have
The preceding relations, expressed in matrix form, are
thereby explaining the above phenomenon.
Quantum mechanical applicationEdit
In quantum mechanical and quantum chemical computations matrix diagonalization is one of the most frequently applied numerical processes. The basic reason is that the time-independent Schrödinger equation is an eigenvalue equation, albeit in most of the physical situations on an infinite dimensional space (a Hilbert space).
A very common approximation is to truncate Hilbert space to finite dimension, after which the Schrödinger equation can be formulated as an eigenvalue problem of a real symmetric, or complex Hermitian, matrix. Formally this approximation is founded on the variational principle, valid for Hamiltonians that are bounded from below.
First-order perturbation theory also leads to matrix eigenvalue problem for degenerate states.