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

数学光头

function f=p_judge(A,alpha)
8 o3 d4 e9 `, R! jalpha=0.05;
  H- Q; x2 k6 ?; D$ I% U9 j8 h% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
9 U' t! D: E. P[mu,sigma]=normfit(A);
- x% B% a) ^- V' R6 X& W" Jp1=normcdf(A,mu,sigma);
[H1,s1]=kstest(A,[A,p1],alpha)
1 O" r0 Z4 ~# \  v, P3 `n=length(A);
3 r3 z) O$ U+ p8 Tif H1==0
disp('该数据源服从正态分布。')
+ p# `) B8 R1 _9 melse
disp('该数据源不服从正态分布。')
end
( I, H4 N+ \6 \( N0 x$ C) L4 a  iphat=gamfit(A,alpha);
# @1 y5 c: K1 F. W- Lp2=gamcdf(A,phat(1),phat(2));
; p& A9 W/ ]5 }8 s. P4 ?) M[H2,s2]=kstest(A,[A,p2],alpha)
7 O9 K9 F4 ?3 ^- i1 m" ^# ~if H2==0
- t# d6 y' j: P4 vdisp('该数据源服从γ分布。')
else
& B& }6 b( K7 W1 y# x- v; W* O7 {disp('该数据源不服从γ分布。')
end
/ J$ Q. e3 p- W/ F3 g. t. k9 R6 tlamda=poissfit(A,alpha);
p3=poisscdf(A,lamda);
5 j. e! D9 J6 k' K" ^8 N! P: O[H3,s3]=kstest(A,[A,p3],alpha)
" s4 u+ O$ g" j' ~if H3==0
, P9 k9 ~/ N$ ~6 U6 `; q; Xdisp('该数据源服从泊松分布。')
else
disp('该数据源不服从泊松分布。')
end
1 N& `- r2 t4 X' F' Kmu=expfit(A,alpha);
p4=expcdf(A,mu);
[H4,s4]=kstest(A,[A,p4],alpha)
if H4==0
disp('该数据源服从指数分布。')
$ ^1 K+ `7 |6 b# u9 Z7 ]else
" @7 |1 u5 B  P# w% gdisp('该数据源不服从指数分布。')
1 E3 |- w  m$ }3 n# I8 c: G5 ?end
" }9 ~, ~+ D% Q# t' A7 Q[phat, pci] = raylfit(A, alpha)
: L$ l3 e6 t7 w6 n8 M! Fp5=raylcdf(A,phat);
[H5,s5]=kstest(A,[A,p5],alpha)
+ B# ~& ]) V2 I1 ~7 m! f$ q. Hif H5==0
1 K) Y+ n* N* c/ M" u1 n/ z6 @disp('该数据源服从rayleigh分布。')
4 x0 S; l) c3 K5 v0 welse
5 s4 ~' m! r. y7 D  hdisp('该数据源不服从rayleigh分布。')
2 V1 \( L/ w9 @! S) p* u( Nend

小盆子。。。

感觉不是很懂得感觉

一场梦12

有点深奥 ，不错

数学光头

lillianxiong 发表于 2012-8-20 17:16
1 b. H" [1 I' M" d; ?( K+ g  {0 Q9 l上面的是什么？

: H2 o+ R$ m6 `' V你是说整个么，是代码

lillianxiong

上面的是什么？