Trasformare un (4D matrice * 1D vettore) il funzionamento in (a matrice 3D * 0D scalare) indipendente operazioni senza loop

Question

C'è un modo per vettorizzare il seguente:

for i = 1:6 
     te = k(:,:,:,i).*(c(i)); 
    end 

sto cercando di moltiplicare una matrice 4D , k, con un vettore, c, suddividendolo in operazioni indipendenti (a matrice 3D * scalare). Ho già altri due loop inevitabili all'interno di un ciclo while in questo file di funzione, e sto facendo del mio meglio per evitare i loop.

Qualsiasi commento su questo sarà molto apprezzato!

-SC

Answer 1

È possibile farlo utilizzando MTIMESX - uno strumento di moltiplicazione di matrici veloce con il supporto multidimensionale da James Tursa, si trovano in scambio di file di Matlab.

E 'semplice come:

C = mtimesx(A,B) 

esegue il calcolo C = A * B

Answer 2

A meno che non mi manca qualcosa, questo è un caso per bsxfun:

te=bsxfun(@times, k, permute(c,[3 4 1 2])); % c is a row vector 

O

te=bsxfun(@times, k, permute(c,[3 4 2 1])); % c is a column vector 

Si presume che la quarta dimensione di k abbia le stesse dimensioni di c. In caso contrario, è possibile utilizzare l'indicizzazione per sottotitoli:

te=bsxfun(@times, k(:,:,:,1:length(c)), permute(c,[3 4 2 1])); % c is a column vector