- 在线时间
- 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)
4 C) O, B7 I7 _4 C5 p0 _alpha=0.05;- `7 @* k+ [! Y% z
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。. y$ m* L# w6 `, f
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
+ m! B+ L. [- J" ?[mu,sigma]=normfit(A);
( b- ]' X! B2 F2 xp1=normcdf(A,mu,sigma);
5 G8 ~1 C% G! ?1 A, T9 p# u9 }[H1,s1]=kstest(A,[A,p1],alpha)2 b' v9 }/ x& l2 u
n=length(A);% {8 I% P0 N+ Z
if H1==0
4 ~$ D' y! g: C- s9 o' e G3 wdisp('该数据源服从正态分布。')
7 V7 P( m% |4 \9 U0 belse: Q3 S* T l3 M
disp('该数据源不服从正态分布。')
; I$ W* u) h- U, lend/ q# g- u1 m: g, w+ s' O
phat=gamfit(A,alpha);. k9 f# o4 p. r* j% u6 M
p2=gamcdf(A,phat(1),phat(2));
# g/ ^( s3 y: S[H2,s2]=kstest(A,[A,p2],alpha)5 u% V0 t0 H+ W+ n6 p
if H2==0+ S! X0 }% R1 z7 E8 c
disp('该数据源服从γ分布。')
# b# g& O. S$ m6 a& z+ x Nelse% }( y3 W' v4 ]% M8 L: I, q
disp('该数据源不服从γ分布。')
1 t# D1 Q0 _3 s3 \end0 M: ?# |6 m# ?5 P2 r
lamda=poissfit(A,alpha);& x5 K A, v$ J
p3=poisscdf(A,lamda);# V9 Q# B8 d3 o) n+ f2 X8 ]
[H3,s3]=kstest(A,[A,p3],alpha), ~% q. b' d/ D' @. E4 z
if H3==0* T Q2 T$ p8 D% D6 n
disp('该数据源服从泊松分布。')
6 n$ E) q' C$ B1 T3 o Lelse
# ^/ X3 W9 ?8 T7 Xdisp('该数据源不服从泊松分布。') ) N# R' x. M2 F/ H; m) N
end/ `9 l, A6 q9 ^& t1 t1 I5 H
mu=expfit(A,alpha);' R N" w9 B$ |8 A% d5 V+ K( H7 ^
p4=expcdf(A,mu);
0 X2 e3 B5 t4 t: `/ m. R- |. K[H4,s4]=kstest(A,[A,p4],alpha)" @5 Y. o' k$ K1 q
if H4==0; }( x# c; [6 F9 N+ }
disp('该数据源服从指数分布。') 0 i$ U5 _% y6 j8 m: p! o
else% x" g. j) p3 t
disp('该数据源不服从指数分布。')
8 _4 T" Z; Q; U/ Send6 ~' j* R) r4 ?5 V4 b
[phat, pci] = raylfit(A, alpha)
. X2 M z7 J+ w! ?( I' d z Cp5=raylcdf(A,phat);* ]* W& ~6 p. z: R
[H5,s5]=kstest(A,[A,p5],alpha)' E/ \( Y3 _) E2 S6 v
if H5==0
1 J) C) h' W1 \( g& [disp('该数据源服从rayleigh分布。') ' i; f" C7 S1 o- d/ Y9 M" q0 @
else
1 y$ ~" F: K* Y7 I$ t9 Ydisp('该数据源不服从rayleigh分布。')8 Z+ Z* H0 a4 |3 ~
end |
zan
|