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

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

[复制链接]
(1)画竖直的errorbar) q3 E; @7 q3 P; I9 ^3 o! z; J. t
例子:$ t: R4 A" e; n0 n
close all;figure;
2 S, b4 ?9 P& G4 ^9 u# e& D    h=bar(volume_month_mean./1e6);       %--volume_month_mean 为12×2的矩阵* n  x& g! n2 D: S
    set(h,'BarWidth',0.9);
+ ~5 T, s$ x# M# p+ D3 E    hold on;7 H; t( _* c0 r( ]/ ?7 z
    set(h(1),'facecolor',[139 35 35]./255)    %--设置bar的颜色( ]& G4 W7 f) G, ]
    set(h(2),'facecolor','k')9 N; Y8 C/ t  p- l$ r/ e3 U
# H0 {  q" M$ I7 ]

, D1 ]8 |# p3 Y0 u; W  J    ngroups = size(volume_month_mean,1);4 T8 @( ~* F  m! a
    nbars = size(volume_month_mean,2);
5 e/ b6 `. r- \3 ^    groupwidth =min(0.8, nbars/(nbars+1.5));, i' v0 g% q  I0 v
0 p" O/ Q8 |! O, ]

! c3 i  w% K  b, N    hold on;" s8 R' x8 W/ ?1 o: f$ Q
    for i = 1:nbars              %--画errorbar
- t# y; r' k, t; `        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);6 {. \/ j; F4 q
        errorbar(x,volume_month_mean(:,i)/1e6,volume_month_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);& W& t2 t# T. ^6 h( U
    end
: O* a4 b8 v. y4 S    ; r+ o% x7 p0 A* p$ \9 S8 ~! ^
    for n=1:12
" r9 V- u: d2 U6 W: A5 a. ~3 ~        if n<=9
! ?+ {. g8 G2 M) u            time_lab{n,1}=['0' num2str(n)];
( B9 v/ `" K9 T6 j        else! b# ~: L9 ~+ i2 s
            time_lab{n,1}=num2str(n);& Q+ W/ D# S/ Z6 h
        end0 p5 ^+ V& ]: y; A( e
    end
# z3 V0 |; c* n" P" G5 c% l2 bset(gca,'XTickLabel',time_lab,'fontsize',14,'linewidth',2)
% E1 q7 \( H8 s. [    ylim([-0.5 2.5])
1 v% j, n9 r8 w  _# m4 U% |5 n    xlim([0 13])
! f% s- j2 C# e! a- |7 }8 z' ^    set(gca,'ytick',-0.5:0.5:2.5)
  j# ~# v, V3 |8 k   
8 b) N2 }5 n9 G% K- o8 V7 p7 d4 \8 _. [" {5 g8 u3 y

7 J5 V" I0 `6 k8 h. x; `1 L1 O(2)画水平errorbar2 B3 M7 r& D7 `$ W' Q6 M& x
例子:$ F  q8 |+ U+ b1 J; w% r8 `& N
close all;figure;3 H( |: H9 k0 q
    h=bar(volume_mean./1e6);           %--volume_mean是4*2的矩阵
8 Z' q+ Z+ l4 ]# f: U    set(h,'BarWidth',0.9);; V( e: s: p  r' I! d( I. t
    hold on;% N/ W/ Q4 W( T6 p$ \: f/ ]
    set(h(1),'facecolor',[139 35 35]./255)7 A+ q* P8 ~6 \9 n) i, @
    set(h(2),'facecolor','k')- Y3 M% m8 }0 I

) E8 n- R9 l+ a3 N1 \" w% H  j0 ~$ u, n6 P, v
    ngroups = size(volume_mean,1);+ X8 [; X, D: p/ t2 N
    nbars = size(volume_mean,2);) P5 ^9 r4 h; ?' M5 K
    groupwidth =min(0.8, nbars/(nbars+1.5));; [- ^5 J4 Q( T' x% g* N# A3 v0 k

: P6 ~) Y. r& D2 Z1 P: ^
6 q$ q% m& `0 K5 [# h; P$ l    hold on;
: G  |9 z  i* ?' N  J' w" C    for i = 1:nbars
& T4 l; S5 y3 R0 |8 S2 U        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);, t7 ~4 T- {9 R* Y. {
        errorbar(x,volume_mean(:,i)/1e6,volume_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);
8 }$ q) e, a7 e& G0 G- h+ m7 K$ w    end( m  o' r# c9 G" @. U
6 ~$ x6 Y  w* O* S6 f8 o) L# w7 B

6 M/ x' Q4 e) H    view(-90,90)4 L8 Y! p- j& y; Y# T
    set(gca,'xdir','reverse','ydir','reverse','XTickLabel',{'Winter','Spring','Summer','Autumn'},'fontsize',14,'linewidth',2)
7 x6 N$ y" W0 q    ylim([-0.5 2.5])) r8 x; w0 C8 C: L) D
    set(gca,'ytick',-0.5:0.5:2.5)" o6 C9 t. I9 E$ n/ w- e5 l. a( i( w# m
    ye_xylabel(gca,' ','Volume[Sv]')% Z: H0 s" i. b7 b: w
    legend('A','B','location','NorthEast')
8 {. _! `$ a2 l4 @( w, H% l) T   2 n* n- S/ N  F
; N; o, c" a; r

- U: h( A! V; C                    
* h  U6 `  |; D* O) H1 Q- S- d+ r; E0 A
                                        转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。
回复

举报 使用道具

相关帖子

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