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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
, Q1 L9 o/ p) ^3 Z3 P" C' N: |volume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
7 |; _/ O7 @3 I$ [* |9 Lvolume_std=[0.65,0.17;
" {) b4 a9 U  _3 G7 Y' [* U                     0.35,0.14;
, d8 D& w" G9 M( W7 |8 g                     0.44,0.13];
/ b/ B* F  l9 P1 E" A# O%绘图                     
! l8 z. z( X9 [0 S4 e+ lclose all;figure;
h=bar(volume_mean);
: q! B% f3 e8 L/ P! Nset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
, R! A0 Z) Y: h; I" E' }: _' B8 sset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
) H1 _- ]; q% H4 k# m, jset(h(2),'facecolor','k')        % 第二列数据视图颜色


ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
. [9 D2 k* \0 igroupwidth =min(0.8, nbars/(nbars+1.5));
' W1 l/ c2 J! K4 j! }! i
! i: e) N3 [( T
1 d/ N6 @$ r- H6 {' E5 Z2 k. v% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
4 V- v" d. ?# A0 ?+ u  M# z0 U3 q% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
. H! B! B/ D/ x% v0 n. A4 O% ahold on;
for i = 1:nbars
    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

' s, w9 |% K, I+ h: @$ }. e+ l+ Q+ _
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
! f* f3 S6 }5 M9 [8 B9 X* f5 ~6 qylim([0 1.5])
; R8 C; z3 P9 g& i6 x+ b; vset(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
# j. w) r0 Y9 Y* R/ x: Q7 ]8 h
5 e& H& ]  D) o3 z
9 z9 e0 |1 ^( D+ \! X5 r以上实例可以参考使用。
; P9 H3 W5 {4 T& @
% ~0 R+ r, A" i  z' u, Y5 @; ?
errorbar的局部调整：
  G" S0 D; ?. F7 W. z/ P1.头部宽度调整
% Create errorbar
. r# }7 ]/ o3 G5 ]" d3 WX = 0:pi/10:pi;
0 t& s/ o8 ^0 ~1 H! X3 U6 C  }Y = sin(X) + 1;
E = std(Y) * ones(size(X));
- d. C7 {0 J2 f1 M3 y. ^( Q% @ha = errorbar(X, Y, E);
/ p$ f# o% [2 ^  r% Width of the top and bottom lines of errorbar
xlength = 0.2;
- y/ z& i3 v5 W% Make horizontal lines with 'line'
! P1 z/ {- Y  R& I1 qfor k = 1:length(X)
  ?" q' x! h7 o- [' L" N3 ? x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
4 P2 H& l8 ]6 X9 ` y_b = [Y(k) - E(k), Y(k) - E(k)];
1 j, y: T% |9 ^5 { line(x, y_b);
; \- G' e, Y- L- x# f& Aend
, m3 e* {8 M' G参考：www.52ocean.cn
1 C. X& Z  \5 x+ v- n  O9 t

& T; s! r6 }5 J* ]* q                    
: M( n# p* l& d9 t4 k- w! E" o
; X5 H* p9 b. B" e3 k                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。