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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
5 g/ N8 s- ~. c# k' R7 z实例：
7 I0 l6 c: y1 @8 C0 u  V+ \% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
# z0 g( s0 r, _/ e  K6 R. I: }                        0.70,0.42];                        
9 @" F7 g" k' j9 k4 avolume_std=[0.65,0.17;
: k5 H& L( ^* W/ W$ _* ]                     0.35,0.14;
                     0.44,0.13];
( D% g) t9 u" `/ S%绘图                     
close all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
! [6 ]6 d1 ~. Mhold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
. H# \1 M& R: x& q+ x1 _set(h(2),'facecolor','k')        % 第二列数据视图颜色
' j* s6 C6 w& ^0 k6 @% P+ m) T

ngroups = size(volume_mean,1);
* @& F4 A( i. R) E% z. n7 fnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
5 q) ], S8 @: [1 e# l) j3 p2 b
' e. t: E* q; h! t5 M0 Z8 w
% 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
8 t0 k* ]6 A9 v6 V    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
' C# u6 T) v6 f: X9 v    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end


set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
$ \3 B4 S9 a1 q# R/ r, ]5 a% \ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
) _# K4 n5 _" h' k1 h8 o$ Q" Y* q

以上实例可以参考使用。


6 h9 ]7 C: ^8 @1 Y; M/ W- O' ~errorbar的局部调整：
. D' O' I  n% r/ `4 }1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
5 I+ F/ z% _& H/ \  W/ jY = sin(X) + 1;
: {. y- G; n5 |1 aE = std(Y) * ones(size(X));
  U! m' E" T( Q# Mha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
) v( ]! V5 u/ @1 hxlength = 0.2;
( O: ^3 ]( V$ z  x) g8 f% Make horizontal lines with 'line'
8 }6 i% q/ _/ \- Kfor k = 1:length(X)
2 F) n0 J+ Z; D1 @5 D 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
" U! ^' X, J# c' t
: f2 Z/ {5 e$ v5 S7 F
8 J/ D7 D2 N5 c+ [                    
+ U1 Q: ~! L3 G4 ^" M8 x; W
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。