QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17769|回复: 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)
    8 n  z- p: I% M! ~2 ^3 zalpha=0.05;/ W. O, m' h2 K9 n& ?( f
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。: `0 X1 ^8 i; V5 u) n" a0 G
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    $ V7 r7 A3 _# Z  A, ^6 f* \[mu,sigma]=normfit(A);3 u6 D% s$ ?! u1 b6 B4 N) w" p: P7 ?* Z
    p1=normcdf(A,mu,sigma);$ D1 C. M, ^" h6 X6 n. ~8 o
    [H1,s1]=kstest(A,[A,p1],alpha)' l6 \9 V3 x$ e4 i6 \3 y3 ^  ?+ B
    n=length(A);* [* c" e* E+ w6 q+ a  n3 r5 J
    if H1==0% X3 u& m! D. X. A6 M4 C
    disp('该数据源服从正态分布。')
    $ n1 V; q- ?+ s- I  [. [; Qelse
    1 I2 s) o0 \5 q$ m- I+ F& y3 L7 f' hdisp('该数据源不服从正态分布。') 6 |7 e3 T: [, l- W3 N
    end
    ( ^+ X* n- q4 s+ M" hphat=gamfit(A,alpha);
    ' t4 ~6 M% L" ^5 ?4 rp2=gamcdf(A,phat(1),phat(2));% \- A5 q7 W5 P8 I/ q, R4 }
    [H2,s2]=kstest(A,[A,p2],alpha)
    5 _) p% Z! M% n  n& fif H2==0
      m1 C! v$ [9 N' K+ fdisp('该数据源服从γ分布。')* c; y- v+ g; Q. q6 H
    else
    $ @' J3 ^% c5 D. g* s; c, xdisp('该数据源不服从γ分布。')
    * p& }8 a  a7 G* _  [end
    ! X4 @( a9 F8 B$ A& W- D& a; \lamda=poissfit(A,alpha);
    : ]" s$ @! A1 o% y5 R8 x, sp3=poisscdf(A,lamda);
    ; u4 T, o: i- h8 X* B8 Z" e[H3,s3]=kstest(A,[A,p3],alpha)4 U- J% y1 N" |
    if H3==0
    ( L9 c) ]1 _9 hdisp('该数据源服从泊松分布。') # {7 d7 m6 ^1 ]
    else
    8 k% V0 s, Q3 T7 \! ~9 K4 Idisp('该数据源不服从泊松分布。')
    + S# ]. e3 X0 iend
    7 D" J  q4 v( F2 r1 W* ~' U) Rmu=expfit(A,alpha);
    ' O( C$ w' D; Z) ]3 vp4=expcdf(A,mu);$ g8 U: O; k4 B) U9 |' _/ ^
    [H4,s4]=kstest(A,[A,p4],alpha)
    % R5 E, I0 n: \6 Eif H4==08 Y! T/ L9 S& Z- O
    disp('该数据源服从指数分布。')
    - x4 n' L4 y. W) N, Zelse4 Z2 ^+ e0 s+ I$ y
    disp('该数据源不服从指数分布。') % r3 S. a$ P1 a1 o
    end! i# V, S# x1 [# `
    [phat, pci] = raylfit(A, alpha)
    $ q* o# o+ p  Fp5=raylcdf(A,phat);
    % K) A" J: f& n% b[H5,s5]=kstest(A,[A,p5],alpha)5 ~; h" c, f8 k/ n: q; A) S+ a1 t- e' m
    if H5==0
    2 N* Z$ \" D7 }4 A8 `: Z* t; ]disp('该数据源服从rayleigh分布。')
    5 \2 Z" l# N# d& ]else' I/ w/ s+ m- p
    disp('该数据源不服从rayleigh分布。')# h' s  u, G" x# u
    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
    7 l+ f4 j4 F, z- ~3 J" @* r上面的是什么?
    5 z& Q' r( r. a/ 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-25 01:07 , Processed in 2.700361 second(s), 71 queries .

    回顶部