Discussion :

[azboul]

Bonsoir,
Voila je dois dans le cadre d'un projet, afficher dans un premier temps un tore de revolution dont l'équation paramétrique est :

x = (c+a*cos(v))*cos(u)
y = (c+a*cos(v))sin(u)
z = a*sin(v)

avec u,v dans [0,2pi]..

J'ai donc procédé de la sorte :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
u = 0:0.001:2*pi;
v = 0:0.001:2*pi;
a = 4;r = 2;
 
x = (a + r.*cos(u)).*cos(v);
y = (a + r.*cos(u)).*sin(v);
z = r.*sin(u);
 
plot3(x,y,z)

Seulement le mieux que j'obtienne c'est une courbe, en fait je ne sais pas quelle fonction utiliser j'ai essayé plot3, surf et surface.
Merci d'avance

[Jerome Briot]

Utilise cette contribution du FEX : TubePlot

Code :

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

1
2
3
4
5
6
7
8
9
10
u = 0:0.001:2*pi;
v = 0:0.001:2*pi;
a = 4;r = 2;
 
x = (a + r.*cos(u)).*cos(v);
y = (a + r.*cos(u)).*sin(v);
z = sin(u);
 
tubeplot([x;y;z],r)
axis vis3d

[azboul]

OK merci pour ta réponse, le problème c'est que le résultat ne correspond pas vraiment à un Tore.
Voilà ce que j'obtiens : http://img47.imageshack.us/img47/9250/tubent2.png
comparé à un tore de révolution classique : http://upload.wikimedia.org/wikipedia/fr/1/19/Tore.png
J'ai vérifier mon équation est bonne.

[Jerome Briot]

Je n'avais même pas pris le temps de regardé les équations... désolé

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
u = linspace(0,2*pi,50);
v = linspace(0,2*pi,50);
a = 4;r = 2;
 
[u,v]=meshgrid(u,v);
 
x = (a + r.*cos(u)).*cos(v);
y = (a + r.*cos(u)).*sin(v);
z = r.*sin(u);
 
surf(x,y,z)
 
axis equal vis3d

[azboul]

C'est pas grave, je te remercie bien pour ton aide rapide, ça marche nickel.
Bonne soirée