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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
" N! r) N2 S, A; N8 m8 i                        0.70,0.42];                        
volume_std=[0.65,0.17;
                     0.35,0.14;
                     0.44,0.13];
; _' w/ N: A; I+ Q& u5 j( [$ k%绘图                     
close all;figure;
h=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
/ b4 }. X7 N. m$ g4 w7 lset(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
* U! p3 w, P; k) C5 Zset(h(2),'facecolor','k')        % 第二列数据视图颜色
$ n( j& l( V3 F4 o1 D( w
: d$ K6 Q% D* d
ngroups = size(volume_mean,1);
nbars = size(volume_mean,2);
5 P6 L* `& v- V' s; g/ U/ X( Ugroupwidth =min(0.8, nbars/(nbars+1.5));

0 ?# R0 U) V( O
% 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];
: ^/ O6 ?0 D& e- J5 zhold on;
for i = 1:nbars
) ]; Y1 N; T' W2 \1 ^    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

/ z- N2 B0 U3 Z7 \0 F8 V
" {- {6 M0 G( k! Nset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
xylabel(gca,' ','Volume[Sv]')
+ t6 K5 X' s$ o7 h, m) x" b+ }- Flegend('data1','data2','location','NorthEast')

# q% K2 w5 ^3 x! w: e
/ G3 [2 c: e! B- E& O: A+ |( P/ l以上实例可以参考使用。
* \* V- l. r  Y5 I  s: e; m
) i! v% n' f) i
7 k0 p3 f- d" N5 g% c8 @$ Rerrorbar的局部调整：
1.头部宽度调整
% Create errorbar
X = 0:pi/10:pi;
/ T7 |: h/ k# P5 ^' I5 fY = sin(X) + 1;
E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
2 g& A$ a: W9 K9 }xlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
* T. @4 W; e) l* o3 ] x = [X(k) - xlength, X(k) + xlength];
y_h = [Y(k) + E(k), Y(k) + E(k)];
2 w1 y$ @1 ^9 K line(x, y_h);
; e- _: D0 a3 k+ {5 Y  t y_b = [Y(k) - E(k), Y(k) - E(k)];
/ t3 C$ B6 n4 o: h line(x, y_b);
2 P. j9 x9 [- Q2 p% nend
参考：www.52ocean.cn
4 Y% y4 Z8 R/ _1 U% f1 J6 P- o0 R
$ A* ^; M4 }* P5 z
                    
) w/ B" ~9 m$ g, i
1 m- ?4 ^& n0 d. [8 }5 K" s                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。