Supongamos que tengo dos matrices Nx2, Mx2 que representan N, M 2d vectores respectivamente. ¿Hay una manera simple y buena de calcular distancias entre cada par de vectores (n, m)?
La manera fácil pero ineficiente es, por supuesto:
d = zeros(N, M);
for i = 1:N,
for j = 1:M,
d(i,j) = norm(n(i,:) - m(j,:));
endfor;
endfor;
La respuesta más cercana que he encontrado es bsxfun
, usada así:
bsxfun(inline("x-y"),[1,2,3,4],[3;4;5;6])
ans =
-2 -1 0 1
-3 -2 -1 0
-4 -3 -2 -1
-5 -4 -3 -2
cartprod
, así que ahora puedo escribir: (1) x = cartprod(n(:,1), m(:,1));
(2) y = cartprod(n(:,2), m(:,2));
(3) d = sqrt((x(:,1)-x(:,2)).^2+(y(:,1)-y(:,2)).^2)
... ¡que funciona mucho más rápido!