Эргодические циклические ресурсные сети. I. колебания и равновесные состояния при малых ресурсах Текст научной статьи по специальности «Математика»

УДК 519.177 ББК 22.18

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

I. КОЛЕБАНИЯ И РАВНОВЕСНЫЕ СОСТОЯНИЯ ПРИ МАЛЫХ РЕСУРСАХ1

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

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

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

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

1. Введение

Работа посвящена исследованию свойств нерегулярных эргодических ресурсных сетей и их функционирования при малых ресурсах. Ресурсная сеть - динамическая пороговая потоковая модель. Регулярные ресурсные сети описаны в работах [3-5, 7],

1 Работа поддержана грантом РФФИ 11-01-00771.

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

там же даны основные определения. Было доказано, что в регулярных сетях существует пороговое значение ресурса Т, единственное для каждой сети: при суммарном ресурсе, не превосходящем Т (такой ресурс в дальнейшем будем называть малым), процессы в ресурсной сети эквивалентны случайным блужданиям на графе, хотя такая сеть является полностью детерминированной моделью. Случайные блуждания и модели рассеяния на графах используются для решения различных задач в широком кругу предметных областей. Исчерпывающий обзор таких моделей, исследование их свойств и описание областей приложения были проведены в работе [10].

Если суммарный ресурс выше порогового значения, некоторые вершины переключаются на другое правило функционирования. Это правило сходно с тем, как «ходят» вершины в модели, названной «играми выстреливания фишек» (chip-firing games) [9, 11]. Chip-firing game - графовая целочисленная пороговая модель, в которой в вершинах находятся стопки фишек. Вершина может сделать ход, выстреливая по одной фишке вдоль каждого исходящего ребра, только в том случае, если она имеет достаточно фишек. Выстреливание может последовательным и параллельным, как, например, описано в работе [12].

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

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

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

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

Сеть представляет собой ориентированный граф G = (V, Е), IV = „, с матрицей пропускных способностей Я = (Гу)„х„.

Q(t) = (д1(:), ..., q„(t)) - состояние сети в момент :; qi(t) > 0 -количество ресурса в вершине vi. Q - предельное состояние.

„ „

Г„ = ^ г}1 и г1ои: = ^ ^ - входная и выходная пропускные

}=1 }=1

способности вершины vi соответственно.

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

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

Г

™ qi (:) единиц ресурса, если qi (:) < г°ш (правило 2).

„ОШ 11

Матрице пропускных способностей R соответствует стохастическая матрица R', получаемая из R нормированием по строкам:

R' =

( Г '11 Г12

Г°Ш ГОШ

Гп1 Гп2

оШ Г V п оШ Г п

1 п

Сеть, все вершины которой функционируют по правилу 2, описывается цепью Маркова с соответствующей стохастической матрицей [2]. Закон перераспределения ресурса описывается рекуррентной формулой:

+ і) = ат.

Сеть будем называть однородной , если пропускные способности всех ее ребер одинаковы; эйлеровой, если для каждой

г

1 Определение однородной ресурсной сети следует отличать от определения однородной цепи Маркова; функционирование всех

вершины r,in = r,ou‘, i = 1, ..., n; несимметричной, если в ней существуют вершины с разными входными и выходными способностями.

W - суммарный ресурс сети. Сеть функционирует с выполнением закона сохранения: W = const.

2.2. ЭРГОДИЧЕСКИЕ ЦИКЛИЧЕСКИЕ СЕТИ

При малых ресурсах сеть описывается соответствующей цепью Маркова. Следуя классификации цепей Маркова, данной в [6], эргодической сетью назовем сеть, граф которой сильно связен; эргодической компонентой неэргодической сети назовем сильно связную компоненту, из которой нет выходящих ребер. Если дополнительно к этому НОД длин всех циклов графа равен единице, сеть будем называть регулярной. Эргодическую нерегулярную сеть, НОД длин всех циклов которой строго больше единицы, будем называть циклической.

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

В эргодической цепи Маркова выделяют циклические классы состояний [6]. С течением времени текущее состояние движется по циклическим классам в определенном порядке, причем каждые d шагов система оказывается в одном и том же циклическом классе. Количество циклических классов d называют пе-

классов ресурсных сетей при малых ресурсах описываются однородными цепями Маркова.

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

Стохастическая матрица Я циклической сети импримитив-на с индексом импримитивности d [2]. Это означает, что Я имеет d различных собственных чисел, равных по модулю 1.

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

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

3.1. ФУНКЦИОНИРОВАНИЕ ОДНОРОДНОГО ПРОСТОГО ЦИКЛА ПРИ W = 1

Пример 1. Пусть сеть представлена двудольным графом (рис. 1).

Рис. 1. Двудольный граф с пятью вершинами, d = 2 Ее матрица пропускных способностей:

(1) R =

(0 2 0 0 0Л

2 0 0 2 0

0 0 0 2 2

0 2 2 0 0

V 0 0 2 0 0,

, = 16.

Соответствующая ей стохастическая матрица имеет вид ( 0 1 0 0 0 ^

0,5 0 0 0,5 0

R' = 0 0 0 0,5 0,5

0 0,5 0,5 0 0

0 0 10 0

ч У

На рис. 2 показан процесс перераспределения ресурса, соответствующий начальному состоянию Q(0) = (1, 0, 0, 0, 0).

Ф)

Рис. 2. Колебания при Q(0) = (1, 0, 0, 0, 0)

Протокол функционирования сети, соответствующий рис. 2, представлен в таблице 1.

г

Таблица 1. Протокол функционирования 2-циклической сети с

матрицей (1) и ^і(0) = (1, 0, 0, 0, 0)

V! V2 Vз V4 V5

0 1 0 0 0 0

1 0 1 0 0 0

2 0,5 0 0 0,5 0

3 0 0,75 0,25 0 0

4 0,375 0 0 0,5 0,125

5 0 0,625 0,375 0 0

6 0,313 0 0 0,5 0,188

7 0 0,563 0,438 0 0

8 0,281 0 0 0,5 0,219

9 0 0,531 0,469 0 0

10 0,266 0 0 0,5 0,234

11 0 0,516 0,484 0 0

12 0,258 0 0 0,5 0,242

13 0 0,508 0,492 0 0

14 0,254 0 0 0,5 0,246

15 0 0,504 0,496 0 0

16 0,252 0 0 0,5 0,248

17 0 0,502 0,498 0 0

18 0,251 0 0 0,5 0,249

19 0 0,501 0,499 0 0

20 0,25 0 0 0,5 0,25

21 0 0,5 0,5 0 0

22 0,25 0 0 0,5 0,25

23 0 0,5 0,5 0 0

Из рис. 2 и таблицы 1 видно, что на четных тактах ресурс содержится в вершинах ^ и ^, на нечетных - в V2 и v3. В этой сети й = 2, циклическими классами являются множества V V4, V5} и ^2, Vз}.

В такой сети отсутствует равновесное предельное состояние. Протокол функционирования показывает, что в пределе получается два вектора, последовательно сменяющие друг друга: дГ = (0,25; 0; 0; 0,5; 0,25) и Й* = (0; 0,5; 0,5; 0; 0).

Помещая ресурс W = 1 в начальном состоянии поочередно в

61*

1 и

61* ^

2 , одинаковые для всех начальных состояний:

0(0) = (1, 0, 0, 0, 0), 0(0) = (0, 1, 0, 0, 0), 6(0) = (0, 0, 1, 0, 0), 6(0) = (0, 0, 0, 1, 0), 6(0) = (0, 0, 0, 0, 1).

Поскольку й = 2, последовательность степеней стохастической матрицы Я' состоит из двух сходящихся подпоследовательностей:

Я™ =

' 0 0,5 0,5 0 0 '

0,25 0 0 0,5 0,25

0,25 0 0 0,5 0,25

0 0,5 0,5 0 0

V 0 0,5 0,5 0 0 ,

' 0,25 0 0 0,5 0,25'

0 0,5 0,5 0 0

0 0,5 0,5 0 0

0,25 0 0 0,5 0,25

ч 0,25 0 0 0,5 0,25у

я? =

Матрицы Я1'ш и Я2'ш состоят из чередующихся векторов Q11* и Q21*, причем Я2'ш представляет собой инверсию Я/”: в ней каждый вектор Q1 заменен вектором Q2 , и наоборот. Строки

1, 4, 5, соответствующие вершинам первого циклического класса, равны; то же самое относится и к строкам 2, 3, соответствующим вершинам второго циклического класса.

Оба вектора Q1 и Q2 являются собственными векторами

п ” п ҐЛ 1* ҐЛ 1* ҐЛ 1*

матрицы Я2 , а матрица Я1 переводит вектор Ql в Q2 , Q2 -

в Q11 . То же самое происходит и при умножении этих векторов

на матрицу Я':

Ql1*я = Ql*, Q!*Я,= й1*.

Умножение на стохастическую матрицу Я' переводит матрицу Я1'” в Я2'” и наоборот.

Заметим, что сумма двух предельных матриц дает матрицу с одинаковыми строками:

' 0,25 0,5 0,5 0,5 0,25'

0,25 0,5 0,5 0,5 0,25

0,25 0,5 0,5 0,5 0,25

0,25 0,5 0,5 0,5 0,25

ч 0,25 0,5 0,5 0,5 0,25 ,

я1” + я2” =

Вектор, образующий строки этой матрицы, является левым собственным вектором обеих матриц Я/” и Я2'“, а также стохастической матрицы Я'; более того, если в матрице Я1'ш + Я2'“ нормировать строки, окажется, что она равна предельной матрице А, к которой суммируется по Чезаро последовательность Я'к

[5]:

(2)

1 к

А = Нт—^ Я]

к к ]=1

, я1ш+я;ш

А = —------—

' 0,125 0,25 0,25 0,25 0,125'

0,125 0,25 0,25 0,25 0,125

0,125 0,25 0,25 0,25 0,125

0,125 0,25 0,25 0,25 0,125

ч 0,125 0,25 0,25 0,25 0,125,

Матрица Я' неразложима по построению. Неразложимая неотрицательная матрица имеет единственный неотрицательный собственный вектор (с точностью до скалярного множителя) [2]. Собственный вектор стохастической матрицы сети определяется как сумма векторов 61 и 62 . Для доказательства этого факта достаточно сложить равенства 61 Я' = 621 и 621* Я' = 611*.

Вектор абсолютных предельных вероятностей находится по формуле:

6"=2(6Г + 6Г).

2

3.2. ПРЕДЕЛЬНЫЕ ВЕКТОРЫ И ЦИКЛЫ D-ЦИКЛИЧЕСКОЙ СЕТИ

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

В [6] доказано, что предельная матрица А, к которой суммируется по Чезаро последовательность R*, (формула (2)) имеет вид:

A = 1-а,

где 1 - вектор-столбец, состоящий из единиц, а а является единственным неотрицательным левым собственным вектором матрицы R', соответствующим собственному числу Л = 1:

(3) aR = а.

Также выполняется

AR = RA = А.

Пусть R' - стохастическая матрица эргодической

d-циклической ресурсной сети. Рассмотрим последовательность ее степеней: R', ..., R'd, R'd+1, ... Она состоит из d сходящихся подпоследовательностей:

1) R, R'd-R, Rad-R', R'3d-R, ...

2) R'2 R"d R|2 Ri^d ri2 Ri3d ri2

d) R"d R"d R^d R&d Rid Ri3d Rid

Таким образом, видно, что все пределы R1'“, ..., RJX выра-

П |0

жаются через одну предельную матрицу Rd :

Rl« = R° R ;

R° = R°R2;

R ш = R шRd~1

Сумма матриц R1'“, ..., RJX строго положительна. В [1] доказано, что предельная матрица А (2) определяется по формуле

(4) A = — (E + R +... + Rd-1)R° = -Rd0(E + R' +... + Rd_1) .

d d

Для любого начального распределения единичного ресурса выполнится:

(5) а = 61 (0)А .

Формула (5) означает, что каждая строка матрицы А соответствует предельным вероятностям цепи Маркова в смысле среднего арифметического. Далее при описании ресурсных сетей вместо символа а, обозначающего предельный вероятностный вектор, будем использовать обозначение 6 - вектор предельного состояния при единичном ресурсе: 6 = 6 (0)А; матрица А состоит из п одинаковых вектор-строк, равных 6 .

Рассмотрим процесс функционирования сети по правилу 2. Последовательность векторов состояний делится на й сходящихся подпоследовательностей, так же как и последовательность степеней стохастической матрицы.

6(кй + 1) = 6(0)Я'Ы Я',

6(кй + 2) = 6(0)Я'Ы Я'2,

6(кс!) = 6(0)Я'Ы Я'й = 6(0)Я'(Ш)й к = 0, 1, ...

Переходя к пределу, получим:

(6) 6* = 6(0)Я'ГЯ', 6* = 6(т'ГЯ'2, . 6* = 6(0)Яй°.

Непосредственно из выражений (6) следует, что

6*+1 = 6г Я', \ =1, ., й - 1, 6*1 = 6\ Я\

Векторы 61, •••, 6 й будем называть предельными векто-

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

Таким образом, доказано следующее предложение.

Пр е д л о ж е н и е 1. В эргодической й-циклической ресурсной сети, функционирующей по правилу 2, существует предельный цикл длины й, состоящий из предельных векторов 61, ..., 6 а, определяемых по формуле (6).

Пр е д л о ж е н и е 2. Для эргодической й-циклической ресурсной сети, функционирующей по правилу 2, единственным положительным собственным вектором матрицы Я' при Ж = 1

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

векторов и определяемый по формуле

п* п*

(7) еі* = .

d • W

Доказательство. Суммируя выражения (6), имеем:

а*+...+а*=а(оЯ (е+я' +...+я*-1).

Преобразуя правую часть по формуле (4), получим:

а*+...+а*=*а(0) а .

Но любой вектор 2(0) при условии, что сеть с нулевого такта функционирует по правилу 2, можно представить в виде а(0)=Жа:(0). Отсюда

а* +... + а* = dWQl(0)а или **

+...+а = а*(0) а .

dW

а* + + а*

Из формулы (5) имеем: а1* = ——ё . То есть матрица

А состоит из п вектор-строк, определяемых формулой (7). А из (3) следует, что 0>1'Я' = а*. □

При умножении на матрицу Я', каждый вектор а г переходит в вектор а г+1 (сложение индексов происходит по модулю ё). Но существует матрица, которая оставляет все эти векторы неподвижными.

Пр е д л о ж е н и е 3. Предельные векторы эргодической ё-циклической сети а * ..., а ё являются собственными векторами матрицы Яа™, соответствующими собственному числу Л = 1 кратности ё.

Доказательство. Запишем равенства (6) в виде

(8) а*=а(0)ях, і = і, ..., а.

Умножим обе части (8) на матрицу Ясправа. а*Яё" = а(0)Яё" ЯгЯ;, но по формуле (4) это можно переписать в виде

аХ"=а(0)(кх°°Я,г'=а(0)Яё°Яг

По формуле (8) правая часть равенства равна а г. Тогда

6*яО° = 6*. □

Замечание. В предельном случае, когда НОД всех циклов сети равен единице и, соответственно, сеть представляет собой один циклический класс (является регулярной), доказанные утверждения остаются верными. В сети, функционирующей по правилу 2, существует единственный вектор предельного состояния 6 . Он является собственным вектором стохастической матрицы Я (предложение 1), для него верно равенство (7): 6 = 6 /Ж из предложения 2, и он является собственным вектором единственной предельной матрицы Я , соответствующим простому собственному числу Л = 1 (предложение 3). Справедливость всех этих утверждений следует непосредственно из свойств регулярных цепей Маркова. Для ресурсных сетей они сформулированы в [3].

3.2. ДОСТИЖЕНИЕ ГЛОБАЛЬНОГО РАВНОВЕСИЯ ПРИ МАЛЫХ РЕСУРСАХ

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

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

Л е м м а 1. В эргодической й-циклической ресурсной сети, функционирующей по правилу 2, предельные векторы 6 1, •••, 6 й при постоянном суммарном ресурсе в циклических классах не зависят от распределения ресурса по вершинам внутри классов в начальном состоянии.

Доказательство. Матрица Я импримитивна с индексом импримитивности й. В [2] доказано, что ее степень Я'й разлагается на й примитивных матриц, которые имеют одно и то же максимальное характеристическое число. Это означает, что существует такая нумерация вершин, при которой Я'й блочно-диагональна. Блоки соответствуют циклическим классам.

' Лп 0 . .. 0 "

= 0 и’ё Л22 . .. 0

V 0 0 . Л ё • • )

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

Л/1 = Нт (л'ё Т также является блочной с блоками (л-) .

k^0

Каждая пара блоков и (Лиё)° имеет один и тот же единственный положительный предельный вектор щ, соответствующий максимальному собственному числу Л = 1, - = 1, ..., ё.

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

Любой вектор начального состояния 0(0) можно разложить в сумму векторов Qi(0), каждый из которых соответствует одному циклическому классу. Тогда поскольку матрицы регулярны, при Ж = 1 для любого Qi(0) имеют место равенства

й (0)Лг= 0* = пг.

Для произвольного Ж при условии, что сеть с первого такта функционирует по правилу 2, для любого 0г(0) выполнится:

(ОЛ00 = 0* = , где Ж - суммарный ресурс в --м цик-

лическом классе.

Прямая сумма векторов 0 ё даст единственный вектор 0 ё. Остальные векторы предельного цикла получаются из 0 ё умножением на соответствующие степени матрицы Л' и тоже определяются единственным образом. □

Л е м м а 2. В эргодической ё-циклической ресурсной сети, функционирующей по правилу 2, при суммарном ресурсе Ж = 1, сосредоточенном в начальном состоянии в т-м циклическом классе, предельные векторы 0 ь ..., являются строками предельных матриц Л/” каждая из которых состоит из ё

различных строк, причем на строках, соответствующих номерам вершин из i-го циклического класса, находится вектор Q k, где k є [l, d],

(9) k = m + j + i - 3 (mod d) + l, i, j = l, ..., d.

Доказательство. Рассмотрим начальное состояние, в котором ресурс находится в первом циклическом классе (m = І). Не нарушая общности, предположим, что этому классу принадлежит вершина vl. Пусть Q(0) = (І, 0, ..., 0). Обозначим этот вектор начального состояния через Ql(0).

Имеем:

Q* = Ql(0)Rl"= R )і,

Q* = Ql(0)R2" = (r2“ )l,

Q* = R )l,

где (R/'“)l - первая строка матрицы R/".

По лемме І, вместо вершины vl может быть взята любая вершина из первого циклического класса. Таким образом, в матрицах Rj в строках, соответствующих номерам вершин из первого циклического класса, стоят векторы Q j.

Для каждого вектора Q i+l (i = І, ..., d- І) справедливо Q i+l = Q i R, Q l = Q d R ; для предельных матриц выполняется Rj+ І'ш = R/“ R'. Это означает, что каждая матрица R/” в строках, соответствующих номерам вершин из i-го циклического класса, содержит векторы с номером j + i - І при j + i - І < d, и с номером j + i - І - d при j + i - І > d. Индекс k вектора Qk , удовлетворяющий этим условиям, находится по формуле k = j + i - 2 (mod d) + l.

Обобщая эту формулу на произвольный класс с номером m, получим

k = m + j + i - 3 (mod d) + l. □

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

Л е м м а 3. В эргодической d-циклической ресурсной сети, функционирующей по правилу 2, предельные векторы Q l, ..., Q а совпадают, если в каждом из циклических классов в

начальном состоянии находится одинаковое количество ресурса Ж/ё и для каждой вершины выполнено qi(0) < г°и‘. Вектор предельного состояния равен Ж01 , где 01 - любая строка матрицы А, определяемой по формуле (4).

Доказательство. Рассмотрим ё векторов начального состояния с единичным ресурсом, сосредоточенным в некоторой вершине каждого из циклических классов: 0т(0), где т - номер циклического класса, т = 1, ..., ё.

Обозначим их предельные векторы через 0* , где ] - номер

предельного вектора при заданном начальном состоянии. Получим следующие равенства:

0*т = 0т (0)Л0 Л7 = й(0) Л 7+т = 0*,

где 0* = 0* - предельные векторы при ресурсе, находящемся в

первом классе, индекс k е [1, ё] вычисляется по формуле (9) при т = 1.

Сложим векторы 0г(0) и найдем первый предельный вектор:

(й(0) +...+0ё (0)) Лё0 Л = 0* +...+0* = 0* +...+0*.

Для второго предельного вектора имеем

(01(0) +... + 0ё (0)) Лё0 л'2 = 0*2 +... + 0*2 = 0* +...+0* + 0;.

В общем случае т-й предельный вектор будет суммой век-

/"Л * /"Л *

торов 0\ + ... + 0ё с циклически переставленными слагаемыми.

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

Если ресурс в каждом классе пропорционально увеличить в Ж/ё раз, предельные векторы будут равны при условии, что все вершины сети функционируют по правилу 2.П

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

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

Т е о р е м а 1. В эргодической ё-циклической ресурсной сети, функционирующей по правилу 2, начиная с некоторого момента $ предельные векторы 0 ь ..., 0 ё совпадают тогда и только тогда, когда в каждом из циклических классов при $ = $ находится одинаковое количество ресурса Ж/ё. Вектор предельного состояния равен Ж01 , где 01 - любая строка матрицы А, определяемой по формуле (4).

Доказательство. 1. Пусть при $ = $ вся сеть функционирует по правилу 2 и в каждом циклическом классе находится одинаковое количество ресурса. Приняв 0(0 за новое начальное состояние, получим условие леммы 3.

2. Пусть в сети, функционирующей по правилу 2, векторы 0 ь ..., 0 ё совпадают. По определению для предельного цикла выполняются соотношения: 0 г+1 = 0 х Л', - = 1, ..., ё- 1,

0 1 = 0 аЛ'. Тогда имеем: 0 х = 0 х Л', - = 1, ..., ё . Умножение на матрицу Л переводит весь суммарный ресурс каждого класса в следующий циклический класс. Совпадение векторов происходит только при одинаковых суммарных ресурсах в каждом циклическом классе. □

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

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

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

0* = 01 = 02* = (0,5, 1, 1, 1, 0,5).

Рис. 3. Функционирование сети, заданной графом с матрицей (1), Q(0) = (4, 0, 0, 0, 0)

Таблица 2. Протокол работы сети

t VI *2 *3 *4 *5

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

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

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

3 1,000 1,500 0,500 1,000 0,000

4 0,750 1,500 0,500 1,000 0,250

5 0,750 1,250 0,750 1,000 0,250

20 0,501 1,002 0,998 1,000 0,499

21 0,501 1,001 0,999 1,000 0,499

22 0,500 1,000 1,000 1,000 0,500

23 0,500 1,000 1,000 1,000 0,500

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

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

В качестве потенциального приложения ресурсных сетей рассматриваются модели децентрализованного управления в многоагентных системах, в том числе задача нахождения консенсуса [1]. Эта задача может быть поставлена и решена в тех классах сетей, в которых предельное состояние существует и зависит от начального. Задача управления распределением ресурса, и, в частности, достижения консенсуса - объект дальнейшего исследования ресурсных сетей.

Литература

1. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Сходимость и устойчивость в задачах согласования характеристик (обзор базовых результатов) // Управление большими системами. -

2010. - Специальный выпуск 30.1 «Сетевые модели в управлении». - С. 470-505.

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

3. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. I. Процессы стабилизации при малых ресурсах // Автоматика и телемеханика. - 2011. - №4. - С. 133-143.

4. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. II. Потоки при больших ресурсах и их стабилизация // Автоматика и телемеханика. - 2012. - №6. - С. 103-118.

5. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. III. Исследование предельных состояний // Автоматика и телемеханика. - 2012. - №7. - С. 67-77.

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

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

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

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. - 353 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.

I. OSCILLATIONS AND EQUILIBRIUM AT LOW RESOURCE

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

Abstract: We study operation of ergodic nonregular resource networks when total amount of resource is small. We show that for an arbitrary initial state undamped oscillations occur between several limit vectors. The number of different limit vectors is equal to the number of cyclic classes of a network. The global equilibrium can be achieved when all limit vectors coincide. We derive conditions on initial states leading to the global equilibrium, and obtain formulae for the limit powers of stochastic matrices and for limit vectors.

Keywords: resource network, stochastic matrix, ergodic Markov chain, equilibrium, limit state.

Статья представлена к публикации членом редакционной коллегии М.В. Губко

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

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