Введение 3
Постановка задачи 5
Обзор литературы 6
Глава 1. Аналитическая модель НСО 8
1.1. Модель НСО 8
1.2. Модель НСО с абсолютными приоритетами 12
Глава 2. Описание алгоритма 16
2.1. Формирование списка возможных состояний 16
2.2. Заполнение матрицы коэффициентов 20
2.3. Реализация решения 21
2.4. Реализация имитационного моделирования 23
Глава 3. Практический пример 25
Заключение 33
Список литературы
Теория массового обслуживания (ТМО) является относительно новым разделом теории вероятности. Целью изучения ТМО является выбор структуры системы обслуживания и процесса обслуживания на основе изучения потоков требований на обслуживание, поступающих в систему и выходящие из неё, длительности ожидания и длины очередей.
Большинство авторов используют модели ТМО в предположении, что очередь заявок бесконечна, существует стационарный режим. Это позволяет оценить эффективность системы обслуживания на длительном интервале времени, получить характеристики «в среднем». Какой в среднем процент заявок будет теряться, какая средняя длительность ожидания обслуживания, какая в среднем нагрузка на один канал и т.д.
В то же время большой практический и теоретический интерес представляют модели нестационарных систем обслуживания (НСО). Изучению таких систем посвящено сравнительно мало работ.
Под моделью НСО понимается модель системы обслуживания с переменными во времени вероятностями состояний, каждое из которых определяет вероятность числа находящихся в НСО заявок и числа получивших обслуживание заявок. Число заявок, поступающих на обслуживание в рассматриваемых моделях, конечно. Стационарный режим не выделяется и основным предметом исследования является переходный процесс. Соответственно и важными характеристиками системы являются совсем другие: например среднее время обработки всех или некоторого числа заявок или вероятность не уложиться в требуемые сроки, вероятность превышения критической длины очереди, вероятность отказа в обслуживании критического числа заявок и т.д. Для иллюстрации важности моделирования НСО можно рассмотреть работу приемной комиссии университета. Заявления на поступление, приходят концентрированно в течение небольшого промежутка времени. Кроме того нас не интересует сколько из них мы можем потерять ведь важно обработать их все. Поэтому одной из важных характеристик будет являться количество людей (каналов обслуживания) необходимых для обработки заявлений в срок.
Имеется большое число реальных ситуаций, когда в потоке требований могут содержаться требования нескольких типов, причем требования первого типа обслуживаются вне всякой очереди, если только в очереди нет требований того же типа. По отношению к требованиям третьего и следующих типов правом преимущества пользуются требования второго типа и т.д. В качестве хорошего примера можно привести обслуживание на телеграфе, где срочные телеграммы передаются раньше обычных, даже когда простые телеграммы сданы ранее срочных.
Далее будут рассматриваться системы с так называемыми абсолютными приоритетами, в которых обслуживание заявки меньшего приоритета в случае прихода более «важной» заявки прерывается, и возобновляется лишь после обработки всех более приоритетных заявок.
Вероятности нахождения системы в том или ином состоянии есть функции от времени, которые являются решением обратного уравнения Чепмена – Колмогорова. Однако его составление даже для одноканальной НСО с абсолютными приоритетами представляется проблематичным.
В данной работе представлен алгоритм формирования нижнетреугольной матрицы коэффициентов линейной стационарной системы обыкновенных дифференциальных уравнений (ОДУ) для вышеуказанных систем без вывода общего уравнения системы. Для проверки корректности составления матрицы, проводится сравнение решения системы ОДУ с указанной матрицей коэффициентов, и результатами имитационного моделирования. Так же приведены графики, отображающие наиболее интересные характеристики исследуемой системы.
Постановка задачи
Рассматривается одноканальная нестационарная система массового обслуживания с абсолютными приоритетами. На вход системы поступают N заявок K типов (приоритетов) (K ≤ N). Известно сколько заявок каждого типа может прийти. Распределения временных интервалов между моментами поступления заявок описывается экспоненциальными законами с интенсивностями, зависящими от номера типа заявок. Закон распределения времён обслуживания – экспоненциальный с интенсивностями, так же зависящими от номера типа заявок. Вычислительную систему представим в виде одноканальной нестационарной системы обслуживания без потерь.
Необходимо сформулировать правила формирования списка состояний, а так же способ генерации нижнетреугольной матрицы коэффициентов системы однородных дифференциальных уравнений, описывающих вышеуказанную систему.
В данной работе описан алгоритм формирования списка состояний для НСО с абсолютными приоритетами таким образом, чтобы матрица A коэффициентов уравнения Чепмена – Колмогорова, описывающего данную систему, имела нижнетреугольный вид.
Правильность выполнения поставленной задачи показана, путем сравнения решения уравнения с результатами имитационного моделирования.
Для решения уравнения Чепмена – Колмогорова был реализован метод матричной экспоненты, имеющий очевидные недостатки. Одним из них является накапливающаяся величина погрешности округления, растущая вместе с увеличением рассматриваемого отрезка времени. Так же, в связи с разреженностью матрицы A, выбранный метод выполняет множество ненужных операций, что значительно увеличивает время работы программы. Чтобы максимально избежать описанных проблем, планируется использовать численно-аналитический метод, описанный в статье [1].
В дальнейшем планируется так же рассмотрение НСО с абсолютными приоритетами с несколькими каналами обслуживания, а так же с распределениями временных интервалов между моментами поступления и обработки заявок подчиняющимися гиперэкспоненциальному распределению, а так же распределениям Эрланга и Кокса.
Рост количества входных заявок увеличивает рассматриваемую систему на порядки, делая её крайне неудобной для исследования. В связи с этим возникает потребность в объединении состояний в группы по некоторому принципу, тем самым аппроксимируя исходную громоздкую систему. Выбор таких принципов так же входит в планы дальнейшего рассмотрения НСО с приоритетами.
1. Бубнов В. П., Еремин А. С., Сергеев С. А. Особенности программной реализации численно-аналитического метода расчета моделей нестационарных систем обслуживания // Труды СПИИРАН 2015. Вып. 38. С. 218-232.
2. Бубнов В. П., Сафонов В.И. Разработка динамических моделей нестационарных систем обслуживания. СПб.: Издательство «Лань». 1999. 64 с.
3. Бубнов В. П., Тырва А. В., Еремин А. С. Комплекс моделей нестационарных систем обслуживания с распределениями фазового типа // Труды СПИИРАН. 2014. Вып. 37. С. 61-71.
4. Вентцель Е.С. Теория вероятностей: Учебник для вузов. 6-е изд. стер. М.: Высш. шк., 1999. 576 c.
5. Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. Изд. 5-е, испр. М.: Издательство ЛКИ, 2011. 400с.
6. Кемени Д. Дж., Снелл Дж. Л. Конечные цепи Маркова. М.: Наука, 1970, 271 с.
7. Клейнрок Л. Теория массового обслуживания. Пер. с англ./ Пер. И. И. Грушко; ред. В. И. Нейман. М.:Машиностроение, 1979. 432с.
8. Матвеев В. Ф., Ушаков В. Г. Системы массового обслуживания. М.: Изд-во МГУ, 1984. 240 c.
9. Сафонов В. И. Моделирование систем обслуживания с изменяющейся интенсивностью поступления запросов. Деп. в войсковой части 11520 № В 216 // УПИМ. 1987. Вып. 3, сер, Б.
10. Хинчин А. Я. Работы по математической теории массового обслуживания / Под редакцией Б. В. Гнеденко. М.: Физматгиз, 1963. 236 с.