QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17736|回复: 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)5 @# ^, b) d  B: h  W! L' D
    alpha=0.05;
    ' T9 d( ?  F5 J( g' ~- A( _% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。9 u: c1 N; b+ O% `
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';7 V; x$ L: w$ j3 i% t4 ]% Y' \
    [mu,sigma]=normfit(A);" |0 p) e% [- C( n) h
    p1=normcdf(A,mu,sigma);
    * C* C& {- r! H0 }" _4 N[H1,s1]=kstest(A,[A,p1],alpha)4 z( k1 Q8 b. j: u5 E; m; t9 h
    n=length(A);) L3 V5 Z0 l) b2 \) {9 u
    if H1==0
    3 y0 o3 D' a, N" j+ l4 ?! c4 {disp('该数据源服从正态分布。')
    $ |! T7 A( A6 ]$ v1 qelse- r6 C  J  B, C* U: Q' y
    disp('该数据源不服从正态分布。') + Y; `5 g; `' e& k5 k
    end$ F. L# D" @0 N3 e
    phat=gamfit(A,alpha);# F* G% F3 P: t
    p2=gamcdf(A,phat(1),phat(2));
    - t+ J, f% b$ q" l6 e[H2,s2]=kstest(A,[A,p2],alpha)
    ) l5 T; W) f: w1 w) X2 oif H2==0
    7 m5 _8 ?: X8 Y- [$ ]  Cdisp('该数据源服从γ分布。')1 U, _6 Q5 e- h% Y+ z3 u; c
    else# B* H! f! @3 a$ N
    disp('该数据源不服从γ分布。'), t2 H2 J: F' T# l
    end0 J% Q  C) u1 `$ J+ S
    lamda=poissfit(A,alpha);
    & x$ a" w9 h+ }p3=poisscdf(A,lamda);2 y$ B  i2 I3 i, E
    [H3,s3]=kstest(A,[A,p3],alpha), g3 s  s# x+ p  @* v6 K8 O
    if H3==0
    6 k8 I4 x  w$ Z2 V7 M: X: Jdisp('该数据源服从泊松分布。')
    : \; w; C( b" i2 {# yelse
    : S- W# f2 x4 D: l$ \disp('该数据源不服从泊松分布。') 9 x: J+ e# ]% Q: W) ]
    end
    3 [' M- X5 N0 l, j% }9 h8 O7 Vmu=expfit(A,alpha);
    4 U0 w. \% h9 G( lp4=expcdf(A,mu);0 O' L' U& N! ]0 M* @
    [H4,s4]=kstest(A,[A,p4],alpha)
    3 d- F2 O) S4 o$ |if H4==0; m( d2 w9 p/ h
    disp('该数据源服从指数分布。')
    9 W4 p, G/ E: }. zelse$ k  D, G9 m/ G# l( F1 w6 Q
    disp('该数据源不服从指数分布。') 3 f; s4 @. t7 X2 }) O: M* M
    end$ }; C. f: K$ {) p( ^/ Q
    [phat, pci] = raylfit(A, alpha)
    ' N4 f/ z! }# q  sp5=raylcdf(A,phat);* Y' c0 U( A0 _0 p: Z
    [H5,s5]=kstest(A,[A,p5],alpha)
    ! P" ^% Z! @; }  Gif H5==0
    - R) H6 j% A' d" hdisp('该数据源服从rayleigh分布。')
    % _, w6 g3 I7 G& `1 [else& A* V9 h# }3 H2 g" M3 w3 N: \3 Z7 J
    disp('该数据源不服从rayleigh分布。')
    6 x. p9 U/ b: y9 w" ~; xend
    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
    1 M6 b& m( v- P. e4 G* n上面的是什么?

    & k1 X1 Z& w: 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-13 20:10 , Processed in 0.446220 second(s), 72 queries .

    回顶部