Discussion :

[piploum]

Bonjour la communauté,

A partir d'un vecteur de départ Mat=[0.16;0.06;0.04;0.03;0.56;25.2;...], je cherche à connaître les effectifs de Mat par tranche [0.01,0.1,1,10,100]
Et donc tracer un histogramme (données de Mat en Abscisse selon une échelle logarithmique, et effectif en ordonnée).
J'ai cru comprendre que je pourrais m'en sortir avec hist ou histc et/ou bar? plus semilogx mais j'ai beau faire des essais, je n'arrive jamais avoir l'histogramme que je désire...
un simple

Code :

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

histogram(log(Mat))

s'en rapproche...
Des conseils?

merci d'avance

[phryte]

Bonjour,

Un exemple :

Code :

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

1
2
3
4
5
6
Mat=[0.16;0.06;0.04;0.03;0.56;25.2;1.5;12;3;0.5];
[n, xout] = hist(Mat);
bar(Mat,n,'barwidth',3);
set(gca, 'XScale', 'log')
set(gca,'FontSize',15, 'XLim',[0.01,100]);
grid

[piploum]

Merci phryte
ça marche sur cet exemple mais pas sur ma matrice qui est beaucoup plus grande: taille: 21000*1
puis vu que n et xout sont de taille 10,

Code :

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

bar(Mat,n)

rale:

Code :

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

1
2
3
4
5
Error using bar (line 172)
X must be same length as Y.
 
Error in Mafonction(line 76)
bar(Mat,n);

[phryte]

Bonjour,
Un autre exemple :

Code :

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

1
2
3
4
5
6
7
8
9
clear
nr=10000;
Mat=10*rand(1,nr);
[n,xout]=hist(Mat);
%n=n/max(n);% Normalisation
bar(xout,n)
set(gca, 'XScale', 'log')
set(gca,'FontSize',15, 'XLim',[0.01,100]);
grid

[le fab]

salut

je ferais plutôt qqch comme :

Code :

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

1
2
histogram(Mat,[1e-2 1e-1 1e0 1e1 1e2])
set(gca, 'XScale', 'log')

Fabien

[Jerome Briot]

Une légère amélioration au dernier code donné par Fabien :

Code :

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

histogram(Mat, logspace(-2, 2, 5))

ou encore mieux :

Code :

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

histogram(Mat, 10.^(-2:2))

Et pour obtenir directement les bornes min et max en fonction des valeurs contenues dans Mat :

Code :

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

1
2
3
4
e(1) = floor(log10(min(Mat(:))));
e(2) = ceil(log10(max(Mat(:))));
 
histogram(Mat, 10.^(e(1):e(2)))