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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
' P; U; ?5 Z% Y! m# S* z: r  [实例：
% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
+ z" N# P, \  y% ~: {: B8 g                        0.70,0.42];                        
volume_std=[0.65,0.17;
                     0.35,0.14;
* k( _* G4 |: x% @- G                     0.44,0.13];
%绘图                     
$ E, W9 r+ ^8 o5 Z0 l1 yclose 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')        % 第二列数据视图颜色


" f& _5 @$ C1 xngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
( }% f9 P6 J  `% W+ G0 r

+ @3 G1 T* d2 y' _) F( h" M% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
) t' _, {6 _" t6 s" d%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
& _+ k; v0 d* t4 \) \6 S" ?% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
* _" {6 l" H# a6 \* uhold on;
for i = 1:nbars
/ e, z& F& [# D" v& m    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 M' r1 V$ z, b
+ r; R$ T* m, y9 F" ?
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
1 R" x0 A# B: i! u6 @# Iylim([0 1.5])
5 g2 n, X# }$ i* t- V# cset(gca,'ytick',0:0.5:1.5)
$ `- n5 c, m" Fxylabel(gca,' ','Volume[Sv]')
, P) e' N; \& Z3 [! `1 Elegend('data1','data2','location','NorthEast')


以上实例可以参考使用。

1 i& z" _4 C; E/ ]; S; `  Z: \
, b6 X3 q% H: W6 T. H. ^# [errorbar的局部调整：
1 o6 q9 f1 q- r7 }4 F1.头部宽度调整
/ d5 @! J1 @& {" X% Create errorbar
8 C5 E. e1 a& X3 lX = 0:pi/10:pi;
% m( q4 C* @' P5 a: }' a( H! SY = sin(X) + 1;
E = std(Y) * ones(size(X));
) U0 N. i4 X9 L3 @. uha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
, W$ V& ^& g1 b1 l' _% `5 `xlength = 0.2;
% Make horizontal lines with 'line'
) Q' Z! W" b: q$ mfor k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
; T; X/ C' ^' {: c: W. x y_h = [Y(k) + E(k), Y(k) + E(k)];
) a7 }9 Y9 w* _ line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
/ T9 M" e& l8 D) y8 ]+ { line(x, y_b);
end
参考：www.52ocean.cn
5 B4 q$ ?9 ?: T2 j/ M
& J) Z. E$ n% Z
1 y$ L4 B6 K5 k3 D                    
/ M+ Z3 Y3 |5 B' f0 g8 ^: }# W
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。