QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17673|回复: 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)0 j) _( |4 Y8 D& ]  S- Z
    alpha=0.05;
    ' T' c' D/ n: F/ s  A1 Z% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
    9 ~* }2 ]/ _  H0 e! X# `A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    , o* G& [) Q" ]: j; Q[mu,sigma]=normfit(A);
    4 W/ [+ J: Q8 A% {% Mp1=normcdf(A,mu,sigma);
    ) }% [( U4 ^4 V  a[H1,s1]=kstest(A,[A,p1],alpha)
    % ^2 j: `( x1 Nn=length(A);
    ) j. ~) v" x: S# {' ^6 lif H1==09 k/ o- ]0 p! E9 Z4 ~
    disp('该数据源服从正态分布。')2 F5 |) q3 I, a
    else* ?4 l- Y9 \* X1 X
    disp('该数据源不服从正态分布。') % E& ?4 ^7 A7 c" o: T! D
    end5 _. n5 {! t! j2 U, g6 t
    phat=gamfit(A,alpha);/ }1 S. u- l3 t. }
    p2=gamcdf(A,phat(1),phat(2));$ V( ?% f& J  d, N" \
    [H2,s2]=kstest(A,[A,p2],alpha)7 {+ z8 R+ O; z. N9 n4 S$ e
    if H2==04 b& L* R+ E7 }( I* x5 K
    disp('该数据源服从γ分布。')
    . C; l6 P7 n. w- D) ]! Melse
    9 d; J+ ~6 e% Q6 F3 J8 j# v0 zdisp('该数据源不服从γ分布。')
    . k7 c( O# J" v% B4 ]end
    6 j7 S0 w& W- F) u6 Z0 X/ W- Clamda=poissfit(A,alpha);
    7 L3 B# Y+ \4 `p3=poisscdf(A,lamda);9 b8 k7 ?( @+ j
    [H3,s3]=kstest(A,[A,p3],alpha)1 f$ K* }* Z6 r( j8 W
    if H3==0: m2 \( O# C: ?+ K+ E
    disp('该数据源服从泊松分布。')
    8 C! r' B6 y8 K8 }* a1 kelse
    ) u8 t- d2 K0 u! u* s6 H/ gdisp('该数据源不服从泊松分布。')
    1 _3 E' N: ~6 j/ D9 l6 jend
    ; ~6 Q' A4 J" q7 A' h, Omu=expfit(A,alpha);, i8 U/ `& L% l- k
    p4=expcdf(A,mu);* h9 G$ f6 F0 \. Q4 U$ b0 J
    [H4,s4]=kstest(A,[A,p4],alpha)
    9 V, Y2 u# m' |; I; c' ?2 f" \% bif H4==0' U( R/ k# A  D. S* K
    disp('该数据源服从指数分布。') . _/ J' s3 }2 ~' \; s$ u' s
    else& v& b  u6 c7 ?; R. p3 U) C: D9 v
    disp('该数据源不服从指数分布。') * @; P! e1 }; p# {: L' @
    end( `3 T# F. V1 U& r5 A' R  X& l
    [phat, pci] = raylfit(A, alpha)
    2 O3 |: _0 F: ep5=raylcdf(A,phat);
    , C9 |" {: W+ h$ N[H5,s5]=kstest(A,[A,p5],alpha)4 S' ]9 s5 P! j6 N) H
    if H5==07 {2 W- n1 F/ y; C( \! e; u
    disp('该数据源服从rayleigh分布。') 8 r3 Q( J- d" w2 e9 ^) X  G
    else, \7 {; t/ g1 H: s1 _0 T4 s: [
    disp('该数据源不服从rayleigh分布。')
    4 t5 n/ x/ N: g& _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 % z" R4 O* t2 Y7 y: E
    上面的是什么?

    ' ~. Z% @0 d2 l& \2 j  [你是说整个么,是代码
    回复

    使用道具 举报

    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-6-14 06:25 , Processed in 0.457444 second(s), 72 queries .

    回顶部