对一组数据怎么进行分布拟合，求高手知道matlab程序

数学光头

function f=p_judge(A,alpha)
1 Z$ [$ z8 A* p) y0 a# V% m- N) palpha=0.05;
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
1 t% r1 x- _; d+ i[mu,sigma]=normfit(A);
p1=normcdf(A,mu,sigma);
) w! ^! E& u- \[H1,s1]=kstest(A,[A,p1],alpha)
) P4 h, L7 k4 P$ h: {  I- dn=length(A);
  `* V$ e7 ~' Q1 Tif H1==0
% Y# l. [  m# Zdisp('该数据源服从正态分布。')
else
3 R% r6 U$ c6 U/ P4 H* V5 udisp('该数据源不服从正态分布。')
; q' K9 k" T! p% {$ W# Q* send
5 X% I6 i1 O6 u2 X1 hphat=gamfit(A,alpha);
8 O( Y6 q2 d6 H- ~& W6 G% I1 z7 dp2=gamcdf(A,phat(1),phat(2));
  @7 E6 @' a* `[H2,s2]=kstest(A,[A,p2],alpha)
! H) O; q2 S6 v. {1 _if H2==0
disp('该数据源服从γ分布。')
else
  M. X& Y0 U8 S! G- ~8 Rdisp('该数据源不服从γ分布。')
end
' \6 ?. _" H4 W, Q& wlamda=poissfit(A,alpha);
p3=poisscdf(A,lamda);
- k+ U1 d: y1 t- w1 N7 s[H3,s3]=kstest(A,[A,p3],alpha)
9 J2 q5 D3 R4 U2 a- j% G' `/ ^if H3==0
% \4 z  [/ o! S6 a* D  Gdisp('该数据源服从泊松分布。')
else
1 r' V5 n! L$ \$ Sdisp('该数据源不服从泊松分布。')
end
mu=expfit(A,alpha);
& o' g- E% g; _# `2 Mp4=expcdf(A,mu);
7 q2 Q" O" e, o! r[H4,s4]=kstest(A,[A,p4],alpha)
: f" u+ O! b( t; z; j% F0 P- l& zif H4==0
9 Z: r1 c  g4 b% z% P; xdisp('该数据源服从指数分布。')
else
disp('该数据源不服从指数分布。')
end
[phat, pci] = raylfit(A, alpha)
/ D* ?7 ?$ r0 }4 M) `p5=raylcdf(A,phat);
[H5,s5]=kstest(A,[A,p5],alpha)
if H5==0
disp('该数据源服从rayleigh分布。')
else
disp('该数据源不服从rayleigh分布。')
: |' B% i- N7 Y  _) d& \end

lillianxiong

上面的是什么？

数学光头

lillianxiong 发表于 2012-8-20 17:16
8 v% V0 C7 v: T4 t上面的是什么？

你是说整个么，是代码

一场梦12

有点深奥 ，不错

小盆子。。。

感觉不是很懂得感觉