微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

R中的特征值分解

如何解决R中的特征值分解

请原谅这个简单的问题,但如何在 R 中进行特征值分解?公式为 A=VDV^(-1) 其中 A 是方阵,V 是包含 A 的特征向量的矩阵,D 是包含 A 的不同特征值的对角矩阵。感谢您的帮助。

下面是一个可重现的例子:

##create the matrix
matrixa <- cbind(c(0.589,0.202),c(0.033,0.869))
##This is what I tried but it doesn't seem right
(eigen(matrixa)[[2]])*(eigen(matrixa)$values)*(solve(eigen(matrixa)[[2]]))

解决方法

你可以试试%*% + diag

with(
  eigen(matrixa),vectors %*% diag(values) %*% solve(vectors)
)

给出

      [,1]  [,2]
[1,] 0.589 0.033
[2,] 0.202 0.869

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。