nuage de points d'une sphere sur matlab

Version imprimable

29/09/2008, 17h34
lion12

nuage de points d'une sphere sur matlab

salut
mon probleme est comme suite:
je veut en premiere lieu programmé les equation des formes classique (sphere,cube ,parllelepede .....)
1- je veut recipuré les points x,y,z
2-tracé en 3d la forme souhité
merci infiniement
29/09/2008, 17h41
mr_samurai

Salut,

tu bloque sur un passage particulier ?
t'as vu les fonctions SPHERE et CYLINDER ?

++
30/09/2008, 10h23
Jerome Briot

A quoi cela va-t-il servir ?
Quel est le contexte de l'étude ?
30/09/2008, 23h40
lion12

récipuration d'un fichier en format .txt ou .xyz et par la suite la recostruction
de la forme CAO
01/10/2008, 09h25
Jerome Briot

Et donc... c'est quoi le problème au juste ?
01/10/2008, 12h00
lion12

le probleme est dans le syntaxe du programme
j'ai essayé mais ca marche pas
01/10/2008, 12h02
Jerome Briot

Montre nous le code et dis nous clairement ce que tu cherches à faire
01/10/2008, 14h09
lion12
voila par exemple le programme d'une sphere
Code:

1 2 3 4 5 6 7 A=[]; for x=-20:1:20 for y=-20:1:20 z=sqrt(x.^2+y.^2).^0.5; A=[A;x,y,z]; end end
mais quand je récupère le nuage de points je trouve que ça ne ressemble pas
à une sphère
01/10/2008, 14h36
Caro-Line

J'ai de gros doutes sur ton équation...
Même si mes souvenirs sont rouillés.

Déjà là :

Code:

sqrt(x.^2+y.^2).^0.5

Tu appliques 2 fois la racine carrée : une fois avec sqrt et une fois avec .^0.5.

De plus tu ne peux pas y arriver avec juste un vecteur x et un vecteur y, puisque pour une sphère, pour un couple (x,y) donné tu peux avoir 2 z différents.
Exemple : pour une sphère de rayon 1 centrée en (0,0,0), si x,y=0,0 on a z=-1 et z=1.

Et dans ton équation d'ailleurs il semblerait que le rayon ne soit pas pris en compte.
Mais bon tout ça est plus un problème de mathématiques au départ que de MATLAB.

As-tu regarder la fonction SPHERE ce serait plus simple ?

Salut,

j'ai 3 remarques :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
 
%-- 1 
% ici tu fait (x^2 + y^2).^(1/4)
sqrt(x.^2+y.^2).^0.5 
 
%-- 2 
% Pour un sphère :
Z = sqrt( R^2 -X^2 - Y^2);
 
%-- 2
% Il faut  rajouter -Z sinon tu obtiens un hémisphère

Essayes ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
clc
 
xx = -20:0.5:20;
yy = -20:0.5:20;
R = 20;
 
[X Y] = meshgrid(xx,yy);
 
Z = sqrt(R^2 - X.^2 - Y.^2);
 
Z(imag(Z)~= 0) = nan;
 
surf(X,Y,Z)
axis equal
hold on
surf(X,Y,-Z)

Citation:

Envoyé par mr_samurai

Salut,

j'ai 3 remarques :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
 
%-- 1 
% ici tu fait (x^2 + y^2).^(1/4)
sqrt(x.^2+y.^2).^0.5 
 
%-- 2 
% Pour un sphère :
Z = sqrt( R^2 -X^2 - Y^2);
 
%-- 2
% Il faut  rajouter -Z sinon tu obtiens un hémisphère

Essayes ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
clc
 
xx = -20:0.5:20;
yy = -20:0.5:20;
R = 20;
 
[X Y] = meshgrid(xx,yy);
 
Z = sqrt(R^2 - X.^2 - Y.^2);
 
Z(imag(Z)~= 0) = nan;
 
surf(X,Y,Z)
axis equal
hold on
surf(X,Y,-Z)

mais mon frere votre programme donne un demi sphere

02/10/2008, 00h19
mr_samurai

Citation:

Envoyé par lion12

mais mon frere votre programme donne un demi sphere et non complet

Combine le code avec les remarques pour faire ton propre programme. Même si le code me parais assez explicite :aie: .

Tu bloques sur un point particulier?

++
02/10/2008, 07h41
Caro-Line

Et puis modifies les limites des axes aussi :roll:

Code:

set(gca,'zlim',[-20,20])

Ceci dit, elle a des trous cette sphère Mr Samurai :mouarf:

Et donc sinon la fonction SPHERE ne te plait pas ?
Parce qu'au moins elle elle fait des jolies sphères.
02/10/2008, 10h17
mr_samurai

Citation:

Envoyé par Caro-Line

Ceci dit, elle a des trous cette sphère Mr Samurai :mouarf:

Tu as raison :aie: , mais si le but est de récupérer un nuage de point, cela devrait suffire.

Sinon, SPHERE est bien la solution optimale ;) .

++
02/10/2008, 13h54
lion12

MERCI POUR L'aide
02/10/2008, 14h30
Caro-Line
Ce que tu veux n'est toujours pas clair...
Qu'est-ce qui ne te va pas dans ce que tu obtiens ?
A quoi veux-tu que tes données ressemblent ?

Je le redemande : pourquoi n'utilises-tu pas la fonction SPHERE puisque ce que tu veux c'est juste les points ?
Exemple :
Code:

1 2 3 4 [x,y,z]=sphere(40);%donne une sphere de rayon 1 centree en (0,0,0) %donc on multiplie par 20 pour avoir un rayon 20 x=20*x;y=20*y;z=20*z; surf(x,y,z);
02/10/2008, 14h33
mr_samurai
Complément au code de Caroline :
Code:

1 2 Sortie = [x(:) y(:) z(:)];
++
02/10/2008, 14h37
Caro-Line

@samurai : hum, t'es sur ?
Parce que SPHERE retourne quand même 3 matrices et non pas des vecteurs, et si on fait comme tu dis,

Code:

surf(Sortie)

ne donne plus du tout la même chose....
Ou alors j'ai rien compris :aie:

Mais c'est vrai qu'il a l'air de vouloir une matrice N*3, donc SPHERE ne conviendrait pas.

Code:

1 2 3 4 5 6 7 [x,y,z]=sphere(40); x=20*x; y=20*y; z=20*z; Sortie = [x(:) y(:) z(:)]; scatter3(Sortie(:,1), Sortie(:,2), Sortie(:,3))
chaque ligne de Sortie correspond à un point de la sphère.

++

02/10/2008, 18h53
lion12

MERCI POUR L'aide mes freres
le probleme est resolu