Discussion :

[Casimirs]

Bonjour à tous,

Je viens de créer un petit programme permettant de déterminer l'intégrale (en utilisant la méthode des trapèzes et après avoir crée un maillage 3D régulier à partir d'un fichier texte (x,y,z,E)) d'une fonction E(x,y,z) par rapport à la variable "y" seulement. Le résultat est une fonction E_TF(x,z) que je n'arrive malheureusement pas à visualiser en utilisant la fonction surf (je souhaiterais en effet obtenir une carte 2D donnant la valeur de E_TF dans le plan (x-z) par une échelle de couleur). La dernière partie du programme (en gras) est la partie "visualisation" de cette fonction qui ne marche pas et je pense que cela est dû à une mauvaise utilisation des matrices qui stockent les données crées après l'utilisation de griddata3...

Quelqu'un peut t-il me dire comment modifier la dernière partie de ce programme (en gras) pour pouvoir visualiser E_TF(X,Z) ?

Voici en détail le code commenté que j'ai écrit :

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
%%LECTURE DES DONNEES
load text.txt % Le fichier text (voir pièce-jointe)
x=text(:,1); % x : extrait la première colonne de text.txt
y=text(:,2); % y : ----------seconde------------                      
z=text(:,3); % z  : ----------troisième-----------
E=text(:,4); % E : -----------quatrième---------

%%CREATION DU MAILLAGE 3D ET INTERPOLATION
xlin=linspace(min(x),max(x),5); % 5 points régulièrement espacés / x
ylin=linspace(min(y),max(y),5); % ------------------------------/y
zlin=linspace(min(z),max(z),5); %------------------------------/z
[X,Y,Z]=meshgrid(xlin,ylin,zlin); % maillage 3D régulier et mise en forme 
EE=griddata3(x,y,z,E,X,Y,Z); % Interpolation (linéaire) des données sur le %nouveau maillage 3D régulier

%%INTEGRATION SUIVANT y
delta=ylin(2)-ylin(1); % donne l'écart entre les points sur l'axe y
E_TF=delta*trapz(EE); % Intégration par la méthode des trapèzes /y

%%VISUALISATION DE E_TF(x,z) EN UTILISANT SURF(*,*,*)
surf(X(1,:,1),Z(1,:,1),E_TF) %% NE MARCHE PAS
pcolor(X(1,:,1),Z(1,:,1),E_TF);
set(get(gca,'Children'), 'EdgeColor', 'none')
colorbar

Merci pour votre aide !!

[Jean Dumoncel]

Bonjour,

Comme ceci :

Code :

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

1
2
3
4
5
surfc(xlin,zlin,squeeze(E_TF)) 
hold on
pcolor(xlin,zlin,squeeze(E_TF));
set(get(gca,'Children'), 'EdgeColor', 'none')
colorbar

[Casimirs]

Merci Magelan pour cette super réponse ! Je ne connaissais pas cette fonction squeeze. Cependant j'ai dû modifier un peu l'ordre des arguments dans cette fonction (par rapport à ce que tu proposais).
Finalement pour visualiser E_TF(x,z) dans le plan (x-z), j'ai dû écrire :

Code :

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

1
2
3
4
5
surf(zlin,xlin,squeeze(E_TF)) 
hold on
pcolor(zlin,xlin,squeeze(E_TF));
set(get(gca,'Children'), 'EdgeColor', 'none')
colorbar

C'est-à-dire passer zlin en premier argument et xlin en second argument dans les fonctions surf et pcolor.

Merci aussi d'avoir ajouté les balises de code ! C'est plus propre !!

Casimir