我想使用chi2gof命令來分析一組大樣本數(shù)據(jù)的分布情況，程序如下：function f=p_chi2gof(x,alpha)

% 本程序用于判別所給數(shù)據(jù)源的概率分布情況。

h1 = chi2gof(x,'cdf',@(z)normcdf(z,mean(x),std(x)),'nparams',2);

if h1==0

h1

disp('該數(shù)據(jù)源服從正態(tài)分布。')

else

disp('該數(shù)據(jù)源不服從正態(tài)分布。')

end

lamda=poissfit(x,alpha);

h2 = chi2gof(x,'cdf',@(z)poisscdf(z,lamda),'nparams',1);

h2;

if h2==0

disp('該數(shù)據(jù)源服從泊松分布。')

else

disp('該數(shù)據(jù)源不服從泊松分布。')

end

[ahat,bhat]=unifit(x,alpha);

h3= chi2gof(x,'cdf',@(z)unifcdf(z,ahat,bhat),'nparams',2);

if h3==0

disp('該數(shù)據(jù)源服從均勻分布。')

else

disp('該數(shù)據(jù)源不服從均勻分布。')

end

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

h4 = chi2gof(x,'cdf',@(z)raylcdf(z,phat),'nparams',1);

if h4==0

disp('該數(shù)據(jù)源服從rayleigh分布。')

else

disp('該數(shù)據(jù)源不服從rayleigh分布。')

end

mu=expfit(x,alpha);

h5= chi2gof(x,'cdf',@(z)gamcdf(z,mu),'nparams',1);

if h5==0

disp('該數(shù)據(jù)源服從指數(shù)分布。')

else

disp('該數(shù)據(jù)源不服從指數(shù)分布。')

end

phat=gamfit(x,alpha);

h6= chi2gof(x,'cdf',@(z)gamcdf(z,phat(1),phat(2)),'nparams',2);

if h6==0

disp('該數(shù)據(jù)源服從γ分布。')

else

disp('該數(shù)據(jù)源不服從γ分布。')

end

之后，我用各種隨機(jī)數(shù)生成命令分別生成了相應(yīng)分布類型的隨機(jī)數(shù)數(shù)組

norm=normrnd(2,5,[1 1000]);

beta=betarnd(2,6,[1 1000]);

logn=lognrnd(3,6,[1 1000]);

exp=exprnd(10,[1 1000]);

poiss=poissrnd(8,[1 1000]);

rayl=raylrnd(30,[1 1000]);

gam=gamrnd(6,9,[1 1000]);

unif=unifrnd(1,100,[1 1000]);

再利用函數(shù)p_chi2gof()來檢驗(yàn)這些數(shù)組相應(yīng)的分布形式

出現(xiàn)了幾個(gè)問題：

檢驗(yàn)指數(shù)分布隨機(jī)數(shù)據(jù)的結(jié)果為：該數(shù)據(jù)源不服從指數(shù)分布，該數(shù)據(jù)源服從γ分布。

檢驗(yàn)泊松分布隨機(jī)數(shù)據(jù)的結(jié)果為：該數(shù)據(jù)源不服從泊松分布(且不服從任意其他分布)

檢驗(yàn)瑞利分布隨機(jī)數(shù)據(jù)的結(jié)果為：該數(shù)據(jù)源服從rayleigh分布，該數(shù)據(jù)源服從γ分布。

會(huì)出現(xiàn)這樣的問題回事什么原因？我在程序上有哪些不正確的地方嗎？

謝謝^_^

總結(jié)

以上是生活随笔為你收集整理的matlab三角定位分析程序,Matlab chi2gof 分布分析程序求教的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。