Discussion :

[gps15000]

Bonjour,

je possede une matrice xyz composée de 3 colonnes x, y et z.
Je souhaiterai tracer tous mes points et obtenir une surface composée des points de ma matrice.

J'ai trouvé cette exemple sur le site de matlab.

Code :

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

1
2
3
4
5
6
7
[X,Y] = meshgrid(-3:.25:3);
Z = peaks(X,Y);
[XI,YI] = meshgrid(-3:.125:3);
ZI = interp2(X,Y,Z,XI,YI);
mesh(X,Y,Z), hold, mesh(XI,YI,ZI+15)
hold off
axis([-3 3 -3 3 -5 20])

J'ai voulu remplacé par mes variables.
A la premiere ligne, j'ai mis

Code :

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

[X,Y] = meshgrid(x,y);

Matlab me renvoie un probleme de mémoire.

Code :

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

1
2
3
4
5
6
7
   ??? Maximum variable size allowed by the program is exceeded.
 
Error in ==> meshgrid at 44
    xx = xx(ones(ny, 1),:);
 
Error in ==> test7 at 3
[X,Y] = meshgrid(x,y);

Mes matrices x et y sont de dimension:164953*1 chacune
Je suis un peu bloqué

[Jerome Briot]

Je suis un peu bloqué

MATLAB aussi

=> :Comment déterminer approximativement la quantité de mémoire nécessaire pour le stockage d'une variable ?

Et aussi :

• http://www.mathworks.com/support/tec...1100/1107.html
• http://www.mathworks.com/support/bugreports/300319

[gps15000]

j'ai reduit la taille de ma matrice. meshgrid continue de bloquer avec la taille.

Existe t il une autre methode pour tracer une surface définie par les points d'une matrice 3d?

merci

[Jean Dumoncel]

Bonjour,

164953 données dans chaque direction x et y, ça fait à peu près 30 milliards de points . AS-tu réellement besoin d'autant de points? Pour commencer si tu faisais des tests avec quelques millions (voir centaine de millier) de points pour voir le résultat?

[gps15000]

j'étais passé a 22640 points dans ma réduction.

j'ai décidé de réécrire mon code et de spécifier les paramètres manuellement pour meshgrid:

Code :

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

1
2
3
4
[X,Y] = meshgrid(200:0.1:400, -30:0.1,1); %X et Y sont de dimension 31*2001
Z=z; % Z est de dimension 22644x1
[XI,YI] = meshgrid(200:0.05:400, -30:0.05,1); % XI et YI sont de dimension 31x4001
ZI = interp2(X,Y,Z,XI,YI);

quand j'execute mon code.

Code :

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

1
2
??? Error using ==> interp2
When Z is a vector, X and Y must also be vectors.

X et Y doivent donc etre des vecteurs mais générés par meshgrid.
Meshgrid ne genere que des array d'apres la doc?

[Jean Dumoncel]

Il y a quelques points qui ne sont pas clairs :

• dans ton titre, tu parles de tracer une courbe, et dans ton post tu parles d'une surface??
• Si tu avais uniquement 22640 points qui définissent ta surface, meshgrid aurait dû marcher, à moins que tu ais très très peu de ram sur ton ordi...
• pourquoi veux-tu passer par une interpolation de tes données?