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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
! V* T$ D4 e. d, a- L% 数据
volume_mean=[0.73,0.45;
, O( Z7 T% f0 X6 f( e                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
                     0.35,0.14;
                     0.44,0.13];
%绘图                     
3 h: }1 T! @3 U+ U0 S" |) Gclose all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
  @" a) ]6 ^6 \2 j+ l. @; C; A4 mset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色
2 Y/ G1 G4 \5 x+ n5 z$ h9 X

ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
- C- g  Q+ X5 |" T9 ?4 q3 n, M4 c
9 a) w: x+ \5 ]$ }& u
( B( t* E- O1 M% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
  W: Z9 H5 m, v5 y& K! c% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
for i = 1:nbars
) p, K) V! k! @4 A6 S- z7 O6 A    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
. D8 F' Q3 t+ x1 pend
3 B' n$ v( J* W- D

) a& \: T+ `" E/ A( D& d( Qset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
' ^& V; |. u& r$ P( Y3 G2 |ylim([0 1.5])
" g: J7 r: s6 O$ Q9 I1 Wset(gca,'ytick',0:0.5:1.5)
8 O) j1 w9 S/ S1 u0 Sxylabel(gca,' ','Volume[Sv]')
% x0 m5 ~6 e% P/ k- E) G# t7 Ylegend('data1','data2','location','NorthEast')


! d7 U) w( u% K- a% [1 ~以上实例可以参考使用。


6 A- {/ d. F& s7 X4 lerrorbar的局部调整：
1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
Y = sin(X) + 1;
3 }1 p9 F( p# x6 g/ N8 pE = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
' z4 C, r' G" |9 e y_h = [Y(k) + E(k), Y(k) + E(k)];
, D5 R, _" e( u line(x, y_h);
+ ~0 y) D) G0 g; Q" e y_b = [Y(k) - E(k), Y(k) - E(k)];
4 L+ E" s" M* t5 w line(x, y_b);
end
5 d9 t6 H2 F: G6 ?' V% d& y参考：www.52ocean.cn

  v  h5 M2 \" P: X1 k: z. \
1 j2 S6 n/ i6 {3 Z6 ?1 Q                    
7 `9 E4 y! o1 s2 ~
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。