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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
1 C9 \/ {$ m% e1 V% Y! J7 A2 ^- {4 Q' B& Q实例：
% 数据
. v* G' a7 D3 K6 Qvolume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
volume_std=[0.65,0.17;
9 H- l2 p2 X6 C" Z) T# N                     0.35,0.14;
                     0.44,0.13];
%绘图                     
close all;figure;
h=bar(volume_mean);
7 M/ ]4 G8 Z6 i5 Q# C4 D% `6 k# q$ Wset(h,'BarWidth',0.9);        % 柱状图的粗细
( E2 ?' |5 c! f- [hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
+ b; J  d6 A* L$ d* a* gset(h(2),'facecolor','k')        % 第二列数据视图颜色


ngroups = size(volume_mean,1);
( l0 w0 t& P" jnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));
! f& h: k: S( ]/ t7 y0 s3 h' C
: v. Y% S9 g, s" D5 J1 d  F. I
% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
. a$ {2 S  T: Q6 t5 X( \' x3 @%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
' V  C. Y* E7 B! k  b' b% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
hold on;
for i = 1:nbars
: l/ O1 D0 h- ?8 C    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
: Q* T4 r( ~9 B  a# P' g4 J) @    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end
5 d. s/ \1 N" y

set(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
. A( e* @& o+ W! X1 E% ixylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
7 Z, U1 }$ p" L. D7 u
* D0 _! L9 `" X' A' S: b) O
以上实例可以参考使用。
% O* w+ m  H8 x. \! F* S

errorbar的局部调整：
, e2 e' M- e) F; T6 I$ ~1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
) v' a  Z/ C) E% H% X4 TY = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
% A& g% K9 \* o% O7 y6 axlength = 0.2;
% }$ F  A3 b6 a! j. j6 n$ {% Make horizontal lines with 'line'
( a1 r0 ]" [0 x; Z  c+ Wfor k = 1:length(X)
  w! @" P5 D+ @! Q% H5 Q, w x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
4 O: u, U0 k- d: D) W  E$ A4 C" { line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
" T7 O2 N" t( `1 i line(x, y_b);
end
参考：www.52ocean.cn
" J4 c- A7 \3 U. B! l' }" K. S/ q, Q
" U7 |8 o, |. u! A' l
) i& W3 B- e* H2 g. U! R. h                    
2 g3 S) {. B  g7 V! n' Y! c
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。