QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17775|回复: 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)+ H$ Y5 B5 Q5 ^
    alpha=0.05;* Y! d: g$ A% x" }8 \0 `8 E
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。- d% e3 F- v( v3 |
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    9 I& K/ h5 D2 C[mu,sigma]=normfit(A);( o6 C+ q9 T$ b0 _1 P
    p1=normcdf(A,mu,sigma);
    6 ^7 J9 V$ Q; x; ]" ?. n[H1,s1]=kstest(A,[A,p1],alpha)" Z7 t5 f2 c5 {) x
    n=length(A);
    1 n8 a9 U" G. D$ \$ pif H1==0" B' x6 s0 n" D, C* V: |3 W
    disp('该数据源服从正态分布。')
      Y  w6 o; w: y6 j" X; Uelse
    . R! g: G3 L/ V: L8 P: y; _/ ^1 }disp('该数据源不服从正态分布。')
    + l, x/ H$ ~. send
    ( F0 z8 V. @( t  Z$ i/ J/ i9 bphat=gamfit(A,alpha);
    6 M% w+ |& l2 lp2=gamcdf(A,phat(1),phat(2));; |' j4 J3 H) f
    [H2,s2]=kstest(A,[A,p2],alpha)
    2 e5 {" I! Z$ f1 W: qif H2==0& ]/ {' t. w9 x! A2 t$ c/ U  I1 p
    disp('该数据源服从γ分布。'), R$ C6 i, p2 b, ]
    else
    & w7 n2 u) S  ~+ {disp('该数据源不服从γ分布。')& ?9 z+ s" Y$ F8 R4 b+ u! K+ N6 G
    end2 h4 d. F, V! i7 n7 `/ m
    lamda=poissfit(A,alpha);
    2 s: `; l. o& i! K1 q/ Tp3=poisscdf(A,lamda);
    & B. W, D6 X2 l0 p1 Y[H3,s3]=kstest(A,[A,p3],alpha)) @" y7 R9 `- |  @
    if H3==0
    9 Q0 u! z5 W- [1 f9 Pdisp('该数据源服从泊松分布。')
    9 x: `3 ?  ~7 U2 m& Felse
    5 H7 c5 w; X1 fdisp('该数据源不服从泊松分布。') ) B) `7 k* c% ~/ j
    end4 ?+ X, h2 D9 [+ x" C4 E  g8 q
    mu=expfit(A,alpha);
    + k* R& ~5 H  `) j6 m6 D) Q1 Gp4=expcdf(A,mu);$ V* c4 W/ q# z' w5 L$ x$ c
    [H4,s4]=kstest(A,[A,p4],alpha)0 J  ]/ `  w- K% J+ m7 q
    if H4==0: w2 b! G2 O5 p% u# H7 m
    disp('该数据源服从指数分布。')
    ! L) }8 t5 S5 F: t, l  Y' Z" G" ^  |6 Nelse4 [: u$ M( [+ T8 P" l
    disp('该数据源不服从指数分布。') 6 r2 Z0 \" Y; e7 r" M
    end7 F) K- d8 h3 h7 Y0 I
    [phat, pci] = raylfit(A, alpha)! N) y1 t, K8 Z+ t+ c. e, f8 B! `
    p5=raylcdf(A,phat);
    " f& ]: Z+ [. E* `# e" ~- [4 ?' n[H5,s5]=kstest(A,[A,p5],alpha)
    : s3 S, J) V9 b; R* b) L! ]if H5==09 A% u  V; J8 X" j( j* n
    disp('该数据源服从rayleigh分布。') $ Q) v! K; A7 O( d' t
    else) s& H8 m8 X7 N- n1 L# L
    disp('该数据源不服从rayleigh分布。')1 N. F1 F9 }7 d" l6 r$ P1 o
    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
    + A- x  @% G; W2 w上面的是什么?

    * M) \, G. P, m2 p5 {9 j& z你是说整个么,是代码
    回复

    使用道具 举报

    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-9-27 16:06 , Processed in 2.325730 second(s), 72 queries .

    回顶部