2024年4月15日发(作者：)

function f=p_judge(A,alpha)

A=A(:);

%正态分布检验

[mu,sigma]=normfit(A);

p1=normcdf(A,mu,sigma);

[H1,s1]=kstest(A,[A,p1],alpha)

n=length(A);

if H1==0

disp('该数据源服从正态分布。')

else

disp('该数据源不服从正态分布。')

end

%γ分布检验

phat=gamfit(A,alpha);

p2=gamcdf(A,phat(1),phat(2));

[H2,s2]=kstest(A,[A,p2],alpha)

if H2==0

disp('该数据源服从γ分布。')

else

disp('该数据源不服从γ分布。')

end

%泊松分布检验

lamda=poissfit(A,alpha);

p3=poisscdf(A,lamda);

[H3,s3]=kstest(A,[A,p3],alpha)

if H3==0

disp('该数据源服从泊松分布。')

else

disp('该数据源不服从泊松分布。')

end

%指数分布检验

mu=expfit(A,alpha);

p4=expcdf(A,mu);

[H4,s4]=kstest(A,[A,p4],alpha)

if H4==0

disp('该数据源服从指数分布。')

else

disp('该数据源不服从指数分布。')

end

%rayleigh分布检验

[phat, pci] = raylfit(A, alpha)

p5=raylcdf(A,phat);

[H5,s5]=kstest(A,[A,p5],alpha)

if H5==0

disp('该数据源服从rayleigh分布。')

else

disp('该数据源不服从rayleigh分布。')

end

%对数正态分布检验

para=lognfit(A);

p6=logncdf(A,para(1),para(2));

[H6,s6]=kstest(A,[A,p6],alpha)

n=length(A);

if H6==0

disp('该数据源服从对数正态分布。')

else

disp('该数据源不服从对数正态分布。')

end