매트랩 1000 5 × 5 행렬의 공분산 계산 :I 이런 1000 5 × 5 행렬 (</strong><strong>내지 Xm)가
각각 $ (x_ij) 미터을 $로부터 인출 점 추정치 분포. 나는 각각의 $ x {ij} $의 공분산 cov
을 계산하고 싶습니다. 여기서 i = 1..n, j = 1..n은 빨간색 화살표 방향입니다.
예를 들어, $ X_m $의 분산은 'var (X, 0,3)'이며, 분산의 5x5 행렬을 제공합니다. 같은 방식으로 공분산을 계산할 수 있습니까? 대답에서
시도가
지금까지이 작업을 완료했습니다 당신이 (명령 창에서 edit cov
) cov
보면
for m=1:1000
Xm_new(m,:)=reshape(Xm(:,:,m)',25,1);
end
cov(Xm_new)
spy(Xm_new) gives me this unusual looking sparse matrix:
아 매우 좋습니다! 고맙습니다! 그러나 각 변수 x_ij의 1000 회의 관측 (3 차원에서 서로 겹쳐진 온 탑)에 기반한 5x5 공분산 행렬을 찾고 있습니다. 예를 들어 Xm (1,1, :)은 같은 변수에 대한 관찰 값입니다.이 함수가 무엇을하고 있습니까? – HCAI
5-by-5 입력에 대해 'cov' 함수를 1000 번 호출하고 5x5 출력을 저장하는 것과 같습니다 :'Xm_out (:, :, 1) = cov (Xm (:, :, 1)), Xm_out (:, :, 2) = cov (Xm (:, :, 2))', ...,'Xm_out (:, :, 999))','Xm_out (:, :, 1000) = cov (Xm (:, :, 1000))'로 시작한다. – horchler
아 아 ...각 x_ij에 대한 공분산을 찾으려고합니다. 즉, x12는 열 벡터 x1과 x2 사이의 공분산이고, 빨간색 화살표의 방향입니다. 올바른 질문을하지 않았습니까? – HCAI