数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
6 p7 Q! {. G& _6 N
alpha=0.05;
2 Y4 h b. `$ f7 X, X
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
* j, h2 N: a& v! O. y, I8 V. L
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
) A8 C ~ f6 n$ g% T' d, I
[mu,sigma]=normfit(A);
" O! C# f6 C n9 ^% z% Z
p1=normcdf(A,mu,sigma);
- Q& L( `5 b" l) W. P! s
[H1,s1]=kstest(A,[A,p1],alpha)
4 F, ]- n) b" ?8 `4 Z" ]
n=length(A);
9 S1 \8 @0 c/ P- O% [1 d% l
if H1==0
6 @; a/ M7 B% J* e2 G5 K
disp('该数据源服从正态分布。')
' ?( x0 |9 h+ C4 \% r$ a4 {9 l3 ^
else
% C) c- U! H+ O, S# l- N/ J4 _
disp('该数据源不服从正态分布。')
+ Z* ^2 P* m7 Y' y( H8 ~
end
2 f0 E' `7 l) q% C$ [9 K ]
phat=gamfit(A,alpha);
/ E$ y/ e' F3 `3 q
p2=gamcdf(A,phat(1),phat(2));
& m; N1 i! q3 ~; }+ J" g) c- o& D
[H2,s2]=kstest(A,[A,p2],alpha)
8 l. }4 U9 ^/ s1 C+ R& f
if H2==0
) G" J- W$ Y. S" K. j1 Z
disp('该数据源服从γ分布。')
1 e8 t% }" N: X9 H
else
. _& I& w ^/ E/ V. I/ s
disp('该数据源不服从γ分布。')
- S [! D L; f% G
end
P# X9 k: w5 C. Q
lamda=poissfit(A,alpha);
! U; F$ F( @9 [6 R% G( g2 c f
p3=poisscdf(A,lamda);
9 X% Z) a0 n) T' g3 s( b- W) X
[H3,s3]=kstest(A,[A,p3],alpha)
+ r# Y+ a; x: F \8 x
if H3==0
6 W4 z9 y6 q; J
disp('该数据源服从泊松分布。')
Q4 U3 V* D4 X. N/ W3 \
else
. s! U; I5 l% v& _
disp('该数据源不服从泊松分布。')
% A4 l* ~: ?5 H1 p" k1 [# w9 ~9 ?
end
7 @, {" L: P6 C5 }
mu=expfit(A,alpha);
7 d- m+ |) f% \8 X
p4=expcdf(A,mu);
; z9 C; _# @7 a
[H4,s4]=kstest(A,[A,p4],alpha)
* o$ R8 S9 T1 s
if H4==0
V: N8 q, ]6 t& e# p
disp('该数据源服从指数分布。')
: i# x: } b5 a/ M2 E
else
O% m) F2 M; |6 }( O! J
disp('该数据源不服从指数分布。')
; w! u& Z( f4 C: P, C# `
end
: ~" L3 p3 N1 I' o# G
[phat, pci] = raylfit(A, alpha)
3 i+ v# L* e7 b" ^/ k% H" z) ?0 I
p5=raylcdf(A,phat);
, t. Y6 ]' B& k: }7 q4 t
[H5,s5]=kstest(A,[A,p5],alpha)
# N6 I9 q$ p& c2 t
if H5==0
& x: x. t* K! g" z+ \. C
disp('该数据源服从rayleigh分布。')
; L: M5 K/ |( q9 z* R, N0 h7 z
else
7 T& }8 s! z1 T1 v( e
disp('该数据源不服从rayleigh分布。')
6 f5 }1 A0 t% R. A$ r! t
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
7 }! ^1 n& ^2 c% E h7 |8 W& G
上面的是什么？
: @# P3 j* z W7 i7 v
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
7 }! ^1 n& ^2 c% E h7 |8 W& G