如何在没有循环的MATLAB中计算矩阵所有列的外积的总和？

Question

我试图让标题说明一切。

我有一个矩阵，我需要计算列组合的外部产品的总和。我知道如何使用for循环。但是，有没有矢量化的方法呢？

testData = rand(x);   
answer = zeros(x);
for i=1:x
  for j=1:x
    answer = answer + testData(:,i)'*testData(:,j);
  end
end

Answer 1

答案在于这个问题Fast Outer Product

answer = testData.'*testData;

N.B。这给出了行的外部产品的总和。