Марковские сети массового обслуживания

До сих пор рассматривались марковские системы, в которых каждое требование проходило одну операцию обслуживания. Такие системы можно называть однофазными. Далее рассматриваются системы с многофазным обслуживанием, в которых требование получает обслуживание более чем о водном приборе (узле). Таким образом, можно говорить о сети узлов, каждый из которых представляет собой СМО (некоторые узлы могут иметь несколько обслуживающих приборов) с накопителем для образования очереди. Требования поступают в систему в различных точках, ждут в очередях обслуживания и, покинув один узел, поступают в другой для дальнейшего обслуживания.

При исследовании сетей возникает много новых аспектов. Например, важной становится топологическая структура сети, так как она определяет возможные переходы между узлами. Требуется также каким-нибудь способом описать пути отдельных требований. Большое значение имеет описание природы вероятностных потоков с помощью основных вероятностных процессов; например, в случае последовательной цепочки СМО, при которой требование, покидающее i-й узел, сразу поступает на обслуживание в -й узел; очевидно, что промежутки времени между последовательными уходами требований из предыдущего узла равны промежуткам времени между последовательными поступлениями в следующий узел.

Рассмотрим простейшую последовательную систему с двумя узлами, показанную на рис. 1.

Рис. 1. Система с двумя последовательными узлами

Каждый овал на этом рисунке изображает СМО, состоящую из очереди и обслуживающего прибора; внутри каждого овала указан номер узла (важно не путать такие схемы физических сетей с абстрактными диаграммами интенсивностей переходов, которые были использованы ранее). Предположим, что входящий поток является пуассоновским с интенсивностью , причем каждое требование поступает сначала в узел 1; предположим также, что этот узел содержит единственный обслуживающий прибор, время обслуживания которого распределено по показательному закону с интенсивностью . Таким образом, узел 1 представляет собой в точности СМО типа M/M/1. Предположим, далее, что узел 2 также состоит из единственного обслуживающего прибора с показательным временем обслуживания с интенсивностью . Основная задача состоит в вычислении распределения промежутков времени между последовательными требованиями, поступающими в узел 2; это естественно, эквивалентно задаче вычисления распределения промежутков времени между последовательными требованиями, уходящими из узла 1.

Пусть означает плотность распределения вероятностей промежутков между последовательными требованиями на выходе узла 1 и, как обычно, пусть означает ее преобразование Лапласа. Перейдем к вычислению в момент, когда требование покидает узел 1. Возможно одно из двух событий: либо в очереди имеется второе требование, готовое немедленно поступить на обслуживание в узел 1, либо требования нет (накопитель пуст). В первом случае промежуток времени, через которое это следующее требование покинет узел 1, распределен точно также, как и время обслуживания, и в этом случае получаем

.

С другой стороны, если при уходе рассматриваемого первого требования узел оказывается пустым, то приходится ожидать в течение двух промежутков времени: первый промежуток – время до поступления следующего требования и второй – время обслуживания этого требования. Так как эти два промежутка распределены независимо, то плотность распределения вероятностей их суммы равна свертке плотностей распределения суммируемых величин. Соответственно преобразование Лапласа плотности распределения суммы равно произведению преобразований исходных плотностей распределения, и следовательно,

,

где для преобразования Лапласа плотности распределения промежутков между поступающими требованиями уже известно явное выражение. Так как время обслуживания является показательно распределенной случайной величиной, то можно записать ; кроме того, как было показано ранее, вероятность того, что требование покинет систему пустой, а именно, равна . Это позволяет записать следующее безусловное преобразование Лапласа для плотности распределения промежутков времени между уходящими требованиями:

.

Используя проделанные ранее вычисления, получаем

.

Простые алгебраические преобразования дают

(1)

и, следовательно, распределение промежутков времени между уходящими требованиями

.

В самом деле =

Делая подстановку , получаем

.

Таким образом, приходим к замечательному выводу о том, что промежутки времени между уходящими требованиями, так же как и промежутки времени между поступающими требованиями, распределены по показательному закону с тем же самым параметром. Иначе говоря (с случае стационарной СМО), входящий поток, протекая через обслуживающий прибор с показательным распределением времени обслуживания, порождает выходящий пуассоновский поток. Этот основополагающий результат обычно называют теоремой Берке.

Фактически теорема Берке говорит больше, а именно, что исходящий поток стационарной СМО типа M/M/m с пуассоновским входящим потоком с параметром и показательным распределением времени обслуживания с параметром в каждом из m приборов является пуассоновским с тем же параметром .

Доказано также, что система M/M/m является единственной системой с обслуживанием в порядке поступления, обладающей таким свойством. Возвращаясь к рис. 1, видим, что в узел 2 поступает независимый пуассоновский поток, и, следовательно, этот узел также представляет собой систему M/M/1, что позволяет рассматривать его независимо от узла 1. Тем самым теорема Берке говорит о том, можно соединить последовательно¹ несколько узлов, состоящих из многолинейных СМО (с показательным распределением времени обслуживания для каждого прибора), и при этом будет сохраняться описанное свойство разложения на отдельные узлы.

Этот вопрос для произвольной сети массового обслуживания исследовал Джексон. Он рассматривал сеть, содержащую N узлов, причем каждый i-й узел состоит из обслуживающих приборов с показательным временем обслуживания с параметром ; в каждый i-й узел извне поступает пуассоновский поток требований с интенсивностью . Таким образом, при получаем обычную систему M/M/m. Покидая i-й узел, требование с вероятностью поступает в j-й узел, причем в формулировке задачи допускается, что . С другой стороны, вероятность того, что после обслуживания в i-м узле требование покинет сеть (и никогда не вернется обратно), равна . Необходимо вычислить полную интенсивность потока требований в заданный узел. Для этого нужно просуммировать (пуассоновские) потоки, поступающие извне, и потоки требований (необязательно пуассоновские), поступающие от других узлов сети. Обозначая через полную интенсивность потока, входящего в i-й узел, легко показать, что это множество параметров должно удовлетворять следующей системе уравнений:

. (2)

Для того, чтобы все узлы сети описывались эргодическими цепями Маркова, для всех i должно выполняться требование . Нужно снова предостеречь читателя от путаницы между понятием узла, используемом в данном обсуждении, и понятием состояния системы, которое ранее изображалось узлом графа. Самое удивительное состоит в том, что Джексон доказал, что каждый узел (скажем i-й) ведет себя в сети так, как если бы он был независимой системой M/M/m с входящим пуассоновским потоком с параметром . В общем случае полный входящий поток не является пуассоновским. Состояние рассматриваемой системы с N узлом описывается вектором , где означает число требований в i-м узле (включая обслуживаемые требования). Обозначим через стационарную вероятность этого состояния. Аналогично, через обозначим маргинальное (частное) распределение вероятностей того, что в состоянии равновесия в i-м узле будет находиться требований. Джексону удалось доказать, что совместное распределение по всем узлам разлагается в произведение маргинальных распределений, т. е. что

(3)

и что представляют собой стационарные вероятности для классической системы M/M/m. Этот последний результат обычно называют теоремой Джексона.

Модификация сети Джексона была исследована Гордоном и Ньюэллом. Исследованная ими модификация представляет собой замкнутую марковскую сеть в том смысле, что в системе рассматривается конечное и фиксированное число K требований, которые задерживаются в сети так, что ни одно из них не может покинуть сеть и ни одно другое требование не может поступить в нее; это соответствует теореме Джексона, в которой и для всех i (интересный пример систем этого класса, известных под названием циклических, был рассмотрен Кенигсбергом; в циклической сети узлы последовательно связаны друг с другом, причем последний связан с первым). В общем случае, рассмотренном Гордоном и Ньюэллом, не приходится надеяться на решение в виде произведения, так как между компонентами вектора состояний имеется следующая зависимость:

.

Так как рассматриваемая модель является частным случаем модели Джексона, то предполагается, что рассматриваемый марковский процесс с дискретным множеством состояний неприводим, и, следовательно, существует единственное распределение вероятностей . Однако в этой модели множество состояний систем конечно: в частности, нетрудно видеть, что число различных состояний системы равно числу способов размещения K требований по N узлам.

1 При этом не допускаются переходы, ведущие назад, так как это может нарушить пуассоновскую природу потока, выходящего после обратной связи. Последующее обсуждение работ Джексона показывает, что, несмотря на это, даже в сети с обратной связью отдельные узлы ведут себя так, как если бы на их вход поступал чистый пуассоновский поток, хотя на самом деле это не соответствует действительности.

72