- 在线时间
- 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)$ g4 j" ^2 t! C% O
alpha=0.05;
: K3 {0 ?5 F* B( g% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
' p8 ^. ]( Z; Q& g' \% L7 xA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
- _3 D. @; A+ U7 A. d3 w5 x) v[mu,sigma]=normfit(A);
3 R3 I- p, V1 f9 sp1=normcdf(A,mu,sigma);( @: ^% D: W' N. e
[H1,s1]=kstest(A,[A,p1],alpha)# ^6 k7 g G8 V3 b O4 G
n=length(A);# D& W5 b0 F8 ?# D1 z
if H1==0: q/ H6 W4 I3 V" _
disp('该数据源服从正态分布。')
* d8 F* \9 C2 s9 [9 f( lelse/ B; L8 J9 g+ V6 O
disp('该数据源不服从正态分布。')
/ ^8 ? K$ T1 z$ rend' r+ ?4 P, {- K2 q0 F
phat=gamfit(A,alpha);
0 M$ m o' L1 x9 Hp2=gamcdf(A,phat(1),phat(2));1 C! ~) b9 l9 C' H8 ?8 H, @
[H2,s2]=kstest(A,[A,p2],alpha)
8 F* w* y$ n; M0 u( k; uif H2==07 J& D' u; u2 Q6 }7 ]
disp('该数据源服从γ分布。')* d+ w6 D' t( N# @! e) ~
else; _3 H. U8 i1 q, ~2 y: D
disp('该数据源不服从γ分布。')0 B1 x( U W+ c" B& [8 {) [, [/ J
end
6 v2 R! i+ b+ @, b- Nlamda=poissfit(A,alpha);- R% y) R1 c% ~1 m1 b
p3=poisscdf(A,lamda);
6 T3 O2 t) D+ V" N% f" D& E[H3,s3]=kstest(A,[A,p3],alpha)# n {8 L R) k
if H3==03 N0 j7 j) W L8 t T8 Z& t
disp('该数据源服从泊松分布。') 3 Y$ p: m$ i. R
else
/ }' Q+ \/ |9 N: Jdisp('该数据源不服从泊松分布。') , i+ n6 ]/ q1 L/ x) i9 s( ~% a+ N) q& p
end
8 j# b. E& k+ imu=expfit(A,alpha);
# k: z0 d& ^7 G" g5 h+ up4=expcdf(A,mu);
8 f z% K5 P* k2 Z[H4,s4]=kstest(A,[A,p4],alpha)
5 Q2 J+ K6 U( Y# ` Uif H4==0
- H/ G: F; P$ H8 o t1 pdisp('该数据源服从指数分布。') / [0 P! S) u" m2 _% C4 C
else
4 _% E! ?! o) T; z2 n8 Q/ mdisp('该数据源不服从指数分布。') 0 i2 x1 o' M- {4 l: s
end) o k5 p1 y6 F
[phat, pci] = raylfit(A, alpha): F: C& K+ Z! M+ d, l
p5=raylcdf(A,phat);
- P6 b3 W" \" `; U+ y+ C& }+ d[H5,s5]=kstest(A,[A,p5],alpha)
1 K/ N+ \6 t5 L0 z8 }' Pif H5==0% y) d: s7 k, K( O9 _0 Q% c1 b, D
disp('该数据源服从rayleigh分布。')
# b$ E9 ]- P* n6 G2 P9 ^2 R- T9 Selse
. N* N( G" z/ q) m; k4 Udisp('该数据源不服从rayleigh分布。')1 X" s( a( O3 Z" |. i ]6 M$ g2 a! S* [
end |
zan
|