Loopless affectation de sous-matrice dans Matlab

Question

I ont une matrice de taille FD-by-N et un vecteur de longueur AN d'entiers aléatoires dans la plage [1,a]. Je veux créer une matrice M de taille D * a telle que chaque colonne M(:,i) a le vecteur F(:,i) à partir de l'index (A(i)-1)*D+1 à (A(i)-1)*D+D.

Exemple:

F = [1 2 3 10 
    4 5 6 22] 

A = [3 2 1 2] 

a = 4 

M = [0 0 3 0 
    0 0 6 0 
    0 2 0 10 
    0 5 0 22 
    1 0 0 0 
    4 0 0 0 
    0 0 0 0 
    0 0 0 0] 

je peux le faire avec une simple boucle

 for i = 1 : N 
      M((A(i)-1)*D+1:(A(i)-1)*D+D,i) = F(:,i); 
     end 

mais pour les grandes N cela pourrait prendre un certain temps. Je cherche un moyen de le faire sans boucle.

Answer 1

Vous pouvez utiliser bsxfun pour une approche basée sur linear-indexing -

[D,N] = size(F); %// Get size of F 

start_idx = (A-1)*D+1 + [0:N-1]*D*a;   %// column start linear indices 
all_idx = bsxfun(@plus,start_idx,[0:D-1]'); %//'# all linear indices 

out = zeros(D*a,N); %// Initialize output array with zeros 
out(all_idx) = F; %// Insert values from F into output array 

run échantillon -

F = 
    1  2  3 10 
    4  5  6 22 
A = 
    3  2  1  2 
a = 
    4 
out = 
    0  0  3  0 
    0  0  6  0 
    0  2  0 10 
    0  5  0 22 
    1  0  0  0 
    4  0  0  0 
    0  0  0  0 
    0  0  0  0