数学建模社区-数学中国
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
[打印本页]
作者:
数学光头
时间:
2012-8-20 16:34
标题:
对一组数据怎么进行分布拟合，求高手知道matlab程序
function f=p_judge(A,alpha)
" ]. G: n) W: `0 e. _' [
alpha=0.05;
# v: C! I8 W! Q. U8 W$ L0 Z
% 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
0 M* r+ v7 o* j: s) t6 D
A=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]';
! ^. ?' [& e! h( V0 T
[mu,sigma]=normfit(A);
% I/ P; b6 e. H$ s2 Q0 p
p1=normcdf(A,mu,sigma);
2 O# t4 Q4 v3 D$ b
[H1,s1]=kstest(A,[A,p1],alpha)
1 t$ B! H2 y( A9 ?/ n& L# Y5 @" E
n=length(A);
0 _0 ~$ h. ~) G; H+ _4 X
if H1==0
E' l+ v" c8 K% ?' q
disp('该数据源服从正态分布。')
* c) s' ]8 J E
else
$ X m0 w; E+ i3 H) G/ g }2 P
disp('该数据源不服从正态分布。')
/ A7 v' s: ] U) h3 S( \
end
/ b5 Z( n3 c+ Z% [' W7 i8 |
phat=gamfit(A,alpha);
2 K+ ^ ]5 O$ x# X9 v [! ~2 |4 K
p2=gamcdf(A,phat(1),phat(2));
7 T. b X" \/ L: o; Y+ Z
[H2,s2]=kstest(A,[A,p2],alpha)
0 l6 P- O; ]8 G* g. V
if H2==0
5 a9 d& Y8 t& H3 z6 [
disp('该数据源服从γ分布。')
" H3 H1 V, H( _" r
else
4 V% s& x$ Z, G3 K5 O- S9 B+ q7 t
disp('该数据源不服从γ分布。')
! E% q) ]+ D$ v0 d4 @1 v2 S
end
; f" X4 L- g6 a7 p' U) w3 o3 Y
lamda=poissfit(A,alpha);
* T9 v% |! y5 R% F. d4 {$ P
p3=poisscdf(A,lamda);
& J* m9 v3 o+ E# N B
[H3,s3]=kstest(A,[A,p3],alpha)
, @. k- ~ e) F& p2 Z
if H3==0
9 G$ b$ c, O0 \$ f
disp('该数据源服从泊松分布。')
, ?& D* g# p3 n: {) m Y
else
; t9 |3 k* b1 j" N6 {
disp('该数据源不服从泊松分布。')
6 b! P) ?- Q1 H& L; L
end
$ F# {+ m6 R4 v9 c0 Z9 l2 n& E
mu=expfit(A,alpha);
# [" z. G- x- ]% \# W- l
p4=expcdf(A,mu);
' I' Y! L" Z& m. c" O
[H4,s4]=kstest(A,[A,p4],alpha)
b3 \1 W2 F% j& \& [+ a6 p
if H4==0
5 s* v' A3 w7 V" q+ o
disp('该数据源服从指数分布。')
4 C! J6 k3 ?4 w8 R
else
& u9 d9 ?) c! P E
disp('该数据源不服从指数分布。')
4 ~4 C# k4 U3 p
end
( C$ z: z: a2 j1 b" {, }
[phat, pci] = raylfit(A, alpha)
2 \$ }+ @0 {1 q1 G4 g* r
p5=raylcdf(A,phat);
7 |/ `7 Q% s: c' R3 T
[H5,s5]=kstest(A,[A,p5],alpha)
8 K/ W7 G; S% e H, v. r
if H5==0
`0 b$ a/ q& p; e# b: @5 ~/ E
disp('该数据源服从rayleigh分布。')
' f* D; f& _# k; _+ G. c, d
else
6 _/ H& W7 c# D5 @1 m2 x
disp('该数据源不服从rayleigh分布。')
/ I; b8 h2 y: P. y
end
作者:
lillianxiong
时间:
2012-8-20 17:16
上面的是什么？
作者:
数学光头
时间:
2012-8-20 18:17
lillianxiong 发表于 2012-8-20 17:16
6 u+ H5 E# t: n: n3 ?% h
上面的是什么？
( s9 U: b+ H0 m C
你是说整个么，是代码
作者:
一场梦12
时间:
2012-8-20 19:32
有点深奥 ，不错
作者:
小盆子。。。
时间:
2013-8-23 13:13
感觉不是很懂得感觉
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
6 u+ H5 E# t: n: n3 ?% h