Positive-definite matrix
From Wikipedia, the free encyclopedia
In linear algebra, a positive-definite matrix is a Hermitian matrix which in many ways is analogous to a positive real number. The notion is closely related to a positive-definite symmetric bilinear form (or a sesquilinear form in the complex case).
Contents |
[edit] Equivalent formulations
Let M be an n × n Hermitian matrix. Denote the transpose of a vector a by aT, and the conjugate transpose by a * .
A matrix M can be tested for positive definiteness by looking for one of the following properties:
1. | For all non-zero vectors z ∈ Cn,
Note that the quantity z * Mz is always real. |
2. | All eigenvalues λi of M are positive. Recall that any Hermitian M, by the spectral theorem, may be regarded as a real diagonal matrix D that has been re-expressed in some new coordinate system (i.e., M = P − 1DP for some unitary matrix P whose rows are orthonormal eigenvectors of M, forming a basis). So this characterization means that M is positive definite if and only if the diagonal elements of D (the eigenvalues) are all positive. In other words, in the basis consisting of its eigenvectors of M, the action of M is component-wise multiplication with a (fixed) element in Cn with positive entries. |
3. | The sesquilinear form
defines an inner product on Cn. (In fact, every inner product on Cn arises in this fashion from a Hermitian positive definite matrix.) |
4. | M is the Gram matrix of some collection of linearly independent vectors
for some k. More precisely, M arises by defining each entry The vectors xi may optionally be restricted to fall in Cn. In other words, M is of the form A*A where A is not necessarily square but must be injective in general. |
5. | All the following matrices (the leading principal minors) have a positive determinant (the Sylvester criterion):
|
These properties are equivalent: if a matrix satisfies one property, it satisfies them all.
For real symmetric matrices, these properties can be simplified by replacing with , and "conjugate transpose" with "transpose."
[edit] Quadratic forms
Echoing condition 3 above, one can also formulate positive-definiteness in terms of quadratic forms. Let K be the field R or C, and V be a vector space over K. A Hermitian form
is a bilinear map such that B(x, y) is always the complex conjugate of B(y, x). Such a function B is called positive definite if B(x, x) > 0 for every nonzero x in V.
[edit] Further properties
Every positive definite matrix is invertible and its inverse is also positive definite. If M is positive definite and r > 0 is a real number, then rM is positive definite. If M and N are positive definite, then the sum M + N and the products MNM and NMN are also positive definite; and if MN = NM, then MN is also positive definite. Every positive definite matrix M, has at least one square root matrix N such that N2 = M. In fact, M may have infinitely many square roots, but exactly one positive definite square root.
[edit] Negative-definite, semidefinite and indefinite matrices
The n × n Hermitian matrix M is said to be negative-definite if
for all non-zero (or, equivalently, all non-zero ). It is called positive-semidefinite if
for all (or ) and negative-semidefinite if
for all (or ).
Equivalently, a matrix is negative-definite if all its eigenvalues are negative, it is positive-semidefinite if they are all greater than or equal to zero, and it is negative-semidefinite if they are all less than or equal to zero.
A matrix M is positive-semidefinite if and only if it arises as the Gram matrix of some set of vectors. In contrast to the positive-definite case, these vectors need not be linearly independent.
For any matrix A, the matrix A*A is positive semidefinite, and rank(A) = rank(A*A). Reversely, any positive semidefinite matrix M can be written as M = A*A; this is the Cholesky decomposition.
A Hermitian matrix which is neither positive- nor negative-semidefinite is called indefinite.
[edit] Non-Hermitian matrices
A real matrix M may have the property that xTMx > 0 for all nonzero real vectors x without being symmetric. The matrix
provides an example. In general, we have xTMx > 0 for all real nonzero vectors x if and only if the symmetric part, (M + MT) / 2, is positive definite.
The situation for complex matrices may be different, depending on how one generalizes the inequality z*Mz > 0. If z*Mz is real for all complex vectors z, then the matrix M is necessarily Hermitian. So, if we require that z*Mz be real and positive, then M is automatically Hermitian. On the other hand, we have that Re(z*Mz) > 0 for all complex nonzero vectors z if and only if the Hermitian part, (M + M*) / 2, is positive definite.
In summary, the distinguishing feature between the real and complex case is that, a bounded positive operator on a complex Hilbert space is necessarily Hermitian, or self adjoint. The general claim can be argued using the polarization identity. That is no longer true in the real case.
There is no agreement in the literature on the proper definition of positive-definite for non-Hermitian matrices.
[edit] See also
[edit] References
- Roger A. Horn and Charles R. Johnson. Matrix Analysis, Chapter 7. Cambridge University Press, 1985. ISBN 0-521-30586-1 (hardback), ISBN 0-521-38632-2 (paperback).