QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17712|回复: 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), u2 U' g# }9 U9 k- f
    alpha=0.05;
    2 O( r  g  e' w* _7 q+ R' j1 R2 d" D9 R% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    2 I4 h9 X3 o  R4 k5 n1 QA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    : M! a% z1 ?% h[mu,sigma]=normfit(A);. c- v! u6 x, w  M+ Q  X% \! e0 w
    p1=normcdf(A,mu,sigma);, q) L- L7 X; W) H
    [H1,s1]=kstest(A,[A,p1],alpha)( x  W. O8 X; Y% z
    n=length(A);! `/ L% w4 P& m7 f$ W7 r5 u: X3 k
    if H1==0* j' `9 j- P3 x- ^# a( l
    disp('该数据源服从正态分布。')
    6 ~4 d- y' G2 \( aelse+ u) t, V2 G" E  W+ a% A; @* j
    disp('该数据源不服从正态分布。')
    9 Q" M4 T' x4 c2 X  i* u+ `end
    2 I" C1 T5 F2 l6 ~  h! f, iphat=gamfit(A,alpha);
    " l1 I  _6 k! L9 J$ H- ip2=gamcdf(A,phat(1),phat(2));& I. r2 e+ o  E4 B! f2 P' f
    [H2,s2]=kstest(A,[A,p2],alpha)0 M9 ^: y& ?1 Q
    if H2==0+ S1 j" Q. q5 @
    disp('该数据源服从γ分布。'); d' C4 I, ^" }4 {  r  H* `
    else/ J" j( ^; S3 n
    disp('该数据源不服从γ分布。')4 d$ }4 J! c" y" Z# m
    end1 c2 |0 F: S! E: k) T# R: Q, ?
    lamda=poissfit(A,alpha);
    " z1 ?, F7 U+ ?! Ip3=poisscdf(A,lamda);
    2 D1 K+ s9 `. Y- D/ w! g: z4 K[H3,s3]=kstest(A,[A,p3],alpha)
    $ x; A" E7 n7 i8 P- E& A$ Yif H3==0
    $ `; f6 u/ R8 w. I% m& H* p) L+ ldisp('该数据源服从泊松分布。')
    ( Y/ @' d! V- velse
    % c5 G; \4 y6 s( k6 Idisp('该数据源不服从泊松分布。')
    , w, X( B4 G4 @0 H* Q7 L7 ^; }end* l* y+ o# _2 S) H8 J
    mu=expfit(A,alpha);
    - C* t& u2 ]# ?- ~) mp4=expcdf(A,mu);
    / s: w( C6 S* o7 \; N; B6 Z[H4,s4]=kstest(A,[A,p4],alpha)
    / s) o  N" P! B5 \if H4==0  ~$ C+ {' l0 j/ e  E+ }
    disp('该数据源服从指数分布。') 7 x7 v* Z* e3 W
    else
    0 @' c: a2 I, n7 _disp('该数据源不服从指数分布。') 5 U& B2 y% p: ~
    end
    , {7 [. j9 t0 _3 |5 L[phat, pci] = raylfit(A, alpha)/ |' e0 r* d) r; v8 J& `
    p5=raylcdf(A,phat);7 S+ ~; a: }2 |& ~" l) C
    [H5,s5]=kstest(A,[A,p5],alpha)
    : A% ^  c  R  ^1 Jif H5==04 }2 k4 M+ S7 Z! b* G
    disp('该数据源服从rayleigh分布。') 3 d( g8 V) F  P2 e9 m
    else
    ! |0 p( c! C9 r" Mdisp('该数据源不服从rayleigh分布。')
    4 E, T/ \, t0 i* pend
    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 + e# }1 i6 P- x3 f
    上面的是什么?

    * g: \5 p4 k# C6 D你是说整个么,是代码
    回复

    使用道具 举报

    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-8-2 05:40 , Processed in 0.450758 second(s), 74 queries .

    回顶部