数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
: V; M2 m" a& K) k9 @0 \3 b) o- l
alpha=0.05;
1 }) d% T7 z k9 [
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
$ e( h& J& F0 \1 R6 |" i6 k" Q
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
* R3 v9 }. b. ?
[mu,sigma]=normfit(A);
7 c" X# ?5 M. U! h6 Z* u
p1=normcdf(A,mu,sigma);
; Q) f" C) `% M" ]) ~7 i& E
[H1,s1]=kstest(A,[A,p1],alpha)
$ u2 l: \: W! c2 C( b/ b( g
n=length(A);
$ ]% y. c. o& I* W0 W. n
if H1==0
$ {$ {" N- H4 [8 p) N" a
disp('该数据源服从正态分布。')
0 I' o* C: m# A, t. c
else
6 e- c& _- }# K0 ~" a$ f# t
disp('该数据源不服从正态分布。')
/ ] a# I$ g3 g' |6 a
end
* v0 G& L1 o3 c
phat=gamfit(A,alpha);
# l. {+ j! w+ M( |0 O9 ^. C9 A- S
p2=gamcdf(A,phat(1),phat(2));
" ~2 K G3 H: p3 Q8 J0 `: F- ~
[H2,s2]=kstest(A,[A,p2],alpha)
& q! p5 ^1 }' U5 I/ W. r$ l
if H2==0
2 ]% ~# N3 T; W% s: G
disp('该数据源服从γ分布。')
0 W! Z @, e" F3 h y
else
- |' G% ~. s! t
disp('该数据源不服从γ分布。')
$ [7 l K" Y+ q E+ o: @: y. I" p
end
6 X) {+ U9 d( Z+ l* f; D6 E8 I
lamda=poissfit(A,alpha);
# b4 E: ^3 b( V ~: P6 @1 u
p3=poisscdf(A,lamda);
* n0 z( c6 V" C7 o
[H3,s3]=kstest(A,[A,p3],alpha)
4 ^) Q! ]3 E a8 ~1 \$ F9 ^
if H3==0
: ~7 J. P/ N9 `; ]: ~
disp('该数据源服从泊松分布。')
; \& S! h8 I. `/ g1 g+ L
else
: L/ g" |7 z7 t' P
disp('该数据源不服从泊松分布。')
0 ^) R& R$ i# C C. Q1 @
end
! c5 b- F$ K, v
mu=expfit(A,alpha);
; b" I# _( X; Y$ l6 T2 `
p4=expcdf(A,mu);
7 t# g& |, y( S! C
[H4,s4]=kstest(A,[A,p4],alpha)
( _( d! O8 [8 [
if H4==0
9 e+ j5 }" c6 M. o; N
disp('该数据源服从指数分布。')
" M" J8 w8 p% F* r5 [( l. o }; S
else
6 S( [" x$ R- M+ k. d* E
disp('该数据源不服从指数分布。')
/ ^3 a* Z4 Z; e# |( t
end
9 y8 o6 j4 @5 M' M" L
[phat, pci] = raylfit(A, alpha)
+ w: e" c. b% H" r
p5=raylcdf(A,phat);
+ z5 ]0 T+ k: j, U6 M
[H5,s5]=kstest(A,[A,p5],alpha)
; u K& ~3 J% U
if H5==0
6 a$ C8 x0 ] n# N/ \7 K E5 t$ ^) U
disp('该数据源服从rayleigh分布。')
+ n6 b. |2 B7 [( b! w
else
; H$ b% r7 E9 A( @$ J
disp('该数据源不服从rayleigh分布。')
4 Q) d; y7 R' V& P
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
5 j. d1 F7 @$ @( R9 c
上面的是什么？
' Y4 j: Y$ j/ Q h0 v0 o* i% v! J/ }
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
5 j. d1 F7 @$ @( R9 c