How to do a quick division between a column array and a matrix?

1 次查看（过去 30 天）

Benson Gou 2019-1-28

评论： Benson Gou 2019-1-28

Hello, All,

I want to divid a column array b (N x 1) by each column in a big sparse matrix A (N x M). My code is as follows:

dividCol = zeros(N,M);

for i = 1 : M

nonCol = find(abs(A(:,i))>0);

dividCol(nonCol,i) = b(nonCol)./A(nonCol,i);

end

But it is very slow. I do know if there exist a faster approach to write this code.

Thanks a lot in advance.

Benson

Omer Yasin Birey 2019-1-28

编辑：Omer Yasin Birey 2019-1-28

Hi Benson, try this

% b = randi(25,25,1);
% A = randi(25,25,5);
[row,col] = find(abs(A(:,:))>0);
dividCol = bsxfun(@rdivide,b(row,1),A(row,col));

Omer Yasin Birey 2019-1-28

Yes I checked the code and it does every operation twice. I had to use loop and updated the code. Try it instead.

% b = randi(25,7,1);
% A = randi(25,7,14);
for i = 1:M
nonCol = find(abs(A(:,i))>0);
dividCol = bsxfun(@rdivide,b(nonCol),A(nonCol,i));
end

Benson Gou 2019-1-28

Omer,

Thanks for your prompt reply.

Is it possible to do this without using a loop? because loop normally takes more time.

Benson

在 Help Center 和 File Exchange 中查找有关 Denoising and Compression 的更多信息

Find the treasures in MATLAB Central and discover how the community can help you!

Translated by