【MATLAB】一组含几个bar的柱状图的绘制以及errorbar的添加

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
2 X- V; Y# n$ i/ ^& B8 g% 数据
1 |' ]  t$ N8 W/ [" h! b6 d( A3 d& T2 Nvolume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
; k7 v% c. h% A0 m% L, Yvolume_std=[0.65,0.17;
                     0.35,0.14;
6 n5 J! S% ~: _  e* H6 T                     0.44,0.13];
%绘图                     
+ e4 U( m& p, W7 ^& Xclose all;figure;
$ l7 y1 Z# a5 V; f1 }1 v9 Uh=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
& X0 |3 o$ e7 l: ~set(h(2),'facecolor','k')        % 第二列数据视图颜色
$ a1 g) H) I  u' x
8 z! q7 E& i6 h4 W, s4 q
ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
: p  E. @9 a/ w9 L' E+ Dgroupwidth =min(0.8, nbars/(nbars+1.5));

* b0 ]" R5 x9 k$ a7 \  P
% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
. Y( M7 V& A+ k( \( z6 P/ Q+ E1 pfor i = 1:nbars
5 H- P# W1 A/ P) f# C/ j2 o0 ^    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
/ [, ]+ j2 b6 \3 G9 {    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end

2 R8 h: I& b0 W8 J$ c& C
, T1 v5 t7 T+ p% gset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
( @+ m' H5 m+ @8 x1 T1 @2 O; }5 lylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
+ k/ v" O+ q! f' F8 K& |xylabel(gca,' ','Volume[Sv]')
8 t# w# s/ D% X- z- [% @legend('data1','data2','location','NorthEast')
0 n* M0 _; l7 a

以上实例可以参考使用。


errorbar的局部调整：
7 ?( X8 ^+ {  x% O, p0 |1.头部宽度调整
% Create errorbar
. Q$ t; J( C" G4 }2 cX = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
; R/ U, n$ b3 y# e& ]6 _) m( }xlength = 0.2;
4 R, N+ a4 j! A# X% Make horizontal lines with 'line'
for k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
0 w- O9 c. Q, {  N8 H; b y_h = [Y(k) + E(k), Y(k) + E(k)];
9 {: m6 h5 c8 p4 }* q' O' s line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
end
  c! v& K. K% S参考：www.52ocean.cn
# _8 f1 ]6 N, a  K+ c1 V! k: Q

                    

1 l2 L/ K7 W$ W$ A8 o  d                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。