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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
/ g& o9 n: q3 w. x  J实例：
% 数据
volume_mean=[0.73,0.45;
& g3 t* h- g( A$ d. K! S                        0.42,0.43;
                        0.70,0.42];                        
# C9 [( f4 b7 n0 ^$ S# d" G4 E% g% hvolume_std=[0.65,0.17;
) H+ s9 t2 ~5 Y0 X8 @# l4 \; G                     0.35,0.14;
0 y7 \5 y0 S9 i& [1 y% v                     0.44,0.13];
%绘图                     
close all;figure;
- j7 j+ y! h( W& Y( E1 Y6 Wh=bar(volume_mean);
8 M3 r9 m3 R3 c" m5 F& H3 s. rset(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色

" T1 z7 N9 [  Q) q
ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
. M$ T' K5 ^) [4 i1 {  Z$ h" }groupwidth =min(0.8, nbars/(nbars+1.5));
* X" L+ X' y0 ]7 w+ ^

' {; s8 q% D$ x( T1 @1 ?% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
2 \; l( M5 P" `2 F! |%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
7 R$ a4 x- [' S1 J) T% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
0 J( c0 r' r( chold 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);
end
7 |8 n9 F7 `7 v/ z1 J2 K
" h# k3 Y3 {, q  ~0 m2 K$ A* J
set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
, @5 H5 ~9 l8 x5 k3 n3 m2 hylim([0 1.5])
2 V+ s$ \5 J- Xset(gca,'ytick',0:0.5:1.5)
0 v) q& v3 S0 e8 Y" D  b; axylabel(gca,' ','Volume[Sv]')
5 p% O  S$ P/ H; |* y1 p  mlegend('data1','data2','location','NorthEast')

( V, @- F. M2 p  T$ {
  T. c: x' w% V  w) G以上实例可以参考使用。
! B9 h+ t; j3 F0 G

errorbar的局部调整：
1.头部宽度调整
- H, Z  k# l! F2 M' u6 X; M! ^% Create errorbar
X = 0:pi/10:pi;
Y = sin(X) + 1;
! `7 \% J  j+ j! f. }: L+ Y8 a' _E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
xlength = 0.2;
/ @7 C( g0 q* O" I" g9 w1 c% Make horizontal lines with 'line'
& k3 q8 E& x/ g4 nfor k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
! I: D+ p% C5 b% i3 d( } line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
end
4 {: P+ g0 f2 z  ]参考：www.52ocean.cn
* \4 U8 m) s/ }4 W

; H+ U/ {) a6 u# n7 i, t) O: q                    
6 P% }; x2 C( i; F& s- Z  H# e5 N
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。