Компьютерное моделирование. Тарасов В.Н - 117 стр.

UptoLike

7
eks:=eks/power(sig,4)-3;
Label16.Caption:='ek: '+ floattostrf(eks,ffGeneral,3,3);
max:=x[1];
min:=x[1];
for i:=2 to n do begin
if x[i]>max then max:=x[i] else
if x[i]<min then min:=x[i];
end;
ddx:=(max-min)/20;
dx:=min;
j:=1;
for i:=1 to n do begin
if x[i]=max then inc(col[20]) else
j:=(trunc(x[i]*100000)div trunc(ddx*100000))+1;
inc(col[j]);
end;
{while dx<max do begin
for i:=1 to n do
if (x[i]>=dx)and(x[i]<=dx+ddx) then col[j]:=col[j]+1;
j:=j+1;
dx:=dx+ddx;
end;
dx:=min;
s:='';
memo1.lines.Clear;
for i:=1 to 20 do begin
s:=floattostrf(dx,ffGeneral,3,3)+'--'+floattostrf(dx+ddx,ffGeneral,3,3)+' '+inttostr(col[i]);
memo1.Lines.Add(s);
dx:=dx+ddx;
end; }
for i:=1 to 20 do
pjk[i]:=col[i]/n;
with series2 do clear;
for i:=1 to 20 do chart2.Series[0].AddY(pjk[i]);
sum:=0;
sum1:=0;
for i:=1 to 20 do sum:=sum+pjk[i];
for i:=1 to 20 do sum1:=sum1+col[i];
{label6.Caption:=(floattostr(sum));
label7.Caption:=(floattostr(sum1));}
dx:=0;
i:=1;
while dx<=1 do begin
pj[i]:=(exp((-1)*lamz*dx)-exp((-1)*lamz*(dx+0.05)));
dx:=dx+0.05;
i:=i+1;
end;
xi2:=0;
for i:=1 to 20 do xi2:=xi2+sqr(n*pjk[i]-n*pj[i])/(n*pj[i]);
{ label6.Caption:='Согласно критерию Пирсона:';
 eks:=eks/power(sig,4)-3;
 Label16.Caption:='ek: '+ floattostrf(eks,ffGeneral,3,3);
 max:=x[1];
 min:=x[1];
for i:=2 to n do begin
 if x[i]>max then max:=x[i] else
 if x[i]=dx)and(x[i]<=dx+ddx) then col[j]:=col[j]+1;
 j:=j+1;
 dx:=dx+ddx;
end;
 dx:=min;
 s:='';
 memo1.lines.Clear;
for i:=1 to 20 do begin
 s:=floattostrf(dx,ffGeneral,3,3)+'--'+floattostrf(dx+ddx,ffGeneral,3,3)+' '+inttostr(col[i]);
 memo1.Lines.Add(s);
 dx:=dx+ddx;
             end; }
for i:=1 to 20 do
 pjk[i]:=col[i]/n;
 with series2 do clear;
for i:=1 to 20 do chart2.Series[0].AddY(pjk[i]);
 sum:=0;
 sum1:=0;
for i:=1 to 20 do sum:=sum+pjk[i];
for i:=1 to 20 do sum1:=sum1+col[i];
 {label6.Caption:=(floattostr(sum));
 label7.Caption:=(floattostr(sum1));}
 dx:=0;
 i:=1;
while dx<=1 do begin
 pj[i]:=(exp((-1)*lamz*dx)-exp((-1)*lamz*(dx+0.05)));
 dx:=dx+0.05;
 i:=i+1;
end;
  xi2:=0;
for i:=1 to 20 do xi2:=xi2+sqr(n*pjk[i]-n*pj[i])/(n*pj[i]);
 { label6.Caption:='Согласно критерию Пирсона:';


                                                                                                 7