- •Основы теории массового обслуживания Методические указания
- •Основы теории массового обслуживания
- •Оглавление
- •Лабораторная работа №1 Введение в Mathcad. Переменные, функции, графика
- •1.1 Интерфейс пользователя
- •1.2 Области рабочего документа
- •1.3 Определение переменных
- •1.6 Работа с функциями
- •1.7 Выделение выражения
- •1.8 Построение двумерных графиков в декартовой системе координат
- •1.4 Определение дискретного аргумента
- •1.5 Ввод текста
- •1.9 Построение графиков в трехмерной системе координат
- •1.10 Построение нескольких графиков в одном графическом регионе
- •1.11 Форматирование графиков
- •1.12 Решение уравнений
- •Лабораторная работа №2 Введение в Mathad. Матричные операции, программирование функций
- •2.1 Состав панели программирования
- •2.2 Программирование в системе Mathcad
- •2.3 Работа с векторами и матрицами
- •2.3.1 Создание вектора или матрицы:
- •2.3.2 Перемножение двух матриц:
- •2.3.3 Среднее и дисперсия:
- •2.4 Генерирование случайных чисел
- •Лабораторная работа №3 Марковские цепи. Определение и построение
- •3.1 Определение Последовательность случайных величин образует дискретную цепь Маркова, если для всех n и всех возможных случайных величин выполняется равенство:
- •3.2 Стохастическая матрица
- •3.3 Неприводимая и однородная цепь Маркова
- •3.4 Эргодическая цепь Маркова
- •3.5 Стохастическая маршрутизация в сетях с коммутацией пакетов
- •Лабораторная работа №4 Марковские цепи. Исследование эргодических свойств
- •4.1 Обозначения и расчетные формулы
- •4.2 Функция для расчета траектории движения пакета по сети
- •Лабораторная работа №5 Система массового обслуживания g/g/1. Формирование управляющих случайных последовательностей
- •5.1 Модель системы массового обслуживания
- •Решение системы уравнений
- •Система m/m/1
- •5.7.2 Гамма – распределение
- •5.7.3 Логнормальное распределение
- •5.7.4 Распределение хи - квадрат
- •Распределение Эрланга
- •Распределение Вейбулла
- •Статистические характеристики
- •Лабораторная работа №6 Система массового обслуживания g/g/1. Исследование зависимостей параметров от типа функций распределения управляющих последовательностей
- •Полное описание модели и полученных в результате моделирования характеристик смотри в прилагающейся к лабораторной работе Mathcad – программе «Система массового обслуживания».
- •Лабораторная работа № 7 Система массового обслуживания m/g/1. Формула Хинчина –Поллячека
- •7.1 Характеристики m/g/1
- •7.2 Характеристики m/d/1
- •7.3 Характеристики m/м/1
- •Литература
- •Основы теории массового обслуживания
Лабораторная работа №5 Система массового обслуживания g/g/1. Формирование управляющих случайных последовательностей
Цель работы: Моделирование управляющих случайных последовательностей с различными функциями распределения.
Подготовка к лабораторной работе:
-
Изучить встроенные функции системы Mathcad.
-
Изучить различные законы распределения случайных величин.
-
Повторить обозначения систем массового обслуживания.
-
Повторить понятия входного потока и времени обслуживания.
-
Повторить понятия математическое ожидание и дисперсия.
-
Изучить формулы для расчета математического ожидания и дисперсии для заданных распределений.
Краткая теория:
5.1 Модель системы массового обслуживания
В предлагаемой лабораторной работе исследуется модель системы массового обслуживания (СМО), состоящей из очереди и прибора.
Р
исунок
5.1 Модель СМО.
На вход системы поступает поток заявок с интенсивностью . Обслуживание осуществляется в единственном приборе с интенсивностью . Подразумевается, что очередь имеет неограниченную емкость (т.е. обслуживание без отказов). Любое требование, поступающее в систему, в то время когда прибор занят, ждет пока не обслужатся все стоящие перед ним в очереди требования. Входными параметрами для исследования СМО являются последовательности {n} – последовательность промежутков времени между поступлениями требований и {n} – последовательность промежутков времени обслуживания.
Задача данной лабораторной работы заключается в том, чтобы правильно сформировать управляющие последовательности в зависимости от вида СМО.
Для сравнения систем массового обслуживания различного типа необходимо, чтобы аналитические характеристики (математическое ожидание и дисперсия) распределений случайных величин ({n} и {n}), подаваемых на вход этих СМО, совпадали. Таким образом, параметры распределений определяются решением системы уравнений:
(
5.1)
-
Решение системы уравнений
Для решения системы уравнений в Mathcad можно воспользоваться встроенной функцией Find. Это можно сделать следующим образом:
-
Присвоить начальные значения всем искомым переменным;
-
Ниже ввести ключевое слово Given (Дано);
-
Перейти на строку ниже и ввести систему уравнений (использовать жирный знак равенства из панели отношений);
-
Ниже выполнить вызов функции Find(x,y).
Пример:

-
Система m/m/1
Распределение времени между поступлениями – показательное:
, (5.2)
со средним:
(5.3)
и дисперсией:
. (5.4)
Распределение
времени обслуживания – показательное
.
Среднее и дисперсия
вычисляются по формулам 5.3 - 5.4 (с
параметром ).
Для генерации показательно распределенных
случайных чисел используется встроенная
функция Mathcad
– rexp(N,),
где N
– размерность
вектора
– параметр
распределения.
-
Система M/G/1
Распределение времени между поступлениями – показательное (среднее и дисперсия смотри выше). Распределение времени обслуживания – общего вида, в данной лабораторной работе это любое распределение из таблицы 5.1. Характеристики распределений общего вида описаны ниже.
-
Система G/M/1
Распределение времени между поступлениями – общего вида. Распределение времени обслуживания – показательное.
-
Система G/G/1
Распределение времени между поступлениями заявок и времени обслуживания – общего вида.
-
Распределения общего вида и их характеристики
5.7.1 Равномерное распределение
Функция распределения
случайной величины x,
равномерно распределенной
на интервале
выглядит следующим образом:
. (5.5)
Плотность вероятности:
, (5.6)
со средним:
(5.7)
и дисперсией:
. (5.8)
Генерация вектора
из N
случайных чисел, равномерно распределенных
на интервале
,
производится с помощью встроенной
функции Mathcad
– runif(N,a,b).
Параметры
распределения (границы интервала)
находятся решением системы уравнений
(5.1).
