2014-10-11 267 views

回答

1
> B <- matrix(1:4, 2) 
> eig <- eigen(B) 
> eig$vectors[,which.max(eig$values)] 
[1] -0.5657675 -0.8245648 

或像这样

> n <- 3 
> A <- matrix(round(runif(n*n),2),nrow=n) 
> A 
    [,1] [,2] [,3] 
[1,] 0.54 0.90 0.82 
[2,] 0.09 0.42 0.95 
[3,] 0.17 0.75 0.69 
> evv.A <- eigen(A) 
> evv.A 
$values 
[1] 1.6230202 0.3095823 -0.2826025 

$vectors 
      [,1]  [,2]  [,3] 
[1,] -0.7455955 -0.97996498 0.3069212 
[2,] -0.4464598 0.18600090 -0.7788676 
[3,] -0.4947332 0.07122006 0.5469594 

> 
> k <- which(abs(evv.A$values)==max(abs(evv.A$values))) 
> evv.A$vectors[,k] 
[1] -0.7455955 -0.4464598 -0.4947332 
> k 
[1] 1 
+0

感谢您的答复。你的第二个建议将取值中第一大的特征值并显示其特征向量。根据你的例子,我需要选择两个特征值1.6230202 0.3095823,并且必须显示它的特征向量。任何方法做到这一点? – 2014-10-11 04:58:11

+0

做这个'k < - 1:2',然后'evv.A $ vectors [,k]'语句应该给出你想要的矩阵。用'?eigen'阅读特征的帮助。 – Bhas 2014-10-11 05:32:08

+0

如果有帮助,请接受解决方案作为答案 – RUser 2014-10-11 05:58:28

相关问题