유클리드 거리 계산을 위해 행렬 더하기 및 곱셈을 조작하는 방법은 무엇입니까?

Question

난 다음 가지고 D는 N1의 엑스 D1 행렬

A = [1 2 3; 4 5 6; 7 8 9]; 
B = [10 11 12; 13 14 15]; 

[N1, D1] = size(A); 
[N2, D2] = size(B); 

A_sq = sum(A.^2, 2); 
B_sq = sum(B.^2, 2)'; 
D = A_sq(:,ones(1,N2)) + B_sq(ones(1,N1),:) - 2.*(A*B'); 

.

나는 하나의 단계 D에 대한 표현을 쓰고 싶은, 즉,이 같은 (이 그림의 목적이지만, 위의 코드와 같은 유클리드 거리 계산한다) :

D = sum(A - B).^2; 

I을 어떤 조언을 주셔서 감사합니다.

Answer 1

당신은 통계 도구 상자를 가지고 있다면 그냥하지 pdist2, 사용할 수 있습니다

D = pdist2(A,B).^2 

아니면 bsxfun 및 permute 사용하여 수동으로 작업을 수행 할 수 있습니다 귀하의 예를 들어

D = permute((sum(bsxfun(@minus, A, permute(B, [3 2 1])).^2,2)), [1 3 2]); 

이

행렬을

A = [1 2 3; 4 5 6; 7 8 9]; 
B = [10 11 12; 13 14 15]; 

상기 g ives

D = 
    243 432 
    108 243 
    27 108