QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17823|回复: 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)
    : Q0 t+ i  A+ W: o  I/ W; nalpha=0.05;
    ) _8 Z6 S8 f" J+ E2 a4 e+ p- X6 |) E% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。5 g" ~5 q0 T1 T8 c
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';9 L" r7 E% D% P6 W" d4 O9 z
    [mu,sigma]=normfit(A);: n# E7 Y. E0 _# x9 R" v" _8 w
    p1=normcdf(A,mu,sigma);
    " ]9 J; y! l9 d7 q[H1,s1]=kstest(A,[A,p1],alpha)
    & \5 Q0 u; n3 @. Z' g! B+ Y3 rn=length(A);
    2 `0 m: d6 g4 W) a% R' ~' b, |if H1==0
    0 Q, o$ j3 D$ L2 P& ~2 t! K8 mdisp('该数据源服从正态分布。')
    $ C/ g1 m  \' nelse
    ( c* I; Q. l* I( t3 ?2 f% f  p$ _" adisp('该数据源不服从正态分布。')
    1 T& Q2 P8 n2 }2 N8 I( J" @end) s  H' {4 V* M8 y
    phat=gamfit(A,alpha);' `7 m0 J- B2 Q: P
    p2=gamcdf(A,phat(1),phat(2));
    " I. F0 Z& i2 r1 y[H2,s2]=kstest(A,[A,p2],alpha)" s1 F2 D% N$ p
    if H2==0
    7 {: M/ o# m8 P& a) h" v- E* [( Qdisp('该数据源服从γ分布。')3 f8 g' {$ p3 O1 D7 X
    else
    1 h4 X+ {% r+ X! a2 A- S' c% ]& M7 Ydisp('该数据源不服从γ分布。')4 R; v% \8 L/ w! M* o
    end
    $ a3 J. P! X' q) f9 P+ S/ Qlamda=poissfit(A,alpha);
    $ q, {3 P) u5 O$ _- z1 n  ]  qp3=poisscdf(A,lamda);* }- [4 r( t; p7 T+ e
    [H3,s3]=kstest(A,[A,p3],alpha)
    7 @3 ^, U& _% `% D" [if H3==0
    # n# A% T) s+ L, X3 _: u8 B! w% Ldisp('该数据源服从泊松分布。') 7 s' V! k) `$ t& C/ O
    else
    ( ?$ M9 Z! C& j5 qdisp('该数据源不服从泊松分布。') - ^8 M. O  O  m
    end
    ! N' |# L* E( t2 z0 I! U  Gmu=expfit(A,alpha);& {1 J3 _. I, y; q/ T
    p4=expcdf(A,mu);
    3 p' l" u& d3 A: n5 ?: \, s[H4,s4]=kstest(A,[A,p4],alpha)- `3 Q4 s2 t5 j1 y! g) s
    if H4==0
      j7 O  \9 G4 M& r' j  {9 udisp('该数据源服从指数分布。') ' l# g7 t# V6 |/ a/ l8 j4 t" g3 m
    else
    , s; q4 D( G! z8 Z, ]: Adisp('该数据源不服从指数分布。')
    9 P+ V% v  g; pend
    4 F- c3 n8 n5 f- B4 {0 [& S[phat, pci] = raylfit(A, alpha)5 U5 O6 E4 w$ k$ p" Y9 k; @. i/ S
    p5=raylcdf(A,phat);/ b' O4 g: h# q! o+ o( B) k% O* b3 F
    [H5,s5]=kstest(A,[A,p5],alpha)/ ~8 ?, Y3 t% i8 S
    if H5==0
    $ S- @! |$ `. C0 `0 Mdisp('该数据源服从rayleigh分布。') " t# C- ?% u! c
    else  h! l# K# i( f; p  D
    disp('该数据源不服从rayleigh分布。')
    ; f6 ]8 n4 w  O; y3 W+ H! Hend
    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
    + j# a; l1 B' ?上面的是什么?
    5 Z9 o1 c' |& R
    你是说整个么,是代码
    回复

    使用道具 举报

    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题讨论群组

    回复

    使用道具 举报

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

    fastpost qq
    收缩
    • 电话咨询

    • 04714969085

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-11-1 01:39 , Processed in 0.532756 second(s), 72 queries .

    回顶部