Matlab equivalente a la transmisión de Numpy?

Question

Estoy tratando de encontrar algún modo de sustraer un vector de tamaño 3 de cada columna de una matriz de 3 * (un gran número) en Matlab. Por supuesto que podría usar un bucle, pero estoy tratando de encontrar una solución más eficiente, un poco como la transmisión numpy. Ah, y no puedo usar repmat porque simplemente no tengo suficiente memoria para usarlo (ya que crea otra matriz 3 * (un número grande)) ...

¿Esto es posible?

Answer 1

Los bucles ya no están mal en MATLAB gracias a las optimizaciones del compilador como just-in-time acceleration (JITA). etc. La mayoría de las veces, he notado que una solución con bucles en las versiones actuales de MATLAB es mucho más rápido que complicado (aunque, genial: D) one-liners.

bsxfun podría hacer el truco, pero en mi experiencia, que tiende a tener problemas de memoria, así pero menos que repmat.

Así que la sintaxis sería:

AA = bsxfun(@minus,A,b) donde b es el vector y A es su gran matriz

pero insto a crear perfiles para la versión descabellada y luego decide! Probablemente, debido a restricciones de memoria, es posible que no tenga una opción :)

Answer 2

No sé si esto acelerará el código, pero la sustracción de un escalar de un vector no tiene problemas de memoria. Dado que el tamaño de su matriz es tan asimétrico, la sobrecarga de un bucle for en la dimensión corta es insignificante.

Así que tal vez

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

por supuesto, usted puede hacer esto en su lugar, pero yo no sabía si quería conservar la matriz original.

Answer 3

En realidad, parece que http://www.frontiernet.net/~dmschwarz/genops.html (sobrecarga del operador con archivos mex) hace el truco también, a pesar de que aún no lo he probado.

Answer 4

Las otras respuestas están un poco desactualizadas - Matlab R2016b parece tener added broadcasting as a standard feature. Un ejemplo de esa publicación de blog que coincide con la pregunta:

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