Discussion :

[haydens]

Bonsoir tout le monde

Voila, j'ai petit probleme. J'ai créé un prgramme qui me calcul la solution de l'équation de la chaleur 2D + temps.

Les solutions sont stockées dans le fichier "solution.dat". Chaque colonne correspond la solution au temps t_i.
Le fichier coordonnées.dat contient les coordonnée de chaque points du maillage.

Voila mon programme

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
 
%lecture du fichier
a=load('solution.dat','-ascii');
b=load('coordonnees.dat','-ascii');
 
%adapation afin de creer le graphe
ti = 0:0.01:1;
[XI,YI] = meshgrid(ti,ti);
 
for i=1:13
	x=b(:,1);
	y=b(:,2);
	z=a(:,i);
	ZI = griddata(x,y,z,XI,YI);
 
	%affichage
	surf(XI,YI,ZI);
	view(0,90);
	%pause
end

Qd une image apparait je peux faire apparaitre en plus une barre verticale a coté correspondant a une echelle des couleurs en cliquant sur l'icone adequate. J'aimerais avoir une commande qui fait directement apparaitre cette barre.

Il y a aussi les trait du maillage sur l'image, peut on les enlever ?

Et le plus important je voudrai fixer les couleurs.
Pour le t=0 la valeur max=30 et min =0 il fixe donc rouge pour30 et bleu pour 0 mais pour t=20 la valeur max=30 et min=20 et la il fixe le bleu pour le 20 et rouge pour 30 je voudrais qu'il reprenne la couleur d'avant en fait et que donc la valeur 20 soit en orange/rouge.

Je ne sais pas si j'ai été tres claire.

Je vous remercie d'avance !

[Jerome Briot]

Qd une image apparait je peux faire apparaitre en plus une barre verticale a coté correspondant a une echelle des couleurs en cliquant sur l'icone adequate. J'aimerais avoir une commande qui fait directement apparaitre cette barre.

Code :

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

doc colorbar

Il y a aussi les trait du maillage sur l'image, peut on les enlever ?

Mettre la propriété EdgeColor à None dans l'appel à SURF

Et le plus important je voudrai fixer les couleurs.

Voir la propriété CLim de l'objet Axes sur lequel est tracé la surface.

[haydens]

merci beaucoup

Mais j'ai aps reussi pour la propriété CLim de l'objet axes.

J'ai ajouter

Code :

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

set(gca, 'CLim', [0;30]);

mais ça marche pas tres bien...

Merci

[haydens]

Code :

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

1
2
3
4
5
6
7
	%affichage
	surf(XI,YI,ZI,'edgecolor','none');
    %surf(XI,YI,ZI);
	view(0,90);
    colorbar
    handlecolorbar=colorbar('CLim', [0 30]);
    get(handlecolorbar,'CLim')

Ne marche pas non plus

[Jerome Briot]

Essaie ceci :

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
 
%adapation afin de creer le graphe
ti = 0:0.01:1;
[XI,YI] = meshgrid(ti,ti);
 
figure
s = surf(nan(4),'edgecolor','none');
set(gca,'clim',[min(a(:)) max(a(:))]);
colorbar
view(0,90);
 
for i=1:13
 
	ZI = griddata(b(:,1),b(:,2),a(:,i),XI,YI);
 
	%affichage
	set(s,'xdata',XI,'ydata',YI,'zdata',ZI);
	pause(.5)
 
end

[haydens]

impeccable merci infiniment