Discussion :

[pasqual]

Bonjour tout le monde,

Je possède une matrice 3d (p(x,y,z) qui contient pour chaque trio (x,y,z) une valeur de pression. J'aimerais afficher le résultat de la manière suivante mais je ne m'en sors pas.



Je parviens a afficher mais seulement en coupe de la manière suivante

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
25
26
27
28
29
30
31
32
33
 
clc;
close all;
lx=8;
ly=6;
lz=4;
 
n1=1;
n2=1;
n3=1;
 
x=0:0.05:lx;
y=0:0.05:ly;
z=0:0.05:lz;
s1=size(x);
s2=size(y);
s3=size(z);
 
a=s1(2);
b=s2(2);
c=s3(2);
 
p=zeros(a,b,c);
 
for i=1:a
    for j=1:b
        for k=1:c
           p(i,j,k)=cos(n1*pi*x(i)/lx)*cos(n2*pi*y(j)/ly)*cos(n3*pi*z(k)/lz); 
        end
    end
end
 
imagesc(p(:,:,1));

J'ai farfouiller un peu la doc, mais pas trouver d'exemples clairs.

Quelqu'un a une idée?

Merci d'avance

[Jerome Briot]

Essaie ceci dans un premier temps :

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
25
26
27
28
29
30
31
32
33
34
35
36
clc;
close all;
lx=8;
ly=6;
lz=4;
 
n1=1;
n2=1;
n3=1;
 
x=0:0.05:lx;
y=0:0.05:ly;
z=0:0.05:lz;
s1=size(x);
s2=size(y);
s3=size(z);
 
a=s1(2);
b=s2(2);
c=s3(2);
 
p=zeros(a,b,c);
 
for i=1:a
    for j=1:b
        for k=1:c
           p(i,j,k)=cos(n1*pi*x(i)/lx)*cos(n2*pi*y(j)/ly)*cos(n3*pi*z(k)/lz); 
        end
    end
end
 
s=slice(p,[1 b],[1 a],[1 c]);
set(s,'edgecolor','none')
 
box on
axis equal tight vis3d

Ensuite, il faudra certainement modifier quelques propriétés de l'objet Axes pour peaufiner les proportions de l'affichage.

[pasqual]

Merci bien,

c'est exactement ce qu'il me fallait

A plus