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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
5 J( x; `2 R" }* h& p: o  J% 数据
volume_mean=[0.73,0.45;
$ N8 D6 s1 C- Z9 a                        0.42,0.43;
+ C1 a# ~' x% I                        0.70,0.42];                        
9 J- B# n3 \# S* P! T8 z- {; ?( ovolume_std=[0.65,0.17;
                     0.35,0.14;
                     0.44,0.13];
0 L) a2 N3 y8 M; w%绘图                     
3 T9 ~+ a* I) Y3 I9 M, mclose all;figure;
% k1 G0 H, L' C, b0 b: t7 Sh=bar(volume_mean);
9 l! v* }: W( ]3 n1 Bset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
8 Z9 w2 V3 P7 X! p! bset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
' N4 I; k; k7 Z6 p0 Fset(h(2),'facecolor','k')        % 第二列数据视图颜色
, _2 A& Z8 ]& o7 L0 _3 k! b
9 \  ^6 n2 h+ O% z" v2 r
3 B2 E& P) P4 R. Xngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
. e" f" Y) r* g8 K
4 m3 r* m. M4 g! H
% 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];
. D; M0 q* J/ S& }, U9 K) Whold on;
for i = 1:nbars
( B2 }1 g6 {1 W. @1 `7 m    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
. m- R$ a9 ]* e7 X9 a    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end
) i& M5 O2 G) {7 W
4 X; C8 w0 n; b1 Z% O4 P8 P1 y
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
" ]( q% C: E- A% i+ Wset(gca,'ytick',0:0.5:1.5)
6 X( E5 h0 j% j4 nxylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')


' a  g8 Z4 [, C0 s& m以上实例可以参考使用。


errorbar的局部调整：
8 V7 L: Y; z7 w0 S8 j; q8 K0 |; b1.头部宽度调整
9 r9 ]- Y; l9 \2 X- X; M2 _% Create errorbar
5 ^2 t7 o/ W, _( ^8 LX = 0:pi/10:pi;
7 q3 O! k2 D- ?" ^: x& I0 T4 HY = sin(X) + 1;
" {/ F  V& Y# Y/ C8 [" h4 ^' tE = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
8 R6 g$ A/ S; |( M% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
, e. r3 i6 D& x% k. \for k = 1:length(X)
1 I) I  w( W4 w% X! j0 m" j' b x = [X(k) - xlength, X(k) + xlength];
# u1 e/ Z7 ~( v$ [2 @  }0 i 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
, C  B8 T& G) W参考：www.52ocean.cn


                    
# ~: c+ r. ~/ J7 B6 M6 l4 n1 T
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。