Асимптотический анализ системы массового обслуживания map/gi/∞ с высокоинтенсивным входящим потоком

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА 2015 Управление, вычислительная техника и информатика № 3

УДК 519.872

DOI 10.17223/19988605/32/6

А.Н. Моисеев

АСИМПТОТИЧЕСКИЙ АНАЛИЗ СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ MAP/GI/ю С ВЫСОКОИНТЕНСИВНЫМ ВХОДЯЩИМ ПОТОКОМ

Результаты получены в рамках выполнения государственного задания Министерства образования и науки Российской Федерации № 1.511.2014/К.

Представлен анализ системы массового обслуживания с неограниченным числом приборов, входящим MAP-потоком и произвольным распределением времени обслуживания. Показано, что в условиях растущей интенсивности входящего потока стационарное распределение вероятностей числа заявок в системе является асимптотически гауссовским. Получены параметры соответствующей гауссовской аппроксимации. Кроме того, получена аппроксимация более высокого (третьего) порядка, установлены области применимости обеих аппроксимаций.

Ключевые слова: система массового обслуживания; MAP-поток; неограниченное число приборов; асимптотический анализ.

Системы массового обслуживания являются адекватными моделями современных систем передачи и обработки информации [1], социально-экономических систем [2, 3], систем управления транспортными потоками [4] и т.д. В научной литературе, как правило, исследуются системы с пуассоновским входящим потоком, однако имеются исследования [5], которые показывают, что, например, потоки сообщений в современных сетях передачи данных не являются пуассоновскими, а более адекватной моделью таких потоков является так называемый MAP-поток (Markovian Arrival Process) [6-8]. Таким образом, исследование моделей массового обслуживания с входящим MAP-потоком является актуальной научной проблемой.

В настоящей работе представлен анализ системы массового обслуживания с входящим MAP-потоком, неограниченным числом приборов и произвольным обслуживанием, которая в нотации Кендалла [9] обозначается как MAP/GI/да. Исследование проведено в асимптотическом условии растущей интенсивности входящего потока, которое относится к классу широко используемых условий «heavy traffic» [10-12]. В результате получена гауссовская аппроксимация распределения вероятностей числа заявок в системе. Этот результат согласуется с известными результатами, связанными с диффузионной аппроксимацией [10, 13, 14], но предоставляет более удобную форму для использования. Кроме того, изложенная в настоящей работе методика исследования позволяет получить аппроксимацию более высокого порядка точности, которая обеспечивает достаточно малую погрешность и при небольшой интенсивности входящего потока.

С другими, близкими по тематике исследованиями в области анализа систем массового обслуживания и случайных потоков можно ознакомиться в работах [15-18].

1. Постановка задачи

Рассмотрим систему массового обслуживания с неограниченным числом приборов, на вход которой поступает высокоинтенсивный MAP-поток, время обслуживания заявки в системе является случайной величиной с функцией распределения B(t). Входящий MAP-поток задан представлением [8] (ND0, NDi), где D0 и Di - квадратные матрицы порядка M, а скалярный мультипликатор N используется для того, чтобы выделить асимптотическую составляющую высокой интенсивности этого потока. Дело в том, что интенсивность MAP-потока, заданного таким образом, составляет AN, где величина А определяется выражением

Х = 0D1e.

56

Здесь e - вектор-столбец, состоящий из единиц, 0 - вектор-строка, элементы которой задают закон стационарного распределения вероятностей цепи Маркова m(t), управляющей MAP-потоком. Этот вектор удовлетворяет следующей системе уравнений:

№„ + D! ) = 0,

[0е = 1.

Таким образом, интенсивность описанного потока неограниченно возрастает при N ^ да. Потоки такого рода называют высокоинтенсивными [18].

Обозначим i(t) - количество заявок (или, что то же самое, количество занятых приборов) в системе в момент времени t. Ставится задача нахождения распределения вероятностей процесса i(t) в произвольный момент времени в стационарном режиме функционирования системы.

2. Метод динамического просеивания. Уравнения Колмогорова

Для составления уравнений Колмогорова, которые являются основным механизмом исследования случайных процессов, необходимо провести так называемую марковизацию исследуемого процесса i(t),

т. е. преобразовать его таким образом, чтобы полученный в результате преобразований процесс был марковским. Стандартным инструментом такой марковизации является метод дополнительных компонент [9, 19, 20]. Для рассматриваемого процесса i(t) этими дополнительными компонентами будут состояние m(t) цепи Маркова, управляющей входящим MAP-потоком, и значения остаточного времени обслуживания для каждой заявки, находящейся в системе в момент t. Таким образом, получаем многомерный марковский процесс с переменным числом компонент.

Прямое исследование такого процесса не представляется возможным ввиду переменного и потенциально неограниченного количества его компонент. В связи с этим в настоящей работе применяется метод динамического просеивания (просеянного потока) [21], который позволяет решить указанную задачу. Основная идея этого подхода заключается в том, чтобы заменить исследование многомерного процесса с переменным числом компонент на исследование двумерного марковского процесса {n(t), m(t)}, где n(t) - считающий процесс просеянного потока, который конструируется следующим образом. Зафиксируем некоторый момент времени T. Будем считать, что заявка, поступившая в систему в момент времени t < T, с вероятностью S(t) = 1 - B(T - t) формирует событие просеянного потока, а с вероятностью 1 - S(t) эта заявка не рассматривается. Через n(t) обозначим число событий просеянного потока, наступивших до момента времени t. Тогда, если в начальный момент времени t0 < T система свободна, получаем, что для момента времени T имеет место равенство

P{i(T) = k} = P{n(T) = k)}

для любых целых k > 0. Таким образом, найдя распределение вероятностей значений процесса n(t) в момент времени t и полагая t = T, в силу произвольности выбора T получим искомое стационарное распределение вероятностей процесса i(t).

Для распределения вероятностей P(n, m, t) = P{n(t) = n, m(t) = m} значений марковского процесса {n(t), m(t)} можно составить следующую систему дифференциальных уравнений Колмогорова:

N= | [P(n,/,t) (D ),m + P(n,I,t) (A ) m 0 - S(t)) + P(n -1,/,t) (A )m S(t)]

для n > 0, m = 1, M . Для частичной характеристической функции

H(u,m,t) = 2ejunP(n,m,t).

n=0

где j

=v-r,

эта система перепишется в виде

N dH t) =Щн [, t )(D0 ) + H (u,l, t )[ ) (1 - S (t)) + H (u,I, t )(DX \mS (t )eJu

/=1

Приводя подобные и используя обозначения D = D0 + Db H(u, t) = {H(u, 1, t), ..., H(u, M, t)}, получим матричное дифференциальное уравнение

N= H(u,t)[D + Dx (eju - 1)S(t)]

57

с начальным условием

H

Так как прямое решение уравнения (3) не представляется возможным, для решения задачи (3)-(4) воспользуемся методом асимптотического анализа [21] в условии растущей интенсивности входящего потока, т.е. при N ^ да.

3. Асимптотический анализ

3.1. Асимптотика первого порядка Выполним в выражениях (3)-(4) следующие замены:

N = Б, u = sw , H(u, t) = Fi(w,t, в).

Тогда задача (3)-(4) примет вид

5Ft( w, t, s) ' dt

Fi( w, t, s)[D + Di (ejw - l)s (t)],

Fi(w, to, e) = 0.

Относительно асимптотического решения F1(w, t) = lim F1(w, t, в) этой задачи имеет место следующее

утверждение.

Теорема 1. Асимптотическое решение задачи (5)-(6) имеет вид

F1(w, t) = 0 expj jwXjS(x)dx l,

где 0 и X определяются выражениями (2) и (1) соответственно.

Доказательство выполним в два этапа.

Этап 1. Положим в (5) e ^ 0, получим

F1(w, t)D = 0.

Сравнивая это равенство с первым уравнением (2), можем сделать вывод, что F1(w,t) может быть представлена в виде

F1

где Ф1(w, t) - некоторая скалярная функция, в силу (6) удовлетворяющая условию Ф^, t0) = 1.

Этап 2. Домножим

501

dt

= Ф1( w, t) jwXS (t),

решение которого с учетом начального условия дает

Ф^,t) = expj jw XJS(x)dxj.

Подставляя данное выражение в (8), получаем (7). Теорема доказана.

3.2. Асимптотика второго порядка

Обозначим через H2(u, t) многомерную функцию, удовлетворяющую выражению

H(u,t) = H2(u,t)expj juNXJS(x)dxj .

Подставляя это выражение в (3)-(4), получим следующую задачу Коши:

, N Ш^t} + H2 (u,t)juXS(t) = H2 (u,t)[d + D1 (eju - 1))t)]

H 2(u, t0) = 0.

t

58

Выполним в (10) следующие замены:

1 2

— = 8 , u = Sw , H 2 (u, t) = F2(w,t,s). С использованием этих обозначений (10) перепишется в виде

2 d¥2(w,t,s) + f2(w,t,s) jswXS(t) = F2(w,t,s)[d + D1(Sw - 1)(t) dt

S

F2 (w, to, S) = 0.

Относительно асимптотического решения F2(w, t) = limF2(w,t, s) этой задачи имеет место следующее

S^0

утверждение.

Теорема 2. Асимптотическое решение задачи (12) имеет вид

F2( w, t) = 0 exp-

(jw)2

AJ S (x)dx + k| S 2(x)dx

где

к = 2g

а вектор-строка g удовлетворяет линейному матричному уравнению

gD = 0

здесь I - единичная матрица.

Доказательство этой теоремы выполним в три этапа.

Этап 1. В (12) выполним предельный переход s ^ 0, получим следующую систему уравнений:

[F2(w, t )D = 0,

[F2(w, t0) = 0.

Отсюда с учетом (2) следует, что F2(w, t) можно записать в виде

F2

где Ф2^, t) - некоторая скалярная функция, удовлетворяющая условию Ф2^, t0) = 1.

Этап 2. С учетом (16) функцию F2(w, t, s) представим в виде разложения

F2

где g - некоторая вектор-строка, а O(s2) - вектор-строка, состоящая из величин порядка малости s2 и выше. Подставим (17) и разложение esw = 1 + jsw +O(s2) в первую строку (12), получим

0 • jswXS(t) = 0D + 0D1 • jswS(t) + gD • jswS(t) + O(s 2).

Выполняя в этом равенстве предельный переход s ^ 0, получаем матричное уравнение относительно вектора g

gD = 0(AI - D1),

которое совпадает с (15).

Этап 3. Умножая обе части уравнения (12) на вектор е, используя (17) и разложение

ejw = 1 + jsw + + O(s3),

2

в результате несложных преобразований с учетом обозначения (14) получаем следующее дифференциальное уравнение относительно функции Ф2^, t):

=ф 2( w, ,Ш- [ (t) + kS 2(t)]l.

dt *' " [ 2

Решение этого уравнения при имеющемся начальном условии Ф2^, t0) = 1 дает выражение

Ф 2( w, t) = exp-

(jw)2

2

AJ S (x)dx + k| S 2(x)dx

подставляя которое в (16), получаем (13). Теорема доказана.

2

0

t

0

59

3.3. Гауссовская аппроксимация стационарного распределения вероятностей числа заявок в системе

Выполнив в выражении (13) замены, обратные к (11) и (9), получим следующее выражение для частичной характеристической функции H(u, t) числа событий, наступивших в просеянном потоке до момента времени t:

H(u, t) = 0 expJ juNIJ S(x)dx +

(ju )2

N1J S (x)dx + NkJ S 2 (x)dx

t0 t0

2

0

Умножим обе части этого равенства на е. Тогда, возвращаясь к исходному исследуемому процессу i(t) числа заявок в системе, подставим t = T, а в качестве начального момента (когда система свободна) выберем t0 = -да, получим следующее выражение для характеристической функции h(u) числа заявок в системе в стационарном режиме:

h(u) = expJ juNlb + j) [N1b + NkP] l,

где b = J(1 - B(t))t есть среднее время обслуживания, а Р = j(1 - B(t))2dt.

00

Таким образом, стационарное распределение числа заявок в системе MAP/GI/да является асимптотически гауссовским с математическим ожиданием Nkb и дисперсией [Nb + NkP].

3.4. Асимптотика третьего порядка

Применение гауссовской аппроксимации (19) достаточно удобно, однако, как будет показано в

п. 4, допустимая погрешность достигается при больших значениях N (30 и больше). Для меньших значений интенсивности входящего потока можно получить более точную аппроксимацию, используя асимптотику третьего порядка.

Итак, введем функцию H3(u, t), удовлетворяющую выражению

H 2 (u, t) = H 3(u, t) exp-

(ju )2

2

N1J S (x)dx + NkJ S 2(x)dx

подставим (20) в (10), получим следующую задачу Коши:

J N dt

H3(u t0) = 0.

В (21) выполним замены

- dH3

2

1S

N = s , u =SW , H3(u,t) = F3(w,t,e),

получим

3 dF3

dt

F3( w, t0, e) = 0.

Обозначим F3 (w, t) = limF3 (w,t, e) . Тогда имеет место следующая теорема.

E—>0

Теорема 3. Асимптотическое решение F3(w,t) задачи (22) имеет вид

F3(w, t) = 0 exp-

(jwf

1J S (x)dx + 3k J S 2(x)dx + p J S3 (x)dx

tn t* t*

где

0

0

E

6

60

ц = 6g2(Di - T!)e - 3Kg,

g определяется выражением (15), а вектор-строка g2 удовлетворяет линейному матричному уравнению

g 2 D = 0

К

--X|I -Di

+ Xg.

Доказательство этой теоремы выполним в три этапа.

Этап 1. Выполним в (22) предельный переход s ^ 0, получим

^(w, t )D = 0,

(^(w, to) = 0.

Отсюда с учетом (2) делаем вывод, что F3(w, t) можно представить в виде

F3(w, t) = 0Ф3(^, t),

где Ф3(^, t) - некоторая скалярная функция.

Этап 2. Решение задачи (22) будем искать в виде разложения

F3

0 + gi • jBwS

+ O(s3),

где g1 и g2 - некоторые вектор-строки. Подставим разложения (27) и (18) в (22). После несложных преобразований получаем равенство

jswXS

- 0 + Xg1

2 1

= jswS (t )[0D1 + S (t )[0D1 + g1D] + (/sw)2 S 2(t )[g1D1 + g 2D] + O(s 2).

Приравнивая выражения при одинаковых степенях s, получаем, во-первых, X0 = 0Dt + gtD, откуда следует gjD = 0(XI - Dj) , что совпадает с (15). С учетом этого вторые слагаемые левой и правой частей (28) уничтожаются, и мы получаем равенство

К

—0 + Xg = gD1 + g 2D,

откуда следует выражение (25).

Этап 3. Умножим каждую часть (22) на е, тогда с использованием обозначения f(w, t, s) = F3(w, t, s)e получаем

+ f (w,t,s)-|jswXS(t) +-(jsw)[xS(t) + kS2(t)]| = ] (w,t,s)D1e(ejsw - 1)(t),

f(^ to, s) = 1.

Подставляя сюда (27) и разложение

jsw , (jsw)2 (jsw)3 4.

e] = 1 + jsw +42—>— + —>— + O(s 4),

2 6

получаем

3 ao3

s

dt

+ Ф 3(w, t )|1 + ge

jsw +

(jsw)2

S (t) + g 2e(jsw)2 S 2(t) Nx

jsw +

(jsw)2

S (t) + - (jsw)2 S 2(t )| =

= Ф 3 (w, t )|0 + g

jsw +

(jsw)2

2

S(t) + g2 (jsw)2 S2(t) Id^

jsw +

(jsw)2 , (jsw)3

2

6

S (t) + O(s 4).

После несложных преобразований с учетом

= Ф 3 (w, t )| [ (t) + 3kS 2 (t) + MS3 (t) |

2

2

61

решение которого при имеющемся начальном условии Ф3(^, t0) = 1 дает

Ф 3( w, t) = exp-

СjwY

6

Я/ S (x)dx + 3kJ S 2 (x)dx + pj S 3(x)dx

to to to

l0 l0 l0

Подставляя это выражение в (26), получаем (23). Теорема доказана.

Делая обратные замены и аналогично п. 3.3 возвращаясь к характеристической функции числа заявок в системе в стационарном режиме, получим

h(u) = exp] juNkb + (ju)[Nb + NkP] + j)Nb + 3NkP + Npy],

-

где у = j(l - B(t))3 dt.

o

К сожалению, пользоваться аппроксимацией (29) не так удобно, как гауссовской (19) (для получения ряда распределения придется выполнить обратное преобразование Фурье от этого выражения), однако, как будет показано ниже, выражение (29) обеспечивает меньшую погрешность и более широкую область применимости асимптотических результатов.

4. Численные результаты

Демонстрацию точности аппроксимаций (19) и (29) проведем на конкретном численном примере. Пусть входящий MAP-поток задан представлением (ND0, NDi), где

'- 6,5o o,5o o,97 > '5,00 0,00 0,03N

Do = o,36 - 1,9o o,96 , D1 = 0,04 0,50 0,04

v o,19 o,54 - o,85 j v 0,01 0,06 0,05 j

Время обслуживания имеет гамма-распределение с параметрами масштаба и формы, равными 0,5. Пример составлен таким образом, что параметр X приблизительно равен 1, среднее время обслуживания также равно 1, и, таким образом, интенсивность входящего потока превышает интенсивность обслуживания в N раз.

0,5 0,4 -0,3 0,2 0,1

0

0 12 3 4

10

20

30

40

50

Рис. 1. Сравнение законов распределения вероятностей для эмпирического распределения

62

На рис. 1 представлены графики законов распределения вероятностей для эмпирического распределения, составленного на основе результатов имитационного моделирования, и распределений, полученных на основе аппроксимаций (19) и (29). Из графиков видно, что при достаточно больших значениях N (N > 10) распределение по форме близко к нормальному.

Для числовой оценки точности аппроксимаций воспользуемся расстоянием Колмогорова [22]. В таблице приведены значения расстояния Колмогорова для различных значений величины N. Если в качестве допустимой погрешности выбрать расстояние Колмогорова, равное 0,03, то из таблицы видно, что третья асимптотика дает достаточно точную аппроксимацию уже при N = 10, в то время как гауссовская аппроксимация дает такую же точность только при N = 30.

Расстояние Колмогорова между эмпирическим распределением и аппроксимациями второго и третьего порядка для различных значений интенсивности входящего потока N

Интенсивность входящего потока N Расстояние Колмогорова

для аппроксимации второго порядка (гауссовской) для аппроксимации третьего порядка

1 0,242 0,186

5 0,099 0,048

10 0,060 0,027

30 0,029 0,013

100 0,017 0,007

Заключение

В работе представлено исследование системы массового обслуживания с неограниченным числом приборов, входящим MAP-потоком и произвольным распределением времени обслуживания. Исследование выполнено в асимптотическом условии растущей интенсивности входящего потока. Показано, что для достаточно больших значений параметра N интенсивности входящего потока стационарное распределение вероятностей числа заявок в системе аппроксимируется гауссовским (нормальным) распределением, параметры которого получены в работе. Для повышения степени точности выполнен асимптотический анализ третьего порядка и получены параметры аппроксимации третьего порядка для стационарного распределения вероятностей числа заявок в системе.

Числовые примеры показывают, что гауссовская аппроксимация дает достаточно малую погрешность при N > 30, где N - параметр, определяющий, во сколько раз интенсивность входящего потока превышает интенсивность обслуживания. Аппроксимация третьего порядка расширяет область применимости до значений N > 10, что существенно увеличивает точность полученных результатов. Однако при использовании полученных в работе аппроксимационных формул следует учитывать, что гауссовская аппроксимация более проста в применении. Таким образом, при N > 30 целесообразно использовать более удобную гауссовскую аппроксимацию, а для меньших значений N следует воспользоваться аппроксимацией третьего порядка.

Результаты могут быть использованы для анализа систем массового обслуживания вида MAP/GI/да, а методика исследования - для выполнения исследования многофазных систем и сетей массового обслуживания с входящим высокоинтенсивным MAP-потоком.

ЛИТЕРАТУРА

1. Грачев В.В., Моисеев А.Н., Назаров А.А., Ямпольский В.З. Многофазная модель массового обслуживания системы распреде-

ленной обработки данных // Доклады ТУСУР. 2012. № 2(26), ч. 2. С. 248-251.

2. Маталыцкий М.А., Статкевич С.Э. НМ-сети как новые стохастические модели прогнозирования доходов различных объек-

тов // Вестник ГрГУ. Сер. 5. Экономика. 2009. № 1. С. 107-115.

3. Гарайшина И.Р. Исследование математических моделей процессов государственного пенсионного страхования : дис. ...

канд. физ.-мат. наук. Томск : ТГУ, 2005. 148 с.

4. Гасников А.В. и др. Введение в математическое моделирование транспортных потоков : учеб. пособие / под ред. А.В. Гасни-

кова. М. : МФТИ, 2010. 362 с.

63

5. Heyman D.P., Lucantoni D. Modelling multiple IP traffic streams with rate limits // IEEE/ACM Trans. on Networking. 2003. V. 11,

No. 6. P. 948-958.

6. NeutsM.F. A versatile Markovian arrival process // J. of Appl. Prob. 1979. V. 16. P. 764-779.

7. Lucantoni D.M. New results on the single server queue with a batch Markovian arrival process // Stochastic Models. 1991. V. 7.

P. 1-46.

8. Breuer L., Baum D. An introduction to queueing theory and matrix-analytic methods. Springer, 2005. 271 p.

9. Kendall D.G. Stochastic processes occurring in the theory of queues and their analysis by the method of the imbedded Markov chains

// Ann. Math. Statist. 1953. V. 24(3). P. 338-354.

10. Whitt W. On the heavy-traffic limit theorem for GI/G/да queues // J. Adv. Appl. Prob. 1982. V. 14. P. 171-190.

11. Kingman J.F.C. On queues in heavy traffic // J. of the Royal Statistical Society. 1962. Series B 24(2). P. 383-392.

12. IglehartD.L., Whitt W. Multiple channel queues in heavy traffic // Adv. Appl. Prob. 1970. V. 2. P. 150-172.

13. IglehartD.L. Limit diffusion approximations for the many server queue and the repairman problem // J. of Appl. Prob. 1965. V. 2. P. 429-441.

14. Боровков А.А. О предельных законах для процессов обслуживания в многоканальных системах // Сибирский математический журнал. 1967. Т. 8, № 5. С. 983-1004.

15. Горбатенко А.Е., Назаров А.А. Исследование МАР-потока в условиях растущей интенсивности // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2008. № 3(4). C. 66-70.

16. Лопухова С.В. Исследование ММР-потока асимптотическим методом m-го порядка // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2008. № 3(4). C. 71-76.

17. Назаров А.А., Семенова И.А. Исследование системы MMP|GI|<» методом просеянного потока // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 4(17). C. 74-84.

18. Моисеев А.Н., Назаров А.А. Исследование системы массового обслуживания HIGI|GI|<» // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 2(23). С. 75-83.

19. Cox D.R. The analysis of non-Markovian stochastic processes by the inclusion of supplementary variables // Proc. Camb. Phil. Soc. 1955. V. 51. P. 433-441.

20. Henderson W. Alternative approaches to the analysis of M/G/1 and G/M/1 queue // J. Oper. Res. Soc. Japan. 1972. V. 15. P. 92-101.

21. Назаров А.А., Моисеева С.П. Метод асимптотического анализа в теории массового обслуживания. Томск : Изд-во НТЛ,

2006. 112 с.

22. Рыков В.В., Иткин В.Ю. Математическая статистика и планирование эксперимента : учеб. пособие. М. : МАКС Пресс,

2010. 308 с.

Моисеев Александр Николаевич, канд. техн. наук, доцент. E-mail: moiseev.tsu@gmail.com Томский государственный университет

Поступила в редакцию 20 апреля 2015 г.

Moiseev Alexander N. (Tomsk State University, Russian Federation).

Asymptotic analysis of queueing system MAP/GI/<x> with high-rate arrivals

Keywords: queueing system; Markovian arrival process; infinite-server system; asymptotic analysis.

DOI 10.17223/19988605/32/6

In the paper, the infinite-server queueing system with the Markovian arrival process and with the general form of service time distribution is considered. The analysis is performed under the asymptotic condition of arrivals’ rate infinite growth. The high intensive Markovian arrival process is defined by representation (ND0, ND@ where N is the scalar that characterizes the high rate of arrivals. Let D0 and D1 be square matrices. A rate of the arrivals is equal to NX. Here X = 0D^, 0 is a row vector of stationary probability distribution of the embedded Markov chain of the MAP, e is the column vector with entries all equal to 1. The service time for each customer is a random variable with the cumulative distribution function B(t).

The approximations of the 2nd and 3rd orders for the characteristic function of number of customers in the stationary regime are obtained. The approximation of the second order is the Gaussian distribution with mean NXb and variance [NXb + Мф]. Here

b = J(1 - B(t ))dt is the average service time, P = J(1 - B(t ))2 dt , к is defined by the expression к = 2g(Di - XI)e, and the row 0 0

vector g satisfies the linear matrix equation gD = 0(XI - D1), where D = D0 + D1 and I is the identity matrix.

The approximation of the third order is equal to

(ju d

h(u) = expt juNXb + ^-[NXb + NkP] + ^-[ЖЬ + 3NkP + Npy]

(ju )3

_____ да 3

where j = V— 1 , y = j(1 — B(t)) dt, p = 6g2(D1 - XI)e - 3Kg, the row vector g2 satisfies the following linear matrix equation

64

g 2 D = 0

к

--X\l - Dx

+ A,g.

The numerical analysis of the applicability of the obtained asymptotic results is also performed. It is shown that the asymptotic approximation of the second order (Gaussian approximation) provides enough accuracy when the arrivals rate is greater than rate of the servicing by 30 times or more (N> 30). On other hand, the third order approximation provides the similar accuracy for values N> 10. Therefore, it is reasonable to use the obtained Gaussian approximation as a more comfortable form for usage when N > 30, and we need take the approximation of the third order when a value of the parameter N is less than 30.

REFERENCES

1. Grachev, V.V., Moiseev, A.N., Nazarov, A.A. & Yampol'skiy, V.Z. (2012) Multi-stage queueing model of distributed data processing

system. Doklady TUSUR - Proceedings of TUSUR. 2(26/2). pp. 248-251. (In Russian).

2. Matalytskiy, M.A. & Statkevich, S.E. (2009) NM-seti kak novye stokhasticheskie modeli prognozirovaniya dokhodov razlichnykh

ob"ek-tov [НМ-networks as new stochastic models for forecasting of revenues of different objects]. Vestnik GrGU. Ser. 5. Ekonomika. 1. pp. 107-115.

3. Garayshina, I.R. (2005) Issledovanie matematicheskikh modeley protsessov gosudarstvennogo pensionnogo strakhovaniya [Investiga-

tion of mathematical models of the State pension insurance processes]. Physics and Mathematics Cand. Dis. Tomsk: Tomsk State University.

4. Gasnikov, A.V. et al. (2010) Vvedenie v matematicheskoe modelirovanie transportnykh potokov [Introduction into mathematical

modeling of traffic flows]. Moscow: MFTI.

5. Heyman, D.P. & Lucantoni, D. (2003) Modelling multiple IP traffic streams with rate limits. IEEE/ACM Trans. on Networking.

11 (6). pp. 948-958.

6. Neuts, M.F. (1979) A versatile Markovian arrival process. Journal of AppliedProb. 16. pp. 764-779.

7. Lucantoni, D.M. (1991) New results on the single server queue with a batch Markovian arrival process. Stochastic Models. 7. pp. 1-

46. DOI: 10.1080/15326349108807174

8. Breuer, L. & Baum, D. (2005) An introduction to queueing theory and matrix-analytic methods. Springer.

9. Kendall, D.G. (1953) Stochastic processes occurring in the theory of queues and their analysis by the method of the imbedded Mar-

kov chains. The Annals of Mathematical Statistics. 24(3). pp. 338-354.

10. Whitt, W. (1982) On the heavy-traffic limit theorem for GI/G/да queues. Journal of Advances in Applied Probability. 14. p. 171-190. DOI: 10.2307/1426738

11. Kingman, J.F.C. (1962) On queues in heavy traffic. Journal of the Royal Statistical Society. B 24(2). pp. 383-392.

12. Iglehart, D.L. & Whitt, W. (1970) Multiple channel queues in heavy traffic. Journal of Advances in Applied Probability. 2. pp. 150172.

13. Iglehart, D.L. (1965) Limit diffusion approximations for the many server queue and the repairman problem. Journal of Applied Probability. 2. pp. 429-441. DOI: 10.2307/3212203

14. Borovkov, A.A. (1967) O predel'nykh zakonakh dlya protsessov obsluzhivaniya v mnogokanal'nykh sistemakh [On limit laws for service processes in multi-channel systems]. Sibirskiy matematicheskiy zhurnal - Siberian Mathematical Journal. 8 (5). pp. 9831004.

15. Gorbatenko, A.E. & Nazarov, A.A. (2008) Issledovanie MAR-potoka v usloviyakh rastushchey intensivnosti [The research of MAP in conditions of growing intensity]. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i in-formatika - Tomsk State University Journal of Control and Computer Science. 3(4). pp. 66-70.

16. Lopukhova, S.V. (2008) Issledovanie MMR-potoka asimptoticheskim metodom m-go poryadka [Research of MMP-process by the asymptotical analysis method of the m-th order]. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 3(4). pp. 71-76.

17. Nazarov, A.A. & Semenova, I.A. (2011) Investigation of non-Markovian system of mass service with unlimited number of servers and input MMP-flow by the method of sifted flow. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(17). pp. 74-84. (In Russian).

18. Moiseev, A.N. & Nazarov, A.A. (2013) Investigation of queueing system HIGI|GI|<». Vestnik Tomskogo gosudarstvennogo universi-teta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(23). pp. 75-83. (In Russian).

19. Cox, D.R. (1955) The analysis of non-Markovian stochastic processes by the inclusion of supplementary variables. Mathematical Proceedings of the Cambridge Philosophical Society. 51. pp. 433-441. DOI: 10.1017/S0305004100030437

20. Henderson, W. (1972) Alternative approaches to the analysis of M/G/1 and G/M/1 queue. J. Oper. Res. Soc. 15. pp. 92-101.

21. Nazarov, A.A. & Moiseeva, S.P. (2006) Metod asimptoticheskogo analiza v teorii massovogo obsluzhivaniya [Asymptotic analysis method in queueing theory]. Tomsk: NTL.

22. Rykov, V.V. & Itkin, V.Yu. (2010) Matematicheskaya statistika i planirovanie eksperimenta [Mathematical statistics and design of experiments]. Moscow: MAKS Press.

65