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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
7 t# v& \! }$ w" |+ _volume_mean=[0.73,0.45;
" R, X2 N6 T& x1 J. e+ E                        0.42,0.43;
                        0.70,0.42];                        
( c5 h+ r+ }% f! [+ hvolume_std=[0.65,0.17;
1 y2 m9 J$ G  @( o, Y4 F1 D                     0.35,0.14;
                     0.44,0.13];
%绘图                     
close all;figure;
h=bar(volume_mean);
% i  m5 x5 w! h# ^. Wset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色
3 L3 h9 y. o8 S. j! U
3 J  S  z9 _% j& x* ~4 i
ngroups = size(volume_mean,1);
. J& \0 @3 N# A9 o" _) \nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));

) h* _) m7 {# e
% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
- H4 c8 @1 J6 e! O' V%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
( S" v6 x! c/ i% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
+ H1 \  Q* |' Chold on;
for i = 1:nbars
  ]7 q) x% x2 r+ l3 d    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);
6 L' x; Y: Y% i6 bend

: S; A' x" x& T! @0 ]- F2 n
% \/ X  u. J$ i* r2 E, e; Oset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
, R4 I* e, r  e  T5 S; Qxylabel(gca,' ','Volume[Sv]')
' Y3 F7 s, Y& c- l# L2 @8 ^legend('data1','data2','location','NorthEast')

( E: k) ^( L* U1 s  }: n+ Y: T
以上实例可以参考使用。
2 S% a- h" L! U8 U& K
+ n' H7 j! d4 G' f3 S) I3 ?4 L( q
) ~9 L0 y/ U0 M0 berrorbar的局部调整：
! L& n* _% U) E1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
7 g& |+ @* X4 _, LY = sin(X) + 1;
% g' y5 o; x% M% g8 cE = std(Y) * ones(size(X));
& I: {0 m9 j) s0 I* eha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
( Q; `, n. I" G1 t% Make horizontal lines with 'line'
for k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
. m7 P! f. i- X y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
; g+ E9 ~$ U: R$ y8 g y_b = [Y(k) - E(k), Y(k) - E(k)];
5 H2 M1 F" }/ u% \ line(x, y_b);
end
$ _( \3 z  k" D7 {" G参考：www.52ocean.cn

: V9 C" f" l% p  L2 B5 C
1 G1 A8 z0 ^) d% j" A                    
( a! ~% q, T0 m7 E* }- k
# f  F$ c# B4 y. E/ ]                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。