- 在线时间
- 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)
: v9 Y5 I3 |/ C0 g( O9 qalpha=0.05;
$ v& n/ s2 G) G% g5 N' V% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
/ ~9 a6 y% f$ W5 ?: J- [A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
; U% s# K$ m; {[mu,sigma]=normfit(A);
- V; L c/ \8 b# S, v% @% Gp1=normcdf(A,mu,sigma);2 J6 w! O: L: `
[H1,s1]=kstest(A,[A,p1],alpha)
0 f9 @! `7 R. b$ `$ _! A4 d! en=length(A);' Z* Y% G: `! T) R* g. v
if H1==0( O, Y3 K5 i( q# z
disp('该数据源服从正态分布。')' [. G F! l7 o) X2 [; T5 b
else
' ?1 w3 z- U4 \: fdisp('该数据源不服从正态分布。')
. E" F: v& G; O& w& }end3 H4 y& D. u% _( b
phat=gamfit(A,alpha);
0 H3 X4 q4 W' B8 cp2=gamcdf(A,phat(1),phat(2));; `$ Y2 W9 o3 \, b
[H2,s2]=kstest(A,[A,p2],alpha)
/ j% i6 w# |9 w" o0 e2 G wif H2==0
% p9 G5 x2 @7 N; Q7 cdisp('该数据源服从γ分布。')
. o+ ]- S: |! ]else3 v0 u# N/ Z) F. m8 E8 V
disp('该数据源不服从γ分布。')4 q2 i7 l, J, I! Y8 M) J; O/ W
end: g% ~' B8 P" V) M4 [1 o
lamda=poissfit(A,alpha);
\. }9 n, ]# P' q- A! P$ gp3=poisscdf(A,lamda);6 p* O/ n* o3 Q7 L8 S5 n
[H3,s3]=kstest(A,[A,p3],alpha) H5 e" O; Y/ w' h9 ]4 R
if H3==0
1 t; `$ Y; f$ A' Ndisp('该数据源服从泊松分布。') , n, U) ^! t& {$ [9 h1 v
else" J$ }; C( y5 M5 j" D
disp('该数据源不服从泊松分布。')
n3 c; O. e g& L6 d( G. lend
1 T" Y; k- G8 Z) C* m8 _5 Imu=expfit(A,alpha);( p, J' u- |+ l* M$ G
p4=expcdf(A,mu);, O/ g6 W* w, _# @8 z8 H
[H4,s4]=kstest(A,[A,p4],alpha)
( p9 I5 W' S1 @: T) m- Aif H4==0
" z# ?! k5 z- O: v" o4 v# bdisp('该数据源服从指数分布。') " Y. a4 l2 d- y8 A3 e5 N; O
else
$ l, J8 {3 w8 W) {( J) E) ydisp('该数据源不服从指数分布。')
* B& {5 M( e, O. _" D, ]" send
( s7 \: q$ v Y' }# m8 X[phat, pci] = raylfit(A, alpha)
5 W3 x$ x: ?9 l: i( [# C5 Mp5=raylcdf(A,phat);
# z1 l% v6 J# ], a# n) U) h[H5,s5]=kstest(A,[A,p5],alpha)
1 i. i" t* w6 b5 Lif H5==01 J! q1 ?) ]: {# F& p" x" y
disp('该数据源服从rayleigh分布。') 9 w- i$ I' n6 q- o
else
: Y: e8 E+ o4 C% `+ O) V1 Fdisp('该数据源不服从rayleigh分布。')
8 P3 X, k+ B7 N ~/ C- a3 |end |
zan
|