- 在线时间
- 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)" v( [% G; ]$ s7 h, V
alpha=0.05;# n/ V; Q$ J8 n6 O2 ]
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
& W$ X4 _/ Y% M4 B4 w9 `4 kA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';/ \) v8 P d/ y0 D( K# r1 o
[mu,sigma]=normfit(A);
$ T2 {6 N' _ Qp1=normcdf(A,mu,sigma);6 j4 b0 g& ~0 y) z
[H1,s1]=kstest(A,[A,p1],alpha)9 p. Q" b- |$ ?) V( ?
n=length(A);
6 q9 B0 x" E$ [. |$ U9 x9 b' tif H1==0
# s( o. j; C g) j1 x {disp('该数据源服从正态分布。')
) `# Z9 { t/ j# ]else
) |; K4 r/ c2 S2 V1 xdisp('该数据源不服从正态分布。') ; g& t. y0 f* d
end4 R* B) `2 z) ?$ d [' ~( r
phat=gamfit(A,alpha);
" c5 G7 ^0 H1 n% e# D! J) g4 Rp2=gamcdf(A,phat(1),phat(2));
2 o, y4 t! G* G4 V[H2,s2]=kstest(A,[A,p2],alpha)
( W" C4 E( B4 v; J; d( ^if H2==0
+ z+ M2 q, w- v: M; ]. t( Idisp('该数据源服从γ分布。')" W( l7 B1 O2 x! m5 O1 J1 @# F
else! S# {8 X- L3 G! v% h/ q3 b
disp('该数据源不服从γ分布。')6 ]6 F! H# r7 U
end
" J# g) P1 |3 Flamda=poissfit(A,alpha);2 l, J% h& h0 r% d2 D
p3=poisscdf(A,lamda);
. z( c1 ]% k1 G( \$ r/ S6 \[H3,s3]=kstest(A,[A,p3],alpha)
1 x4 g2 D" @0 e0 V0 a! `if H3==0
p( |. h6 j* @ Pdisp('该数据源服从泊松分布。') " x T3 J1 M* A& K% p" L3 A$ k) `4 `
else) _* ?$ I2 s7 ?. }. o3 R2 ?
disp('该数据源不服从泊松分布。')
9 T- u6 t1 w" {) G* t- wend
5 V7 i6 b0 g) Z, b% W8 E" Smu=expfit(A,alpha);7 Z% F! W0 Y6 `7 w6 K
p4=expcdf(A,mu);
% W# p k" }. ]% I3 U' }[H4,s4]=kstest(A,[A,p4],alpha)
) t8 n S* r: B: o; V4 O- d L1 m7 x% vif H4==0
$ \8 S6 p% z9 I/ q! Xdisp('该数据源服从指数分布。') 8 e2 M4 _& e' |2 b
else
* w9 z% q5 _ m) n+ Y/ o: Q/ o( K/ Jdisp('该数据源不服从指数分布。') / U/ s+ S2 G8 A7 \3 [. X
end
0 a( L7 m& ]' t9 _# n5 `[phat, pci] = raylfit(A, alpha)0 J5 v+ m g# J
p5=raylcdf(A,phat);- W! M4 W& v+ U/ L' x8 w
[H5,s5]=kstest(A,[A,p5],alpha)
: e l% B' n& Z" p% U5 ?) }if H5==0
_5 f5 H5 K; g; N, l5 F4 Udisp('该数据源服从rayleigh分布。') 6 I/ g; R; s% V7 }& a4 s
else0 R: r0 `: G' A9 p; R u
disp('该数据源不服从rayleigh分布。')0 a0 Z+ z4 z7 j( c
end |
zan
|