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

数学光头

function f=p_judge(A,alpha)
& o0 L9 v2 U+ S, l) r, z6 Ialpha=0.05;
" \4 `! ]% D) v: P: C5 |$ C' [% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
[mu,sigma]=normfit(A);
# N5 q. o( G$ F. |" A3 X( N4 N. w7 rp1=normcdf(A,mu,sigma);
[H1,s1]=kstest(A,[A,p1],alpha)
n=length(A);
if H1==0
disp('该数据源服从正态分布。')
; W' }( A% B' v7 D/ I: xelse
6 L1 @* [  _9 S. ]disp('该数据源不服从正态分布。')
end
phat=gamfit(A,alpha);
8 }9 O1 k4 K0 ?' d2 Y  up2=gamcdf(A,phat(1),phat(2));
4 ?4 U( B$ g) K5 K4 b# r[H2,s2]=kstest(A,[A,p2],alpha)
if H2==0
disp('该数据源服从γ分布。')
% P& v# T9 v# w1 R  K2 v) p- nelse
disp('该数据源不服从γ分布。')
# J. w: E# ?0 ^- P3 M0 O' `end
8 s4 ~- {+ R" V  |6 c5 dlamda=poissfit(A,alpha);
% T5 J9 }! K# R' p0 w' Y5 ep3=poisscdf(A,lamda);
[H3,s3]=kstest(A,[A,p3],alpha)
: q4 A0 V2 ~6 L. Lif H3==0
disp('该数据源服从泊松分布。')
- D) T/ W% H- j  lelse
, s: j$ n/ H3 L8 n5 V" ydisp('该数据源不服从泊松分布。')
end
mu=expfit(A,alpha);
p4=expcdf(A,mu);
2 j& X  N! X/ Y8 ~" [% R[H4,s4]=kstest(A,[A,p4],alpha)
if H4==0
disp('该数据源服从指数分布。')
$ o9 w& `& S2 \1 ~2 D- _) h: eelse
4 c" i* |* @  s- ydisp('该数据源不服从指数分布。')
3 x( Y+ v# T1 o9 X! Q7 }+ X- mend
2 P+ `' K6 N6 x( I[phat, pci] = raylfit(A, alpha)
p5=raylcdf(A,phat);
1 b+ f" y- T: a[H5,s5]=kstest(A,[A,p5],alpha)
# a4 D3 E0 |5 C. z; vif H5==0
disp('该数据源服从rayleigh分布。')
else
disp('该数据源不服从rayleigh分布。')
9 ^' R' ?5 a: e4 H* |end

lillianxiong

上面的是什么？

数学光头

lillianxiong 发表于 2012-8-20 17:16
* ?$ D  |# _8 E7 H1 w$ L" ]( ~上面的是什么？

! P; W6 `4 w8 Y! S% S你是说整个么，是代码

一场梦12

有点深奥 ，不错

小盆子。。。

感觉不是很懂得感觉