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

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

[复制链接]
(1)画竖直的errorbar1 [5 f7 |$ B0 u; \
例子:2 n+ S  g5 i  R
close all;figure;
7 D: W! b1 ^- N( U    h=bar(volume_month_mean./1e6);       %--volume_month_mean 为12×2的矩阵+ u& v" {8 D1 U8 N
    set(h,'BarWidth',0.9);
! C$ }& ~9 P$ e- P    hold on;
% Q2 _) H) Y" Y# K    set(h(1),'facecolor',[139 35 35]./255)    %--设置bar的颜色
; C) m( s7 Z8 {: k) ?: s    set(h(2),'facecolor','k')# r4 J  [/ C" `8 M3 P; C" y! ?

6 k; h/ i7 ~0 M7 y- P; w
2 c; V6 |0 E, Q4 ?    ngroups = size(volume_month_mean,1);& A5 v. Z' z/ W
    nbars = size(volume_month_mean,2);
  o0 _: |' L# d& Z/ q! L" B+ ]    groupwidth =min(0.8, nbars/(nbars+1.5));$ V  [9 X9 h) L2 D& E

- D; I: \9 _: u3 ]. G, O1 c0 ]; I  ^# v
    hold on;! y) y% x. `6 S4 M
    for i = 1:nbars              %--画errorbar
6 V# M' s2 T' N- a! F* C: H        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);0 k) Q1 x/ `6 w8 ~- J$ @
        errorbar(x,volume_month_mean(:,i)/1e6,volume_month_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);
8 D: w( _3 g) a* u' F( }    end+ E* n: R7 C* U+ L: t
    3 u$ @. C6 X2 O6 u3 H
    for n=1:121 j( I- Q5 I2 H: M7 A6 J% p
        if n<=9
: x' I6 }1 K! s% }            time_lab{n,1}=['0' num2str(n)];: G5 f7 g7 X, p) H
        else
/ ^0 _4 l2 v% I  H. p2 y3 W* I            time_lab{n,1}=num2str(n);$ k, V7 ?7 M3 d7 c/ M; ]& g' t
        end- Y; d/ `/ ~7 \* |! `
    end5 |) d& \/ o7 [' t$ f) v
set(gca,'XTickLabel',time_lab,'fontsize',14,'linewidth',2)# z7 ?5 {+ t, v( \: A% L1 f* X5 m
    ylim([-0.5 2.5])/ }+ s: I( M& q4 D
    xlim([0 13])* O# `# q- s* K& ?4 F
    set(gca,'ytick',-0.5:0.5:2.5)
9 X' @9 ?& u4 S: H$ @   3 c3 }% a7 u4 e
1 }6 G! x7 w# }, L# n

( m, }  s2 `5 [8 a+ M(2)画水平errorbar
9 @" t( Y4 x) z5 F例子:4 F: t. w" {! V# A  K! x( K
close all;figure;9 c* |" V" J/ `4 i4 W+ t+ Q" q
    h=bar(volume_mean./1e6);           %--volume_mean是4*2的矩阵
' I! O1 \0 O# P! t. k    set(h,'BarWidth',0.9);/ i; i! t; Z# l& Y
    hold on;1 z' V) ?6 o- n) i5 a
    set(h(1),'facecolor',[139 35 35]./255)
# G4 o7 D8 f! |" [( b% @( k8 M    set(h(2),'facecolor','k')6 k( T( `0 G5 x$ H0 m
# \/ x6 W9 f; P. i* i' q

! l: o4 L" K/ d! H% p5 i1 J6 ?  w# \    ngroups = size(volume_mean,1);  w3 V: \9 l. Y
    nbars = size(volume_mean,2);9 x5 a5 }  e2 O: K/ N, y; N
    groupwidth =min(0.8, nbars/(nbars+1.5));1 D! d8 I: h/ ]' X$ s9 D6 r. ^
5 J/ {6 q0 b1 k
2 l. f" d+ s# `0 @) x' [) p' H
    hold on;
' ^% J8 L# Q0 Y6 V, j    for i = 1:nbars( e$ w% d  F9 @$ E* l
        x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
; m: p/ K1 X! V        errorbar(x,volume_mean(:,i)/1e6,volume_std(:,i)/1e6,'o','color',[.5 .5 .5],'linewidth',2);6 v8 e) \6 s6 ^, X& {
    end
; R2 e  X! P  @$ E$ H
% P) {+ {" _& a. m+ T+ P1 o  ^
# `0 S# _2 ~0 S8 u4 [    view(-90,90)
' A* {0 _0 O) U5 `, U    set(gca,'xdir','reverse','ydir','reverse','XTickLabel',{'Winter','Spring','Summer','Autumn'},'fontsize',14,'linewidth',2), c/ _3 O% Q* i8 J
    ylim([-0.5 2.5])# n7 s6 k/ U/ L( ^. h
    set(gca,'ytick',-0.5:0.5:2.5)+ G" J1 L* ]' Z- {2 `& ]
    ye_xylabel(gca,' ','Volume[Sv]')
1 R3 o  h: B) ~3 W+ e/ ~    legend('A','B','location','NorthEast')
3 {, l8 V2 C: z   4 b7 f3 x* l2 m5 G3 J
, g. b; W7 I0 j* g

! C  e  E& ^8 I& V9 O; m                    
5 u0 x# d  w, O( d, |- T# \: i2 s  Z4 l
                                        转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。
回复

举报 使用道具

相关帖子

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