Discussion :

[med.doc]

Bonjour à tous,

je cherche de l'aide pour implémenter un algo qui prend en entrée une matrice de dimension L*C et qui pour chaque colonne contenant des valeurs me retourne une matrice de à priori autant de colonnes C mais en nombre de lignes toutes les combinaisons possibles.

Par exemple, si en entrée :

Code :

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

1
2
3
4
5
M =
 
   11   21   31
   10   20   30
    9   19   29

en sortie :

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
V =
 
   11   21   31
   10   21   31
    9   21   31
   11   20   31
   10   20   31
    9   20   31
   11   19   31
   10   19   31
    9   19   31
   11   21   30
   10   21   30
    9   21   30
   11   20   30
   10   20   30
    9   20   30
   11   19   30
   10   19   30
    9   19   30
   11   21   29
   10   21   29
    9   21   29
   11   20   29
   10   20   29
    9   20   29
   11   19   29
   10   19   29
    9   19   29

J'ai essayer avec la fonction

Code :

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

ndgrid

le problème est qu'elle prend en paramètre que des vecteurs ... et son utilisation me paraît un peu restrictive.
Je vous remercie d'avance pour votre aide.

[Jerome Briot]

Une solution :

Code :

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

1
2
3
4
5
6
7
8
9
M = [11   21   31
    10   20   30
    9   19   29];
 
A = repmat(M(:,1), 9, 1);
B = kron(M(:,2), ones(3,3));
C = repmat(M(:,3), 1, 9).';
 
V = [A(:) B(:) C(:)];

[med.doc]

merci beaucoup pour votre réponse. Cela marche parfaitement sauf que là vous supposez connaître les dimensions de M(3x3), mais si M est de dimension différente, par exemple 13x5 ? Je souhaite le faire dans le cas général, à savoir un algorithme qui ne dépendrai que de nL (nombre de lignes de ma matrice d'entrée) et nC (nombre de colonnes).

Merci encore pour votre aide

[Jerome Briot]

Je suis loin d'être un spécialiste, il y a peut être une solution plus simple, plus efficace et surtout plus robuste :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
nL = size(M,1);
nC = size(M,2);
 
V = zeros(nL^nC,nC);
 
for n = 0:nC-1
 
    t = kron(M(:,n+1), ones(nL^n,1));
 
    t = repmat(t,1, nL^(nC-n-1));
 
    V(:,n+1) = t(:);
 
end

A tester…

[med.doc]

je m'incline devant votre génie !
merci beaucoup pour votre aide.