- 在线时间
- 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 a0 @# D" L+ ]9 W* O
alpha=0.05;6 ]9 C( A4 E6 h" `% }+ q2 Q* l
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
5 H4 u& @0 a4 U! U! I4 Q2 IA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
8 @7 o) k$ z0 b4 B[mu,sigma]=normfit(A);. A5 c$ x: ^# k7 w+ c+ G7 m- [
p1=normcdf(A,mu,sigma);0 ^ W( O! j' l/ A) @+ `+ w
[H1,s1]=kstest(A,[A,p1],alpha)
; w1 s. ~( f6 A, Mn=length(A);
2 I. ?( |- n' \# U. pif H1==0
2 I5 i; S% V& D" A* [0 b3 bdisp('该数据源服从正态分布。')8 S0 Q7 i) s9 F
else0 {0 \. N; \$ ^3 P
disp('该数据源不服从正态分布。') . B6 K. `: S( u% ]" ]
end8 T" c( s5 F, a3 E% Z& j7 u
phat=gamfit(A,alpha);
5 b3 T/ a% s% S7 O# e2 Vp2=gamcdf(A,phat(1),phat(2));- v- b, ?. `9 v- K- D
[H2,s2]=kstest(A,[A,p2],alpha)
8 j2 ^; H- m3 t+ H8 wif H2==0$ P% ?7 |7 A* P. A2 I6 v* |/ l
disp('该数据源服从γ分布。')! A8 J& T. y0 d
else5 ~/ O% o+ @5 H
disp('该数据源不服从γ分布。')) V( Q G i6 i* V
end- D O' o- X' H9 Y D' E% u5 q
lamda=poissfit(A,alpha);$ x7 y5 v: o2 ]) U
p3=poisscdf(A,lamda);
5 h) w0 ~8 G4 |: p5 [2 m% B1 }[H3,s3]=kstest(A,[A,p3],alpha) B1 A- r- z& I9 |; k$ H
if H3==0& a, N! X6 k2 ~: ?+ Z# \
disp('该数据源服从泊松分布。') " n, g1 V$ ]/ v/ Z! ^
else5 U* ?5 X8 a) I# x6 U* R+ F
disp('该数据源不服从泊松分布。')
4 X3 T8 t5 U b( q2 }7 ]end
# j3 u& p) ]: r: Cmu=expfit(A,alpha);/ W" J* g- a# C
p4=expcdf(A,mu);
/ a; z6 d6 p9 b[H4,s4]=kstest(A,[A,p4],alpha)+ K0 y. r3 g- x1 N" q5 ]( G
if H4==0, v& V1 S+ \! A. Z
disp('该数据源服从指数分布。')
9 I4 U: N, J, Z- H% l, e) {else
0 A% { x% r8 O1 z0 kdisp('该数据源不服从指数分布。')
n6 `+ P/ J" n* D( y O& Xend
# j K, A8 H- \. H$ b2 c! y; L5 V[phat, pci] = raylfit(A, alpha)- l+ y4 X4 I1 B: o
p5=raylcdf(A,phat);7 n" y; b0 d" [4 T8 p& ]
[H5,s5]=kstest(A,[A,p5],alpha)
2 Q8 L; P4 n6 i( H' e- xif H5==0& i* Y2 V% t1 ^$ i7 |' o6 c. T
disp('该数据源服从rayleigh分布。') " p) }7 D* Y+ [4 o
else
& P$ A: o* m; s& m# m; P% Vdisp('该数据源不服从rayleigh分布。')
1 p) |; V" x bend |
zan
|