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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
! |; a; K& t# z实例：
% 数据
* t% q5 W- @9 j6 Y3 Z9 n" ~volume_mean=[0.73,0.45;
& [& k# D: F, M                        0.42,0.43;
$ R% o: Y% B( ^1 |                        0.70,0.42];                        
  x0 Y1 I. f3 |" Xvolume_std=[0.65,0.17;
                     0.35,0.14;
4 _( _/ K! M7 p6 e8 h% t                     0.44,0.13];
%绘图                     
+ ?/ G# |9 [- wclose all;figure;
h=bar(volume_mean);
' J. K! h- C. S$ J% g' r% A" ~/ Sset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
. S* A* j0 Y: T4 vset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
/ L" C* w& j$ Dset(h(2),'facecolor','k')        % 第二列数据视图颜色
# V$ W7 f4 ~0 f) f1 @1 s

. N9 _5 j/ O. p; K+ l0 Ungroups = size(volume_mean,1);
nbars = size(volume_mean,2);
) P' ]: ]6 k8 w4 ~& Bgroupwidth =min(0.8, nbars/(nbars+1.5));
+ Q2 g( U* Q0 x3 w# g) x

% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
; i  e1 J5 V3 Q- z( \- p3 q% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
8 M. p$ K0 b2 s3 l$ k( p+ @' o& ~hold on;
6 {* N  l7 @0 Xfor i = 1:nbars
    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
) b" t' f+ S# @: G    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
/ c' x3 S+ _$ r  C) c$ X8 G& t" ^end


set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
7 c; I$ h. y2 T" X7 u. H/ W* C3 ~ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
& G# ^' V- ]8 }! N0 m- T; p

! D# A- z6 w) i) ^* W' X  a以上实例可以参考使用。

. u- j' u0 M9 G2 c. J" _
errorbar的局部调整：
1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
0 Y7 L' s3 y3 y2 x( `" g: K2 R% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
. T, H! j" u: q9 v x = [X(k) - xlength, X(k) + xlength];
* _. n. P" G, b' N% x y_h = [Y(k) + E(k), Y(k) + E(k)];
2 S, {, ], \3 y; }* B, v line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
4 }% [$ w- u5 G5 \9 b line(x, y_b);
* ^9 p( B. _, E3 \8 send
: c  E2 ]$ P# J0 Q5 G5 a参考：www.52ocean.cn
- R4 K1 Y& v% s8 ~, X

( _6 m9 C0 R4 e: _                    
$ U$ b1 U3 {& p/ n6 n
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。