收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

【MATLAB】如何画带errorbar的直方图

[复制链接]
(1)画竖直的errorbar8 ?; E/ ~+ x3 B
例子:8 N# i5 q) ~# i
close all;figure;
- i# K4 n: _9 s    h=bar(volume_month_mean./1e6);       %--volume_month_mean 为12×2的矩阵: P2 s$ i) V. p& O7 q/ v
    set(h,'BarWidth',0.9);
6 d8 G# a! L+ L2 g0 h    hold on;
6 s! i/ Z6 F! e  p4 m; {    set(h(1),'facecolor',[139 35 35]./255)    %--设置bar的颜色
$ C/ }$ W8 J2 ~/ T: V) ]    set(h(2),'facecolor','k')1 P4 f  }) I, ], `
0 {6 [: o6 O, D& p( l: `

  @, t4 l) E0 J- |( S2 [+ o5 s    ngroups = size(volume_month_mean,1);: ]1 k4 N) _+ l) ^( p! v2 e
    nbars = size(volume_month_mean,2);
  Q  l: M- K& d% ^    groupwidth =min(0.8, nbars/(nbars+1.5));+ n+ g1 Z, Z; ~4 _/ P9 M" S
% W1 N3 k9 U" \3 U) |, D
* `% y; v1 Y" }- A, i1 |6 o9 V
    hold on;
* s- Z/ u! _: u3 @+ q    for i = 1:nbars              %--画errorbar# O) i8 C( u6 r; G: r
        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);+ u0 m6 c( Z( y! ]
        errorbar(x,volume_month_mean(:,i)/1e6,volume_month_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);! O" K# \+ H8 A" a  p
    end+ `+ D( x) m/ t7 R8 L% F, N5 M
    8 ]& B$ |- g( O
    for n=1:12
8 \* d7 p! W. a        if n<=9
' R! Z) A5 T% g            time_lab{n,1}=['0' num2str(n)];  @9 |! \# H; }0 T7 h" S( ?
        else+ o/ E& Z2 M1 k' G1 _! \3 R' M
            time_lab{n,1}=num2str(n);
% X, i; Z/ Q7 ?- l' h/ ?        end
/ D* Y8 \# @0 [# D    end
# Q5 Q/ R. |( H. B# D5 xset(gca,'XTickLabel',time_lab,'fontsize',14,'linewidth',2)
  o# Y2 y6 E3 e; s  _# T5 X9 W    ylim([-0.5 2.5])8 {  G6 `6 P: K' f5 c0 y7 m: m3 s7 F
    xlim([0 13])- f3 h0 `  c* m1 Z& w# k% _5 K
    set(gca,'ytick',-0.5:0.5:2.5)
- ]2 N* q: o$ b* X; Y! G) a   
7 X* S: U# I6 _- _0 I, r
* d7 J/ r0 v+ W) e/ S) d+ G5 Y5 u1 n( F/ U
(2)画水平errorbar
$ i4 h3 O9 }" n% H例子:3 ]' H, q. E  R& `! S* [! a2 E
close all;figure;
# C2 `) T- O+ ?2 O, s6 G/ T/ W    h=bar(volume_mean./1e6);           %--volume_mean是4*2的矩阵
& m* N1 U: N( [2 J    set(h,'BarWidth',0.9);
+ M8 c+ s& }; G' g    hold on;) _5 Q' G% u/ ?5 U5 c
    set(h(1),'facecolor',[139 35 35]./255)8 [# j# [9 M# n/ d% n' V$ p+ S
    set(h(2),'facecolor','k')
6 U8 E; W: R; N; ^* c
0 [' Y* T- Y4 W# @7 x% M! c' p5 r
    ngroups = size(volume_mean,1);) z0 M% j1 ?5 F: j
    nbars = size(volume_mean,2);
2 Z0 O4 r" e5 [3 X2 r9 b    groupwidth =min(0.8, nbars/(nbars+1.5));
0 z) n& a6 A3 T" ?! G# H& i( I; \) }, u/ j2 z$ f3 V! T

! q! H, o; J/ V# J/ U. G  \    hold on;, L( F: d4 @7 {
    for i = 1:nbars
% q, N. k- }- i        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);3 M* E- g7 e: x7 E& \3 h
        errorbar(x,volume_mean(:,i)/1e6,volume_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);
1 J3 T7 ?0 s  z- [/ D4 t    end
( f! D7 I( j( L! [( G5 r: @) X, W9 }4 M6 e! i

; O' o1 [7 c. O8 I4 R: @) N( P    view(-90,90)* P9 x5 S& D) q3 b
    set(gca,'xdir','reverse','ydir','reverse','XTickLabel',{'Winter','Spring','Summer','Autumn'},'fontsize',14,'linewidth',2). `, h& o. f9 U: ~+ d/ V
    ylim([-0.5 2.5])% O0 t% n3 [$ t1 a& `
    set(gca,'ytick',-0.5:0.5:2.5)$ z$ K2 @: g6 j9 z; {& g# w0 b  m
    ye_xylabel(gca,' ','Volume[Sv]')$ ]0 j! ^8 H! ?' K" A1 E( w; F
    legend('A','B','location','NorthEast')/ [$ N7 t1 Q. S9 b& ]0 x" ]
   * F" v0 Y& n. s
2 d+ O/ `! @  o2 k/ H4 [) d2 b
. B0 X: K+ `7 j2 U" |  P
                    
$ |4 h. M- _: b
% b: z' F  p+ e$ x" z9 Q' [1 a                                        转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
黄金品质
活跃在2021-12-8
快速回复 返回顶部 返回列表