
- •Методические материалы
- •1. Основы имитационного моделирования1
- •1.1. Общие положения
- •1.2. Событийный подход к построению алгоритма модели
- •1.3. Подход сканирования активностей
- •1.4. Учет законов распределения случайных параметров
- •1.5. Возможные показатели, оцениваемые с помощью имитационной модели. Проверка гипотез
- •1.6. Схемы проведения имитационных экспериментов. Статистическая обработка результатов экспериментов
- •1.6.1. Метод повторений
- •1.6.2. Метод подынтервалов
- •1.6.3. Метод циклов
- •2. Пример построения блок-схемы имитационной модели
- •3. Способы понижения погрешностей выборочных оценок расчетных параметров
- •3.1. Дополняющие выборки
- •3.2. Общие потоки случайных чисел
- •3.3. Управляющие переменные
- •3.4. Стратифицированные выборки
- •3.5. Значимые выборки
- •4. Некоторые типы систем массового обслуживания (для проведения имитационных экспериментов)
- •4.1. Модели массового обслуживания с пуассоновскими потоками
- •4.1.1. Обобщенная стационарная модель смо
- •4.1.2. Модель смо (m/m/1):(gd//)
- •4.1.3. Модель смо (m/m/1):(gd/n/)
- •4.1.4. Модель смо (m/m/с):(gd//)
- •4.1.5. Модель смо (m/m/с):(gd/n/), cn
- •4.1.6. Модель смо (m/m/):(gd//)
- •4.1.7. Модель смо (m/m/r):(gd/k/k)
- •4.2. Модель смо (m/g/1):(gd//). Формула Поллачека-Хинчина
- •Приложение
- •Рекомендуемая литература
4.1.3. Модель смо (m/m/1):(gd/n/)
Эта модель отличается от предыдущей ограничением на емкость системы: система вмещает не более N клиентов, соответственно, максимальная длина очереди равна N–1.
Таким образом, формально можем положить:
n
=
0, n N.
n = , n = 0, 1, …
рn
=
(4.12)
0, n > N.
Соответственно, уравнение для p0 имеет вид p0(1 + + … + N) = 1. Таким образом:
р0
=
(4.13)
,
=1,
В этой модели значение параметра = / не обязательно должно быть меньше 1, поскольку поступление клиентов в систему ограничивается ее конечной емкостью. Таким образом, в данном случае интенсивность поступления клиентов в систему естественно характеризовать величиной
эфф = (1 – рN).
Следует ожидать, что эфф будет меньше . {Почему?}
С учетом (4.3), (4.12) и (4.13) среднее число клиентов в системе равно:
Ls
=
(4.14)
,
=1,
{Студенты
должны самостоятельно проделывать
вывод последней формулы. Первую формулу
из (4.14) можно получить с помощью суммы
для арифметико-геометрической прогрессии:
12.
Возможны иные подходы.}
Зная значения Ls и эфф можно получить значения других функциональных характеристик.
4.1.4. Модель смо (m/m/с):(gd//)
Эта модель обобщает модель из п. 4.1.2 по параметру с (число каналов обслуживания).
Пусть интенсивность поступления клиентов в систему, интенсивность обслуживания клиентов одним каналом.
Так как нет ограничений на число клиентов в системе, то эфф = .
Наличие параллельных каналов обслуживания приводит к возрастанию фактической интенсивности обслуживания (здесь также лучше говорить об интенсивности выходного потока): если nс, то интенсивность равна n; если n>с, то интенсивность равна c. Используя принятые выше обозначения, можно записать:
n
n
=
c, n>с,
n = , n = 0, 1, … .
Соответственно,
pn
=
(4.15)
,
n>с,
Если предположить, что /c<1, то можно получить следующую формулу для p0:
.
(4.16)
Используя (4.15) и (4.16), можно рассчитать функциональные характеристики. Например, среднее число клиентов в очереди:
Lq
=
.
Остальные характеристики определяются в соответствии с общими формулами {проделать самостоятельно}.
4.1.5. Модель смо (m/m/с):(gd/n/), cn
Эта модель обобщает модель из п. 4.1.3 по параметру с (число каналов обслуживания).
Максимальная длина очереди в системе равна N–c.
Пусть интенсивность поступления клиентов в систему, интенсивность обслуживания клиентов одним каналом.
Приведем некоторые расчетные соотношения для параметров и характеристик системы.
n
=
0, n>N .
n
n
=
c, cnN .
pn
=
,
cnN
.
р0
=
.
Lq
=
.
эфф = (1 – рN).
4.1.6. Модель смо (m/m/):(gd//)
В этой модели неограниченны емкости источника заявок и системы и, кроме того, полагается неограниченным количество каналов обслуживания. На практике такие модели соответствуют ситуации, когда клиент, одновременно выступает в роли канала (само)обслуживания.
Полагается, что интенсивность поступления клиентов в систему постоянна, а также постоянна интенсивность обслуживания клиентов одним каналом.
В обозначениях общей модели можно записать:
n = ; n = n, n = 0, 1, …,
Следовательно:
.
Из равенства
следует, что
.
Откуда получаем:
,
n = 0, 1, ….
Таким образом вероятности состояний системы в этой модели совпадают с вероятностями распределения Пуассона с математическим ожиданием Ls = . Кроме того в этой системе, очевидно Lq = Wq = 0.