QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17671|回复: 4
打印 上一主题 下一主题

对一组数据怎么进行分布拟合,求高手知道matlab程序

[复制链接]
字体大小: 正常 放大

5

主题

4

听众

163

积分

升级  31.5%

  • TA的每日心情
    无聊
    2014-3-17 18:33
  • 签到天数: 54 天

    [LV.5]常住居民I

    群组学术交流B

    跳转到指定楼层
    1#
    发表于 2012-8-20 16:34 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    function f=p_judge(A,alpha)" v( [% G; ]$ s7 h, V
    alpha=0.05;# n/ V; Q$ J8 n6 O2 ]
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    & W$ X4 _/ Y% M4 B4 w9 `4 kA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';/ \) v8 P  d/ y0 D( K# r1 o
    [mu,sigma]=normfit(A);
    $ T2 {6 N' _  Qp1=normcdf(A,mu,sigma);6 j4 b0 g& ~0 y) z
    [H1,s1]=kstest(A,[A,p1],alpha)9 p. Q" b- |$ ?) V( ?
    n=length(A);
    6 q9 B0 x" E$ [. |$ U9 x9 b' tif H1==0
    # s( o. j; C  g) j1 x  {disp('该数据源服从正态分布。')
    ) `# Z9 {  t/ j# ]else
    ) |; K4 r/ c2 S2 V1 xdisp('该数据源不服从正态分布。') ; g& t. y0 f* d
    end4 R* B) `2 z) ?$ d  [' ~( r
    phat=gamfit(A,alpha);
    " c5 G7 ^0 H1 n% e# D! J) g4 Rp2=gamcdf(A,phat(1),phat(2));
    2 o, y4 t! G* G4 V[H2,s2]=kstest(A,[A,p2],alpha)
    ( W" C4 E( B4 v; J; d( ^if H2==0
    + z+ M2 q, w- v: M; ]. t( Idisp('该数据源服从γ分布。')" W( l7 B1 O2 x! m5 O1 J1 @# F
    else! S# {8 X- L3 G! v% h/ q3 b
    disp('该数据源不服从γ分布。')6 ]6 F! H# r7 U
    end
    " J# g) P1 |3 Flamda=poissfit(A,alpha);2 l, J% h& h0 r% d2 D
    p3=poisscdf(A,lamda);
    . z( c1 ]% k1 G( \$ r/ S6 \[H3,s3]=kstest(A,[A,p3],alpha)
    1 x4 g2 D" @0 e0 V0 a! `if H3==0
      p( |. h6 j* @  Pdisp('该数据源服从泊松分布。') " x  T3 J1 M* A& K% p" L3 A$ k) `4 `
    else) _* ?$ I2 s7 ?. }. o3 R2 ?
    disp('该数据源不服从泊松分布。')
    9 T- u6 t1 w" {) G* t- wend
    5 V7 i6 b0 g) Z, b% W8 E" Smu=expfit(A,alpha);7 Z% F! W0 Y6 `7 w6 K
    p4=expcdf(A,mu);
    % W# p  k" }. ]% I3 U' }[H4,s4]=kstest(A,[A,p4],alpha)
    ) t8 n  S* r: B: o; V4 O- d  L1 m7 x% vif H4==0
    $ \8 S6 p% z9 I/ q! Xdisp('该数据源服从指数分布。') 8 e2 M4 _& e' |2 b
    else
    * w9 z% q5 _  m) n+ Y/ o: Q/ o( K/ Jdisp('该数据源不服从指数分布。') / U/ s+ S2 G8 A7 \3 [. X
    end
    0 a( L7 m& ]' t9 _# n5 `[phat, pci] = raylfit(A, alpha)0 J5 v+ m  g# J
    p5=raylcdf(A,phat);- W! M4 W& v+ U/ L' x8 w
    [H5,s5]=kstest(A,[A,p5],alpha)
    : e  l% B' n& Z" p% U5 ?) }if H5==0
      _5 f5 H5 K; g; N, l5 F4 Udisp('该数据源服从rayleigh分布。') 6 I/ g; R; s% V7 }& a4 s
    else0 R: r0 `: G' A9 p; R  u
    disp('该数据源不服从rayleigh分布。')0 a0 Z+ z4 z7 j( c
    end
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持1 反对反对0 微信微信

    0

    主题

    5

    听众

    57

    积分

    升级  54.74%

  • TA的每日心情
    难过
    2012-9-7 12:42
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    群组学术交流A

    回复

    使用道具 举报

    5

    主题

    4

    听众

    163

    积分

    升级  31.5%

  • TA的每日心情
    无聊
    2014-3-17 18:33
  • 签到天数: 54 天

    [LV.5]常住居民I

    群组学术交流B

    lillianxiong 发表于 2012-8-20 17:16
    7 T2 d0 j& `- g2 {上面的是什么?
    # r6 v) _/ U6 e! J
    你是说整个么,是代码
    回复

    使用道具 举报

    1

    主题

    5

    听众

    26

    积分

    升级  22.11%

  • TA的每日心情
    开心
    2012-8-25 08:50
  • 签到天数: 5 天

    [LV.2]偶尔看看I

    自我介绍
    软件工程系,即将参加数学建模比赛

    群组学术交流A

    回复

    使用道具 举报

    5

    主题

    7

    听众

    145

    积分

    升级  22.5%

  • TA的每日心情
    郁闷
    2013-9-16 12:49
  • 签到天数: 28 天

    [LV.4]偶尔看看III

    自我介绍
    西政一学生

    群组2013认证赛B题讨论群组

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-13 19:36 , Processed in 0.433298 second(s), 72 queries .

    回顶部