УДК 519.872
ФИЗИЧЕСКАЯ ОЧЕРЕДЬ В ПРИКЛАДНОЙ ТЕОРИИ МАССОВОГО ОБСЛУЖИВАНИЯ
© 2016 А.С. Титовцев
Казанский национальный исследовательский технологический университет
Статья поступила в редакцию 09.09.2016
В статье рассмотрена проблема различных типов очередей, возникающих в системах массового обслуживания. Введено понятия физической очереди и получены общие математические формулы первых и вторых моментов основных дискретных и непрерывных случайных величин, характеризующих поведение физических очередей для различных моделей систем массового обслуживания смешанного типа.
Ключевые слова: физическая очередь, качество обслуживания, система массового обслуживания
Основные определения. На практике при эксплуатации различного рода технических объектов, работающих по принципу систем массового обслуживания (СМО), весьма часто приходится иметь дело не только с привычной и хорошо изученной математической очередью, но и с другими типами очередей. В частности, можно рассмотреть и другую числовую характеристику СМО, которая, например, в стандартном отчёте системы имитационного моделирования GPSS, носит название «очередь без нулевых входов». В данном случае под нулевым входом понимается такой приход заявки в систему, при котором в многоканальном устройстве имеется, по крайней мере, один свободный канал, и в этом случае вновь поступившая в систему заявка обслуживается немедленно. Таким образом, в отличие от обычной математической очереди, «очередь без нулевых входов» означает такую ситуацию, при которой на момент прибытия в систему новой заявки все каналы в обслуживающем устройстве могут быть заняты, но при этом очередь, как таковая, может и отсутствовать. В последнем случае поступившая заявка не встречает перед собой других заявок, а находится непосредственно перед обслуживающим устройством, в котором заняты все каналы обслуживания. Определённая таким образом «очередь без нулевых входов» должна рассчитываться с учётом только тех заявок, которые действительно ожидали начала обслуживания и без учёта заявок, которым не пришлось ждать, поскольку на момент их прибытия в систему хотя бы один обслуживающий канал был свободен. Средняя длина очереди без нулевых входов, очевидно, должна быть больше средней длины общеизвестной и более привычной математической очереди. При этом минимальная средняя длина такой очереди - это нуль, как и у обычной математической очереди, то есть в среднем в такой очереди, как и в математической очереди, может находиться любое число заявок.
Титовцев Антон Сергеевич, кандидат технических наук. E-mail: notna6683@mail.ru
Таким образом, если обычная математическая очередь рассчитывается как среднее по всем заявкам, побывавшим в системе, то очередь без нулевых входов должна рассчитываться как средняя величина за вычетом тех требований, которые были обслужены немедленно, поскольку попали в систему в тот момент, когда, по крайней мере, один из каналов обслуживания был свободен. В соответствии с терминологией, принятой в системе имитационного моделирования ОРББ, будем называть нулевым входом такой приход заявки в систему, при котором имеется, по крайней мере, один свободный канал обслуживания, и заявка в этом случае обслуживается немедленно. Назовём физической очередью очередь, рассчитанную без учёта нулевых входов. В случае занятости всех обслуживающих каналов вновь прибывшая в систему очередная заявка будет вынуждена ожидать обслуживания в накопителе, при этом минимальное количество заявок в накопителе равно нулю.
Как и привычная математическая очередь, рассчитываемая с учётом нулевых входов, физическая очередь также характеризуются моментами числа требований, ожидающих обслуживания, и соответственно, моментами времени ожидания обслуживания вновь прибывшей заявкой в соответствующей очереди. Идея расчёта этих моментов заключается в том, что в этом случае при расчёте параметров очереди следует переопределить условие нормировки с учётом только тех состояний СМО, при которых отсутствуют нулевые входы. Говоря другими словами, нужно учесть при расчётах только те состояния системы, в которых присутствует физическая очередь. Таким образом, алгоритм расчёта параметров физической очереди полностью идентичен алгоритму расчёта математической очереди за исключением вероятностей возможных состояний СМО, которые теперь должны рассчитываться исходя из нового условия нормировки, исключающего состояния, соответствующие нулевым входам.
В данной работе в рамках модели наиболее общей системы массового обслуживания с произвольным количеством источников и ограничений по длине очереди [5] получены общие математические выражения для первых и вторых моментов числа требований, ожидающих обслуживания в физической очереди. Также рассчитаны соответствующие функции распределения, и как следствие, первый и второй моменты времени ожидания начала обслуживания заявкой в очереди указанного типа.
Вероятностные характеристики и моменты числа требований в физической очереди. В работе авторов [5] была представлена наиболее общая математическая модель открытой
] ] 1=0 1=1
ограничения длины очереди (объёма накопителя) для заявок )-й компоненты;
й /г А
1=0 1=1 1=2 где — интенсивность потока заявок )-й компоненты; ц — интенсивность обслуживания заявок одним обслуживающим устройством;
/г к И д
*о ^=1 ра ъ =1л;- = А; к. =—
1=0 1=1 1=2 Ц
многоканальной СМО, имеющей т обслуживающих устройств одинаковой производительности с экспоненциально распределенным временем обслуживания. Входной поток требований в этом случае является суперпозицией произвольного числа компонент к, каждая из которых представляет собой пуассоновский поток заявок, обслуживаемых в порядке очереди. Для каждого типа требований, поступающих в систему из )-го источника, действует своё ограничение на длину
очереди 8 ,, при этом в0 < < в2 < ■ • ■ < вк . В настоящей работе мы также будем использовать обозначения, принятые в работе [5]:
где р] — приведенная интенсивность потока заявок )-й компоненты.
Для расчёта параметров физической очереди на базе разработанной и представленной в работе [5] универсальной математической модели необходимо переписать условие нормировки таким образом, чтобы при суммировании в этом условии учитывались лишь те состояния системы, которые соответствуют физической очереди. В результате имеем
Ер = 1
(1)
- условие нормировки для расчёта физической очереди.
Подставив соответствующие выражения для вероятностей стационарных состояний Р в представленное выше условие нормировки, и решая полученное уравнение относительно Р0, получим следующее выражение величины Р0, используемое для расчета физической очереди:
^0 =
Т)ш пт I, я-1 ( 7) \
от! от! ' -
1=0 V т У
Я
т - Я Е , Я = т
( / Я ^ \ 1 - ~
, V т У
V у у
Я Ф т
состояния системы Р0, то есть вероятность полного отсутствия заявок в системе, которая в этом расчёте не нужна, а некоторый коэффициент пропорциональности, который может быть и большим единицы. Физический смысл при этом имеют лишь величины вероятностей соответствующих состояний системы Ц, которые, конечно, для всех значений 1 > т остаются меньшими, чем единица. При этом формулы для остальных вероятностных характеристик, предназначенных для расчёта параметров физической очереди, будут выглядеть точно так же, как и при расчёте характеристик обычной математической очереди с заменой Р,, на Р0. Имеем
Р =
(/г., Л
7+1
т
и я ЛЕ>
П
8=0
V т
я:
т!
т + е < 1 < т + е 1,
0 < 1 < к -1
- вероятности всевозможных состояний СМО в стационарном режиме;
= П
я,
т!
(
т
1 -
Я.
т
Е N
т - Я, Е, Я, = т
, Я Ф т
- базовая вероятность для расчёта характеристик СМО;
Подчеркнём, что в данном случае Ро представляет собой уже не вероятность нулевого
Р =П
т+ек 1 1
\(я„ ¡с
8=1
— Р
I 0
т!
У
У
Е
т
X
- вероятность полной занятости системы;
1 к
1'ц- =— ^аЛРЫ
Л0 1=1
- вероятность ожидания обслуживания вновь прибывшей заявкой в физической очереди.
В этом случае относительная пропускная способность системы д = Р№, абсолютная пропускная способность, очевидно, А = АЛ
Моменты числа заявок, ожидающих обслуживания в физической очереди рассчитываются по той же схеме, что и для математической очереди, однако в соответствующие формулы подставляются вероятности с тильдой, полученные с учётом соответствующего условия нормировки. Подробный вывод формул для моментов длины очереди приводится в работе [5]. В результате имеем
- среднее число заявок, находящихся в физической очереди
1 =1
Л
т-Кг
Ф т
к г>
+у£ ^кр ■
¿-1 Гви
12=1
тЕ.
г Р +
Л
т + Л
т—ЛАт—Л.
'г- + 2 Е, ,
т — Л
Е +
2Л
_I
т - Л
- + 2е
Р ,Лфуп
т+ е, ' /
(Е-Щ^^ + е^^, Л=т
+е:Р .
п т+8\
Приведём также варианты этих зависимостей для нескольких частных случаев, представляющих определённый практический интерес. Для комбинированной модели массового обслуживания, представленной в работе [4], имеем
^0 =
Г
лГ
1+-
т\ V
Л1 [1—("
т — Л
Е Л = т
Ч 1Е1 р
Кт ) т — р
т
Л ф т
Р,=
Л
-Р„, т<1<т+Е,
т) т\
1—т—Е,
т ) т\
¿>т + Е,
В данной модели принята следующая система обозначений:
А А А
р0 =—; р =—; Р =—;
МММ
о Ло О Л1
Л0 = — = Ро +Рх + Р; Л = — = Рх + Р-М М
В этом случае первый и второй моменты числа заявок, ожидающих обслуживания в физической очереди имеют вид
г=2 т
- второй начальный момент числа требований в физической очереди
I =
Л
т-Л Е, +1
'—{Рвх-ЕА+Е), Л*т
рт, л=т
+р
12 = -
1
т т — р
Л
р„
—ЕЛЕ, +
т - Л 2т
т + Л т - Л
1-Р -
ГВ1
Р.
т — Л
(Ех +1)( 2Ех +1)
+Р
Е2
6 1
Л Фт
Р К =т
1 В1'1Н
2 Е,
т + р
Рт
т — р\ т — р) где р = Р, ~ вероятность полной загрузки
накопителя, р — приведенная интенсивность потока заявок, для которых не действуют ограничения на длину очереди.
Для частной модели СМО с полным набором накопителей, предложенной в работе [6], вероятности стационарных состояний системы при расчёте физической очереди принимают вид
Р0 =
щг ^ ¡^ ^ ^ Н
т\
т! ж=1}=1 т
1-т Р Т>т
^=1 т т!
»г +1 < г <т+к.
и тогда моменты числа заявок, ожидающих обслуживания в физической очереди
т+к 1-т Т> Т}Ш
,=т+1 г=1 т т!
Для модели Х. Такаджи [7] вероятности стационарных состояний системы для расчёта физической имеют вид:
при этом вероятности стационарных состояний
1=1
1=1
+
Р0 =
пт пт
К0 , К0
т\ т\
Я
1-1 Я
т - Я
, Я Ф т
Е, Я = т
Р = I ^
т
К
т\
Р0, т<г <т + Е1
и моменты числа заявок, ожидающих обслуживания в физической очереди
/ =
Я,
т-Я,
/2 = -
я,
т - Я1
т + Я,
т - Я1
3-Я -
1 в\
-Е\Е +
2т т - Я,
Р
Д Ф да
(Е1 +1)( 2Е1 +1)
6
Точно такая же система формул будет, очевидно, действительной и для однокомпонентной модели СМО с очередью конечной длины (по классификации Дж. Кендалла - модель М/М/т/Е). В этом случае в вышеуказанных формулах следует лишь положить Я0 = Я, = р, откуда в соответствии с [2, 3] имеем
/ =
Р
(
т-р1 Ех + \
Л
1-Е^-Р
Р = т
I 2 Лря + Л)
(т-р )2 т-Р1 V 1 т-Р1 У р
Е1(Е1 + 1)(2Е1 +1)
6 '
В этих формулах
Л
1-1А
(да -1)! (да - р1)
и,™-1
рфт.
р = т
рхФ т
(т-1)!
- вероятность ожидания заявкой обслуживания, то есть вероятность того, что поступающее требование найдет все каналы занятыми
Р =
Р
т\т т
р
Е О'
(т -1)!
Руфт
Р =т
- вероятность отказа, который заявка получает, если приходит в момент, когда заняты все т обслуживающих каналов и все Е1 мест в очереди.
Наконец, для частной модели, впервые изученной в работе Дж. Коэна [1], имеем следующие вероятности стационарных состояний системы для расчёта параметров физической очереди:
^0 =
Щ
т
т! т - р
РМ*-
т
Т)т
г>м.
т!
(2)
В этом случае первый и второй моменты числа заявок, ожидающих обслуживания в физической очереди
Р
I =-
-Г,,
/2 =
т - р
_р(т + р)~
, ч2 1
(т-р)
Для наиболее простой однокомпонентной модели с неограниченной очередью (модель М/М/т) в формулах (2) следует принять Я0 = р и в этом случае отсюда имеем
I =■
Р
т - р
р = р( т + р)
(т-р)2
ставить соответствующие вероятности с тильдой, полученные с учётом нового условия нормировки (1). Кроме того, с учётом этого обстоятельства следует переопределить относительную пропускную способность системы Ц . В данном случае для универсальной модели СМО, представленной работе [5], имеем следующие конечные формулы для моментов времени ожидания заявки в физической очереди:
^ 1 V \К1Рв,-ЕР т-Я,
Время ожидания начала обслуживания заявкой в физической очереди. Алгоритм расчёта моментов времени ожидания обслуживания заявкой в физической очереди вполне идентичен алгоритму расчёта моментов времени ожидания заявки в обычной математической очереди, приведенному в работе [5]. В нем вместо обычных вероятностных характеристик лишь следует под-
А
2 (Рп-ЕР
\ ВI I т-
)
/л(т -Я1)
1 +
т
Р_
п
(К }
V я У
[е(е +1)(е+2 )-
1
¿=1
)pm+e¡ r ^
m2 и
mn-( m — R )
— Si—i (S—i + l)(s—i + 2)],R =
m
Для комбинированной модели массового обслуживания, приведенной в работе [4], при этом имеем
~ -1 tw-j
R
m-R
Г, (j'¡;] Е\Ё'т+Е1 ) +
Р = I w Я
т 1 Fl
m 1
(
Е1+-
m
m — p
1_ A:
/u^m-R^
2 ,, ( E\PÉ+E, )
m — H ^ 1 '
m — R
'1
Ф m
F -1 ~ 1 -P i? = /w
3m и
K-
2 \E+
¿ rFLS ^
m V m — p
и(m—p)
Для частной модели СМО с полным набором накопителей [6], имеем
7 __1 т+Ь
ш А Ат/л Д } '
Для модели Х. Такаджи [7]:
1 = - Rl (Р -ЕР ) = -■
P = I
lw ~
A m-R 1
А [¿u(m-R) -RiR+l^lR
2
R1
m — R Ф m
( p; RPm+Ei ) "
E2-1 -
-i-PBvR=m
3mp
; I
Ajum
tn+R '
~ -LA
£•,+1
P Ф m
2À
Pi = m
- дисперсия времени ожидания заявки в физической очереди
[(£.+!)(£.+ 2) -2 _ 3/1, А °'=\ (£,+!)(£,+5)
-2
"iL — tW ■
р1Фт
12À
Pi = m
Наконец, для модели Дж. Коэна [1]:
2 PRIS
i P p _L
A FlS A '
A m-p A
t2 =
Äju(m — p)
Эти же формулы будут справедливыми и для классической модели М/М/т/Е, если в них
положить Л0 = Л = р, в этом случае имеем
Для наиболее простой однокомпонентной модели с неограниченной очередью (модель M/M/m) при условии R = Р отсюда, очевидно, следует
~ Р ' /г2 = ^
= Т7-w 2 1 \1'
X\m-р) 2 р [m-р)
Следует отметить, что результаты, полученные по этому алгоритму, полностью совпадают с результатами имитационного моделирования СМО вышеуказанных типов в системе имитационного моделирования GPSS World.
СПИСОК ЛИТЕРАТУРЫ:
1. Cohen, J. W. Certain delay problems for a full availability trunk group loaded by two sources // Communication News. 1956. Vol. 16. No. 3. P. 105-113.
2. Кирпичников, А.П. Прикладная теория массового обслуживания. - Казань: Изд-во Казанского ун-та, 2008. 118 с.
3. Кирпичников, А.П. Методы прикладной теории массового обслуживания. - Казань: Изд-во Казанского ун-та, 2011. 200 с.
4. Kirpichnikov, A.P. Open systems of multicomponent flows differentiated service / A.P. Kirpichnikov, A.S. Titovtsev // Ciencia e Técnica Vitivinícola. 2014. Vol. 29 No. 7. P. 108-122.
5. Kirpichnikov, A. Mathematical model of a queuing system with arbitrary quantity of sources and size-limited queue / A. Kirpichnikov, A Titovtsev // International Journal of Pure and Applied Mathematics. 2016. Vol. 106, No. 2. P. 649-661.
6. Kirpichnikov, A. Mathematical model of open queuing system with full set of memories / A. Kirpichnikov, A. Titovtsev // International Journal of Pure and Applied Mathematics. 2016. Vol. 107, No. 1. P. 139-143.
7. Takagi, H. Explicit delay distribution in first-come first-served M/M/m/K and M/M/m/K/n queues and mixed loss-delay system // International Journal of Pure and Applied Mathematics. 2007. Vol. 40, No. 2. P. 185-200.
PHYSICAL QUEUE IN THE APPLIED QUEUING THEORY
© 2016 A.S. Titovtsev Kazan National Research Technological University
The problem of various types of queues arising in queuing systems is considered in the paper. The concept of a physical queue is introduced; universal mathematical formulas of the first and second moments of main discrete and continuous random variables characterizing behavior of physical queues for various models of mixed type queuing systems are obtained.
Key words: physical queue, quality of service (QoS), queuing system
Anton Titovtsev, Candidate of Technical Sciences. E-mail: notna6683@mail.ru
x