Discussion :

[Jlien]

Bonjour a tous

Je veux cree un maillage avec des quqdrilateres sur une sphere centree a l'origine et de rayon 1
j'ai fait le code mais juste la derniere etape
j'ai les deux angles

Code :

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

1
2
theta =0:pi/3:pi ;
phi=0:pi/3:7*pi/3;

Soit

Code :

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

1
2
3
x=cos(Phi).*sin(Theta);
y=sin(Theta).*sin(Phi);
z=cos(Theta);

Je veux calculer les trois variables x, y, z seulement aux points (theta,phi) de l'emsenble qui contient seulement 14 points
{(0,Pi/3), (0,Pi), , (0,5*Pi/3),

(Pi/3,0), (Pi/3,2*Pi/3), (Pi/3,4*Pi/3), (Pi/3,2*Pi),

(2*Pi/3,Pi/3), (2*Pi/3,Pi), (2*Pi/3,5*Pi/3), (2*Pi/3,7*Pi/3),

(PI,2*Pi/3), (Pi,4*Pi/3), (Pi,2*Pi) }

puis afficher les noeuds obtenus sur une sphere
Comment peut on extraire cet ensemble des deux vecteurs theta et phi
Merci d'avance

[Jerome Briot]

Comme ceci ?

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
17
18
19
20
21
22
23
24
TP = [0, pi/3
    0, pi
    0, 5*pi/3
    pi/3, 0
    pi/3, 2*pi/3
    pi/3, 4*pi/3
    pi/3, 2*pi
    2*pi/3, pi/3
    2*pi/3, pi
    2*pi/3, 5*pi/3
    2*pi/3, 7*pi/3
    pi, 2*pi/3
    pi, 4*pi/3
    pi, 2*pi ];
 
x = cos(TP(:,2)).*sin(TP(:,1));
y = sin(TP(:,1)).*sin(TP(:,2));
z = cos(TP(:,1));
 
figure
sphere;
hold on
plot3(x, y, z, 'r*');
axis equal

[Jlien]

Bonsoir a tous

Comment a partir par example de 4 points donnees creer un quqdrilatere sur une sphere.
Peut on utiliser la liste des points donnees pour creer un maillage en forme quqdrilatere d'une sphere
Merci

[Jlien]

Bonjour
Comment tracer les lignes entre points pour former un maillage sous forme quadrilateral de la sphere
Merci d'avance

Comme ceci ?

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
17
18
19
20
21
22
23
24
TP = [0, pi/3
    0, pi
    0, 5*pi/3
    pi/3, 0
    pi/3, 2*pi/3
    pi/3, 4*pi/3
    pi/3, 2*pi
    2*pi/3, pi/3
    2*pi/3, pi
    2*pi/3, 5*pi/3
    2*pi/3, 7*pi/3
    pi, 2*pi/3
    pi, 4*pi/3
    pi, 2*pi ];
 
x = cos(TP(:,2)).*sin(TP(:,1));
y = sin(TP(:,1)).*sin(TP(:,2));
z = cos(TP(:,1));
 
figure
sphere;
hold on
plot3(x, y, z, 'r*');
axis equal