2017-12-11 10 views
0

나는 줄리아에서 임의의 스파 스 대칭 긍정적 인 명확한 매트릭스를 만들려고합니다. 내가 시도줄리아의 스파 스 대칭 긍정적 인 한정 매트릭스

한 가지 방법이었다

n = 10 
density = 0.1 
X = sprand(n,n,density) 
X = 0.5*(X+X') + n*eye(n) 

@test isposdef(X) == true 

이 작동하지만 그것은 당신에게 내가 원하는하지 않는 대각선 지배 매트릭스를 제공합니다. 항목은 모두 동일한 임의의 분포에서 나온 것이어야합니다.

+0

샘플링 할 정확한 무작위 분포를 지정해야한다고 생각합니다 (sym pos def 행렬 세트는 압축되지 않으므로 균일하지 않을 수 있습니다). –

+0

관련 메모에서 특정 언어로 어떻게 할 수 있는지 걱정하기 전에 수학적으로 정확하게 원하는 것을 지정할 수 있어야합니다. –

답변

0

벡터의 바깥 쪽 합계에 더하여 하나의 ID (예 : 원할 경우 작음)를 복수로 사용할 수 있습니다 (예 : M = eps*I + X*X' + Y*Y' + .). 바깥 쪽 제품이 양의 값을 가지므로 eps> 0이면 양의 값이됩니다.

모든 요소가 동일한 분포를 갖는 명확한 매트릭스를 가질 수 있습니까? 예를 들어 M [0,0] * M [1,1]> M [0,1] * M [1,0]과 같이 보장 할 수 없기 때문에 나는 생각하지 못했습니다. 긍정적 인 명확한 것.