- 在线时间
- 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), u2 U' g# }9 U9 k- f
alpha=0.05;
2 O( r g e' w* _7 q+ R' j1 R2 d" D9 R% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
2 I4 h9 X3 o R4 k5 n1 QA=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
: M! a% z1 ?% h[mu,sigma]=normfit(A);. c- v! u6 x, w M+ Q X% \! e0 w
p1=normcdf(A,mu,sigma);, q) L- L7 X; W) H
[H1,s1]=kstest(A,[A,p1],alpha)( x W. O8 X; Y% z
n=length(A);! `/ L% w4 P& m7 f$ W7 r5 u: X3 k
if H1==0* j' `9 j- P3 x- ^# a( l
disp('该数据源服从正态分布。')
6 ~4 d- y' G2 \( aelse+ u) t, V2 G" E W+ a% A; @* j
disp('该数据源不服从正态分布。')
9 Q" M4 T' x4 c2 X i* u+ `end
2 I" C1 T5 F2 l6 ~ h! f, iphat=gamfit(A,alpha);
" l1 I _6 k! L9 J$ H- ip2=gamcdf(A,phat(1),phat(2));& I. r2 e+ o E4 B! f2 P' f
[H2,s2]=kstest(A,[A,p2],alpha)0 M9 ^: y& ?1 Q
if H2==0+ S1 j" Q. q5 @
disp('该数据源服从γ分布。'); d' C4 I, ^" }4 { r H* `
else/ J" j( ^; S3 n
disp('该数据源不服从γ分布。')4 d$ }4 J! c" y" Z# m
end1 c2 |0 F: S! E: k) T# R: Q, ?
lamda=poissfit(A,alpha);
" z1 ?, F7 U+ ?! Ip3=poisscdf(A,lamda);
2 D1 K+ s9 `. Y- D/ w! g: z4 K[H3,s3]=kstest(A,[A,p3],alpha)
$ x; A" E7 n7 i8 P- E& A$ Yif H3==0
$ `; f6 u/ R8 w. I% m& H* p) L+ ldisp('该数据源服从泊松分布。')
( Y/ @' d! V- velse
% c5 G; \4 y6 s( k6 Idisp('该数据源不服从泊松分布。')
, w, X( B4 G4 @0 H* Q7 L7 ^; }end* l* y+ o# _2 S) H8 J
mu=expfit(A,alpha);
- C* t& u2 ]# ?- ~) mp4=expcdf(A,mu);
/ s: w( C6 S* o7 \; N; B6 Z[H4,s4]=kstest(A,[A,p4],alpha)
/ s) o N" P! B5 \if H4==0 ~$ C+ {' l0 j/ e E+ }
disp('该数据源服从指数分布。') 7 x7 v* Z* e3 W
else
0 @' c: a2 I, n7 _disp('该数据源不服从指数分布。') 5 U& B2 y% p: ~
end
, {7 [. j9 t0 _3 |5 L[phat, pci] = raylfit(A, alpha)/ |' e0 r* d) r; v8 J& `
p5=raylcdf(A,phat);7 S+ ~; a: }2 |& ~" l) C
[H5,s5]=kstest(A,[A,p5],alpha)
: A% ^ c R ^1 Jif H5==04 }2 k4 M+ S7 Z! b* G
disp('该数据源服从rayleigh分布。') 3 d( g8 V) F P2 e9 m
else
! |0 p( c! C9 r" Mdisp('该数据源不服从rayleigh分布。')
4 E, T/ \, t0 i* pend |
zan
|