这个用矩阵相乘怎么表示？

清水小鱼

P为某矩阵
r（i），i=1···N为矩阵R的N个列向量

或者帮写个求和的matlab程序更好，感谢

sogooda

r用二维数组表示,r的第i列即为ri。在matlab中可以用r(:,i)表示ri，然后就可以作了，没啥难的啊。
s=0;
for i=1:N
s=s+r(:,i)'*P*r(:,i)*r(:,i)*r(:,i)';
end

清水小鱼

非常感谢楼上。谢谢:@) :@)