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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
* w7 P6 L/ d% }0 l% ?volume_mean=[0.73,0.45;
                        0.42,0.43;
7 B% g# ]! B5 \                        0.70,0.42];                        
  |4 T. H% N( j- L5 W( H$ v2 nvolume_std=[0.65,0.17;
& Y* i( S! B; T6 x1 d                     0.35,0.14;
4 \5 \' i$ O" F& O                     0.44,0.13];
%绘图                     
% E2 v# F  i0 k, r1 m! U( oclose all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色
$ L9 Z$ p* S7 |+ C/ \% ^( @

4 |7 `8 E) _6 R9 J4 {9 [ngroups = size(volume_mean,1);
5 R; r. x. M0 Q' o  Nnbars = size(volume_mean,2);
- R( ^( r6 ], A  H' l0 |groupwidth =min(0.8, nbars/(nbars+1.5));

& s9 ~3 y* D0 n) r# h
% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
% X7 B, K; i2 J- U%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
. j) r' M5 @: ?- {% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
for i = 1:nbars
6 J5 u+ {. t/ E3 X7 W    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);
end

4 C0 R4 b5 e7 U7 F8 d2 U. v5 V; N" y
5 }6 u) Q) b; }  y0 G9 {8 t9 ?set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
- d! V, F; Z7 o" W) Xylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
( D1 d* m+ m' A' t1 q1 [# glegend('data1','data2','location','NorthEast')
  q) R5 |2 M' p. e
! ]! P: D- C, u' Z& q0 Q/ E/ W
1 N+ b' G: V6 ?' F1 k以上实例可以参考使用。
" D; A% E3 V+ P7 L6 ?. s
6 {6 M3 v; C& a! s; a
7 f& J) @! T/ `& q0 R9 Kerrorbar的局部调整：
1.头部宽度调整
1 v/ W. ?$ P' O9 b3 b$ h4 m% Create errorbar
X = 0:pi/10:pi;
1 @9 R7 A! L5 ~1 K+ P: a( ]Y = sin(X) + 1;
E = std(Y) * ones(size(X));
( `% E( @; P" H" c9 cha = errorbar(X, Y, E);
/ G0 P) v2 v; K2 x- P6 a% Width of the top and bottom lines of errorbar
% @( ^# B2 A1 I& dxlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
* j, u# ~3 a% j% C0 [/ x, _ x = [X(k) - xlength, X(k) + xlength];
( o9 b" w$ U$ j3 z4 O+ z/ y y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
+ p& P, M: K/ R5 x) |) ~ line(x, y_b);
end
  u9 K  P7 S+ B: j5 h/ s参考：www.52ocean.cn
% N: m" u- N( N( _2 J3 P1 }

* [" o: g/ Z& `+ B  [9 M                    
( S2 H) v  h; u
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。