2014-09-04 3 views
-1

R을 사용하는 방법을 아는 사람이라면 chr 요소로 구성된 일반 매트릭스의 고유 값의 특성 다항식을 계산하고 표시 할 수 있습니까? 는 일반 매트릭스의 고유 값 찾기 R

m <- matrix(c('a','b','c','d','e','f','g','h','i','l','m','n'),4,4) 

내가

가 사전에 감사합니다 매우 큰 행렬에이 방법을 적용 할 필요가 있음을 고려하시기 바랍니다 즉 말을

+0

chr 요소로 구성된 제네릭 행렬의 고유 값의 특성 다항식이 정확히 무엇인지에 대한 정보를 포함하면 도움이 될 수 있습니다. 여기에있는 대부분의 사람들은 프로그래머이고 문자 공간에서 비공식적 인 수학에 익숙하지 않다 ... –

+1

나는 생각한다 * (??) OP는 기호 계산을 찾고 있을지도 모른다. (이것은 R에서 불가능하다. –

+0

제안 해 주셔서 감사합니다. 그러나 간단한 정의가 행렬 대수에 익숙하지 않은 사람 (또는 문자 공간에서 일반 수학 방법을 일반적인 수학 방법이라고 생각하는 경우)을 도움이되지 않을 수 있다고 생각합니다. 물론 물론 줄 수는 없습니다. 여기에이 주제에 대한 수업이 있습니다. 어쩌면 전문 사용자 만이 대답 할 수 있습니다. – Shatz

답변

1

당신은 Ryacas 패키지와 함께이 작업을 수행 할 수 있지만, 당신은 할 수 있습니다 먼저 시스템에 Yacas를 설치하기 위해 필요한 농구를 뛰어 넘으십시오.

library("Ryacas") 
m <- matrix(letters[1:16],4,byrow=TRUE) 
yrow <- function(x) paste0("{",paste(x,collapse=","),"}") 
yrow(m[1,]) ## "{a,b,c,d}" 
ymat <- function(x) yrow(apply(x,1,yrow)) 
cheqstr <- function(x) { 
    paste0("Expand(CharacteristicEquation(", 
      ymat(x),",x),x)") 
} 
yacas(cheqstr(m)) 
## (a-x)*(f-x)*(k-x)*(p-x)-(a-x)*(f-x)*l*o+(a-x)*h*j*o-d*e*j*o- 
## (a-x)*g*j*(p-x)+(a-x)*g*l*n-(a-x)*h*(k-x)*n+d*e*(k-x)*n+c*e*j*(p-x)- 
## c*e*l*n+c*h*i*n-d*g*i*n-b*e*(k-x)*(p-x)+b*e*l*o-b*h*i*o+d*(f-x)*i*o+ 
## b*g*i*(p-x)-b*g*l*m+b*h*(k-x)*m-d*(f-x)*(k-x)*m-c*(f-x)*i*(p-x)+ 
## c*(f-x)*l*m-c*h*j*m+d*g*j*m;