- 在线时间
- 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* w; W6 _& n0 {" f+ _" kalpha=0.05;. P6 b- Q; W. M W. m
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
1 C2 L T0 D4 `+ [7 ]) jA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';1 y; H3 I" ]. E7 z$ u( O
[mu,sigma]=normfit(A); c( H( }! P* _2 O* D9 p8 Z
p1=normcdf(A,mu,sigma);
/ v# w9 v9 f# h8 J: H% n' C8 U[H1,s1]=kstest(A,[A,p1],alpha)
1 Y6 j3 e$ N& f- o- Jn=length(A);, x8 U; _" d0 |2 ]. x' _
if H1==0* Y; M! i5 D T& ^& h3 [
disp('该数据源服从正态分布。')
{ w" s2 L+ yelse7 A/ q r9 G) y
disp('该数据源不服从正态分布。') 9 r: l/ U% W( M% h* K9 f
end4 y1 B: a% L$ I |
phat=gamfit(A,alpha);
$ ~( n% j- a: H% A# M) \p2=gamcdf(A,phat(1),phat(2));" h0 p- y$ a5 [! R2 f
[H2,s2]=kstest(A,[A,p2],alpha)
; X2 Q* r2 J. c9 |8 V, ^8 B9 Nif H2==07 }& h6 f# A: K/ N9 h
disp('该数据源服从γ分布。')
1 c6 G4 i6 \- E( Aelse
' k* T; k C4 r: S& T0 E) O1 Pdisp('该数据源不服从γ分布。')% ~% @) ]1 i! P D5 F
end( [9 n; _1 P: h! {" T& P
lamda=poissfit(A,alpha);- Q! Y9 R( \" F a$ I% Y
p3=poisscdf(A,lamda);
1 S& ^) N' ]$ H0 ^+ l[H3,s3]=kstest(A,[A,p3],alpha)8 I' Q5 w5 F$ \- z$ z" V
if H3==0
0 o0 u% q( |; S. fdisp('该数据源服从泊松分布。')
) \; ?' h. s% @3 Nelse0 z' X+ j: O4 o% C# x) B% [
disp('该数据源不服从泊松分布。')
2 q+ N `+ t) {6 f. J$ ^1 q- V& _end
( i! y( ^9 o6 E7 Smu=expfit(A,alpha);
5 g! M9 }8 _4 e1 W* \" ?& Hp4=expcdf(A,mu);
6 {2 x) p& C. Z/ e[H4,s4]=kstest(A,[A,p4],alpha). Z) W- L( \' P) _
if H4==0
, ]3 l/ | p" ~: Bdisp('该数据源服从指数分布。') 6 l& O" o3 X$ J0 r. c
else; `& @' \" g+ G) {' z6 i" D; h
disp('该数据源不服从指数分布。') 4 a1 b$ e6 ~) `) W( F/ F
end+ t y# I4 p) l) x7 b( o! D
[phat, pci] = raylfit(A, alpha) N9 L$ O* l9 c6 |" v
p5=raylcdf(A,phat);; F" }2 l4 |1 c# K3 e5 m0 R( I
[H5,s5]=kstest(A,[A,p5],alpha)
7 y1 F. s& G/ _+ Iif H5==03 h0 K/ U2 l1 x$ p' K6 g9 O
disp('该数据源服从rayleigh分布。')
. c3 ^0 w2 V* g& p1 E' |else; x8 {4 I& d. }/ Q. V
disp('该数据源不服从rayleigh分布。')# q: Z8 }* v$ ?8 c0 x
end |
zan
|