CC BY
44
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2014 Управление, вычислительная техника и информатика № 1 (26)
УДК 519.872
Л.А. Жидкова, С.П. Моисеева
ИССЛЕДОВАНИЕ ЧИСЛА ЗАНЯТЫХ ПРИБОРОВ В СИСТЕМЕ MMPP|M|^
C ПОВТОРНЫМИ ОБРАЩЕНИЯМИ
Работа выполнена в рамках государственного задания Минобрнауки РФ на проведение научных исследований в Томском государственном университете на 2012-2013 гг. (№ 8.4055.2011).
Рассматривается система массового обслуживания MMPP|Ma> с повторными обращениями в систему.
Найдены аналитические выражения для нахождения первого и второго моментов, характеризующих
число занятых приборов в системе, а также их асимптотическое приближение.
Ключевые слова: система массового обслуживания с повторными обращениями; марковский модулированный поток; метод асимптотического анализа.
Одним из перспективных направлений современной экономической науки является привлечение для описания и количественного выражения экономических реалий методов их математического описания и моделирования. В связи с этим возникает задача формализации той высокой неопределенности, которая свойственна реальным экономическим процессам.
В качестве математических моделей социально-экономических и демографических процессов часто используют системы массового обслуживания (СМО) с неограниченным числом обслуживающих приборов [1-3]. Как правило, в таких системах число потенциальных клиентов (страховых и торговых компаний, пенсионных фондов) считается неограниченным. В реальных технических системах число обслуживающих приборов конечно, но в случае, когда вероятностью потери заявки можно пренебречь, такие системы можно аппроксимировать бесконечнолинейными СМО. Исследование таких систем с пуассоновским входящим и произвольным временем обслуживания можно встретить в работах В.В. Рыкова, П.П. Бочарова, А.В. Печинкина и других авторов [4-7].
Однако применение пуассоновского потока для расчета характеристик качества обслуживания в реальных системах дает большую погрешность. Обоснованность адекватности применения марковского модулированного пуассоновского потока для описания информационных потоков в мультисер-висных сетях связи и телекоммуникационных системах следует из работ W.E. Leland, M.S. Taqqu, W. Willinger, V. Paxson, C. Lindemann, M. Lohmann и др. [8-9]. Системам с непуассоновскими входящими потоками (марковскими модулированными, полумарковскими) посвящены статьи D. Baum и L. Breuer [10, 11], M. Parulekar и A.M. Makowski [6], A.K. Jayawardene и O. Kella [12] и многих других.
Основными методами исследования СМО с неограниченным числом приборов, как правило, являются метод вложенных цепей Маркова и метод дополнительной переменной, в последние двадцать лет стали развиваться матрично-аналитические методы [13]. В тех случаях, когда не удается найти характеристики системы в явном виде, применяют асимптотические методы [14-16].
Одной из модификаций СМО с неограниченным числом приборов являются системы массового обслуживания с повторными обращениями, которые применяются для описания математических моделей, например, страховых или торговых компаний. Исследования таких систем с пуассоновским входящим потоком приведены в работах [2, 17], где методом производящих функций определены такие процессы, как число занятых приборов в системе и характеристики потоков, обращающихся в систему. Такие модели используют для описания процесса зависимости величины капитала коммерческой организации от стимулирующих программ, например предоставление скидки при приобретении товара; кроме того, подобные системы предлагаются в качестве математических моделей распределительных вычислительных сетей [1, 18].
Для аналогичных систем с произвольным временем обслуживания в работе [19] предложен метод предельной декомпозиции, позволяющий свести исследование бесконечнолинейной системы массового обслуживания к исследованию совокупности однолинейных систем. К сожалению, данный метод не удается применить для исследования систем не пуассоновским входящим потоком (ММРР, МАР, рекуррентный поток), для таких потоков необходимо освоение новых методов исследования.
В данной статье приводится исследование СМО с повторными обращениями, на вход которой поступает поток ММРР|М|да. С помощью метода начальных моментов найдены основные характеристики системы, такие как среднее число занятых приборов в системе и второй момент числа занятых приборов. Кроме того, для более детального исследования процессов, указанных выше, предложен метод асимптотического анализа.
Рассмотрим систему массового обслуживания с неограниченным числом приборов (рис. 1), на вход которой поступает марковский модулированный поток (ММРР), управляемый цепью Маркова к(і), заданной матрицей инфинитезимальных характеристик Q = ||^-|| и матрицей условных интенсивностей Л [16].
Продолжительность обслуживания заявки имеет экспоненциальное распределение с параметром ц [20]. Поступившая заявка занимает любой из свободных приборов, завершив обслуживание на котором, с вероятностью 1-г покидает систему или с вероятностью г возвращается в неё для повторного обслуживания [2,17].
Обозначим /(?) число занятых приборов в момент времени ?, п(?) - число повторных заявок, обратившихся за время ?, к(?) - состояние управляющей цепи Маркова. Полученный двумерный случайный процесс {к(0, /(?)} является марковским. Для его распределения вероятностей Р(к,/,п,?) = Р{к(?) = к,/(?) = /,п(?) = п} получаем систему дифференциальных уравнений Колмогорова
1. Постановка задачи
1-7*
-Г
і(і)
Рис.1. Система массового обслуживания ММРР|М|а> с повторным обращением к блоку
Введем частичные характеристические функции [21] вида
Н(к,и,м>,і) = ££ЄшЄппР(к,і,п,і) .
і п
Учитывая, что
дH (k, ы, і )
ды
= ] Х іе]ШР (к, і, і),
і=0
из (1) имеем следующую систему уравнений:
дН(£,и,?) + ^^_ ^дН(Ат,^^к( _ ^^^?)+ ^ (u,?).
ді ' ! ды к
Запишем данную систему в виде матричного уравнения
дН (ы, і)
уфк
где
ді
]ц(1-г)(е-]ы -1))^ = Н(ы,^[л(е]ы -1)
Н(ы,і) = [Н(1,ы,і),Н(2,ы,і),...,Н(К,ы,і)] :
0 . . 0 " Ч11 Чп . . Ч1К
Л = 0 Х 2 . . 0 , Q = Ч21 Ч22 . . Ч2К
0 0 . . Х К ] [Які Як 2 . . Якк ]
(2)
2. Метод начальных моментов
Используя уравнение (2), определим начальные моменты для числа занятых приборов в системе, для этого рассмотрим нулевой момент времени, в который имеем число повторных обращений w = 0 (так как число обслуженных заявок равно нулю), тогда (2) перепишем в следующем виде:
ди
Чтобы определить момент первого порядка, продифференцируем уравнение (3) по и:
ч _ ,и дН (и) ./ _ и чд2 Н (и) дН (и) г / ,и \
.(1 _+л.(1 _г)( _ =^1[л( _О
обозначив
дН ( и )
(3)
+ ]егыН(ы), (4)
ды
= 7ті:
ы=0
получаем следующую систему линейных уравнений:
(1 _ г) = m1Q + ИЛ.
Здесь И = Н(0) - вектор стационарного распределения вероятностей управляющей цепи Маркова [16], которое определяется системой уравнений RQ = 0 и удовлетворяет условию нормировки КЕ = 1. Решая данную систему, получаем выражение, определяющее вектор ш1 :
Ш1 = ИЛ[. (1 _ г) I _ Q]_1. (5)
Тогда значение среднего числа занятых приборов в системе определяется следующим выражением:
М {і (і)}= . 1 . ЯЛЕ.
Ц (1 - г)
Для определения момента второго порядка продифференцируем по и выражение (4), получим
- ]Ц (1 - г )е"
-]ы дН (ы )
ды
-2ц(1 - г)е
- ]ы
д2 Н (ы)
ды
,1 - г)(е-- - .ДИ
д 2Н (ы)
Л( -1) + Q] + 2]е]ы Л - е]ыН (ы).
ды
Тогда, обозначив
д 2Н (ы)
= ] т2 ,
ы=0
получим
т2 [Q - 2ц(1 - г)!] + т1 [2Л + ц(1 - г)I] + ЯЛ = 0,
(6)
где I - единичная диагональная матрица. Подставляя в (5) выражение (6) и решая данное уравнение относительно величины ш 2, имеем
т.
: КЛ{[>.(1 _г)I{ [2Л + .(1 _ г)1 ] _ I} _2.(1 _ г) 1}_ .
Тогда второй момент числа занятых приборов в системе имеет вид
М {/2 (?)} = ш2Е = RЛ{[Q-.(1 _ г)]_* [2Л + .(1 _ г)] _ I}{Q _ 2.(1 _ г)1}_1 Е .
3. Метод асимптотического анализа
Для более полного исследования применим метод асимптотического анализа, заключающийся в нахождении аппроксимации характеристической функции числа занятых приборов в системе ММРР|М| да при определенных условиях. Для нашей системы будем рассматривать условие растущего времени обслуживания [16], т.е. высокой загрузки системы.
3.1. Асимптотика первого порядка
Обозначим
ц = є, ы = гу, Н (ы ) = Ц (у, 8). Перепишем (3) с учетом введенных обозначений:
(1 - г) (]8у -1)^1^ 8) = Е (у, 8) [Л(е]8у -1) + Q
Тогда при є ^ 0 имеем уравнение
решение которого имеет вид
0 = Е (у) Q, *1 (У ) = ЯФі ( У).
(7)
(8)
Домножая уравнение (8) на единичный вектор Е и выполняя все необходимые преобразования, имеем
](1 - г)(е-]8у -і)^1^ 8) Е = *1 (у, є)(е]8у - і)ЛЕ .
Раскладывая в этом уравнении экспоненты в ряд, а также разделив обе части уравнения на є и полагая є ^ 0, получаем
, . дЕ (у) , .
(1 - г)у—^Е = ]у*1 (у)ЛЕ . (10)
ду
Подставляя (9) в систему (10) и учитывая условие ЯЕ = 1, получаем дифференциальное уравнение первого порядка для нахождения функции Ф1(у):
Е=—]— Ф1(у)ЛЕ. ду (1 - г) 1
Откуда, учитывая начальное условие Ф1 (0) = 1, имеем
]у
Ф1 (y ) = exp {—- к- j, к- = R • A • E.
Следовательно, учитывая (9),
Fi( y) = R exp Ki j.
В силу замены (7) и предыдущего равенства можно записать следующее асимптотическое равенство:
jy
Н (и ) = F1 ^ 8) * F1 (y) = R • exp {-—- К- j .
Поэтому для характеристической функции стационарного процесса i(t) запишем
M {eMt}} = Н (u )E * exP \j Ki j = exP {j Kj.
Данное равенство будем называть асимптотикой первого порядка для СМО ММРР|М|да с повторными обращениями.
3.2. Асимптотика второго порядка Выполним следующую замену:
Н(и) = Н (и)expU-.
Тогда уравнение (3) для Н2 (и) будет иметь вид
j-ц(1 - r)(.-( - l)дн2^u) = H2 (u)[q + Л(.ju -1) + к- (.- -l)1
(11)
ды
Обозначим
ц = є2, u = єу, H2 (u) = F2 (y, є) (12)
и из(11) получим
]є(1 - г)(е-]єу -1); є) = *2 (у, є)^ + Л( -1) + Кі (е- - і)т'
(13)
тогда при е ^ 0 имеем уравнение
о = Р2 (у) Q.
Следовательно, решение (у, е) в уравнении (13) можно записать в виде разложения
г2 (У, е) = Ф2 (У) [ К + ]£У*2 ] + О (е2). (14)
Подставляя данное выражение в (13) и раскладывая в ней экспоненты в ряд, получаем уравнение для нахождения функции Г2
И (Л _ к^ ) + ^ = 0.
Чтобы найти функцию Ф2 (у), домножим обе части уравнения (13) на единичный вектор Е,
далее, раскладывая в данном выражении экспоненты в ряд и подставляя (14), получаем дифференциальное уравнение первого порядка
(1 _ г )“;“■ = ] 2уФ2 (у)[К • Л • Е + Г2 (Л _ )Е] .
ду
Решение данного уравнения, удовлетворяющее начальному условию Ф2 (0) = 1, имеет вид
Ф2 (у) = ехр<
(]у)
2
2 (1 _ г ) 2
где к 2 = к1 + ^ (Л _ к1I) Е .
Тогда в силу замен (12) запишем асимптотическое равенство для Н2 (и):
Н2 (и) = Г2 (У, е) * Г2 (у) = К • ехР2((1у)г) К2 ^ И • ехр)-. ^ . (15)
Это равенство будем называть асимптотикой второго порядка.
Учитывая (11) и (15), имеем следующее выражение для Н(и):
Н(и)= Н2(и)• ехри^к.} = К.ехр\*-^ + -(’^ к2
I ^ | 1 / ч | 5
1 _г . \ 11 _г . 2(1 _ г) . I
поэтому для характеристической функции числа занятых приборов /(?) в системе получим
]и к , (р )2 к1
М {е]ы,( )} = Н (ы )Е = ехр
- г ц 2 (1 - г) ц I
Из последнего выражения следует, что стационарное распределение вероятностей числа занятых приборов в системе ММРР|М|да можно аппроксимировать гауссовским распределением со следующими параметрами:
к2
=м {і (і )}=ціГТ) ■ °2 = м {(і (і)-" )2}=
-г)
4. Область применения асимптотических результатов
Выясним, насколько результаты, полученные с помощью асимптотического анализа, близки к аналитическим результатам.
Пример 1.
Используем следующие значения параметров:
"-0,5 0,3 0,2 ' "1 0 0 '
Q = 0,2 -0,6 0,4 , Л = 0 0,6 0
[ 0,5 0,3 -0,8 0 0 0,8
Для заданных значений, изменяя параметр ., получаем следующие результаты (табл. 1).
Т а б л и ц а 1
Сравнение асимптотических и аналитических результатов
1 0,5 0,25 0,1 0,05
М Б М Б М Б М Б М Б
Асимптотические результаты 4,067 4,244 8,133 8,488 16,267 16,976 40,667 42,440 81,333 84,881
Аналитические результаты 4,067 4,210 8,133 8,450 16,267 16,936 40,667 42,399 81,333 84,839
Вычислим относительную погрешность между величиной дисперсии, полученной методом моментов, и асимптотическим методом (табл. 2).
Т а б л и ц а 2
Относительная погрешность
1 0,5 0,25 0,1 0,05
А 0,82% 0,45% 0,24% 0,10% 0,05%
Следует отметить, что при данных входных параметрах максимальная относительная погрешность составляет менее 1%. Отметим, что при уменьшении параметров времени обслуживания результаты, полученные асимптотическим методом, приближаются к аналитическим результатам, что позволяет применять асимптотический алгоритм для широкого круга задач.
Пример 2.
Рассмотрим еще один пример с большой интенсивностью входящего потока (табл. 3). Используем следующие значения параметров:
"-0,5 0,3 0,2 ' "5 0 0" ( 5 0 0 1 ' 1 0 0 1 (11
Q = 0,2 -0,6 0,4 , Л = 0 8 0 , л = 0 8 0 , Е = 0 1 0 , I = 1
[ 0,5 0,3 -0,8 0 0 10 ,0 0 10V ч 0 0 1V л
Т а б л и ц а 3
Сравнение асимптотических и аналитических результатов
1 0,5 0,25 0,1 0,05
М Б М Б М Б М Б М Б
Асимптотические результаты 36,667 59,383 73,333 118,765 146,667 237,531 336,667 593,827 733,330 1187,654
Аналитические результаты 36,667 55,354 73,333 114,343 146,667 232,882 336,667 589,031 733,330 1182,807
Вычислим относительную погрешность между величиной дисперсии, полученной методом моментов, и асимптотическим методом.
Т а б л и ц а 4
Относительная погрешность
1 0,5 0,25 0,1 0,05
А 7,2S% 3,S7% 2,00% 0,si% 0,41%
Заметим, что в данном примере относительная погрешность, приведенная в табл. 4, оказалась в разы больше погрешности в примере 1; это говорит о том, что на асимптотические результаты влияет не только продолжительность времени обслуживания, но и входящие параметры, а именно интенсивность входящего потока.
Заключение
В настоящей работе была построена математическая модель обслуживания заявок в системе ММРР|М| да, а также определены аналитические выражения для нахождения первого и второго моментов, характеризующих число занятых приборов в системе и их асимптотического приближения. Проведен численный анализ полученных результатов.
ЛИТЕРАТУРА
1. Моисеева С.П., Захорольная И.А. Математическая модель параллельного обслуживания кратных заявок с повтор-
ными обращениями // Автометрия. 2011. Т. 47, № 6. С. 51-58.
2. Морозова А.С., Моисеева С.П., Назаров А.А. Исследование экономико-математической модели влияния ценовой
скидки для постоянных клиентов на прибыль коммерческой организации // Вестник Томского государственного университета. 2006. № 293. С. 49-52.
3. Носова М.Г. Автономная немарковская система массового обслуживания и ее применение в задачах демографии :
дис. ... канд. физ.-мат. наук. Томск, 2010. 204 с.
4. Кёнинг Д., Рыков В., Штоян Д. Теория массового обслуживания. М. : Московский институт нефтехимической и
газовой промышленности, 1979. 112 с.
5. Бочаров П.П., Печинкин А.В. Теория массового обслуживания. М. : Изд-во РУДН, 1995. 520 с.
6. Parulekar M., Makowski A.M. Tail probabilities for M/G/да input processes (I): Preliminary asymptotics // Queueing Sys-
tems. 1997. V. 27, is. 3-4. P. 271-296.
7. Baltzer J.C. On the fluid limit of the M/G/да queue queueing systems // Theory and applications. August 2007. V. 56, is. 3-4.
Р. 255-265.
8. Leland W.E., Willinger W., Taqqu M.S., Wilson D.V. On the self-similar nature of Ethernet traffic, ACM SIGCOMM //
Computer Communication Review. 1995. V. 25. P. 202-213.
9. Klemm A., Lindemann C., Lohmann M., Modeling I.P. Traffic Using the Batch Markovian Arrival Process (extendend ver-
sion) // Performance Evaluation. 2003. V. 54. P. 149-173.
10. Baum D. The infinite server queue with Markov additive arrivals in space // Proceedings of the international conference «Probabilistic analysis of rare events». Riga, Latvia, 1999. P. 136-142.
11. Breuer L., Baum D. The Inhomogeneous BMAP/G/infinity queue // Proceedings 11th GI/ITG Conference on measuring, modelling and evaluation of computer and communication systems (MMB 2001). Aachen, Germany, 2001. P. 209-223.
12. Jayawardene A.K., Kella O. M/G/да with alternating renewal breakdowns // Queueing Systems. 1996. V. 22, is. 1-2. P. 79-95.
13. НазаровА.А., ТерпуговА.Ф. Теория массового обслуживания. Томск : Изд-во НТЛ, 2005. 228 c.
14. Iglegart D.L. Limit diffusion approximations for the many server queue and the repairman problem // J. Appl. Prob. 1965. Ко. 2. P. 429-441.
15 Reynolds J.F. Some results for the bulk-arrival infinite-server Poisson queue. // Oper. Res. 1968. V. 16. P. 186-189.
16. Назаров А.А., Моисеева С.П. Метод асимптотического анализа в теории массового обслуживания. Томск : Изд-во НТЛ, 2006. 112 с.
17. Жидкова Л.А., Моисеева С.П. Математическая модель потоков покупателей двухпродуктовой торговой компании в виде системы массового обслуживания с повторными обращениями к блокам // Известия Томского политехнического университета. 2013. Т. 322, № 6. C. 5-9.
18. Жидкова Л.А., Моисеева С.П. Исследование систем параллельного обслуживания кратных заявок простейшего потока // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 4(17). С. 49-55.
19. Моисеева С.П., Ананина И.А., Назаров А.А. Исследование потоков в системе M|GI|o> с повторными обращениями методом предельной декомпозиции // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2009. № 3(8). С. 56-66.
20. Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. М. : КомКнига, 2005. 408 с.
21. Хинчин А.Я. Работы по математической теории массового обслуживания. М. : Физматгиз, 1963. 236 с.
Жидкова Любовь Александровна. E-mail: zhidkovala@mail.ru Моисеева Светлана Петровна. E-mail: smoiseeva@mail.ru
Томский государственный университет Поступила в редакцию 28 октября 2013 г.
Zhidkova Lyubov A., Moiseeva Svetlana P. (Tomsk State University, Russian Federation).
Investigation of the queueing system MMPP|M|<x> with repeated service.
Keywords: Queueing system with repeated service; Markov modulated process; method of asymptotic analysis.
In this paper the investigation of the queueing infinity-server system with the Markov modulated arrival process and repeated service is presented. The Markov modulated Poisson process is defined by the infinitesimal generator Q = | |qj|| and conditional rates
Л . Each incoming customer occupies any of the vacant server for stochastic time interval distributed according to the exponential law with the parameter ц. After service the customer leaves system with probability r-1 and comes back for repeated service with probability r.
The goal of the research is to obtain the distribution parameters of the number of busy servers in stationary mode. The system of differential equations for partial characteristic functions of the number of busy servers is obtained in the following matrix form:
£5M+„(1_ r)((u - ^ = H(u,t)[a( -1) + Q] .
The formulas for the average and variance of the number of busy servers are found by the moments method. The detailed research is performed by the asymptotic analysis method when the service time increases. It is shown that the asymptotic characteristic function of the number of busy servers in stationary mode can be approximated by the Gaussian distribution with the following parameters:
"=M!'(' )!=^, = M {()-a f)= •
where values к1, к2 are defined by the expressions
к1 = R • Л • E ,
к2 = к1 + f2 (Л - K1I) E .
Here the vector f2 is the decision of the linear equations system
R (Л - K1I) + f2Q = 0 .
The applicability range of the asymptotic approximation is determined by numerical calculations. It is shown that asymptotic results depend on the service time and arrivals rate.
REFERENCES
1. Moiseeva S.P., Zakhorol'naya I.A. Matematicheskaya model' parallel'nogo obsluzhivaniya kratnykh zayavok s povtornymi
obrashcheniyami. Avtometriya. 2011. V. 47. No. 6. P. 51-58. [Moiseeva S.P., Zakhorol'naya I.A. Mathematical model of parallel retrial queueing of multiple requests. Optoelectronics, Instrumentation and Data Processing. V. 47. No. 6. P. 567572. (2011).]
2. Morozova A.S., Moiseeva S.P., Nazarov A.A. Issledovanie ekonomiko-matematicheskoy modeli vliyaniya tsenovoy skidki
dlya postoyannykh klientov na pribyl' kommercheskoy organizatsii. Vestnik Tomskogo gosudarstvennogo universiteta. 2006. No. 293. P. 49-52. [Morozova A.S., Moiseeva C.P., Nazarov A.A. Investigation of the economic-mathematical model of discount for patrons influence on income of trading company. Tomsk State University Journal. No. 293. P. 4952. (2006).]
3. Nosova M.G. Avtonomnaya nemarkovskaya sistema massovogo obsluzhivaniya i ee primenenie v zadachakh demografii :
dis. ... kand. fiz.-mat. nauk. Tomsk, 2010. 204 p. [Nosova M.G. Autonomous non-markovian Queueing system and appli-cationin tasks of demography: PhD. Thesis. Tomsk State University , Tomsk. (2010).]
4. Kening D., Rykov V., Shtoyan D. Teoriya massovogo obsluzhivaniya. M. : Moskovskiy institut neftekhimicheskoy i
gazovoy promyshlennosti, 1979. 112 p. [Konig D., Rykov V., Stoyan D. Queueing Theory. М.: Moscow Institute of Oil and Gas. (1979).]
5. Bocharov P.P., Pechinkin A.V. Teoriya massovogo obsluzhivaniya. M. : Izd-vo RUDN, 1995. 520 p. [Bocharov P.P. Pechinkin A.V. Queueing Theory. M.: Peoples’ Friendship University of Russia (PFUR). (1995).]
6. Parulekar M., Makowski A. M. Tail probabilities for M/G/» input processes (I): Preliminary asymptotics. Queueing Sys-
tems. V. 27. Issue 3-4. P. 271-296. (1997).
7. Baltzer J.C. On the fluid limit of the M/G/» queue queueing systems. Theory and applications. V. 56. Issue 3-4. P. 255-
265. (August 2007).
8. Leland W.E., Willinger W., Taqqu M.S., Wilson D.V. On the self-similar nature of Ethernet traffic", ACM SIGCOMM.
Computer Communication Review. V. 25. P. 202-213. (1995).
9. Klemm A., Lindemann C., Lohmann M., Modeling IP Traffic Using the Batch Markovian Arrival Process (extendend ver-
sion). Performance Evaluation. V. 54. P. 149-173. (2003).
10. Baum D. The infinite server queue with Markov additive arrivals in space. Proceedings of the international conference “Probabilistic analysis of rare events”. Riga, Latvia. P. 136-142. (1999).
11. Breuer L., Baum D. The Inhomogeneous BMAP/G/infinity queue. Proceedings 11th GI/ITG Conference on measuring, modelling and evaluation of computer and communication systems (MMB 2001). Aachen, Germany. P. 209-223. (2001).
12. Jayawardene A.K., Kella O. M/G/» with alternating renewal breakdowns. Queueing Systems. V. 22. Issue 1-2. P. 79-95. (1996).
13. Nazarov A.A., Terpugov A.F. Teoriya massovogo obsluzhivaniya. Tomsk : Izd-vo NTL, 2005. 228 c. [Nazarov A.A., Terpugov A.F. Queueing Theory. Tomsk: NTL. (2005).]
14. Iglegart D. L. Limit diffusion approximations for the many server queue and the repairman problem. J. Appl. Prob. No. 2. P. 429-441. (1965).
15 Reynolds J.F. Some results for the bulk-arrival infinite-server Poisson queue. Oper. Res. V. 16. (1968).
16. Nazarov A.A., Moiseeva S.P. Metod asimptoticheskogo analiza v teorii massovogo obsluzhivaniya. Tomsk : Izd-vo NTL, 2006. 112 s. [Nazarov A.A., Moiseeva C.P. Method of asymptotic analyze on Queueing Theory. Tomsk: NTL. (2006).]
17. Zhidkova L.A., Moiseeva S.P. Matematicheskaya model' potokov pokupateley dvukhproduktovoy torgovoy kompanii v vide sistemy massovogo obsluzhivaniya s povtornymi obrashcheniyami k blokam. Izvestiya Tomskogo politekhnicheskogo universiteta. 2013. V. 322, no. 6. P. 5-9. [Zhidkova L.A. Moiseeva S.P. Mathematical model of consumer traffic in two -commoditty commercial company in the form of queueing system with repeated block access. The Bulletin of Tomsk Polytechnic University. V. 322. No. 6. P. 5-9. (2013).]
18. Zhidkova L.A., Moiseeva S.P. Issledovanie sistem parallel'nogo obsluzhivaniya kratnykh zayavok prosteyshego potoka. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel’naya tekhnika i informatika. 2011. No. 4(17). P. 49-55. [Zhidkova L.A., Moiseeva S.P. Investigation of the parallel service system with multiple claims of the Poisson process. Tomsk State University Journal of Control and Computer Science. No. 4(17). P. 49-55. (2011).]
19. Moiseeva S.P., Ananina I.A., Nazarov A.A. Issledovanie potokov v sisteme M|GI|» s povtornymi obrashcheniyami metodom predel'noy dekompozitsii. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel’naya tekhnika i informatika. 2009. No. 3(8). P. 56-66. [Moiseeva S.P., Ananina I.A., Nazarov A.A. Research of streams in system M|GI|» with repeated references the method of limiting decomposition. Tomsk State University Journal of Control and Computer Science. No. 3(8). P. 56-66. (2009).]
20. Gnedenko B.V., Kovalenko I.N. Vvedenie v teoriyu massovogo obsluzhivaniya. M. : KomKniga, 2005. 408 p. [Gnedenko B.P.,Kovalenko I.N. Introduction in Queueing Theory. M.: KomKniga. (2005).]
21. Khinchin A.Ya. Raboty po matematicheskoy teorii massovogo obsluzhivaniya. M. : Fizmatgiz, 1963. 236 s. [Hinchin A.Y. Papers about mathematical theory of Queueing. M.: Fizmatgiz. (1963).]
б2
