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

数学光头

function f=p_judge(A,alpha)
3 S4 h, u$ `, I' b0 J5 w* Kalpha=0.05;
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
- r7 |0 h7 S0 ^! ]) i: _A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
8 i9 d8 S5 Z5 j# |! Y, E0 ^[mu,sigma]=normfit(A);
p1=normcdf(A,mu,sigma);
[H1,s1]=kstest(A,[A,p1],alpha)
/ O: Q# q3 m. a7 z, @n=length(A);
if H1==0
1 z0 G9 _( J/ R6 @. N6 T" Adisp('该数据源服从正态分布。')
; _* m8 ]# H& kelse
disp('该数据源不服从正态分布。')
" K9 Z/ F' r2 F6 \$ I( Dend
phat=gamfit(A,alpha);
p2=gamcdf(A,phat(1),phat(2));
[H2,s2]=kstest(A,[A,p2],alpha)
! G" s+ ^1 G& R! Aif H2==0
disp('该数据源服从γ分布。')
else
5 h/ N. P6 p; m6 o" Wdisp('该数据源不服从γ分布。')
3 b9 N) k8 P+ j' ~" w' Rend
lamda=poissfit(A,alpha);
p3=poisscdf(A,lamda);
$ Z" q# Y$ d4 A4 b) x1 i5 H1 @  H[H3,s3]=kstest(A,[A,p3],alpha)
) z$ Y+ A& o5 c; d! uif H3==0
$ l  A' L/ t# m; Kdisp('该数据源服从泊松分布。')
else
1 ~: D) v$ }, c4 n! W7 V! }disp('该数据源不服从泊松分布。')
end
mu=expfit(A,alpha);
p4=expcdf(A,mu);
+ ~9 d) c; F+ G. f0 }# }  M" ^[H4,s4]=kstest(A,[A,p4],alpha)
if H4==0
disp('该数据源服从指数分布。')
else
$ Z9 m. A; E* ^disp('该数据源不服从指数分布。')
% V: h4 W8 p5 W6 ]3 z( h" Wend
7 U5 R) q1 `, P- v4 L/ p% I, \7 i[phat, pci] = raylfit(A, alpha)
# F% F: g# a0 K3 p( r( ap5=raylcdf(A,phat);
[H5,s5]=kstest(A,[A,p5],alpha)
if H5==0
* Z- h# v- k5 o8 j% Y; D1 Rdisp('该数据源服从rayleigh分布。')
else
disp('该数据源不服从rayleigh分布。')
end

小盆子。。。

感觉不是很懂得感觉

一场梦12

有点深奥 ，不错

数学光头

lillianxiong 发表于 2012-8-20 17:16
上面的是什么？

你是说整个么，是代码

lillianxiong

上面的是什么？