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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
+ V* ]& d1 ^& T$ O1 r实例：
% 数据
8 R% a$ q( {9 G/ s  cvolume_mean=[0.73,0.45;
" c, F1 i2 p, A, v+ i& V                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
2 r. s. }3 y4 L6 t                     0.35,0.14;
  G4 y, O7 k1 v                     0.44,0.13];
%绘图                     
close all;figure;
h=bar(volume_mean);
! u: G  i& y: @1 B  o9 |; fset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色

$ D3 }0 }% q+ ^" d; T
; W* N8 h2 P4 _8 F; v6 t/ P( wngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
/ ~9 z' W( G  k2 u

% B0 A' @* E/ ^% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
7 ^! }" p% B* M8 \; o$ u: P2 W* ~2 L% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
& f5 y. X# A& W4 E7 {+ M8 H: ehold on;
% N) I$ g7 u, u( F# T. u, Vfor 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);
% U0 |: Q& \, i2 N" Send


set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
% a- q1 a3 R) Vylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
# B9 `: z% _& C8 Y; [legend('data1','data2','location','NorthEast')


以上实例可以参考使用。
8 p$ s  E: e  V. m" H
9 e- k- W, r: L% C. q% S( t" x* N
errorbar的局部调整：
* }; O% V# }/ V) @5 w1.头部宽度调整
% Create errorbar
, d( Q/ Q" x0 o; ~: NX = 0:pi/10:pi;
: u5 a# ~6 h9 h. e+ Q. n# LY = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
4 j+ u' g& ^9 d x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
9 `8 h" M# g& ]8 a; M6 L line(x, y_h);
! j% ?) Z. K% p$ L0 S4 p y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
+ t0 M2 F( K3 b% Kend
参考：www.52ocean.cn
3 W1 v; M$ g! g6 d4 N
/ p/ y+ Y5 M( B- \4 G+ f* B
                    

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