- 在线时间
- 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)
8 n z- p: I% M! ~2 ^3 zalpha=0.05;/ W. O, m' h2 K9 n& ?( f
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。: `0 X1 ^8 i; V5 u) n" a0 G
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
$ V7 r7 A3 _# Z A, ^6 f* \[mu,sigma]=normfit(A);3 u6 D% s$ ?! u1 b6 B4 N) w" p: P7 ?* Z
p1=normcdf(A,mu,sigma);$ D1 C. M, ^" h6 X6 n. ~8 o
[H1,s1]=kstest(A,[A,p1],alpha)' l6 \9 V3 x$ e4 i6 \3 y3 ^ ?+ B
n=length(A);* [* c" e* E+ w6 q+ a n3 r5 J
if H1==0% X3 u& m! D. X. A6 M4 C
disp('该数据源服从正态分布。')
$ n1 V; q- ?+ s- I [. [; Qelse
1 I2 s) o0 \5 q$ m- I+ F& y3 L7 f' hdisp('该数据源不服从正态分布。') 6 |7 e3 T: [, l- W3 N
end
( ^+ X* n- q4 s+ M" hphat=gamfit(A,alpha);
' t4 ~6 M% L" ^5 ?4 rp2=gamcdf(A,phat(1),phat(2));% \- A5 q7 W5 P8 I/ q, R4 }
[H2,s2]=kstest(A,[A,p2],alpha)
5 _) p% Z! M% n n& fif H2==0
m1 C! v$ [9 N' K+ fdisp('该数据源服从γ分布。')* c; y- v+ g; Q. q6 H
else
$ @' J3 ^% c5 D. g* s; c, xdisp('该数据源不服从γ分布。')
* p& }8 a a7 G* _ [end
! X4 @( a9 F8 B$ A& W- D& a; \lamda=poissfit(A,alpha);
: ]" s$ @! A1 o% y5 R8 x, sp3=poisscdf(A,lamda);
; u4 T, o: i- h8 X* B8 Z" e[H3,s3]=kstest(A,[A,p3],alpha)4 U- J% y1 N" |
if H3==0
( L9 c) ]1 _9 hdisp('该数据源服从泊松分布。') # {7 d7 m6 ^1 ]
else
8 k% V0 s, Q3 T7 \! ~9 K4 Idisp('该数据源不服从泊松分布。')
+ S# ]. e3 X0 iend
7 D" J q4 v( F2 r1 W* ~' U) Rmu=expfit(A,alpha);
' O( C$ w' D; Z) ]3 vp4=expcdf(A,mu);$ g8 U: O; k4 B) U9 |' _/ ^
[H4,s4]=kstest(A,[A,p4],alpha)
% R5 E, I0 n: \6 Eif H4==08 Y! T/ L9 S& Z- O
disp('该数据源服从指数分布。')
- x4 n' L4 y. W) N, Zelse4 Z2 ^+ e0 s+ I$ y
disp('该数据源不服从指数分布。') % r3 S. a$ P1 a1 o
end! i# V, S# x1 [# `
[phat, pci] = raylfit(A, alpha)
$ q* o# o+ p Fp5=raylcdf(A,phat);
% K) A" J: f& n% b[H5,s5]=kstest(A,[A,p5],alpha)5 ~; h" c, f8 k/ n: q; A) S+ a1 t- e' m
if H5==0
2 N* Z$ \" D7 }4 A8 `: Z* t; ]disp('该数据源服从rayleigh分布。')
5 \2 Z" l# N# d& ]else' I/ w/ s+ m- p
disp('该数据源不服从rayleigh分布。')# h' s u, G" x# u
end |
zan
|