CC BY
49
УДК 519.177 ББК 22.18
РАСПРЕДЕЛЕНИЕ РЕСУРСА МЕЖДУ АТТРАКТОРАМИ В РЕГУЛЯРНЫХ НЕСИММЕТРИЧНЫХ РЕСУРСНЫХ СЕТЯХ1
Жилякова Л. Ю.2
(ФГБУН Институт проблем управления им. В.А. Трапезникова РАН, Москва)
Рассматриваются несимметричные регулярные ресурсные сети с несколькими вершинами-аттракторами. Для них показано, что распределение ресурса сверх порогового значения АЖ = Ж - Т происходит по тому же закону, что и в соответствующей поглощающей сети (полученной из несимметричной удалением выходных ребер аттракторов). Однако существуют поправки, зависящие от пропускных способностей выходных ребер аттракторов и начального распределения ресурса. Поправки оцениваются сверху; находятся начальные состояния, при которых они максимальны. Определяются начальные состояния, при которых поправки отсутствуют.
Ключевые слова: ресурсная сеть, графовая динамическая пороговая модель, вершины-аттракторы.
1. Введение
Ресурсная сеть - графовая динамическая параллельная модель с пороговым переключением правил функционирования, предложенная в [6] и получившая развитие в ряде работ (в том числе [3, 4]). Она представляет собой ориентированный взве-
1 Работа выполнена при частичной финансовой поддержке РФФИ (ппроекты № 14-01-00422а, 15-07-02488а).
2 Людмила Юрьевна Жилякова, доктор физико-математических наук (zhilyakova. ludmila@gmail. com).
шенный граф, в котором вершины обмениваются ресурсами в дискретном времени. Вершины могут хранить неограниченное количество ресурса; веса ребер обозначают их пропускные способности: максимальное количество ресурса, которое может пройти через ребро за один такт времени. В отличие от классической потоковой модели Форда-Фалкерсона [7] и ее динамических модификаций [8], в ресурсной сети отсутствуют источники и стоки, и, соответственно, нет направленности движения ресурса. Его распространение происходит по всем возможным путям одновременно. При малом ресурсе, когда все вершины сети функционируют по одному правилу и отдают весь ресурс на каждом такте, модель представляет собой рассеяние на графе [12, 14]. Модели рассеяния хорошо изучены, - для графов с постоянной топологией они описываются однородными цепями Маркова, - и поэтому ресурсная сеть представляет гораздо больший исследовательский интерес при функционировании с большим суммарным ресурсом. Это означает, что некоторые вершины сети изменяют правило функционирования и накапливают излишки ресурса, которые они не в состоянии отдать. Такие «бутылочные горлышки» сети, суммарные входы которых больше суммарных выходов, были названы вершинами-аттракторами [3]. Эти вершины начинают функционировать сходно со стреляющими вершинами в пороговой модели «игра выстреливания фишек» [11, 18] и «абелева куча песка» [15, 20]. Модель кучи песка примечательна тем, что при своей простоте находит применение в моделировании сложных природных и социальных процессов, названных «самоорганизованной критичностью» [9, 10, 13].
В отличие от игры выстреливания, где фишки целочислен-ны, ресурс в ресурсной сети является бесконечно делимым. Благодаря этому свойству разнообразные переходные процессы сходятся к предельным состояниям, которые в целочисленных пороговых моделях могут не достигаться.
Для каждой сети понятия малый и большой ресурс строго определены. Малое и большое значения ресурса лежат по разные стороны от порогового значения Т, которое определяется через суммарные выходные пропускные способности вершин и
координаты вектора предельных вероятностей однородной цепи Маркова, соответствующей ресурсной сети, в случае если граф регулярен. Свойство регулярности [5] означает, что граф сильно связен (т.е. состоит из одной эргодической компоненты) и НОД длин всех его циклов равен 1 (заметим, что в разных источниках регулярность определяется по-разному и данное определение отличается от определений, используемых в [2, 19]). Если граф циклический (НОД длин всех его циклов отличен от 1), роль предельного вероятностного вектора играет предел по Чеза-ро [5]. В общем случае пороговое значение не превосходит суммарной пропускной способности сети. В эйлеровых сетях, в которых каждая вершина имеет одинаковые входную и выходную пропускные способности, Т совпадает с суммой пропускных способностей всех ребер сети. В остальных сетях выполняется строгое неравенство.
Задача исследования в данной работе формулируется следующим образом. Пусть в сети есть несколько аттракторов и суммарный ресурс в ней велик. Доказано, что предельное состояние в ней существует. Компоненты вектора предельного состояния, соответствующие неаттрактивным вершинам, определяются однозначно; формулы для них найдены в [3]. Необходимо найти распределение излишков ресурса сверх Т между аттракторами. Пропорция, в которой эти излишки распределяются, очевидно, зависит от начального состояния и от некоторых параметров сети.
Функционирование сетей при большом ресурсе описывается неоднородными цепями Маркова, исследование которых гораздо сложнее. Из существования предельного состояния при больших ресурсах [3] следует, что эти неоднородные цепи Маркова являются сильно эргодическими [1, 16, 17]. Однако для нахождения компонент вектора предельного состояния, соответствующих вершинам-аттракторам, требуются дополнительные приемы. Для этого несимметричным сетям ставятся в соответствие поглощающие сети, исследованные в [4]. Граф поглощающей сети имеет переходную компоненту и стоковые вершины. Порогового значения в такой сети не существует.
Поглощающая сеть, соответствующая регулярной несимметричной сети, - это сеть, в которой удалены все ребра, выходящие из аттракторов. Таким образом, это сеть со стоковыми вершинами, соответствующими аттракторам исходной сети. Для поглощающих сетей получены формулы, позволяющие аналитически находить предельные стохастические матрицы неоднородной цепи Маркова, и, соответственно, по любому начальному состоянию определять предельное состояние сети - в данном случае пропорцию распределения ресурса между стоками.
2. Основные обозначения и характеристики ресурсных сетей
Ресурсная сеть представляет собой взвешенный ориентированный граф G = (V, E), |V = n. Веса, соответствующие пропускным способностям ребер, задаются матрицей R = (rij)nxn.
Q(t) = (q\(f), ..., qn(t)) - состояние сети в момент t; qi(t) - количество ресурса в вершине vi в момент t; Q* - предельное состояние.
r'n и rCMt - суммарные входная и выходная пропускные способности вершины vi.
Правила функционирования сети. В момент t вершина vi по ребру eim отдает в смежную ей вершину vm:
• rim единиц ресурса, если qi(t) > r°ut (правило 1);
r
• -mq, (t) единиц ресурса, если q (t) < rout (правило 2).
rout
R' - стохастическая матрица, полученная из матрицы R нормированием по строкам: R = D lR, где D = diag(r1out,...,r"0ut).
W = const - суммарный ресурс сети.
Z(t) - множество вершин, для которых q{ (t) < r°ut; Z+(t) - множество вершин, для которых q{ (t) > r°ut.
Т - пороговое значение ресурса, такое, что при W < Т все вершины, начиная с некоторого t', переходят в зону Z-(t); при W > T зона Z+(t) непуста, начиная с некоторого t".
T = min ,
ie{1,...,n} q.
где Q1*- вектор предельного состояния при W = 1; Q - вектор предельного состояния при W = Т.
Вершина Vj несимметричной регулярной сети является аттрактором тогда и только тогда, когда для нее выполняется
r°ut
(1) j = arg mn .
¿e{1,--n} q1
Введем такую нумерацию вершин сети, при которой аттракторы имеют номера от 1 до l. В [3] доказано, что
Q = (-°ut,..., -Г
, ql+1,...,qn) . Вектор Q пропорционален вектору
Q1* с коэффициентом Т: Q = QQ* T .
При этом компоненты с номерами от l + 1 до n не изменяются при любом W > T, и все излишки распределяются только среди аттракторов. Вектор предельного состояния для W> T имеет вид
i
Q = (-°ut + Аql..,-r + Aq*,ql+l,..,qn) , где £Aq* = W- T.
i =1
Поток в ресурсной сети.
Ресурс, выходящий из вершины Vi по ребру ej в момент t, приходит в вершину Vj в момент t + 1; между моментами t и t + 1 он находится в ребре ej. Этот ресурс назовем потоком fj(t). Общий поток сети описывается матрицей F(t) = (f(O)nxn. Величиной потока будем называть сумму
n n
fsum
= H fj (t).
i=1 j=1
Суммарные исходящий и входящий потоки вершины Vi:
n n
£fj (t) = f°ut(t); (t) = fin(t +1). Будем полагать,
j=1 i=1
что j (0) = 0 .
F°"'(t) = (f°ut (t), ..., C' (t)) - вектор исходящего потока;
^"(0 = (/Г(1) , •••, /ПП(0 ) - вектор входящего потока. Р - матрица предельного потока.
п п
/1т = XX/у - суммарный предельный поток; Р"* и Р}иГ* -i=\ ]=\
векторы двух предельных потоков.
Доказано, что при любом значении Ж> Т суммарный предельный поток в сети равен Т.
3. Свойства поглощающих сетей
Исследованию поглощающих сетей посвящена работа [4]. Кратко перечислим полученные в ней результаты, которые понадобятся при описании регулярных сетей с несколькими аттракторами.
Пусть поглощающая сеть имеет I стоков с номерами от 1 до I, тогда ее матрица пропускных способностей состоит из блоков:
(2) Я =
где (О) х I - диагональная матрица; (О\)г х"-1) - нулевая матрица; матрица (Яг)"-I) х I состоит из пропускных способностей ребер, ведущих из переходной компоненты в стоки; матрица (Я2)" -1) х (" - I) - из пропускных способностей ребер, соединяющих вершины внутри переходной компоненты.
Стохастическая матрица Я' будет иметь вид
Я'=
В 0\
Я\ Я2
Г Е 0\ ^
1 Я\ Я2J
Даже если матрица О имеет на диагонали нулевые элементы (т.е. одна или несколько стоковых вершин не имеют петель), блок Е\ будет единичной матрицей.
Свойство 1 (Теорема 2 [4]). В поглощающей сети порогового значения ресурса Т не существует.
Свойство 2 (Лемма 1 [4]). Пусть поглощающая сеть с I стоками представлена матрицей (2). Тогда - предел степеней стохастической матрицы Я - определяется по формуле
(3) Д'ш =
Е
(Е2 - )
О 1
О
2 /
где Е2 - единичная матрица размера (п - I) х (п - I).
Свойство 3 (Теорема 3 [4]). В поглощающей сети с I стоками матрица остается неизменной при любых изменениях диагональных элементов.
Свойство 4 (Теорема 4 [4]). Пусть поглощающая сеть с I стоками представлена матрицей (2). Тогда для любого ресурса Ж и любого начального состояния Q(0) = (#1(0), ..., #п(0)) предельное состояние рассчитывается по формуле ( = ((0)Д'Ш,
где - предельная матрица, определяемая по формуле (3).
Следствие. Вектор предельного состояния (* не зависит от наличия или отсутствия петель в вершинах.
Свойство 5 (Теорема 5 [4]). Пусть поглощающая сеть с I стоками представлена матрицей (2). Тогда элементы ¡-й строки матрицы (/' > I) соответственно равны компонентам вектора предельного состояния при начальном состоянии (¡(0) = (0, ..., 0, 1, 0, ..., 0), где единице равна ¡-я компонента: ( „т Л
(4) Д'ш =
а
1+1
в*.
(е\, ..., ег - первые I вектор-столбцов единичной матрицы
(Е)ПХП).
т
е
4. Распределение ресурса между аттракторами в несимметричной регулярной сети
Поведение регулярных несимметричных сетей с несколькими аттракторами сложнее поведения поглощающих сетей. Так, в поглощающей сети ресурс, помещенный в сток в начальном состоянии, остается в нем в любой момент времени. В регулярных сетях это условие не выполняется. Для иллюстрации этого утверждения рассмотрим пример функционирования сети с двумя потенциальными аттракторами.
Пример 1. Несимметричная сеть (" = 5) с двумя аттракторами VI и V2 задана матрицей
(5) Я =
Основные ее характеристики следующие.
• Кортеж р = ((г/п;г°и'), ..., (г^п;г™')) для этой сети имеет
вид р = ((10, 5), (10, 5), (5, 10), (5, 10), (5, 5)).
• Предельное состояние при Ж = 1: 01* = (0,250, 0,250, 0,1(6), 0,1(6), 1(6)).
• Пороговое значение ресурса Т = 20.
• 0 = (5, 5, 3,(3), 3,(3), 3,(3)).
Рассмотрим начальное состояние 0(0) = (0, 5, 20, 0, 0). Суммарный ресурс 25 больше порогового значения, значит, излишки ресурса должны распределиться между аттракторами. Аттрактор V2 имеет ресурс в начальном состоянии, аттрактор V: - нет.
Предельное состояние для такого вектора 0(0) будет 0* = (10, 5, 3,(3), 3,(3), 3,(3)). Как видно, весь излишек оказался в первой вершине (поскольку она связана с источником, имеющим большой ресурс, ребром с большой пропускной способно-
Г \ \ \ \ Л
\ \ \ \ \
6 \ \ \ \
\ 6 \ \ \
V \ ! ! \ V
стью). Вершине V2 не удалось накопить больше того, что она имела в начале.
Даже с ресурсом, превышающим выходную пропускную способность, потенциальный аттрактор бывает неспособным перейти в зону 2+*. Для той же сети перераспределим ресурс 25 в начальном состоянии: 0(0) = (0, 8, 17, 0, 0). Предельное состояние в этом случае будет точно таким же: 00 = (10, 5, 3,(3), 3,(3), 3,(3)).
Более того, если ресурс 8 переложить в вершину V4 (источник, связанный ребром с большой пропускной способностью с V2), так, чтобы начальное состояние было 0(0) = (0, 0, 17, 8, 0), то предельное состояние будет: 0 = (9,75, 5,25, 3,(3), 3,(3), 3,(3)). То есть в этом случае второй аттрактор получает излишки сверх своей пропускной способности, чего не наблюдалось, когда мы помещали ресурс непосредственно в него.
Таким образом, потенциальный аттрактор может притягивать ресурс из источников, но не всегда может удержать ресурс, находящийся в нем в начальном состоянии.
В данном исследовании будем рассматривать начальные состояния, в которых ресурс находится в источниках.
4.1. ПОГЛОЩАЮЩАЯ СЕТЬ, СООТВЕТСТВУЮЩАЯ НЕСИММЕТРИЧНОЙ, И ПОПРАВКА НА РЕГУЛЯРНОСТЬ 5Ж
Исследуем распределение ресурса между потенциальными аттракторами при Ж > Т. В [3] получена формула предельного состояния при больших ресурсах. В ней определены все компоненты вектора предельного состояния, кроме компонент, соответствующих вершинам-аттракторам. Если аттрактор в сети один, предельное состояние единственно. Если потенциальных аттракторов несколько, ресурс между ними в предельном состоянии может распределяться по-разному для разных начальных состояний. Опишем зависимость предельного состояния в аттракторах от начального состояния сети, как это было сделано для поглощающих сетей.
Рассмотрим несимметричную сеть с двумя аттракторами, заданную матрицей (5).
Пусть в начальном состоянии весь суммарный ресурс Ж > T находится в одном из источников. Найдем закономерность распределения ресурса Ж - Т между аттракторами при изменении Ж Рассмотрим два начальных состояния и соответствующие им предельные состояния:
(6) ((0) = (0, 0, 50, 0, 0) ^ ( = (29,3, 10,7, 3,(3), 3,(3), 3,(3)).
(7) ((0) = (0, 0, 90, 0, 0) ^ ( = (59,3, 20,7, 3,(3), 3,(3), 3,(3)).
Суммарный ресурс вершин, находящихся в зоне в предельном состоянии равен 10. При ((0) = (0, 0, 50, 0, 0) и ((0) = (0, 0, 90, 0, 0) в зоне Т+* распределяется ресурс, равный 40 и 80 соответственно. При суммарном приросте 40 прирост в вершинах составил соответственно 30 и 10, т.е. 0,75ДЖ и 0,25 Д Ж Прирост ресурса распределился в пропорции 3:1.
При любом дальнейшем увеличении ресурса пропорция 3:1 в распределении ресурса ДЖ между аттракторами сохраняется.
В силу симметричности аттракторов, если в начальном состоянии весь ресурс находился во втором источнике (вершине в предельном состоянии его прирост в аттракторах распределится в пропорции 0,25ДЖ и 0,75ДЖ (соответственно, 1:3). А при помещении ресурса в начальном состоянии в нейтральную вершину прирост в аттракторах распределится как 0,5 ДЖ, 0,5 ДЖ, т.е. в пропорции 1:1.
Для объяснения этого факта построим поглощающую сеть, соответствующую данной несимметричной сети. То есть превратим аттракторы в стоки, заменив все элементы матрицы Я, соответствующие входным пропускным способностям аттракторов, нулями. Матрицу пропускной способности такой сети обозначим через Яаь$огь:
Я,
аЬзогЬ
(0 0 0 0 01
0 0 0 0 0
6 1 1 1 1
1 6 1 1 1
V 1 1 1 1 1,
Предел степеней соответствующей стохастической матрицы обозначим через К^ь^ь. Найдем эту матрицу, воспользовавшись формулой (4). Для этого вычислим векторы предельного состояния при единичном ресурсе, поочередно находящемся в каждой из переходных вершин:
0(0) = (0, 0, 1, 0, 0) ^ 0 = (0,75, 0,25, 0, 0, 0); 0(0) = (0, 0, 0, 1, 0) ^ 0 = (0,25, 0,75, 0, 0, 0); 0(0) = (0, 0, 0, 0, 1) ^ 0 = (0,5, 0,5, 0, 0, 0). Тогда
(8)
R »
R absorb
1 0 0 0 01
0 1 0 0 0
0,75 0,25 0 0 0
0,25 0,75 0 0 0
0,5 0,5 0 0 0,
Можно заметить, что именно в таких пропорциях распределяется ресурс AW из каждой не-аттрактивной вершины между аттракторами.
Возникает вопрос, зависит ли эта пропорция от первых двух строк матрицы R или выходные пропускные способности аттракторов могут быть любыми до тех пор, пока они остаются аттракторами. При изменении выходных пропускных способностей аттракторов меняется пороговое значение Т. Остается ли при больших W пропорция, заданная Rr^bsorb. ресурса?
Рассмотрим следующие примеры.
Пример 2. Сеть задана матрицей 0,5 0,51
для излишков
(9) R =
f 0,5 0,5 0,5
0,5 0,5 0,5
6 1 1
1 6 1
. 1 1 1
0,5 1 1 1
0,5 1 1 1
Для этой сети Т = 10; Q1* = (0,250, 0,250, 0,1(6), 0,1(6), 1(6));
V
в =(2,5, 2,5, 1,(6), 1,(6), 1,(6)).
Функционирование сети при разных начальных состояниях дает следующие предельные состояния:
((0) = (0, 0, 85, 0, 0) ^ = (60,5, 19,5, 1,(6), 1,(6), 1,(6));
((0) = (0, 0, 45, 0, 0) ^ ( = (30,5, 9,5, 1,(6), 1,(6), 1,(6)).
Ресурс, равный 5, остается в зоне Тт*, остальной распределяется между аттракторами. Разность ресурса для этих двух начальных состояний ДЖ = 40 разделилась в пропорции 30:10, т.е. 3:1, как и раньше.
Это соотношение выполняется не только для однородных выходных пропускных способностей аттракторов. Но чтобы обе вершины оставались потенциальными аттракторами, их выходные пропускные способности нельзя изменять произвольно. Не приводят к потере аттрактивности любые изменения в пропускных способностях петель аттракторов, а также одинаковые изменения в пропускных способностях ребер, связывающих аттракторы с нейтральной вершиной \5 и друг с другом.
Пример 3. Сеть задана матрицей
Я =
Для нее имеем: Т = 23;
= (0,261, 0,478, 0,116, 0,116, 0,029); в = (6, 11, 2,(6), 2,(6), 0,(6)).
Как видно, вектор предельных вероятностей изменился; в той же пропорции изменились компоненты вектора в . В зоне Т * в предельном состоянии оказывается ресурс, равный 6.
Рассмотрим две пары векторов (начальное состояние ^ предельное состояние): ((0) ^ (*.
((0) = (0, 0, 86, 0, 0) ^ (* = (56,910, 23,090, 2,(6), 2,(6), 0,(6)). ((0) = (0, 0, 46, 0, 0) ^ = (26,910, 13,090, 2,(6), 2,(6), 0,(6)).
(3 11 1 01
1 8 1 0
6 11 1 1
1 61 1 1
V 1 11 1 1,
Разность ресурса AW = 40 вновь разделилась в пропорции 3:1, несмотря на вышеописанные изменения векторов предельного состояния.
Эта пропорция имеет место только для разности между двумя различными векторами. В пределах одного вектора предельного состояния она не выполняется. Рассмотрим первую пару векторов:
Q(0) = (0, 0, 86, 0, 0) ^ Q* = (56,910, 23,090, 2,(6), 2,(6), 0,(6)).
rout = 6, r°* = 11.
Aqx = q* - r°out = 50,910, Aq2 = q* - r2°ut = 12,090.
Aq : Aq = 50,910 : 12,090 > 4.
Вернемся к матрице (5) и проверим выполнение этой пропорции для нее. Из примера 1 имеем:
Q(0) = (0, 0, 100, 0, 0) ^ Q*= (66,8, 23,2, 3,(3), 3,(3), 3,(3)).
ТТ " „out г „out г
Для этой матрицы r = 5, r = 5.
Aq = q* - r°ut = 61.8, Aq = q* - rf = 18,2.
Как видно, пропорция 3:1 вновь не выполняется.
Вычислим разность Aqx - 3 Aq2:
Aq -3Aq2 = 61,8 - 18,2 • 3 = 61,8 - 54,6 = 7,2.
Проверим, чему равна эта разность при других значениях W.
Из (6), (7) следует, что при W = 50 и W = 90 также выполняется равенство (10) Aq - 3Aq2 = 7,2.
Для достаточно больших значений W в сети, заданной матрицей (5), значение разности (10) остается неизменным. Обозначим его через SW. Начиная с некоторого Wmin > T, SW = const.
Таким образом, для поглощающей сети 5W = 0. Как только поглощающая сеть превращается в регулярную, (т.е. у стоков появляются выходные ребра), величина 5W становится положительной. Будем называть SWпоправкой на регулярность.
Опишем процесс нахождения 5W для произвольной сети с l аттракторами и произвольного начального состояния Q(0), в котором ресурс, больший Wmin, сосредоточен в не-аттрактивных
вершинах. Пусть Я - матрица пропускных способностей сети с аттракторами, имеющими номера от 1 до I. Пусть (* - вектор предельного состояния при суммарном ресурсе Ж > 2гшт (количество ресурса 2гШт заведомо больше Жтт). Построим вектор Д( = С - в . Его последние п - I компонент равны нулю, а первые I компонент содержат ресурс, накопленный аттракторами сверх их пропускных способностей. Я,ьхогь - матрица поглощающей сети, первые I строк которой равны нулю, а остальные п - I строк совпадают со строками матрицы Я. Рассмотрим вектор предельного состояния этой поглощающей сети как функцию от количества ресурса ( = (*(Ж). Будем увеличивать Ж от нуля до тех пор, пока какая-нибудь компонента не совпадет с компонентой вектора Д(. Обозначим это количество ресурса через Ж. Тогда все компоненты разности Д( - (*(Ж) неотрицательны по построению. Их сумма и составляет поправку на регулярность 5Ж. Это такое количество ресурса, которое нарушает пропорциональность, имеющую место для поглощающей сети, распределяясь между аттракторами. Эта величина заведомо положительна. Используя введенные обозначения, можно дать формальное определение.
Поправкой на регулярность ЗЖ называется суммарное количество ресурса в предельном состоянии в одном или нескольких аттракторах, нарушающего пропорцию, заданную матрицей
ЯaЬsorь'■ ЗЖ = ((* - в - (*(Ж))-1, где 1 - вектор-столбец из п единиц.
Для того чтобы найти Жтт и объяснить, откуда возникает вклад ресурса ЗЖ, нарушающий пропорцию, имеющую место для стоков, рассмотрим функционирование сети при ресурсах Ж > Т, близких к пороговому значению, т.е. таких, что Ж - Т < ЗЖ.
Пример 4. Сеть задана матрицей (5). Начиная с Ж = Т, будем увеличивать суммарный ресурс в сети, помещая его в начальном состоянии в третью вершину. Для этой сети Т = 20. Выпишем начальные и соответствующие им предельные состояния:
0(0) = (0, 0, 20, 0, 0) ^ Q = (5, 5, 3,(3), 3,(3), 3,(3)). 0(0) = (0, 0, 25, 0, 0) ^ Q = (10, 5, 3,(3), 3,(3), 3,(3)). 0(0) = (0, 0, 27, 0, 0) ^ Q = (12, 5, 3,(3), 3,(3), 3,(3)). 0(0) = (0, 0, 27,2, 0, 0) ^ Q = (12,2, 5, 3,(3), 3,(3), 3,(3)). 0(0) = (0, 0, 27,21, 0, 0) ^
^ Q = (12,2075, 5,0025, 3,(3), 3,(3), 3,(3)). Из этих пар векторов видно, что при значениях W- T < 7,2 ресурс из источника не распределяется пропорционально, а накапливается в одном аттракторе. Ресурс второго аттрактора остается в предельном состоянии равным r= 5 вплоть до
W = 27,2, и, соответственно, SW = 7,2.
При суммарном ресурсе W > 27,2 второй аттрактор тоже попадает в зону Z+*, и начинает выполняться соотношение (q* -5-7,2): (q* -5) = 1:3
В общем виде для сети c двумя симметричными аттракторами V1 и V2, запитанными на источники V3 и V4 соответственно, можно утверждать следующее. Если ресурс в начальном состоянии находится в вершине V3 и при этом его величина удовлетворяет неравенству W > T + SW, в предельном состоянии выполнится соотношение
(q* - r- - SW) : (q* - ^) = ralorb 31 : ralorb32, где матрица R'absorb задана формулой (8).
То есть, если известно значение SW, предельное состояние при заданном начальном находится однозначно.
Лемма 1. Пусть в регулярной несимметричной сети с двумя аттракторами v\ и V2, ресурс W> T в начальном состоянии
Л Л 'со 'со
находится в источнике v¡, для которого r ъ ъ , = max rъ ь ., .
J ke{3,...n}
Тогда вектор предельного состояния равен
(П)
. \Ü + (W-T, О, О, ..., 0) при W<T + SW, Q ={Q + {SW + ralorbjl^W, r:ZorbJ2-AW,0,...,0) при W>T + SW;
где АЖ = Ж - (Т + ЗЖ); - предельная матрица погло-
щающей сети, соответствующей данной регулярной сети.
Д о к а з а т е л ь с т в о. Дисбаланс в аттракторах происходит по следующей причине. Пока аттрактор V2 добирает свой ресурс до выходной пропускной способности, он функционирует по правилу 2, т.е. на каждом такте отдает весь свой ресурс. Аттрактор VI переходит на функционирование по правилу 1 и начинает отдавать по полной выходной пропускной способности, получая больше этого значения. Излишек на каждом такте остается в V!. Таким образом, ЗЖ состоит из ресурса, накопленного одним аттрактором, пока он работает по правилу 1, в то время как другой аттрактор работает по правилу 2. Если значение ресурса суммарного ресурса Ж < Т + ЗЖ, то по определению ЗЖ ресурс аттрактора V2 асимптотически стремится к г^"', и V2 не переходит в зону £+({), и верхняя строка (11) доказана.
Если Ж > Т + ЗЖ, оба аттрактора оказываются в зоне Z+(t). В некоторый момент времени выполнится: - 1) < г^"' и > г. Тогда д^') = г°"' + ЗЖ . Покажем, что при переходе второго аттрактора в зону ресурс начинает распределяться между аттракторами в пропорции, заданной Я'аЬтгЬ. Выходной поток в аттракторах при их переходе в зону стабилизировался: он равен г""' в первом аттракторе и г""' во втором - в каждое ребро они отдают по полной пропускной способности. Накопление в аттракторах происходит, пока их входной поток больше выходного. Оно складывается из разницы стабильной и нестабильной компонент, поскольку в сбалансированном (предельном) потоке для аттракторов выполняется: у,п* = уош* = гО"', ут* = уош* = г. Если рассматривать только нестабильную компоненту потока, ей соответствует матрица, в которой выходные пропускные способности аттракторов равны нулю, а это и есть матрица ЯаЬтгЬ.
Таким образом, ресурс сверх Т + ЗЖ распределяется между аттракторами в пропорции, заданной матрицей . □
Лемма легко обобщается на случай произвольного количества аттракторов. Пусть в сети с I аттракторами Vl, ..., V;, ресурс
Q =
W > T в начальном состоянии весь сосредоточен в источнике Vj, j > l. Обозначим через SWi - поправки в каждом аттракторе при любом W > T - достаточно большом, чтобы часть ресурса распределялась в предельном состоянии в соответствии с матрицей
Rabsorb. Ясно, что 5Wi = const, поскольку после перехода последнего из аттракторов в зону Z+(t) выходной поток в них стабилизируется, и ресурс начинает прибывать в них с пропорции, заданной матрицей Я'^огь •
Теорема 1. Пусть в регулярной несимметричной сети с l аттракторами vi, ..., vi ресурс W> T в начальном состоянии находится в источнике Vj, j > l. Тогда вектор предельного состояния равен: (12)
б + (0,..., W-T, 0,..., 0) при W<T + (SWm-SW_m), & + (ЩЛ,...,ЩЛ, 0,...,0) при T + (SWm-SW_m)<W<T + SW, в + (Щ +r'2orbjl-AW, ...,SWl+r2orbjl-AW, 0,...,0)
прuW > T + SW;
где вектор (0, W- T, 0, 0) имеет ненулевую компоненту с номером m < l; m находится из условия m = arg max ra^b km;
ke{1,...l}
0 < SWiЛ < SWi - поправки при T + (SWm - SW-m) < W< T + SW;
SWm = max SWk - поправка, соответствующая аттрак-
ke{1,..J}
тору Vm;
SW = max SW, - вторая по величине поправка
Ще{Щ,...Щ }\SWm
при ресурсе, сконцентрированном в источнике Vj;
i
SW = YJSW1 ; AW = W-(T + SW);
i=1
причем существует аттрактор Vi, для которого SWi = 0, i < l. Д о к а з а т е л ь с т в о.
Докажем каждую строку из формулы (12) отдельно. 1. Справедливость формулы для случая
W < T + (SWm - SW-m) вытекает из доказательства леммы 1. Аттрактор Vm первым переходит на правило 1 и начинает наби-
рать излишки ресурса, составляющие поправку ЗЖт" = (Ж - Т). Если второй (по величине ЗЖ,) аттрактор не может перейти на правило 1, все излишки ресурса в предельном состоянии окажутся в вершине ут. Переход второго аттрактора на правило 1 происходит, когда аттрактор Ут имеет ресурс, равный (ЗЖт - ЗЖ-т). Отсюда непосредственно вытекает условие выполнения равенства ( = (% + (0,..., Ж - Т, 0,..., 0).
2. При ресурсе, превосходящем значение Т менее чем на суммарную поправку ЗЖ, при функционировании сети часть аттракторов может не перейти в зону 2+(0, соответствующие поправки ЗЖ,Л = ЗЖг'(Ж) для них будут равны нулю, и компоненты предельного состояния будут равны соответствующим компонентам вектора ( = (гЦ"*,...,,%+1(0), ..., %(0)). По мере увеличения суммарного ресурса сверх величины Т + (ЗЖт - ЗЖ-т) аттракторы один за другим смогут перейти в зону 7+{(). Если ввести такую нумерацию аттракторов, что они будут упорядочены по убыванию величины ЗЖ,, т.е. ЗЖ\ := ЗЖт, ЗЖ2 := ЗЖ-т, и т.д., то можно вычислить значение ресурса Ж, при котором аттрактор с номером к будет иметь в предельном состоянии ЗЖк(Ж) > 0. Для этого необходимо, чтобы все аттракторы с номерами, меньшими к, тоже имели положительные поправки. Разобьем множество значений суммарного ресурса Ж > Т на интервалы и рассмотрим зону для каждого из них.
1) Т< Ж< Т + (5Ж\ - ЗЖ2): только аттрактор У\ функционирует по правилу 1. Введем обозначение: Т + (ЗЖ\ - ЗЖ2) = Т\. Далее границы полуинтервалов будем обозначать через Т,.
2) Т\ < Ж< Т + (ЗЖ\ - ЗЖ2) + 2(ЗЖ - ЗЖз) = Т + ЗЖ\ + ЗЖ2 -- 2ЗЖз = Т2: аттракторы У\ и У2 функционируют по правилу \
3) Т2 < Ж < Т + (ЗЖ\ - ЗЖ2) + 2(ЗЖ - ЗЖз) + 3(СЖз - ЗЖ4) = = Т + ЗЖ\ + ЗЖ2 + ЗЖз - ЗЗЖ4 = Тз: аттракторы У\, У2 и Уз функционируют по правилу \.
I - 2) Т-з < Ж< Т + ЗЖ\ + ... + ЗЖ-2 - (I - \)ЗЖ-\ = Т-2: аттракторы V!, У2, ..., У1-\ функционируют по правилу \
l - 1) Ti-2 < W< T + SWi + ... + SWi-V. аттракторы vi, ..., v-i функционируют по правилу 1, но поправка в аттракторе v-i еще не достигла своего максимального значения.
i) W> T + SWi + ... + SWi-i = T + SW все аттракторы переходят в зону Z+(t), и поправки достигли максимальных значений. При этом, в аттракторе vi поправки нет: SWi = 0, так как больше нет аттракторов, которым нужно накапливать ресурс.
3. Третья строка формулы (i2) доказана в пункте i). При W > T + SW все аттракторы переходят в зону Z+(t), и выходной поток в них стабилизируется за конечное число тактов. Далее в сети начинается перераспределение ресурса в соответствии с
матриЦей Reabsorb . D
Замечание i. Если вершина-аттрактор, для которой выпол-
'сО
няется условие m = arg max rabsorb km, не единственна, то ресурс
W - T в верхней строке формулы (12) распределится в равных долях между всеми этими аттракторами.
Замечание 2. Проиллюстрируем графически формулы в пунктах 1) ^ i).
На рис. 1 представим, что пунктирные линии - смещающаяся вниз ось абсцисс. Они задают уровень, ресурс над которым виден. Чем ось ниже, тем больше вершин-аттракторов имеют положительные поправки. При непрерывной («главной») оси абсцисс все вершины имеют максимальные поправки. При дальнейшем увеличении ресурса излишки будут распределяться в пропорции, заданной матрицей R'^bsorb.
Как и во всех гистограммах количество ресурса рассчитывается как площадь частей столбиков, выделенных тем или иным цветом. Ширину каждого столбца считаем единичной. На рис. 2 количество ресурса над соответствующей осью показано цветом. Этот рисунок иллюстрирует возникновение каждого слагаемого в суммах из выражений 1) ^ i).
Рис. 1. Расположение величин ЗЖ, в порядке убывания
Рис. 2. Для перехода аттрактора в зону 7+($), ресурс должен превысить сумму ресурсов над пунктирной чертой
Следствие. Непосредственно из доказательства теоремы вытекает, что если в несимметричной сети имеется I потенци-
альных аттракторов, а ресурс в начальном состоянии сосредоточен в вершине-источнике V], ] > I, то значения суммарного ресурса при Ж > Т можно разбить на I интервалов, для каждого из которых поправки ЗЗЖ, будут иметь 7 вершин (7 = 1, ..., /). Границы интервалов задаются формулами из пп. 1) ^ I).
При ресурсе Ж > Т + ЗЖ все вершины имеют максимально возможную поправку, ресурс сверх Т + ЗЖ распределяется в соответствии с матрицей Я'^огЬ.
Значения ЗЖ, можно найти экспериментально для любой сети. Чтобы получить аналитические зависимости и выявить закономерности, рассмотрим изменение ЗЖ, при изменении строк матрицы пропускных способностей, соответствующих аттракторам. Если все элементы в них равны нулю, сеть превращается в поглощающую.
4.2. ЗАВИСИМОСТЬ Ш ОТ ВЫХОДНЫХПРОПУСКНЫХ СПОСОБНОСТЕЙ АТТРАКТОРОВ
Вновь рассмотрим сеть с двумя потенциальными аттракторами. Здесь поправка ЗЖ накапливается только в одной вершине, поэтому можно записать ЗЖ1 = ЗЖ и далее оперировать величиной ЗЖ.
Проследим зависимость ЗЖ от выходных пропускных способностей вершин-аттракторов: зж = зж г*), 7 = 1, 2, считая, что первые две строки матрицы Я однородны, как в матрицах (5) и (9).
Пример 5. В матрице ^ г г г г г^ г г г г г Я =61111 16 111 х 11111 ,
будем изменять пропускные способности г от нуля с некоторым шагом. Данные экспериментов представлены в таблице 1.
Таблица 1. Характеристики предельного состояния
при разных выходных пропускных способностях аттракторов
№ r yGUt yGUt T SW
1. 0 0 0 0
2. 0,05 0,25 1 (0,25, 0,25, 0,1(6), 0,1(6), 0,1(6)) 0,875
3. 0,1 0,5 2 (0,5, 0,5, 0,(3), 0,(3), 0,(3)) 1,75
4. 0,15 0,75 3 (0,75, 0,75, 0,5, 0,5, 0,5) 2,625
5. 0,2 1 4 (1, 1, 0,(6), 0,(6), 0,(6)) 3,5
6. 0,25 1,25 5 (1,25, 1,25, 0,8(3), 0,8(3), 0,8(3)) 4,375
7. 0,3 1,5 6 (1,5, 1,5, 1, 1, 1) 5
8. 0,35 1,75 7 (1,75, 1,75, 1,1(6), 1,1(6), 1,1(6)) 5,375
9. 0,4 2 8 (2, 2, 1,(3), 1,(3), 1,(3)) 5,75
10. 0,45 2,25 9 (2,25, 2,25, 1,5, 1,5, 1,5) 6,125
11. 0,5 2,5 10 (2,5, 2,5, 1,(6), 1,(6), 1,(6)) 6,5
12. 0,55 2,75 11 (2,75, 2,75, 1,8(3), 1,8(3), 1,8(3)) 6,875
13. 0,6 3 12 (3, 3, 2, 2, 2) 7,2
14. 0,65 3,25 13 (3,25, 3,25, 2,1(6), 2,1(6), 2,1(6)) 7,28
15. 0,7 3,5 14 (3,5, 3,5, 2,(3), 2,(3), 2,(3)) 7,3
16. 0,75 3,75 15 (3,75, 3,75, 2,5, 2,5, 2,5) 7,31
17. 0,8 4 16 (4, 4, 2,(6), 2,(6), 2,(6)) 7,32
18. 0,85 4,25 17 (4,25, 4,25, 2,8(3), 2,8(3), 2,8(3)) 7,335
19. 0,9 4,5 18 (4,5, 4,5, 3, 3, 3) 7,35
20. 0,95 4,75 19 (4,75, 4,75, 3,1(6), 3,1(6), 3,1(6)) 7,365
21. 1 5 20 (5, 5, 3,(3), 3,(3), 3,(3)) 7,2
22. 1,05 5,25 21 (5,25, 5,25, 3,5, 3,5, 3,5) 6,945
23. 1,1 5,5 22 (5,5, 5,5, 3,(6), 3,(6), 3,(6)) 6,67
24. 1,15 5,75 23 (5,75, 5,75, 3,8(3), 3,8(3), 3,8(3)) 6,4
№ г уОиг Т 02 Ж
25. 1,2 6 24 (6, 6, 4, 4, 4) 6,13
26. 1,25 6,25 25 (6,25, 6,25, 4,1(6), 4,1(6), 4,1(6)) 6,155
27. 1,3 6,5 26 (6,5, 6,5, 4,(3), 4,(3), 4,(3)) 6,16
28. 1,35 6,75 27 (6,75, 6,75, 4,5, 4,5, 4,5) 6,04
29. 1,4 7 28 (7, 7, 4,(6), 4,(6), 4,(6)) 5,78
30. 1,45 7,25 29 (7,25, 7,25, 4,8(3), 4,8(3), 4,8(3)) 5,515
31. 1,5 7,5 30 (7,5, 7,5, 5, 5, 5) 5,18
Таблица конечна. Далее пропускная способность увеличиваться не может, так как матрица
Я =
Г 1,5 1,5 6 1 1
1,5 1,5 1 6 1
1,5 1,5 1 1 1
1,5 1,5 1 1 1
1,51 1,5 1 1 1
задает сеть стремя аттракторами. Это вершины VI, V2 и V5 - для всех трех вершин выполняется критерий аттрактивности (формула (1)). В самом деле, для вершин VI, V2 и V5 выполняется:
оМ Г2
1* <?2
7,5
= 30;
г
ОМ
5_
д д2 0,25 ' с^* 1/6 при этом из таблицы 1 следует, что для этой сети Т = 30.
Представим данные, полученные в таблице в виде графика (рис. 3). Из него видно, что изменения 3!¥ с ростом пропускных способностей аттракторов происходят немонотонно. Сначала эта величина возрастает, причем присутствуют участки линейного роста; при пропускных способностях, равных приблизительно 0,95, достигает максимума, и затем начинает убывать. Однако при пропускных способностях, равных 1,25 и 1,3, монотонность вновь нарушается.
5
= 30;
ОШ
Рис. 3. Зависимость величины ЗЖот выходных пропускных способностей аттракторов (и порогового значения Т)
4.3. НЕЛИНЕЙНЫЕ ПЕРЕХОДЫ ПРИ ФУНКЦИОНИРОВАНИИ СЕТИ И НЕОДНОРОДНАЯ ЦЕПЬ МАРКОВА
Функционирование несимметричной сети с ресурсом, большим порогового значения, описывается неоднородной цепью Маркова.
Для того чтобы вся сеть при большом ресурсе функционировала по правилу 2 и при этом ее состояния на каждом такте совпадали с состояниями сети, работающей по двум правилам, необходимо и достаточно на каждом такте изменять пропускные способности петель вершин, имеющих избыток ресурса, так чтобы весь избыток уходил в петлю и возвращался на следующем такте. Таким образом, каждая вершина будет отдавать весь свой ресурс.
На каждом такте ^ строится матрица пропускных способностей Я(0, отличающаяся от матрицы Я только диагональными элементами вершин, которые функционируют по правилу 1. То есть
г «)=
Г, если цг ^) < г°,
Ч ) "Е Г ' если ^) > Г °
1 *г
<
Нормируя матрицы R(t), получим стохастические матрицы R'(t), которые задают неоднородную цепь Маркова. Вектор
m
Q(t +1) = б(0)П R'(t )
t=0
совпадает с вектором состояния сети, заданной матрицей R, с начальным состоянием 0(0).
Если в начальном состоянии весь ресурс сосредоточен в источнике с номером j, матрица R(t) будет отличаться от исходной матрицы R только j-й строкой. В процессе функционирования некоторые вершины могут оказываться в зоне Z+(t), и тогда соответствующие им строки матрицы R(t) будут отличаться от строк матрицы R. За конечное число шагов зона Z+(t) стабилизируется, и в ней остаются только потенциальные аттракторы. Все остальные строки матрицы R(t) и соответствующей ей матрицы R'(t) совпадают с соответствующими строками матриц R и R'.
Опишем процессы, происходящие в сети, когда одна из вершин изменяет правило функционирования. Снова для простоты рассмотрим сеть с двумя потенциальными аттракторами. Пусть W > T + 5W, тогда при начальном состоянии, в котором весь ресурс сосредоточен в одном источнике, существует, как минимум, три точки смены правила функционирования: первый и второй аттракторы по очереди переходят на правило 1, источник с ресурсом переходит на правило 2.
Рассмотрим этапы, на которые делится функционирование сети.
0. В начальном состоянии в зоне Z+(0) находится единственная вершина - источник с ресурсом W > T.
1. Смена правила в первой вершине-аттракторе. Поскольку в начальном состоянии этот аттрактор не имел ресурса, существует как минимум один такт времени, в который эта вершина функционирует по правилу 2. Через конечное число тактов эта вершина переходит на правило 1. Зона Z+(t) содержит две вершины: аттрактор и источник.
2. Второй аттрактор переходит на правило 1. Зона 2+(0 содержит три вершины: два аттрактора и источник. Выходной поток из аттракторов стабилизируется.
3. Источник покидает зону 2+(Х). В зоне 2+(Х) остаются два аттрактора.
На промежутке между переходами 1 и 2 происходит формирование поправки 5Ж. На промежутке между переходами 2 и 3 начинается формирование разности ресурса, соответствующее Я'™ЬтгЬ . Если переход 2 отсутствует, второй аттрактор не получает ресурс, больший своей выходной пропускной способности.
После перехода 3 зона 2+(Х) стабилизируется.
Строки стохастических матриц Я'(0, соответствующие вершинам, функционирующим по правилу 1, изменяются на каждом такте. Рассмотрим матрицы, задающие изменение вектора исходящего потока. Они изменяются только при переходах 1-3, в остальные моменты времени они стационарны. Представим матрицу стохастическую матрицу Я в виде: Я = (г/,г2' ,...,Г' ), где гк - столбцы Я\ Пусть V] - вершина-источник содержащая ресурс в начальном состоянии, VI, V2 - аттракторы, причем, пусть, для определенности первый переходит в зону 2+(Х) аттрактор VI. Тогда изменение исходящего потока на последовательных тактах в соответствии с этапами 0-3 можно записать следующим образом:
0. ^(г +1) = ^(0(гДг2\г,',..., в1 ,•••, г'), х < ¿1;
1. ^(г +1) = ¥°и(0^,г2\г,',..., е},..., гп'), ¿1 < X < Х2;
2. ^(г +1) = ^(о(е„е2,г,',..., е1,..., г„'), ¿2 < X < Хз;
3. ^(г +1) = Р°ш(г)(е1,е2,г,',..., г/,..., гп' ), X > х2.
Обозначим матрицы, на которые умножается вектор через Р], Ру, Рщ и Р12 соответственно. Для каждого из четырех этапов входящий поток вычисляется по формуле
г'п (г +1) = Р"* (г)я '.
Из сказанного выше следует, что ЗЖ вычисляется как первая компонента вектора
т = [р" (г+1) - тои' (г+1))= (г) я ' - ^ (ор; )=
í=í1 í=í1
= £ Роггг (Г)(я'-р;).
г=г1
Матрица (К - Ру) имеет два ненулевых столбца с номерами 1 и
Преобразуем это выражение далее, воспользовавшись зависимостью векторов исходящего потока на двух последовательных тактах. В рамках второго этапа ¥оиг (г +1) = ¥оиг(г)р ., следовательно, (г) = ¥°г (г1)(р1; )г-г1.
Тогда
Вектор 5\¥ имеет ненулевые компоненты с номерами 1 и соответствующие аттрактору и источнику; = (<5V) .
Значения 5Жг для сети с I аттракторами (/ = 1, ..., I) определяются аналогично. Количество промежуточных переходов будет I + 1 - по одному на каждый аттрактор и один для источника.
4.4. ПРЕДЕЛЬНОЕ СОСТОЯНИЕ СЕТИ ПРИ¡Шт(?) > Т
Заметим, что значения (/ = 1, ..., I) максимальны, когда ресурс в начальном состоянии сосредоточен только в одной вершине-источнике, т.е. асимметрия его распределения максимальна. При других начальных состояниях 5Ж будут меньше. В предельном случае 55Ж = 0. Так, например, эксперименты показывают, что если во все вершины в начальном состоянии поместить ресурс, больший предельного значения при Ж = Т, то весь излишек будет распределяться в пропорции, заданной матрицей Я^огЪ. Докажем, что это верно для любой несимметричной регулярной сети.
Теорема 2. Пусть регулярная несимметричная ресурсная сеть с l потенциальными аттракторами (l > 1) имеет суммарный ресурс W> T и начальное ее состояние Q(0) представимо в виде Q(0) = Q + AQ(0), где Q - предельное состояние при
W = T, AQ(0) - произвольный вектор, задающий распределение излишков ресурса W- T. Тогда предельное состояние представимо в виде суммы двух векторов:
(13) Q* = Q + AQ(0)RZsorb.
Доказательство. Когда в регулярной сети все вершины имеют ресурс, не меньший соответствующей компоненты вектора Q, поток в сети состоит из двух частей: Fout (t) = F + AFout (t) , где стабильная компонента F равна предельному потоку при любом ресурсе W > T, а нестабильная компонента AF°ut(t), уменьшающаяся со временем, отвечает за перераспределение излишков ресурса AQ(t). Причем, поскольку аттракторы при W = T отдают по полной пропускной способности и больше отдать уже не могут, то соответствующие им компоненты вектора AF°ut(t) равны нулю - ресурс, соответствующий AQ(t), в аттракторы входит, но не может выйти. Таким образом, излишки распределяются по закону:
AQ(t + 1) = AQ(t) Rabsorb, где
В предельном состоянии:
( Е 01
IR R2 J
AQ* = AQ(0)
Е Oi
(e2 - R) 1R O2
= AQ(0) Ratorb. □
4.5. НАЧАЛЬНЫЕ СОСТОЯНИЯ, НЕ СОЗДАЮЩИЕ ПОПРАВКУ
Теорема 2 содержит достаточное условие того, что ресурс в аттракторах распределится в той же пропорции, что и в стоках соответствующей поглощающей сети. Однако в ней описан узкий класс начальных состояний, приводящих к такому предельному состоянию. Предельное состояние описывается формулой (13) для еще одного класса начальных состояний.
Теорема 3. Пусть в регулярной несимметричной ресурсной сети с I аттракторами с номерами 1, ..., I и к источниками с номерами 1+1, ..., 1+к начальное состояние задано вектором:
где величины
Щ = Е Щ/ (] = /+1, ..., /+к)
г^{\,...},1+к+\,...п1} Гг
- необходимые минимальные излишки в каждом источнике; гП - суммарная входная пропускная способность вершин из Z~(0) по ребрам из вершин-источников, ^¿(0) - дефициты1 всех вершин из 2^(0): 7 = 1, ..., 1,1+к+1, ..., п, АQj > 0, ] = 1+1, ., 1+к. Тогда предельное состояние представимо в виде
е=е+Аб • я ' ^,
где АQ - вектор длины п с ненулевыми компонентами АQj (] = 1+1, ..., 1+к).
Доказательство. В теореме утверждается, что если все источники имеют ресурс и его количество достаточно велико, то накопления поправки ЗЖ в аттракторах не происходит. Чтобы доказать это утверждение и оценить количество ресурса, которое должно находиться в каждом из источников, обратимся к описанию смены правил вершинами в сети, данному в разделе 4.3. Накопление величины поправки происходит, когда при функционировании сети случается асинхронный переход аттракторов из зоны 2 в зону 2+(Х): пока в аттракторах из 2+(Х) происходит накопление, аттракторы, находящиеся в 2"(Х), продолжают отдавать весь свой ресурс.
При достаточно больших ресурсах во всех источниках аттракторы переходят на правило 1 одновременно. Следовательно, накопления поправки ресурса не происходит и ЗЖ = 0.
1 Под дефицитом вершины V, мы понимаем количество ресурса, недостающее до значения .
Определим понятие «достаточно большого ресурса».
Для каждой вершины V,, которая в состоянии 0(0) имеет дефицит ресурса, обозначим через г'"- ее суммарную входную пропускную способность из вершин-источников. Дефицит вершины V, обозначим, как и прежде, &(0). Тогда, чтобы покрыть дефицит каждой вершины до величины, соответствующей компоненте 0 , источник V] должен отдать
I
йг (0)
г
ге{1,...Лг+к+1,...,п} Г'
ресурса, где 1 пробегает множество индексов вершин, имеющих дефицит в начальном состоянии.
Отсюда АЖ7 = I ^^ ги .
" } ¿—I 1П- }'
1е{1,...1,1+к+1,...,п} Г'
Соответствующая компонента вектора А0 вычисляется по формуле: А0,. = (0) - - АЖ/. □
Проиллюстрирует теорему примером, в котором предельное состояние может быть получено из начального аналитически.
Пример 6. Сеть задана матрицей (5). Пусть ресурс Ж> Т присутствует в обоих источниках, причем, в каждом источнике его достаточно много: 0(0) = (0, 0, 20, 40, 0). Тогда предельное состояние для этой сети будет = (20, 30, 3,(3), 3,(3), 3,(3)).
Опишем алгоритм нахождения 0*. Все его компоненты, кроме двух первых, совпадают с компонентами вектора 0 . Для расчета первых двух компонент найдем дефициты вершин в начальном состоянии и пропорцию восполнения этих дефицитов из обоих источников, содержащих излишки ресурса. Каждый источник должен отдать ресурс в нейтральную вершину в количестве
э! ,
2 э
и в вершины-аттракторы по 5 единиц - чтобы они получили гО' и гО' единиц ресурса. Итого каждый из двух источников дол-
жен отдать 20/3 ресурса для насыщения сети и оставить себе 3,(3) = 10/3.
лш 20 ~ 10
А^3,4 = у , ~3,4 = у .
Излишек вершины У3 составит
АQз = 20 - 20 - 10 = 20 - 10 = 10, Q 3 3 ,
в вершине У4 остается излишек АQ4 = 40 - 10 = 30. Эти излишки вершины должны распределить между аттракторами в пропорциях 3:1 и 1:3 соответственно (рис. 4).
Рис. 4. Распределение излишков ресурса между аттракторами
Имеем: д! = 5 + 7,5 + 7,5 = 20; д * = 5 + 2,5 + 22,5 = 30. Эти значения соответствуют первым двум компонентам вектора предельного состояния.
4.6. ОЦЕНКА ПОПРАВКИ ЗЖ
Теоремы 2-3 описывают начальные состояния, при которых для каждого из аттракторов ЗЖ, = 0. Своего максимума величины ЗЖ, достигают, если ресурс в начальном состоянии сосредоточен в одном источнике.
Теорема 4. Пусть в регулярной несимметричной ресурсной сети с I аттракторами (I > 2), имеющими номера 1, ..., I, суммарный ресурс W> Т + АЖ/(0) (АЖ/(0) > 5Щ в начальном состоянии находится в вершине-источнике V. Тогда предельное состояние представимо в виде
6=6+аос^ чт, т2, т, о,0),
где А<2 = 6(0) — ^ 0,...,0, ^ , ненулевая компонента
имеет номер
Причем для аттрактора с номером 10 = а^шт г'аЬтгЬ ]к
выполняется 5\¥, = 0, а остальные значения ЗЖ, находятся по
10
формуле
'ад
г
(14) Щ = (д* — гГ) — (д* — С') М]г
СО
r
absorb jl0
Доказательство. Теорема 4. двойственна теореме 1. Вер-
7 • (СО <>
шина-аттрактор с номером /0 = arg min r absorb m , последней из
j 0 j
аттракторов перейдет на правило 1. Излишек SW\ в каждом аттракторе накапливается, пока хотя бы один аттрактор функционирует по правилу 2. Таким образом, вершина, перешедшая на правило 1 последней, не имеет излишка.
Формула (14) - следствие из теоремы 1. Когда все аттракторы функционируют по правилу 1, ресурс перераспределяется между ними в пропорции, задаваемой матрицей R}™bsorb . Излишки в каждой вершине-аттракторе выражаются через компоненту вектора предельного состояния, соответствующую аттрактору, для которого SW^ = 0 . □
Теорема объясняет результаты, полученные в приме pax 2-4. Пусть даны два вектора начальных состояний:
öi(0) = (ft(0),q2(0), ...,qM+AWM,...,ql+k+AWl+k,...,qnm, ß2(0) = (<7i(0),<fe(0), ...,ql+1 +2AWl+l,...,ql+k +2AWl+k,...,qn(Q)\
АЖг - любые заведомо большие величины (например, можно считать, что АЖ > Т). Предельные состояния для них:
01* = +Аая' :ЫогЬ + (т, т,..., т, 0, ..., 0),
02* = 2 + А&Я'Ъо-ъ + (Ж 5Ж2, ¿Ж, 0, 0).
Разность между ними равна
А А ~
02 " 01 = й + (А02 ' :ЫогЪ .
Таким образом, разность двух предельных состояний не зависит от поправок в аттракторах.
Компоненты вектора предельного состояния, соответствующие потенциальным аттракторам, можно вычислить через потоки.
Теорема 5. Пусть регулярная несимметричная ресурсная сеть имеет потенциальные аттракторы с номерами 1, ..., I. Пусть суммарный ресурс Ж> Т + АЖ/(0) в начальном состоянии находится в не-аттрактивных вершинах. Тогда для каждого из потенциальных аттракторов выполнится соотношение
(15) !(/Т(0 - гГ)=я* - С.
'Я (')>г,Ш
Доказательство. Левая часть равенства (15) характеризует ресурс, накопленный сверх суммарной выходной пропускной способности вершины на протяжении всего функционирования сети - каждое слагаемое равно разности входного и выходного
* о"' |—|
потоков, а это и есть предельная разность я - г . и
Следствие. Если в начальном состоянии Q(0) ресурс сосредоточен в одной неаттрактивной вершине, и его величина Ж > Т такова, что в предельном состоянии по крайней мере один из потенциальный аттрактор имеет ресурс, равный своей суммарной выходной пропускной способности: я* = г™", т < I, а при любом его увеличении Ят > г'0"', то для всех остальных аттракторов, имевших я* > г'0"', выполнится:
!(т (') - г0"' )=Щ .
'Я (')>г°"'
Доказательство вытекает непосредственно из формулы (15). □
Отсюда следует, что ¿¿Ж, = я* - г°и.
Пример 7. Пусть сеть задана матрицей (5) и суммарный ресурс Ж = Т + 5Ж\ = 27,2 в начальном состоянии находится в вершине Уз: Q(0) = (0, 0, 27,2, 0, 0). Предельное состояние для этой сети будет Q* = (12,2, 5, 3,(3), 3,(3), 3,(3)).
5Ж\ = 7,2 может быть найдена как сумма разностей двух графиков в дискретных точках (рис. 5):
Рис. 5. Разность входящего и исходящего потоков
Из рисунка видно, что вклад в величину ¿¿Жг дают только первые несколько тактов функционирования сети.
5. Заключение
В работе исследовано распределение излишков ресурса между несколькими аттракторами в регулярных несимметричных сетях.
Получены результаты, позволяющие находить пропорцию распределения ресурса Ж - Т при определенных начальных состояниях.
Показано, что в общем случае эта пропорция повторяет пропорцию индуцированной поглощающей сети за исключением некоторых аддитивных поправок, названых «поправками на регулярность». Эти поправки оценены сверху. Таким образом, чем больше в сети суммарный ресурс, тем меньше относительная погрешность, с которой может быть вычислено предельное состояние.
Рассмотрены переходные процессы при распределении ресурса, когда происходит смена правил функционирования вершин. Описана динамика изменения ресурса в зависимости от весов исходящих из аттракторов ребер.
Однако задача аналитического нахождения предельного состояния в аттрактивных вершинах при произвольном начальном состоянии не решена. Статья содержит много примеров, чтобы проблема была описана максимально полно, и автор надеется, что то, что не получилось у него, может получиться у проницательного читателя.
Литература
1. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Сходимость и устойчивость в задачах согласования характеристик (обзор базовых результатов) // Управление большими системами. -2010. - №30.1 «Сетевые модели в управлении». -С. 470-505.
2. ГАНТМАХЕР Ф.Р. Теория матриц. - М.: Физматлит. 2004. - 560 с.
3. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. III. Исследование предельных состояний // Автоматика и телемеханика. - 2012. - №7. - С. 67-77.
4. ЖИЛЯКОВА Л.Ю. Управление предельными состояниями в поглощающих ресурсных сетях // Проблемы управления. -2013. - №3. - С. 51-59.
5. КЕМЕНИ ДЖ., СНЕЛЛ ДЖ. Конечные цепи Маркова. - М.: Наука, 1970. - 272 с.
6. КУЗНЕЦОВ О.П. Однородные ресурсные сети. I. Полные графы // Автоматика и телемеханика. - 2009. - №11. -С. 136-147.
7. ФОРД Л.Р., ФАЛКЕРСОН Д Р. Потоки в сетях. Пер. с англ. - М.: Мир, 1996. - 334 с.
8. AHUJA R.K., MAGNANTI T.L., ORLIN J.B. Network Flows: Theory, Algorithms, and Applications. - Prentice Hall, United States, 1993. - 864 p.
9. BAK P. How Nature Works: The Science of Self-Organized Criticality. - New York: Copernicus, 1996. - 212 p.
10. BAK P., TANG C., and WIESENFELD K. Self-organized criticality. // Physical Review. - 1988. - A 38. - P. 364-374.
11. BJORNER A., LOVASZ L. Chip-firing game on directed graphs // J. Algebraic Combinatorics. - 1992. - No. 1. -P.305-328.
12. 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.
13. DHAR D. Self-organized critical state of sandpile automation models // Physical Review Letters - 1990. - No. 64. -P.1613-1616.
14. IVASHKEVICH E.V., PRIEZZHEV V.B. Introduction to the sandpile model // Physica. - 1998. - A 254. - P. 97-116.
15. HAJNAL J. The ergodic properties of non-homogeneous finite Markov chains // Proc. Cambridge Philos. Soc. - 1956. -Vol. 52. - P. 67-77.
16. HAJNAL J. Weak ergodicity in non-homogeneous Markov chains // Proc. Cambridge Philos. Soc. - 1958. - Vol. 54. -P. 233-246.
17. 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.
18. PRISNER E. Parallel Chip Firing on Digraphs // Complex Systems. - 1994. - No. 8. - P. 367-383.
19. SENETA E. Non-negative Matrices and Markov Chains. -Springer, 2006. - 279 p.
20. SPEER E.R. Asymmetric abelian sandpile models // Journal of Statistical Physics. - April 1993. - Vol. 71, Issue 1-2. -P. 61-74.
ALLOCATION OF RESOURCE AMONG ATTRACTOR-VERTICES IN NONSYMMETRIC REGULAR RESOURCE NETWORKS
Ludmila Zhilyakova, Institute of Control Sciences of RAS, (Moscow, Profsoyuznaya st., 65), Dr.Sci., Leader Scientist, zhilyako-va.ludmila@gmail.com.
We study nonsymmetric regular resource networks with several attractor-vertices. In these networks there exists a threshold value of total resource W = T, such that for any initial distribution the limit state is uniquely determined, but when W > T extra resource AW = W- T allocates among attractor-vertices depending on the initial state. We prove that this allocation obeys the same law as the allocation in corresponding absorbing network (derived from asymmetric by deletion output edges of attractors). However, there are adjustments dependent on the graph characteristics and the initial distribution of resources. The upper bounds of these adjustments are estimated. The initial states are determined that lead to precise limit states with no adjustments.
Keywords: resource network, graph dynamic threshold model, attractor-vertices, limit state.
Статья представлена к публикации членом редакционной коллегии П.Ю. Чеботарев.
Поступила в редакцию 31.08.2015.
Опубликована 31.03.2016.
