这是我的问题的简化版本。假设我有以下矩阵:
[0 4x
.6 0]
我想求解x,因为我希望我的第一个(最大)特征值等于1。有没有一种简单的方法可以在R中解决这个问题?
尝试
f <- function(x){
m <- matrix(c(0, 4*x, .6, 0), 2, 2, byrow=TRUE)
Re(eigen(m)$values[1]) - 1
}
uniroot(f, c(-1, 1))$root