Теория массового обслуживания и математическая статистика

                            sum_otkaz[f]:=sum_otkaz[f]+1;

                            goto 1; {увеличиваем сетчик отказов}

                        end;

                    end;

1:              end

            else

            begin

                sum_obr[f]:=sum_obr[f]+1; {увеличиваем счетчик  обращений}

                x[i,j]:=1; {иначе элемент равен единице}

            end;

        end

        else if (x[i,j]>l[i]*(0.005))  then  {прверяем  следующее условие}

        begin x[i,j]:=0; {считаем, что компьютер не пытается подключиться к сети, и ставим в ячейку 0.}

            if x[i,j-1]=1 then

            begin

                {если в предыдущем столбце этой строки

                стоит 1 (компьютер в сети), проверяем,  отключится

                ли он в текущий момент времени  следующим образом}

                if (x[i,j]<=mu[i]*(0.005)) then x[i,j]:=0

                    {читаем, что компьютер выходит из сети и освобождает

                    канал, и в ячейку ставим 0. В  противном случае

                    ставим 1(компьютер продолжает в  сети работать).}

                else x[i,j]:=1;

            end;

        end;

        write(file_,x[i,j]:7);

    end;

    writeln(file_);

end;

 

 

  {-----------------оценка  вероятности----------------}

 

    textcolor(5);

    writeln(file_,'obr = ',sum_obr[f]);

    writeln(file_,'otk = ',sum_otkaz[f]);

 

    p[f]:=sum_otkaz[f]/sum_obr[f];

    writeln(file_,'p[',f,']=  ',p[f]:4:2);

    writeln(file_,'==============================================');

    writeln(file_);

 

end;

 

  otkaz:=0;

  obraw:=0;

  for f:=1 to 10 do begin

    otkaz:=otkaz+sum_otkaz[f];

    obraw:=obraw+sum_obr[f];

  end;

 

  textcolor(4);

  p_sum:=otkaz/obraw; {оцениваем  вероятность отказа}

  write(file_,' summa obraweniy = ',obraw,'; ');

  writeln(file_,'  summa otkazov   = ',otkaz);

  writeln(file_,' p = ',p_sum:5:2);

 

  {-------------- doveritelynuy interval -----------------------}

 

  {Находим границы  доверительного интервала}

  Ga:=(p_sum-U*sqrt((p_sum*(abs(1-p_sum)))/(otkaz)));

  Gb:=(p_sum+U*sqrt((p_sum*(abs(1-p_sum)))/(otkaz)));

  writeln(file_,' doveritelynuy interval = [',Ga:3:2,';',Gb:3:2,']');

 

  textcolor(15);

 

  write(#13#10,' summa obraweniy = ',obraw,'; ');

  writeln(#13#10,'  summa otkazov   = ',otkaz);

  writeln(#13#10,' p = ',p_sum:5:2);

  writeln(#13#10,' doveritelynuy interval = [',Ga:3:2,';',Gb:3:2,']');

  writeln(#13#10,'          Results are saved at a.txt -->>> "Enter"');

 

  readln;

  close(file_);

end.

 

Вывод:

В результате работы программы подсчитываются пары чисел (sum_obraw[f],sum_otkaz[f]), f = 1..10, где sum_obraw[f] – счётчик попыток подключения (обращений) и sum_otkaz[f]- счётчик отказов. Далее программа оцениваем вероятность отказа  p_sum:=otkaz/obraw. Требовалось построить доверительный интервал. Он задается формулой

,

Квантиль нормального распределения  = =U_0,975=1,96.

Вероятность отказа рассчитывалась двумя  методами:

- при решении аналитическим методом вероятность отказа равна 0.30.

- при решении методом Монте-Карло вероятность отказа равна 0.33.

summa obraweniy = 602;   summa otkazov   = 201

p =  0.33

doveritelynuy interval = [0.27;0.40]

Значение, полученное аналитическим  способом попадает в доверительный  интервал.