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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
2 @5 r/ q* `2 R( Z. |0 a                     0.35,0.14;
# E6 O( r3 k" Y2 N' r' M0 Y                     0.44,0.13];
$ x6 E/ h4 l6 L3 X%绘图                     
close all;figure;
3 L8 y& P1 E7 E+ J* wh=bar(volume_mean);
" w6 ~& Y: F, X- Bset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
5 x+ {+ I, y* c* t; G! aset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色

3 y; F" |5 |( n1 l
  C- y7 v4 ?0 I; ~7 W" T) R( b9 yngroups = size(volume_mean,1);
8 g( W0 ~8 w! T/ ?. E  R* |8 t$ jnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
( @" _* O  n6 w5 ]- w; O

4 Q- P9 v9 `2 U; H5 z% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
3 ^' K6 ^/ E" {3 c( `%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];
9 [% J- ?, }: x8 t, Khold on;
for i = 1:nbars
4 R2 D, G) m& a/ A6 E    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);
0 M7 l- J2 J: d, oend

# `- p3 f- A- \8 k+ v
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
& @6 x6 C8 d  h! gylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
9 f. W/ v! y% g0 B0 nxylabel(gca,' ','Volume[Sv]')
- f) E) r& ]! I7 J/ K% E: Wlegend('data1','data2','location','NorthEast')


' s  @7 Q% M9 G* x0 V$ g1 b: `以上实例可以参考使用。
2 ]: [: \% N- E: d

3 G* y0 v5 A8 j8 j# i" W. Ierrorbar的局部调整：
& I" {# @/ @4 S0 f6 ]1 f1 F- M1.头部宽度调整
- {8 e; F* W" G1 V/ k% Create errorbar
X = 0:pi/10:pi;
6 D# K7 E$ B, d- [, |Y = sin(X) + 1;
9 r3 F+ v% W$ s  }5 n0 H" u/ PE = std(Y) * ones(size(X));
: b/ r$ y& X9 f# h. Q0 a3 Tha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
2 @8 e6 p4 `: @2 P0 x, u. ]$ R% Make horizontal lines with 'line'
for k = 1:length(X)
- y0 s" L, k9 A' H4 T8 J x = [X(k) - xlength, X(k) + xlength];
% W' R+ _. z  _( K6 E# g y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
, J0 R  {1 C4 X2 gend
参考：www.52ocean.cn
- d/ E6 ]" ]  s3 f4 Y$ i

- M- e5 u: K5 {* b4 N                    

                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。