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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
  L' x: M9 Y/ J) l0 Z5 M) @volume_mean=[0.73,0.45;
% N0 p9 \! l" x8 N$ j                        0.42,0.43;
' f% z! O2 l8 o) y* S; ]: X1 O6 L. h                        0.70,0.42];                        
2 ^/ [5 c/ W+ d0 ?, d+ Lvolume_std=[0.65,0.17;
                     0.35,0.14;
/ i& v+ ^% n7 F5 I! v: O7 K                     0.44,0.13];
%绘图                     
close all;figure;
2 x3 u  s- U) |7 X4 l& [1 `5 j: Ph=bar(volume_mean);
  _0 E0 M6 l3 J) T8 f1 H4 wset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
' \! }/ Z: V9 j3 `) b+ lset(h(2),'facecolor','k')        % 第二列数据视图颜色


$ Q* y' U7 D" H( f2 I& n5 x0 n; Engroups = size(volume_mean,1);
! p+ i9 o0 I6 l4 i& rnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
* H7 i2 }  d% ^2 L$ |' K+ F

/ ?# w7 v! u6 }) h8 U; n% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
  T$ E+ V, j4 e0 P' P5 m1 K+ j6 y$ Y%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
- O4 i2 l8 C0 [% w% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
  H8 v6 h/ ]& h, d+ Lfor i = 1:nbars
' B0 N5 ^2 ?- S3 D! q- J$ |    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
0 S5 _! t  q1 V: Z' [# R, k    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end
4 @6 q+ W; S! @8 X

8 n* k6 [$ r) U0 W% yset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
. d/ b" t' M4 {- d. f& W5 bxylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')

- s/ i4 W0 X+ B5 D  s6 ?
4 {# \7 Y% w, ]4 l1 w以上实例可以参考使用。
1 `! {. l0 F8 t) A) k

errorbar的局部调整：
1.头部宽度调整
% Create errorbar
" a% Q4 E2 O4 oX = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
7 p$ q6 x* V; Oha = errorbar(X, Y, E);
# h- \* o/ B9 y7 E$ P$ V' G9 M% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
% z( x1 T& C( Dfor k = 1:length(X)
; b5 I6 @+ m/ \+ \5 j+ L 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)];
line(x, y_b);
end
参考：www.52ocean.cn
* Q8 `% u, O4 z& D' G
; d, ?5 u- S$ r. C% c6 C
                    

                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。