2014-05-25 2 views
1

[0,2]에 고유 값을 갖는 작고 잘 조절 된 헤르 미안 행렬 L이 있습니다. 행렬의 최소 고유치의 동일 역으로 역의 두 번째 규범한다고 때문에 이상하다왜 내가 matlab에 잘못된 매트릭스 표준을 얻고 있습니까?

>> norm(inv(L)) 

ans = 

    2.0788 

>> min(eig(L)) 

ans = 

    0.5000 

: L의 역의 규범을 계산하는 동안 나는 이상한 결과를 얻고있다.

기계 연산에 의해 발생하는 오류에 대해 알고 있지만, 작고, 허름하고 잘 설명 된 예제에서는 무시할 수있을 것으로 예상했습니다.

여기 리눅스 민트 16 (페트라)의 매트릭스 https://www.dropbox.com/s/nh1wegrnn53wb6w/matrix.mat

I는 MATLAB 8.2.0.701 (R2013b)를 사용하고있다.

+1

729x729는 "작은"매트릭스입니까? – Daniel

+0

@ 대니얼, 네, 그렇습니다. 그리고이 하나는 드문 드문 한 매트릭스입니다. – Moonwalker

+0

당신은'[0,1]에 고유치가있는 잘 조절 된 헤르 미안 행렬 L을 말합니다. 그러나 이것은 제가 얻은 것이 아닙니다 :'e = eig (L); max (e) -> 1.3789, min (e) -> 0.5000' – Nasser

답변

5

매트릭스가 잘 조절되었음을 지적 했으므로 수치적인 문제는 아닙니다. 역의

초 표준 (즉, 양의 확실한) 행렬

이 행렬은 포지티브 고유과 에르 미트 경우에만 참의 최소 고유치의 동일한 역되어야한다. 위키 피 디아에서 : 행렬 A의 스펙트럼 표준은 A의 가장 큰 특이 값, 즉 양의 semidefinite 행렬의 가장 큰 고유 값의 제곱근입니다. A * A

여기에서 역함수를 다음과 같이 계산합니다.

[v,d] = eig(L'*L); 
1.0/sqrt(min(diag(d))) = 2.0788539 
norm(inv(L)) = 2.0788539 

예상대로.

+3

그러면 매트릭스가 실제로 헤미메이트가 될 수 없다는 것을 보여줍니다. 즉, 실제 엔트리가 있으면 매트릭스입니다. 아니면 확실하지 않습니다. 대칭 양의 확정 (spd) 행렬의 경우 원래의 추론이 정확하므로 고유 값 또한 고유 값입니다. – LutzL

+0

@ LutzL, 예, 행렬은 실제로 헤르메이트가 아닙니다. 도와 주셔서 감사합니다! – Moonwalker

+0

LutzL에게 감사드립니다. 귀하의 의견을 반영하도록 작은 편집을했습니다. – y300

관련 문제