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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
6 J( L3 c& [) l实例：
0 U; _3 B+ n+ y% s  }' c/ E% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
7 ^4 @# j: Y: u5 a( lvolume_std=[0.65,0.17;
7 b% o1 P3 a1 @6 t3 @8 c6 h4 u/ g4 y                     0.35,0.14;
+ d, X8 u: r9 f, C                     0.44,0.13];
%绘图                     
8 ]1 P5 n+ `/ z: E1 eclose all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
" W" ^! o3 c' S/ o* y5 v, C/ Vset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
0 ~4 `) T$ P# S+ O* ?3 q+ zset(h(2),'facecolor','k')        % 第二列数据视图颜色


ngroups = size(volume_mean,1);
% R7 @. j5 w! z; H+ {( W; I& a2 Ynbars = size(volume_mean,2);
/ t# q8 Z2 O- y3 e% ]! hgroupwidth =min(0.8, nbars/(nbars+1.5));
) M) Z8 d4 C/ J. V4 m
6 B7 W8 y) V& O* k2 `6 u! ~# m8 x" @
  ]0 G! v# _* k8 j& ?* J% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
: ^  f4 r% [8 Z8 @% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
" s0 s2 j; k3 J/ w; \( I$ @/ d. O% Qhold on;
* U! l$ f' f1 {0 k; f4 t! }( Mfor i = 1:nbars
    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);
( |5 {4 l& h' v7 Kend
$ Z5 W$ ^- b% L2 ?

set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
! \3 K5 @& g0 b) _& dxylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
: G8 H% f6 ^; m% s/ l

以上实例可以参考使用。
. j% f9 X) J- p2 w+ l9 D1 c& U' q) p: L

errorbar的局部调整：
1.头部宽度调整
3 g4 O- m$ I, p$ E  `( S. W% Create errorbar
X = 0:pi/10:pi;
  C& }+ ?# ^# t2 ?% UY = sin(X) + 1;
$ m7 v/ v: b$ _E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
+ a# L% ]4 {( Q" A( j+ V, L  ?2 z1 i% Make horizontal lines with 'line'
; c4 C% ~) Y3 d- b4 v7 {for k = 1:length(X)
- n4 ?2 P% }( c1 S8 h x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
- k- `, d4 g# L, K- D/ b y_b = [Y(k) - E(k), Y(k) - E(k)];
7 B( W7 ^& I* t1 D line(x, y_b);
: W+ Z- [" T1 g( m8 h9 Zend
参考：www.52ocean.cn

; L! D1 I, W/ v: d) b
                    
# f- z; y7 h7 C: x
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。