Matlab有效的稀疏矩陣乘法

Question

我有一個稀疏矩陣，它只有三個對角線中的元素。例如。

我也有在那裏我希望每個元素由對應的元件列向量的每一行中乘以稀疏矩陣的每行中的列向量。有沒有在MATLAB中做到這一點的有效方法？如果稀疏矩陣稱爲A和列向量B，我只試過

A.*repmat(B,[1,9]) 

這顯然是低效的。

Answer 1

這裏有一種方法：

C = bsxfun(@times, A, B) 

據docs，結果矩陣C是稀疏的：

二元運算產生稀疏的結果，如果兩個操作數都很少，如果兩個全部結果都滿。 對於混合操作數，除非操作保留稀疏度，否則結果將爲滿。如果S稀疏並且F滿，則S + F，S * F和F \ S滿，，而S. * F和S & F稀疏。在某些情況下，即使矩陣只有零個元素，結果也可能很稀疏。