QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17767|回复: 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 l$ E* A  S" J& M% {1 E, K! u3 [
    alpha=0.05;1 S7 L% e$ Q& z6 c8 S. f
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。" X) @5 u! ^0 ?3 m- ?
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';* P1 N  D& {3 G; k, I
    [mu,sigma]=normfit(A);% I9 M" ?$ H+ X1 h; ?
    p1=normcdf(A,mu,sigma);; N/ v/ O3 u% L8 h' b/ i; K/ x
    [H1,s1]=kstest(A,[A,p1],alpha)
    " k$ n( e9 `0 i* h5 i) Pn=length(A);
    + x0 }" v4 f3 Aif H1==0( V; b  J6 `8 l, _: u
    disp('该数据源服从正态分布。')6 ]4 N7 G; M" ~5 w  w1 m
    else. h' O) `7 C/ l  B6 t7 m
    disp('该数据源不服从正态分布。')
    : T3 {" `+ R4 n4 d: n5 {3 D$ d# Yend' p& u& Y& x( c4 s5 g' q% r
    phat=gamfit(A,alpha);
    1 m0 |2 W! z& N8 C0 ^; Op2=gamcdf(A,phat(1),phat(2));" V0 B6 T# ^7 ~& E$ Q# f8 P
    [H2,s2]=kstest(A,[A,p2],alpha)
    ! i, o9 G8 j5 N& Vif H2==0
    ) u. a( X9 n9 r" gdisp('该数据源服从γ分布。')
    3 V9 G' X9 ~9 A* S- j# p6 Qelse; i; ~5 p  V' X, V3 f  ~5 q
    disp('该数据源不服从γ分布。')
    , u7 z2 `' T6 C4 pend
    / b4 l+ m9 `6 h8 d0 Xlamda=poissfit(A,alpha);2 f5 K6 F  d8 M4 [; o# ~
    p3=poisscdf(A,lamda);
    6 v) g& Q/ z1 W0 F' }[H3,s3]=kstest(A,[A,p3],alpha): w6 T0 t3 M7 w; l* O6 ?
    if H3==04 U% D; g, v- i  V2 k" o. y
    disp('该数据源服从泊松分布。')
    1 K" |" S" U& d& |" z0 aelse; o5 [/ d. W4 c1 Q4 x
    disp('该数据源不服从泊松分布。') 4 B5 O) t  T( ^- G" B
    end' z$ O0 u( i5 f" a. |% G  O
    mu=expfit(A,alpha);1 {9 N! L; }( C0 N" H
    p4=expcdf(A,mu);
    % R. i; _+ t* K( I. }# L: A& p0 d[H4,s4]=kstest(A,[A,p4],alpha)
    2 {9 i* P: I7 S3 Xif H4==0- q( E' E! h- A9 f: k2 ^' W. c* e
    disp('该数据源服从指数分布。') , o; Q* U/ Y9 X1 L
    else, X* P% G6 f8 W! b& s8 _
    disp('该数据源不服从指数分布。') " C" p" r, p! {) L4 X+ s
    end1 C! P; w- l$ v
    [phat, pci] = raylfit(A, alpha)
    $ P; g9 H% l% \9 [* O7 V% _p5=raylcdf(A,phat);
    1 e) _5 Z1 P2 I# F[H5,s5]=kstest(A,[A,p5],alpha)( u  F3 D. E) t3 T# g' }
    if H5==0
    1 M/ @& d+ L8 p' r6 g, ^" Xdisp('该数据源服从rayleigh分布。')
    3 ]& h! C4 W. \else7 a  O+ q$ |0 E
    disp('该数据源不服从rayleigh分布。'). t6 J# i: [, T$ {' 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 ; B1 O& \8 Q/ p. P
    上面的是什么?

    * U' |; g+ T$ p' O) O$ @你是说整个么,是代码
    回复

    使用道具 举报

    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 21:10 , Processed in 0.411575 second(s), 72 queries .

    回顶部