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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
  o' Y; T7 z; ]! X; `. y) Qvolume_mean=[0.73,0.45;
                        0.42,0.43;
/ |8 j7 n1 E) i! Z6 e7 w9 q                        0.70,0.42];                        
volume_std=[0.65,0.17;
                     0.35,0.14;
' C1 f- }+ ~% e  O                     0.44,0.13];
/ s- ?0 Y& H' W/ ^* E  s9 }" [%绘图                     
. T; e& z  n/ _1 l' s0 q/ N! ?% g- ~close all;figure;
h=bar(volume_mean);
, x  Q) T/ p$ x4 e* `3 xset(h,'BarWidth',0.9);        % 柱状图的粗细
/ M/ M! l1 X: D3 w2 U  z! R- I) D( Shold on;
0 P# j; @# A$ f% ?' a& J! \set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色
4 U" m7 ^1 Y4 y2 Z! X! {
& a  S' T# y9 ?+ b- Y) S3 {
ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
6 k, q+ r4 m2 u/ @  W
6 c: W& w( X1 Z6 m3 v$ U  t1 k3 O
8 P) p/ F: ?' V' s$ _6 G' F; n% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
# C+ o8 P9 U: n% Z* F2 u+ w/ [%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
  L# R8 J# \( c( q2 `% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
for i = 1:nbars
+ R( R% e7 Y3 Z) b6 D3 M; A( C    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
  |+ D" z, a" M9 F9 w0 h7 z1 g8 Q: }

set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
$ w) c0 H, c# d- Rset(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')


以上实例可以参考使用。
- M- Q% R% G3 B; x

* _6 [8 q8 \" ?6 @1 k9 P, Ferrorbar的局部调整：
( n3 q4 G% ?% W1.头部宽度调整
% Create errorbar
: g* l+ y6 G% a4 @* u( V  c+ YX = 0:pi/10:pi;
Y = sin(X) + 1;
E = std(Y) * ones(size(X));
5 `& T/ w6 F8 Bha = 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)
2 _. S3 v4 o: p6 x( h x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
: i# I9 E) ?% F3 p7 s( F- N, K line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
. O1 {' Q# b; ^. M6 W line(x, y_b);
end
# P0 |. q% }5 W6 \  b参考：www.52ocean.cn


2 p  J+ j8 Q& f8 g0 @% Q2 b$ v6 f& s                    
2 ?( Q7 X+ {* k$ i- Z
7 Z* l! \( q% k. A: c  ^* X                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。