CC BY
55
УДК 519.177 ББК 22.18
РЕСУРСНАЯ СЕТЬ С ОГРАНИЧЕНИЯМИ НА ЕМКОСТЬ АТТРАКТОРОВ1
Жилякова Л. Ю.2
(ФГБУН Институт проблем управления им. В.А. Трапезникова РАН, Москва)
Рассматриваются регулярные несимметричные ресурсные сети, в которых вершины-аттракторы, аккумулирующие ресурс сверх порогового значения, имеют ограничения на емкость. В этом случае излишки ресурса перераспределяются в некоторое множество вершин, названных вторичными аттракторами. Изучается функционирование сети с ограничениями на емкость аттракторов при разных значениях суммарного ресурса и определяются пороговые значения, при которых сеть изменяет поведение.
Ключевые слова: ресурсная сеть, графовая динамическая пороговая модель, аттракторы, ограничения на емкость.
1. Введение
Ресурсная сеть, впервые предложенная в [8], - графовая динамическая потоковая модель, в которой в дискретном времени происходит перераспределение ресурса между вершинами. Сеть представлена ориентированным взвешенным графом; веса ребер обозначают их пропускные способности.
Вершины распределяют ресурс в исходящие ребра по двум разным правилам с пороговым переключением. Переключение
1 Работа выполнена при частичной финансовой поддержке РФФИ (проекты № 14-01-00422а, 15-07-02488а).
2 Людмила Юрьевна Жилякова, доктор физико-математических наук ^НПуакоуа. ludmila@gmail. сот)
происходит в зависимости от количества ресурса, содержащегося в вершине. Если это количество превосходит сумму весов всех исходящих из нее ребер, в каждую смежную вершину передастся ресурс, равный пропускной способности соответствующего ребра. В противном случае вершина отдаст весь свой ресурс, разделив его пропорционально пропускным способностям. Модель параллельна: на каждом такте дискретного времени ^ вершины, имеющие ресурс, отдают его одновременно.
Для такой модели получены результаты, позволяющие для каждой топологии сети и пропускных способностей ее ребер, для любого суммарного ресурса и его начального распределения по вершинам определить предельное состояние и предельный поток. Определена формула, которая позволяет для произвольной сети вычислить агрегированный показатель -пороговое значение ресурса Т [6]. При суммарном ресурсе в сети, не превосходящем Т, ее функционирование, начиная с некоторого конечного такта, описывается однородной цепью Маркова [4]. Более интересные результаты получены для ресурсов выше порогового значения. В этом случае процесс перераспределения ресурса описывается неоднородной цепью Маркова и в сетях с разными топологиями возникают различные эффекты.
При малых ресурсах (меньших порогового значения Т) функционирование сети полностью определяется стохастической матрицей, полученной из матрицы пропускных способностей нормированием строк [4]. Предельное состояние для регулярных сетей (которые исследуются в настоящей работе) всегда существует и единственно; вектор предельного состояния пропорционален вектору предельных вероятностей однородной цепи Маркова. На сети с малым ресурсом переносятся все результаты, полученные для случайных блужданий и процессов рассеяния на графах [3, 7, 9, 14, 17]. Модели, основанные на однородных цепях Маркова, используются во многих предметных областях, в частности, для решения задач о нахождении консенсуса в многоагентных системах и социальных сетях [1, 2], при исследовании распространения эпидемий, в моделировании регуляторных сетей экспрессии генов и многих других.
Разнообразные прикладные приложения, основанные на модели рассеяния на графах, описаны в книге [14].
При больших ресурсах ресурсная сеть наследует часть свойств целочисленных пороговых моделей, таких как, в частности, игра выстреливания фишек, в которой вершины графа обмениваются фишками с выполнением закона сохранения. Вершина может выстрелить, если ее суммарное количество фишек больше количества выходных ребер. Тогда она посылает по одной фишке вдоль каждого ребра. Классическая игра представляет собой последовательные выстреливания вершин [12, 13]; очередность вершин выбирается случайным образом. Были исследованы и параллельные игры выстреливаний [18]. Эта математическая модель нашла применение при описании процессов, получивших название «самоорганизованная критичность», таких, в частности, как сход лавин [10, 11, 15, 16].
Ресурсная сеть при больших ресурсах, объединяя в себе свойства моделей рассеяния и пороговых моделей, демонстрирует проявление новых свойств. К ним относится способность вершин отдавать ресурс независимо от его количества (таким образом, обеспечивается устойчивость и сходимость, которых нет в целочисленных моделях); способность некоторых вершин накапливать излишки ресурса (чего нельзя достичь в однородных марковских процессах), при этом все остальные вершины имеют ресурс, пропорциональный компонентам вектора предельных вероятностей соответствующей однородной цепи Маркова [4-6].
2. Основные определения, обозначения и предварительные результаты
В настоящем разделе будут конспективно перечислены основные определения и предварительные результаты, полученные в [8, 4-6] и др.
Ресурсная сеть представляет собой ориентированный граф О = (V, Е), IV = п. Ребра ву имеют неотрицательные веса ггу, которые обозначают их пропускные способности; Я = (гу)пхп -матрица пропускных способностей.
Вершины сети содержат ресурс, который перераспределяется в дискретном времени: qi(t) > 0 - количество ресурса в вершине v, в момент времени t. Состояние сети в момент t описывается вектором Q(t) = (qi(t), ..., qn(t)). Предельное состояние сети будем обозначать через Q*.
Обозначим через r'" и r""1 сумму всех пропускных способностей входящих и исходящих ребер соответственно. Будем называть эти величины входной и выходной пропускными способностями вершины v,.
Правила функционирования сети. В момент t вершина v, по ребру e,m отдает в смежную ей вершину vm:
rim единиц ресурса, если qi(t) > r°ut (правило 1);
r
-Jim-q (t) единиц ресурса, если q (t) < rout (правило 2).
rout
Суммарный ресурс сети обозначается через W. При функционировании выполняется закон сохранения: W = const.
Множество вершин, для которых q,(t) < r°ut, называется зоной Z (t). Вершины из Z (t) функционируют по правилу 2. Z+(t) - множество вершин, ресурс которых больше их выходной пропускной способности, они функционируют по правилу 1.
Т - пороговое значение ресурса, такое что при W < Т все вершины, начиная с некоторого t', переходят в зону Z-(t); при W > T зона Z+(t) не пуста и не изменяется, начиная с некоторого t". Для несимметричных сетей существование и единственность Т доказаны в теореме 3 [4], в дальнейших работах было показано, что Т существует и единственно в любой эргодической сети (сети без стоковых вершин).
Для предельного состояния Q обозначим эти зоны через Z-* и Z+*.
Доказано, что при W < Т сеть любой топологии эквивалентна однородной цепи Маркова, заданной стохастической матрицей R', полученной из матрицы пропускных способностей R нормированием по строкам (этот результат для регулярных сетей был получен в [4] и впоследствии перенесен на остальные классы сетей). Процессы перераспределения ресурса в регуляр-
ных сетях при малом ресурсе всегда сходятся. Предельное состояние существует и единственно (не зависит от начального распределения ресурса по вершинам).
Вектор предельного состояния сети при Ж = 1 совпадает с вектором предельных вероятностей соответствующей цепи Маркова. Обозначим этот вектор через 01*.
При Ж > Т в некоторых вершинах начинают накапливаться излишки ресурса.
Введем обозначение: (1) Д/- = г<" - г°и.
По знаку Дг, вершины делятся на три класса:
1) вершины-приемники, для которых Дг, > 0;
2) вершины-источники, для которых Дг, < 0;
3) нейтральные вершины, для которых Дг, = 0.
Вершины, способные при Ж > Т оказаться в зоне 2+*, называются аттракторами. В роли аттракторов могут выступать некоторые приемники и некоторые нейтральные вершины (критерий аттрактивности будет приведен ниже после введения всех требуемых обозначений), причем первые способны притягивать ресурс и поэтому названы активными аттракторами, вторые могут лишь сохранить ресурс, если он был достаточно большим в начальном состоянии. Такие аттракторы называются пассивными.
Поток в ресурсной сети
Ресурс, выходящий из вершины V, по ребру в- в момент приходит в вершину V]- в момент ^ + 1; между моментами ^ и ^ + 1 он находится в ребре в--. Этот ресурс назовем потоком /-((). Общий поток сети описывается матрицей = (/](0)пхп.
Величиной потока будем называть сумму
п п
fsum () .
1=1
Введем следующие обозначения.
п п
I/ (0 = /Г (0; I /у (0 = 17 ((+1).
1=1 -=1
Будем полагать, что /'п (0) = 0 .
F°ut(t) = (f°ut(t), ... , fOut(t) )T - вектор исходящего потока; Fin(t) = (f'" (t) , ... , f'n (t) ) - вектор входящего потока; F - матрица предельного потока;
n n
fmm - SS fij - суммарный предельный поток; Fn* и F°ut* -'-1 j-1
векторы двух предельных потоков. Входное и выходное соседство.
Назовем входными соседями вершины vi множество V'"1 — [vk\rki > 0}; выходными соседями вершины vi будем называть множество V'0ut — [vk\rik > 0} .
3. Описание сети с ограничениями на емкость аттракторов
3.1. МОТИВАЦИЯ ИССЛЕДОВАНИЯ
В работах [5, 6] исследовано поведение регулярных несимметричных сетей с ресурсом, превышающим пороговое значение Т. Напомним, что регулярной называется сеть, представленная сильно связным графом, у которого длины всех циклов взаимно просты, а несимметричной - сеть, в которой присутствует хотя бы пара вершин, для которых А г, (формула (1)) отлична от нуля. То есть в несимметричной сети необходимо существует хотя бы один источник и хотя бы один приемник. При этом хотя бы один приемник будет в такой сети активным аттрактором.
Несимметричные сети с несколькими активными аттракторами - как правило, искусственные конструкции. В них для нескольких вершин должен выполняться критерий аттрактивно-сти ([6]):
r"ut
j = arg min -¡-¡т,
ie{1,..."} q1
1*
где q - элементы вектора предельных вероятностей однородной цепи Маркова Q1*.
В противном случае, даже если в сети много приемников, активным аттрактором из них может быть только один.
При любом начальном распределении суммарного ресурса Ж > Т в регулярной несимметричной сети с одним аттрактором будет единственное предельное состояние, в котором аттрактор забирает все излишки, сколь большими они бы ни были. (Если в сети два и более аттрактора, излишки ресурса распределятся между ними.) Распределение ресурса между аттракторами зависит от начального состояния. В остальных вершинах предельные значения ресурса постоянны при сколь угодно большом количестве суммарного ресурса и его начальном распределении по вершинам - не-аттрактивные вершины в предельном состоянии всегда содержат ресурс, строго меньший своей выходной пропускной способности, даже если это вершины-приемники [5].
Новая модель, представленная в настоящей статье, предполагает, что аттракторы должны «делиться» излишками ресурса с другими вершинами сети. Это достигается тем, что в аттракторах вводится ограничение на емкость - максимальное количество ресурса, которое может получить аттрактор V, сверх величины г°и' (величине г°и равна соответствующая компонента вектора предельного состояния при Ж = Т).
Если ресурс в сети достаточно велик, то после насыщения аттракторов он начнет распределяться в другие вершины. Новые вершины, принимающие излишки ресурса, назовем вторичными аттракторами. Ясно, что процесс ограничения можно продолжать и дальше, вводя в сети аттракторы все более высоких рангов. «Изначальные» аттракторы (аттракторы сети без ограничений) будем называть первичными. Поскольку процесс ограничения первичных, вторичных и т.д. аттракторов одинаков, рассмотрим лишь переход от первичных аттракторов к вторичным.
3.2. ПОСТАНОВКА ЗАДАЧИ
Как уже было отмечено, вектор предельного состояния при Ж = 1 является вектором предельных вероятностей цепи Маркова, заданной графом сети. Для него введено обозначение 01*.
Вектор предельного состояния при Ж = Т обозначим
через Q .
Введем такую нумерацию вершин сети, при которой первичные аттракторы имеют номера от 1 до I. В [5] было доказано,
что б = (СгГ,) , причем, ~ < ,] = I + 1, ..., п.
Из единственности предельного состояния следует, что
вектор б пропорционален вектору б1* с коэффициентом Т:
б = б1* Т. При этом компоненты с номерами от I + 1 до п не
изменяются при любом Ж > Т, и все излишки распределяются только среди аттракторов (Следствие 5 из теоремы 1 [5]). Вектор предельного состояния для Ж > Т имеет вид
I
б* = (гГ + ля1...,гГ + Л<ЪЙ+1,..,&), где ^Лд* = Ж- Т.
1=1
Предположим теперь, что для всех аттракторов введено ограничение на емкость р, такое что каждый аттрактор может получить лишь ресурс, равный + р и дальнейшее накопление в этой вершине становится невозможным.
Ясно, что при этом появятся новые вершины, накапливающие излишки.
Перечислим основные задачи исследования получившейся модели:
1. Исследование динамики сети и переходных процессов при разном суммарном ресурсе.
2. Нахождение пороговых значений количества суммарного ресурса, при которых происходят качественные изменения функционирования сети - насыщение первичных аттракторов и начало накопления во вторичных. Описание функционирования в каждом из интервалов между двумя соседними пороговыми значениями.
3. Исследование устойчивости: существование предельных потоков и предельных состояний в каждом из интервалов.
4. Нахождение критерия «вторичной аттрактивности».
5. Нахождение векторов предельных состояний и потоков.
В данной работе мы остановимся на первых двух пунктах
(и частично третьем), но сначала опишем алгоритм ограничения емкости аттракторов и обоснуем его использование.
4. Алгоритм функционирования сети
Ограничение на емкость аттракторов формулируется просто и достаточно естественно: если представить вершины в виде некоторых резервуаров ресурса, то у аттракторов на высоте p от значения г°и ставится заслонка и нового ресурса вершина принять не может.
Сделаем одно важное замечание относительно начальных состояний. Если в начальном состоянии в первичном аттракторе содержится ресурс, превосходящий его ограничение, то будем полагать, что такой аттрактор ограничить нельзя. А поскольку заранее неизвестно, какие вершины в сети являются первичными аттракторами, будем рассматривать такие начальные состояния 2(0), в которых ресурс выше значений г™1 + р может содержаться только в вершинах-источниках или в тех нейтральных вершинах, которые гарантированно не будут первичными аттракторами.
4.1. ФУНКЦИОНИРОВАНИЕ СЕТИ ПОСЛЕ НАСЫЩЕНИЯ ПЕРВИЧНЫХ АТТРАКТОРОВ
Представим фазу функционирования сети с ограничениями, на которой происходит накопление ресурса в аттракторах. На этой фазе входной поток в аттракторы превосходит выходной. При наличии ограничений в вершинах-аттракторах на каждом такте должно проверяться условие, не превышены ли их емкости. Если емкость вершины превышена, она должна принять не весь ресурс, а лишь его часть. При этом нужно решить, что делать с оставшимся ресурсом.
При моделировании физического объекта логично было бы оставить излишки ресурса в ребрах. Такое решение хорошо еще и тем, что на каждом последующем такте входящий в первичный аттрактор поток будет строго равен выходящему из нее потоку, потому что остатки ресурса во входящих ребрах пропорционально уменьшили бы их пропускную способность, частично «закупорив» их собой. Тем самым в аттракторе будет поддерживаться стабильный ресурс г™' + р. Однако в модели выполняется закон сохранения. На каждом такте времени весь ресурс должен быть в вершинах, и этот выход не является приемлемым. Поэтому перед каждым тактом вершины анализируют
свой потенциальный входящий поток. Если в некоторой вершине на некотором ненулевом такте ^ входной поток избыточен, т.е. при принятии всего потока выполнится неравенство qi (г) > г°и + р, она оповещает свое входное соседство, и все вершины Ук е V™ по соответствующим ребрам отдают не весь планируемый поток, а пропорционально убавляют его на некоторую часть с коэффициентом с() е (0, 1). Суммарный входной поток уменьшается с этим же коэффициентом:
Лп-п™(') = Лп (') - с-(г) ■ лп(г) = (1 - с-(г))• Лп(г).
При этом, поскольку вершина функционирует по правилу 1, на каждом такте она отдает ресурс, равный своей выходной пропускной способности гсми. Чтобы удержать в ней стационарный ресурс, равный г'Г' + р, дивергенция потока (разность между входным и выходным потоком) в этой вершине должна быть равна 0, и, следовательно, входной поток тоже должен быть равен г°ш .
Имеем:
(1 - с-(г))• /Г (г) = г°и;
откуда коэффициент пропорциональности
(2) с.(г) = 1 = /Ш-Г = ДШ 1 Лп (г) Лп (') Лп (')
где Ддг(0 - потенциальный избыток ресурса в аттракторе V;.
Коэффициент с() зависит только от входного потока в вершину-аттрактор V*. Неотданный ресурс остается в этих вершинах на такте (.
В следующем разделе формула для коэффициента с() претерпит некоторые изменения.
4.2. МАРКИРОВКА ПЕРВИЧНЫХ АТТРАКТОРОВ
Для того чтобы функционирование сети с ограничениями на аттракторы было корректным, первичные аттракторы необходимо пометить до начала процесса перераспределения ресурса. Это условие, на первый взгляд, выглядит искусственным и усложняет модель. Однако на практике оказывается, что в сети, которая о себе ничего не знает (до начала ее работы неизвестно,
какие вершины являются первичными аттракторами), не всегда возможно корректно поставить ограничения на емкость.
Рассмотрим две иллюстрации этого утверждения.
1. Представим себе, что в сети имеется более одного первичного аттрактора, и рассмотрим некоторую пару из множества первичных аттракторов. Пусть количество ресурса настолько велико, что начиная с некоторого такта ^ оба они набрали бы максимально допустимое количество ресурса, и при этом в момент ^ в некоторых источниках еще содержался бы избыток ресурса, который они отдают в аттракторы. Поскольку (по определению первичных аттракторов) при таких процессах их входной поток превосходит выходной, на каждом такте в обоих аттракторах появляются потенциальные излишки, которые должны остаться в тех вершинах, откуда они пришли. Представим теперь, что эти аттракторы связаны двусторонней парой ребер. Тогда первый оставляет свои излишки у второго, а второй - у первого. При этом ни первый, ни второй не могут их оставить у себя, потому что имеют максимально возможное количество ресурса + р. Эти излишки должны быть перераспределены между другими вершинами. Но другие вершины уже один раз получили излишки назад и пересчитали свой ресурс. Процесс пересчета ресурса у входящих соседей придется повторить. Избежать такого итеративного распределения излишков на каждом такте можно лишь пометив аттракторы. Тогда будет достаточно не возвращать ресурс из первичных аттракторов в смежные с ними первичные аттракторы, а распределять излишек между не-аттрактивными вершинами из их входного соседства. Формула (2) для расчета коэффициента
пропорциональности с() в этом случае изменится: (,) = ,
г V / п ?
I /кг(( )
к=1+1
если помеченные первичные аттракторы имеют номера от 1 до I.
Весь потенциальный избыток ресурса в аттракторе будет пропорционально делиться лишь между не-аттрактивными вершинами.
2. Если аттракторы в сети не помечены, они определяются во время функционирования следующим образом. Предполага-
ется, что если вершина в начальном состоянии имела малый ресурс, а за время функционирования ей удалось набрать ресурс выше предельной емкости, то эта вершина - активный аттрактор - по-другому определить аттрактор, не зная вектора предельного состояния, невозможно [6]. В большинстве случаев это предположение верно. Однако в некоторых случаях динамика сети может быть такой, что определенные вершины могут сначала набирать, а затем терять ресурс (т.е. на протяжении значительного времени демонстрировать поведение аттракторов, не являясь таковыми). Рассмотрим пример, иллюстрирующий такую динамику.
Пример 1. Пусть сеть с пятью вершинами задана матрицей пропускных способностей:
Я
Это полный граф, два ребра которых имеют пропускные способности, отличные от остальных.
Напомним, что в матрице Я 1-я строка соответствует пропускным способностям ребер, выходящих из вершины V*, а 7-й столбец - пропускным способностям ребер, входящих в вершину V,.
В этой сети два приемника - вершины V: и V2, один источник - vз и две нейтральных вершины V4 и V5. Аттрактор этой сети - вершина V2, несмотря на то, что у вершины V! входная пропускная способность значительно выше V «паразитирует» на Vl).
Рассмотрим функционирование этой сети без ограничения на аттракторы. Пусть вектор начального состояния имеет вид 2(0) = (0, 0, 0, 0, 100), т.е. ресурс 100 находится в нейтральной вершине.
(1 3 1 1 1^1
1 11 1 1
100 11 1 1
1 11 1 1
V 1 11 1 V
Вектор предельного состояния б* = (6,266; 84,104; 3,210; 3,210; 3,210).
Из рис. 1 видно, что прежде чем уйти в зону Z~(í), вершина У1 более ста сорока тактов держится в зоне вместе с вершиной У2, при этом больше семидесяти тактов количество ресурса в ней возрастает. До некоторого такта « 75 обе эти вершины неотличимы по поведению и годятся на роль аттракторов. Если накладывать ограничения на емкость, не зная заранее, сколько в сети аттракторов и каковы они, то ограничения придется наложить на обе вершины, хотя вершина является вторичным аттрактором этой сети, и именно она должна получить все излишки ресурса.
Единственным первичным аттрактором является вершина У2 - в конечном счете, несмотря на немонотонную динамику сети, она притянула к себе все излишки, и именно ее одну нужно ограничивать. Если поставить на нее ограничение р = 10, получим динамику, представленную на рис. 2.
Рис. 1. Функционирование сети без ограничений на аттракторы при большом ресурсе
Рис. 2. Функционирование сети с ограничением на емкость аттрактора р = 10 (вершина У2 имеет ресурс, равный 15; излишки ресурса - в вершине VI)
Предельное состояние при функционировании сети с ограничением будет О* = (74,557; 15; 3,438; 3,503; 3,503).
Однако если предварительно не иметь информации о том, какая из вершин VI и V2 является первичным аттрактором, нет никакой возможности узнать, какую из них нужно ограничивать, когда ресурс в них достигнет величины, равной + р, -для вершины VI эта величина равна 17, для вершины V2, соответственно, 15.
Единственный выход, позволяющий избежать некорректного задания ограничений, как и в предыдущем случае, состоит в том, чтобы, используя критерий аттрактивности, определять первичные аттракторы сети и ставить на них ограничения до начала функционирования.
Замечание. Как видно из сравнения двух векторов предельного состояния, в сети с ограничениями изменились все компоненты вектора О*. Вектор предельного состояния в сети с Ж = 100 без ограничений: О* = (6,266; 84,104; 3,210; 3,210; 3,210), с ограничениями: О* = (74,557; 15; 3,438; 3,503; 3,503). Ресурсы в не-аттрактивных вершинах увеличились. Объяснение
причин увеличения компонент вектора предельного состояния и, как следствие, суммарного потока в сети с ограничениями будет дано в следующем разделе.
5. Исследование потока при ограничениях на аттракторы
В сети без ограничений на аттракторы при Ж > Т для вычисления компонент вектора предельного потока количество избыточного ресурса (сверх Т) не имеет значения. При Ж = Т + А Ж для любого сколь угодно большого значения АЖ величина суммарного предельного потока в регулярной сети неизменна и равна Т (теорема 1 [5]).
Пусть по-прежнему нумерация вершин сети такова, что первичные аттракторы имеют номера от 1 до I. В [5] доказано, что при Ж = Т предельное состояние описывается вектором
(3) £ = (гГ,.,гГ\ ),
независимо от начального состояния. Этим же вектором описываются входной и выходной предельные потоки при любом Ж = Т + АЖ, АЖ > 0:
Рт" = = (гГ,..., гГ, дм,..., ~); £т = Т.
Для доказательства утверждений о величине потока в сетях с ограничениями сформулируем лемму о свойстве активных аттракторов сети без ограничений. Напомним, что активным аттрактором называется вершина-приемник, способная притягивать излишки ресурса при Ж > Т.
Лемма 1. В несимметричной регулярной сети до каждого активного аттрактора существует путь хотя бы из одной неаттрактивной вершины.
Доказательство. Рассмотрим произвольный активный аттрактор VI. Поскольку Vг - приемник, для него выполняется неравенство г'" > г°и'. Если его входное соседство Vсодержит не-аттрактивную вершину, лемма доказана. Если V'" не содержит не-аттрактивных вершин, то все входные соседи VI сами являются аттракторами (активными или пассивными), т.е. это приемники или нейтральные вершины. А это означает, что для
вершин, принадлежащих V-", должны существовать свои входные соседи, которые будут компенсировать ресурс, отдаваемый ими в активный аттрактор Vi, и служить для них источниками. Таким образом, существует по крайней мере один источник, из которого имеется путь в Vi. □
Сформулируем утверждения об оценке величины суммарного потока в сетях с ограничением при разных количествах суммарного ресурса.
Утверждение 1. Пусть в регулярной ресурсной сети с ограничением на аттракторы r°ut + p, i = 1, ..., l, суммарный ресурс равен W = T + AW, AW> 0, и в начальном состоянии ни один аттрактор не имеет ресурс, больший своего ограничения. Тогда:
1) если AW < lp, предельный поток существует и f*um = T;
2) если AW > lp, выполняется неравенство:
(4) fsum (t)|T .
Доказательство. 1) Пусть AW < lp. По определению функционирования сети с ограничениями, ограничение на аттрактор включается лишь тогда, когда его ресурс достиг величины r°ut + p. До этого момента вершина-аттрактор функционирует так же, как если бы ограничения не было.
1.1) Если аттрактор в сети единственен, то условие AW < lp редуцируется до AW < p, предельное состояние описывается вектором Q = (rlout + AW,q2,...,qn) ; предельный поток существует и fS*um = T .
1.2) Рассмотрим сеть с двумя аттракторами V1 и V2. Докажем, что при t ^ да каждый аттрактор будет функционировать без ограничений. Если ни в какой момент ресурс в V1 и V2 не достигает ограничения, утверждение доказано. Предельное состояние описывается вектором
Q = (r°out + AW, ro0ut +AW2, ~3,..., qn ),
AW1 < p, AW2 < p, AW1 + AW2 = AW; предельный поток существует и fs*um = T .
Пусть в некоторый момент времени t ресурс в одном из аттракторов (не ограничивая общности, положим, что это V1)
достиг ограничения г-^" + р. Тогда этот аттрактор возвращает излишки в соседние не-аттрактивные вершины. Но поскольку сеть регулярна, из этих не-аттрактивных вершин существуют пути в вершину V2, которая по определению первичного аттрактора притягивает ресурс. При этом вершина V2 может, как максимум, набрать ресурс, равный г™1 + р, потому что общий излишек ресурса в сети не превосходит величины 2р. Если А Ж = 2р, то вектор предельного состояния будет иметь вид ё = (г^ии + р, г+ р, ~3,..., ), и аттракторы достигнут верхней границы своей емкости, но не будут принимать, и, следовательно, отдавать излишки.
1.3) Для I аттракторов рассуждения аналогичны. Ресурс А Ж < 1р может какое-то время перераспределяться между I аттракторами. Но поскольку хотя бы один аттрактор в сети будет не насыщенным, он будет притягивать к себе ресурс, пока в сети не установится сбалансированный поток /*ит = Т.
2. Докажем что при АЖ > 1р выполняется неравенство (4).
Процессы стабилизации потоков в регулярной несимметричной сети без ограничений при Ж>Т описаны в теореме 1 [5]. Весь ресурс АЖ в сети без ограничений должен распределиться между вершинами-аттракторами. Суммарный поток при этом равен Т. При АЖ > 1р в сети без ограничений гарантированно существует хотя бы один аттрактор V;, в котором ресурс будет больше ограничения на емкость г°и + р (принцип Дирихле).
Тогда при наличии ограничения в этом аттракторе появится избыток входного ресурса, который он будет возвращать своим входным соседям. По лемме 1 в регулярной несимметричной сети в этот аттрактор существует путь от не-аттрактивной вершины. Если эта вершина находится во входном соседстве аттрактора, это означает, что к ней будет возвращаться часть ресурса из вершины V;. А поскольку она с какого-то момента будет функционировать по правилу 2 и отдавать весь свой ресурс, то выходящий поток из нее будет больше, чем в сети без ограничений за счет возврата от аттрактора, и, соответственно, возрастет суммарный поток в сети. Если во входном соседстве Vг не существует не-аттрактивных вершин, то ресурс будет возвращаться к не-аттрактивной вершине вдоль пути, описанного в лемме 1 (в обратном направлении). И в этой вершине за счет
этого возврата поток будет больше, чем в сети без ограничений. Таким образом, при АЖ> /р суммарный поток в сети строго больше значения Т при ^ ^ <х>.П
Легко доказать, что верно обратное утверждение. Следствие. Пусть в регулярной ресурсной сети с ограничением на аттракторы г°"* + р, 1 = 1, ..., /, суммарный ресурс равен Ж = Т + АЖ, АЖ > 0, и в начальном состоянии ни один аттрактор не имеет ресурс, больший своего ограничения. Тогда:
1) если /1 = Т, то АЖ < /р;
2) если /шп (0ЦШ> Т , то АЖ > /р.
Доказательство следует из доказательства утверждения 1.
Замечание. Величинар берется одной и той же для всех аттракторов, потому что в сети ни количество аттракторов, ни сами аттракторы не известны заранее. Они помечаются перед началом функционирования, но уже после распределения ресурса по вершинам. Поэтому разные ограничения можно задавать лишь случайным образом, однако это обобщение не усложнит модель.
Покажем, что в сети с ограничениями существует еще одна интегральная характеристика - второе пороговое значение ресурса Т1.
Утверждение 2. В регулярной ресурсной сети с ограничением на аттракторы + р, 1 = 1, ..., /, существует дополнительное пороговое значение (помимо Т) Т11 такое, что в зависимости от количества ресурса в сети существует четыре режима функционирования с тремя пороговыми переключениями, задающимися константами Т, /р и Т11:
1) при Ж < Тс некоторого конечного момента ^ вся сеть функционирует по правилу 2 и представляет собой однородную регулярную цепь Маркова;
2) при Т < Ж < Т + /р суммарный поток в сети не изменяется с ростом Ж, предельный поток существует, Е'"* = р°ии*т = ^ ,
/* = Т •
J 5ип '
3) при T + lp < W< J11 + lp суммарный поток возрастает с ростом W, суммарный предельный поток существует и
f* = W - lv;
J sum " lF з
4) при W > Tl + lp суммарный поток не изменяется с ростом W, суммарный предельный поток существует и fSUm = T11.
Доказательство.
1) Очевидно, при малом ресурсе функционирование сети с ограничениями эквивалентно функционированию сети без ограничений, для которых пункт 1) выполняется.
2) Доказательство этого пункта следует из утверждения 1.
3) Рассмотрим поток в сети при W > T + lp как функцию от W: fSUm = fsum(W, t). При W = T + lp во всех аттракторах наступает насыщение. Предельный поток существует. Будем постепенно увеличивать W. При появлении избытка сверх T + lp хотя бы в одном из аттракторов «срабатывает» ограничение. Тогда он принимает из вершин-источников не весь ресурс, а его часть; ресурс в источниках, смежных с этим аттрактором (или связанных с ним путем (лемма 1)), а, следовательно, и суммарный поток в сети начинает возрастать. Поскольку сеть регулярна, ограничения срабатывают во всех аттракторах, сколь бы ни было мало значение W - (T + lp). Так как все вершины, кроме аттрактивных, функционируют по правилу 2, они отдают весь имеющийся у них ресурс на каждом такте t, и он становится потоком. Таким образом, суммарный поток в сети выше, чем при T < W < T + lp, и он растет с ростом величины W.
Ресурс в аттракторах стабилизируется на значениях r°ut + v. Все остальные вершины функционируют по правилу 2, следовательно, весь их ресурс является потоком. Тогда независимо от существования предельного потока в сети, суммарный предельный поток существует и fsum (W) = W - lp . Это весь ресурс вершин, за исключением накопленных излишков в первичных аттракторах. При росте W это равенство выполняется до тех пор, пока в предельном состоянии ресурс некоторой вершины, не бывшей аттрактором, не достигнет пороговой величины rout . Тогда эта вершина при дальнейшем увеличении W будет накапливать ресурс (таких вершин может быть несколько). Эти вершины являются вторичными аттракторами. Значение W, при
котором для некоторых не-аттрактивных вершин выполнится Ч* = г°и*, является пороговым. Обозначим его через Т11.
4) при Ж > Т11 вторичные аттракторы функционируют по правилу 1 и притягивают к себе все излишки ресурса. Суммарный предельный поток при этом: /**ип = Т11 .□
Замечание 1. В сетях без ограничений при любом, сколь угодно большом ресурсе, суммарный предельный поток постоянен и равен пороговому значению Т. Именно при этом значении аттрактор получает в предельном состоянии ресурс, равный своей выходной пропускной способности. Аналогично, в сети с ограничениями на аттракторы таким пороговым значением является величина Т11. При Ж > Т11 весь ресурс сверх Т11 поглотится вторичными аттракторами.
Замечание 2. Существование Т11 доказано в утверждении 2. Вопрос о его единственности пока остается открытым, так же как и вопрос о существовании предельных потоков и состояний сети при Ж > Т + /р.
Замечание 3. При р = 0 режим 2) отсутствует, и при возрастании суммарного ресурса вплоть до значения Т11 суммарный предельный поток все время возрастает без промежуточного плато. При р > 0 суммарный предельный поток стабилизируется при возрастании Ж до насыщения емкостей аттракторов.
6. Заключение
В работе предложена модификация графовой динамической пороговой модели «ресурсная сеть». При большом суммарном ресурсе в исходной модели существовали вершины, накапливающие в себе основную часть ресурса при его перераспределении. Как правило, в произвольной сети такая вершина единственна. Для того чтобы изменить это положение вещей и распределить излишки более равномерно, в модели вводятся ограничения на емкости вершин-аттракторов. Тогда при их насыщении излишки ресурса начинают переходить в другое множество вершин. Эти вершины названы вторичными аттракторами.
Настоящая работа посвящена описанию и всестороннему анализу предложенной модели, а также исследованию ее устойчивости. В сети с ограничениями на емкости аттракторов для любого суммарного ресурса суммарный предельный поток определяется однозначно. Доказано существование второго порогового значения суммарного ресурса Т11, при переходе через которое сеть изменяет поведение. Показано, что в сети с ограничениями для значения суммарного ресурса существует четыре промежутка, на которых сеть демонстрирует различное поведение: (0, Т], (Т, Т + 1р], (Т + 1р, Т11 + 1р], (Т11 + 1р, да).
В ходе дальнейших исследований планируется доказать, что в сети существуют предельные состояния и потоки (не только суммарный поток стационарен); доказать единственность Т11 и найти формулу для его нахождения, найти критерий вторичной аттрактивности, найти векторы предельных потоков и состояний.
Автор выражает глубокую благодарность д.ф.-м.н. Р.П. Агаеву за конструктивные замечания и полезные обсуждения в ходе подготовки статьи.
Литература
1. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Метод проекции в задаче о консенсусе и регуляризованный предел степеней стохастической матрицы // Автоматика и телемеханика. -2011. - №12. - С. 38-59.
2. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Сходимость и устойчивость в задачах согласования характеристик (обзор базовых результатов) // Управление большими системами. -2010. - №30.1 - С. 470-505.
3. ГАНТМАХЕР Ф.Р. Теория матриц. - М.: Физматлит. 2004. - 560 с.
4. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. I. Процессы стабилизации при малых ресурсах // Автоматика и телемеханика. - 2011. - №4. - С. 133-143.
5. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. II. Потоки при больших ресурсах и их стабилизация // Автоматика и телемеханика. - 2012. - №6. - С. 103-118.
6. ЖИЛЯКОВА Л.Ю. Несимметричные ресурсные сети. III. Исследование предельных состояний // Автоматика и телемеханика. - 2012. - №7. - С. 67-77.
7. КЕМЕНИ ДЖ., СНЕЛЛ ДЖ. Конечные цепи Маркова. - М.: Наука, 1970. - 272 с.
8. КУЗНЕЦОВ О.П. Однородные ресурсные сети. I. Полные графы // Автоматика и телемеханика. - 2009. - №11. -С. 136-147.
9. РОБЕРТС Ф.С. Дискретные математические модели с приложениями к социальным, биологическим и экологическим задачам. - М. Наука, 1986. - 496 с.
10. BAK P., TANG C., WIESENFELD K. Self-organized criticali-ty // Phys. Rev. A. - 1988. - Vol. 38, No. 1. - P. 364-374.
11. BAK P., CHEN K. Self-organized criticality // Scientif. Am. -1991. - No. 264. - P. 46-53.
12. BJORNER A., LOVASZ L., SHOR P. Chip-firing games on graphs // Europ. J. Comb. - 1991. - No. 12. - P. 283-291.
13. BJORNER A., LOVASZ L. Chip-firing game on directed graphs // J. Algebraic Combinatorics. - 1992. - No. 1. -P. 305-328.
14. 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.
15. DHAR D. The abelian sandpile and related models // Physica A: Statist.l Mechan. Appl. - 1 February, 1999. - Vol. 263, No. 1-4. - P. 4-25.
16. DHAR D., SADHU T., CHANDRA S. Pattern formation in growing sandpiles // Eur. Phys. Lett. - 2009. - Vol. 85, No. 4. -48002. 2009. arXiv:0808.1732 [cond-mat.stat-mech]
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.
RESOURCE NETWORK WITH LIMITED CAPACITY OF ATTRACTOR VERTICES
Ludmila Zhilyakova, Institute of Control Sciences of RAS, (Moscow, Profsoyuznaya st., 65), Doctor of Science, Leader Researcher, zhilyakova.ludmila@gmail.com.
Abstract: An asymmetric regular resource network with limited capacity of attractor vertices is considered. Provided attractors' capacity is unlimited they accumulate all the resource of a network, except some small volumes left in the rest vertices. On the other hand, attractors of limited capacity are saturated, and then another set of vertices starts attracting the rest resource. We refer to this phenomenon as to secondary attractivity, and such vertices, respectively, are called secondary attractors. We study the flows in a network with different total amount resource and search the threshold values which lead to a new behavior of a network.
Keywords: resource network, graph dynamic threshold model, attractors, capacity limitation.
Статья представлена к публикации членом редакционной коллегии П.Ю. Чеботарёвым
Поступила в редакцию 20.08.2015.
Опубликована 30.11.2015.
