数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
+ W8 y/ a. h: K7 g% \& e, c) H
alpha=0.05;
7 f6 p9 F" N! _
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
' u3 w8 ^9 \# U7 o5 q+ m
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
& e s/ }8 F6 N5 }
[mu,sigma]=normfit(A);
7 t6 |& W3 z, A+ W. H. h U/ K
p1=normcdf(A,mu,sigma);
2 l! T) v8 ~4 G8 T8 z9 @+ d
[H1,s1]=kstest(A,[A,p1],alpha)
0 q! O- a8 G8 U7 b
n=length(A);
- M y( E5 v4 c+ c9 M
if H1==0
3 N- p$ h' U; W! F* ]
disp('该数据源服从正态分布。')
4 s# I5 Q, u0 X, m, W" n
else
, b2 R3 k$ |6 |$ j" t$ ~1 I }2 k
disp('该数据源不服从正态分布。')
' N. t% d5 }* H" n' s& T
end
% Z) l; y2 V* k3 q6 @
phat=gamfit(A,alpha);
+ P+ d# v7 G# u& ?/ s$ v4 G
p2=gamcdf(A,phat(1),phat(2));
" f2 E4 ^" m/ h% L. Y4 }4 b/ b: y8 S
[H2,s2]=kstest(A,[A,p2],alpha)
$ V7 P2 s4 z. g# U& d2 h
if H2==0
* Z7 y& \3 J! V( G) b
disp('该数据源服从γ分布。')
0 R- T3 y7 M }- [9 T
else
( Q1 u( ?& h" s- j( g
disp('该数据源不服从γ分布。')
: c, a4 C! a) i1 ]" d5 y: z
end
$ Q5 y7 _6 Y( q2 Q. N
lamda=poissfit(A,alpha);
# v2 g5 w! S4 S/ J5 R* U
p3=poisscdf(A,lamda);
3 C4 _+ r4 j' r6 g, p) E5 L
[H3,s3]=kstest(A,[A,p3],alpha)
' r Y) O$ c N+ O: Z* g
if H3==0
# {0 B3 r4 e( G9 i+ J* @# h/ b8 S
disp('该数据源服从泊松分布。')
{" i8 R J! } B
else
" ~. @' i$ N3 `# u5 W: d
disp('该数据源不服从泊松分布。')
0 d9 Z+ W! ~& |) u/ Z
end
; Z. [+ \) O$ b8 z7 o
mu=expfit(A,alpha);
$ J$ B( o. v/ U5 H/ J
p4=expcdf(A,mu);
; E1 Q3 h# |* _! P2 f$ B
[H4,s4]=kstest(A,[A,p4],alpha)
6 `* P/ p2 h1 }0 m1 k( @1 y( _ e% ]' n
if H4==0
L; r6 Q U# Z4 N8 Y7 t3 l
disp('该数据源服从指数分布。')
8 F6 e% f) V2 V8 o, e
else
9 O, M, u3 x. d. F* n- q6 e
disp('该数据源不服从指数分布。')
# P1 U3 E& x3 f" a
end
9 S s, n% T3 h" ?
[phat, pci] = raylfit(A, alpha)
+ h, H( }( v# K- p- V( M# v' b
p5=raylcdf(A,phat);
- G) e2 L) D9 U/ j( O
[H5,s5]=kstest(A,[A,p5],alpha)
% T( ]& @0 a4 ^2 S# O
if H5==0
& Q: P5 `2 v0 f# U3 ~3 K* i
disp('该数据源服从rayleigh分布。')
5 d# j* u4 J' x( m
else
P% Y9 j( f. P. K2 ?
disp('该数据源不服从rayleigh分布。')
6 I" B' N4 }3 Y; b& V
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
! [9 g( @" ?5 ?& W5 j
上面的是什么？
2 X/ R. n2 _+ ^9 l. \
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
