QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17636|回复: 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)
    1 x7 ?5 @0 ~% o) \. oalpha=0.05;/ H2 \2 U8 V- e4 A, h7 R( b
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。5 U" p/ Q% q( G: ~
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';8 S/ o7 f  i$ {/ q
    [mu,sigma]=normfit(A);
    7 m) Z& @6 y% ?$ ~6 y5 L2 Yp1=normcdf(A,mu,sigma);
    : L) s! Y' s/ H! L[H1,s1]=kstest(A,[A,p1],alpha)6 {. @5 u$ S0 `
    n=length(A);
    ( J( z0 `8 i. n7 s: l5 Iif H1==0) h3 N. w* r+ n
    disp('该数据源服从正态分布。')3 W0 ^/ e" I0 W! ?! X+ \5 k
    else2 Z! X4 d9 i% j6 _# N; I! i
    disp('该数据源不服从正态分布。') , P7 {2 G! m- k( X& N
    end
    & L0 ~5 Z- H  F" Kphat=gamfit(A,alpha);9 W; \* n$ L6 e) Y6 q
    p2=gamcdf(A,phat(1),phat(2));
    ) w6 E. D$ z" t7 |1 w$ A; K4 g[H2,s2]=kstest(A,[A,p2],alpha)
    * C9 K7 y3 G( p- S; |, J% ~" kif H2==0
    % `( g% b9 l: a3 d, c5 T* ]5 w9 fdisp('该数据源服从γ分布。')( W3 C/ k$ A4 ^& C+ U; u2 Q
    else
    9 K. S- D" i5 u, g( ]disp('该数据源不服从γ分布。')8 s9 D: E7 y4 l" o
    end
    " W/ `5 t7 l% D' ~2 C- _9 plamda=poissfit(A,alpha);
    " S. x( q+ d6 X( @( tp3=poisscdf(A,lamda);
    7 g/ U8 \7 `1 A2 `- e# P' X2 C[H3,s3]=kstest(A,[A,p3],alpha)+ I: C1 N* Y% g* V$ d) e
    if H3==0
    9 a5 C  R) Z/ t" ^% C" T; Adisp('该数据源服从泊松分布。') 3 j  N9 e; @% ~6 h5 H
    else, O* d0 ~" u) L* k) D
    disp('该数据源不服从泊松分布。')
    ' U8 K- a3 @8 tend
    ' ^$ Z4 I" c) f  k6 [mu=expfit(A,alpha);
      W% ~2 S. W8 j4 @/ ]p4=expcdf(A,mu);
    # }# V7 b4 L3 v1 @0 v6 H$ C& ?[H4,s4]=kstest(A,[A,p4],alpha)
    ) Z: I% f1 g% k2 L1 Q2 Hif H4==09 M- n, o! H. f$ l- Q0 N9 K
    disp('该数据源服从指数分布。') & }( K4 d$ m; E' N0 q
    else# X! K# T) |) n( Y) V
    disp('该数据源不服从指数分布。')   n! k& K$ q" |" u& t+ Y. K3 S; x
    end
    5 P2 }  k$ S  E* |3 z: T[phat, pci] = raylfit(A, alpha)
    * z5 n0 X" K  h/ Up5=raylcdf(A,phat);
    % J2 P6 j$ q: r  L7 V6 n7 ?[H5,s5]=kstest(A,[A,p5],alpha)
    - B6 H+ b' V2 ^) i# W# m4 Aif H5==0/ g6 x: r" i; J  r
    disp('该数据源服从rayleigh分布。')
      z0 l$ M5 s! M+ e3 g( }& s, qelse/ Y8 c, h9 a/ o. U' Q& x+ N6 C' t
    disp('该数据源不服从rayleigh分布。')  d& F% d  _( x6 ]/ ^
    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
      W. ^! Y( u! B1 l" k上面的是什么?

    9 E& Q+ x7 x7 j) Q1 P你是说整个么,是代码
    回复

    使用道具 举报

    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-4-17 21:25 , Processed in 2.485216 second(s), 72 queries .

    回顶部