모든 쌍의 행렬을 계산하십시오.

Question

동일한 행 수의 2 행렬을 취하여 스칼라 값을 구하는 함수가 f입니다. 이제는 두 개의 행렬 목록을 취하고 모든 쌍에 f을 호출하는 새 함수를 만들 가능성을 찾고 있습니다.

나는이 루프의 무어 효율적인 구현이 필요합니다

% X = cell of matrices 
% Y = cell of matrices 
for k=1:length(X) 
    for l=1:length(Y) 
     M(k,l) = f(X{k},Y{l}); 
    end 
end 

제곱 거리

을의 평균 수 예를 f를 들어

(그것은 X와 Y는 세포 인 요구 사항은 아닙니다)

f = @(X,Y) mean(mean(bsxfun(@plus,dot(X,X,1)',dot(Y,Y,1))-2*(X'*Y))); 

그러나 질문이 없습니다 f, 그것은 훨씬 더 복잡한 문제의 예일뿐입니다.

Answer 1

먼저 M을 사전 할당해야합니다. 그럼 당신은 X와 Y의 모든 요소 쌍의 목록 생성 meshgrid를 사용하여 아주 쉽게 하나 개의 루프를 잘라 수 있습니다 : 당신이 그렇게 제대로 벡터화되어 f을 구현할 수있는 경우에, 그러나

[K, L] = meshgrid(1:length(X), 1:length(Y)); 

M = zeros(size(K)); %//This preallocation alone should give you a significant speed up 
for ii = 1:numel(K) 
    M(ii) = f(X{K(ii)},Y{L(ii)}); 
end 

을, 당신은 수 있습니다 두 개의 3D 행렬, 즉 X와 Y 쌍의 전체 목록을 전달할 수 있고 전혀 루프없이 수행 할 수 있습니다. 그러나 이것은 당신의 f이하는 일에 전적으로 의존합니다. 내가 이것을 할 수있는 더 좋은 방법이 생각하지만 cellfun에게

fy = @(x,Y) cellfun(@(y) f(x,y), Y); 
fx = @(X,Y) cellfun(@(x) fy(x,Y), X, 'UniformOutput', false); 
newF = @(X,Y) cell2mat(fx(X,Y)')'; 

을 사용하고 작동 가능한 솔루션에

Answer 2

.