CC BY
14ASIMPTOTIC ANALYSIS OF TWO-PHASE RETRIAL QUEUING SYSTEM M/M/l UNDER DELAY CONDITION OF ORBIT
A. A. Nazarov, A. A. Anisimova
Tomsk State University, 634050, Tomsk, Russia
In this paper, we consider a M/M/l retrial queuing system with two phases of essential service. Customers arrive at the system according to a Poisson process with arrival rate A. Service times follow an exponential distribution with parameters ¡¿i and //2 for the first and the second phase respectively. If a customer finds the server free, it is served immediately; otherwise the customer joins the orbit. Customers in an orbit come back to be served after intervals of time exponentially generated with parameters o\ and 02 for the first and the second phase respectively. If the server is free, it will receive the service; otherwise it will be re-inserted into the orbit. Each phase includes one server and one orbit. The aim of this paper is to obtain the mean values of the number of customers in the orbits using the asymptotic analysis method under delay condition of orbit. The results of asymptotic analysis and the results of simulation are compared.
Key words: retrial queuing system, two-phase system, Poisson process, delays condition of
orbit.
References
1. Art ale jo J. R., Gomez-Corral A. Retrial Queuing Systems: A Computational Approach. Berlin, Springer, 2008.
2. Falin G. I., Templet on .J. G. C. Retrial queues. London, Chapman & Hall, 1997.
3 Gnedenko B. V., Kovalenko I. N. Vvedenie. v te.oriyu massovogo obsluzhivaniya [Introduction to Queuing Theory]. Moscow, Nauka, Gl. red. fiz.-mat.lit., 1987.
4. Krishna C. M. and Lee Y. H. A study of two-phase, service.. Operations Research Letters, 1990, Vol. 9. P. 91 97.
5. Doshi, B. T. Analysis of a two phase, queuing system with general service times /7 Operations Research Letters. 1991. Vol. 10. P. 265 272.
6. Krishnakumar B. and Arivudainambi D. An M/G/l retrial queuing system with two phases of services and preemptive, resume. /7 Annals of Operations Research. 2002. Vol. 113. P. 61 79.
7. Choudhurv G. A single server queuing system with two-phases of service and vacations. /7 QTQM. 2008. Vol. 5 (1). P. 33 49.
8. Nazarov A. A., Sudvko E. A. Method of asymptotic, semi invariants for studying a mathematical model of a random access network /7 Problems of information transmission. 2010. Vol. 46. N 1. P. 86 102.
9. Anisimova A. A. Imitatsionnoe. modelirovanie. dvukhfaznoi RQ-siste.my [Simulation modeling of two-phase retrial queuing system] /7 Nauchnoe tvorchestvo molodezhi. Matematika. Informatika: materialv XX Vserossiiskoi nauchno-prakticheskoi konferentsii (28 29 aprelva 2016 g.), sost . Y. A. Naumkina. Tomsk, Izd-vo Tom. un-ta, 2016, Part 1, P. 81 85.
АСИМПТОТИЧЕСКИЙ АНАЛИЗ ПЕРВОГО ПОРЯДКА ДВУХФАЗНОЙ ШЭ-СИСТЕМЫ М/М/1 В УСЛОВИИ БОЛЬШОЙ ЗАДЕРЖКИ В ИСТОЧНИКАХ ПОВТОРНЫХ
ВЫЗОВОВ
А. А. Назаров, А. А. Анисимова
Национальный исследовательский Томский государственный университет
634050, Томск, Россия
УДК 519.872
В данной статье рассматривается двухфазная система массовшх) обслуживания еледующмх) вида. На вход системы поступает простейший входящий ноток заявок, которые последовательно проходят две фазы. Каждая фаза содержит один обслуживающий прибор и один источник повторных вызовов (ИПВ), куда отправляется заявка, если застает прибор занятым. Приборы на обеих фазах обрабатывают заявку в течение экспоненциально раепределенших) промежутка времени. Находясь в ИПВ, заявка через экспоненциально распределенные промежутки времени делает попытки вновь встать на обслуживание. После обслуживания на второй фазе заявка считается обработанной и покидает систему. Целью данной работы является исследование числа заявок в ИПВ на обеих фазах. Исследование проводится методом асимнтотическохх) анализа в условии большой задержки в ИПВ. Получены семиинварианты нервохх) порядка и проведено сравнение с результатами имитационной модели.
Ключевые слова: система массовохх) обслуживания, RQ-система, источник повторных вызовов.
Введение. С середины XX века все большую роль стали играть телекоммуникационные системы: компьютерные и телефонные сети, системы передачи данных, радио, телевидение, мобильная связь и другие. Такие системы подходят иод определение систем массового обслуживания |1|: они предназначены дня удовлетворения массовых запросов на выполнение каких-либо услуг. Однако, в отличие от классических систем массового обслуживания, дня телекоммуникационных систем характерна ситуация, при которой заявка, заставшая обслуживающий прибор занятым, не встает в очередь, а уходит в источник повторных вызовов (ИПВ), откуда через некоторые промежутки времени предпринимает попытки вновь обратиться за обслуживанием. Такие модели описываются в виде систем массового обслуживания с повторными вызовами или RQ-систем (Retrial queuing system) |2,3|.
В данной работе исследуется RQ-система М/М/1 с двумя фазами, каждая из которых содержит один обслуживающий прибор и ИПВ. Впервые двухфазные RQ-системы были рассмотрены С. М. Krishna и Y. Н. Lee |4|, Также их анализом занимались В. Т. Doshi |5|, В. Krishnakumar |6|, G. Choudhury |7| и другие.
1. Постановка задачи. Рассмотрим RQ-систему, на вход которой поступает простейший входящий ноток интенсивности А. Каждая заявка проходит последовательно 2 фазы.
Рис. 1. Модель двухфазной RQ-систомы
Продолжительность обслуживания на обеих фазах имеет экспоненциальное распределение с параметрами ц\ и /12 соответственно. Если на очередной фазе заявка застает прибор занятым, она отправляется в соответствующий ИПВ, откуда через случайные моменты времени, имеющие экспоненциальное распределение с параметрами о\ и ег2, делает попытки вновь обратиться за обслуживанием. После завершения обслуживания на второй фазе заявка считается обработанной и покидает систему (рис. 1).
Обозначим через = 1,2} число заявок в ИПВ, соответствующее ^'-ой фазе, в
момент времени Процессы { = 1,2} не являются марковскими, поэтому введем
дополнительную переменную:
{kj(t)
прибор па j-й фазе свободен, прибор на j-й фазе занят.
Пусть Pkik2(ii,i2,t) = P{ki(t) = ki,k2(t) = k2,ii(t) = H,i2(t) = i2} — вероятность того, что в системе в момент времени t при условии, что приборы па обеих фазах находятся в состояниях к\ и к2 соответственно, в ИПВ па первой и второй фазе содержится i\ и г2 заявок. В данной статье мы будем рассматривать систему в стационарном режиме, при котором lim Pkik2(ii,i2,t) = Ркifc2 (¿1,г2).
г—Юо
Очевидно, что процесс {ij(t),j = 1,2, kj(t),j = 1,2} является марковским. Дня получения распределений вероятностей Рк\к2(н, ¿2) составим систему уравнений Колмогорова:
'-(A-Hi<7i +г202)-Роо(®1,»2) +At2-Poi(«l,«2) = О,
АРю(«1 - 1,«з) + ЛРоо(гьг2) + oi(ii + l)Poo(«i + М2) ~{\ + fJ,i + i2a2)P10(i1,i2)+ +ß2Pii(h,i2) = О,
ßiP0i(ii,i2) + /xiPii(»i,»2 - 1) + 0-2(^2 + 1)Роо(гъг2) - (А + ц2 + «i<Ti)Poi(»i,i2) = О, APn(«i - 1,г2) + APoi(н,г2) + <Ti(ii + l)P0i(«i + M2) + ^(¿г + l)Pio(®i,®2 + 1)--(Л + ßl + ß2)Pll{i!,i2) = 0.
(1)
В системе (1) перейдем к частичным характеристическим функциям:
Нк1к2(иъ и2) = У^У^ехр{j(uiii + и2г2)}Рк1к2(гг, г2), ¿1 ¿2
где 3 = л/— 1 — мнимая единица. Тогда система уравнений (1) для характеристических функций перепишется в виде:
\ тт ( ч 1 9Яоо(«1,«2) 1 дн00(иии2) , ,
-\Нт(иъи2) - -<71-д---:<72-■=--1"А»2-П01 («1,142) = О,
3 дщ 3 ди2
Хе^Н10(иъи2) + ХН00(и1,и2) + - (А + ^)Н10(иъи2)-
1 аяю(и1,и2) . гт ( Ч п
--Г02-я--К/Х2Яц(мьм2) = О,
3 ди2
ßiHw(uuu2) + /ле^ЯпКиа) + - (Л + ß2) Н01(щ,и2)-
3 ои2
1 0#oi(ui,U2) п
—:0"1-д-= О,
J дщ
Л jui тт ( \ , i гг / X . 1 _7-и10Я01(и1,и2) . 1 _7Чи0Яю(и1,И2)
(2)
3 дщ з ди2
"(А + 1^1+ /¿2)Яц («1,1x2) = 0.
Систему (2) будем решать методом асимптотического анализа [8] в условии большой задержки в ИПВ, характеризующегося тем, что (Т\,сг2 —> 0,
2. Асимптотика первого порядка, В соответствии с условием большой задержки в ИПВ, обозначим а& = 7^<7, а —> 0. В системе (2) выполним замены:
а = е, ик = етк, Нк1к2{иъи2) = Рк\к 2(уоъь)2£). (3)
Получим:
/
\ 771 / ч 1 dFoo(wx,w2,£) 1 dFoo(wx,w2,£)
—XF0O(wi,W2,£) - -71-----Г72---h /i2-Foi(wi,iü2,e) = о,
J CWi 3 OW2
Xe^F10(Wl,w2,e) + XFm{wuw2,e) + I^-i«*0^1'«*,*)
i öwi
/X , t \ 1 0Fio(«;i,iU2,e) , ^ / ч n
-(A + /xi)Fio(i(;i,W2,e) - -72-ä--ß2F-il(wi,w2)e) = 0,
3 ow 2
p f ^ , ( \ , 1 -7И02 öi?oo(wbw2,e) (4)
3 ow2
/л 1 * p / 4 1 9F0I(wi,W;2,£) n
- (A + /i2) Foi^i/u^e) - J7i-g^-= 0>
A e^FnK^e) + АЗДи^е) +
3 ow\
+-726 36W2-ä--(A + //1 + /¿2) Fn(it;i,ii;2,e) = 0.
3 ow2
Теорема 1, Если существуют пределы lim FkikzipiWi, s) — ^шЦ,^), то
¿ЪшОад) = (5)
где {Ик1к2,к1,к2 = 1,2} — стационарные вероятности значений цепи Маркова — 1,2}, определяемые уравнением (9), а величины а\ и а2 определяются из системы (13),
Доказательство. Пусть Рк\к2 (^1,^2) = Я/люФ (^1,^2)- Тогда после предельного перехода е —> 0 система (4) запишется в виде:
Л с ж/ \ 1 о дФ^!^) 1 п дФ^^) ,
-ХКооФ^г^) - -71Д00-я---:72-Коо-~-102+
3 ди)! з д
+ц2Я 01Ф(г01,ги2) = О,
АДюФ^ь^г) + ХВооФ^иьъ) + - (А + ц^КюФ^г^)-
+ = О,
3 дт2
^ ж. 1 <9Ф(г1>1,г<;2)
/Xliíl0®(tüi,tü2) + ¡JLlRl^{WUW2) + T72#00
(6)
3 dw2
- (А + ß2) ЯогФ(W1,w2) - =
ХППФМ + АД01Ф(ВД) + -ъПогЩ^ + ~ъЪ0дФМ
j dwi j dw2
— (A + /ii + ¡i2) RI^(wi,W2) — 0.
Чтобы избавиться ОТ производных ——и —{^liw2) положим Ф (Wl WA =
OW1 ow2
ej(w1a1+w2a2)^ 0ТКуда ——(^ь^г) ф (wi,w2) = jak. Тогда после сокращения системы (6)
" dwk / на Ф (wi,w2) получаем:
r — Ai?oo - 7i°iPoo - 72O2-R00 + V2R01 = 0,
A-Rio + АДоо + JiaiRoo - (А + (J,i)Rw - 72^10 + V2R11 = 0, ^
AiiДю + Ati-Rn + 72Й2-Коо - (А + ¡л2) Roi - 71^1-Roi = 0, k XR11 + XR01 + 7iaii?oi + ^2a2RiQ — (A + + fi2) Rn = 0.
С учетом условия нормировки Rqq + Rqi + Rio + Ru = 1, запишем систему (7) в матричном виде:
(
RM = 0, , .
RE = 1. (8)
Здесь R = ( Roo Roi Rio Ru ). E единичный вектор-столбец, a матрица M пред-ставима в виде:
М =
( - (Л + 71О1 + 72а2) А + 71 Й1 72а2 0 \
¿¿2 О -(А + ц2 + 7101) А + 71Й1
О -(/¿1+72^2) 72 а2
\ О /Х2 -(//1+Д2) У
Введя обозначения М = (М |.Е), г; = (Ет0 |1), перепишем систему (8) в виде:
ДМ = V,
откуда уравнение для стационарного распределения значений цепи Маркова {А^ = 1,2} имеет вид:
Д = {уМ)Т {ММТ)~1.
(9)
Здесь Д = /{(«] ,а2). Для нахождения «1 и а2 просуммируем уравнения системы (4) и после сокращений получим:
А (е-76™1 - 1) Гю(«л,ш2,е) + /л (е^2 - 1) ^и^.ша.е) + А (е^1 - 1)
+21 _ 1) д*Ьо(м1,М2,е) + 72 /
+72 _ ^ + 71 ^ _ ^ ^01(^1 №,е) = (10)
3 4 ' 3
Поделив уравнение (10) на е при е —> 0 и учитывая, что
дги 1
lim-(e±jswь-l)=±jwk,
е—»0 £ 4 '
получим уравнение:
дт 1
+72го2
&Р00(г<;1,г<;2)
+ 72гу2
а^0(г<;1,гу2)
+ 7^1
9-^01 (^1,г(;2)
ди>2 ди)2 ди> 1
Приравнивая коэффициенты при и ги2, получим систему:
А^-Г1о(гУ1,гу2) + А^Цги^гиг) - ъ-бт:--Ъ-~-= 0,
= 0. (11)
ди)1
■ V < \ дРоо^!^) дР10 (И)Ъ1Ю2) _ МгпКМъЩ) ~ 72-^--72-~-= 0.
дтг
(12)
_ дт2 д'ш2
Поделив систему (12) на 3 и используя представление (5) для функций ^ша^ъгиг); получим систему уравнений для нахождения а\ и а2:
АДю(а1,а2) + АДц(а1,а2) - 7101-^00(01,02) - 7101-^01(01,02) = 0, /¿1-йц(а1 ,а2) - 72^2-^00(01,02) - 7202-^10(01,02) = 0.
В силу замены (3) и равенства (5) можно записать:
Таблица 1
Сравнение асимптотических семиинвариант с результатами имитационного моделирования
а Асимптотические семиинварианты Результаты имитационного моделирования «1 _ --7771 а «2 _ --7772 а
а\ /а «2 /О" 7711 7712
0,100 0,010 0,001 5 50 500 0,828 8,284 82,823 5,479 49,728 481,340 1,190 8,175 79,010 0,479 0,272 18,660 0,362 0,109 3,813
Нк1к2{иъи2) = Ршя(г^г^е) » Як1к2ехр ¡^(щаг + и2а2)
Отсюда безусловная характеристическая функция
Мез{и 1Ч+и2г2)
процессов
= 1,2} имеет вид:
Мез{и,п+и2ъ) = и ехр и(и1а1 + и2а2)
к1 к2 ^
следовательно, двумерный процесс {ij(t),j = 1,2} числа заявок в ИПВ обеих фаз имеет асимптотически нормальное распределение. Величины ^ и которые определяют асимптотические средние значения числа заявок в источниках повторных вызовов, будем называть асимптотическими семиинвариантами первого порядка,
3. Результаты имитационного моделирования, В данном разделе проводится сравнение полученных асимптотических семиинвариант первого порядка с результатами имитационного моделирования для математических ожиданий числа заявок в ИПВ при
а —У 0, Оценки математических ожиданий определяются по формуле гЩ = — гДе
пк ~ общее число заявок в ИПВ на к-й фазе, Т.— период времени, в течение которого в ИПВ на к-й фазе находилось г заявок, и Т полное время моделирования. Моделирование проводилось при числе заявок во входящем потоке, равном 100000,
Покажем отклонение значений семиинвариантов первого порядка, полученных с помощью метода асимптотического анализа от эмпирических оценок математических ожиданий числа заявок в ИПВ при заданных значениях параметров: Л = 0,5, ¡1\ = [12 = 1, 71 = 1, 72 = 5 (табл. 1),
Из таблицы 1 видно, что при уменьшении а повышается точность аппроксимации, но при о = 0,001 результаты моделирования оказались хуже, чем при больших значениях а. Это связано с тем, что системе не хватило 100000 заявок во входящем потоке для сбора всей необходимой информации, поэтому число заявок в этом случае следует увеличить.
Заключение, В настоящей работе исследована двухфазная КС^-система М|М|1 и найдены моменты первого порядка для числа заявок в ИПВ, Результаты имитационного моделирования подтверждают теоретические значения числовых характеристик.
■
Список литературы
1. Гнеденко Б. В. Введение в теорию массового обслуживания. 2-е изд., перераб. и доп. // Б. В. Гнеденко, И.Н. Коваленко. М.: Наука. Гл. ред. физ.-мат. лит., 1987.
2. Artalejo J.R., Gomez-Corral A. Retrial Queueing Systems: A Computational Approach. Berlin: Springer, 2008.
3. Falin G.I., Templeton J.G.C. Retrial queues. London: Chapman & Hall, 1997.
4. Krishna С. M. and Lee Y. H. A study of two-phase service // Operations Research Letters. 1990. Vol. 9. P. 91-97.
5. Doshi B.T. Analysis of a two phase queueing system with general service times // Operations Research Letters. 1991. Vol. 10. P. 265-272.
6. B. Krishnakumar and D. Arivudainambi. An M/G/l retrial queueing system with two phases of services and preemptive resume // Annals of Operations Research. 2002. Vol. 113. P. 61-79.
7. Choudhury G. A single server queueing system with two-phases of service and vacations // QTQM. 2008. Vol. о (1). P. 33-49.
8. Nazarov A. A., Sudyko E. A. Method of asymptotic semiinvariants for studying a mathematical model of a random access network // Problems of information transmission. 2010. Vol. 46, N 1. P. 86102.
используемых при исследованиях математических моделей телекоммуникационных систем в теории телетрафика и массового обслуживания. Опубликовал (в соавторстве) 6 монографий и 2 методических пособия. Им опубликовано более четырехсот научных работ, из которых более 150 входят в РИНЦ и более 70 — в Scopus.
Anatoly Nazarov received M. D. degree in Mathematics at Tomsk State University in 1970, Ph.D. degree in physical-mathematical sciences in 1976. His thesis was on optimal control of queuing systems. He worked at Kemerovo State University from 1976 to 1980. He was Associate Professor of Probability Theory and Mathematical Statistics Department at the Faculty of Applied Mathematics and Cybernetics at Tomsk State University from 1980 to 1985. Since 1985, he has been a Professor. He received Doctor of Science degree in 1985. He has been the Head of the Department of Probability theory and Mathematical Statistics since 1999. He is a well-known specialist in teletraffic theory and queuing theory. He is also one of the developers of the method of asymptotic analysis in queuing theory, the dynamic screening method and the limit decomposition method, which are widely used in studying mathematical models of telecommunication systems in teletraffic theory and queuing theory. He published 6 monographies jointly with other professors. He published over
Nazarov А. A., Dr. Sc.
(Engineering), Professor, Head of the Department, e-mail:
Анатолий Назаров в 1970 г. окончил механико-математический факультет Томского государственного университета по специальности „Математика". В 1976 г. защитил диссертацию на соискание ученой степени кандидата физико-математических наук, посвященную вопросам оптимального управления в системах массового обслуживания. С 1976 г. по 1980 г. работал в Кемеровском государственном университете. С 1980 г. — доцент кафедры теории вероятностей и математической статистики факультета прикладной математики и кибернетики Томского государственного университета. В 1985 г. защитил диссертацию на соискание ученой степени доктора технических наук. С 1988 г. — профессор кафедры теории вероятностей и математической статистики. С 1999 г. — заведующий кафедрой теории вероятностей и математической статистики. Является известным специалистом в области теории телетрафика и массового обслуживания, а также одним из основных разработчиков методов асимптотического анализа в теории массового обслуживания, метода просеянного потока и метода предельной декомпозиции, широко
400 papers; about 150 papers are included in Russian Science Citation Index and about 70 papers — in Scopus.
Anisimova А. А., В A
(Economics), Student, e-mail: [email protected].
Анна Анисимова в 2016 году окончила бакалавриат факультета прикладной математики и кибернетики Томского государственного университета по специально-
сти „Экономика". В данный момент обучается в магистратуре по направлению „Прикладная математика и информатика". Ею опубликованы статьи на тему имитационного моделирования и теории массового обслуживания.
Anna Anisimova received her Bachelor degree in Economics at Tomsk State University in 2016. Now she is a Master's Degree student at the Faculty of Applied Mathematics and Cybernetics. She published papers in the areas of simulation modeling and queuing theory.
A
Дата поступления — 14-01.2017
