R语言检查矩阵是否正定
省流摘要
install.packages("matrixcalc")
library(matrixcalc)
A <- diag( 1, 3 ) #创建一个对角阵
is.positive.definite( A ) #返回TRUE证明该矩阵是正定矩阵
B <- matrix( c( 2, -1, 2, -1, 2, -1, 2, -1, 2 ), nrow=3, byrow=TRUE )
is.positive.definite( B ) #返回FALSE证明该矩阵不是正定矩阵
#参考链接:https://www.rdocumentation.org/packages/matrixcalc/versions/1.0-6/topics/is.positive.definite