CC BY
16
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2016 Управление, вычислительная техника и информатика № 4 (37)
УДК 519.872
DOI: 10.17223/19988605/37/6
А.А. Назаров, Е.А. Фёдорова
МОДИФИКАЦИЯ МЕТОДА АСИМПТОТИЧЕСКОГО АНАЛИЗА В УСЛОВИИ БОЛЬШОЙ ЗАГРУЗКИ НА ПРИМЕРЕ RQ-СИСТЕМЫ M|M|1
Работа выполнена в рамках проекта Министерства образования и науки РФ №1.511.2014/К.
В работе предлагаются три метода модификации асимптотического анализа в условии большой загрузки RQ-систем. Описание методов проведено на примере RQ-системы М|М|1. Предложены «промежуточная» асимптотика, имеющая вид гипергамма распределения, модификация «промежуточной» асимптотики в виде сдвига аргумента на параметр 5 и полусумма распределений модифицированной и второй асимптотики. Представлены результаты численного анализа предложенных модификаций.
Ключевые слова: RQ-система; метод асимптотического анализа; большая загрузка; гипергамма-распределение.
Retrial Queueing System (системы массового обслуживания с повторными вызовами, или RQ-системы) - математические модели, широко применяемые для анализа и оптимизации различных телекоммуникационных систем, сетей мобильной связи, call-центров и других технических и экономических систем [1-6]. Характерной чертой таких моделей является наличие повторных обращений заявок к прибору после неудачной попытки обслуживания спустя некоторое случайное время. Такие ситуации могут быть вызваны не только отсутствием свободных серверов в моменты поступления заявок в систему, но техническими причинами.
Возникновение моделей RQ-систем прежде всего связывают с работами американских ученых R.I. Wilkinson [1] и J.W. Cohen [2] середины ХХ в., которые были посвящены практическим задачам, возникающим в телефонных сетях, и описанию влияния эффекта повторных вызовов на производительность технических систем. Первые подходы к математическому описанию RQ-систем были предприняты G. Gosztony [3] и A. Elldin [4].
Наиболее полное и детальное описание RQ-систем и их сравнение с классическими СМО было отражено в монографиях J.R. Artalejo, A. Gomez-Corral, G.I. Falin и J.G.C. Templeton [7, 8]. Ими получены аналитические результаты для RQ-систем М|М|1, M|GI|1, М|М|С и других систем с пуассоновским входящим потоком, а также рассмотрены разнообразные методы для исследования таких систем и проведено детальное сравнение RQ-систем с их классическими аналогами [9].
На сегодняшний день исследованию RQ-систем посвящено большое количество работ [7], однако аналитических формул получено немного (лишь для систем с простейшим входящим потоком) [8]. В основном возникающие задачи решаются численными методами или с помощью имитационного моделирования [10-12], и как следствие, результаты таких исследований имеют узкое применение.
Ранее в работах [13, 14] был предложен метод асимптотического анализа (первого и второго порядков) для исследования RQ-систем при условии большой загрузки. Предложенным методом были исследованы системы с различными типами входящего потока и дисциплиной обслуживания (в том числе RQ-система ММРР^|1 [15]). Было показано, что асимптотика первого порядка имеет вид гамма-распределения вне зависимости от типа системы, но область ее применения достаточно узка. Асимптотика второго порядка позволила расширить область применения результатов в 4 раза, однако асимптотическая формула не обладает единообразием и требуется ее достаточно трудоемкий вывод для каждой системы.
В связи с этим в данной статье предлагаются способы модификации метода асимптотического анализа в условии большой загрузки RQ-систем, позволяющие увеличить точность асимптотического метода. Описание методов приведено на примере простейшей RQ-системы М|М|1 для возможности
сравнения асимптотических результатов с известным точным. В дальнейшем предложенные модификации могут быть применены для исследования различных Я^-систем, в том числе с непуассоновским входящим потоком.
1. Математическое описание ЯО-системы М|М|1
Рассмотрим однолинейную Я^-систему (рис. 1), на вход которой поступает простейший поток заявок с параметром X, время обслуживания каждой заявки распределено по экспоненциальному закону с параметром ц. Если поступившая заявка застает прибор свободным, то она занимает его для обслуживания. Если прибор занят, то заявка переходит в источник повторных вызовов (ИПВ), где осуществляет случайную задержку, продолжительность которой имеет экспоненциальное распределение с параметром о. Из ИПВ после случайной задержки заявка вновь обращается к обслуживающему прибору с повторной попыткой его захвата. Если прибор свободен, то заявка из ИПВ занимает его для обслуживания, в противном случае она мгновенно возвращается в источник повторных вызовов для реализации следующей задержки.
* (' ) =
Рис. 1. Я^-система М|М|1
Пусть /(0 - число заявок в ИПВ, а к(?) определяет состояние прибора следующим образом:
[0, если прибор свободен, [1, если прибор занят.
Обозначим Р{к(0 = к, /(0 = /} = Р(к,/^) вероятность того, что прибор в момент времени ^ находится в состоянии к и в источнике повторных вызовов находится / заявок.
Причем процесс {к(0, /(0} изменения состояний данной системы во времени является марковским. Тогда ставится задача нахождения распределения вероятностей числа заявок в источнике повторных вызовов такой системы.
Для распределения вероятностей Р(к,/^) состояний рассматриваемой Я^-системы составим систему дифференциальных уравнений Колмогорова:
dP(0, t)
dt
dP(1, i, t)
dt
= i,t) - (к + ia)P(0,i,t),
= kP(0,i, t) + (i + 1)a • P(0,i + 1, t) - (к + ц)P(1,i,t) + kP(1,i -1,t).
(1)
Перейдем к стационарным частичным характеристическим функциям H(к, u) = £eJulP(k, i), где
i
P(k, i) = lim P(k , i,t), а J - мнимая единица.
Введем параметр p = —, характеризующий загрузку системы. Тогда система (1) в стационарном
ц
режиме для характеристических функций перепишется в следующем виде:
ц du
H (1, u )-pH (0, u ) = -a J ^^
pH(0,u)+(p( -1)-1)(1,u) = jüe-JudH(u.
ц du
2. Метод асимптотического анализа в условии большой загрузки
В работах [13-15] для решения системы уравнений (2) был предложен метод асимптотического анализа в условиях большой загрузки, т.е. при р Т 1, который заключается в следующем.
Вводятся обозначения в = 1 - р (тогда в условии большой загрузки в I 0), и = вм, Н(0,и) = в^0(м,в), Н(1,и) = ^(м,в) .
Таким образом, система (2) принимает вид
в) - (1 -в)в^ (м, в) + j = 0,
ц См
(1 - в)в^(м, в) + ((1 - в)(ем -1)- ), в) - j 1 е-м ^^^ = 0.
Ц СМ
(3)
Далее используются следующие разложения:
^ (м, в) = ^ (м) + в- Л (М + О(в 2), (4)
где О (в ) - бесконечно малая величина порядка в2 .
Допредельная характеристическая функция Н (и) = Н (1, и) + Н (0, и) в условиях большой загрузки может быть приближенно определена равенством
Н (и) « ^
( и >
1 -Р
+ О (в 2).
Обозначим к1(и) асимптотическую характеристическую функцию первого порядка, которая равна
Ми) = ^
( и ^
1 -Р
В [13] было доказано, что ЛДи) удовлетворяет следующей теореме.
Теорема 1. Асимптотическая характеристическая функция числа заявок в ИПВ в Я^-системе М|М|1 имеет вид характеристической функции гамма-распределения
/ \"а и
К (и) = 1 - ]—
(5)
с параметрами формы а = —+1 и масштаба р = 1 - р .
а
В качестве критерия допустимости асимптотических результатов использовалось условие: расстояние Колмогорова Д < 0,05 .
С помощью численного анализа было показано, что метод асимптотического анализа первого порядка в условиях большой загрузки может быть применим при значениях загрузки р > 0,95, что является достаточно узкой областью применимости. Поэтому также была получена асимптотическая характеристическая функция второго порядка К2(и):
К2(и) = ^ I — | + (1 -р)
^т-РГ л I 1-р
Для получения этой асимптотики второго порядка вместо разложений функций (4) необходимо было рассмотреть разложения
^(м,в) = ^(м) + в • /к(м) + в2 • (м) + О(в3) и, подставив их в систему (3), решить систему пяти дифференциальных уравнений.
Тогда допредельная характеристическая функция может быть определяется равенством
Н (и) = К2(и) + О(в3). Была доказана следующая теорема о виде функции К2(и) [14].
Теорема 2. Асимптотическая характеристическая функция второго порядка распределения вероятностей числа заявок в ИПВ в Я^-системе М|М|1 имеет вид
И2(и) = 1 1" 1
1 "Р
ц+а а
1+(1 "Р)
1 1П1 - 1
1 "Р а
1"Р
]и
"|Ц+111 "Р"
1 "Р
а
1"
_]и_
1 "Р
Численный анализ показал, что область применения асимптотики второго порядка р > 0,8, т.е. увеличивается в 4 раза (по сравнению с (5)). Однако для более сложных Я^-систем, например ММРР|01|1, приходилось решать систему из семи и более дифференциальных уравнений, что достаточно трудоемко.
3. Гипергамма-асимптотика
Заметив, что в ходе доказательства теоремы 1 была найдена и функция -^о(^), было предложено провести численный анализ «промежуточной» асимптотики И (и):
и (и) = ^
( и 1
1 "Р.
+ (1 "РЖ,
( и 1
1"Р
В работе [13] показано выполнение следующих равенств:
ц
Ро(м) = С(1 " ]м>)~а ,
<
ад = с (1"
Имеем условие нормировки И(0) = Н0(0) + Н1(0) = 1. Так как Н1(и) « Р1(и), а Н0(и) « (1"р)(и),
1
то получаем И(0) = (1" р) (0) + ^(0) = (1" р)С + С = 1. Отсюда С =
2 "р
Тогда имеем
И*(и) = 11"
2"Р1 1 "Р
—1 1
а 1"Р Ь и
+--I 1" ]-
2 " Р I 1" Р
Заметим, что слагаемые функции И (и) являются характеристическими функциями гамма-распределения:
ц(и) = |1"^ , Г0(и) = I 1 , где а = —+1 и р = 1 "Р .
V
р
]и
X
р
ц
а
Таким образом, получили, что «промежуточная» асимптотическая характеристическая функция распределения числа заявок в ИПВ имеет вид взвешенной суммы двух характеристических функций гамма-распределения
И*(и) = дГ1 (и) + (1" д)Г0 (и), (6)
где д =
1
2"Р
Будем называть такое распределение гипергамма-распределением (по аналогии с гиперэкспоненциальным).
В табл. 1 приведены Д1, Ар, Д2 - расстояния Колмогорова между точным распределением и первой асимптотикой, гипергамма-распределением и второй асимптотикой соответственно.
Из табл. 1 видно, что гипергамма-асимптотика значительно точнее асимптотики первого порядка, но уступает асимптотике второго порядка.
2
2
и
а
Т а б л и ц а 1
Расстояние Колмогорова между точным и асимптотическими распределениями
Значения параметров системы Р = 0,7 Р = 0,8 Р = 0,9 Р = 0,95
а = 10 Д1 0,251 0,170 0,086 0,043
АР 0,097 0,060 0,028 0,014
Д2 0,035 0,022 0,009 0,005
а = 2 Д1 0,266 0,179 0,091 0,046
Ар 0,157 0,104 0,050 0,025
Д2 0,088 0,047 0,016 0,005
а = 1 А1 0,304 0,208 0,106 0,053
АР 0,220 0,147 0,074 0,037
А2 0,124 0,052 0,013 0,003
Однако вывод асимптотической формулы второго порядка для более сложных Я^-систем затруднителен, тогда как «промежуточная» асимптотика не требует дополнительных выкладок, так как ее форма одинакова для Я^-систем с различными типами входящего потока и дисциплины обслуживания, а параметры гипергамма-распределения совпадают с параметрами первой асимптотики. Поэтому целесообразность ее использования для большинства практических задач очевидна.
4. Модификация гипергамма-асимптотики
Известно, что Н1(0) = В1 = р, Н0(0) = В0 = 1 - р . Однако в полученной «промежуточной» асимптотике (6)
д^(0) = д * Я1 и (1 - д)Г>(0) = 1 - д * В>.
Перейдем от характеристических функций к функциям распределения вероятностей, обозначив их соответственно Г0 (х) и Г1 (х). Тогда дГ1 (да) = д * В1 и (1 - д)Г0 (да) = 1 - д * В0 .
В связи с этим предлагается модифицировать асимптотику следующим образом:
>1( х) = д [Г 1( х + 5) -Г 1(5) ],
х) = ж + (1 -д) [Г 0( х + 5)-Г 0(5) ],
где х) и х) - модифицированные частичные асимптотические функции распределения, а параметры 5 и п могут быть найдены из условия нормировки:
1>1(да) = д[1 -Г1(5) 1^0 (да) = ж + (1 - д)
Из первого уравнения (7) имеем
= р~ , (7)
1 - Гэ(5)]= 1 -р.
г1(5)=1 -р=д-р. д д
Подставляя введенные ранее обозначения, нетрудно получить
ГД5) = (1 -р)2. (8)
Выразим величину п из второго уравнения (11):
(1 -р)2 + (1 - р)Г0 (5)
ж = •
2-р
Подставляя (8), получим следующее выражение:
ж=Г,(5) + (1-р)Г0(5) _д (5) + (1 _д)Г0(5). (9)
2-р
Модифицированная асимптотическая функция распределения записывается как сумма
^ (х) = х) + х). Учитывая полученные выражения (11)—(13), имеем следующее равенство:
^ (х) = дГ1 (х + 5) + (1 - д)Г 0 (х + 5). Проведем численный анализ модификации асимптотики и сравним ее с результатами, полученными ранее асимптотики (табл. 2, где ДМ - расстояние Колмогорова между точным распределением и модифицированной асимптотикой).
Т а б л и ц а 2
Расстояние Колмогорова между точным, асимптотическими и модифицированным распределениями
Значения параметров системы Р = 0,5 Р = 0,6 Р = 0,7 Р = 0,8 Р = 0,9 Р = 0,95
о = 10 Д1 0,378 0,321 0,251 0,170 0,086 0,043
Др 0,180 0,137 0,097 0,060 0,028 0,014
ДМ 0,065 0,059 0,050 0,038 0,022 0,012
Д2 0,042 0,041 0,035 0,022 0,009 0,005
о = 2 Д1 0,432 0,342 0,266 0,179 0,091 0,046
Др 0,271 0,213 0,157 0,104 0,050 0,025
ДМ 0,079 0,077 0,068 0,050 0,031 0,017
Д2 0,192 0,139 0,088 0,047 0,016 0,005
о = 1 Д1 0,487 0,394 0,304 0,208 0,106 0,053
Др 0,364 0,291 0,220 0,147 0,074 0,037
ДМ 0,093 0,087 0,074 0,052 0,030 0,016
Д2 0,313 0,214 0,124 0,052 0,013 0,003
Данные табл. 2 показывают, что «модифицированное» распределение существенно улучшает результаты «промежуточной» асимптотики. При этом для построения модифицированного распределения необходимо лишь найти величину сдвига 5 из условия (8).
Рассмотрим комбинацию «модифицированного» распределения и второй асимптотики
Рт(2) = 2 (Рш + Р(2)).
Будем называть такое распределение «модификация-2». Проведем численный анализ модификации асимптотики и сравним ее с результатами, полученными ранее.
Т а б л и ц а 3
Расстояние Колмогорова между точным, асимптотическими и модифицированными распределениями
Значения параметров системы Р = 0,5 Р = 0,6 Р = 0,7 Р = 0,8 Р = 0,9 Р = 0,95
о = 10 Д1 0,378 0,321 0,251 0,170 0,086 0,043
Др 0,180 0,137 0,097 0,060 0,028 0,014
ДМ 0,065 0,059 0,050 0,038 0,022 0,012
Д2 0,042 0,041 0,035 0,022 0,009 0,005
ДМ2 0,034 0,024 0,018 0,013 0,009 0,005
о = 2 Д1 0,432 0,342 0,266 0,179 0,091 0,046
Др 0,271 0,213 0,157 0,104 0,051 0,025
ДМ 0,079 0,077 0,068 0,050 0,031 0,017
Д2 0,192 0,139 0,088 0,047 0,016 0,005
ДМ2 0,111 0,073 0,039 0,019 0,013 0,008
о = 1 Д1 0,487 0,394 0,304 0,208 0,106 0,053
Др 0,364 0,291 0,220 0,147 0,074 0,037
ДМ 0,093 0,087 0,074 0,052 0,030 0,016
Д2 0,313 0,214 0,124 0,052 0,013 0,003
ДМ2 0,150 0,087 0,046 0,022 0,012 0,007
В результате численного анализа был замечен интересный результат: «модификация-2» точнее второй асимптотики. При этом предложенная модификация применима для значений загрузки р>0,7, а при больших значениях параметра задержки о - для значений загрузки р>0,5.
Заключение
Таким образом, в работе предложены три модификации метода асимптотического анализа в условии большой загрузки Я^-систем: «промежуточная» асимптотика, имеющая вид гипер гамма-рас -
пределения; модификация «промежуточной» асимптотики в виде сдвига распределения на параметр 5; полусумма распределений модифицированной и второй асимптотики. В статье также представлены результаты численного сравнения предложенных асимптотик с известным точным, позволяющие сделать заключения об области применения каждого метода.
Обобщая вышесказанное, можно сделать вывод, что для практических задач, не требующих высокой точности аппроксимации, целесообразно использование результатов «промежуточной» асимптотики в силу ее простоты построения (имеет единую форму вне зависимости от типа входящего потока) и достаточно широкой области применимости (р > 0,8), а для задач, требующих более высокой точности аппроксимации, целесообразно применять предлагаемые в работе модификации 1 и 2.
ЛИТЕРАТУРА
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. Cohen J.W. Basic problems of telephone trafic and the influence of repeated calls // Philips Telecommunication Review. 1957. V. 18,
No. 2. P. 49-100.
3. Gosztony G. Repeated call attempts and their efect on trafic engineering // Budavox Telecommunication Review. 1976. V. 2. P. 16-26.
4. Elldin A., Lind G. Elementary Telephone Trafic Theory. Ericsson Public Telecommunications, 1971.
5. Kuznetsov D.Yu., Nazarov A.A. Analysis of non-Markovian models of communication networks with adaptive protocols of multiple
random access // Avtomatika i Telemekhanika. 2001. V. 5. P. 124-146.
6. Nazarov A.A., Tsoj S.A. Common approach to studies of Markov models for data transmission networks controlled by the static ran-
dom multiple access protocols // Avtomatika i Vychislitel'naya Tekhnika. 2004. V. 4. P. 73-85.
7. Artalejo J.R., Gomez-Corral A. Retrial Queueing Systems. A Computational Approach. Springer, 2008.
8. Falin G.I., Templeton J.G.C. Retrial queues. London : Chapman & Hall, 1997.
9. Artalejo J.R., Falin G.I. Standard and retrial queueing systems: A comparative analysis // Revista Matematica Complutense. 2002.
V. 15. P. 101-129.
10. Neuts M.F., Rao B.M. Numerical investigation of a multiserver retrial model // Queueing Systems. 2002. V. 7, No. 2. P. 169-189.
11. Ridder A. Fast simulation of retrial queues // Third Workshop on Rare Event Simulation and Related Combinatorial Optimization Problems. Pisa, 2000. P. 1-5.
12. Kim C.S., Mushko V.V., Dudin A. Computation of the steady state distribution for multi-server retrial queues with phase type service process // Annals of Operations Research. 2012. V. 201, No. 1. P. 307-323.
13. Moiseeva E., Nazarov A. Asymptotic Analysis of RQ-systems M/M/1 on Heavy Load Condition // Proceedings of the IV International Conference Problems of Cybernetics and Informatics. Baku, Azerbaijan, 2012. P. 164-166.
14. Назаров А.А., Фёдорова Е.А. Метод асимптотического анализа в условии большой загрузки 2-го порядка на примере исследования RQ-системы M|M|1 // Информационные технологии и математическое моделирование (ИТММ-2013) : материалы XII всерос. науч.-практ. конф. с междунар. участием им. А.Ф. Терпугова : в 2 ч. Томск : Изд-во Том. ун-та, 2013. Ч. 2. С. 59-65.
15. Моисеева Е.А., Назаров А. А. Исследование RQ-системы MMPP|GI|1 методом асимптотического анализа // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 4 (25). С. 84-94.
Назаров Анатолий Андреевич, д-р техн. наук, профессор. E-mail: [email protected] Фёдорова Екатерина Александровна, канд. физ.-мат. наук. E-mail: [email protected] Томский государственный университет
Поступила в редакцию 14 августа 2016 г.
Nazarov Anatoly A., Fedorova Ekaterina A. (Tomsk State University, Russian Federation).
Modification of the asymptotic analysis method under heavy load condition on the example of the study of the retrial queueing system M|M|1.
Key words: retrial queueing system; asymptotic analysis method; heavy load; hypergamma distribution. DOI: 10.17223/19988605/37/6
In the paper, three methods of the asymptotic analysis modification are offered for retrial queueing systems under heavy load condition as mathematical models of telecommunication systems. The characteristic feature of these models is existence of repeated attempts of calls to get service after a random time. There are a large number of papers devoted to retrial queues researching, but mainly problems are solved by numerical methods and simulation. Analytical results are obtained only for systems with Poisson input process.
Earlier we proposed the method of asymptotic analysis (first and second order) for retrial queueing systems researching under heavy load condition [13, 14]. It has been shown that the first order asymptotics has the form of gamma distributions for all types of the system, but its application range is quite narrow. The asymptotics of second order has allowed to extend the application range of the results by 4 times, but the asymptotic formula has not unified form and it needs time-taking rather laborious formula derivations for each system. Therefore, we suggest ways to modify the method of the asymptotic analysis under heavy load condition to increase the method accuracy.
The method is described on the example of the simplest retrial queue M|M|1 in order to compare the results with the known exact distribution. The system of differential Kolmogorov equations for Markov process {k(t), i(t)} is composed, then it is rewritten for characteristic functions in steady state.
X
The parameter p = — is the load rate of the system. The principle of asymptotic analysis method under heavy load condition (where
p 11) is briefly described. Asymptotic characteristic functions of the first and second order are presented.
Note that during the derivation of the first order asymptotic formula the function F0(w) was obtained. We propose to study "intermediate" asymptotics h (u), which has the form of a weighted sum of two gamma distribution characteristic functions:
h (u) = qri(u) + (1 - q)ro(u) ,
i i \-a+1 i
where r1 (u) = i 1 - , r0 (u) = il - , q = — , « = ^ +1 u P = 1 -p .
We suggest to call this distribution as hypergamma distribution (like a hyperexponential distribution).
The numerical analysis shows that hypergamma asymptotics is more accurately than the first order asymptotics, but it is worth than the second order asymptotics. However, the derivation of the second order asymptotic formula for more complex retrial queues is difficult, whereas "intermediate" asymptotics does not require additional calculations, because its shape is the same for all retrial queues and hypergamma distribution parameters are obtained in the first order asymptotic formula theorem.
Then a modification of hypergamma asymptotics of shifting the argument to the parameter 5 is proposed:
F(x) = qf 1 (x + 8) + (1 - q)f 0(x + 8),
where 5 is found from the condition f 1(8) = (1 - p)2 .
The numerical analysis shows that "modified" distribution significantly improves results of hypergamma asymptotics.
In the paper conclusion a combination of "modified" distribution and the second asymptotics Pm(2) = 2(Pm + P(2)). is considered.
As a result of numerical analysis, it was observed that this modification is more accurately than the second order asymptotics and can be applied for load rate p>0.7 and even p>0.5 for large values of delay parameter a.
In summary, it can be concluded that for practical tasks which does not require high accuracy of approximation, it is useful to apply the "intermediate" asymptotics results because of its simplicity of construction and rather broad area of applicability (p> 0.8), and for applications requiring higher accuracy of approximation, it is advisable to apply the proposed modifications 1 and 2. In the future, the modifications can be applied to the study of various retrial queueing systems including with non Poisson incoming process.
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. Cohen, J.W. (1957) Basic problems of telephone traffic and the influence of repeated calls. Philips Telecommunication Review. 18(2).
pp. 49-100.
3. Gosztony, G. (1976) Repeated call attempts and their efect on trafic engineering. Budavox Telecommunication Review. 2. pp. 16-26.
4. Elldin, A. & Lind, G. (1971) Elementary Telephone Traffic Theory. Ericsson Public Telecommunications.
5. Kuznetsov, D.Yu., Nazarov, A.A. (2001) Analysis of non-Markovian models of communication networks with adaptive protocols of
multiple random access. Automation and Remote Control. 62(5). pp. 124-146. DOI: 10.1023/A:1010231008994.
6. Nazarov, A.A. & Tsoj, S.A. (2004) Common approach to studies of Markov models for data transmission networks controlled by the
static random multiple access protocols. Avtomatika i Vychislitel'naya Tekhnika - Automatic Control And Computer Sciences. 4. pp. 73-85. (In Russian).
7. Artalejo, J.R. & Gomez-Corral, A. (2008) Retrial Queueing Systems. A Computational Approach. Springer.
8. Falin, G.I. & Templeton, J.G.C. (1997) Retrial queues. London: Chapman & Hall.
9. Artalejo, J.R. & Falin, G.I. (2002) Standard and retrial queueing systems: A comparative analysis. Revista Matematica Complutense.
15. pp. 101-129. DOI: 10.5209/rev_REMA.2002.v15.n1.16950
10. Neuts, M.F. & Rao, B.M. (2002) Numerical investigation of a multiserver retrial model. Queueing Systems. 7(2). pp. 169-189. DOI: 10.1007/BF01158473
11. Ridder, A. (2000) Fast simulation of retrial queues. Third Workshop on Rare Event Simulation and Related Combinatorial Optimization Problems. Pisa. pp. 1-5.
12. Kim, C.S., Mushko, V.V. & Dudin, A. (2012) Computation of the steady state distribution for multi-server retrial queues with phase type service process. Annals of Operations Research. 201(1). pp. 307-323. DOI: 10.1007/s10479-012-1254-7
13. Moiseeva, E. & Nazarov, A. (2012) Asymptotic Analysis of RQ-systems M/M/1 on Heavy Load Condition. Problems of Cybernetics and Informatics. Proc. of the Fourth International Conference. Baku, Azerbaijan. pp. 164-166.
14. Nazarov, A. & Fedorova, E. (2013) [The second order asymptotic analysis method under heavy load condition on example of the retrial queue M|M|1 research]. Informatsionnye tekhnologii i matematicheskoe modelirovanie (ITMM-2013) [Information technology and mathematical modeling (ITMM 2013)]. Proc. of the 12th All-Russian Conference with International Participation. Vol. 2. Tomsk. pp. 59-65. (In Russian).
15. Moiseeva, E. & Nazarov, A. (2013) Researching of Retrial Queueing system MMPP|GI|1 by using asymptotic analysis method on heavy load condition. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(25). pp. 84-94. (In Russian).
