CC BY
13
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2019 Управление, вычислительная техника и информатика № 49
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 519.872
DOI: 10.17223/19988605/49/4
А.А. Назаров, Я.Е. Измайлова
АСИМПТОТИЧЕСКИЙ АНАЛИЗ СИСТЕМЫ С ПОВТОРНЫМИ ВЫЗОВАМИ, ВЫТЕСНЕНИЕМ ЗАЯВОК И ФАЗОВЫМ ДООБСЛУЖИВАНИЕМ
Проводится исследование системы с повторными вызовами, в которой входящий поток является простейшим, с вытеснением заявок и фазовым дообслуживанием методом асимптотического анализа. Использовалось предельное условие растущей задержки заявок на орбите. Применяя предлагаемый метод для данной системы с повторными вызовами (RQ-системы), показано, что характеристическая функция числа заявок на орбите является гауссовской. Найдены пропускная способность системы и параметры распределения вероятностей числа заявок на орбите.
Ключевые слова: вытеснение заявок; зоны орбиты; фазовое дообслуживание заявок; RQ-система; большая задержка.
В связи с тем, что в современном мире наблюдается ежегодное увеличение объема и размерности информационных данных, проблема проектирования и анализа телекоммуникационных систем является одним из приоритетных направлений развития любой страны. Система с повторными вызовами (RQ-система) - это модель массового обслуживания, которая характеризуется тем, что заявка, пришедшая в систему, если обнаружит прибор занятым, повторит попытку обращения для обслуживания через некоторое случайное время, в течение которого находится на орбите в зоне ожидания. В литературе показано, что такие модели являются адекватными математическими моделями телекоммуникационных сетей связи. Между повторами заявки (клиенты, требования) находятся на орбите. Первые системы с повторными вызовами были рассмотрены R.I. Wilkinson [1]. Наиболее полное и глубокое исследование разнообразных процессов в RQ-системах было приведено в работах Г.И. Фалина [2].
Одной из важных характеристик системы массового обслуживания является процесс обслуживания. В данной работе предлагается обслуживание в два этапа. На первом этапе прибор обрабатывает входящие заявки (первичная обработка), и после обслуживания на первой фазе заявки отправляются на вторую (фактическое обслуживание). Задачи, возникающие в системах управления и требующие двухфазного обслуживания, были предложены B.T. Doshi [3]. K.C. Madan [4] были рассмотрены системы массового обслуживания M/G/1, где заявка после обслуживания на первой фазе может покинуть систему или перейти на вторую фазу обслуживания. В работе B. Krishna Kumar [5] рассматривается система, в которой заявки, прибывшие в систему, вытесняют заявки, находящиеся только на первой фазе, и при этом обслуживающиеся в тот момент требования фактически покидают систему без обслуживания.
Как правило, рассматриваются системы, в которых все требования однотипные, хотя в реальных системах некоторые должны претендовать на первоочередное обслуживание, т.е. на приоритет. Системам с приоритетами посвящено большое количество исследований, например работы K. Avrachenkov [6], G. Ayyappan, A. Muthu Ganapathi [7], П.П. Бочарова, О.И. Павловой [8], I.M. Atencia [9], C. Kim [10], A.N. Dudin [11]. Литература по RQ-системам очень обширна. Некоторые обзоры моделей можно найти у J.R. Artalejo [12-13]. В данной работе будет рассмотрена RQ-система с повторными вызовами, фазовым дообслуживанием и вытеснением заявок.
1. Математическая модель и постановка задачи
Рассмотрим RQ- систему МЮ1/1 с вытеснением заявок и обслуживанием заново (рис. 1).
Орбита
О Q :а \
В(х)
Рис. 1. RQ-система M/GI/1 с вытеснением заявок Fig. 1. Retrial queueing system M/GI/1 with exclusion customers
Описание модели приведено в работе [14]. Показано, что пропускная способность RQ-системы с вытеснением заявок и обслуживанием их заново имеет вид Б = В'(0), где В(х) - функция распределения времени обслуживания. Откуда для времени обслуживания, распределенного по закону Эр-ланга, получаем Б = 0 . То есть стационарного режима в такой системе не существует.
Поэтому рассмотрим RQ-систему с вытеснением заявок и фазовым дообслуживанием (рис. 2).
На вход системы с повторными вызовами поступает простейший поток заявок, имеющий интенсивность X. Если пришедшая заявка застает прибор свободным, то она занимает его для обслуживания на первой фазе. Время обслуживания на первой фазе имеет экспоненциальное распределение с параметром Ц1. После успешного окончания обслуживания на первой фазе заявка переходит на вторую фазу обслуживания. Время обслуживания на второй фазе имеет экспоненциальное распределение с параметром Ц2. Если в момент прихода заявки прибор оказывается занятым, то заявка, которая пришла, вытесняет обслуживаемую заявку и встает на прибор. Заявка, вытесненная мгновенно, переходит на орбиту в зону ожидания. Орбита условно разделяется на две зоны. Если заявка обслуживалась на первой фазе, то она уходит в первую зону орбиты, если на второй фазе, то во вторую. В зонах орбиты заявки осуществляют случайную задержку. Продолжительность таких задержек имеет экспоненциальное распределение с параметром о. С орбиты после случайной задержки заявка вновь обращается к прибору с повторной попыткой обслуживания. С первой зоны орбиты обращается для обслуживания на первую фазу, со второй зоны - на вторую, таким образом, происходит фазовое до-обслуживание заявок, обслуживание которых было прервано.
Так как орбита разделяется на две зоны, то будем полагать:
стх = ухст - интенсивность обращения заявок с первой зоны орбиты;
ст2 = у2ст - интенсивность обращения заявок со второй зоны.
Обозначим \ (г) - число заявок в первой зоне орбиты; /2 (г) - число заявок на орбите во второй зоне, к (г) отвечает за состояние прибора: к (г) = 0, если прибор свободен; к (г) = 1, если прибор обслуживает заявку на первой фазе; к (г) = 2, если прибор обслуживает заявку на второй фазе.
Перед нами стоит задача нахождения двумерного распределения вероятностей числа заявок на орбите и состояний прибора.
Орбита
Яч- Чч-
ff ■ П ГТ ГТ
>а1 а, а, а.
Рис. 2. RQ-система M/E2/1 с вытеснением заявок и дообслуживанием Fig. 2. Retrial queueing system M/E2/1 with exclusion customers and phase follow-up
2. Уравнения Колмогорова
Рассмотрим трехмерный процесс |к(г0,^(0, 120:)} .
Обозначим стационарные вероятности Р{к(:) = к,',({) = 'т,'2(:) = '2} = Рк('1,'2), к = 0,1,2. Введем функции следующего вида:
да да ....
Ик(щ,и2) = Е Е е^1е"2'2Рк(ц,¡2), к = 0,1,2, ¡1 =0 '2 = 0
где у = V-т - мнимая единица, которые являются частичными характеристическими функциями.
Можно показать, что система уравнений для частичных характеристических функций имеет вид:
,„, ч . дИ0(иТ, и2) . 8И0(иТ, и2) -Ш0(щ, Щ) + -+ -^ 1 + И2 И2 ("т , Щ) = 0,
ощ ои2
-ц1И1(и1,и2) + ХИ0(иТ,и2) + \е]и2И2(ит,и2) + \е]и1 И1(и1,и2) - ХИ1(и1,и2) + уст2 и_
си 2
]Щ дИ0(и1,и2) Ш2 - и дИ 2 (и1, и2 ) _ п
уСТ1е _ уст1е _ = 0, дит дит
-ХИ2(ит,и2)- ц2И2(ит,и2) + цТИТ(ит,и2) + уст, дИ2) -
ди т
}и2 дИ0(и1,и2) р}"1-и2 дИ1(и1,и2) _ П
уст2е - уст2е - = 0-ди 2 ди 2
Данная система трудно решается аналитически, поэтому для исследования используем метод асимптотического анализа [15] в предельном условии большой задержки (ст ^ 0).
3. Асимптотический анализ
Введем обозначения:
кт - предельное среднее значение числа заявок на орбиты в первой зоне; к2 - предельное среднее значение числа заявок на орбите во второй зоне; Як , к = 0,1,2 - распределение вероятностей состояний прибора.
Теорема 1. Пусть ц ), '2 ^) - число заявок в зонах орбиты RQ-системы с вытеснением заявок и дообслуживанием. Тогда в предельном условии ст^ 0 для вероятностей состояний прибора Я = Р{к(г) = к} выполняются равенства
Я = 1----, Я = —, Я = —, (1)
Мт И2 Мт И2
а для характеристических функций выполняется предельное равенство
Ит М ехр{ ущст',^) + уи2ст'2(:)} = ехр{ у^к, + уи2к2},
где параметры к , к2 имеют вид:
Л 2 Л 2
к = Л М 2 , _^Ит_. (2)
1 У1 (М1И2 -ЧМт +И2Ж 2 У2 (ИтИ2 -Х(М1 +М2))
Из полученных в теореме 1 равенств (2) допредельные средние значения числа заявок в зонах
к
на орбите можно аппроксимировать величинами — , к = 1,2 .
ст
В работе [14] дано определение пропускной способности.
Следствие. Из первого равенства в (1) теоремы 1 получаем значение пропускной способности RQ-системы с вытеснением заявок и фазовым дообслуживанием в виде:
5 = > 0.
M1 + M2
Поэтому стационарный режим в такой системе существует для любых X < S .
Теорема 2. Пусть ix (t), /2 (t) - число заявок в зонах орбиты RQ-системы с вытеснением заявок и фазовым дообслуживанием. Тогда имеет место предельное равенство
lim M exp | jui + л } = exp {Ц1»2 Qn + С/^2)! „22 + },
где к , к2 определены выше, величины Qn, Q12, Q22 - решение неоднородной системы уравнений ßl 1 ( Y1R + Y1R - Y1Rz0 + Y1Rz1 + Y1Rz1 - Y1Rz2 ) + Q12 (—Y2R1 — Y2Rz0 + Y2Rz2 — Y2R1z1 + Y2R1z2 ) =
=1XR1 + 1 YlK1R0 + 1R2YlK1 + 1 RlY2K2 + YlK1R0zi -XR1Z1 +YlK1R2zi — Y2K2R1z2> ß22 (Y2R0 + Y2R1 — Y2RУо + Y2RУ2 + Y2Ry — Y2R^1) +ß12 (—YlR2 — YiRoУ0 + YiRoУ1 — YiR2У2 + YlR2У1) = = 1XR2 + 1Y 2K2 R0 + 1R2 YlK1 + 1 RlY 2K2 — XR2 У1 +Y 2K2 R0 У2 +Y 2 K2 R^2 — YlK1R2 Уь
Q11 (—Y1R2 — Y1R0P0 + Y^A + Y1R2P1 — Y1R2P2 ) + Q12 (YR + Y1R2 + Y2r + Y2R — Y2R0P0 +
+Y 2 R0 P2 — Y 2 R1P1 + Y 2 R1P2 — Y1R S0 + YR S1 — Y1R2S2 + Y1R2S1) +
+Q22 (—Y2R1 — Y2R0s0 + Y2R0s2 + Y2R1s2 — Y2R1s1) = —R2Y1K1 — R1Y2K2 + (3)
+Y1K1R0P1 — ЩП + Y1K1R2P1 — Y 2K2 R1P2 — XR2S1 + Y 2K2 -^0^2 + Y 2 K2 R1S2 — Y1K1R2 S1-Величины zo, zi, Z2; yo, yi, yi, po, pi, pi, so, si, S2 являются решением систем уравнений (4)-(7) соответственно.
Для zo, zi, Z2:
— (X + Y1K1 + Y2k2 ) z0 + (X + Y1K1) z1 + Y2k2 z2 = —Y1K1, — (M1 + Y2k2 ) z1 + (М + Y2K2 ) z2 = X + Y2K2>
M2 z0 + (X + Y1K) z1 — (X + M2 + Y1K1) z2 = —Y1K1 • (4)
Для yo, yi, У2:
— (X + Y1K +Y2k2 ) У0 + (X + Y1K) У1 +Y2K2y2 =—Y2k2. —(M1 +Y2k2 ) У + (М +Y2k2 ) У2 =—Y2K2.
Ц2У0 + (X + Y1K1) У1 — (X + M2 + Y1K1) У2 = X + Y1K1 • (5)
Для po, pi, P2:
—(X + Y1K +Y2k2 ) P0 + (X + Y1K) P1 + Y2K2P2 = —Y2k2. — (M1 +Y2k2 ) P1 + (М +Y2k2 ) P2 = —Y2K2.
M2 P0 + (X + Y1K1) P1 — (X + M2 + Y1K1) P2 = X + Y1K1 • (6)
Для So, Si, S2:
— (X + Y1K1 +Y2k2 ) s0 +(X + Y1K1) s1 +Y2K2S2 = —Y1K1, — (M1 +Y2k2 ) s1 +(М +Y2K2 ) s2 = X + Y2K2>
M2s0 + (X + Y1K1) S — (X + M2 + Y1K1) s2 = —Y1K1 • (7)
Запишем характеристическую функцию распределения вероятностей числа заявок на орбите в следующем виде:
H и2) = exp { А k1 + k2 + Qn + Q22 + jj Q12 }•
I ст ст 2ст 2ct ct i
Из вышеизложенного получаем, что двумерное распределение вероятностей числа заявок на орбите является в пределе гауссовским распределением.
Заключение
В работе была исследована система с повторными вызовами, фазовым дообслуживанием и вытеснением заявок методом асимптотического анализа. Показано, что предельная характеристическая функция числа заявок на орбите имеет вид характеристической функции гауссовского распределения. Получено значение пропускной способности данной системы и сделан вывод о том, что при значениях параметра интенсивности входящего потока Х< S для RQ-системы с вытеснением заявок и фазовым дообслуживанием стационарный режим существует, в то время как для RQ-системы с вытеснением заявок и обслуживанием заново стационарного режима не существует. Считаем целесообразным исследовать RQ-систему с вытеснением заявок и фазовым дообслуживанием общего типа, т.е. когда присутствует N фаз.
ЛИТЕРАТУРА
1. Wilkinson R.I. Theories for toll traffic engineering in the USA // The Bell System Technical Journal. 1956. V. 35, No. 2. P. 421-
507.
2. Falin G.L., Templeton J.G.C. Retrial queues. London : Chapman & Hall, 1997. 328 р.
3. Doshi B.T. Analysis of a two phase queueing system with general service times // Operations Research Letters. 1991. V. 10, No. 5.
P. 265-272.
4. Madan K.C. An M/G/1 queue with second optional service // Queueing Systems. Theory and Applications. 2000. V. 3, No. 4.
P. 37-46.
5. Krishna K.B., Vijaykumar A., Arivudainambi D. An M/G/1 Retrial queueing system with two phase service and preemptive
resume // Annals of Operations Research. 2002. V. 113, No. 1-4. P. 61-79. D0I:10.1023/A:1020901710087.
6. Avrachenkov K., Dudin A., Klimenok V. Queueing model MMAP/M 2/1 with two orbits // Lecture Notes in Comput. Sci. 2010.
V. 6235. P. 107-118.
7. Ayyappan G., Muthu G.A., Gopal S. M/M/1 Retrial queuing system with loss and feedback under pre-emptive priority service //
Int. J. of Computer Applications. 2010. V. 2, No. 6. P. 27-34.
8. Bocharov P.P., Pavlova O.I., Puzikova D.A. M| G| 1 |r retrial queueing systems with priority of primary customers // Mathematical
and computer Modelling. 1999. V. 30, No. 3-4. P. 89 -98.
9. Atencia I.M. A Geo/G/1 retrial queueing system with priority services // European J. of Operational Research. 2017. Vol. 256,
No. 1. P. 178-186.
10. Kim C. Priority tandem queueing system with retrials and reservation of channels as a model of call center // Computers & Industrial Engineering. 2016. V. 96. P. 61-71.
11. Dudin A.N., Kim C., Dudin S., Dudina O. Priority retrial queueing model operating in random environment with varying number and reservation of servers // Applied Mathematics and Computation. 2015. V. 269. P. 674-690.
12. Artalejo J.R. Aceessible bibliography on retrial queues // Mathematical and Computer Modelling. 1999. V. 30. P. 1-6.
13. Artalejo J.R. Accessible bibliography on retrial queues: Progress in 2000-2009 // Mathematical and Computer Modelling. 2010. V. 51, No. 9-10. P. 1071-1081.
14. Назаров А.А., Черникова Я.Е. Исследование RQ-системы M / GI / 1 с вытеснением в условии большой задержки // Известия Томского политехнического университета. 2013. Т. 323, № 5. С. 16-20.
15. Назаров А.А., Моисеева С.П. Методы асимптотического анализа в теории массового обслуживания. Томск : Изд-во НТЛ, 2006. 112 c.
Поступила в редакцию 17 мая 2019 г.
Nazarov A.A., Izmailova Y.E. (2019) ASYMPTOTIC ANALYSIS RETRIAL QUEUEING SYSTEM WITH EXCLUSION CUSTOMERS AND PHASE FOLLOW-UP. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie vychislitelnaja tehnika i informatika [Tomsk State University Journal of Control and Computer Science]. 49. pp. 29-34
DOI: 10.17223/19988605/49/4
In this paper, we consider a retrial queueing system, which has a two-phase serviced and displacement of customers. Frequently such tasks arise in control systems. The input of the system is a simplest request flow with intensity X. If the received customer finds the server free, it takes it for service in the first phase. The service time in the first phase has an exponential distribution with a parameter Ц1. After successful completion of service in the first phase, the request goes to the second phase of service. Service time in the second phase has an exponential distribution with a parameter Ц2. If the server is busy, then the incoming request pushes the previous request out of service and takes its place, while the pushed out request goes to the orbit. The orbit is conventionally divided into two zones. If the customer was serviced in the first phase, it goes to the first zone of the orbit, but if it was serviced in the second
phase, then it comes in the second phase of the orbit. In the orbit zones, the requests are subject to a random delay. The duration of such delays has an exponential distribution with a parameter a. After a random delay the customer re-accesses the server from an orbit to the re-attempt service. From the first zone of the orbit, the customer is appealed for service to the first phase, from the second zone, one is appealed to the second phase, thus wise, the phase maintenance of customers whose service has been interrupted occurs. The throughput of RQ-system with the displacement of customers and phase follow-up. By throughput we mean the maximum average number of requests that the system can serve per unit of time. In this RQ-system it has the form S = > 0 , so the stationary mode in such a system exists for any X < S. Further, the observation of the system is carried out using the method of asymptotic analysis. It was shown that the limit characteristic function of the number of customers in orbit has the form of a characteristic function of the Gaussian distribution. Numerical experiments were carried out. The asymptotic probability distribution of the number of customers in orbit was obtained for the selected parameter values.
Keywords: exclusion of customers; zones in the orbit; phase follow-up customers; retrial queueing system; the long delay.
NAZAROVAnatoly Andreevich (Doctor of Technical Sciences, Professor, Head of the department of probability Theory and mathematical statistics, National Research Tomsk State University, Tomsk, Russian Federation). E-mail: nazarov.tsu@gmail.com
IZMAYLOVA Yana Evgenevna (Candidate of Physics and Mathematics, Associate Professor, National Research Tomsk State University, Tomsk, Russian Federation). E-mail: evgenevna.92@mail.ru
REFERENCES
1. Wilkinson, R.I. (1956) Theories for toll traffic engineering in the USA. The Bell System Technical Journal. 35(2). pp. 421-507.
DOI: 10.1002/j.1538-7305.1956.tb02388.x
2. Falin, G.L. & Templeton, J.G.C. (1997) Retrial queues. London: Chapman & Hall.
3. Doshi, B.T. (1991) Analysis of a two phase queueing system with general service times. Operations Research Letters. 10(5).
pp. 265-272. DOI: 10.1016/0167-6377(91)90012-E
4. Madan, K.C. (2000) An M/G/1 queue with second optional service. Queueing Systems. Theory and Applications. 3(4). pp. 37-46.
DOI: 10.1023/A:1019144716929
5. Krishna, K.B. & Vijaykumar, A. & Arivudainambi, D. (2002) An M/G/1 Retrial queueing system with two phase service and
preemptive resume. Annals of Operations Research. 113(1-4). pp. 61-79. DOI: 10.1023/A:1020901710087
6. Avrachenkov, K. & Dudin, A. & Klimenok, V. (2010) Queueing model MMAP/M 2/1 with two orbits. Lecture Notes in Comput-
er Science. 6235. pp. 107-118.
7. Ayyappan, G. & Muthu Ganapathi, A. & Sekar, G. (2010) M/M/1 Retrial queuing system with loss and feedback under pre-
emptive priority service. International Journal of Computer Applications. 2(6). pp. 27-34. DOI: 10.5120/672-945
8. Bocharov, P.P., Pavlova, O.I. & Puzikova, D.A. (1999) M| G| 1 |r retrial queueing systems with priority of primary customers.
Mathematical and computer Modelling. 30(3-4). pp. 89-98. DOI: 10.1016/S0895-7177(99)00134-X
9. Atencia, I.M. (2017) A Geo/G/1 retrial queueing system with priority services. European Journal of Operational Research.
256(1). pp. 178-186. DOI: 10.1016/j.ejor.2016.07.011
10. Kim, C. (2016) Priority tandem queueing system with retrials and reservation of channels as a model of call center. Computers & Industrial Engineering. 96. pp. 61-71. DOI: 10.1016/j.cie.2016.03.012
11. Dudin, A.N., Che Soong Kim, Dudin, S. & Dudina, O. (2015) Priority retrial queueing model operating in random environment with varying number and reservation of servers. Applied Mathematics and Computation. 269. pp. 674-690. DOI: 10.1016/j.amc.2015.08.005
12. Artalejo, J.R. (1999) Accessible bibliography on retrial queues. Mathematical and Computer Modelling. 30. pp. 1-6. DOI: 10.1016/S0895-7177(99)00128-4
13. Artalejo, J.R. (2010) Accessible bibliography on retrial queues: Progress in 2000-2009. Mathematical and Computer Modelling. 51(9-10). pp. 1071-1081. DOI: 10.1016/j.mcm.2009.12.011
14. Nazarov, A.A. & Chernikova, Y.E. (2013) Investigation of preemptive RQ-system M| GI| 1 under extensive delay. Izvestiya Tomskogo politekhnicheskogo universiteta - Bulletin of the Tomsk Polytechnic University. vol. 323. no 5. pp. 16-20.
15. Nazarov, A.A. & Moiseeva, S.P. (2006) Metody asimptoticheskogo analiza v teorii massovogo obsluzhivaniya [The asymptotical analysis method in queueing theory]. Tomsk: NTL.
