QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 17708|回复: 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)$ k4 c9 \% Q0 v
    alpha=0.05;0 O0 [+ B. }. L2 K- ]1 v+ J$ M) ]
    % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。/ K" ]% E, J& _, K
    A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
    2 W$ o: l: ]" o; T' {5 }' ?( e; L[mu,sigma]=normfit(A);1 s; @: H; M! E2 a
    p1=normcdf(A,mu,sigma);( H$ A9 W. {4 s
    [H1,s1]=kstest(A,[A,p1],alpha)$ k: K/ j- b% }$ @% N; S; j
    n=length(A);
    4 J& e9 s/ _1 P  ]( t9 F% N, ?  Kif H1==0- ~( T, ?8 S. s
    disp('该数据源服从正态分布。')
    ; l/ i. n# V7 t' Y+ Zelse; b# \, O, z' O2 h; B
    disp('该数据源不服从正态分布。')
    6 Y( A% m& _+ T+ P! v( A8 Gend7 m% W2 Q2 T/ ?
    phat=gamfit(A,alpha);
    # n) P* D3 F0 H' _p2=gamcdf(A,phat(1),phat(2));0 y8 h$ [! \! ?. V  a, T* b# J& V
    [H2,s2]=kstest(A,[A,p2],alpha)8 e( l$ J, b( T
    if H2==0
    . ^) q3 A$ b1 _1 q8 ~. L" Bdisp('该数据源服从γ分布。')
    5 q- \4 u6 I  Selse5 X+ q  J; `' w
    disp('该数据源不服从γ分布。')+ c. Y3 n  E5 F) k
    end* Q7 y, ^* p+ }- T
    lamda=poissfit(A,alpha);
    * u4 Q; A2 D; j$ z) w; |) @p3=poisscdf(A,lamda);
    9 s8 v2 z, m' f( h9 o[H3,s3]=kstest(A,[A,p3],alpha)
    ( |/ F4 W8 x% J( p/ K7 }" xif H3==0
    3 t2 Z( B  d. Udisp('该数据源服从泊松分布。')
    ( o( D7 _9 A8 B3 P( @else: _" W# M: ?0 p0 q7 l
    disp('该数据源不服从泊松分布。')
    ) C4 R- U" l( D3 K( Tend
    # n: ~2 G+ G, e! x# m1 y1 n" I' wmu=expfit(A,alpha);, E  D2 g4 F. s
    p4=expcdf(A,mu);
    # u- D& r, `. o1 [$ o[H4,s4]=kstest(A,[A,p4],alpha), A5 e  V. T5 F- U2 q* x; s, L# h
    if H4==0% E8 {  J8 x+ l0 g8 \4 f
    disp('该数据源服从指数分布。')
    # M9 P* }  ]; y4 X7 Gelse
    7 e7 w$ X8 l- Odisp('该数据源不服从指数分布。') + m: u, `7 G  u* L
    end& L- l: o. k$ ?, l  I0 \
    [phat, pci] = raylfit(A, alpha)
    - `' q/ _. T& e/ Yp5=raylcdf(A,phat);
    1 f) l* N; h( v! M! t% A$ J) D[H5,s5]=kstest(A,[A,p5],alpha)
    * S1 g' g8 |7 s7 q; O( }8 vif H5==0
    5 j' t; K) t- V  q& A7 R. `" C' idisp('该数据源服从rayleigh分布。')
    ' S2 ~0 P7 c* Y% @: E# _! belse! t' e2 c( n  \) D. f4 S! ]
    disp('该数据源不服从rayleigh分布。')
    % j- ?& ~7 f' b% D6 l& W; d$ xend
    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
    6 w. p1 g& ?% z# ^* h! D9 G上面的是什么?

    5 H+ d# X5 ~3 S5 C* ]9 E1 m你是说整个么,是代码
    回复

    使用道具 举报

    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-7-28 09:25 , Processed in 0.338109 second(s), 72 queries .

    回顶部