Discussion :

[Clem44120]

Bonjour à tous,

Je me permets de poster un message sur le forum car j'avoue que j'ai un peu de mal avec l'utilisation de Matlab.
Je cherche à modéliser une forme (écaille).
J'ai donc trouvé une fonction qui correspondait à peu près à la forme recherchée seulement je suis bloqué pour la tracer sous matlab.
La fonction en question est : Z = ( -X.^2 / (0.5 - abs(Y) ).^2 )+ 1 ; Quand je la rentre sur matlab cela ne fonctionne pas.
J'ai utilisé la commande [X,Y]=meshgrid(x,y) en ayant au préalable défini x et y.
J'ai ensuite utilisé mesh(X,Y,Z).

Je ne sais pas si mon message est très clair..
Merci d'avance.

[Jerome Briot]

cela ne fonctionne pas.

Peut-on en savoir plus ?

[Clem44120]

Voici le petit code que j'ai rentré :

Code :

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

1
2
3
4
5
6
clear all
x=0:1:100;
y=0:1:100;
[X,Y]=meshgrid(x,y);
Z = ( -X.^2 / (0.5 - abs(Y) ).^2 )+ 1 ;
mesh ( X, Y, Z)

Et en réponse :

Warning: Matrix is singular to working precision.

Sur la ligne 5 donc sur :

Code :

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

Z = ( -X.^2 / (0.5 - abs(Y) ).^2 )+ 1 ;

[phryte]

Bonjour,

Peut-être :

Code :

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

Z = ( -X.^2 ./ (0.5 - abs(Y) ).^2 )+ 1 ;

[Clem44120]

En effet... Erreur bête.
Je n'obtiens pas le résultat escompté mais au moins j'ai un tracé !
Merci beaucoup.

[Clem44120]

[Clem44120]

[Clem44120]

Je suis censé obtenir la 1ère des deux courbes que j'ai posté.. (La personne l'ayant obtenu avait utilisé la même fonction que moi)
J'ai utilisé pas mal de commandes (help surf, help mesh...)
Je ne vois pas comment l'obtenir (j'obtiens la deuxieme capture d'écran).
Si jamais vous avez des idées..

[Jerome Briot]

En ajustant légèrement les paramètres d'affichage, MATLAB trace la même surface que Wolfram Alpha.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
figure
 
x = -100:4:100;
y = -100:4:100;
[X,Y] = meshgrid(x,y);
 
Z = -X.^2 ./ (0.5 - abs(Y)).^2 + 1;
 
surf(X,Y,Z)
axis tight
 
zlim([-100 max(Z(:))])

Voir aussi en jouant sur les limites des axes avec les fonctions xlim, ylim et zlim

[Clem44120]

Merci beaucoup pour les réponses.
Je vais me pencher là dessus.. Impossible d'obtenir la figure verte, je ne sais pas comment il a pu faire.
Bonne journée