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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
# u/ C! ~; W  S: ^  u) n实例：
; h7 G: y. t, t0 G9 P' d' m$ ?3 j% 数据
volume_mean=[0.73,0.45;
% }4 f" E2 U# N# P; U                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
/ ]& B! j; u5 C8 i( Y/ A7 v$ G                     0.35,0.14;
" y' e! n% S0 r' W# {                     0.44,0.13];
1 B  ^. ~& t/ D9 |7 _/ S%绘图                     
/ y- H# T* a4 A+ Aclose all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
' e5 n  w# V7 \7 {4 Whold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
' y: R, g  @, U! U1 l- w" P  n8 B8 Mset(h(2),'facecolor','k')        % 第二列数据视图颜色
5 P# S/ M; B; o1 Q# W* ^4 [6 W0 a

/ K0 r  H6 `" _' R* R/ Bngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
, v0 w/ P( p0 p. P0 A
3 U  r2 f+ U) q4 C9 a
! H  D( L" C; Y/ d/ c9 j3 U. ?% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
+ F+ \( R& K6 b$ r) c! T6 y% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
; W  P1 }0 ]5 n8 o2 b4 qhold on;
for i = 1:nbars
) O1 p1 t/ {( R% o    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

( {& u4 n$ a9 A. P
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
7 K9 o& k7 r- }' e8 mxylabel(gca,' ','Volume[Sv]')
: A% @6 a* M" K' _; S" x" xlegend('data1','data2','location','NorthEast')
; q. G/ j& A9 r/ S/ R" F! B1 {2 H

以上实例可以参考使用。

  w2 j) {" i# [2 H# ~
errorbar的局部调整：
3 J) a3 N" e5 s+ ?/ P6 @& p* M4 o1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
5 I: V6 k6 c8 i# S  ~0 sY = sin(X) + 1;
, ^; {3 ]( q# ~: x2 `+ A5 LE = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
$ r; ?/ g( r8 V" \4 D% Width of the top and bottom lines of errorbar
. Z% P8 I. ?$ c7 v7 p$ u$ J5 Exlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
4 W4 x/ W2 e$ l" ^ line(x, y_b);
end
% ]0 n  g! e6 O' N- x参考：www.52ocean.cn

  b( R" v$ l) |! ^& J' z7 F/ [5 T
                    
/ Q$ z6 n0 V- l5 \+ }
+ u9 K6 E8 K, E3 ?                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。