В.Н. Тарасов, Н.Ф. Бахарева
НЕПРЕРЫВНЫЕ ДИФФУЗИОННЫЕ МОДЕЛИ МАССОВОГО ОБСЛУЖИВАНИЯ И МЕТОДИКА РАСЧЕТА ИХ ХАРАКТЕРИСТИК
В статье предлагается методика расчета характеристик стохастических сетей на основе новых диффузионных моделей массового обслуживания. Рассмотрены системы массового обслуживания с бесконечной очередью, с конечной очередью и потерями, с переменными параметрами поступления и обслуживания заявок.
Для расчета характеристик систем массового обслуживания (СМО) в предположении, что времена между заявками во входном потоке и обслуживания распределены по произвольному закону, не существует точных методов. Использование же традиционных методов, основанных на пуассоновских входных потоках и экспонен-циальности времени обслуживания, не всегда оправдано. Существующие аппроксимацион-ные методы, использующие два первых момента распределений входного потока и времени обслуживания, имеют большую погрешность в случае отклонения коэффициентов вариаций этих распределений от единицы.
Для расчета характеристик сетевых моделей систем сеть декомпозируется на отдельные узлы (рисунок 1).
В этой статье излагается новый метод декомпозиции разомкнутых однородных сетевых моделей систем управления на отдельные узлы на уровне двух моментов распределений параметров потоков. С целью определения параметров выходного потока, а также основных харак-
теристик отдельного узла вводится двумерная диффузионная аппроксимация процессов функционирования узла. Указанная методика также обобщается на сетевые модели с неоднородными потоками заявок.
Для определения среднего значения и дисперсии распределения времени между соседними заявками на входе и выходе каждого узла сети необходимо выразить среднее и дисперсию времени в выходном потоке отдельного узла через параметры входного потока и времени обслуживания. Для решения этой задачи выведены точные формулы:
?вых Рот°, (!)
Ошх = р'0 0 + р0(1 - р0ХП)2, (2)
гдетвыхБвых, Бц - соответственно средние и дисперсии времени между заявками выходного потока и времени обслуживания, т'х, Б'х - среднее и дисперсия остаточного времени (времени простоя узла), а рО - вероятность того, что обслуженная заявка оставляет узел пустым. Для системы без потерь средние времена между за-
Выход из сети
Рисунок 1. Структура 1-й СМО сети (А - узел композиции, В - узел декомпозиции)
явками в выходном и входном потоках совпадают. Уравнения (2) решаются совместно со следующими уравнениями преобразования дисперсии времени между заявками:
Бъ=(11/ Ах) А +(12/ Ах) А (3)
- для суммы двух независимых потоков с интен-сивностямиЯ1;Я2 и дисперсиями Э Б2 (узел композиции потоков) и
Бр = -Р
Б +
1 - Р
(4)
- для потока с интенсивностью 1 и дисперсией времени между заявками Б, в котором заявки с вероятностью р уходят из потока (узел декомпозиции потоков). Из уравнений (2)-(4) следуют основные уравнения декомпозиции сети, так называемые уравнения баланса для средних и дисперсий времени между заявками на входе и выходе /-го узла сети (7=1,...,«)
1 =т:
+ Х Рл1! 1=1
Авхл
1
\3
А
+ 1
1=1
Рр1! 1
3
1
Рр
Ав
1 - Рд Р12
(5)
(6)
где тш иБ0, - среднее и дисперсия времени между заявками в потоке, идущем от внешнего источника на вход 7-го узла, 1 - интенсивность
потока на входе и выходе 7-го узла, Авхг-, Авыхг- -дисперсии времени между заявками во входном и выходном потоках 7-го узла, Рд, - вероятность передачи заявки от у'-го узла к 7-му, а п - количество узлов в сети. Для решения системы (6) совместно с (5) и (2) предложена итерационная процедура, где в качестве первого приближения в уравнении (2) используется замена Р0 = Р0,Т'Х=ТХ, Бх = Бх . Здесь Р0 =1 -Р -вероятность простоя узла,Тх иБх - среднее и дисперсия времени между заявками во входном потоке. Это приближение в случае экспоненциальной сети не вносит погрешности, а уравнения (6) становятся линейными относительно искомых дисперсий. На последующих этапах значения дисперсий входных и выходных потоков уточняются с использованием уравнений (2) и (6), одновременно вычисляются основные характеристики отдельных узлов и всей сети в целом.
При приближенном расчете замкнутых сетевых моделей систем в выражениях (5) и (6) первые слагаемые исчезают. В этом случае в итерационной процедуре значения интенсивностей 1 подбираются так, чтобы с заданной точностью е выполнялось условие
^М, = N,
1=1
2
Рисунок 2. Аппроксимация СМО 01/0/1/ш с потерями, а также с переменными параметрами двумерным диффузионным процессом (хр х2) (а) и выделенная область (б) (индекс к определяет уровень по оси хр а 1 - по х2). Траектория 1 - для СМО 01 / О /1/ ¥ при т , траектория 2 - для СМО с потерями, траектория 3 - для СМО с переменными
параметрами.
где через N обозначено среднее количество заявок в 7-м узле, а N - количество заявок, постоянно циркулирующих в замкнутой сети.
Для определения параметров р0 ,Т'Х,, входящих в формулы (1) и (2), а также основных характеристик узла процесс его функционирования на периоде занятости аппроксимируется двумерным диффузионным процессом (х1, х2). Процесс х^) аппроксимирует число заявок N1(t), поступивших в узел к моменту времени t, а х() - число заявок N2(t), покинувших узел к тому же времени. Текущее значение N числа заявок в узле определяется разностью N = [х ]-[х2 ], где [х]-целая часть от х (рисунок 2а, траектория 1).
Потребуем, чтобы компоненты двумерного процесса (х1, х2) в моменты времени первого прохождения целочисленного уровня (моменты поступления и ухода заявок) имели средние значения и дисперсии, совпадающие со средними и дисперсиями компонент дискретного процесса N N2). Тогда коэффициенты сносаа1 =т-1 и диффузии в, = Цт"3 процесса х. можно выразить через средние значениям и дисперсии Б. интервалов времени между скачками дискретного процесса N.(7=1,2). В области Б, определенной условием N > 0 и приведенной на рисунке 2а, плотность распределения ю((, х1, х2) векторного диффузионного процесса удовлетворяет уравнению Колмогорова
Эю = Д, ^ £
Ь, Э2ю
2 Эх.
Эю
„ - а, 2 Эх
(7)
дельного узла. Указанные распределения можно представить через решение уравнения (7), т. е. плотность распределения! ((, х1, Х2 ). Для этого в статье с использованием аппарата теории функций Грина(/, Хр х2 | к, у2) получены рекуррентные формулы:
Так как период занятости начинается с уровня х1=1, то начальным условием для (7) бу-детю(о,х1,х2)=5( -1)5(х2). Рассматривая функционирование узла на периоде занятости, добавим к уравнению (7) граничное условие поглощения®^ = 0. Граница Г, определенная условием N ] = 0, имеет ступенчатый характер, а попадание ординаты процесса на границу физически означает завершение периода занятости. Вследствие сложного характера границы, решение уравнения (7) будем искать в виде совокупности решений в подобластях ВК = (х1 < к +1, х2 < к), (к = 1,2,...) и «сшивая» их на границахх1 = к, (к = 2,3,...). Введем обозначения ¥к (у1) для распределения ординаты процесса х1 в момент достижения уровня х2=к (границы Г) ирк (у2) для распределения ординаты процесса х2 в момент прохождения процессом (х1, х2) уровня х1=к+1 (границы области Бк). Эти распределения позволяют определить все основные характеристики функционирования от-
юк ((, х1, х2) = 1 <Рк -1 к Х1,х2 1 k,)йУ2, (8)
Рк(У2)=
¥
Р- 1(у2Яр(У2 юлур^=
о
= вр(У2 I 0))
(9)
бр(У2|У2)=
1 а а
ехрЬ1 + У2- У2 +1)] X
ж^ЬЬ
Ь1 Ь2
х ^ К (2л/Ь^) - ^ К )];
1 . (у;-у2 +1)2 „ 1 . (у2 + у2 +1)2
1 2Ь1 2Ь2
; Ь = Ж +
2Ь
22 а1 + а2 .
г=—+—; у2 е[0, ¥]; 2Ь1 2Ь2 2
К1(.) - функция Макдональда;
У(У1> =
¥
= Р -1( у2% (Уl, у2)лУ2 (У1( =
о г
= (у1|0)), (10)
1 + у о а а„
(у1|у2) = • ехр[-^(1- у1) + -2(1 + у2 )]х
Ьз =
х Ц ь К1 (^^^З7^) К1 )];
(1 - У1)2 + (1+у2 )2; д = (^+(1+У2 )2.
2Ь
2Ь
4 '
2Ь
21 У2 е [0, ¥].
2Ь
Если ту и - среднее и дисперсия распределения у(у1 )=£ук (у1), то искомые параметры
к=1
двумерного диффузионного приближения р0 Хх, В, входящие в уравнение (2) и формулы для основных характеристик узла, можно выразить через эти величиныт¥, Ву и извест-
2
2
1
2
ные параметры входного потока: среднееf и дисперсию Di времени между заявками. Решение этой задачи связано с использованием функции плотности времени первого достижения марковским процессом заданного уровня g(t|y^). Окончательные формулы имеют следующий вид:
fl=fl- my, (11)
Di = Dxmv + f1 Dy
(12)
узла и
P0 = 1/E kPk PK = JVk (yi )dy
(13)
м . , Хоб = Eim)
(14)
- для интенсивности поступления потока обобщенных заявок на вход 7 - го узла (7=1,.п);
f$ = e
м (i(m) / N
(m)
m=1
м
V
О об
m (m) = 1/ ( )
Vi /f(m)
/ m
\ ' ^ j
(15)
m=1
Dm=e d^i+((t)
/ io
i06 = Et)
m=1
DOб = E (( / 1 )
(16)
об 0i
m=1
м
p об = e p (m ]if) /1
(17)
для среднего и дисперсии времени простоя
- для вероятности того, что обслуженная заявка оставляет узел пустым. Подставив последние формулы в (1) и (2) получим окончательные формулы для определения параметров выходного потока.
В случае сетевых моделей с неоднородными потоками, в которых заявки различаются маршрутами и законами поступления и обслуживания, параметры потоков различных классов (т=1,...,М) усредняются с целью приведения неоднородного потока к однородному. Эти параметры будут описывать так называемую обобщенную заявку. При этом соблюдается условие, чтобы однородный поток создавал такую же нагрузку в каждом узле сети, что и неоднородный поток. Для определения параметров потоков обобщенных заявок получены следующие формулы:
- для значений обобщенной матрицы вероятностей передач.
Итерационная процедура расчета сети с неоднородными потоками будет такая же, что и для однородного потока с приведенными параметрами, т. е. заключается в решении уравнений (5) и (6) совместно с (2). Используя полученные характеристики для отдельного узла сети, определяются характеристики сети для каждого типа заявок, аналогично сетям с однородными потоками.
Далее рассмотрим модели систем с ограниченной емкостью накопителя и с потерями (рисунок 2а траектория 2), а также с переменными параметрами законов поступления и обслуживания заявок (рисунок 2а траектория 3). Такие модели отражают процесс прохождения задач в отдельной ЭВМ и в более сложных комплексах - сетях передачи данных, имеющих буферные накопители ограниченной емкости. Последняя модель также отражает работу саморегулирующихся динамических систем.
Для определения характеристик моделей систем, содержащих узлы с потерями и ограничениями на длину очереди, сеть также декомпозируется на отдельные узлы на уровне двух моментов. При этом уравнения баланса (5) и (6) модифицируются: учитываются потоки отказов, идущие от узлов с потерями к другим узлам. Вводится матрицаQ = } вероятностей передач заявок потоков отказов от 7-го узла на вход у-го (7,у'=1,...,п), а в системы (5) и (6) добавляются соответствующие слагаемые, учитывающие потоки отказов. Тогда уравнения баланса примут вид:
Ibx fo
e РЛвых + e qji1
Ji JOTK
(18)
/ __, \3
- для среднего и дисперсии времени обслуживания;
Dbx =
Л,,
Doi+E j=1
j=1
i p
7ВЫ11 Ji
J=1
\3
1b
Pni
1 - P»
D вых + - Ji
1 P
JB^IX1 Ji
+e
J=1
/1 q
л]'откЧ Ji
1'
1b
DJrTK +
1 - If.
]°TK 1 q ^JqtkH Ji
(19)
- для интенсивности и дисперсии потока от внешнего источника;
Для определения среднегоТотк и дисперсии БТотк времени между заявками в потоке отказов получены формулы:
^отк = ?вх?вых /(вых — *вх )' (20)
0
1
+
q
гдетвх, итвых - средние времена между последующими заявками во входном и выходном потоках, и
^тотк отк^отк / Тц ,
(21)
где Тц - среднее время цикла занятости, отк -дисперсия числа потерянных заявок на цикле занятости. Численное определениеотк выполняется аналогично вычислению распределения^ (У2): вычисляются условные распределения числа и квадрата числа потерянных заявок.
Основные характеристики для узла с потерями определяются с учетом того, что уравнение Колмогорова решается в области Б, определенной двухсторонними условиями N > 0 и N < т, т - максимально допустимое число заявок в узле (рисунок 2а траектория 2). Граничные условия для уравнения (7) физически эквивалентны условиям поглощения на границе Г1, гдею|Г1 = 0 и отражения на границе Г2 где gradю|T2 = 0. '
В задаче анализа систем с переменными параметрами законов поступления и обслуживания, зависящими от состояния системы (длины очереди к ресурсу), вышеуказанные распределения ук (у1) и (рк (у2) определяются отдельно для каждой из подобластей0'к(к = 1,...,т;г = 1,2,...). ОбластиОгк физически соответствует период функционирования узла, на котором в узле находится ровно к заявок (индекс к определяет уровень по оси х 1 - уровень по х2). Этот период характеризуется своими интенсивностями поступления и обслуживания заявок, а также дисперсиями этих времен, зависящими от состояния системы. Такой подход означает решение уравнения Колмогорова с переменными коэффициентами сноса и диффузии, зависящими от длины очереди в узле (рисунок 2б).
Приведены результаты машинных экспериментов по оценке точности метода двумерного диффузионного приближения. Точность методов исследована для диапазонов изменения
а) С = 1
W
20
15
10
1 = 0,5
1 = 0,1
12 3 4 5
б) С = 1
W А
20
15
10
1 = 0,9
1 = 0,5
1 = 0,1
1 2 3 4 5 С„
5
5
С
л
Рисунок 3. Зависимость времени задержки в узле при различных значениях интенсивности 1 входного трафика (время обслуживания нормированное): а) - от коэффициента вариации входного потока; б) - от коэффициента вариации времени обслуживания.
коэффициента загрузки от 0,1 до 0,9 и коэффициентов вариаций распределений длин интервалов между заявками во входом потоке и времени обслуживания от 0 до 5. Значения среднего числа заявок в узле сравнивались с результатами имитационного моделирования.
В качестве одного из параметров моделирования задавалось количество циклов занятости, которое в зависимости от загрузки изменялось от 1000 до 20 000. Для отдельных узлов сети результаты двумерного диффузионного приближения совпадают с результатами имитационного моделирования в пределах 2-5% и в некоторых крайних случаях отклоняются до 10%. При расчете сетевых моделей добавляются ошибки декомпозиции, поэтому погрешность получалась несколько выше - до 15%.
Эмпирические распределения и гистограммы реальных потоков, полученные в результате измерений вычислительной нагрузки и потоков, циркулирующих в системе, также свидетельствуют об отличии потоков от пуассоновс-ких, в особенности в системах с коммутацией
пакетов, в которых в режиме насыщения доминируют регулярные потоки. Исследованиями установлена нестационарность характеристик информационных потоков и неприемлемость аппроксимации гистограмм и эмпирических распределений экспоненциальным распределением. На основные показатели производительности систем, такие как время задержки сообщения, вероятность переполнения буферной памяти, время занятости канала и других существенное влияние оказывает вид распределения времени между соседними заявками во входном потоке и времени обслуживания. Эти факты отражены на рисунке 3.
Эксперимент также подтвердил, что моменты более высокого порядка на характеристики системы оказывают несущественное влияние. Сказанное подтверждает правомерность использования методов двумерной диффузионной аппроксимации как методов двухмоментной аппроксимации для анализа сетевых моделей сложных систем управления.
Список использованной литературы:
1. Кругликов В.К., Тарасов В.Н. Анализ и расчет сетей массового обслуживания методом двумерной диффузионной аппроксимации. - Изв. АН СССР - Автоматика и телемеханика, 1983. №8. - 74-84с.
2. Кругликов В.К., Тарасов В.Н. Приближенный метод декомпозиции разомкнутой стохастической сети. - Изв. АН СССР - Техническая кибернетика, 1983. №6. - 142-147 с.