数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
) H! Y( J( {+ e/ f. a0 K
alpha=0.05;
7 H/ i8 n4 P+ N4 H5 c
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
1 c8 T7 U% K: n' E
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
5 j: S0 r/ M$ }* s" {4 `
[mu,sigma]=normfit(A);
S6 W. ^: _& M2 f# [
p1=normcdf(A,mu,sigma);
/ w3 \8 {; [7 a1 b. t
[H1,s1]=kstest(A,[A,p1],alpha)
4 w# {" K2 M* a1 v
n=length(A);
! `4 N( B9 v( i' F+ z& [" k! e
if H1==0
Z, G2 ]0 ^0 P2 x" B Z0 K9 h
disp('该数据源服从正态分布。')
6 Z; N& f2 K# [$ W. n
else
5 i' q/ ?% e( }
disp('该数据源不服从正态分布。')
1 h: N0 g+ W1 }& o
end
- S. B2 g% t5 e
phat=gamfit(A,alpha);
7 U1 r4 c( O% M, o
p2=gamcdf(A,phat(1),phat(2));
" g* B) i% s+ Q9 m/ l1 r F7 X( _
[H2,s2]=kstest(A,[A,p2],alpha)
* q/ _) ^+ U) ^2 u& p1 N
if H2==0
; I1 t3 |$ s. c7 r' v( ]
disp('该数据源服从γ分布。')
& ]$ T: @: B7 P2 D% i
else
9 C8 G% r; f$ B; s
disp('该数据源不服从γ分布。')
4 i( ?4 [; X2 N9 S. [- _
end
3 q, f9 \2 |; f. E
lamda=poissfit(A,alpha);
6 R: _3 h0 |& e
p3=poisscdf(A,lamda);
& O" E; A% g/ _( A& a' Q
[H3,s3]=kstest(A,[A,p3],alpha)
6 F9 X$ |3 ?4 ~: ^& h
if H3==0
8 S4 O, z4 h4 C: I* }2 Z
disp('该数据源服从泊松分布。')
! B7 @' o2 s5 e- A3 F1 s, J
else
( \/ M0 M# R% I3 t& i5 d) |
disp('该数据源不服从泊松分布。')
: ^9 ^: @ p. d8 ^, k% I4 t8 o6 H
end
) w6 T" |" J, ?+ D# Y+ ~
mu=expfit(A,alpha);
L9 S2 j' W# r
p4=expcdf(A,mu);
% |' D5 E5 } ?5 R, V2 u' w# `
[H4,s4]=kstest(A,[A,p4],alpha)
& t8 p( }* b* e8 }
if H4==0
+ o \) r. } [# F
disp('该数据源服从指数分布。')
8 `4 z9 h+ E) v- ]) L
else
1 L% _5 C' b8 z8 s4 T% ~
disp('该数据源不服从指数分布。')
% R0 m0 r' z& X! ~8 Y" J" f* T
end
. u/ ~- Z2 p* U" _# o' h* k* y4 d
[phat, pci] = raylfit(A, alpha)
. _$ s. X6 e7 q, @" R4 h3 v
p5=raylcdf(A,phat);
0 {, ^. U. R" j3 J; V; `9 X) X
[H5,s5]=kstest(A,[A,p5],alpha)
' d% \9 p2 r/ p9 r
if H5==0
& W/ P. ]6 I. l0 M r% s* _! e
disp('该数据源服从rayleigh分布。')
4 o/ w) U- X6 i+ M7 f2 k; T3 o
else
( B7 ^, V9 g; n6 c* G
disp('该数据源不服从rayleigh分布。')
( K+ y& J1 S. m) I% ~: {
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
* z8 w3 [- q/ T
上面的是什么？
# _1 p- D( L* z1 {4 _' Y
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
* z8 w3 [- q/ T