Как векторизовать векторно-матричную поэтапную операцию?

У меня есть следующий код:Как векторизовать векторно-матричную поэтапную операцию?

A = rand(N1,N2); 
b = rand(1,N1); 
B = zeros(N1,N2); 

for i=1:N1 
    for j=1:N2 
     B(i,j) = A(i,j)*b(i); 
    end 
end

Вопрос заключается в том, чтобы записать его в векторной форме операции? Что-то вроде B(:,:) = A(:,:).*b(:).

источник

2013-06-06 Alex

Поскольку все переменные являются нулями, в чем смысл? просто 'B = нули (N1, N2);' достаточно – bla

Ваш внутренний цикл тривиально удален. Векторизация внешнего цикла - это то, что поддерживает Octave, но MATLAB этого не делает, хотя, как указано ниже, 'bsxfun' может быть быстрее, чем один явный цикл, и, конечно же, должен быть быстрее, чем вложенные явные циклы, которые вы указали выше. – user57368

@ natan Я почти уверен, что OP поставил их там просто, чтобы проиллюстрировать размеры. У тебя наверняка есть какая-то вера в человека! – paddy

Простое использование для bsxfun:

B = bsxfun(@times, A, b')

источник

2013-06-06 04:19:26 paddy

, увеличивает ли производительность? – Alex

Да, 'bsxfun' поразит ваши мелкие вложенные петли. знак равно – paddy

Вы также можете попробовать:

B = A*.(repmat(b,N2,1))';

Здесь, во-первых, вы производите N2 повторную версию вектора b и умножить его A в поэлементного образом

источник

2013-06-06 04:57:41

Как векторизовать векторно-матричную поэтапную операцию?

ответ

Смежные вопросы