Discussion :

[Casimirs]

Bonjour à tous,

Je souhaiterais visualiser une fonction f(x,y,z) en utilisant l'outil slice de matlab. Voici mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
[x,y,z] = meshgrid(-10:1:10,-10:1:10,-10:1:10);
 
if (x > 0)
v = x.^2+y.^2+z.^2;
else 
v = z;    
end
 
xslice = []; yslice = 0; zslice = [];
slice(x,y,z,v,xslice,yslice,zslice)
colormap hsv

Ma fonction vaut en effet x^2+y^2+z^2 lorsque x >0 et vaut z pour x < =0.

Seule la valeur de z est plotée dans le slide...Quelqu'un peut t-il me dire pourquoi la valeur de ma fonction dans le domaine x >0 n'est t-elle pas plotée ?

merci

[Invité]

Bonjour,

Ma fonction vaut en effet x^2+y^2+z^2 lorsque x >0 et vaut z pour x < =0.

Non ton code ne fait pas du tout cela.
Ici x est un tableau 3D contenant des valeurs positives et négatives, et écrire if x > 0 revient à à dire si TOUS les éléments de x sont positifs. Ce qui est bien évidemment faux, c'est v = z; qui est exécuté indifféremment du signe de x, et v est donc égal à z.

En utilisant l'indexage logique :

Code :

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

1
2
3
xp = x > 0;
v = z;
v(xp) = x(xp).^2+y(xp).^2+z(xp).^2;