数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
7 |* T2 W5 P) ]
alpha=0.05;
! m2 u4 N( j7 D/ K, E
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
( P. [' p4 T3 I) n! [
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
) Z9 `/ u n* U4 X3 b8 V
[mu,sigma]=normfit(A);
1 Y- ]5 a- J& c+ |
p1=normcdf(A,mu,sigma);
, O" i7 {6 L6 E1 T
[H1,s1]=kstest(A,[A,p1],alpha)
3 l9 s: R3 Z; U, k9 }% Z; U
n=length(A);
; n$ c- Y. G2 \0 R y9 |9 S3 V9 D2 v
if H1==0
1 b. q8 [% ]9 g) A8 c* j
disp('该数据源服从正态分布。')
! `0 A/ e, w$ e* R! Z* p' u
else
3 v. u1 ^% r' {2 ^8 H3 h2 Y
disp('该数据源不服从正态分布。')
- L1 b' s* D9 h9 x# Z% l
end
2 h0 Q2 ]: O; e; x! \
phat=gamfit(A,alpha);
8 l/ A4 R, `7 g" R
p2=gamcdf(A,phat(1),phat(2));
3 i4 W& ?1 I) E2 A' F
[H2,s2]=kstest(A,[A,p2],alpha)
. L% I1 n0 W: z5 n% a
if H2==0
; J/ f- V( v8 I' ]8 M
disp('该数据源服从γ分布。')
8 \- g y, M% w4 m& P
else
i$ c9 B# Z+ D; P# ]2 v
disp('该数据源不服从γ分布。')
/ F6 r# U: h, R+ o" U1 r! D" ^( q4 Q9 ]7 R
end
+ ~: K* F* E$ d# ?% W! r
lamda=poissfit(A,alpha);
7 a" V& ~8 W0 {/ f. G+ `" Z, m
p3=poisscdf(A,lamda);
: Y2 q& T0 D( m' ?# t8 J
[H3,s3]=kstest(A,[A,p3],alpha)
" q$ ]6 m& s$ V: I9 d! R: i
if H3==0
/ V4 i! B1 e7 O: _+ Z- h
disp('该数据源服从泊松分布。')
, I; ^' J n/ L7 H" Z. `, ~4 R7 H
else
& c% R; W5 q. k$ ^7 [* c. n' y
disp('该数据源不服从泊松分布。')
. B! Z6 g# M# u9 O
end
, v6 \+ @& v- b& D
mu=expfit(A,alpha);
7 }5 X. X9 `. c- g6 @1 S
p4=expcdf(A,mu);
$ D1 ?) E$ D% l
[H4,s4]=kstest(A,[A,p4],alpha)
/ ]: W% n0 V" A) @
if H4==0
. ]7 ]9 v& p- V
disp('该数据源服从指数分布。')
$ E R( z; L$ O5 k" N4 e# {% ]5 i
else
9 k$ u k, B* l
disp('该数据源不服从指数分布。')
/ K. G. S+ \2 o. M* j2 L0 q# L
end
8 Z8 O. v" h; K4 G" W8 \
[phat, pci] = raylfit(A, alpha)
u# X3 `( A% K/ Z
p5=raylcdf(A,phat);
+ Y8 a8 m; }% Z& _$ T8 F7 N ~
[H5,s5]=kstest(A,[A,p5],alpha)
0 L" u t, J/ ]& }- E! _! c2 H' O
if H5==0
8 e4 G% i" ]. s. D+ j
disp('该数据源服从rayleigh分布。')
4 P) u, T& N% u) s
else
9 }" L% V. H+ ^5 ?/ x3 r
disp('该数据源不服从rayleigh分布。')
+ N! G! G! \9 l7 N
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
$ w% J# m; h8 |5 s8 y1 a# p
上面的是什么？
$ M" G9 ~7 m* \4 c. O4 }& I
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
