|
你的数据服从正态分布。
以前写过一个假设检验的函数,主要是利用K-S检验的一个MATLAB函数kstest。其他分布自己按照格式照猫画虎填写即可:
- function f=p_judge(A,alpha)
- % 本程序用于判别所给数据源在置信率为0.05时的概率分布形式。A的形式为n×1。
- 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
- [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
复制代码 |
评分
-
1
查看全部评分
-
|