Diagonal matrix

In linear algebra, a diagonal matrix is a matrix in which the entries outside the main diagonal are all zero. The term usually refers to square matrices. An example of a 2-by-2 diagonal matrix is ${\begin{bmatrix}3&0\\0&2\end{bmatrix}}$; the following matrix is a 3-by-3 diagonal matrix:${\begin{bmatrix}6&0&0\\0&7&0\\0&0&19\end{bmatrix}}$. An identity matrix of any size, or any multiple of it, will be a diagonal matrix.

Background

As stated above, the off-diagonal entries are zero. That is, the matrix D = (d_i,j) with n columns and n rows is diagonal if

$\forall i,j\in \{1,2,\ldots ,n\},i\neq j\implies d_{i,j}=0$.

However, the main diagonal entries are unrestricted.

Rectangular diagonal matrices

The term diagonal matrix may sometimes refer to a rectangular diagonal matrix, which is an m-by-n matrix with all the entries not of the form d_i,i being zero. For example:

${\begin{bmatrix}1&0&0\\0&4&0\\0&0&-3\\0&0&0\\\end{bmatrix}}$ or ${\begin{bmatrix}1&0&0&0&0\\0&4&0&0&0\\0&0&-3&0&0\end{bmatrix}}$

Symmetric diagonal matrices

The following matrix is a symmetric diagonal matrix:

${\begin{bmatrix}1&0&0\\0&4&0\\0&0&-2\end{bmatrix}}$

If the entries are real numbers or complex numbers, then it is a normal matrix as well.

In the remainder of this article we will consider only square matrices.

Scalar matrix

A square diagonal matrix with all its main diagonal entries equal is a scalar matrix, that is, a scalar multiple λI of the identity matrix I. Its effect on a vector is scalar multiplication by λ. For example, a 3×3 scalar matrix has the form:

${\begin{bmatrix}\lambda &0&0\\0&\lambda &0\\0&0&\lambda \end{bmatrix}}\equiv \lambda {\boldsymbol {I}}_{3}$

The scalar matrices are the center of the algebra of matrices: that is, they are precisely the matrices that commute with all other square matrices of the same size. All other diagonal matrices which are not scalar only commute with other diagonal matrices and not with any matrix unlike scalar matrices. Intuitively, this stems from the fact that scalar matrices are Identity matrices multiplied with scalars.

For an abstract vector space V (rather than the concrete vector space $K^{n}$), or more generally a module M over a ring R, with the endomorphism algebra End(M) (algebra of linear operators on M) replacing the algebra of matrices, the analog of scalar matrices are scalar transformations. Formally, scalar multiplication is a linear map, inducing a map $R\to \operatorname {End} (M),$ (send a scalar λ to the corresponding scalar transformation, multiplication by λ) exhibiting End(M) as a R-algebra. For vector spaces, or more generally free modules $M\cong R^{n}$, for which the endomorphism algebra is isomorphic to a matrix algebra, the scalar transforms are exactly the center of the endomorphism algebra, and similarly invertible transforms are the center of the general linear group GL(V), where they are denoted by Z(V), follow the usual notation for the center.

Matrix operations

The operations of matrix addition and matrix multiplication are especially simple for symmetric diagonal matrices. Write diag(a₁, ..., a_n) for a diagonal matrix whose diagonal entries starting in the upper left corner are a₁, ..., a_n. Then, for addition, we have

diag(a₁, ..., a_n) + diag(b₁, ..., b_n) = diag(a₁ + b₁, ..., a_n + b_n)

and for matrix multiplication,

diag(a₁, ..., a_n) · diag(b₁, ..., b_n) = diag(a₁b₁, ..., a_nb_n).

The diagonal matrix diag(a₁, ..., a_n) is invertible if and only if the entries a₁, ..., a_n are all non-zero. In this case, we have

diag(a₁, ..., a_n)⁻¹ = diag(a₁⁻¹, ..., a_n⁻¹).

In particular, the diagonal matrices form a subring of the ring of all n-by-n matrices.

Multiplying an n-by-n matrix A from the left with diag(a₁, ..., a_n) amounts to multiplying the ith row of A by a_i for all i; multiplying the matrix A from the right with diag(a₁, ..., a_n) amounts to multiplying the ith column of A by a_i for all i.

Operator matrix in eigenbasis

Main articles: Finding the matrix of a transformation and Eigenvalues and eigenvectors

As explained in determining coefficients of operator matrix, there is a special basis, e₁, ..., e_n, for which the matrix $A$ takes the diagonal form. Hence, in the defining equation $A{\vec {e}}_{j}=\sum a_{i,j}{\vec {e}}_{i}$, all coefficients $a_{i,j}$ with i ≠ j are zero, leaving only one term per sum. The surviving diagonal elements, $a_{i,i}$, are known as eigenvalues and designated with $\lambda _{i}$ in the equation, which reduces to $A{\vec {e}}_{i}=\lambda _{i}{\vec {e}}_{i}$. The resulting equation is known as eigenvalue equation and used to derive the characteristic polynomial and, further, eigenvalues and eigenvectors.

In other words, the eigenvalues of diag(λ₁, ..., λ_n) are λ₁, ..., λ_n with associated eigenvectors of e₁, ..., e_n.

Properties

The determinant of diag(a₁, ..., a_n) is the product a₁...a_n.

The adjugate of a diagonal matrix is again diagonal.

A square matrix is diagonal if and only if it is triangular and normal.

Any square diagonal matrix is also a symmetric matrix.

A symmetric diagonal matrix can be defined as a matrix that is both upper- and lower-triangular. The identity matrix I_n and any square zero matrix are diagonal. A one-dimensional matrix is always diagonal.

Applications

Diagonal matrices occur in many areas of linear algebra. Because of the simple description of the matrix operation and eigenvalues/eigenvectors given above, it is typically desirable to represent a given matrix or linear map by a diagonal matrix.

In fact, a given n-by-n matrix A is similar to a diagonal matrix (meaning that there is a matrix X such that X⁻¹AX is diagonal) if and only if it has n linearly independent eigenvectors. Such matrices are said to be diagonalizable.

Over the field of real or complex numbers, more is true. The spectral theorem says that every normal matrix is unitarily similar to a diagonal matrix (if AA^∗ = A^∗A then there exists a unitary matrix U such that UAU^∗ is diagonal). Furthermore, the singular value decomposition implies that for any matrix A, there exist unitary matrices U and V such that UAV^∗ is diagonal with positive entries.

Operator theory

In operator theory, particularly the study of PDEs, operators are particularly easy to understand and PDEs easy to solve if the operator is diagonal with respect to the basis with which one is working; this corresponds to a separable partial differential equation. Therefore, a key technique to understanding operators is a change of coordinates–in the language of operators, an integral transform–which changes the basis to an eigenbasis of eigenfunctions: which makes the equation separable. An important example of this is the Fourier transform, which diagonalizes constant coefficient differentiation operators (or more generally translation invariant operators), such as the Laplacian operator, say, in the heat equation.

Especially easy are multiplication operators, which are defined as multiplication by (the values of) a fixed function–the values of the function at each point correspond to the diagonal entries of a matrix.

See also

• Anti-diagonal matrix
• Banded matrix
• Bidiagonal matrix
• Diagonally dominant matrix
• Jordan normal form
• Multiplication operator
• Tridiagonal matrix
• Toeplitz matrix
• Toral Lie algebra
• Circulant matrix

References

• Roger A. Horn and Charles R. Johnson, Matrix Analysis, Cambridge University Press, 1985. ISBN 0-521-30586-1 (hardback), ISBN 0-521-38632-2 (paperback).