
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ
ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
ЛИПЕЦКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАФЕДРА АВТОМАТИЗИРОВАННЫХ СИСТЕМ УПРАВЛЕНИЯ
Лабораторная работа №4
«Определение стационарного режима при моделировании на GPSS»
по дисциплине
«Моделирование систем»
Вариант 9, 4
|
Студент |
|
|
|
Дырдин Н.В. |
| ||||||||||
|
|
|
подпись, дата |
|
фамилия, инициалы |
| ||||||||||
|
Группа |
|
АС-09 |
|
|
| ||||||||||
|
|
|
|
|
|
| ||||||||||
|
Принял |
|
|
|
|
| ||||||||||
|
|
|
|
|
Гаев Л.В. |
| ||||||||||
|
ученая степень, звание |
|
подпись, дата |
|
фамилия, инициалы |
|
Липецк 2012
Задание кафедры
Определить интервалы времени моделирования, после завершения которых можно говорить о установившемся стационарном режиме. Исследования проводить при различных загруженностях прибора.
Закон распределения дискретной случайной величины (4):
Возможные значения |
a |
2a |
3a |
4a |
Вероятности |
0.4 |
0.45 |
0.08 |
0.07 |
Закон распределения непрерывной случайной величины (9):
xi |
d |
(4d+e)/5 |
(3d+2e)/5 |
(2d+3e)/4 |
(d+4e)/5 |
e |
yi |
0 |
Y |
0 |
0 |
Y |
0 |
2. Подбор параметров
Рассмотрим непрерывный закон распределения.
Мат. ожидание:
M(P) = y(e^2 – d^2)/5
Площадь должна быть равна 1:
2y*(e-d)/5=1, y=5/(2e-2d).
Примем e=300, d=50. Тогда y=0,01. M(x)=175.
Функция плотности распределения:
Функция распределения случайной величины:
Рассмотрим дискретный закон. Его мат. ожидание:
M(X)=0,4a+0,9a+0,24a+0,28a=1,82a
Теперь подберем параметры непрерывного распределения таким образом, чтобы
a |
M(X)н |
M(X)д |
Загр. Прибора |
137 |
175 |
249,34 |
0,7 |
120 |
175 |
218,4 |
0,8 |
107 |
175 |
194,74 |
0,9 |
101 |
175 |
183,82 |
0,95 |
97 |
175 |
176,54 |
0,99 |
Код программы на языке gpss/pc Система с загруженностью 0,7
GEN1 FUNCTION RN1,D4
.4,137/.85,274/.93,411/1,548
ADV1 FUNCTION RN2,C26
0,50/.01,60/.04,70/.09,80/.16,90/.25,100/.34,110/.41,120/.46,130/
.49,140/.5,150/.5,160/.5,170/.5,180/.5,190/.5,200/.51,210/.54,220/
.59,230/.66,240/.75,250/.84,260/.91,270/.96,280/.99,290/1,300
1 GENERATE FN$GEN1
2 QUEUE OCH
3 SEIZE PRB
4 DEPART OCH
5 ADVANCE FN$ADV1
6 RELEASE PRB
7 TERMINATE
8 PLOT QT$OCH,400,0,80000
9 GENERATE 100000,0
10 TERMINATE 1
START 1
При моделировании появится окно, в котором будет графически отображаться значение среднего времени пребывания в очереди.
Из графика видно, что стационарный режим наступает в момент времени после 60 000. Запустим программу еще раз, и в окрестности этого времени прогоним пошагово. Последний раз среднее время пребывания в очереди случается в момент времени 80000, далее система работает в стационарном режиме.
GPSS/PC Report file REPORT.GPS. (V 2, # 37471) 04-19-2012 09:38:29 page 1
START_TIME END_TIME BLOCKS FACILITIES STORAGES FREE_MEMORY
0 80000 9 1 0 17328
LINE LOC BLOCK_TYPE ENTRY_COUNT CURRENT_COUNT RETRY
1 1 GENERATE 324 0 0
2 2 QUEUE 324 1 0
3 3 SEIZE 323 1 0
4 4 DEPART 322 0 0
5 5 ADVANCE 322 0 0
6 6 RELEASE 322 0 0
7 7 TERMINATE 322 0 0
9 8 GENERATE 1 0 0
10 9 TERMINATE 1 0 0
FACILITY ENTRIES UTIL. AVE._TIME AVAILABLE OWNER PEND INTER RETRY DELAY
PRB 323 0.697 172.81 1 324 0 0 0 1
QUEUE MAX CONT. ENTRIES ENTRIES(0) AVE.CONT. AVE.TIME AVE.(-0) RETRY
OCH 3 2 324 180 0.27 67.33 151.49 0