複數矩陣乘法

Question

我想要做這個操作，但是如果你使用for循環，它對於一個大型矩陣來說確實很慢。 有人可以幫助做這個操作嗎？我們知道Y,A,M是矩陣，而mean是一列。

謝謝

Y=zeros(s,y); 
for i=1:y 
    Y(:,i)=A*(M(:,i)-mean); 
end 

Answer 1

Y = A * (M - repmat(mean, 1, y));