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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
1 H* t! I+ j- _% 数据
! {3 y# }+ M8 v2 `+ k7 k7 Jvolume_mean=[0.73,0.45;
                        0.42,0.43;
. B2 ]2 k" K: U6 @3 g                        0.70,0.42];                        
volume_std=[0.65,0.17;
                     0.35,0.14;
2 x1 {/ z' ^, b) w                     0.44,0.13];
& c. p7 f# `5 K7 H%绘图                     
# l) q. `! N% Uclose all;figure;
' _! V3 ?. t6 N; E: R+ S0 r; N& O& Gh=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
# l9 F4 U4 j7 Z% {) t, G6 }hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色


- C: ~) f+ p% sngroups = size(volume_mean,1);
, E: W2 N) ~- r2 Lnbars = size(volume_mean,2);
, c+ f" m% q+ s1 cgroupwidth =min(0.8, nbars/(nbars+1.5));
" l" A& r2 p' y  p) \
1 K; K1 e5 ^* o; i: `# T
% 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
    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);
' U' n# ~! p' j4 ^8 ^' Pend


; m( |$ W" u# G8 b3 `3 O0 s, {set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
, m' m+ s2 x- y, G( k4 wset(gca,'ytick',0:0.5:1.5)
/ I( y  `( ^  q6 Zxylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
2 W3 n5 i, B, _7 V# X

以上实例可以参考使用。
0 M7 T4 i6 X$ T' }4 L

! v# Q6 _) T* N* p1 Rerrorbar的局部调整：
1.头部宽度调整
' q5 B4 {/ a/ `) T: ~) \% Create errorbar
X = 0:pi/10:pi;
. x+ f3 ~* d+ F: U/ Z  ?: @3 @Y = sin(X) + 1;
1 K# f$ T  Y1 z: A1 S- |E = std(Y) * ones(size(X));
  V) O4 E* ^( e" Sha = errorbar(X, Y, E);
* Z- f3 Q9 _0 ^- z9 O2 Y# R5 O) c% Width of the top and bottom lines of errorbar
xlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
. a, X5 @- X& d0 S y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
4 P6 L6 @' e* Y: W- B) U3 H/ j- E line(x, y_b);
8 }# x9 v  O5 u+ jend
参考：www.52ocean.cn

- Z4 D1 d: q) \+ Q
                    

2 x& g9 @; |2 g5 ^3 E                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。