Discussion :

[romwn]

Bonjour à tous,

Je souhaite créer une matrice de telle sorte que je puisse avoir :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
A =
1 1 1
1 2 2
1 3 3
2 1 4
2 2 5
2 3 6
3 1 7
3 2 8
3 3 9

La première colonne représentant les X, la deuxième les Y, la troisième représentant les Z.
J'ai pour cela créé 3 vecteurs :

Code :

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

1
2
3
4
 
x = [1 2 3];
y = [1 2 3];
z = [1 2 3 4 5 6 7 8 9];

J'ai alors tenté le code suivant pour construire ma matrice :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
count = 1;
for i=1:length(x)
for j=1:length(y)
XYZ(i,1)=x(i);
XYZ(i,2)=y(j);
XYZ(i,3)=z(count);
count = count+1;
end
end

J'obtiens alors 9 tableaux en 9x3, donc aucun ne me donne le résultat cherché...

J'ai essayé le même programme en C, et tout se déroule bien, je découvre MATLAB, et ne voit donc pas où se situe mon erreur...

Vous me seriez d'une grande aide !

Romwn

[Gooby]

Bonjour,

Peux tu nous montrer le résultat de ce script en écrivant XYZ dans ta command window?

[romwn]

Après avoir redémarré Matlab j'ai réinitialisé les valeurs à nouveau, je n'obtiens plus le même résultat.
Désormais j'ai bien un tableau en 2D, mais toujours pas le résultat obtenu :

Code :

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

1
2
3
4
5
6
 
XYZ =
 
1	1	3
2	2	6
3	3	9

J'ai cru comprendre que les boucles for n'étaient pas ce qu'il y a de mieux sur Matlab, peut être que la double boucle for est le problème ..?

[Gooby]

Ta variable d'itération varie de 1 à 3 dans ton code, pas étonnant que ta matrice ne comporte que 3 lignes si tu ne trouves pas comment faire, je t'aiderai tout à l'heure

[romwn]

J'ai bien réessayé de faire tourner le même code en C et Java, ça marche bien...
Je vois bien que je fais itérer mon i jusqu'à 3, mais pour chaque valeur de i, je fais itérer j de 1 à 3, donc logiquement je devrais avoir

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
i = 1
   j = 1
   j = 2
   j = 3
i = 2
   j = 1
   j = 2
   j = 3
i = 3
   j = 1
   j = 2
   j = 3

Non ? C'est en tout cas le but recherché pourtant

[Gooby]

Certes, mais ici

Code :

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

1
2
3
4
5
6
7
8
9
count = 1;
for i=1:length(x)
for j=1:length(y)
XYZ(i,1)=x(i);
XYZ(i,2)=y(j);
XYZ(i,3)=z(count);
count = count+1;
end
end

Tu attribues des valeurs aux lignes i=1, i=2 et i=3 de ta matrice. A aucun moment tu n'attribues des valeurs aux lignes suivantes.

Code :

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

1
2
3
4
5
6
7
8
9
count = 1;
for i=1:length(x)
for j=1:length(y)
XYZ(count,1)=x(i);
XYZ(count,2)=y(j);
XYZ(count,3)=z(count);
count = count+1;
end
end

et comme ça? à noter que tu aurais pu utiliser une troisième boucle.

Il y a également d'autres méthodes sans passer par des boucles, en utilisant des concaténations à bon escient par exemple:

Code :

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

1
2
3
 
XYZ= cat(2,[x'; x'; x'], [y'; y'; y'], z);
XYZ(:,2) = sortrows(XYZ(:,2));

[dionisis]

Salut romwn,

Je ne suis pas 100% sur de ce que tu veux faire donc desole si c est a cote de la plaque. J aurais fais quelque chose comme ca:

Code :

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

1
2
3
4
5
n = 3
y = repmat(1:n,1,n)'
x = sort(y)
z = (1:n*n)'
XYZ = [x y z]

[romwn]

Salut gooby, j'ai essayé le premier code, et j'ai obtenu le résultat suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
XYZ = 
1	1	1
2	2	2
3	3	3
1	1	4
2	2	5
3	3	6
1	1	7
2	2	8
3	3	9

J'avais essayé de passer par les concaténations, je vais essayer d'y passer un peu plus de temps..!

dionisis mon sauveur ! C'est (presque) ça ! Au final c'est le résultat recherché, mais la méthode ne peut s'appliquer à mon application

En l'occurence pour mon application, mes X correspondent à un tableau qui va de -200 à +200 avec un pas qui peut varier, mes Y vont de -500 à +4000. Mais les valeurs peuvent changer, il faut donc que je passe par les éléments de X, Y et Z directement et non par un n, de plus on ne peut créer x à partir de y ...

Je vais essayer de passer quelques heures sur la doc, mais j'ai peur qu'on ne puisse créer la table de cette façon ..

[Paul-Ceria]

Le code pour construire A prend une ligne :

Code :

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

1
2
3
4
5
6
 
x=[1 2 3]
y=[1 2 3]
z=[1 2 3 4 5 6 7 8 9]
 
A=[kron(x,ones(1,3))' repmat(y,1,3)' z']

[romwn]

Bonjour Paul-Ceria,

J'ai essayé le code et en effet ça marche, seulement comme je l'expliquait précédemment, mes valeurs de X et Y sont respectivement stockées dans une matrice 1*627 double et 1*1401 double...
mes Z sont normalement stockées dans une matrice 1401*627 double, que j'ai ensuite transformé pour avoir une matrice 1*878427 double en utilisant le code suivant :

Code :

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

z = z(:)';

en essayant

Code :

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

A=[kron(x,ones(1,3))' repmat(y,1,3)' z'];

j'ai obtenu l'erreur suivante :

Code :

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

1
2
Error using horzcat
Dimensions of matrices being concatenated are not consistent.

Seulement à quel moment utilise-t-on horzcat dans l'instruction ..?
L'adaptation vecteur à matrice est il possible ? J'ai essayé de trouver une solution pour transformer une matrice 1*nbElements à un tableau, mais sans succès...

[Paul-Ceria]

Tu nous montres un exemple où z était une matrice 3x3, maintenant tu as une matrice 627*1401. Donc il suffit de remplacer les 3 par le nombre correct de colonnes et de lignes :

Code :

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

A=[kron(x,ones(1,627))' repmat(y,1,1401)' z'];

[romwn]

ça fonctionne !

Pas exactement avec ces valeurs, il faut juste inverser 627 et 1401 de façon à avoir :

Code :

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

1
2
 
A=[kron(x,ones(1,1401))' repmat(y,1,627)' z'];

Je peux enfin avancer dans mon application, merci à tous, et en particulier à Paul-Ceria