Discussion :

[Bebs7]

Je dispose de trois vecteurs : a, b, c.
a contient les abscisses.
b contient les ordonnées.
c contient la valeur d'une fonction au points de coordonnées (a,b).

Je souhaite tracer la fonction qui va de R² dans R : c = f(a,b)

Par exemple (exemple simple):

Code :

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

1
2
3
4
 
a = [ 0 1 2 0 1 2 0 1 2 ] ;
b = [ 0 0 0 1 1 1 2 2 2 ] ;
c = [ 0 1 2 1 2 3 2 3 4 ] ;

Permettrait de représenter c = f(a,b) = a+b

[Jerome Briot]

Tu veux dire une surface plutôt qu'une courbe, non ?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
a = [ 0 1 2 0 1 2 0 1 2 ] ;
b = [ 0 0 0 1 1 1 2 2 2 ] ;
c = [ 0 1 2 1 2 3 2 3 4 ] ;
 
x = reshape(a,3,[]);
y = reshape(b,3,[]);
z = reshape(c,3,[]);
 
figure
surf(x,y,z)
axis equal

[Bebs7]

Merci, c'est exactement ce que je voulais.

PS : J'ai corrigé le titre, il s'agissait bien de surface.

[Bebs7]

Je pensais que c'était bon au début, mais l'exemple que je propose est un peu trop simple... j'ai donné neuf points, donc divisible par 3 ; et les points sont dans un certain ordre "agréable".
En réalité, c'est un peu plus bordélique...

Comment ferais-tu pour l'exemple suivant, avec des listes de 11 éléments :

Code :

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

1
2
3
a = [ 0 1 2 0 1 2 0 1 2 .5 1.2] ;
b = [ 0 0 0 1 1 1 2 2 2 .4 1.9] ;
c = [ 0 1 2 1 2 3 2 3 4 .9 3.1] ;

Je précise que dans mon cas, j'ai quelques dizaines de milliers de points...!

[Jerome Briot]

GRIDDATA est ton ami

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
a = [ 0 1 2 0 1 2 0 1 2 .5 1.2] ;
b = [ 0 0 0 1 1 1 2 2 2 .4 1.9] ;
c = [ 0 1 2 1 2 3 2 3 4 .9 3.1] ;
 
 
[ai,bi] = meshgrid(linspace(min(a),max(a),10),linspace(min(b),max(b),10));
 
ci = griddata(a,b,c,ai,bi);
 
figure
colormap(gray)
surf(ai,bi,ci);
hold on
plot3(a,b,c,'r*')
 
axis equal vis3d

Fais une recherche sur ce forum, ce problème a déjà été traité plusieurs fois...

Je précise que dans mon cas, j'ai quelques dizaines de milliers de points...!

Afin d'alléger le travail, tu peux te demander si tous ces points sont représentatifs ou si certains ne pourraient pas être mis de coté

[Bebs7]

Merci, cette fois c'est définitivement correct.