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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
4 {) n( v6 j- ^$ F) \& m* K4 s% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
$ b, \$ n; q4 H" p' ?) Cvolume_std=[0.65,0.17;
                     0.35,0.14;
                     0.44,0.13];
7 E$ V6 _0 x/ F7 R. D+ @: d/ E0 w%绘图                     
8 T. f$ d5 \+ z5 v$ }# }close all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
" J9 r" v- F  Zhold on;
% m; I8 R8 |3 Tset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色
: V8 {- ~5 J9 ]/ ^: p5 F

ngroups = size(volume_mean,1);
, w8 R' i2 p+ f8 Z9 H! pnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
2 C' P) D% G+ e: k

% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
for i = 1:nbars
% c4 k1 I  _- e; P5 s1 H; R' L( 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
6 O1 `$ Y% t* ~" e4 g$ B

set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
* ~7 ^, C2 R* _1 ~/ Ilegend('data1','data2','location','NorthEast')

5 n9 K5 ]  n* p, d) ?# Q3 z$ w
以上实例可以参考使用。


errorbar的局部调整：
( P4 x9 M# ]5 F" k* \; a1.头部宽度调整
* B  S$ y1 {3 N, v1 v5 S4 ]% Create errorbar
/ T2 W: t. \) w8 S/ EX = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
$ M" }# @; _2 Q% g2 \& uha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
, d$ E4 e0 d; @! R% Make horizontal lines with 'line'
( w, R( @; m0 H4 H: {- r  Y$ O3 Lfor k = 1:length(X)
9 P) x# d: Y0 ?0 o! D, q x = [X(k) - xlength, X(k) + xlength];
4 s0 i4 m* Q. e) T1 ]) c  b y_h = [Y(k) + E(k), Y(k) + E(k)];
; ~! v* \6 ?. l* y line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
end
6 t& m6 ~# F$ [- B参考：www.52ocean.cn
# |2 _' i7 v$ F

& t' X- m; j5 }( q  s                    
$ q3 ?7 T9 x9 [) T" y7 O
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。