Discussion :

[Sachmot]

Bonjour à tous,


Je bloque depuis un moment sur problème et je ne trouve pas la solution.

Contexte:
Je tente de visualiser les tensions exercées sur une surface ronde. Ma tension ne dépends que du rayon et j'en connais déjà les valeur (qui forme une courbe discrète).

en schéma, se que je désirerais :

Problème
Je désire faire faire tourner ma courbe autour de l'axe y. J'ai essayé avec la fonction cylinder, le résultat me donne une rotation autour de x.


code

Code :

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

1
2
3
4
5
x = [0 2.5 5 7.5 10 12.5 15 17.5];
y = [12 11 10 10 10 9 11 30];
 
[X, Y, Z] = cylinder(y,500);
surf(X,Y,Z); xlabel('X'), ylabel('Z'), zlabel('Y');

resultat


Une idée? Merci d'avance en tout cas


Sach

[Jerome Briot]

Il suffit d'utiliser la définition du passage des coordonnées cylindriques en coordonnées cartésiennes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
r = [0 2.5 5 7.5 10 12.5 15 17.5];
z = [12 11 10 10 10 9 11 30];
 
n = 30;
th = linspace(0,2*pi,n);
 
xc = r(:)*cos(th);
yc = r(:)*sin(th);
zc = repmat(z(:),1,n);
 
figure
surf(xc,yc,zc,'facecolor','interp')

Si tu souhaites une répartitions plus progressive des couleurs, il faudra utiliser des points intermédiaires interpolés le long de z.

Par exemple :

Code :

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

1
2
3
4
5
r = [0 2.5 5 7.5 10 12.5 15 17.5];
z = [12 11 10 10 10 9 11 30];
 
ri = linspace(0,17.5,50);
zi = interp1(r,z,ri);

[Sachmot]

Wow merci beaucoup Dut pour la rapidité de ta réponse!

Ca marche du tonnerre ==> [résolut]


Tip top votre forum