matlab矩陣元素單元乘法的替代方法

Question

我正在進行大矩陣乘法。我有一個大矩陣A（至少5000x5000）和一個列向量V（5000x1）。在我的代碼中，每個V將按元素乘以A元素的每列。我做了一個循環

K = zeros(5000, 5000); 
for n=1:5000 
    K(:, n) = V.*A(:, n); 
end 

但它是如此之慢。所以，我創建指定爲V每列一個大矩陣，使得

MV=repmat(V,1,5000); 
K = MV.*A; 

速度快，但浪費了太多的記憶。當矩陣的大小增加時，它會使用太多的內存。是否有任何想法使用更少的內存，但速度很快？

Answer 1

經典bsxfun

K = bsxfun(@times, A, V); 

或者，你可能想看看James Tursa'sMTIMESX（FEX中找到）。