Discussion :

[tarzan_tlse]

Bonjour à toutes et à tous !
Je rencontre actuellement un problème dans l'utilisation de la fonction "pcolor".
L'objectif est de tracer par un dégradé de couleur les effectifs des couples (gradiant de dissipation,durée) dans l'échantillon suivant:

grad. dissip(x) durée (y)
16.7961 3.0000
4.0099 2.0000
14.0942 3.0000
8.3702 1.0000
7.8334 2.0000
7.4838 1.0000
8.9951 3.0000
-20.0976 2.0000
-0.7084 1.0000
-4.9076 1.0000
-1.1695 1.0000
-21.0938 5.0000
-3.5675 1.0000

Le plot 2D obtenu avec "pcolor" est en pièce jointe.

Voici le code:

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
  % Dissipation slopes (column#1) and associated times (column#2)
B= [16.7961    3.0000
    4.0099    2.0000
   14.0942    3.0000
    8.3702    1.0000
    7.8334    2.0000
    7.4838    1.0000
    8.9951    3.0000
  -20.0976    2.0000
   -0.7084    1.0000
   -4.9076    1.0000
   -1.1695    1.0000
  -21.0938    5.0000
   -3.5675    1.0000]
 
n=hist3(B,{min(B(:,1)):5:max(B(:,1)) min(B(:,2)):1:max(B(:,2))});
n1 = n'; 
n1( size(n1,1) + 1 ,size(n1,2) + 1 ) = 0;
 
% % Generate grid for 2-D projected view of intensities
xb = linspace(min(B(:,1)),max(B(:,1)),size(n,1)+1);
yb = linspace(min(B(:,2))-0.5,max(B(:,2))+0.5,size(n,2)+1);
 
% % Make a pseudocolor plot on this grid 
pcolor(xb,yb,n1)
colormap(hot)
grid on
xlabel('Dissipation slope (W)')
ylabel('Slope time (min.)')

Le problème se situe sur la figure pour x=[-6.9 ; 2.6] et y=1, car d'après le code des couleurs il y aurait:
* 1 couple (x,y) vérifiant x=[-6.9 ; -2.149] et y=1
* 3 couples (x,y) vérifiant x=[-2.149 ; 2.6] et y=1

Or, si je reprends le tableau ci-dessus, il y a en fait:
* les couples (-3.5675,1.0000) et (-4.9076,1.0000) vérifiant x=[-6.9 ; -2.149] et y=1, soit 2 couples et non pas 1.
* les couples (-0.7084,1.0000) et (-1.1695,1.0000) vérifiant x=[-2.149 ; 2.6] et y=1, soit 2 couples et non pas 3.

D'où vient l'erreur? Merci pour votre support et votre aide...