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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% f( \! E6 i; ]2 H. t) j5 f8 }6 V% 数据
9 Q! Y" {" E' Z& _volume_mean=[0.73,0.45;
! X$ L5 B1 D: n                        0.42,0.43;
5 ~  y$ G5 x( J( L3 U                        0.70,0.42];                        
' c# K5 Q" o$ }1 D* j6 N7 {5 Mvolume_std=[0.65,0.17;
                     0.35,0.14;
" Z  w$ r3 L/ J! j                     0.44,0.13];
( _/ I) e  t/ s' ~9 M2 N%绘图                     
! g, @4 g8 B! U) D8 E2 Gclose all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
3 L4 y5 S9 v, H) Lset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
& Y$ R' h  Z# h  `6 q+ ~6 jset(h(2),'facecolor','k')        % 第二列数据视图颜色


ngroups = size(volume_mean,1);
( D+ `3 P! `7 Znbars = size(volume_mean,2);
7 L( G$ f% T% S! a# }! ?groupwidth =min(0.8, nbars/(nbars+1.5));


: a; f8 h, v% G1 r8 P  ~$ u0 k6 V: D% 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
' ]! V- r  C5 ]    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
8 K# n& A& l) N" N    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)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
/ g, ]6 N$ f8 b. Rlegend('data1','data2','location','NorthEast')


以上实例可以参考使用。

( `( _1 a5 A6 I+ L" e5 V+ w2 F
errorbar的局部调整：
1.头部宽度调整
+ v% H% q& k$ N5 x2 ?. o- s0 p& A) l% Create errorbar
1 h* G0 E% Z8 D$ J2 X* }+ |X = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
9 D0 X4 B5 ^& n- z5 b% Width of the top and bottom lines of errorbar
xlength = 0.2;
' z7 Q4 @; d! |% f# s% Make horizontal lines with 'line'
for k = 1:length(X)
* O1 Z1 n0 ^1 u/ ^- s- G4 T 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


                    

) ?- r/ ~; D4 E/ ?  w                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。