Discussion :

[funkyplayer]

Bonjour,

Voici un problème relativement simple je pense mais que je n'arrive pas à résoudre après l'utilisation intensive de l'aide Matlab en ligne et des conversations sur ce forum.

J'ai un nuage de points qui ont chacun des coordonnées (x,y). Ces points ont chacun un poids z. J'aimerais utiliser la fonction contour pour représenter la projection dans le plan (x,y) de mes différents poids (avec des couleurs bien sûr).

Dans les exemples que j'ai pu voir précédemment, il y avait une relation z=f(x,y). Ici ce n'est pas le cas, mes "poids" sont totalement indépendants de x et y.
J'ai essayé pas mal de choses mais ça ne fonctionne pas.

Voici par exemple ce que j'ai codé :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
load 'x1.txt'
load 'y1.txt'
load 'pression.txt'
 
x1=-200:1:1500;
y1=-400:1:1000;
[X,Y] = meshgrid(x1,y1);
Z=pression;
contour(X,Y,Z)
xlabel('x');ylabel('y');zlabel('z')

Merci à l'avance de votre aide.

[Jean Dumoncel]

Bonjour,

qu'est ce qui ne fonctionne pas? Il y a un message d'erreur? Le résultat obtenu n'est pas celui que tu attendais?

contour trace les isolignes sur z, est-ce que c'est ce que tu veux faire?

[funkyplayer]

Oui je veux tracer les isolignes sur z.

Bon en fait j'ai trouvé avec le code suivant bien connu :

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
 
 
% Determine the minimum and the maximum x and y values:
xmin = min(x1); ymin = min(y1);
xmax = max(x1); ymax = max(y1); 
 
 
% Define the resolution of the grid:
xres=10;
yres=15;
 
 
% Define the range and spacing of the x- and y-coordinates,
% and then fit them into X and Y
xv = linspace(xmin, xmax, xres);
yv = linspace(ymin, ymax, yres);
[Xinterp,Yinterp] = meshgrid(xv,yv); 
 
 
% Calculate Z in the X-Y interpolation space, which is an 
% evenly spaced grid:
Zinterp = griddata(x1,y1,z1,Xinterp,Yinterp); 
 
 
% Generate the contour plot : 
figure(2)
contour(Xinterp,Yinterp,Zinterp)
xlabel('x(pixels)');ylabel('y(pixels)');

Merci !