Discussion :

[thaomta]

bonjour,

j'ai un problème suivant:

je veux créer une matrice des coordonnées à partir des arrays:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
x = 1:10;
y = 20:50;
z = 100:500;

je peux faire une boucle comme ça:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
matrice = [];
for i = 1:10
   for j = 20:50
      for k = 100:500
            point = [i j k];
            matrice = [matrice; point];
      end
   end
end

mais ce n'est pas optimal pour le calcul car le temps de calcul pour une boucle for...end dans matlab est très long.

Ma question: est-ce qu'il y a des fonctions de Matlab ou des autres manières qui peuvent créer la matrice ci-dessus sans utiliser la boucle for...end ?

je vous remercie pour votre réponse!!!!!

[Jerome Briot]

Comme tu connais les dimensions de la matrice au final, tu peux simplement préallouer la mémoire de la matrice pour accélérer la boucle for-end. Voir la FAQ : Qu'est-ce que la préallocation de mémoire ?

Ensuite pour une solution vectorisée, il faut par exemple utiliser la fonction repmat. Mais cela demande de maîtriser l'ordre de stockage des éléments des matrices en mémoire par MATLAB. Lire ce tutoriel : Introduction à la gestion des matrices sous MATLAB

Bref voici un comparatif des trois solutions :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
x = 1:10;
y = 20:50;
z = 100:500;
 
nx = numel(x);
ny = numel(y);
nz = numel(z);
 
tic
matrice1 = [];
for i = x
    for j = y
        for k = z
            point = [i j k];
            matrice1 = [matrice1; point];
        end
    end
end
toc
 
tic
matrice2 = zeros(nx*ny*nz,3);
u = 1;
for i = x
    for j = y
        for k = z
            matrice2(u,:) = [i j k];
            u = u+1;
        end
    end
end
 
toc
 
tic
 
xx = repmat(x,ny*nz,1);
yy = repmat(y,nz,nx);
zz = repmat(z,1,nx*ny);
matrice3 = [xx(:) yy(:) zz(:)];
 
toc
 
all(matrice1(:)==matrice2(:)) && all(matrice1(:)==matrice3(:))

ce qui donne sur ma machine :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
Elapsed time is 71.831909 seconds.
Elapsed time is 0.259971 seconds.
Elapsed time is 0.010429 seconds.
 
ans =
 
     1

[thaomta]

merci à Dut!


j'ai trouvé l'autre manière pour la créer:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
x = 1:10;
y = 20:50;
z = 100:500;
[x1, y1, z1] = meshgrid(x,y,z);
x2 = reshape(x1,length(x)*length(y)*length(z),1);
y2 = reshape(y1,length(x)*length(y)*length(z),1);
z2 = reshape(z1,length(x)*length(y)*length(z),1);
matrice4=horzcat(x2,y2,z2);

le temps de calcul pour la matrice3 (de Dut) et la matrice4 sur mon ordinateur:

Elapsed time is 0.006817 seconds. (matrice3)
Elapsed time is 0.004366 seconds. (matrice4)

[Jerome Briot]

Sauf que les valeurs dans matrice4 n'apparaissent plus dans l'ordre que tu as spécifié dans ton premier message