QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17716|回复: 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)$ r& l) j/ G. g' ]7 U3 J, ^
    alpha=0.05;; ^: W; K$ L8 k5 @! ]% f0 _
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    - \$ x- u* j4 Z4 x% Q4 a+ {, YA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    / U8 b' y3 \  i% c: _[mu,sigma]=normfit(A);
    ; C8 S" h8 H/ C- D: ]6 |  Jp1=normcdf(A,mu,sigma);
    . L7 t5 g' E- T% |[H1,s1]=kstest(A,[A,p1],alpha)) d' [  p( @1 ]  {# g' C, z. d; i
    n=length(A);
    ' f7 D; @9 @& Dif H1==0
    ' N* X9 E1 E4 W9 qdisp('该数据源服从正态分布。')
    % m! e6 t5 p, N) }else
    0 ^, b8 m: C8 l0 ~/ v: Y2 e' `disp('该数据源不服从正态分布。') # M" H2 e0 C" A3 b9 C0 e" ?8 U
    end: L8 X0 z7 R6 S
    phat=gamfit(A,alpha);
    $ d5 ?4 A  V5 x, ?p2=gamcdf(A,phat(1),phat(2));' O& [1 `4 |6 A: L, c' P
    [H2,s2]=kstest(A,[A,p2],alpha)2 Q: a- G8 Z. ]* H: ]
    if H2==0
    ! D$ {  A) `4 ^( vdisp('该数据源服从γ分布。')9 \8 r# w7 Z  Y
    else& k3 C8 I7 j2 b. O) ~
    disp('该数据源不服从γ分布。'): z! Q' ?7 N& y: m- _! Z5 _  a; \6 ?2 P. [* T
    end
    2 v9 f, I6 p6 z9 plamda=poissfit(A,alpha);" H! D9 S- i* |) P1 H
    p3=poisscdf(A,lamda);+ \* y+ b) Y+ y: t
    [H3,s3]=kstest(A,[A,p3],alpha)
    6 i. |; j$ Q0 F" ~if H3==0  p8 ^4 P, Y# L6 U0 A- c
    disp('该数据源服从泊松分布。') 6 |) d) t* f- F9 _: c, Z
    else
    8 n+ }, X; p* V& P; Ndisp('该数据源不服从泊松分布。') 9 J1 r% p0 x  P* D) e4 F
    end
    6 ^4 F( q8 v3 M0 f  s! Qmu=expfit(A,alpha);: ?7 p6 p8 C5 V' l. U
    p4=expcdf(A,mu);
    ' n, T7 ?! n. R4 K/ e[H4,s4]=kstest(A,[A,p4],alpha)4 E: K; x  e* D/ m) F" R
    if H4==0# f' _6 F% T" _9 K
    disp('该数据源服从指数分布。') 1 G7 k- _8 v( q) \9 Z
    else
    2 x: ^( ?! Z- F$ [9 t7 H5 Edisp('该数据源不服从指数分布。')
      Y1 V( Q8 _3 N* Gend/ `$ K: k1 H1 t& j6 @
    [phat, pci] = raylfit(A, alpha)
    ' B, @- R1 a% E* lp5=raylcdf(A,phat);3 ~$ O$ }6 `4 U/ j
    [H5,s5]=kstest(A,[A,p5],alpha)
    : g! c* W0 g) R" W/ _! H/ Lif H5==0! n; @* T# a. [& Q2 g& ^
    disp('该数据源服从rayleigh分布。') ( U8 u2 u3 E: N( r
    else2 e  x1 m. _+ Y$ ?* w: {
    disp('该数据源不服从rayleigh分布。')7 M/ I' o9 A) f+ _
    end
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持1 反对反对0 微信微信

    5

    主题

    7

    听众

    145

    积分

    升级  22.5%

  • TA的每日心情
    郁闷
    2013-9-16 12:49
  • 签到天数: 28 天

    [LV.4]偶尔看看III

    自我介绍
    西政一学生

    群组2013认证赛B题讨论群组

    回复

    使用道具 举报

    1

    主题

    5

    听众

    26

    积分

    升级  22.11%

  • TA的每日心情
    开心
    2012-8-25 08:50
  • 签到天数: 5 天

    [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 0 d+ p) H6 f9 e* q# Q
    上面的是什么?

    . }7 F  z8 v0 I7 ~你是说整个么,是代码
    回复

    使用道具 举报

    0

    主题

    5

    听众

    57

    积分

    升级  54.74%

  • TA的每日心情
    难过
    2012-9-7 12:42
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    群组学术交流A

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-8-4 15:34 , Processed in 0.462729 second(s), 73 queries .

    回顶部