Matlab/Octaveを使用して2つの行列の対応する列の外積の合計をベクトル化する

Question

2つの行列AとBが次のように列ベクトルで構成されているとします。

A = [a_1,a_2,...,a_N]; 
B = [b_1,b_2,...,b_N]; 

ここでBに対応する列にA内のすべての列の外側の積和の計算をベクトル化するためにどのような方法がありますが、私の非ベクトル化ソリューションです。

S = zeros(size(A,1), size(B,1)); 
for n=1:N 
    S = S + A(:,n)*B(:,n)'; % S = S + a_n * b_n' 
end 

ご協力いただければ幸いです。

Answer 1

あなたはNが何であるかについては明らかではないですが、私が想定しているN =列ベクトルの数 - あなたは、単に*のB」

A = rand(3,4); 
B = rand(3,4); 
N = size(A,2); 
S = zeros(size(A,1), size(B,1)); 
for n=1:N 
    S = S + A(:,n)*B(:,n)'; % S = S + a_n * b_n' 
end 
%Check that you are doing A*B' 
S == A*B' 
>> ans = 

1  1  1 
1  1  1 
1  1  1