CC BY
54
УДК 53.082
Система массового обслуживания с отрицательными заявками, бесконечным накопителем и конечным бункером для вытесненных заявок в непрерывном времени
Р. В. Разумчик
Кафедра теории вероятностей и математической статистики Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, 117198, Москва, Россия
Рассматривается система массового обслуживания, в которую поступают пуассонов-ские потоки обычных и отрицательных заявок. Для обычных заявок имеется накопитель неограниченной ёмкости. Отрицательная заявка при поступлении в систему вытесняет одну обычную заявку из очереди в накопителе и перемещает её в бункер ограниченной ёмкости. Если в момент вытеснения заявки бункер полностью заполнен, то вытесненная заявка теряется.
Если в момент прихода отрицательной заявки накопитель пуст, она покидает систему, не оказывая на неё никакого воздействия. После окончания обслуживания очередной заявки на прибор поступает заявка из накопителя или, если накопитель пуст, из бункера. Длительности обслуживания заявок, как из накопителя, так и из бункера имеют экспоненциальное распределение с одним и тем же параметром. Получены формулы для расчёта совместного стационарного распределения числа заявок в накопителе и бункере и основных вероятностных характеристик системы.
Ключевые слова: система массового обслуживания, отрицательная заявка, бункер.
1. Введение
В последние несколько лет значительное внимание уделяется изучению систем и сетей с отрицательными заявками, что во многом мотивировано их практической ценностью. Впервые введённая в работе [1] отрицательная заявка при поступлении в систему «убивает» (разрушает) одну или несколько обычных заявок, ожидающих в очереди, после чего мгновенно покидает систему вместе с «убитыми» заявками. С помощью подобных систем можно моделировать такие эффекты как потери информации в инфотелекоммуникационных системах в связи с проникновением вирусов, сбоем программного обеспечения или внезапным выходом их строя оборудования в распределённых базах данных и др.
Одной из интересных задач является исследование систем с отрицательными заявками несколько иного типа, по сравнению с тем, который был предложен в [1]. Новый тип предполагает, что отрицательная заявка в момент поступления в систему не «убивает», а вытесняет одну из заявок в основной очереди и помещает её в дополнительную очередь, заявки из которой могут обслуживаются, но только с относительным приоритетом.
В работах [2] и [3] были исследованы подобные системы с бесконечным числом ожиданий мест в обеих очередях и функционирующие как в непрерывном, так и в дискретном времени. Были получены формулы для расчёта совместного стационарного распределения числа заявок в очередях и основных вероятностных характеристик.
Настоящее исследование является развитием работы [2] на случай, когда одна из очередей, а именно дополнительная очередь, имеет ограниченную ёмкость.
Итак, рассматривается однолинейная СМО, в которую поступает пуассонов-ский поток заявок интенсивности Л. Заявки этого потока будем называть положительными. Для положительных заявок имеется накопитель положительных заявок (далее просто накопитель) неограниченной ёмкости.
Статья поступила в редакцию 31 октября 2010 г.
Работа выполнена при финансовой поддержке Российского Фонда Фундаментальных Исследований (проект 09-07-12032).
Помимо положительных заявок, в систему поступает пуассоновский поток отрицательных заявок интенсивности Л-. Отрицательная заявка, поступающая в систему, вытесняет одну (положительную) заявку из конца очереди в накопителе и перемещает её в накопитель для вытесненных заявок или бункер, который имеет ёмкость г (г < то).
Если в момент вытеснения заявки из очереди в накопителе бункер полностью заполнен, то вытесненная заявка теряется.
Если в момент поступления отрицательной заявки в накопителе нет положительных заявок, а на приборе обслуживается заявка, то отрицательная заявка, не прерывая обслуживание на приборе, покидает систему, не оказывая на неё никакого воздействия. То же самое происходит и в случае, когда в момент поступления отрицательной заявки накопитель и обслуживающий прибор пусты.
Выбор заявок на обслуживание производится следующим образом. После окончания обслуживания очередной заявки на прибор становится заявка из накопителя. Если же накопитель пуст, на прибор поступает заявка из бункера. Обслуживание заявок не прерывается новыми как положительными, так и отрицательными заявками.
Длительности обслуживания заявок как из накопителя, так и из бункера имеют экспоненциальное распределение с одним и тем же параметром ц.
2. Система уравнений равновесия
Обозначим через и(¿) число заявок, находящихся в накопителе в момент времени а через — число заявок в бункере. Положим X(¿) = (^(¿),^(£)). Случайный процесс {X(¿), £ ^ 0}, описывающий стохастическое поведение рассматриваемой СМО во времени, является марковским процессом (МП) с непрерывным временем и дискретным счётным множеством состояний. Множество состояний процесса {X(¿), £ ^ 0} является объединением X = и?=0Х^ непересекающихся подмножеств Хг, которые далее мы будем называть слоями. Здесь Хо = {0}, = {(г,3), г ^ 0, 0 < ] < г}. Состояние (0) слоя Хо соответствует полностью свободной системе; состояние (г,^) слоя Х^ означает, что в накопителе находится % заявок, а в бункере ожидают '} заявок, вытесненных из накопителя.
Инфинитезимальная матрица переходов данного процесса является блочной трёхдиагональной вида
Я =
(N0 Л1 0 0 0
Мо Л 0 0
0 М N Л 0
0 0 М N Л
\.
Л
(1)
.
где квадратные матрицы Л, N, М размерности (г + 1) х (г + 1) имеют следующий вид:
(Х 0 0 . .0 0 ( № 0. .. 0 0 \
0 \ 0 . .0 0 0 \- . .. 0 0
Л= 0 0 Л . .0 0 , м = 0 0 ^ . .. 0 0
0 0 0 . . л 0 0 0 0. .. V X-
0 0 0 . .0 X/ 0 0 0. .. 0 ц + А" 7
/-0" + А) 0
I! -(^ + Л)
N1 =
0
0 0
-(/! + Л)
V
0 0
0 0
0 0 0
-(и + Л) 0 I! -(^ + А)/
(4)
и N0 = (-А), Л1 = (Л, 0, 0, ... 0, 0), М0 = (и, 0, 0, ... 0, 0)т.
Матрица N, стоящая на пересечении к-й строки и к-го столбца матрицы Я, описывает переходы МП между состояниями внутри слоя к. Матрица Л, стоящая на пересечении (к - 1)-й строки и к-го столбца матрицы Я, описывает переходы МП из состояний слоя (к - 1) в состояния слоя к. Матрица М, стоящая на пересечении (к + 1)-й строки и к-го столбца матрицы Я, описывает переходы МП из состояний слоя (к + 1) в состояния слоя к.
Обозначим через рт = (ро,рТ,Р\,...) вектор стационарных вероятностей МП {X(£), £ > 0}, где р!Т = (рго,рг1,... ,Ргг) является вектором размерности г + 1. Тогда из системы уравнений равновесия (СУР)
рт Я = 0т
следует, что компоненты вектора р удовлетворяют уравнениям
роМо + рТ Мо = 0, р0 Л1 + рТ N1 + рТ М = 0Т,
рТ- 1Л + рТ N + рТ+1М = 0Т, г > 1,
с условием нормировки
Ро + ^ Рг 1 = 1,
г=о
(5)
(6)
(7)
(8)
(9)
где 1 — вектор-столбец из единиц размерности г + 1.
Существование стационарного распределения для рассматриваемой системы и решение СУР (6)-(9) следует из следующей теоремы [4, стр. 32, теорема 1.7.1.].
Теорема 1. Неприводимый МП {X(£), £ > 0} является положительно возвратным тогда и только тогда, когда минимальное неотрицательное решение К уравнения
Л + ЯМ + Я2М = 0 (10)
имеет спектральный радиус р(Я) < 1, и если существует неотрицательный вектор (р0,р'Т) такой, что выполняется
р0Х0 + рТ Мо = 0,
(11)
ро Л1 + (N1 + ЕМ ) = 0 Т. Тогда решение СУР (6)-(9) имеет вид
р? = р0ТЕг, г > 0,
р 0Т = -р0 Л1(^1 + КМ)
-1
1
Ро =
1 - Л1(^1 + ЕМ)-1(1 - К)-1 1' где I — единичная матрица размерности (г + 1) х (г + 1).
(12)
(13)
(14)
(15)
Покажем, что условия теоремы выполняются.
Нетрудно видеть, что все состояния МП {X(¿), £ ^ 0} с инфинитезимальной матрицей сообщаются между собой, поэтому процесс является неприводимым.
В [5] показано, что достаточным условием того, что максимальное по модулю собственное значение (спектральный радиус) р(Щ имел алгебраическую кратность 1 и все остальные собственные значения были по модулю меньше р(Щ, является неразложимость матрицы .
Однако в данном случае достаточное условие не выполняется, потому что матрица Я является разложимой. В самом деле, если рассмотреть матрицу А:
А=Л+N+М=
(-Х- X- 0. .0 0
0 А- . .0 0
0 0 -Л- . .0 0
0 0 0. . -А- А-
0 0 0. .0 0
(16)
то видно, что она является разложимой, а значит разложимой является и матрица Тем не менее вид матриц Л, N и М позволяет убедиться непосредственной проверкой, что спектральный радиус р(К) имеет алгебраическую кратность 1 и все остальные собственные значения по модулю меньше р(Щ. Покажем это.
Рассматривая уравнение (10) в скалярном виде, нетрудно убедиться, что матрица К является верхней треугольной:
К =
(г0 о г 01 Г02 . .. Г0,г-1 Г0,г\
0 Г11 Г12 . .. Г1,г-1 Г1,г
0 0 Г22 . .. Г2,г-1 Г2,г
0 0 0. .. Гз,г-1 Гз,г
0 0 0. .. 0 ^ т т У
(17)
Следовательно, все собственные значения матрицы К — это её диагональные элементы rjj, 0 ^ у ^ г. Таким образом, матрица Л + НХ + В?М в левой части уравнения (10) является верхней треугольной. Введём производящие функции (ПФ)
Фз (г) = (ЛЬ- + (XX + (МЬ г2, 0 < з < г,
(18)
где через обозначается элемент, стоящий на пересечении j-й строки и j-го столбца соответствующей матрицы. Тогда условие (10) равносильно тому, что
Ъ (Гц) = 0, 0 < з < г. (19)
Если переписать данные уравнения, учитывая вид матриц Л, N и М, то для
j = 0,г — 1 все уравнения будут одинакового вида
rjj +
X — (А + р + \-)rjj + = 0, (20)
а при j = г имеем:
Л — (Л + р + A-)fjj + (р + \-)r2¡j = 0. (21)
Решая эти уравнения, получаем, что минимальное неотрицательное решение уравнения (10) со спектральным радиусом р(К) < 1 существует, причём собственные значения матрицы К определяются следующим образом:
Л + ц + А- - -у/(А + ц + А-)2 - 4\р
2/j,
А + р + А- - у/(Л + р + А-)2 - 4А(и + А-) _ Л
rjj = <......- „Г______________ , (22)
j = 0, г — 1,
3 = г.
2(м + А-) р + Л"
Отсюда видно, что спектральный радиус p(R) = rrr.
Как показано в [4, см. 1.7.11, с. 32], спектральный радиус матрицы R меньше 1 тогда и только тогда, когда
жт Л1 <жт М1, (23)
где ж — единственное положительное решение СУР
жт А = 0т (24)
с условием нормировки
^т1 = 1. (25)
Если матрица А неразложима или разложима, но имеет единственную неразложимую группу индексов, то вектор ж, удовлетворяющий условиям (24) и (25), существует и единственен, и в данном случае имеет вид жт = (0,0,0... 0,1).
Окончательно из (23) получаем, что условие существования стационарного режима выглядит следующим образом: Л < (р + Л-).
Последнее условие теоремы, которое осталось проверить — это то, что существует неотрицательный вектор (р0,рт) вида (14) и (15), удовлетворяющий уравнениям (11) и (12). Проверить это можно непосредственной подстановкой (14) и (15) в СУР. Неотрицательность вектора (pq,Pq ) вытекает из того, что согласно [6] (см. Предложение 1. стр. 366) матрица (I — R)-1 > 0, а матрица N1 + RM имеет отрицательные диагональные элементы (и неотрицательные недиагональные) и, значит, (N1 + RM)-1 < 0.
Таким образом, все условия теоремы 1 выполняются и стационарное распределение задаётся формулами (13)-(15).
Зная стационарное распределение вероятностей состояний, выпишем некоторые вероятностно-временные характеристики системы.
Пусть éi — вектор-столбец размерности (г +1), у которого все элементы равны 0, кроме элемента на г-м месте.
Тогда, вероятность pioss потери положительной заявки вследствие переполнения очереди в бункере имеет вид
Vioss = -у l^Pi er+1 = ~yP0 R(I — R) er+1.
i=1
Среднее число Мпак и дисперсия числа заявок в накопителе рассчиты-
ваются по формулам:
о
мпак = £^1= р^Е(1 - В)-2 1, Впак = 2рТК2(1 - К)-3 1 + Мпак(1 - Мпак). Среднее
число и дисперсия числа заявок в бункере можно найти с
помощью следующих выражений:
оо г г г
Мъип = ^р?^ = Ро (/я^+ь БЪип = (1 Р^+1-Мъип.
г=0 г=1 г=1 г=1
Формула для вычисления среднего числа N заявок в системе имеет вид:
N =1 - ро + МПак + МЪип.
Для среднего времени ш пребывания заявки в накопителе до её ухода либо на прибор, либо в бункер, либо из системы вследствие полностью заполненной очереди в бункере, используя формулу Литтла [7], имеем следующее выражение:
ш = -1 К(1 - К)-21.
3. Заключение
Таким образом, для рассматриваемой системы были получены формулы для расчёта совместного стационарного распределения числа заявок в накопителе и бункере, а также выражения для ряда основных вероятностных характеристик системы, включая вероятность потери заявки из-за переполнения очереди в бункере.
Отдельно отметим, что в формулах расчёта стационарного распределения (13)-(15) участвует матрица К, которая имеет вид (17). Вопрос нахождения матрицы К является хорошо изученным. Для непосредственного расчёта матрицы К можно воспользоваться алгоритмами, предложенными в [7-15].
Для проверки аналитических методов была разработана имитационная модель с помощью программных средств СРЯБ, которая показала хорошее совпадение с результатами численных расчётов.
Литература
1. Gelenbe E., Glynn P., Sigman K. Queues with Negative Arrivals // J. Appl. Probab. — 1991. — Vol. 28. — P. 245-250.
2. Мандзо Р., Касконе Н., Разумчик Р. В. Экспоненциальная система массового обслуживания с отрицательными заявками и бункером для вытесненных заявок // Автоматика и телемеханика. — 2008. — № 9. — С. 103-113. [Mandzo R., Kaskone N., Razumchik R. V. Ehksponencialjnaya sistema massovogo obsluzhivaniya s otricateljnihmi zayavkami i bunkerom dlya vihtesnennihkh zayavok // Avtomatika i telemekhanika. — 2008. — No 9. — S. 103113.]
3. Печинкин А. В., Разумчик Р. В. Система массового обслуживания с отрицательными заявками и бункером для вытесненных заявок в дискретном времени // Автоматика и телемеханика. — 2009. — № 12. — С. 109120. [Pechinkin A. V., Razumchik R. V. Sistema massovogo obsluzhivaniya s otricateljnihmi zayavkami i bunkerom dlya vihtesnennihkh zayavok v diskretnom vremeni // Avtomatika i telemekhanika. — 2009. — No 12. — S. 109-120.]
4. Neuts M. F. Matrix-Geometric Solutions in Stochastic Models: An Algorithmic Approach. — The John Hopkins University Press, 1981.
5. Latouche G, Taylor P. G. Level-Phase Independence for GI\M|1 Type Markov Chains // J. Appl. Probab. — 2000. — Vol. 37, No 4. — Pp. 984-998.
6. Гантмахер Ф. М. Теория матриц. — М., 1966. [Gantmakher F. M. Teoriya matric. — M., 1966.]
7. Бочаров П. П., Печинкин А. В. Теория массового обслуживания. — М.: Изд-во РУДН, 1995. [Bocharov P. P., Pechinkin A. V. Teoriya massovogo obsluzhivaniya. — M.: Izd-vo RUDN, 1995.]
8. Favati P., Meini B. On Functional Iteration Methods for Solving Nonlinear Matrix Equations Arising in Queueing Problems // IMA J. Numer. Anal. — 1999. — No 19. — Pp. 39-49.
9. Guo C.-H. On the Numerical Solution of a Nonlinear Matrix Equation in Markov Chains // Linear Algebra Appl. — 1999. — No 288. — Pp. 175-186.
10. Latouche G. Newton's Iteration for Non-Linear Equations in Markov Chains // IMA J. Numer. Anal. — 1994. — No 14. — Pp. 583-598.
11. Latouche G. Algorithms for Evaluating the Matrix G in Markov Chains of Ph\G\1 Type // Cahiers Centre Etudes Rech. Oper. — 1994. — No 36. — Pp. 251-258.
12. Meini B. New Convergence Results on Functional Iteration Techniques for the Numerical Solution of M|G|1 Type Markov Chains // Numer. Math. — 1997. — No 78. — Pp. 39-58.
13. Neuts M. F. Moment Formulas for the Markov Renewal Branching Process // Adv. in Appl. Probab. — 1976. — No 8. — Pp. 690-711.
14. Ye Q. High Accuracy Algorithms for Solving Nonlinear Matrix Equations in Queueing Models // Advances in Algorithmic Methods for Stochastic Models. — 2000. — Pp. 401-415.
15. Latouche G., Ramaswami V. A Logarithmic Reduction Algorithm for Quasi-Birth-Death Processes // J. Appl. Probab. — 1993. — No 30. — Pp. 650-674.
UDC 53.082
A Queueing System with Negative Claims, Infinite Buffer and Finite Bunker for Superseded Claims
R. V. Razumchik
Department of Probability Theory and Mathematical Statistics Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, 117198, Russia
Study is performed of the queuing system with Poisson incoming flows of ordinary and negative claims. For ordinary claims, there is an infinite-capacity buffer. The arriving negative claim knocks out an ordinary claim queued in the buffer and moves it to a finite-capacity buffer of superseded claims (bunker). If upon the knock out of the ordinary claim bunker is full, knocked out claim is lost. If the buffer is empty, then the negative claim discharges the system without affecting it. After servicing the current claim, the server receives a claim from the buffer or, if the buffer is empty, the bunker. The claims arriving from both the buffer and bunker are served exponentially with the same parameter.
Expressions for the stationary joint probability distribution of number of claims in buffer and bunker and for main probabilistic characteristics of the system were obtained.
Key words and phrases: queueing system, negative claim, bunker.
