QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17766|回复: 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)
    ) h* w; W6 _& n0 {" f+ _" kalpha=0.05;. P6 b- Q; W. M  W. m
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    1 C2 L  T0 D4 `+ [7 ]) jA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';1 y; H3 I" ]. E7 z$ u( O
    [mu,sigma]=normfit(A);  c( H( }! P* _2 O* D9 p8 Z
    p1=normcdf(A,mu,sigma);
    / v# w9 v9 f# h8 J: H% n' C8 U[H1,s1]=kstest(A,[A,p1],alpha)
    1 Y6 j3 e$ N& f- o- Jn=length(A);, x8 U; _" d0 |2 ]. x' _
    if H1==0* Y; M! i5 D  T& ^& h3 [
    disp('该数据源服从正态分布。')
      {  w" s2 L+ yelse7 A/ q  r9 G) y
    disp('该数据源不服从正态分布。') 9 r: l/ U% W( M% h* K9 f
    end4 y1 B: a% L$ I  |
    phat=gamfit(A,alpha);
    $ ~( n% j- a: H% A# M) \p2=gamcdf(A,phat(1),phat(2));" h0 p- y$ a5 [! R2 f
    [H2,s2]=kstest(A,[A,p2],alpha)
    ; X2 Q* r2 J. c9 |8 V, ^8 B9 Nif H2==07 }& h6 f# A: K/ N9 h
    disp('该数据源服从γ分布。')
    1 c6 G4 i6 \- E( Aelse
    ' k* T; k  C4 r: S& T0 E) O1 Pdisp('该数据源不服从γ分布。')% ~% @) ]1 i! P  D5 F
    end( [9 n; _1 P: h! {" T& P
    lamda=poissfit(A,alpha);- Q! Y9 R( \" F  a$ I% Y
    p3=poisscdf(A,lamda);
    1 S& ^) N' ]$ H0 ^+ l[H3,s3]=kstest(A,[A,p3],alpha)8 I' Q5 w5 F$ \- z$ z" V
    if H3==0
    0 o0 u% q( |; S. fdisp('该数据源服从泊松分布。')
    ) \; ?' h. s% @3 Nelse0 z' X+ j: O4 o% C# x) B% [
    disp('该数据源不服从泊松分布。')
    2 q+ N  `+ t) {6 f. J$ ^1 q- V& _end
    ( i! y( ^9 o6 E7 Smu=expfit(A,alpha);
    5 g! M9 }8 _4 e1 W* \" ?& Hp4=expcdf(A,mu);
    6 {2 x) p& C. Z/ e[H4,s4]=kstest(A,[A,p4],alpha). Z) W- L( \' P) _
    if H4==0
    , ]3 l/ |  p" ~: Bdisp('该数据源服从指数分布。') 6 l& O" o3 X$ J0 r. c
    else; `& @' \" g+ G) {' z6 i" D; h
    disp('该数据源不服从指数分布。') 4 a1 b$ e6 ~) `) W( F/ F
    end+ t  y# I4 p) l) x7 b( o! D
    [phat, pci] = raylfit(A, alpha)  N9 L$ O* l9 c6 |" v
    p5=raylcdf(A,phat);; F" }2 l4 |1 c# K3 e5 m0 R( I
    [H5,s5]=kstest(A,[A,p5],alpha)
    7 y1 F. s& G/ _+ Iif H5==03 h0 K/ U2 l1 x$ p' K6 g9 O
    disp('该数据源服从rayleigh分布。')
    . c3 ^0 w2 V* g& p1 E' |else; x8 {4 I& d. }/ Q. V
    disp('该数据源不服从rayleigh分布。')# q: Z8 }* v$ ?8 c0 x
    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 ; L! P: O2 D0 y! _$ H
    上面的是什么?

    5 ]( Q5 P1 B, B: T$ v你是说整个么,是代码
    回复

    使用道具 举报

    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-24 20:44 , Processed in 0.647097 second(s), 71 queries .

    回顶部