- 在线时间
- 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)0 j) _( |4 Y8 D& ] S- Z
alpha=0.05;
' T' c' D/ n: F/ s A1 Z% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
9 ~* }2 ]/ _ H0 e! X# `A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
, o* G& [) Q" ]: j; Q[mu,sigma]=normfit(A);
4 W/ [+ J: Q8 A% {% Mp1=normcdf(A,mu,sigma);
) }% [( U4 ^4 V a[H1,s1]=kstest(A,[A,p1],alpha)
% ^2 j: `( x1 Nn=length(A);
) j. ~) v" x: S# {' ^6 lif H1==09 k/ o- ]0 p! E9 Z4 ~
disp('该数据源服从正态分布。')2 F5 |) q3 I, a
else* ?4 l- Y9 \* X1 X
disp('该数据源不服从正态分布。') % E& ?4 ^7 A7 c" o: T! D
end5 _. n5 {! t! j2 U, g6 t
phat=gamfit(A,alpha);/ }1 S. u- l3 t. }
p2=gamcdf(A,phat(1),phat(2));$ V( ?% f& J d, N" \
[H2,s2]=kstest(A,[A,p2],alpha)7 {+ z8 R+ O; z. N9 n4 S$ e
if H2==04 b& L* R+ E7 }( I* x5 K
disp('该数据源服从γ分布。')
. C; l6 P7 n. w- D) ]! Melse
9 d; J+ ~6 e% Q6 F3 J8 j# v0 zdisp('该数据源不服从γ分布。')
. k7 c( O# J" v% B4 ]end
6 j7 S0 w& W- F) u6 Z0 X/ W- Clamda=poissfit(A,alpha);
7 L3 B# Y+ \4 `p3=poisscdf(A,lamda);9 b8 k7 ?( @+ j
[H3,s3]=kstest(A,[A,p3],alpha)1 f$ K* }* Z6 r( j8 W
if H3==0: m2 \( O# C: ?+ K+ E
disp('该数据源服从泊松分布。')
8 C! r' B6 y8 K8 }* a1 kelse
) u8 t- d2 K0 u! u* s6 H/ gdisp('该数据源不服从泊松分布。')
1 _3 E' N: ~6 j/ D9 l6 jend
; ~6 Q' A4 J" q7 A' h, Omu=expfit(A,alpha);, i8 U/ `& L% l- k
p4=expcdf(A,mu);* h9 G$ f6 F0 \. Q4 U$ b0 J
[H4,s4]=kstest(A,[A,p4],alpha)
9 V, Y2 u# m' |; I; c' ?2 f" \% bif H4==0' U( R/ k# A D. S* K
disp('该数据源服从指数分布。') . _/ J' s3 }2 ~' \; s$ u' s
else& v& b u6 c7 ?; R. p3 U) C: D9 v
disp('该数据源不服从指数分布。') * @; P! e1 }; p# {: L' @
end( `3 T# F. V1 U& r5 A' R X& l
[phat, pci] = raylfit(A, alpha)
2 O3 |: _0 F: ep5=raylcdf(A,phat);
, C9 |" {: W+ h$ N[H5,s5]=kstest(A,[A,p5],alpha)4 S' ]9 s5 P! j6 N) H
if H5==07 {2 W- n1 F/ y; C( \! e; u
disp('该数据源服从rayleigh分布。') 8 r3 Q( J- d" w2 e9 ^) X G
else, \7 {; t/ g1 H: s1 _0 T4 s: [
disp('该数据源不服从rayleigh分布。')
4 t5 n/ x/ N: g& _end |
zan
|