2014-04-04 2 views
-4

제목에서 언급 한 바와 같이. 매트릭스의 고유 값의 곱을 계산하도록 가르쳐 줄 수 있습니까? R 함수 또는 패키지를 사용하는 것이 좋습니다. 감사.행렬의 고유 값 계산 방법 R

+1

아마이 검색 결과 : ['[R] 고유 matrix' (http://stackoverflow.com/search?q= [R] + 고유치 + 매트릭스) – Jaap

+1

어떤 질문을 도움이 될 것입니다 일반적으로 당신이 시도한 것에 대한 예가 적어도 있어야합니다. –

+1

스택 오버플로에 오신 것을 환영합니다. [좋은 질문] (http://stackoverflow.com/help/how-to-ask) 및 [재현 가능한 예] (http://stackoverflow.com/questions/5963269)를 제공하는 방법에 대해 스스로에게 알려주십시오./만드는 방법 - 훌륭한 - r - reproducible - 예) – Jaap

답변

1

이 문제를 해결하기 위해 얼마나 노력 했습니까?
행렬의 고유 값의 곱을 계산하는 것은 매우 기본입니다.
설명서 파일 ?matrix, ?eigen?prod에 대한 설명을 남겨 두겠습니다.

> m <- matrix(rnorm(25, 5, 10), 5, 5) 
> m 
      [,1]  [,2]  [,3]  [,4]  [,5] 
[1,] -6.068987 5.657092 16.106767 11.701708 6.060705 
[2,] 2.574272 7.420419 9.871342 6.049024 -5.930012 
[3,] -10.549236 23.885948 -7.187153 9.084914 -4.843331 
[4,] 17.815534 3.348587 1.925881 -3.900034 10.506535 
[5,] -10.143544 13.401505 11.022406 5.437238 17.646310 

> eigen(m)$values 
[1] 20.720634+0.000000i -15.994616+9.683080i -15.994616-9.683080i 
[4] 9.589576+6.645961i 9.589576-6.645961i 

> prod(eigen(m)$values) 
[1] 986078.9-0i 
0

??eigenvalueeigen()입니다. ?eigen을 통해 도움말 페이지를 읽는 것은 고유 벡터 (eigenvectors) 대신 고유 값을 추출하는 방법을 알려줍니다. 충분히 흥미롭게도 ??product이 아니며prod()을 가리 킵니다. 이제이 기능에 대해 알게되었습니다. 그래서 :

> foo <- matrix(runif(9),3,3) 
> prod(eigen(foo)$values) 
[1] -0.07673157