Discussion :

[Udani]

Bonjour,

J’ai réuss à avoir, à l’aide de matlab, des données sous forme d’un tableau de 3 dimensions ( trois colonnes (x, y, z)).

Mon problème, c’est que j’ai un quatriéme vecteur de valeurs ( par exemple représenté par un vecteur à 1 dimension nommé par exemple « error » ).

J’aimerais faire une représentation graphique qui fait intervenir le vecteur à 3 dimensions (x, y, z) et le vecteur error ensemble. Autrement dit, comment representer graphiquement ( x, y, z, error) ensemble??

Merci de me conseiler,

[duf42]

Bonjour,

Une solution possible serait de tracer tes points en 3D (avec SURF, PLOT3 ou autre) et de leur attribuer à chacun une couleur qui correspond à l'erreur (la 4°dimension de ton problème).

Après, il faut voir si la représentation correspond à ce que tu recherches...

Bonne journée,
Duf

[telecom85]

bonjour;
j'ai pas très très bien compris ton problème.
Mais si tu veux par exemple représenter un cube en 3D, tu peut le faire avec la commande patch mais dans ce cas il faut utiliser 4 variables(x,y,z,t) au lieu de trois.
c'est ça ton problème????

[Udani]

Bonjour duf42, telecom85,
Merci, vos conseils m'ont été très utiles.

Mes données sont sous forme de 4 vecteurs en colonne. J'ai essayé de travailler avec la commande

Code :

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

patch

J'ai reussi à transformer les vecteurs (x,y,z) en matrice à l'aide de la commande

Code :

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

meshgrid ( x, y, z)

Code :

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

1
2
3
4
[x, y, z]=meshgrid(C1, C2, C3);
p1 = patch(isocaps(x,y,z,V,0.1), 'FaceColor', 'interp','EdgeColor','none');
p2 = patch(isosurface(x,y,z,V,0.1) ,'FaceColor','red' ,'EdgeColor','None');
isonormals(x,y,z,V,p2);

Message erreur:

Code :

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

1
2
??? Error using ==> isocaps
V must be a 3D array.

je sais que mon programme n'accepte pas le vecteur V et exige une Matrice V . Mais pourquoi et comment?

Merci de votre aide!

[telecom85]

salut;
bon je n'arrive a comprendre ce que tu veux exactement ,
mais voila un exemple qui te permet de dessiner un cube en (3D), mais pour le faire il faux que tu déclare 4 vecteurs et non pas trois et a l'aide de la commande patch, peut être sa va t'aide .

Code :

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

1
2
3
4
5
6
7
X=[0 0 4 4,0 0 4 4,4 0 0 4,4 4 4 4,0 0 4 4,0 0 0 0];
Y=[0 4 4 0,0 0 0 0,0 0 4 4,0 0 4 4,4 4 4 0,0 4 4 0];
Z=[0 0 0 0,0 4 4 0,4 4 4 4,4 0 0 4,4 0 0 0,0 0 4 4];
s=[1 1 1 1,1 1 1 1,1 1 1 1,1 1 1 1,1 1 1 1,1 1 1 1];
figure
p=patch(X,Y,Z,s);
axis([-10 10 -10 10 -10 10]),title('cube initial');%% affichage du cube

[Jerome Briot]

A ta place, j'utiliserais simplement SURF

Si tes données sont réparties sur une grille régulière pleine, il suffit d'utiliser RESHAPE pour mettre en forme les 4 vecteurs. Sinon, il faut utiliser GRIDDATA.

Fais une recherche sur ce forum, nous avons déjà répondu à ce genre de questions...

[Udani]

Bonsoir,

ce que j'aimerais faire, c'est de travailler avec les données que j'ai rassemblés. Mes données sont sous formes de 4 vecteurs en colonne. J'aimerais représenter (x, y, z) sur 3 axes ( càd en 3D) avec t une 4 éme valeur dont dépend les 3 variables.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
x    y   z    t
6    3   2   2
4    2   4   3
6    3   8   0
2    4   2   5
3    4   3   6
1    3   8   0
3    7   0   2
2    3   3   1
5    9   4   9
2    2   1   3

le problème, c’est que les commandes que vous m'avez donné travaillent avec des matrices, et moi j'ai à ma disposition des vecteurs colonnes.
A+

[Jerome Briot]

Essaie ceci :

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
xyzt = [6    3   2   2
    4    2   4   3
    6    3   8   0
    2    4   2   5
    3    4   3   6
    1    3   8   0
    3    7   0   2
    2    3   3   1
    5    9   4   9
    2    2   1   3]
 
 
[xi,yi] = meshgrid(linspace(min(xyzt(:,1)),max(xyzt(:,1))),linspace(min(xyzt(:,2)),max(xyzt(:,2))));
 
zi = griddata(xyzt(:,1),xyzt(:,2),xyzt(:,3),xi,yi);
ci = griddata(xyzt(:,1),xyzt(:,2),xyzt(:,4),xi,yi);
 
figure
surf(xi,yi,zi,ci,'edgecolor','none')
colorbar

[Udani]

Merci beaucoup pour vos réponses que je vais essayer d'appliquer. Je vous tient au courant.

[Udani]

Bonjour,

désolé, mais j'arrive pas toujours à trouver une solution pour mon problème.
J'ai réussi à voir les données exactes dont j'aimerais faire une représentation graphique. La représentation graphique doit être de dimension 3 pour les variables (x,y,z) et la variable t intervient en tant que la 4 éme variable, par exemple en couleur, ou sous forme d'une boule dont rayon prend la valeur t...

J'ai essayé avec la proposition de Duf et après exécution du code je reçois des erreurs:

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
 
 
function matrice_error
 
load data4d.dat;
x=data4d(:,1);
y=data4d(:,2);
z=data4d(:,3);
t=data4d(:,4);
 
Z=zeros(20,4);
 
Z=[x y z t];
 
 
[xi,yi] = meshgrid(linspace(min(Z(:,1))*100,max(Z(:,1))*100),linspace(min(Z(:,2)*100),max(Z(:,2))*100));
 
zi = griddata(Z(:,1)*100,Z(:,2)*100,Z(:,3)*100,xi,yi);
ci = griddata(Z(:,1)*100,Z(:,2)*100,Z(:,4)*100,xi,yi);
 
figure;
surf(xi,yi,zi,ci,'edgecolor','none');
colorbar;

Error:

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
EDU» matrice_error
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\matrice_error.m at line 16
Number of array elements exceeds 16384
??? Student Edition array size limitation.
 
Error in ==> C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m (linear)
On line 111  ==> tri = tri(t,:);
 
Error in ==> C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m
On line 72  ==>     zi = linear(x,y,z,xi,yi);
 
Error in ==> C:\MATLAB_SE_5.3\work\matrice_error.m
On line 16  ==> zi = griddata(Z(:,1)*100,Z(:,2)*100,Z(:,3)*100,xi,yi);

vous trouverez en piéce jointe les 4 vecteurs sous forme de matrice. voir fichier "data4d.dat"

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
  6.6666667e-002  4.1666667e-002  2.0000000e-002  0.0000000e+000
  6.0301508e-002  3.7688442e-002  2.0904523e-001  0.0000000e+000
  5.6951424e-002  3.5594640e-002  3.0854271e-001  0.0000000e+000
  5.3601340e-002  3.3500838e-002  4.0804020e-001 -3.5527137e-015
  5.0251256e-002  3.1407035e-002  5.0753769e-001 -3.5527137e-015
  4.6901173e-002  2.9313233e-002  6.0703518e-001  0.0000000e+000
  4.0201005e-002  2.5125628e-002  8.0603015e-001 -3.8524909e-001
  4.0201005e-002  2.5125628e-002  8.0603015e-001 -3.8524909e-001
  3.6850921e-002  2.3031826e-002  9.0552764e-001 -7.3045241e-001
  3.3500838e-002  2.0938023e-002  1.0050251e+000 -1.0601348e+000
  2.6800670e-002  1.6750419e-002  1.2040201e+000 -9.9903100e-001
  2.3450586e-002  1.4656616e-002  1.3035176e+000 -7.2097346e-001

Merci de vous conseils.

[Jerome Briot]

Code :

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

1
2
3
4
5
6
EDU» matrice_error
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\matrice_error.m at line 16
Number of array elements exceeds 16384
??? Student Edition array size limitation.

Le problème vient semble-t-il du fait que tu utilises la version Student de MATLAB qui semble limiter la taille des tableaux manipulables.

La représentation graphique doit être de dimension 3 pour les variables (x,y,z) et la variable t intervient en tant que la 4 éme variable, par exemple en couleur, ou sous forme d'une boule dont rayon la valeur de t...

Dans ce cas, utilises SCATTER3

[Udani]

bonjour,

malheureusement, la ou je suis maintenant je peux utiliser que matlab version étudiant. j'ai suivi ton conseil en utilisant scatter3 ( en remplaçant surf, je reçois la même erreur..

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
EDU» matrice_error
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\matrice_error.m at line 16
Number of array elements exceeds 16384
??? Student Edition array size limitation.
 
Error in ==> C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m (linear)
On line 111  ==> tri = tri(t,:);
 
Error in ==> C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m
On line 72  ==>     zi = linear(x,y,z,xi,yi);
 
Error in ==> C:\MATLAB_SE_5.3\work\matrice_error.m
On line 16  ==> zi = griddata(Z(:,1)*100,Z(:,2)*100,Z(:,3)*100,xi,yi);
 
EDU»

si seulement je pourrais en être sûr que le problème est dû à la version de matlab???? je pourrais faire autre chose en attendant de rentrer à la faculté...

[Jerome Briot]

C'est normal que tu obtiennes la même erreur.

Tu peux réduire la taille de ton problème.

Par défaut, LINSPACE produit 100 valeurs régulièrement espacées, mais tu peux spécifier le nombre toi même :

y = linspace(a,b) generates a row vector y of 100 points linearly spaced between and including a and b.

y = linspace(a,b,n) generates a row vector y of n points linearly spaced between and including a and b.

[Udani]

Bonsoir Dut,

Code :

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

Tu peux réduire la taille de ton problème

j'ai changé le code par

Code :

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

1
2
3
4
 
x=linspace(min(Z(:,1))*100,max(Z(:,1))*100,50);
y=linspace(min(Z(:,2)*100),max(Z(:,2))*100,50);
[xi,yi] = meshgrid(x,y);

je reçoi toujours les erreurs suivantes:

Code :

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

1
2
3
4
5
6
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\error_matrix.m at line 55
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\error_matrix.m at line 56

[Jerome Briot]

Ce ne sont pas des erreurs (error en anglais) mais des avertissements (warning en anglais)

Tu as des points pour lesquels il y a plusieurs valeurs.
Dans ce cas, GRIDDATA fais la moyenne de ces valeurs pour n'obtenir au final qu'une seule valeur par point

[Udani]

j'ai encore revu le code et le resultat est dans la piéce jointe. Apparament il y'a toujours quelques chose qui marche pas dans le code.

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
function matrice_error
 
load data4d.dat;
x=data4d(:,1);
y=data4d(:,2);
z=data4d(:,3);
t=data4d(:,4);
 
Z=zeros(20,4);
 
Z=[x y z t];
 
 
[xi,yi] = meshgrid(linspace(min(Z(:,1)),max(Z(:,1)),20),linspace(min(Z(:,2)),max(Z(:,2)),20));
 
zi = griddata(Z(:,1),Z(:,2),Z(:,3),xi,yi);
ci = griddata(Z(:,1),Z(:,2),Z(:,4),xi,yi);
 
figure;
surf(xi,yi,zi,ci,'edgecolor','none');
colorbar;

Code :

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

1
2
3
4
5
6
7
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\matrice_error.m at line 16
Warning: Duplicate x-y data points detected: using average of the z values
> In C:\MATLAB_SE_5.3\toolbox\matlab\polyfun\griddata.m at line 58
  In C:\MATLAB_SE_5.3\work\matrice_error.m at line 17
EDU»

[Udani]

Bonjour,

j'ai fini par avoir une representation graphique( voir piéce-jointe)

Code :

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

1
2
3
4
5
6
7
load data4d.dat;
x=data4d(:,1);
y=data4d(:,2);
z=data4d(:,3);
t=data4d(:,4);
scatter3(x,y,z,200,t,'filled');
colorbar

Maintenant il me faut une interpretation des couleurs!!!
J'aimerais par exemple, selectionner les parametres (x,y,z) pour lesquels la couleurs est bleu. C'est possible

[Jerome Briot]

Maintenant il me faut une interpretation des couleurs!!!

C'est à dire ?

[Udani]

la 4 variable représente en fait, l'erreur absolue entre deux méthodes numériques appliquées pour calculer quelque chose ( un produit financier, une option). D'aprés le resultat graphique, les points (x,y,z) dont la couleur est bleu présentent une erreur absolue plus intéressante. Est t- il possible de localiser avec plus d’éxactitude ces points.