小编Ign*_*cio的帖子

我想知道是否有任何方法可以更有效地执行以下矩阵乘法,例如向量化它:

m=1000;n=500;
a=zeros(n,1);
b=rand(n,1);
A=rand(m,n);
B=rand(m,m);
for i=1:n
    a(i)=b'*(A'*B(i,:)'*B(i,:)*A)*b;
end

提前致谢