Атрибут Matlab для вещания Numpy?

Question

Я пытаюсь найти какой-то способ выставить вектор размером 3 из каждого столбца матрицы 3 * (большого числа) в Matlab. Конечно, я мог бы использовать цикл, но я пытаюсь найти более эффективное решение, немного похожее на широковещательное вещание. О, и я не могу использовать repmat, потому что у меня просто недостаточно памяти для его использования (поскольку он создает еще одну матрицу 3 * (большое число)) ...

Возможно ли это?

Answer 1

Петли не плохо в MATLAB больше благодаря оптимизации компилятора, например just-in-time acceleration (JITA). и т. д. В большинстве случаев я заметил, что решение с петлями в текущих версиях MATLAB: много быстрее, чем сложные (хотя и крутые: D) однострочные.

bsxfun может сделать трюк но в моем опыте, как правило, имеют проблемы с памятью, а также, но в меньшей степени, чем repmat.

Так синтаксис будет:

AA = bsxfun(@minus,A,b) где b вектор и A ваша большая матрица

Но я призываю вас профилировать хитрый вариант, а затем решить! Скорее всего, из-за ограничений памяти у вас может не быть выбора :)

Answer 2

Я не знаю, ускорит ли это код, но вычитание скаляра из вектора не имеет проблем с памятью. Поскольку размер вашей матрицы такой асимметричный, накладные расходы из цикла for на короткий размер незначительны.

Так может быть

matout = matin; 
for j = 1:size(matin, 1) %3 in this case 
    matout(j,:) = matin(j,:) - vec_to_subtract(j); 
end 

конечно, вы могли бы сделать это на месте, но я не знаю, если вы хотите, чтобы сохранить первоначальную матрицу.

Answer 3

На самом деле кажется, что http://www.frontiernet.net/~dmschwarz/genops.html (перегрузка оператора с помощью mex-файлов) тоже делает трюк, хотя я еще не тестировал его.

Answer 4

Другие ответы немного устарели - Matlab R2016b имеет added broadcasting as a standard feature. Пример из этого сообщения в блоге, который соответствует вопросу:

>> A = ones(2) + [1 5]' 
A = 
    2  2 
    6  6