QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17718|回复: 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)$ g4 j" ^2 t! C% O
    alpha=0.05;
    : K3 {0 ?5 F* B( g% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    ' p8 ^. ]( Z; Q& g' \% L7 xA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    - _3 D. @; A+ U7 A. d3 w5 x) v[mu,sigma]=normfit(A);
    3 R3 I- p, V1 f9 sp1=normcdf(A,mu,sigma);( @: ^% D: W' N. e
    [H1,s1]=kstest(A,[A,p1],alpha)# ^6 k7 g  G8 V3 b  O4 G
    n=length(A);# D& W5 b0 F8 ?# D1 z
    if H1==0: q/ H6 W4 I3 V" _
    disp('该数据源服从正态分布。')
    * d8 F* \9 C2 s9 [9 f( lelse/ B; L8 J9 g+ V6 O
    disp('该数据源不服从正态分布。')
    / ^8 ?  K$ T1 z$ rend' r+ ?4 P, {- K2 q0 F
    phat=gamfit(A,alpha);
    0 M$ m  o' L1 x9 Hp2=gamcdf(A,phat(1),phat(2));1 C! ~) b9 l9 C' H8 ?8 H, @
    [H2,s2]=kstest(A,[A,p2],alpha)
    8 F* w* y$ n; M0 u( k; uif H2==07 J& D' u; u2 Q6 }7 ]
    disp('该数据源服从γ分布。')* d+ w6 D' t( N# @! e) ~
    else; _3 H. U8 i1 q, ~2 y: D
    disp('该数据源不服从γ分布。')0 B1 x( U  W+ c" B& [8 {) [, [/ J
    end
    6 v2 R! i+ b+ @, b- Nlamda=poissfit(A,alpha);- R% y) R1 c% ~1 m1 b
    p3=poisscdf(A,lamda);
    6 T3 O2 t) D+ V" N% f" D& E[H3,s3]=kstest(A,[A,p3],alpha)# n  {8 L  R) k
    if H3==03 N0 j7 j) W  L8 t  T8 Z& t
    disp('该数据源服从泊松分布。') 3 Y$ p: m$ i. R
    else
    / }' Q+ \/ |9 N: Jdisp('该数据源不服从泊松分布。') , i+ n6 ]/ q1 L/ x) i9 s( ~% a+ N) q& p
    end
    8 j# b. E& k+ imu=expfit(A,alpha);
    # k: z0 d& ^7 G" g5 h+ up4=expcdf(A,mu);
    8 f  z% K5 P* k2 Z[H4,s4]=kstest(A,[A,p4],alpha)
    5 Q2 J+ K6 U( Y# `  Uif H4==0
    - H/ G: F; P$ H8 o  t1 pdisp('该数据源服从指数分布。') / [0 P! S) u" m2 _% C4 C
    else
    4 _% E! ?! o) T; z2 n8 Q/ mdisp('该数据源不服从指数分布。') 0 i2 x1 o' M- {4 l: s
    end) o  k5 p1 y6 F
    [phat, pci] = raylfit(A, alpha): F: C& K+ Z! M+ d, l
    p5=raylcdf(A,phat);
    - P6 b3 W" \" `; U+ y+ C& }+ d[H5,s5]=kstest(A,[A,p5],alpha)
    1 K/ N+ \6 t5 L0 z8 }' Pif H5==0% y) d: s7 k, K( O9 _0 Q% c1 b, D
    disp('该数据源服从rayleigh分布。')
    # b$ E9 ]- P* n6 G2 P9 ^2 R- T9 Selse
    . N* N( G" z/ q) m; k4 Udisp('该数据源不服从rayleigh分布。')1 X" s( a( O3 Z" |. i  ]6 M$ g2 a! S* [
    end
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持1 反对反对0 微信微信

    5

    主题

    7

    听众

    145

    积分

    升级  22.5%

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

    [LV.4]偶尔看看III

    自我介绍
    西政一学生

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

    回复

    使用道具 举报

    1

    主题

    5

    听众

    26

    积分

    升级  22.11%

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

    [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 ! p3 ^: y2 }4 O9 |9 |% M
    上面的是什么?
    - ^$ w$ d7 O" S6 @% P* A
    你是说整个么,是代码
    回复

    使用道具 举报

    0

    主题

    5

    听众

    57

    积分

    升级  54.74%

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

    [LV.2]偶尔看看I

    群组学术交流A

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-8-5 02:32 , Processed in 0.495794 second(s), 73 queries .

    回顶部