CC BY
87
е— 1 —е
Вестник РУДН, Серия Математика. Информатика. Физика. К2 3-4. 2007. с. 45-52 45
Математическая теория телетрафика
УДК 519.6:621.39
Построение и анализ моделей системы с групповым обслуживанием заявок
К. Е. Самуйлов, Н. В. Першаков, И. А. Гудкова
Кафедра систем телекоммуникаций Российский университет дружбы народов Россия, 117198, Москва, ул. Миклухо-Маклая, 6
Статья посвящена исследованию системы с групповым обслуживанием со случайной длиной поступающей заявки и таймером, ограничивающим время формирования группы. Для простейшей модели типа М|М[к] излагается метод расчёта квантилей функции распределения времени ожидания начала обслуживания. В статье также представлены результаты имитационного моделирования для системы с таймером и примеры численного анализа.
Ключевые словА: система массового обслуживания, групповое обслуживание, таймер формирования группы заявок, время ожидания начала обслуживания.
Введение
Системы массового обслуживания (СМО) с групповым обслуживанием заявок находят применение в таких сферах жизнедеятельности человека, как промышленность, транспорт и связь. Начало исследованиям данного класса СМО в середине 1950-х гг. положили работы зарубежных учёных Н. Бейли и Ф. Даунто-на [1,2]. Бейли впервые рассмотрел систему с пуассоновским входящим потоком и переменным, ограниченным сверху размером группы заявок. Методом производящих функций он получил выражения для среднего значения и дисперсии числа заявок в очереди для цепи Маркова, вложенной по моментам, непосредственно предшествующим моментам окончания обслуживания групп заявок. В качестве дисциплины обслуживания им было рассмотрено распределение х2 (хи-квадрат) с чётным числом степеней свободы. Даунтон, в продолжение исследований Бейли, нашёл связь производящей функции с преобразованием Лапласа-Стилтьеса функции распределения (ФР) времени ожидания начала обслуживания, получил выражения для среднего значения и дисперсии, а также дал некоторые предварительные замечания к методу вычисления квантилей.
В настоящее время существует достаточное число работ (см. например, список источников в [3]), посвящённых исследованию СМО с групповым обслуживанием и предназначенных для решения различных прикладных задач. С целью унификации результатов этих исследований и их применения к анализу показателей качества современных протоколов передачи данных мы модифицируем и расширяем классификацию Башарина-Кендалла [4] для систем с групповым обслуживанием заявок. Анализ простейшей СМО с групповым обслуживанием проведён с использованием матрично-аналитического метода, а для численного анализа системы с таймером применён метод имитационного моделирования.
ф
ф
1. Система массового обслуживания с ограничением на время формирования группы
заявок
Рассмотрим СМО типа GIix |GI[X'T]|1| (R, L), состоящую из одного обслуживающего прибора и накопителя ограниченной ёмкости. На систему поступает рекуррентный стационарный поток заявок с интенсивностью Л, а длина 1\ поступающей заявки является дискретной случайной величиной (СВ) с распределением 7т = Р {1\ = тп}, m = 1, М. Обслуживание заявок производится группами, причём ёмкость накопителя рассчитана на R групп, и предполагается, что обслуживаемые прибором заявки сохраняют за собой место в накопителе и одновременно покидают систему в момент окончания обслуживания. Время формирования каждой группы заявок ограничено таймером 0 ^ T ^ то, и, поэтому, длина X группы является дискретной СВ, максимальное значение которой равно L. Обслуживание группы начинается либо в момент освобождения прибора, если группа сформирована, либо в момент завершения формирования группы, если прибор не занят обслуживанием, а время обслуживания группы заявок распределено по произвольному закону с ФР B(x). На рис. 1 изображена диаграмма, поясняющая процедуру формирования группы, которая считается сформированной либо в момент истечения таймера, либо в момент поступления заявки, длина которой превышает доступное ей место в группе, причём при условии занятости прибора, заявки могут поступать в группу и после истечения таймера.
Заметим, что в частных случаях рассматриваемая СМО была изучена в ряде работ отечественных и зарубежных авторов. Например, в [5] получены вероятностно-временные характеристики (ВВХ) СМО Mi|GK'°]| 1| (то, K) = M|G[K'°] |1|то. В [3] проведено исследование СМО Mi|GK>~l 111 (то, K) = M|G[k]|1|to, кодированной как M|G[K]|1|FS (англ. FS - Fixed Size), в случае, когда размер группы заявок является фиксированным.
Рис. 1. Диаграмма состояний группы заявок
В следующем разделе статьи для анализа простейшей системы типа М\М[к]|1|то применён матрично-аналитический метод [6], а в последнем разделе для численного анализа ВВХ СМО М\х\Д[Х'Т1\1 \{К,Ь) —метод имитационного моделирования [7].
2. Простейшая система с групповым обслуживанием заявок
Рассмотрим СМО типа M|M[к1|1|то, в которой обслуживание группы заявок
распределено экспоненциально с параметром р. Введём марковский процесс (МП)
X (t) как число заявок в очереди в момент времени t ^ 0, диаграмма переходов
состояний которого имеет вид, показанный на рис. 2. Целью анализа является
стационарное распределение pj = lim P {X (t) = j}, j ^ 0 и квантилей ФР
t—
W (ж) времени ожидания начала обслуживания.
ля ял я я
Рис. 2. Диаграмма переходов состояний СМО M|M
Предложение 1. Если K > 1 и р < K, то распределение {pj} МП X (t) вычисляется по формулам
pj
PK-1
j Ri
i=0
р +1-RK '
j = 0, K - 2,
K-2 i
iЕ Ей
г=0 fc=0 p + 1 -
k
+1
-i
1
р +1-RK : У PkR?-K,
, j = K - 1,
j = K, j>K,
(1)
где R
корень уравнения
rK+1 - (р +1) r + р = 0, |r| < 1.
(2)
Доказательство. Для доказательства воспользуемся матрично-аналитическим методом [6], который применим для неприводимых МП с инфинитезимальными матрицами вида
a
g0 g-i g-2 g-3 ■ ■ G-n
g1 fo f-i f-2 ■ ■ f-n+i f f - n 0
g2 fi f0 f-i ■ ■ f-n+2 f-n+i f F-n
g3 f2 fi fo ■ ■ F-n+3 F-n+2 f-n+i f-n
(3)
где с» и е», г ^ —п — прямоугольные матрицы. Нетрудно убедиться, что матрица интенсивностей МП X (¿) может быть приведена к виду (3) в случае п =1 путём разбиения её на перечисленные ниже блоки:
[go]ij = < л
- Л, i = 1, K, j = i i = 1,K -
i = l,K-l, j = i +1, иначе;
e— e
—e e
ф-
-е
г—1 с\К—%
С = (0 , м, 0
г = 1,К;
сг = 0, г > К, сг = К, Уг;
с-!
(ок-1,А), ё1ш с— = К;
Го = - (А + м); Гк = м;
Рг = 0, г = {1,...,К - 1,К + 1,...} ; Р-1 = А.
Вероятности pj, 3 > К можно выразить через вероятность рк, если МП X (¿) положительно-возвратный, а матрицы, соответствующие ведущим главным минорам матрицы, которая получается вычёркиванием первой строки и первого столбца блоков матрицы а, не являются единичными. Оба условия выполнены, и, следуя [6], получаем
Pj = рк ^
]-к
3 >К,
(4)
Е = Иш Фгд,
Фг.1 =
А + м
А + м (1 - Ф-м)
г = 1
, г > К.
(5)
Переходя к пределу в (5), получим уравнение м^к+1 — (А + м) г + А = 0 или (гх + гк_1 + ... + г2 + г — ^ (г — 1) = 0, одним из корней которого и является
исследуемая величина К. Уравнение гк -\-гк~1 _|___|_г2_|_г_Л = о согласно правилу знаков Декарта, имеет ровно один действительный положительный корень. Если Е ^ 1, то не выполнено условие нормировки для распределения стационарных вероятностей pj, и поэтому искомое значение Е С (0; 1) будет единственным действительным положительным корнем уравнения (2). Заметим, что значение
К £ (1+р)к) может быть найдено одним из итерационных методов.
Вероятности pj, 3 ^ К находятся из системы уравнений равновесия
2к-1
(ро,... ,рк-1) со + PjСj-к+l = о, j=к
(ро,... ,рк-1) с-1 - (А + м) Рк + МР2к = 0,
и условия нормировки, что совместно с уравнением (4) доказывает сформулированное предложение. □
г—
А
А
Для расчёта квантилей юв необходимо найти ФР Ш (ж). Эта функция определена ниже в предложении 2, которое мы приводим без доказательства.
Предложение 2. Пусть aj = К (bj + 1) - 3 - 1 — число свободных мест для заявок в формирующейся группе, bj = — число сформированных групп в системе при условии, что в момент поступления заявки в системе находилось 3
ф-
заявок. Если К > 1 и р < К, то ФР Ш (ж) времени ожидания начала обслуживания вычисляется по формулам
W (ж) = <
Е
j=0
Е
j=0
Е
j=0
Г \ai -1 Г nbj fbj-1
Pj I --—e~ dt I £-—e'^dt
L 0+
Pj
0+ x
Pj
0+
(aj - 1)!
цЪ] tbi ~1 \aj taj-1
0+
(bj - 1)!
, aj > 0, bj > 0, ж > 0,
e-^dt e-Aidt
aj = 0, bj > 0, ж > 0, aj > 0, bj = 0, ж > 0,
(6)
1, о, = 0, Ъ, = 0, ж > 0, ч 0, ж < 0.
Вычисление квантилей юп ФР Ш (ж) с использованием формул (6) проводится только численно. Для достижения заданной точности е необходимо найти значения ж 1 и ж2, 0 < Ш (жх) < и)в < Ш (ж2) < 1, такие что |Ш (ж2) — Ш (жх)| < е. Пример результатов вычислений ФР Ш (ж) по формуле (6) показан на рис. 3 для случая К = 5, р-1 = 1, 25 и р £ (0; К). Из графика видно, что для малых и предельных значений р ФР Ш (ж) ведёт себя как распределение с «тяжёлым хвостом», а в сечениях ФР Ш (ж) горизонтальными плоскостями получаем зави-
п
симости квантилей юи от значений величины нагрузки р.
На рис. 4 для наглядности представлена зависимость среднего времени ожидания начала обслуживания, вычисленная с использованием результатов работы [8]. Поведение графика достаточно очевидное, поскольку, чем меньше нагрузка, тем реже поступают заявки и требуется больше времени на формирование группы, и, следовательно, больше значение среднего времени ожидания. При увеличении нагрузки время ожидания убывает до некоторого минимального значения, а затем начинает возрастать, т. к. группы заявок накапливаются в очереди.
30 1.25
3,5 4 4,5 5 5
Рис. 3. ФР времени ожидания начала обслуживания и её 95% квантиль
| И__
S | 20
к S.
|| 15
0 ^
1 5 ю S1
5 0
1,5 2 2,5 3 3,5 Р
4,5
Рис. 4. Среднее время ожидания начала обслуживания
-е
3. Численный анализ модели с таймером
Рассмотрим СМО типа Mix |D[X'T] |1| (Л, L), для численного анализа ВВХ которой была разработана имитационная модель [7], а для вычисления квантили w0'95 использовался метод, предложенный в [9]. Модель была разработана для анализа показателей качества обслуживания протокола управления потоковой передачей (англ. SCTP — Stream Control Transmission Protocol) [10]. Процесс передачи
по протоколу состоит в следующем: поступающие сообщения инкапсулируются в порции данных (заявки), которые объединяются в пакеты (группы заявок). Пакеты передаются по каналу с пропускной способностью 2 Мбит/с, а длины порций данных заданы распределением, указанным в табл. 1, где единицей измерения длины заявки является один байт. Значение таймера Т было выбрано
Таблица 1
Распределение длин порций данных (байт)
1х 60 36 32
7т 0,23 0,41 0,36
равным 12 мс, что соответствует реальным параметрам оборудования компании Cisco Systems, поддерживающего протокол SCTP, а максимальная длина пакета L = 1468 байт, исходя из ограничений на максимальную длину кадра Ethernet.
Как видно из рис. 5, график среднего времени пребывания заявки в системе (задержки передачи сообщения в канал) имеет три экстремальные точки, определяющих четыре промежутка монотонности. Для их объяснения воспользуемся зависимостью среднего времени ожидания начала обслуживания (задержки в очереди на передачу) от нагрузки. При сравнительно низкой нагрузке (на рисунке эта область значений обозначена римской цифрой I) интенсивность поступления заявок мала, время формирования группы определяет таймер и, следовательно, величина задержки близка к значению таймера. С ростом нагрузки время задержки быстро убывает до того момента, когда среднее число заявок в группе, находящейся на обслуживании, становится более одной. При последующем увеличении интенсивности поступления заявок (область II) время формирования группы по-прежнему превышает среднее время обслуживания заявки, которое ещё не достигло максимального значения. Этим обусловливается дальнейшее убывание величины задержки, что в сравнении с предыдущим интервалом происходит медленнее, поскольку число заявок в группе увеличивается, приводя к росту времени их обслуживания. После того как группа достигает максимальной длины L = 1468 (область III), чем больше становится нагрузка, тем меньше на формирование группы влияет таймер, но, как и на только что рассмотренном интервале, время формирования группы превышает время обслуживания, что приводит к дальнейшему более быстрому уменьшению задержки. И, наконец, когда среднее время обслуживания превысит среднее время формирования группы (область IV), заявки начинают накапливаться в накопителе, способствуя увеличению среднего времени ожидания начала обслуживания.
Среднее время пребывания заявки в системе
Среднее время ожидания начала обслуживания
Среднее время обслуживания
Среднее время формирования группы
Среднее число заявок в группе на приборе
Рис. 5. Средние значения вероятностных характеристик СМО Mix |D[X'12] |1| (ж, 1468}
ф
ф
Современные прикладные задачи, и в частности, задачи, связанные с анализом показателей качества функционирования протоколов сетей связи, требуют уточнения имеющихся и разработки новых методов исследований. Так, оказывается, что процесс передачи данных по протоколу ЯСТР — новому решению, предназначенному для обеспечения работы приложений с высокими требованиями к показателям качества обслуживания - может быть описан в виде СМО с ограничением на время формирования группы заявок.
В статье представлена часть результатов исследований системы с групповым обслуживанием с учётом случайной длины заявки и таймера, ограничивающего время формирования группы. Для простейшей модели получен метод расчёта квантилей ФР времени ожидания начала обслуживания и приведены результаты имитационного моделирования для системы с таймером.
В заключении авторы выражают искреннюю признательность профессору кафедры систем телекоммуникаций РУДН Наумову В. А. за помощь и ценные советы в процессе исследований.
1. Bailey N. On Queueing Processes with Bulk Service // Journal of the Royal Statistical Society, Ser. B. - Vol. 16, No 1. - 1954. - Pp. 80-87.
2. Downton F. Waiting Times in Bulk Service Queues // Journal of the Royal Statistical Society, Ser. B. - Vol. 17, No 2. - 1955. - Pp. 256-261.
3. Tadj L. Explicit Solution of a Quorum Queueing System // Stochastic Analysis and Applications. - Vol. 21. - 2003. - Pp. 703-717.
4. Башарин Г. П., Харкевич А. Д., Шнепс М. А. Массовое обслуживание в телефонии. - М.: Наука, 1968.
5. Саати Т. Л. Элементы теории массового обслуживания и ее приложения. — М.: Советское Радио, 1971.
6. Naoumov V. A. Matrix-Analytic Solution for a Class of Markov Processes with Toepliz-like Generator Matrices - Advances in Matrix-Analytic Methods for Stochastic Models // Notable Publications. - New Jersey, 1998.
7. Гудкова И. А., Першаков Н. В., Самуйлов К. Е. Имитационная модель для анализа показателей качества функционирования протокола SCTP // Труды 62-ой Научной сессии, посвященной Дню радио. — Москва: РНТОРЭС, 2007. - С. 231-233.
8. Самуйлов К. Е., Чукарин А. В., Першаков Н. В. Разработка модели функционирования протокола управления потоковой передачей // Вестник РУДН, серия «Прикладная и компьютерная математика». — Т. 4, № 1. — 2005. —
9. Jain R., Chlamtac I. The P2 Algorithm for Dynamic Calculation of Quantiles and Histograms without Storing Observations // Communications of the ACM. — Vol. 28. — 1985. — Pp. 1076-1085.
10. Гайдамака Ю. В., Першаков Н. В., Чукарин А. В. Модель протокола SCTP и ее применение к анализу характеристик сигнального трафика в сетях сотовой подвижной связи // Электросвязь. — № 8. — 2007.
Заключение
Литература
С. 40-47.
ф
ф
е-
-е
UDC 519.6:621.39
Batch Service Queue Model Design and Analysis
K. E. Samouylov, N. V. Pershakov, I. A. Gudkova
Telecommunications Systems Department Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, 117198, Russia
In this paper we introduce an extended Basharin-Kendall queuing notation with queuing system Ghx\GI[X'T]|1| {R, L) in which jobs are served in batches, each job having a random length and the batch construction time limited to a maximum. The simplest fixed-size batch queue model M\M[K]\1\ ra is studied. Expressions are obtained for the distribution function of the waiting time and a brief discussion is included of its quantiles computation. We also develop a simulation model and present an example of its application to estimate quality of service in telecommunication networks.
e— e
—e e
