Discussion :

[thomMonteillet]

Bonjour,

je viens vous voir car j'aimerai vous soumettre une fois de plus un petit problème que je rencontre.

j'aimerai faire un disque comme sur la pièce jointe à partir des coordonnées des points (que je connais) et j'aimerai attribuer à chaque point une valeur afin d'avoir un graph en 3D.

Le soucis est que je ne comprends pas comment faire...

Je pense qu'avec surfc on peut faire ceci mais j'ai regardé l'aide et ça me parait bizarre :

Code :

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

surfc(x,y,z)

je pensais que les coordonnées de mes points devaient etre dans les vecteur X et Y et les valeurs associées dans Z mais ce n'est pas possible puisque X et Y n'ont pas la meme taille !

- pourriez vous me dire s'il vous plait comment faire ce graphique en 2D ?

ps: pour le moment j'ai 3choses :
1) un vecteur X qui contient 13 ordonnées (car j'ai 13points)
2) un vecteur Y qui contient 13 abscisses (car j'ai 13points)
3) un vecteur Z qui contient 13 valeurs d'intensité (car j'ai 13points)

je vous remercie

[Jerome Briot]

Comment se fait-il que les variables X et Y n'aient pas la même taille ?

Pourrais-tu nous fournir des valeurs (dans un fichier attaché) ?

[thomMonteillet]

Bonjour DUT

Comment se fait-il que les variables X et Y n'aient pas la même taille ?

dans mon cas X et Y ont bien la meme taille mais c'est dans l'aide de surfc qu'il y a marqué que ces vecteurs non pas la meme taille et c'est ça que je ne comprends pas.

Pourrais-tu nous fournir des valeurs (dans un fichier attaché) ?

oui, bien sur.

en PJ il y a la definition d'une disque avec 8 rayons (j'ai decoupé le disque en angles de 45°) avec 4 points par rayon (centre compris).

je m'attendais à ce que ceci fonctionne par exemple :

Code :

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

1
2
3
4
5
6
7
8
 
donnees=load('surfaceNEW.txt');
X=donnees(:,1);
Y=donnees(:,2);
Z=donnees(:,3);
figure (1)
plot3(X,Y,Z,'*r')
surfc(X,Y,Z)

bien entendu le graph avec plot3 fonctionne mais avec surfc j'ai ce message d'erreur :

Code :

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

The surface Z must contain more than one row or column.

je comprends ce que ça veut dire mais je ne comprends pas pourquoi je dois avoir une matrice et comment la construire

j'ai testé là aussi en construisant une matrice mais ça ne me donne pas du tout ce que j'attendais mais des surfaces qui s'entrecroisent

[thomMonteillet]

peut etre que le mieux serait d'utiliser path et vertice ?

par contre puis je donner une hauteur et associé à ceci une couleur comme avec surfc ?

[membreComplexe12]

ça m'interesse egalement ce probleme, je suis incapable de le faire avec "surfc" mais avec "patch" j'ai reussi à faires quelques trucs.

j'ai repris ton fichier et j'ai remplacé la 3eme colonne par une colonne de 0 (pour la coordonnée Z) et ensuite j'ai fais :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
clear all;close all;clc;
noeuds=load('surfaceNEW.txt');
nbRayons=6;
nbNodeR=4;
face1=[1 2 5];
face2=[2 3 6 5;
       3 4 7 6];
patch('Vertices',noeuds,'Faces',face1)
patch('Vertices',noeuds,'Faces',face2)

ça me fait bien une portion de disque (par contre j'ai un petit soucis car c'est pas dans le plan...)

[Invité]

Bonsoir,

Une façon simple d'obtenir ce que tu as actuellement :

Code :

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

1
2
[X Y Z] = cylinder(3:-1:0,9);
surfc(X,Y,Z)

Tu pourras t'inspirer de ces données

[thomMonteillet]

merci, en effet ce que tu m'as montré fonctionne très bien, j'ai bien une portion de disque. par contre je n'ai pas le soucis que tu mentionnais, je pense que tu avait conservé le mauvais fichier de départ.

je me suis inspiré de ce que tu m'as montré et voici ce que j'arrive à faire :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
noeuds1=[0 0 0;
         1 0 0;
         0 1 0];
noeuds2=[1 0 0;
         2 0 0;
         0 2 0;
         0 1 0;
         3 0 0
         0 3 0]; 
face1=[1 2 3];
face2=[1 2 3 4;
       6 3 2 5];
patch('Faces',face1,'Vertices',noeuds1,'FaceVertexCData',[1 1 1]','FaceColor','interp','FaceVertexAlphaData',0.1)
patch('Faces',face2,'Vertices',noeuds2,'FaceVertexCData',[1 2 2 1 3 3]','FaceColor','interp','FaceVertexAlphaData',0.1)

c'est pas mal mais je n'arrive pas à mettre :
-> de la transparence
-> et faire que l'amplitude que j'ai donné au noeuds fait que le noeud prenne de l'altitude... (il faut forcement que je modifie la coordonnée suivant "z" pour cela ?)

Bonsoir,

Une façon simple d'obtenir ce que tu as actuellement :

Code :

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

1
2
[X Y Z] = cylinder(3:-1:0,9);
surfc(X,Y,Z)

Tu pourras t'inspirer de ces données

merci winjerome ça ressemble exactement à ce que je souhaite et ça à l'air très simple !
-> il faut juste que j'arrive à comprendre comment tu as fais pour fixer le nombre de portions de la sphere et le nombre d'epaisseur.
-> il faut que j'arrive aussi à attribuer à un element une amplitude, là ça à l'air de ce faire automatiquement...

EDIT :
j'aime bien la technique avec les patch ça me parait assez efficace
par contre j'arrive pas à gérer la transparence je ne comprends pas pourquoi...