Diagonal matrix

Vector-to-matrix diag operator

A diagonal matrix D can be constructed from a vector ${\displaystyle \mathbf{a}={\left[\begin{array}{ccc}{a}_1& \dots & a_n\end{array}\right]}^{\mathsf{\text{T}}}}$ using the ${\displaystyle \mathrm{diag}}$ operator: $${\displaystyle \mathbf{D}=\mathrm{diag}(a_1,\dots ,a_n).}$$

This may be written more compactly as ${\displaystyle \mathbf{D}=\mathrm{diag}(\mathbf{a})}$.

The same operator is also used to represent block diagonal matrices as ${\displaystyle \mathbf{A}=\mathrm{diag}({\mathbf{A}}_1,\dots ,{\mathbf{A}}_n)}$ where each argument A_i is a matrix.

The diag operator may be written as $${\displaystyle \mathrm{diag}(\mathbf{a})=\left(\mathbf{a}{\mathbf{1}}^{\mathsf{\text{T}}}\right)\circ \mathbf{I},}$$ where ${\displaystyle \circ }$ represents the Hadamard product, and 1 is a constant vector with elements 1.

Matrix-to-vector diag operator

The inverse matrix-to-vector diag operator is sometimes denoted by the identically named ${\displaystyle \mathrm{diag}(\mathbf{D})={\left[\begin{array}{ccc}{a}_1& \dots & a_n\end{array}\right]}^{\mathsf{\text{T}}},}$ where the argument is now a matrix, and the result is a vector of its diagonal entries.

The following property holds: $${\displaystyle \mathrm{diag}(\mathbf{A}\mathbf{B})=\sum _j{\left(\mathbf{A}\circ {\mathbf{B}}^{\mathsf{\text{T}}}\right)}_{ij}=\left(\mathbf{A}\circ {\mathbf{B}}^{\mathsf{\text{T}}}\right)\mathbf{1}.}$$

Scalar matrix

A diagonal matrix with equal diagonal entries is a scalar matrix; that is, a scalar multiple λ of the identity matrix I. Its effect on a vector is scalar multiplication by λ. For example, a 3×3 scalar matrix has the form: $${\displaystyle \left[\begin{array}{ccc}\lambda & 0& 0\\ 0& \lambda & 0\\ 0& 0& \lambda \end{array}\right]\equiv \lambda {\mathit{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. By contrast, over a field (like the real numbers), a diagonal matrix with all diagonal elements distinct only commutes with diagonal matrices (its centralizer is the set of diagonal matrices). That is because if a diagonal matrix ${\displaystyle \mathbf{D}=\mathrm{diag}(a_1,\dots ,a_n)}$ has ${\displaystyle a_i\ne a_j,}$ then given a matrix M with ${\displaystyle m_{ij}\ne 0,}$ the (i, j) term of the products are: ${\displaystyle (\mathbf{D}\mathbf{M}{)}_{ij}=a_i{m}_{ij}}$ and ${\displaystyle (\mathbf{M}\mathbf{D}{)}_{ij}=m_{ij}{a}_j,}$ and ${\displaystyle a_j{m}_{ij}\ne m_{ij}{a}_i}$ (since one can divide by m_ij), so they do not commute unless the off-diagonal terms are zero. Diagonal matrices where the diagonal entries are not all equal or all distinct have centralizers intermediate between the whole space and only diagonal matrices.

For an abstract vector space V (rather than the concrete vector space Kⁿ), the analog of scalar matrices are scalar transformations. This is true more generally for a module M over a ring R, with the endomorphism algebra End(M) (algebra of linear operators on M) replacing the algebra of matrices. Formally, scalar multiplication is a linear map, inducing a map ${\displaystyle R\to \mathrm{End}(M),}$ (from a scalar λ to its corresponding scalar transformation, multiplication by λ) exhibiting End(M) as a R-algebra. For vector spaces, the scalar transforms are exactly the center of the endomorphism algebra, and, similarly, scalar invertible transforms are the center of the general linear group GL(V). The former is more generally true free modules ${\displaystyle M\cong R^n,}$ for which the endomorphism algebra is isomorphic to a matrix algebra.

Vector operations

Multiplying a vector by a diagonal matrix multiplies each of the terms by the corresponding diagonal entry. Given a diagonal matrix ${\displaystyle \mathbf{D}=\mathrm{diag}(a_1,\dots ,a_n)}$ and a vector ${\displaystyle \mathbf{v}={\left[\begin{array}{ccc}{x}_1& \cdots & x_n\end{array}\right]}^{\mathsf{\text{T}}}}$, the product is: $${\displaystyle \mathbf{D}\mathbf{v}=\mathrm{diag}(a_1,\dots ,a_n)\left[\begin{array}{c}{x}_1\\ ⋮\\ x_n\end{array}\right]=\left[\begin{array}{c}{a}_1\\ & \ddots \\ & & a_n\end{array}\right]\left[\begin{array}{c}{x}_1\\ ⋮\\ x_n\end{array}\right]=\left[\begin{array}{c}{a}_1{x}_1\\ ⋮\\ a_n{x}_n\end{array}\right].}$$

This can be expressed more compactly by using a vector instead of a diagonal matrix, ${\displaystyle \mathbf{d}={\left[\begin{array}{ccc}{a}_1& \cdots & a_n\end{array}\right]}^{\mathsf{\text{T}}}}$, and taking the Hadamard product of the vectors (entrywise product), denoted ${\displaystyle \mathbf{d}\circ \mathbf{v}}$:

$${\displaystyle \mathbf{D}\mathbf{v}=\mathbf{d}\circ \mathbf{v}=\left[\begin{array}{c}{a}_1\\ ⋮\\ a_n\end{array}\right]\circ \left[\begin{array}{c}{x}_1\\ ⋮\\ x_n\end{array}\right]=\left[\begin{array}{c}{a}_1{x}_1\\ ⋮\\ a_n{x}_n\end{array}\right].}$$

This is mathematically equivalent, but avoids storing all the zero terms of this sparse matrix. This product is thus used in machine learning, such as computing products of derivatives in backpropagation or multiplying IDF weights in TF-IDF, since some BLAS frameworks, which multiply matrices efficiently, do not include Hadamard product capability directly.

Matrix operations

The operations of matrix addition and matrix multiplication are especially simple for 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

$${\displaystyle \mathrm{diag}(a_1,\dots ,a_n)+\mathrm{diag}(b_1,\dots ,b_n)=\mathrm{diag}(a_1+b_1,\dots ,a_n+b_n)}$$

and for matrix multiplication,

$${\displaystyle \mathrm{diag}(a_1,\dots ,a_n)\mathrm{diag}(b_1,\dots ,b_n)=\mathrm{diag}(a_1{b}_1,\dots ,a_n{b}_n).}$$

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

$${\displaystyle \mathrm{diag}(a_1,\dots ,a_n{)}^{-1}=\mathrm{diag}(a_1^{-1},\dots ,a_n^{-1}).}$$

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 i-th 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 i-th column of A by a_i for all i.

Operator matrix in eigenbasis

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 ${\mathbf{A}\mathbf{e}}_j=\sum _i{a}_{i,j}{\mathbf{e}}_i$, all coefficients a_{i, j} with i ≠ j are zero, leaving only one term per sum. The surviving diagonal elements, a_{i, j}, are known as eigenvalues and designated with λ_i in the equation, which reduces to ${\displaystyle {\mathbf{A}\mathbf{e}}_i={\lambda }_i{\mathbf{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.
• Where all matrices are square,
  • A matrix is diagonal if and only if it is triangular and normal.
  • A matrix is diagonal if and only if it is both upper- and lower-triangular.
  • A diagonal matrix is symmetric.
• The identity matrix I_n and zero matrix are diagonal.
• A 1×1 matrix is always diagonal.
• The square of a 2×2 matrix with zero trace is always diagonal.
• The inverse of a diagonal matrix is obtained by taking the reciprocals of the elements on the 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 U^∗AV 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.