Отчет по лабораторной работе №8

Проверка гипотезы о законе распределения

Выполнили: Шубенок Е.М., Корецкая В.С.

Текст исходной программы:

path(path,'e:\shik\laba8\');

clear;

clc;

alpha=0.05;

n=100;

par(1)=2;

par(2)=4;

for i=1:100

x(i)=normrnd(par(1),par(2));

end;

iclass=0;

cint=[0 0 0 0 0 0 0 0 0 0];

k=10;

[cint,ifreq,xmin,xmax,ifail]=g01aef(x,iclass,cint);

v=mynormHI(x,par)

npest=0;

p=normcdf(v,0,1); %[chisq,p,ndf,eval,chisqi,ifail]=g08cgf(ifreq,cint(1:k-1),'n',par,npest);

if p>(1-alpha)

disp('OK')

else

disp('nonOK')

end

v;

vA=chi2inv(0.95,k-1);

if v>vA

disp('myok');

else

disp('myNONok');

end

Текст м-файла mynormHI.m

function v=mynormHI (x,par)

n=100;

iclass=0;

cint=[0 0 0 0 0 0 0 0 0 0];

k=10;

[cint,ifreq,xmin,xmax,ifail]=g01aef(x,iclass,cint);

v=0;

for i=1:k

if i==1

q(i)=normcdf(cint(i),par(1),par(2));

else

if i==k

q(i)=1-normcdf(cint(i),par(1),par(2));

else

q(i)=normcdf(cint(i),par(1),par(2))-normcdf(cint(i-1),par(1),par(2));

end

end

v=v+((ifreq(i)-n*q(i))^2)/(n*q(i));

end

Результат выполнения программы (5 запусков):

v = 78.2107

OK

myok

v = 77.9459

OK

myok

v = 72.8832

OK

myok

v = 72.9150

OK

myok

v = 75.4989

OK

myok