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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
$ G# S" n2 P8 `2 T+ Z, Dvolume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
) g0 j  {% z& L9 ?8 N, @4 r6 ?                     0.35,0.14;
                     0.44,0.13];
) h2 O3 m& e; y9 ^6 g1 k( K%绘图                     
close all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
. i. V5 x( P6 g: M( Nhold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
( B* k/ B$ r$ Q* E( @+ @set(h(2),'facecolor','k')        % 第二列数据视图颜色

+ k- D( D) D5 e( i8 d
ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
2 z, W' r) [7 N3 ~, h. c* M' Bgroupwidth =min(0.8, nbars/(nbars+1.5));
! h3 N9 |- ~  }# n+ u- @

/ y, r# E9 K* e( l" b7 l( ?% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
9 v: A3 c8 n. A0 g+ a2 i' m%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];
hold on;
. a; R2 D% g+ Nfor i = 1:nbars
    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
& g4 ^$ p- e' b    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
8 {& j* R) c- B6 z/ Kend
6 B( [1 a) v1 r" {

+ w" u& u. }" \" V, fset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
& ~% f3 n/ F, \3 Vlegend('data1','data2','location','NorthEast')
7 p4 ?% \# K4 O2 E/ Z

% u) [: t1 u6 Y1 k! q以上实例可以参考使用。
& J- r# j% h2 u

2 v: q! z/ x  R* A1 ?, ^3 ^errorbar的局部调整：
1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
Y = sin(X) + 1;
. N( I  m+ _1 Z- W' c; ME = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% l$ E; J9 B  R7 `2 [  }% Width of the top and bottom lines of errorbar
  Z4 L2 B0 j5 Y' q! U3 Y7 z: cxlength = 0.2;
% Make horizontal lines with 'line'
+ }* E4 T; z  q9 N8 P! xfor k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
* M# _& l/ ?. i# D y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
# d& `8 y7 x2 _' C  M y_b = [Y(k) - E(k), Y(k) - E(k)];
$ V& g2 g1 k/ Y8 h line(x, y_b);
end
/ Y* ]6 f' Q( K9 }参考：www.52ocean.cn
2 v' F3 U5 h) Q
6 `+ V# R; [5 d
1 M5 y$ |* o' `- p) y' [3 t                    
! j) D9 v1 C9 [3 \/ F' Z) p7 r
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。