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

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
实例：
% 数据
volume_mean=[0.73,0.45;
                        0.42,0.43;
                        0.70,0.42];                        
% E% @8 ]6 [5 _$ f2 Q  hvolume_std=[0.65,0.17;
                     0.35,0.14;
) H) d7 {0 i* v  q) [, n                     0.44,0.13];
* c1 P" Z% E: X* d4 H/ P1 W%绘图                     
close all;figure;
4 k. f$ O7 ^5 e3 sh=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
" o3 ?, e# T/ O! u0 Fset(h(2),'facecolor','k')        % 第二列数据视图颜色
3 T! r( Q7 Y( c# r1 j2 ^8 X
# X& h  l. r% P. V$ Q
ngroups = size(volume_mean,1);
6 m3 s" M* D3 C* j3 Hnbars = size(volume_mean,2);
groupwidth =min(0.8, nbars/(nbars+1.5));

3 a) a$ }7 D8 P; Z% s
8 |: A3 G4 X: k# b9 U7 ?0 r% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
2 \3 `# ~7 ?& h) _- Q9 S$ ]5 s. Q% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
% c7 W6 v/ s. y' T+ Qhold on;
for i = 1:nbars
: c8 O* d  h2 j    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);
9 Y7 j1 o5 S1 Lend
" k+ D8 y! j" i* S5 e2 R3 ~

1 I# B( A( b& Iset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
ylim([0 1.5])
set(gca,'ytick',0:0.5:1.5)
/ ?; j/ }3 o5 J/ s, d; ?" [8 Exylabel(gca,' ','Volume[Sv]')
8 Z  F5 }: @8 s2 [legend('data1','data2','location','NorthEast')

- L  p1 }. L1 j. ?% a
4 a0 G9 G  |+ Y8 l/ n5 o" s# }以上实例可以参考使用。


errorbar的局部调整：
1.头部宽度调整
% Create errorbar
9 G7 D1 |  B0 t4 q6 x0 tX = 0:pi/10:pi;
# O1 j0 ^! l) i2 d" KY = sin(X) + 1;
; X2 T. j" W# u5 t+ n9 K: `E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
: R: y, P! `  e) c( W6 s% Width of the top and bottom lines of errorbar
: N# O- {1 k5 o. X1 k. U/ l3 \) zxlength = 0.2;
% Make horizontal lines with 'line'
for k = 1:length(X)
/ L- C5 ]) d' h* k x = [X(k) - xlength, X(k) + xlength];
" L) ^+ o2 x# z5 R  J' y4 a$ } y_h = [Y(k) + E(k), Y(k) + E(k)];
8 |1 F3 f. s, x- j: \  C line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
end
! h8 ]& J* ?2 O3 _2 v4 F( i参考：www.52ocean.cn
/ w% ], q* M2 c# H

                    
5 c, t# Y8 C+ K' t) ?
                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。