1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
\name{is.indefinite}
\alias{is.indefinite}
\title{ Test matrix for positive indefiniteness }
\description{
This function returns TRUE if the argument, a square symmetric real matrix x, is indefinite.
That is, the matrix has both positive and negative eigenvalues.
}
\usage{
is.indefinite(x, tol=1e-8)
}
\arguments{
\item{x}{ a matrix }
\item{tol}{ a numeric tolerance level }
}
\details{
For an indefinite matrix, the matrix should positive and negative eigenvalues. The R function \code{eigen}
is used to compute the eigenvalues. If any of the eigenvalues is absolute value is less than the
given tolerance, that eigenvalue is replaced with zero. If the matrix has both positive and
negative eigenvalues, it is declared to be indefinite.
}
\value{
TRUE or FALSE.
}
\references{
Bellman, R. (1987). \emph{Matrix Analysis}, Second edition, Classics in Applied Mathematics,
Society for Industrial and Applied Mathematics.
}
\author{ Frederick Novomestky \email{fnovomes@poly.edu} }
\seealso{
\code{\link{is.positive.definite}},
\code{\link{is.positive.semi.definite}},
\code{\link{is.negative.definite}},
\code{\link{is.negative.semi.definite}}
}
\examples{
###
### identity matrix is always positive definite
###
I <- diag( 1, 3 )
is.indefinite( I )
###
### positive definite matrix
### eigenvalues are 3.4142136 2.0000000 0.585786
###
A <- matrix( c( 2, -1, 0, -1, 2, -1, 0, -1, 2 ), nrow=3, byrow=TRUE )
is.indefinite( A )
###
### positive semi-defnite matrix
### eigenvalues are 4.732051 1.267949 8.881784e-16
###
B <- matrix( c( 2, -1, 2, -1, 2, -1, 2, -1, 2 ), nrow=3, byrow=TRUE )
is.indefinite( B )
###
### negative definite matrix
### eigenvalues are -0.5857864 -2.0000000 -3.4142136
###
C <- matrix( c( -2, 1, 0, 1, -2, 1, 0, 1, -2 ), nrow=3, byrow=TRUE )
is.indefinite( C )
###
### negative semi-definite matrix
### eigenvalues are 1.894210e-16 -1.267949 -4.732051
###
D <- matrix( c( -2, 1, -2, 1, -2, 1, -2, 1, -2 ), nrow=3, byrow=TRUE )
is.indefinite( D )
###
### indefinite matrix
### eigenvalues are 3.828427 1.000000 -1.828427
###
E <- matrix( c( 1, 2, 0, 2, 1, 2, 0, 2, 1 ), nrow=3, byrow=TRUE )
is.indefinite( E )
}
\keyword{ math }
|
