QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17768|回复: 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 a0 @# D" L+ ]9 W* O
    alpha=0.05;6 ]9 C( A4 E6 h" `% }+ q2 Q* l
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    5 H4 u& @0 a4 U! U! I4 Q2 IA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    8 @7 o) k$ z0 b4 B[mu,sigma]=normfit(A);. A5 c$ x: ^# k7 w+ c+ G7 m- [
    p1=normcdf(A,mu,sigma);0 ^  W( O! j' l/ A) @+ `+ w
    [H1,s1]=kstest(A,[A,p1],alpha)
    ; w1 s. ~( f6 A, Mn=length(A);
    2 I. ?( |- n' \# U. pif H1==0
    2 I5 i; S% V& D" A* [0 b3 bdisp('该数据源服从正态分布。')8 S0 Q7 i) s9 F
    else0 {0 \. N; \$ ^3 P
    disp('该数据源不服从正态分布。') . B6 K. `: S( u% ]" ]
    end8 T" c( s5 F, a3 E% Z& j7 u
    phat=gamfit(A,alpha);
    5 b3 T/ a% s% S7 O# e2 Vp2=gamcdf(A,phat(1),phat(2));- v- b, ?. `9 v- K- D
    [H2,s2]=kstest(A,[A,p2],alpha)
    8 j2 ^; H- m3 t+ H8 wif H2==0$ P% ?7 |7 A* P. A2 I6 v* |/ l
    disp('该数据源服从γ分布。')! A8 J& T. y0 d
    else5 ~/ O% o+ @5 H
    disp('该数据源不服从γ分布。')) V( Q  G  i6 i* V
    end- D  O' o- X' H9 Y  D' E% u5 q
    lamda=poissfit(A,alpha);$ x7 y5 v: o2 ]) U
    p3=poisscdf(A,lamda);
    5 h) w0 ~8 G4 |: p5 [2 m% B1 }[H3,s3]=kstest(A,[A,p3],alpha)  B1 A- r- z& I9 |; k$ H
    if H3==0& a, N! X6 k2 ~: ?+ Z# \
    disp('该数据源服从泊松分布。') " n, g1 V$ ]/ v/ Z! ^
    else5 U* ?5 X8 a) I# x6 U* R+ F
    disp('该数据源不服从泊松分布。')
    4 X3 T8 t5 U  b( q2 }7 ]end
    # j3 u& p) ]: r: Cmu=expfit(A,alpha);/ W" J* g- a# C
    p4=expcdf(A,mu);
    / a; z6 d6 p9 b[H4,s4]=kstest(A,[A,p4],alpha)+ K0 y. r3 g- x1 N" q5 ]( G
    if H4==0, v& V1 S+ \! A. Z
    disp('该数据源服从指数分布。')
    9 I4 U: N, J, Z- H% l, e) {else
    0 A% {  x% r8 O1 z0 kdisp('该数据源不服从指数分布。')
      n6 `+ P/ J" n* D( y  O& Xend
    # j  K, A8 H- \. H$ b2 c! y; L5 V[phat, pci] = raylfit(A, alpha)- l+ y4 X4 I1 B: o
    p5=raylcdf(A,phat);7 n" y; b0 d" [4 T8 p& ]
    [H5,s5]=kstest(A,[A,p5],alpha)
    2 Q8 L; P4 n6 i( H' e- xif H5==0& i* Y2 V% t1 ^$ i7 |' o6 c. T
    disp('该数据源服从rayleigh分布。') " p) }7 D* Y+ [4 o
    else
    & P$ A: o* m; s& m# m; P% Vdisp('该数据源不服从rayleigh分布。')
    1 p) |; V" x  bend
    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 + k- i) f# E- [  [
    上面的是什么?
    + i. y/ @, S2 o4 p! F6 ]
    你是说整个么,是代码
    回复

    使用道具 举报

    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-9-24 22:23 , Processed in 0.378602 second(s), 72 queries .

    回顶部