Эргодические циклические ресурсные сети. Ii. Большие ресурсы Текст научной статьи по специальности «Математика»

УДК 519.177 ББК 22.18

ЭРГОДИЧЕСКИЕ ЦИКЛИЧЕСКИЕ РЕСУРСНЫЕ СЕТИ. II. БОЛЬШИЕ РЕСУРСЫ

Жилякова Л. Ю.1

(ФГБУН Институт проблем управления им. В.А. Трапезникова РАН, Москва)

В работе исследуется функционирование эргодических нерегулярных ресурсных сетей при больших ресурсах. Найдена формула для порогового значения ресурса Т. Показано, что при больших ресурсах циклические сети всегда устойчивы, колебаний в них не происходит. Найдена формула компонент вектора предельного состояния, соответствующих всем неаттрактивным вершинам. Сформулировано условие аттрактивности вершин в циклических сетях.

Ключевые слова: ресурсная сеть, стохастическая матрица, эргодическая цепь, аттракторы, предельное состояние.

1. Введение

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

1 Людмила Юрьевна Жилякова, кандидат физико-математических наук (zhilyakova. ludmila@gmail. com)

остаток на сколь угодно больших тактах t, ресурс в сети превосходит пороговое значение. При ресурсах, меньших порога, сеть сходна с моделями рассеяния на графах [10, 11]; при ресурсах, больших порога, некоторые вершины изменяют правило функционирования, - для них имеется сходство с целочисленными пороговыми моделями, описанными в [8, 9, 12].

Настоящая работа исследует функционирование эргодиче-ских циклических сетей при больших ресурсах. Она является продолжением статьи «Эргодические циклические ресурсные сети. I. Колебания и равновесные состояния при малых ресурсах» [4] При малых ресурсах эргодическая нерегулярная сеть описывается циклической цепью Маркова [5, 7]. При произвольном начальном состоянии в сети возникают периодические колебания между d предельными векторами, названные предельным циклом, где d - количество циклических классов в сети. В том случае, когда все d векторов совпадают, в сети при малых ресурсах возникает равновесие.

В указанной статье были получены следующие результаты.

1. Найдены пределы d сходящихся последовательностей степеней стохастической матрицы сети. Эти предельные матрицы состоят из циклических перестановок d предельных векторов при единичном суммарном ресурсе, находящемся в начальном состоянии целиком в одном из циклических классов.

2. Доказано, что все предельные векторы являются собственными векторами одной и той же предельной матрицы: предела степеней стохастической матрицы, кратных d.

3. Получено необходимое и достаточное условие на начальное состояние, при котором все предельные векторы равны и в сети достигается равновесие.

Настоящая работа посвящена функционированию циклических сетей при ресурсе, большем порогового значения. Как будет показано, при большом суммарном ресурсе свойство цикличности перестает оказывать влияние на процессы перераспределения ресурса, происходящие в сети, и ее поведение становится сходным с поведением регулярных сетей. Однако в отличие от регулярных сетей, в некоторых циклических сетях

становится возможным по свойствам пропускных способностей определить, какие из вершин являются потенциальными аттракторами.

2. Основные определения

2.1. РЕСУРСНЫЕ СЕТИ

Ресурсная сеть (в дальнейшем просто «сеть») представляет собой ориентированный граф G = (V, Е), V = п, с матрицей пропускных способностей R = (ггу)пхп.

Q(t) = ^(0, ..., qn(t)) - состояние сети в момент ^ qi(t) > 0 -количество ресурса в вершине

Состояние сети Q называется предельным, если оно асимптотически достижимо из некоторого начального состояния Q(0).

г-п и г°и - входная и выходная пропускные способности вершины vi соответственно.

Правила функционирования сети. В момент t вершина vi отдает в смежную ей вершину vm:

гт единиц ресурса, если qi ^) > гоии (правило 1);

Г

—^0mU^qi ^) единиц ресурса, если qi (t) < г1°и‘ (правило 2).

i i

Множество вершин, для которых qi(t) < г°и, называется зоной ^({). Вершины из ZГ(t) функционируют по правилу 2. Z+(t) - множество вершин, ресурс которых больше их выходной пропускной способности, они функционируют по правилу 1. (Эти и другие определения описаны более подробно в ряде работ. Среди них, например, [2, 3, 6]).

Т - пороговое значение ресурса, такое, что при W < Т все вершины, начиная с некоторого t', переходят в зону Z-(t); при W > Т зона ^ (t) не пуста, начиная с некоторого ^.

Для предельного состояния Q обозначим эти зоны через и Z+*.

Введем обозначение: Дг = г’п - г°ш . По знаку Аri вершины делятся на три класса:

1) вершины-приемники, для которых Дгг- > 0;

2) вершины-источники, для которых Дгг- < 0;

3) нейтральные вершины, для которых Дгг- = 0.

Вершины, способные при W > Т из некоторого начального состояния перейти в , называются потенциальными аттракторами. Потенциальными аттракторами могут быть некоторые приемники и некоторые нейтральные вершины, причем первые способны притягивать ресурс и поэтому названы активными аттракторами, вторые могут лишь сохранить ресурс, которым обладали в начальном состоянии. Такие аттракторы называются пассивными. Потенциальные аттракторы в регулярных сетях и предельные состояния при больших ресурсах были исследованы в [2, 3].

В d-циклической сети существует d подпоследовательностей векторов состояний, сходящихся к предельным векторам

ГЛ * ГЛ *

Q1 , ., Qd , последовательно сменяющим друг друга при t ^ да. Если при некотором начальном состоянии они равны, в сети достигается равновесное состояние. Если векторы Q1 , ..., Qd равны при любом начальном состоянии, будем говорить, что в сети существует предельное состояние Q :

в: =... = в:=в*.

Поток в ресурсной сети.

Ресурс, выходящий из вершины vi по ребру (V, V) в момент t, приходит в вершину Vj в момент t + 1; между моментами t и t + 1 он находится в ребре ^, v). Этот ресурс назовем потоком /(0. Общий поток сети описывается матрицей F(t) = (/)(0)пхп.

Величиной потока будем называть сумму:

!шт = 1 1/ (t).

i=1 3=1

Введем следующие обозначения.

£ Л (0 = ГГ (О ; £ Л (г) = Г" (г +1). Кроме того, положим

у=1 г =1

по определению /™ (0) = 0 .

Матрицу предельного потока, если он существует, обозначим через Р. Суммарный предельный поток тогда будет

П П

* = ££ г*

J ^-ит J г]

г=1 ] =1

Р°ш(г) = ((г), ..., (г)) - вектор исходящего потока;

Р'П(0 = (1Г (г), ., ГП" (0) - вектор входящего потока.

Пределы Нт Р°иг(г) и lim Ргп (г) , если они существуют, обо-

г

Т~ЛП * Т°иг*

значим через р и Р соответственно.

2.2. ЦИКЛИЧЕСКИЕ СЕТИ И ИХ ТОПОЛГИЯ

Эргодическая сеть - сеть, граф которой сильно связен. Это определение соответствует определению эргодической цепи Маркова, данному в [4]. Эргодическую сеть, наибольший общий делитель всех циклов которой больше единицы, будем называть циклической. Число d, равное НОД длин всех циклов, - период сети. Множество ее вершин делится на d циклических классов. Стохастическая матрица Я такой сети имеет ровно d собственных значений, равных по модулю единице: М = 1^2 = ••• = |Л*| = 1, т.е. является импримитивной с индексом импримитивности d.

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

Более сложный случай - простые циклы. Простой цикл состоит из нескольких циклов, пересекающихся по вершинам, но не пересекающихся по ребрам. Все эти циклы обязаны не быть

взаимно простыми. Если НОД длин этих циклов равен единице, сеть состоит из одного циклического класса; она регулярна и для нее выполняются все результаты, полученные для регулярных сетей.

Общий случай циклической сети - сеть, состоящая из циклов, пересекающихся либо только по вершинам, либо по вершинам и ребрам.

3. Функционирование циклических сетей с различными топологиями при 1/И > T

3.1. ЭЛЕМЕНТАРНЫЕ ЦИКЛЫ

Рассмотрим на примерах функционирование сети, состоящей из однородного и неоднородного элементарного цикла.

Пример 1. Пусть сеть представлена однородным циклом с пятью вершинами и ее матрица пропускных способностей имеет вид:

(1) Я

' 0 5 0 0 0'

0 0 5 0 0

0 0 0 5 0

0 0 0 0 5

V 5 0 0 0 0,

, = 25.

Она однородна и квазисимметрична, т.е. является эйлеровой нерегулярной сетью. Каждая ее вершина - потенциальный аттрактор. Нетрудно убедиться, что насыщение в этой сети наступает, когда все вершины имеют по г°и‘ = 5 единиц ресурса,

и, таким образом, Т = гшт = 25.

При начальном состоянии Q(0) = (30, 0, 0, 0, 0) имеем рис. 1.

Из рис. 1 видно, что за первые пять тактов все вершины получили ресурс, равный своей пропускной способности (вершина 1 оставила себе излишки W - Т), и на этом изменения компонент вектора состояния завершились. Сеть не прекращает функционировать, но ее состояние остается неизменным.

г

Рассмотрим сеть, представляющую собой неоднородный цикл, с матрицей пропускных способностей

<К0

35 ------

Рис. 1. Функционирование элементарного цикла при W > T

(2) R =

' 0 5 0 0 0'

0 0 4 0 0

0 0 0 3 0

0 0 0 0 4

V 5 0 0 0 0,

= 21.

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

Пример 2. Пусть в сети с матрицей (2) начальное состояние Q(0) = (7, 0, 0, 0, 0) (см. рис. 2).

Поскольку суммарный ресурс превосходит пропускные способности ребер, минимальная из которых г33 = 3, он делится на три неравные части, и дальше в сети наблюдается три цикла длины 5: с ресурсом 3, 3 и 1.

Таким образом, если элементарный цикл неоднороден и ресурс в начальном состоянии аккумулирован в одной из вершин,

г

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

W -

W

шт к

1, і

шт Г:, где

У

^ і

и -

целая часть числа.

Рис. 2. Потоки, не превосходящие минимальный разрез Таблица 1. Протокол работы сети_________________________

t VI *2 *3 *4 *5

0 7,000 0,000 0,000 0,000 0,000

1 2,000 5,000 0,000 0,000 0,000

2 0,000 3,000 4,000 0,000 0,000

3 0,000 0,000 4,000 3,000 0,000

4 0,000 0,000 1,000 3,000 3,000

5 3,000 0,000 0,000 1,000 3,000

6 3,000 3,000 0,000 0,000 1,000

7 1,000 3,000 3,000 0,000 0,000

8 0,000 1,000 3,000 3,000 0,000

9 0,000 0,000 1,000 3,000 3,000

10 3,000 0,000 0,000 1,000 3,000

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

рт _ роШ _ (т^п ,...,т|п г,, ) ; Т = П ■ шт Г-,

где п - число вершин; тт г- - минимальная пропускная способность цикла.

При Ж > Т предельный поток существует и равен Т для любого значения суммарного ресурса, а излишки сверх Т скапливаются лишь в вершинах ук, для которых к _ а^тт ги.

1 1

Так, для сети, заданной матрицей (2), Т = 15.

Пример 3. Пусть в сети с матрицей (2) начальное состояние Q(0) = (15, 0, 0, 0, 0). В таком случае в ней за шесть тактов происходит выравнивание ресурса и стабилизация состояний и потока (рис. 3).

В неоднородном цикле при Ж > Т, как уже было сказано, ресурс скапливается в вершинах с минимальной выходной пропускной способностью ук. Такие вершины будут потенциальными аттракторами цикла. Если аттрактор один, предельное состояние единственно. Ресурс Ж - Т в предельном состоянии окажется в аттракторе (рис. 4). Отсюда следует, что в сети с матрицей (2) потенциальным аттрактором является вершина v3.

Излишек ресурса Ж - Т перешел в третью вершину. Предельное состояние описывается вектором Q = (3, 3, 8, 3, 3).

Интересно, что в последних двух примерах вершина VI, обладая суммарным ресурсом сети в начальном состоянии, сначала отдает его весь и только затем получает часть обратно.

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

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

<7(0

\ —1

\ ^

/ / X/ У» *

Шт Ят 1 *тч 1 Ж ■ ™ 1 1 0 1 2 3 4 5 6 7 '

Рис. 3. Ресурс Ж = Т = 15 выравнивается. Вектор предельного состояния Q =(3, 3, 3, 3, 3)

Пример 4. Пусть в сети с матрицей (2) начальное состояние Q(0) = (20, 0, 0, 0, 0).

Пример 5. Рассмотрим неоднородный цикл из семи вершин с двумя потенциальными аттракторами v2 и v5 - они имеют минимальные выходные пропускные способности.

(3) R =

' 0 5 0 0 0 0 0'

0 0 2 0 0 0 0

0 0 0 5 0 0 0

0 0 0 0 4 0 0

0 0 0 0 0 2 0

0 0 0 0 0 0 6

V 7 0 0 0 0 0 0,

, = 31; Т = 2-7 = 14.

Ж = 17 > Т. Начальное состояние Q(0) = (17, 0, 0, 0, 0, 0, 0).

Из рис. 5 видно, что в зону 2+ перешел тот аттрактор, который первым успел задержать ресурс у себя, т.е. вершина v2.

При начальном состоянии Q(0) = (0, 0, 0, 17, 0, 0, 0) предельным состоянием будет Q = (2, 2, 2, 2, 5, 2, 2), т.е. весь излишек останется во втором потенциальном аттракторе v5.

Если ресурс Ж > Т в начальном состоянии находится в нескольких вершинах, его излишек Ж - Т распределяется между аттракторами в пропорции, которая зависит от пропускных способностей всех ребер.

Рис. 5. Цикл с двумя аттракторами. Вектор предельного состояния О* = (2, 5, 2, 2, 2, 2, 2)

г

Так, например, для матрицы (3) при начальном состоянии Q(0) = (8, 0, 0, 9, 0, 0, 0) предельное состояние будет Q = (2, 2, 2, 2, 5, 2, 2), т.е. все излишки вновь окажутся в аттракторе у5. Чтобы вершина у2 тоже получила излишек, ресурс в вершине у1 должен быть больше 8. При начальном состоянии Q(0) = (9, 0, 0, 8, 0, 0, 0) предельное состояние

Q = (2, 3, 2, 2, 4, 2, 2). Однако при перемещении ресурса, равного 8, в начальном состоянии из вершины у4 в у3, т.е. при Q(0) = (9, 0, 8, 0, 0, 0, 0), предельным состоянием вновь будет вектор Q = (2, 2, 2, 2, 5, 2, 2) - все излишки забирает аттрактор у5. Это происходит оттого, что ресурс должен проделать более длинный путь по циклу до возвращения в его начало: вершину у1, служащую источником для аттрактора у2. При Q(0) = (9, 0, 8, 0, 0, 0, 0) вершина у1 пуста три такта (со второго по четвертый), и этого достаточно, чтобы ресурс из аттрактора у2 переместился в аттрактор у5 (рис. 6). Из рисунка видно, что к тому моменту (такт 5), как ресурс вершины у1, становится вновь ненулевым и она могла бы поддерживать баланс в вершине у2, вершина у2 уже имеет ресурс, равный единице, т.е. не имеет излишков.

4(0

10

9 8

б

5

4

3

1

о

0123456789

Рис. 6. Цикл с двумя аттракторами. Вектор предельного состояния О* = (2, 2, 2, 2, 5, 2, 2)

Обобщим все результаты, полученные при рассмотрении элементарных циклов при больших ресурсах.

У т в е р ж д е н и е 1. Для циклической ресурсной сети с n вершинами, представленной элементарным циклом:

1. Пороговое значение ресурса Т вычисляется по формуле: T = n ■ rmin, где rmm = min Гу.

2. При W = T предельный поток и предельное состояние существуют и единственны. Вектор предельного потока равен

rin* T~°ut* ( \

вектору предельного состояния: F = F = Q = (rmin, • • •, Г min).

3. При W > T предельное состояние и поток существуют, причем предельный поток определяется единственным образом: Fn = F°ut = (rmin, ..., rmin). Предельное состояние единственно в том и только в том случае, когда сеть имеет единственный потенциальный аттрактор. В противном случае распределение ресурса сверх порогового значения W - T зависит от начального состояния.

4. Вершина vk является потенциальным аттрактором элементарного цикла в том и только в том случае, если для нее выполняется k = argmin r. Если два и более потенциальных аттракто-

i 1

ра расположены в цикле непосредственно друг за другом, то активным аттрактором является только тот, для которого г"1 > r°ut; остальные аттракторы пассивны.

3.2. ПРОСТЫЕ ЦИКЛЫ И ЦИКЛЫ С ОБЩИМИ РЕБРАМИ Рассмотрим функционирование неоднородной сети, состоящей из двух циклических классов. Циклы в ней пересекаются только по вершинам.

Рис. 7. Двудольный граф с пятью вершинами, d = 2

Пример 6. Пусть сеть за

' 0 1 0 0 0 "

2 0 0 2 0

(3) R = 0 0 0 2 3

0 2 1 0 0

V 0 0 2 0 0,

= 15.

При заведомо большом значении ресурса (большем гшт) W = 20 любое начальное состояние приводит к одному и тому же предельному состоянию Q = (14,5; 2; 1,25; 1,5; 0,75).

Отсюда видно, что потенциальный аттрактор этой сети единственен, им является вершина VI. Пороговое значение ресурса для этой сети находится из того условия, что при W > Т сумма всех компонент вектора предельного состояния, кроме первой, постоянна и равна 5,5, а первая компонента при W = Т равна суммарной выходной пропускной способности вершины, т.е. единице: Т = 6,5. При W = Т предельное состояние описывается вектором Q = (1; 2; 1,25; 1,5; 0,75) для любого начального состояния. Этим же вектором описывается поток при любом W > Т.

При суммарном ресурсе, меньшем 6,5, в сети существует два предельных вектора.

Векторы Q11* и б2*, формирующие предельные матрицы

и R2‘ю, имеют вид:

= (0, 0,616, 0,384, 0, 0), Q1* = (0,308, 0, 0, 0,462, 0,230).

Собственный вектор стохастической матрицы R' ищется по формуле

(4) = б + б1* .

Ql* = (0,154, 0,313, 0,192, 0,231, 0,115).

Проверим, выполняется ли для вершины VI признак аттрак-тивности, введенный для вершин регулярных сетей.

Для аттрактора с номером д должно выполниться условие:

г

к = а^тт-

Вычислим значение

Чг

для всех вершин сети. В качестве

вектора предельного состояния возьмем вектор (4).

'і

1

і*

Ч

омі

4

1*

Ч4

0,154

3

0,231

>6,5.

13

4

1*

Ч2

Ч5

0,313

2 ^ 0,115 ~

17,4.

5

Ч3

0,192

>26.

Как видно, минимум этого отношения действительно достигается на первой вершине, более того, его значение равно 6,5, что совпадает со значением Т, как это и было доказано для регулярных сетей [3].

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

4. Пороговое значение Т, предельные состояния и потоки в циклических сетях при больших ресурсах

В регулярных сетях существование и единственность порогового значения Т вытекали из единственности вектора предельного состояния при малых ресурсах. Все компоненты вектора б при увеличении W вплоть до значения Т увеличивались пропорционально, пока одна или несколько из них не принимали значения, равные своим суммарным выходным пропускным способностям. Это значение суммарного ресурса и было пороговым. Вершины с ресурсом, достигшим выходных пропускных способностей, - потенциальные аттракторы сети. При дальней-

г

Ч

г

г

г

г

шем увеличении ресурса пропорциональность компонент нарушается, и все излишки сверх Т накапливаются в аттракторах.

В циклической сети существует d предельных векторов, в общем случае не равных друг другу. Кроме того, каждый из этих векторов зависит от начального состояния. Единственность порогового значения в такой сети уже неочевидна.

Будем говорить, что если в циклической сети d предельных векторов равны при любом начальном состоянии, в ней существует предельное состояние Q : Q1 = ... = Qd = Q .

Т е о р е м а 1. В эргодической d-циклической ресурсной сети существует единственное пороговое значение Т, такое, что:

при W < Т все вершины за конечное число тактов переходят на правило 2, и в сети имеется d предельных векторов;

при W > Т потоки стабилизируются, и в сети достигается глобальное равновесие при любом начальном состоянии. Предельный поток существует и единственен. Предельное состояние существует; оно не зависит от начального состояния в том и только в том случае, когда сеть имеет один аттрактор.

Д о к а з а т е л ь с т в о. В сети, все вершины которой функционируют по правилу 2, существует предельный цикл из d предельных векторов. Если в начальном состоянии некоторые вершины функционируют по правилу 1, но все они за конечное число тактов N переходят на правило 2, то, взяв Q(N + 1) за новое начальное состояние, вновь получим предельный цикл из d векторов. Верхняя граница суммарного ресурса, при котором все вершины функционируют по правилу 2, равна Т. Т существует для любой сети.

Докажем, что значение Т единственно и не зависит от начального состояния. При W = Т по крайней мере одна вершина имеет ресурс, равный своей выходной пропускной способности. Пусть это вершина Vj. Для нее д* = г°иі в каждом из d предельных векторов.

Матрица R' импримитивна с индексом импримитивности d. Это означает, что ее степень Я1 разлагается в прямую сумму d регулярных матриц, соответствующих циклическим классам Я/ [1]. Предельные векторы сети выражаются через собственные векторы этих матриц. Вершина V. является потенциальным аттрактором в некотором циклическом классе с номером /. Если с первого такта все вершины функционировали по правилу 2,

для него выполняется: б/ (0)Я™ = б* = Wiжi, где

Wi - суммарный ресурс в /-м циклическом классе; п/ - вектор предельных вероятностей соответствующей цепи Маркова. Для сети, описываемой регулярной матрицей Я/, пороговое значение существует и единственно. Обозначим его через Т/. Тогда при W/ = Т/ имеем (д* = г°ии, где (д*). - .-я компонента векто-

ра б* . Но чтобы ресурс вершины V. был равен г°ш на каждом

такте, суммарный ресурс /-го циклического класса должен быть равен Т/ также на каждом такте. Если сеть функционирует по правилу 2, ресурс циркулирует по циклическим классам, не перемешиваясь. Если раскрасить ресурс в каждом классе в свой цвет, ресурсы разных цветов на любом такте будут целиком в одном классе. Пороговое значение Т определяется как максимальный ресурс, при котором все вершины функционируют по правилу 2 при t ^ да. Тогда при t ^ да каждый циклический класс на каждом такте должен содержать Т/ ресурса. А пороговое значение ресурса сети вычисляется по формуле Т = 1Т/.

Поскольку каждая компонента имеет одинаковое количество ресурса, поток в сети стабилизируется, и все предельные векторы одинаковы.

При W > Т как минимум один аттрактор (пусть это вновь будет вершина V.) за конечное число тактов переходит на правило 1. То есть существует такой момент времени t', что при t > t вершина V. функционирует по правилу 1. На каждом такте она отдает ресурс, равный . , и получает ресурс, не меньший . . Тогда в нее можно добавить петлю с пропускной способностью

qjt) - r°ut, по которой на каждом такте будет проходить неотрицательный поток, равный qjt) - r°ut. Функционирование сети при этом не изменится. Но наличие петли влечет за собой потерю цикличности, и сеть превращается в регулярную. В регулярной сети поток существует и единственен. Предельное состояние единственно во всех вершинах, кроме потенциальных аттракторов. Если в сети единственный аттрактор, предельное состояние также единственно. □

Докажем, что для вершин в произвольной d-циклической сети выполняется тот же критерий аттрактивности, что и для вершин в несимметричной регулярной сети.

Т е о р е м а 2. В эргодической d-циклической ресурсной сети вершина Vj является потенциальным аттрактором, если и

out

Г- 1 *

только если j = argmin-1—, где вектор Q определяется по

1 qt

формуле

1 d

(5) Q1 = 1Ш,

d k=1

векторы Q1 (k = 0, ..., d) - предельные векторы при W = 1 и произвольном начальном состоянии.

Д о к а з а т е л ь с т в о. Из теоремы 1 следует, что при W = T предельное состояние сети при любом начальном состоянии единственно, и все ее d предельных векторов равны. Но для W< T выполняется равенство: Q* + ... +Qd = dTQ1*, где вектор Q , вычисляемый по формуле (5), - единственный положительный собственный вектор стохастической матрицы R'. Непосредственно отсюда следует, что предельное состояние сети Q = Q1 = ... = Qd = TQ . Таким образом, функционирование сети начинает описываться регулярной матрицей A, состоящей из n строк Q1 :

A = 1- Q1*.

Для сети, функционирование которой описывается матрицей А, справедливы все результаты, полученные для регулярных сетей. Таким образом, для того, чтобы вершина Vj была аттрак-

тором, необходимо и достаточно выполнение усло-

out

Г-

вия j = argmin-1— [3]. □

- Чг

С л е д с т в и е 1. В эйлеровых d-циклических сетях при больших ресурсах предельное состояние полностью зависит от начального, так как каждая вершина в них - потенциальный аттрактор.

С л е д с т в и е 2. В эргодической d-циклической ресурсной сети пороговое значение Т определяется по формуле

out

Г

T = min^.

- Чг

Таким образом, при W > T для циклических сетей можно сформулировать теоремы о предельном состоянии и потоке, аналогичные теоремам о регулярных сетях. Их доказательства следуют из теорем 1-2 и следствия 2.

Т е о р е м а 3 (о п р е д е л ь н о м с о с т о я н и и). В эргодической d-циклической ресурсной сети при W > T значения неаттрактивных компонент вектора предельного состояния Q = (чг , ..., Чп ) вычисляются по формуле

* 1* ЛТ-Т . .

4t = 4i • T , i ^jk, jk - номера аттракторов,

где:

1) вектор Qг* определяется по формуле (5),

out

2) T = min -^,

- Чг

out

Г-

3) jk определяются из условия: jk = argmin -1— .

- Чг

Оставшийся ресурс распределяется между потенциальными аттракторами.

Проиллюстрируем теорему следующим примером.

Пример 7. Рассмотрим функционирование неоднородной сети, содержащей циклы с общими ребрами.

Рис. 8. Циклы с общими ребрами

В этой сети четыре циклических множества: (уь у5},

V V6}, ^3} и ^4}.

Пусть матрица пропускных способностей сети имеет вид

(0 1

R =

000 3 0 0 0 2 2 0 0 0

00

00

0000

0 0 1 0 0

0 ^ 0 0 0 4 0

= 18.

Тогда при помещении в начальном состоянии ресурса Ж = 1 в любую из шести вершин предельные векторы будут иметь вид е;* = (0,286; 0; 0; 0; 0,714; 0),

Q1* = (0; 0,286; 0; 0; 0; 0,714),

= (0, 0, 1, 0, 0, 0), е4* = (0, 0, 0, 1, 0, 0).

Собственный вектор стохастической матрицы находится как их среднее арифметическое (формула (5)):

е1* = (0,0715; 0,0715; 0,25; 0,25; 0,1785; 0,1785).

Вычислим значение отношения

Чг

для каждой вершины.

1

Ч1

0,0715

14.

ОШ

Г2

1*

Ч2

3

0,0715

42.

Чз

0,25

= 8.

г

г

2

г

=<

rOUt 7 rOUt 4 rOUt i

*------= 28, rL- *---------* 22,4, ^ *------------* 5,6.

q\* 0,25 q\* 0,1785 q'* 0,1785

out

r

Минимум -1-— достигается в вершине v6. qг

Таким образом, v6 - единственный потенциальный аттрактор этой сети. По следствию 2 из теоремы 1 T * 5,6.

При W = T предельное состояние при любом начальном:

Q = (0,4; 0,4; 1,4; 1,4; 1; 1).

Этому же вектору равны предельные входной и выходной потоки.

При W > T, например, W = 20, при любом начальном состоянии предельное состояние равно Q* = (0,4; 0,4; 1,4; 1,4; 1; 15,4).

Излишек ресурса аккумулировался в аттракторе. Остальные компоненты предельного вектора не изменились.

Т е о р е м а 4 (о п р е д е л ь н о м п о т о к е). В эргоди-ческой d-циклической ресурсной сети при W > T предельный поток существует и единственен и определяется по формуле:

rin* rout* ,-Л*^ . Trin* TpOUt* /о1*^ ~

Ji = Ji = q,1; F = Q T = Q,

где

1) вектор Ql* определяется по формуле (5),

out

2) T = min -L**-.

1 qt

З а м е ч а н и е. В регулярных сетях при разделении ресурса на «малый» и «большой», пороговое значение Т, по которому проходит эта граница, относился к малым ресурсам. Это связано с тем, что при W = T сеть описывается регулярной цепью Маркова, как и при меньших значениях суммарного ресурса, и предельное состояние в ней единственно. То есть в регулярных сетях суммарный ресурс разбивается на значения W < T и W > T. В циклических сетях значение Т входит в «большой» ресурс, так как при W > T сеть имеет предельное состояние, в то время как при W< T - лишь предельный цикл из d предельных векторов.

5. Заключение

В работе исследованы процессы перераспределения больших ресурсов в эргодических циклических сетях. Рассмотрены примеры функционирования сетей с большим ресурсом, имеющих различные топологии. Доказано, что циклические сети, как и регулярные, имеют единственное пороговое значение ресурса Т, не зависящее от начального состояния, такое что при W < Т все вершины переходят на правило 2, в сети существует d предельных векторов, в общем случае различных. При W > Т предельные векторы равны для любого начального состояния. Предельный поток единственен, суммарный поток равен Т. Найдены формулы для векторов предельного состояния в неаттрактивных вершинах и потока - во всех вершинах. Сформулирован и доказан критерий аттрактивности вершин.

Из полученных результатов следует, что при больших ресурсах функционирование сети становится регулярным и описывается регулярной матрицей А.

Литература

1. ГАНТМАХЕР Ф.Р. Теория матриц. - М.: Физматлит. 2004. - 560 с.

2. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети.

II. Потоки при больших ресурсах и их стабилизация // Автоматика и телемеханика. - 2012. - №6. - С. 103-118.

3. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети.

III. Исследование предельных состояний // Автоматика и телемеханика. - 2012. - №7. - С. 67-77.

4. ЖИЛЯКОВА Л.Ю. Эргодические циклические ресурсные сети. I. Колебания и равновесные состояния при малых ресурсах // Управление большими системами. - Вып. 43. -М.: ИПУ РАН, 2013. - С. 34 - 54.

5. КЕМЕНИ ДЖ., СНЕЛЛ ДЖ. Конечные цепи Маркова. - М.: Наука, 1970. - 271 с.

6. КУЗНЕЦОВ О.П. Однородные ресурсные сети. I. Полные графы // Автоматика и телемеханика. - 2009. - №11. -С. 136-147.

7. РОБЕРТС Ф.С. Дискретные математические модели с приложениями к социальным, биологическим и экологическим задачам. - М. Наука, 1986. - 496 с.

8. BJORNER A., LOVASZ L., SHOR P. Chip-firing games on graphs // Europ. J. Comb. - 1991. - №12. - P. 283-291.

9. BJORNER A., LOVASZ L. Chip-firing game on directed

graphs // J. Algebraic Combinatorics. - 1992. - №1. -

P. 305-328.

10. BLANCHARD PH., VOLCHENKOV D. Random Walks and Diffusions on Graphs and Databases: An Introduction (Springer Series in Synergetics). - Springer-Verlag - Berlin-Heidelberg, 2011. - 262 p.

11. LOVASZ L., WINKLER P. Mixing of Random Walks and Other Diffusions on a Graph // Surveys in Combinatorics / Ed. P. Rowlinson. - London Math. Soc. Lecture Notes Series 218. -Cambridge Univ. Press, 1995. - P. 119-154.

12. PRISNER E. Parallel Chip Firing on Digraphs // Complex

Systems. - 1994. - №8. - P. 367-383.

ERGODIC CYCLIC RESOURCE NETWORKS. II. HIGH LEVELS OF RESOURCE

Ludmila Zhilyakova, Institute of Control Sciences of RAS, (Moscow, Profsoyuznaya st., 65), cand. sc., senior scientist, zhilyakova. ludmila@gmail .com.

Abstract: We study processes in ergodic nonregular resource networks at high levels of resource. The formula for the threshold value T is derived. It was shown that there are no oscillations if W > T; the network becomes stable from every initial state when resource level is high. The formulae for the vectors of limit state and the limit flow are derived. The criterion of attractivity of vertices is formulated.

Keywords: diffusion on graphs, resource network, attractors, limit state.

Статья представлена к публикации членом редакционной коллегии О.П. Кузнецовым

Поступила в редакцию 06.02.2013.

Опубликована 30.09.2013.