- 在线时间
- 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)
1 x7 ?5 @0 ~% o) \. oalpha=0.05;/ H2 \2 U8 V- e4 A, h7 R( b
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。5 U" p/ Q% q( G: ~
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';8 S/ o7 f i$ {/ q
[mu,sigma]=normfit(A);
7 m) Z& @6 y% ?$ ~6 y5 L2 Yp1=normcdf(A,mu,sigma);
: L) s! Y' s/ H! L[H1,s1]=kstest(A,[A,p1],alpha)6 {. @5 u$ S0 `
n=length(A);
( J( z0 `8 i. n7 s: l5 Iif H1==0) h3 N. w* r+ n
disp('该数据源服从正态分布。')3 W0 ^/ e" I0 W! ?! X+ \5 k
else2 Z! X4 d9 i% j6 _# N; I! i
disp('该数据源不服从正态分布。') , P7 {2 G! m- k( X& N
end
& L0 ~5 Z- H F" Kphat=gamfit(A,alpha);9 W; \* n$ L6 e) Y6 q
p2=gamcdf(A,phat(1),phat(2));
) w6 E. D$ z" t7 |1 w$ A; K4 g[H2,s2]=kstest(A,[A,p2],alpha)
* C9 K7 y3 G( p- S; |, J% ~" kif H2==0
% `( g% b9 l: a3 d, c5 T* ]5 w9 fdisp('该数据源服从γ分布。')( W3 C/ k$ A4 ^& C+ U; u2 Q
else
9 K. S- D" i5 u, g( ]disp('该数据源不服从γ分布。')8 s9 D: E7 y4 l" o
end
" W/ `5 t7 l% D' ~2 C- _9 plamda=poissfit(A,alpha);
" S. x( q+ d6 X( @( tp3=poisscdf(A,lamda);
7 g/ U8 \7 `1 A2 `- e# P' X2 C[H3,s3]=kstest(A,[A,p3],alpha)+ I: C1 N* Y% g* V$ d) e
if H3==0
9 a5 C R) Z/ t" ^% C" T; Adisp('该数据源服从泊松分布。') 3 j N9 e; @% ~6 h5 H
else, O* d0 ~" u) L* k) D
disp('该数据源不服从泊松分布。')
' U8 K- a3 @8 tend
' ^$ Z4 I" c) f k6 [mu=expfit(A,alpha);
W% ~2 S. W8 j4 @/ ]p4=expcdf(A,mu);
# }# V7 b4 L3 v1 @0 v6 H$ C& ?[H4,s4]=kstest(A,[A,p4],alpha)
) Z: I% f1 g% k2 L1 Q2 Hif H4==09 M- n, o! H. f$ l- Q0 N9 K
disp('该数据源服从指数分布。') & }( K4 d$ m; E' N0 q
else# X! K# T) |) n( Y) V
disp('该数据源不服从指数分布。') n! k& K$ q" |" u& t+ Y. K3 S; x
end
5 P2 } k$ S E* |3 z: T[phat, pci] = raylfit(A, alpha)
* z5 n0 X" K h/ Up5=raylcdf(A,phat);
% J2 P6 j$ q: r L7 V6 n7 ?[H5,s5]=kstest(A,[A,p5],alpha)
- B6 H+ b' V2 ^) i# W# m4 Aif H5==0/ g6 x: r" i; J r
disp('该数据源服从rayleigh分布。')
z0 l$ M5 s! M+ e3 g( }& s, qelse/ Y8 c, h9 a/ o. U' Q& x+ N6 C' t
disp('该数据源不服从rayleigh分布。') d& F% d _( x6 ]/ ^
end |
zan
|