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).
Equivalent formulations
Let ''M'' be an ''n'' × ''n''
Hermitian matrix. Denote the transpose of a vector
by
, and the
conjugate transpose by
.
The matrix ''M'' is ''positive definite'' if and only if it satisfies any of the following properties:
{| cellspacing="0" cellpadding="2"
|-
|valign="top"| '1.' || For all non-zero vectors ''z'' ∈ 'C'
''n'',
:
.
Note that the quantity
is always real.
|-
|valign="top"| '2.' || All
eigenvalues
of
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.,
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 the eigenvectors of ''M'', the action of ''M'' is component-wise multiplication with a (fixed) element in 'C'
''n'' with positive entries.
|-
|valign="top"| '3.' || The
sesquilinear form
:
defines an
inner product on 'C'
''n''. (In fact, every inner product on 'C'
''n'' arises in this fashion from a Hermitian positive definite matrix.)
|-
|valign="top"| '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 ''x
i'' may optionally be restricted to fall in 'C'
''n''. In other words, ''M'' is of the form ''A
★ A'' where ''A'' is not necessarily square but must be injective in general.
|-
|valign="top"| '5.' || All the following matrices have a positive
determinant (the
Sylvester criterion):
★ the upper left 1-by-1 corner of
★ the upper left 2-by-2 corner of
★ the upper left 3-by-3 corner of
★ ...
★
itself
In other words, all of the leading
principal minors are positive.
For positive semidefinite matrices, all principal minors have to be non-negative.
The leading principal minors alone do not imply positive semidefiniteness, as can be seen from the example
:
|}
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."
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''.
Negative-definite, semidefinite and indefinite matrices
The ''n'' × ''n'' Hermitian matrix
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
, the matrix ''A''
★ ''A'' is positive semidefinite, and rank(
) = rank(''A''
★ ''A''). Reversely, any positive semidefinite matrix
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'.
Further properties
If
is positive semi-definite, one sometimes writes
and if
is positive-definite one writes
. This may be confusing, as sometimes
nonnegative matrices are also denoted in this way. Our notion comes from
functional analysis where positive definite matrices define
positive operators.
For positive semi-definite matrices
we write
if
, i.e.
is positive semi-definite. Equivalently for
.
{| cellspacing="0" cellpadding="2"
|-
|valign="top"| '1.' ||
Every positive definite matrix is
invertible and its inverse is also positive definite. If
then
.
|-
|valign="top"| '2.' || If
is positive definite and
is a real number, then
is positive definite.
If
and
are positive definite, then the sum
and the products
and
are also positive definite. If
, then
is also positive definite.
|-
|valign="top"| '3.' || If
then the diagonal entries
are real and positive. As a consequence
. Furthermore
:
.
|-
|valign="top"| '4.' || A matrix
is positive definite, if and only if there is a positive definite matrix
with
. One writes
. This matrix
is unique (but only under the assumption
). If
then
.
|-
|valign="top"| '5.' || If
then
. (Here
denotes
Kronecker product.)
|-
|valign="top"| '6.' || For matrices
write
for the entry-wise product of
and
, i.e. the matrix whose
entry is
. Then
is the Hadamard product of
and
. If
then
and if
are
real matrices, the following inequality, due to Oppenheim, holds:
|-
|valign="top"| '7.' || Let
and
hermitian. If
(
) then
(
).
|-
|valign="top"| '8.' || If
are real matrices then
.
|-
|valign="top"| '9.' || If
is real, then there is a
such that
where
is the
identity matrix.
|}
Non-Hermitian matrices
A real matrix ''M'' may have the property that ''x''
T''Mx'' > 0 for all nonzero real vectors ''x'' without being symmetric. The matrix
:
satisfies this property, because for all real vectors
such that
,
:
In general, we have ''x''
T''Mx'' > 0 for all real nonzero vectors ''x'' if and only if the symmetric part, (''M'' + ''M''
T) / 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.
See also
★
Square root of a matrix
★
Schur complement
★
Positive definite kernel
★
Positive-definite function
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).
★ Rajendra Bhatia. ''Positive definite matrices,''. Princeton Series in Applied Mathematics, 2007. ISBN 978-0691129181.