Автор: Пользователь скрыл имя, 19 Декабря 2010 в 15:54, лабораторная работа
Объектом исследования данной работы является система противоракетной обороны (ПРО). Будем рассматривать систему ПРО как систему массового обслуживания, т.е. как -канальную систему, на вход которой поступают ракеты противника. Ракеты, пролетающие в пределах полосы налета, могут быть обстреляны любым из каналов данной системы ПРО.
Ширина полосы налета определяется возможностями обстрела всеми n каналами любой цели в пределах полосы налета. Предполагается, что если ракета летит выше пределов полосы налета (слева или справа), то эти ракеты не могут быть обстреляны ни одним из n каналов данной системы ПРО.
Техническое задание 3
Введение 4
Разработка имитационной модели 8
Математическое описание имитационной модели. 8
Описание блок-схемы алгоритма. 9
Анализ полученных результатов 13
Разработка аналитической модели. 14
Математическое описание аналитической модели 14
Расчёт параметров СМО. 16
Анализ полученных результатов моделирования. 18
Выводы 19
Список используемой литературы 20
Приложения 21
Приложение 1 (листинг программы для имитационного моделирования) 21
Приложение 2 (листинг программы для аналитического моделирования) 23
Приложение 3. Таблицы результатов моделирования. 24
Приложение 4. Графики, полученные в результате моделирования. 25
function faktorial(a:integer):integer;
var ii,fak:integer;
begin
fak:=1;
if a=0 then fak:=1;
for ii:=1 to a do fak:=fak*ii;
faktorial:=fak;
end;
var
p: array [0..n] of real;
i,h,ns:integer;
c:real;
begin
h:=n div ll;
ns:=h;
c:=0;
if n mod ll <>0 then ns:=ns+1;
for i:=0 to ns do p[i]:=0;
for i:=0 to ns-1 do c:=c+exp(i*ln(l))/(faktorial(
c:=c+exp(ns*ln(l))/(faktorial(
p[0]:=1/c;
for i:=1 to ns-1 do p[i]:=(p[i-1]*l)/(i*ll*m);
{exp(i*ln(l))/(faktorial(i)*
p[ns]:=p[ns-1]*l/(n*m);
{exp(ns*ln(l))/(faktorial(ns)*
for i:=0 to ns do writeln('p[',i,']=',p[i]:4:5);
writeln('n=',ns);
writeln('Pobsl=',(1-p[ns]):3:
writeln('l0=',(1-p[ns]):l:3:6)
writeln('ksr=',(l/m)*p[ns]:3:
readln;
end.
Результаты имитационного моделирования:
Число
каналов обслуживания, |
Номер реализации | ||||
1 | 0,443 | 0,560 | 1,680 | ||
2 | 0,449 | 0,554 | 1,663 | ||
3 | 0,448 | 0,552 | 1,656 | ||
4 | 0,446 | 0,555 | 1,665 | ||
5 | 0,446 | 0,557 | 1,672 | ||
6 | 0,446 | 0,555 | 1,666 | ||
7 | 0,443 | 0,560 | 1,680 | ||
8 | 0,446 | 0,580 | 1,673 | ||
9 | 0,448 | 0,556 | 1,667 | ||
10 | 0,445 | 0,555 | 1,664 | ||
Средние значения | 0,446 | 0,5584 | 1,6686 | ||
1 | 0,652 | 0,554 | 2,216 | ||
2 | 0,649 | 0,557 | 2,229 | ||
3 | 0,647 | 0,560 | 2,241 | ||
4 | 0,649 | 0,558 | 2,233 | ||
5 | 0,650 | 0,556 | 2,225 | ||
6 | 0,656 | 0,553 | 2,214 | ||
7 | 0,650 | 0,557 | 2,226 | ||
8 | 0,650 | 0,554 | 2,218 | ||
9 | 0,648 | 0,557 | 2,226 | ||
10 | 0,653 | 0,556 | 2,226 | ||
Средние значения | 0,6504 | 0,5562 | 2,2254 | ||
1 | 0,706 | 0,512 | 2,560 | ||
2 | 0,710 | 0,514 | 2,568 | ||
3 | 0,707 | 0,509 | 2,547 | ||
4 | 0,707 | 0,511 | 2,557 | ||
5 | 0,714 | 0,509 | 2,546 | ||
6 | 0,705 | 0,511 | 2,556 | ||
7 | 0,708 | 0,514 | 2,569 | ||
8 | 0,709 | 0,507 | 2,537 | ||
9 | 0,708 | 0,510 | 2,550 | ||
10 | 0,710 | 0,510 | 2,551 | ||
Средние значения | 0,7084 | 0,5107 | 2,5541 | ||
1 | 0,747 | 0,461 | 2,764 | ||
2 | 0,743 | 0,463 | 2,775 | ||
3 | 0,745 | 0,461 | 2,769 | ||
4 | 0,748 | 0,462 | 2,775 | ||
0,746 | 0,461 | 2,764 | |||
n=6 | 6 |
0,744 | 0,465 | 2,79 | |
7 | 0,745 | 0,464 | 2,781 | ||
8 | 0,743 | 0,461 | 2,766 | ||
9 | 0,743 | 0,467 | 2,802 | ||
10 | 0,742 | 0,466 | 2,798 | ||
средние значения | 0,7446 | 0,4631 | 2,7784 |
Сравнение результатов имитационного и аналитического моделирования:
Число каналов обслуживания, n | Искомый параметр | Аналитическая модель | Имитационная модель |
Робс | 0,4454 | 0,446 | |
n=3 | πз.к. | 0,556 | 0,5584 |
к̅ | 1,66 | 1,6686 | |
Робс | 0,6487 | 0,6504 | |
n=4 | πз.к. | 0,5561 | 0,5562 |
к̅ | 2,221 | 2,2254 | |
Робс | 0,702 | 0,7084 | |
n=5 | πз.к. | 0,5102 | 0,5107 |
к̅ | 2,55 | 2,5541 | |
Робс | 0,7435 | 0,7446 | |
n=6 | πз.к. | 0,4625 | 0,4631 |
к̅ | 2,772 | 2,7784 |
График 1. Зависимость вероятности обслуживания от числа приборов для аналитической и имитационной моделей.
График 2. Зависимость вероятности занятости канала от числа приборов для аналитической и имитационной моделей
График 3. Зависимость среднего числа занятых каналов от числа приборов для аналитической и имитационной моделей.