МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ
ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
ЛИПЕЦКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАФЕДРА АВТОМАТИЗИРОВАННЫХ СИСТЕМ УПРАВЛЕНИЯ
Лабораторная работа №4
«Определение стационарного режима при моделировании на GPSS»
по дисциплине
«Моделирование систем»
Вариант 7, 11
|
|
Студент |
|
|
|
Ретивых С.Ю. |
| ||||||||||
|
|
|
|
подпись, дата |
|
фамилия, инициалы |
| ||||||||||
|
|
Группа |
|
АС-09 |
|
|
| ||||||||||
|
|
|
|
|
|
|
| ||||||||||
|
|
Принял |
|
|
|
|
| ||||||||||
|
|
|
|
|
|
Гаев Л.В. |
| ||||||||||
|
|
ученая степень, звание |
|
подпись, дата |
|
фамилия, инициалы |
| ||||||||||
Липецк 2012
Задание кафедры
Определить интервалы времени моделирования, после завершения которых можно говорить о установившемся стационарном режиме. Исследования проводить при различных загруженностях прибора.
Закон распределения дискретной случайной величины (11):
|
Возможные значения |
a |
a + 1 |
7a |
8a |
|
Вероятности |
0.39 |
0.02 |
0.53 |
0.06 |
Закон распределения непрерывной случайной величины (7):
|
xi |
b |
(3b+c)/4 |
(b+c)/2 |
(b+3c)/4 |
c |
|
yi |
Y |
0 |
Y |
0 |
Y |
2. Подбор параметров
Рассмотрим непрерывный закон распределения.
Поскольку область симметрична, матожидание будет посередине и равно (b+c)/2.
Площадь, ограничиваемая данной фигурой, должна равняться 1, следовательно
y*(c-b)/2=1, y=2/(c-b).
Примем b=100, c=200. Тогда y=0,02. M(x)=150.
Функция плотности распределения:
Функция распределения случайной величины:
Рис. 1 – график функции распределения
Рассмотрим дискретный закон. Его матожидание:
M(X)=0,4a+0,9a+0,24a+0,28a=1,82a
Теперь подберем параметры непрерывного распределения таким образом, чтобы
|
a |
M(X)н |
M(X)д |
Загр. Прибора |
|
45 |
150 |
224,5 |
0,7 |
|
40 |
150 |
199,6 |
0,8 |
|
36 |
150 |
179,68 |
0,9 |
|
34 |
150 |
169,92 |
0,95 |
|
32 |
150 |
159,96 |
0,99 |
Код программы на языке gpss/pc Система с загруженностью 0,7
DRP FUNCTION RN1,D4
.39,45./.41,46/.94,315/1,360
NRP FUNCTION RN2,C21
0,100/.09,105/.16,110/.21,115/.24,120/.25,125/.26,130/.29,135/
.34,140/.41,145/.5,150/.59,155/.66,160/.71,165/.74,170/.75,175/
.76,180/.79,185/.84,190/.91,195/1,200
1 GENERATE FN$DRP
2 QUEUE OCH
3 SEIZE PRB
4 DEPART OCH
5 ADVANCE FN$NRP
6 RELEASE PRB
7 TERMINATE
8 PLOT QT$OCH,400,0,100000
9 GENERATE 100000,0
10 TERMINATE 1
11 START 1
При моделировании появится окно, в котором будет графически отображаться значение среднего времени пребывания в очереди.
Из графика видно, что стационарный режим наступает в момент времени около 50 000. Запустим программу еще раз, и в окрестности этого времени прогоним пошагово. Последний раз среднее время пребывания в очереди случается в момент времени 45750, далее система работает в стационарном режиме.
GPSS/PC Report file REPORT.GPS. (V 2, # 37471) 04-18-2012 19:00:17 page 1
START_TIME END_TIME BLOCKS FACILITIES STORAGES FREE_MEMORY
0 45750 8 1 0 16864
LINE LOC BLOCK_TYPE ENTRY_COUNT CURRENT_COUNT RETRY
1 1 GENERATE 230 0 0
2 2 QUEUE 230 1 0
3 3 SEIZE 229 0 0
4 4 DEPART 229 0 0
5 5 ADVANCE 229 1 0
6 6 RELEASE 228 0 0
7 7 TERMINATE 228 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 229 0.728 145.52 1 230 0 0 0 1
QUEUE MAX CONT. ENTRIES ENTRIES(0) AVE.CONT. AVE.TIME AVE.(-0) RETRY
OCH 8 1 230 91 0.68 135.04 223.45 0