- 在线时间
- 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)+ H$ Y5 B5 Q5 ^
alpha=0.05;* Y! d: g$ A% x" }8 \0 `8 E
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。- d% e3 F- v( v3 |
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
9 I& K/ h5 D2 C[mu,sigma]=normfit(A);( o6 C+ q9 T$ b0 _1 P
p1=normcdf(A,mu,sigma);
6 ^7 J9 V$ Q; x; ]" ?. n[H1,s1]=kstest(A,[A,p1],alpha)" Z7 t5 f2 c5 {) x
n=length(A);
1 n8 a9 U" G. D$ \$ pif H1==0" B' x6 s0 n" D, C* V: |3 W
disp('该数据源服从正态分布。')
Y w6 o; w: y6 j" X; Uelse
. R! g: G3 L/ V: L8 P: y; _/ ^1 }disp('该数据源不服从正态分布。')
+ l, x/ H$ ~. send
( F0 z8 V. @( t Z$ i/ J/ i9 bphat=gamfit(A,alpha);
6 M% w+ |& l2 lp2=gamcdf(A,phat(1),phat(2));; |' j4 J3 H) f
[H2,s2]=kstest(A,[A,p2],alpha)
2 e5 {" I! Z$ f1 W: qif H2==0& ]/ {' t. w9 x! A2 t$ c/ U I1 p
disp('该数据源服从γ分布。'), R$ C6 i, p2 b, ]
else
& w7 n2 u) S ~+ {disp('该数据源不服从γ分布。')& ?9 z+ s" Y$ F8 R4 b+ u! K+ N6 G
end2 h4 d. F, V! i7 n7 `/ m
lamda=poissfit(A,alpha);
2 s: `; l. o& i! K1 q/ Tp3=poisscdf(A,lamda);
& B. W, D6 X2 l0 p1 Y[H3,s3]=kstest(A,[A,p3],alpha)) @" y7 R9 `- | @
if H3==0
9 Q0 u! z5 W- [1 f9 Pdisp('该数据源服从泊松分布。')
9 x: `3 ? ~7 U2 m& Felse
5 H7 c5 w; X1 fdisp('该数据源不服从泊松分布。') ) B) `7 k* c% ~/ j
end4 ?+ X, h2 D9 [+ x" C4 E g8 q
mu=expfit(A,alpha);
+ k* R& ~5 H `) j6 m6 D) Q1 Gp4=expcdf(A,mu);$ V* c4 W/ q# z' w5 L$ x$ c
[H4,s4]=kstest(A,[A,p4],alpha)0 J ]/ ` w- K% J+ m7 q
if H4==0: w2 b! G2 O5 p% u# H7 m
disp('该数据源服从指数分布。')
! L) }8 t5 S5 F: t, l Y' Z" G" ^ |6 Nelse4 [: u$ M( [+ T8 P" l
disp('该数据源不服从指数分布。') 6 r2 Z0 \" Y; e7 r" M
end7 F) K- d8 h3 h7 Y0 I
[phat, pci] = raylfit(A, alpha)! N) y1 t, K8 Z+ t+ c. e, f8 B! `
p5=raylcdf(A,phat);
" f& ]: Z+ [. E* `# e" ~- [4 ?' n[H5,s5]=kstest(A,[A,p5],alpha)
: s3 S, J) V9 b; R* b) L! ]if H5==09 A% u V; J8 X" j( j* n
disp('该数据源服从rayleigh分布。') $ Q) v! K; A7 O( d' t
else) s& H8 m8 X7 N- n1 L# L
disp('该数据源不服从rayleigh分布。')1 N. F1 F9 }7 d" l6 r$ P1 o
end |
zan
|