QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17737|回复: 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)
    & o0 L9 v2 U+ S, l) r, z6 Ialpha=0.05;
    " \4 `! ]% D) v: P: C5 |$ C' [% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。, e! `2 W8 Q  q7 {2 L( ]
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';/ a: J. G' D/ N, e8 j% p
    [mu,sigma]=normfit(A);
    # N5 q. o( G$ F. |" A3 X( N4 N. w7 rp1=normcdf(A,mu,sigma);; V: @& |+ A* M/ P
    [H1,s1]=kstest(A,[A,p1],alpha)+ K) l% a0 S9 o6 c( N* ?2 _3 r
    n=length(A);0 z4 x( @, N: t- ]2 @
    if H1==01 `- J8 [3 Z0 b, E; i- W: T( v7 j: y
    disp('该数据源服从正态分布。')
    ; W' }( A% B' v7 D/ I: xelse
    6 L1 @* [  _9 S. ]disp('该数据源不服从正态分布。')   ~/ W1 W8 w, n9 _. m+ P7 s( |
    end- w! N1 @7 b4 @2 ^+ ?) @2 K7 J
    phat=gamfit(A,alpha);
    8 }9 O1 k4 K0 ?' d2 Y  up2=gamcdf(A,phat(1),phat(2));
    4 ?4 U( B$ g) K5 K4 b# r[H2,s2]=kstest(A,[A,p2],alpha)7 X# C, c% i1 t1 J" E; o* W4 ]+ k
    if H2==0: j- `4 Q6 X% Y; U( b
    disp('该数据源服从γ分布。')
    % P& v# T9 v# w1 R  K2 v) p- nelse3 x* |+ _4 Z9 M
    disp('该数据源不服从γ分布。')
    # J. w: E# ?0 ^- P3 M0 O' `end
    8 s4 ~- {+ R" V  |6 c5 dlamda=poissfit(A,alpha);
    % T5 J9 }! K# R' p0 w' Y5 ep3=poisscdf(A,lamda);! U' H6 Y% C9 t
    [H3,s3]=kstest(A,[A,p3],alpha)
    : q4 A0 V2 ~6 L. Lif H3==0% K6 j" H; ^) q6 x  o
    disp('该数据源服从泊松分布。')
    - D) T/ W% H- j  lelse
    , s: j$ n/ H3 L8 n5 V" ydisp('该数据源不服从泊松分布。') ) P- Z1 Q/ x! [
    end3 D2 B+ D1 t! ~# @
    mu=expfit(A,alpha);! b4 k  I& x9 P6 h" E1 `5 ~
    p4=expcdf(A,mu);
    2 j& X  N! X/ Y8 ~" [% R[H4,s4]=kstest(A,[A,p4],alpha)& [) i+ Y* q# j0 W; u' Z. z2 c" i& z
    if H4==03 U& r* A6 {5 ?. c; T8 U- s
    disp('该数据源服从指数分布。')
    $ o9 w& `& S2 \1 ~2 D- _) h: eelse
    4 c" i* |* @  s- ydisp('该数据源不服从指数分布。')
    3 x( Y+ v# T1 o9 X! Q7 }+ X- mend
    2 P+ `' K6 N6 x( I[phat, pci] = raylfit(A, alpha)2 p  G7 C! c( L7 {8 z
    p5=raylcdf(A,phat);
    1 b+ f" y- T: a[H5,s5]=kstest(A,[A,p5],alpha)
    # a4 D3 E0 |5 C. z; vif H5==0- K  f0 a& _/ Q. ]4 x! d, N4 i
    disp('该数据源服从rayleigh分布。') 7 p: _, r4 w/ @9 G/ l8 Y
    else/ j  n& t2 j& P# w( U9 x' X: [
    disp('该数据源不服从rayleigh分布。')
    9 ^' R' ?5 a: e4 H* |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
    * ?$ D  |# _8 E7 H1 w$ L" ]( ~上面的是什么?

    ! P; W6 `4 w8 Y! S% S你是说整个么,是代码
    回复

    使用道具 举报

    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-13 22:09 , Processed in 0.647055 second(s), 71 queries .

    回顶部