数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
1 c% |( M @1 a
alpha=0.05;
; `4 S7 L8 X% ^8 Z' D: ^0 y. t
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
: i( H& Y0 ?2 X$ s7 q8 e% D6 s
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
& @0 s# ]3 Z0 I7 D: [
[mu,sigma]=normfit(A);
0 |, [" t/ I: |! Z! d+ `
p1=normcdf(A,mu,sigma);
* {0 S7 z" W# B$ {9 \( ` Q
[H1,s1]=kstest(A,[A,p1],alpha)
) r/ _/ n0 j0 ~+ y/ L2 n& P
n=length(A);
7 ` b6 n$ J4 p% M
if H1==0
& `# P: @$ ]1 b* @* j* A/ f7 a6 Z! l
disp('该数据源服从正态分布。')
+ N1 X* U/ k/ j* Y9 p
else
/ @, C" H- R7 A3 R
disp('该数据源不服从正态分布。')
: a4 ]* ^5 y4 f: ~9 q
end
2 H) Y" N1 L. H" ~; R4 B- `
phat=gamfit(A,alpha);
( m! S& z' M Z0 g, u! t: k% L
p2=gamcdf(A,phat(1),phat(2));
$ T d/ Z* [% j0 F Q6 {, Z4 w
[H2,s2]=kstest(A,[A,p2],alpha)
- V" N" a+ E/ U( L ~
if H2==0
U- C$ E3 J8 J( @) x
disp('该数据源服从γ分布。')
7 j6 @) I1 @! r; k
else
: w" ~5 ~0 X5 q$ ?2 _; w
disp('该数据源不服从γ分布。')
! N( N6 i% |9 d }; D
end
# \$ u: V2 o( p" v+ H+ u" Y
lamda=poissfit(A,alpha);
* o' a/ ?* z7 Y! T
p3=poisscdf(A,lamda);
& x; n8 g) D+ K, F& H' _
[H3,s3]=kstest(A,[A,p3],alpha)
0 @( \# Z$ X' D) j( D7 F
if H3==0
) s' K5 F' c6 {4 O' }: U, t: E
disp('该数据源服从泊松分布。')
+ Y- l. U8 u/ ?+ I& W! W* F
else
' a( E" b' g% n! {" I* o) r9 I
disp('该数据源不服从泊松分布。')
3 t8 j7 B2 b* Z* _ n
end
# T1 ?) ]' C- B/ R3 X
mu=expfit(A,alpha);
& `! ]% b& E3 _# F9 Q
p4=expcdf(A,mu);
+ M' L. e/ i% [6 \# E1 g! o5 b
[H4,s4]=kstest(A,[A,p4],alpha)
5 I q b0 h+ B- q
if H4==0
$ k1 W8 M2 L% U; n
disp('该数据源服从指数分布。')
& M9 u- s/ C+ O& |
else
) h3 L0 G; x9 i/ j: P
disp('该数据源不服从指数分布。')
$ W/ N3 `- I) P3 {4 ?' r; h
end
1 ~7 e c+ R% C" Z7 |3 ^/ q% U+ j
[phat, pci] = raylfit(A, alpha)
k! N8 m" w, Y$ S' E
p5=raylcdf(A,phat);
9 h S6 f7 f( y' w* _4 L9 j
[H5,s5]=kstest(A,[A,p5],alpha)
! [# V! T# R6 G& D; h
if H5==0
9 V$ J( v# C+ ~9 I w+ F
disp('该数据源服从rayleigh分布。')
" ^5 }, ?- A; t/ p# h
else
; \4 @1 [4 [ D' R8 o8 z& e
disp('该数据源不服从rayleigh分布。')
% _; @) j& y8 k* M) H
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
4 {. f; d$ D. |9 {% a8 N8 o7 Y
上面的是什么？
' }9 z! @2 M Q: j8 K' U6 ?
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
4 {. f; d$ D. |9 {% a8 N8 o7 Y