- 在线时间
- 34 小时
- 最后登录
- 2014-5-13
- 注册时间
- 2012-2-7
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 370 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 163
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 96
- 主题
- 5
- 精华
- 0
- 分享
- 0
- 好友
- 37
升级   31.5% TA的每日心情 | 无聊 2014-3-17 18:33 |
|---|
签到天数: 54 天 [LV.5]常住居民I
 群组: 学术交流B |
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
|