ИССЛЕДОВАНИЕ ПОЛНЫХ ОДНОРОДНЫХ РЕСУРСНЫХ СЕТЕЙ С «ЖАДНЫМИ» ВЕРШИНАМИ1
Жилякова Л. Ю.2, Чаплинская Н. В.3
(ФГБУН Институт проблем управления им. В.А. Трапезникова РАН, Москва)
Описана модификация графовой динамической модели «ресурсная сеть» - «ресурсная сеть с жадными вершинами». В этой модели вершины графа на каждом такте дискретного времени обмениваются ресурсом по ребрам, обладающим ограниченной пропускной способностью, причем сначала передают имеющийся ресурс себе в петлю, а оставшийся ресурс распределяют в смежные вершины по правилам «стандартной» ресурсной сети. Это два правила с пороговым переключением: если ресурс вершины превышает суммарную пропускную способность всех ее исходящих ребер, она отдает по полной пропускной способности в каждое ребро, в противном случае отдает весь ресурс, деля его пропорционально пропускным способностям исходящих ребер. Исследуется процесс функционирования полной однородной ресурсной сети с «жадными» вершинами при разной величине суммарного ресурса и различных начальных состояниях. Описаны возможные состояния сети; выявлено нехарактерное для стандартной модели состояние - остановка сети. Найдены два пороговых значения суммарного ресурса, разделяющих зоны различного поведения сети: первое разделяет зоны «недостаточного» и «достаточного» ресурса, второе -зоны «достаточного малого» и «достаточного большого» ресурса. В каждой зоне описано функционирование сети и исследованы предельные состояния и потоки. Для всех характерных ситуаций приведены примеры, демонстрирующие согласованность аналитических результатов с численными экспериментами.
Ключевые слова: ресурсная сеть, модель «жадных» вершин, графовая динамическая пороговая модель.
1. Введение
Ресурсная сеть - графовая динамическая потоковая модель, в которой в дискретном времени происходит перераспределение ресурса между вершинами. Впервые модель была предложена
1 Работа выполнена при частичной финансовой поддержке РФФИ, гранты № 20-07-00190А, 19-07-00525А.
2 Людмила Юрьевна Жилякова, д.ф.-м.н., в.н.с., (zhilyakova.ludmila@gmail.com).
3Надежда Васильевна Чаплинская, техник (nadya1462@gmail.com).
в [11] и получила дальнейшее развитие в ряде работ, результаты которых аккумулированы в [9, 10]. Ресурсная сеть представляет собой ориентированный взвешенный граф. Веса ребер обозначают их пропускные способности. Вершины распределяют ресурс по двум правилам с пороговым переключением, которое происходит в зависимости от количества ресурса, находящегося в данный момент в вершине. Если это количество превосходит суммарный вес исходящих из вершины ребер, то в каждую смежную вершину передается ресурс, равный пропускной способности соответствующего ребра. Иначе вершина отдает весь свой ресурс смежным вершинам пропорционально пропускным способностям соединяющих с ней ребер. На каждом такте дискретного времени £ вершины, имеющие ресурс, отдают его одновременно.
Для описанной модели получены результаты, позволяющие для каждой топологии сети и пропускных способностей ее ребер, для любого суммарного ресурса и его начального распределения по вершинам определить предельный поток и предельное состояние вершин сети. Также определена формула, позволяющая для произвольной сети вычислить пороговое значение ресурса Т, при переходе через которое сеть изменяет свое поведение (см. [9]).
При суммарном ресурсе Ш, не превосходящем Т, т.е. при малом ресурсе, все вершины сети функционируют по одному правилу и отдают весь ресурс на каждом такте. Динамика сети в этом случае описывается моделью рассеяния на графах [16, 21]. Модели рассеяния хорошо изучены - для графов с постоянной топологией они описываются однородными цепями Маркова. Аппарат однородных цепей Маркова применяется при моделировании процессов во многих предметных областях, в частности при решении задач о нахождении консенсуса в многоагентных системах [4, 5], социальных сетях [7] и многих других.
При суммарном ресурсе выше порогового значения, т.е. при большом ресурсе, процесс функционирования сети описывается неоднородной цепью Маркова. Некоторые вершины сети в этом случае накапливают излишки ресурса (такие вершины были названы аттракторами [9]). Вершины, имеющие большой ресурс,
ведут себя как «стреляющие вершины» в пороговой модели «игра выстреливания фишек» (chip-firing game) [15, 19, 22]. Эта игра на графе стала широко известна благодаря тому, что нашла применение в описании схода лавин разной природы - она легла в основу математической модели «абелева куча песка» [20, 24], которая, в свою очередь, применяется в моделировании сложных природных и социальных процессов, названных «самоорганизованной критичностью» [13, 14, 17, 18, 23].
За время с момента появления модели «ресурсная сеть» в 2009 году, был создан и исследован ряд ее модификаций. В монографии [10] введена модель с ограничениями на ёмкость вершин-аттракторов. С помощью этого ограничения некоторые не-аттрактивные вершины сети получили возможность также накопить некоторые излишки ресурса. Для таких сетей было введено понятие «вторичной аттрактивности».
В работах [1,2] были рассмотрены ресурсные сети на графах с нестандартной достижимостью [8]. В [1] введены и описаны ресурсные сети с вентильной достижимостью, в [2] - с магнитной достижимостью. В статье [3] предложены двухресурсные сети и, как и в самой первой работе по ресурсным сетям [11], рассмотрен частный случай: полные однородные сети. Динамические ресурсные сети (сети, веса дуг которых могут принимать более одного значения) предложены и описаны в [12].
Еще один научный коллектив исследовал применение ресурсных сетей с нестандартной достижимостью при построении моделей телекоммуникационных систем [6].
Настоящая статья посвящена новой модификации ресурсных сетей - ресурсной сети с «жадными» вершинами. Отличие этой модели заключается в том, что каждая вершина передает ресурс сначала в свою петлю (если она есть) и только затем остатки распределяет в другие вершины. Как и в других случаях, мы начинаем исследование с частного случая: полной однородной сети с петлями.
2. Ресурсные сети - основные понятия
Формальное описание ресурсной сети с «жадными» вершинами наследует все понятия обычной ресурсной сети. Дадим основные определения и опишем функционирование «стандартной» сети.
Ресурсной сетью называется связный ориентированный граф С = (У,Е) с вершинами уг € V и ребрами е^ = ) € Е; IV| = п. Ребра графа имеют неотрицательные веса г^, обозначающие их пропускные способности; К = (г^)пхп - матрица пропускных способностей. Вершинам графа приписываются неотрицательные числа дг(Ь), изменяющиеся в дискретном времени £ и называемые ресурсами. Вершины могут хранить неограниченное количество ресурса.
Состоянием Q(t) сети в момент времени £ называется вектор <3(1) = (д\ (I), ..., дп(£)), содержащий значения ресурсов в вершинах в момент 1
Состояние * = (д*,..., д*п) называется предельным, если для любого е > 0 существует 1£ такое, что для всех 1>1е | д* -дг(¿) | < е, г = 1,2,...,п.
Обозначив суммарный ресурс всех вершин через Ш, можно
п
записать закон сохранения, выполняющийся в сети: V .
г=1
п п
ГТ = г]г и г°иЬ = ^ г^ - входная и выходная пропуск-3 = 1 3=1
ные способности вершины соответственно. Пропускная способность петли, если она существует, входит в обе суммы. Суммар-
п п
ная пропускная способность сети: г8ит = ^ ^ г^.
г=1з=\
В каждый момент £ вершины ьг, г = 1, ...,п, передают ресурс по выходящим ребрам. Стандартная модель ресурсной сети характеризуется следующими правилами:
В момент £ вершина уг отдает в ребро е^, соединяющее ее с вершиной Vj:
• гу единиц ресурса, если д^) > г°иЬ (правило 1);
• зётЧ^) единиц ресурса, если ^ (¿) ^ г^ (правило 2).
г
Смысл этих правил заключается в следующем. Правило 1 работает тогда, когда вершина содержит ресурса больше, чем может отдать, т.е. по каждому ребру отдается ресурс, равный про-
п
пускной способности этого ребра, а всего отдается г([и = ^ г^
3 = 1
ресурса. По правилу 2 вершина отдает весь свой ресурс, причем он распределяется пропорционально пропускным способностям ребер. Если ресурс в вершине равен выходной пропускной способности вершины: qí(t) = г^, то применение правил 1 и 2 даст одинаковые результаты.
Множество вершин, для которых qí(t) ^ г([и1, называется зоной 2-(£); множество вершин, для которых > г™1, называется зоной 2 + (¿).
Распределение ресурса происходит параллельно: на каждом такте все вершины, имеющие ресурс, отдают его во все исходящие ребра по правилу 1 или 2, в зависимости от его количества; вершины, смежные с ними, получают этот ресурс по входящим ребрам на следующем такте.
Важной характеристикой ресурсной сети является пороговое значение ресурса Т: при Ш ^ Т все вершины, начиная с некоторого переходят в зону 2-(£) (функционируют по правилу 2); при Ш > Т зона 2+(£) непуста и не изменяется, начиная с некоторого Ь" (существуют вершины, функционирующие по правилу 1). В [9] доказано, что в любой сети, кроме поглощающей, пороговое значение существует и единственно, а также предложены формулы, позволяющие найти значение Т для всех классов сетей, в которых оно существует.
Поток ресурса
Ресурс, выходящий из вершины Ví по ребру е^ в момент ¿, приходит в вершину Vj в момент £ + 1. Между моментами £ и £ + 1 ресурс находится в ребре е^. Этот ресурс называет-
ся потоком (¿). Общий поток Р(Ь) описывается матрицей
п
Величина = ^ ¡^ (I) - ресурс, выходящий из вер-
3=1
шины Уг в момент
п
Величина Дп(1 + 1) = ^ }гз(^) - ресурс, входящий в вер-
г=1
шину Уу в момент £ + 1. Полагаем, что /™(0) = 0.
3. Модель сети с «жадными» вершинами
3.1. ОПИСАНИЕ И ОСНОВНЫЕ ОПРЕДЕЛЕНИЯ
Рассмотрим ресурсную сеть с наличием хотя бы одной петли. Ее функционирование определяется следующим образом: Вершина г>г, имеющая петлю, на такте Ь отдает в петлю:
• Гц ресурса, если д^) > Гц,
• весь свой ресурс ^(¿), если дг(1) ^ Гц.
При этом в первом случае оставшийся ресурс распределяется между другими вершинами по правилам стандартной ресурсной сети.
Интерпретировать такую модель можно как ресурсную сеть, у которой есть «запасливые» (или «жадные») вершины, т.е. вершины, передающие имеющийся ресурс в первую очередь самим себе.
Определение 1. Будем говорить, что сеть остановилась на такте I, если У г = 1,п д^) ^ Гц.
Это определение говорит о том, что если количество ресурса во всех вершинах не превосходит пропускных способностей петель, то ресурс не перераспределяется между вершинами. Заметим, что остановка сети является частным случаем ее стабилизации.
Определение 2. Назовем вершину Ví насыщенной в момент времени ^ если д^) ^ г и, перенасыщенной, если д^) > г и, и ненасыщенной, если д,,(Ь) < г и.
Утверждение 1. Вершина, насыщенная в момент времени I, в момент времени £ + 1 также останется насыщенной.
Доказательство. Вершина такова, что ^ ги. На такте £ + 1 ресурс равный г^ пойдет в петлю, а оставшийся ресурс (%(£) — гы) распределится между остальными вершинами по правилам стандартной ресурсной сети. Таким образом, на такте £ +1 имеем д^ + 1) ^ г и.
Для перенасыщенной вершины введем понятие свободного ресурса.
Определение 3. Свободный ресурс перенасыщенной вершины Ут, -это ресурс, не задействованный петлей: = gí (Ь)—Гп. Пусть А - множество индексов насыщенных вершин:
(1) А ={к е {1, ...,п} | дк(г) > гиу
Суммарный свободный ресурс в сети на такте Ь обозначим
Дд(£), ДЯ(1) = £ Д^).
íeА
3.2. СВОЙСТВА ПОЛНЫХ ОДНОРОДНЫХ СЕТЕЙ С «ЖАДНЫМИ ВЕРШИНАМИ»
В настоящем исследовании ограничимся исследованием полных однородных ресурсных сетей.
Ресурсная сеть называется полной, если г^ = 0 = 1,п. Однородность сети означает, что г^ = г Уг,] = 1,п.
Утверждение 2. В полной однородной сети с «жадными» вершинами суммарное количество ресурса, полученное вершинами V^ Vj от остальных вершин ук, к е {1, ...,п}\{1, ]}, одинаково на каждом такте I.
Доказательство. Мы рассматриваем полную однородную сеть, а значит, и вершина Ví, и вершина Vj обязательно имеют ребра, соединяющие их с каждой из оставшихся вершин Ьк, к е {1,...,n}\{i,j }.
Суммарный ресурс, получаемый вершиной ы от вершин ,
к Е {1,..., и}\(г, ]}, на такте £ равен
,{1..п}\{ч} = у^ А дк ) ] * ^ п - 1 ,
к^А\{г,]}
где А - множество насыщенных вершин (формула (1)). Аналогично для вершины Vj:
Л1..п}\{г,Л = у^ АЯк ^) ] з ^ п - 1 .
к^А\{г,]}
Таким образом, /{^М = .
3.3. ПОРОГОВЫЕ ЗНАЧЕНИЯ РЕСУРСА
В стандартной ресурсной сети существует одно пороговое значение Т, разделяющее ресурс на малый, когда все вершины на каждом такте отдают весь свой ресурс, и большой, когда излишки ресурса накапливаются в одной или нескольких вершинах (см. раздел 2).
Поведение сети с «жадными» вершинами сложнее. Определение 1 говорит о том, что существуют ситуации, в которых вершины перестают обмениваться ресурсом (сеть останавливается). В стандартных эргодических сетях такая ситуация невозможна.
В ресурсных сетях с «жадными» вершинами существуют два пороговых значения суммарного ресурса Т1 и Т2, которые разделяют зоны различного поведения сети. Если Ш ^ Т сеть останавливается - за конечное число тактов, или асимптотически, в зависимости от начального состояния.
Определение 4. Суммарное значение ресурса, при котором функционирование сети останавливается, будем называть недостаточным.
Максимальное значение ресурса, при котором вершины могут не обмениваться ресурсом, равно т. Если в начальном состоянии этот ресурс распределить равномерно по п вершинам, сеть
«остановится» на нулевом такте. Будем полагать, что для некоторых начальных состояний Т\ = гп. В следующем разделе будет доказано, что это значение единственно для любого начального состояния сети.
Если Ш > гп, на каждом такте будет существовать хотя бы одна вершина с ресурсом, превосходящим пропускную способность петли г. Это означает, что такая сеть никогда не остановится.
Определение 5. Суммарное значение ресурса, при котором функционирование сети не останавливается, будем называть достаточным.
Достаточный ресурс, по аналогии со стандартной сетью, делится порогом Т2 на малый (Т\ < Ш ^ Т2) и большой ресурс ^ > Т2). Будет показано, что величина Т2 совпадает с пороговым значением Т = гп2 стандартной полной однородной сети.
Рис. 1. Схема пороговых значений суммарного ресурса в сети с «жадными» вершинами
3.4. НЕДОСТАТОЧНЫЙ РЕСУРС
Рассмотрим функционирование сети при Ш ^ гп. В этом случае общее количество ресурса не превосходит сумму пропускных способностей петель.
Нетрудно убедиться, что в этом случае предельное состояние зависит от начального распределения ресурса. Например, если ресурс настолько мал, что в каждой вершине не превосходит г, то предельное состояние совпадает с начальным. Сформулируем
ряд утверждений о возможных поведениях сети при недостаточном ресурсе.
Пусть ф(0) - вектор начального состояния сети с п вершинами. Упорядочим вектор Q(0) по убыванию его координат и в соответствии с этим перенумеруем вершины. Тогда вершина VI хранит наибольшее количество ресурса, а вершина уп - наименьшее.
Утверждение 3. (Тривиальный случай.) Если ^(0) ^ г, то сеть остановится на такте £ = 0. Начальное состояние в данном случае будет предельным.
Доказательство. д1 (0) ^ г ^ Уг = 2,п ^(0) ^ д1(0) ^ г, а, значит, по определению, сеть остановилась на такте £ = 0.
Утверждение 4. (Достаточное условие бесконечного функционирования сети.) Пусть Ш ^ гп. Если д1(0) > г и д2(0) ^ г, то сеть будет функционировать бесконечно, останавливаясь в пределе.
Доказательство утверждения приведено в приложении.
Замечание 1. Для каждой насыщенной вершины сходимость к предельному состоянию, в котором она остановится, имеет характер затухающих колебаний.
Рассмотрим частный случай, когда ресурс сосредоточен только в двух первых вершинах. Для этого случая можно найти вектор предельного состояния в явном виде.
Утверждение 5. Пусть в полной однородной сети с «жадными» вершинами Ш ^ гп; д1(0) >г и д2(0) ^ г, д,,(0) = 0, г = 3, п. Тогда ее предельное состояние задается вектором
Ш — 2г № — 2г-
п — 2 , п — 2
Доказательство. По утверждению 2, суммарный ресурс, полученный вершинами VI и г>2 от остальных вершин Ví, г = 3,п, одинаков. Обозначим его через -Ыге3¿(0).
По условию -Ыгеst(0) = 0.
Докажем, что У£ wrest(t) = 0.
В самом деле, если •ыге3¿(0) = 0, это означает, что все вершины Ví, г = 3,п, получили ресурс, больший г, - это напрямую
/ УУ — 2г УУ — 2г \
(2) Я = [г, г, -— , ..., -— .
V п—2 п—2/
следует из утверждения 2. Вершины ы, г = 3,п, получают попарно одинаковый ресурс, а значит, их входной поток одинаков, в силу полноты сети и их симметричности относительно вершин у1 и у2. Но если для некоторого значения £ д^) > г, г = 3,п, то неравенство Ш ^ гп не может быть выполненным: в силу утверждения 1 вершины г>1 и г>2 также останутся насыщенными, и суммарный ресурс превысит гп. Таким образом, Ш ,Гез^) = 0. Тогда динамика ресурса в первых двух вершинах будет следующей:
д1(1) = г + Я1(2) = г +
= г + если - четное, и
= г +
2(0) -
п - 1
(Ц(0) -
(п - 1)2
(Ц(0) -
(п - 1)
(12(0) -
(п - 1)
Я2(1) = Г + я2(2) = г +
2( ) = + 2( ) = +
91(0) - г
п - 1
(12(0) -
(п - 1)2
(12(0) -
(п - 1) '
(Ц(0) -
(п - 1) '
если £ - нечетное.
Пределом последовательности состояний вершин 1 и 2 при £ ^ те будет г. По закону сохранения общего ресурса Ш в сети на остальные вершины в пределе останется ( Ш - 2г) ресурса,
• о--Ш-2г
т.е. на каждую вершину ы, г = 3,п, придется п_2 ресурса:
Ш -2г
Ш < гп ^-< г.
п - 2
Таким образом, предельное состояние описывается вектором (2).
Пример 1. В данном примере и в последующих будем рассматривать функционирование однородной ресурсной сети с пятью вершинами и единичными пропускными способностями:
г = 1.
Рис. 2. Однородная ресурсная сеть с пятью вершинами, г = 1
Начальное состояние задано вектором: ф(0) = (1,7; 1; 0,9; 0,6; 0,4). Суммарный ресурс ^^ — 4, 6 < тп — 5, что удовлетворяет условию «недостаточности» ресурса. При этом ^(0) = 1, 7 > г и ?2(0) = 1 = г, а, значит, согласно утверждению 4, сеть будет функционировать бесконечно, останавливаясь в пределе.
Ресурс насыщенных вершин г>1, ь2 сходится к предельному состоянию в виде затухающих колебаний (рис. 3, таблица 1). Вершина г>э также становится насыщенной на втором такте, поэтому и сходимость ее ресурса к предельному состоянию также имеет характер затухающих колебаний. Вершинам г>4, г>5 для насыщения не хватает начального ресурса; в этих вершинах ресурс увеличивается монотонно и в пределе не превосходит пропускной способности петли.•
Рассмотрим функционирование сети в нетривиальном случае (91(0) > г) с произвольным начальным состоянием.
Утверждение 6. Пусть Ш ^ гп; д1 (0) > г. Тогда если д2(0) ^ гп~пП"-/0, то сеть остановится на первом такте с со-
1,8
0,2
О--1-1-1-1-1-1-1-1-1
123456789
Рис. 3. Динамика ресурса в сети с начальным состоянием Q(0) = (1,7; 1; 0,9; 0,6; 0,4)
Таблица 1. Протокол работы сети с начальным состоянием ф(0) = (1,7; 1; 0,9; 0,6; 0,4)_
¿г VI У2 Уз У4 У5
0 1,7000 1,0000 0,9000 0,6000 0,4000
1 1,0000 1,1750 1,0750 0,7750 0,5750
2 1,0625 1,0188 1,0437 0,8375 0,6375
3 1,0156 1,0266 1,0203 0,8687 0,6687
10 1,0002 1,0002 1,0002 0,8998 0,6998
11 1,0001 1,0001 1,0001 0,8999 0,6999
12 1,0000 1,0000 1,0000 0,8999 0,6999
13 1,0000 1,0000 1,0000 0,9000 0,7000
14 1,0000 1,0000 1,0000 0,9000 0,7000
стоянием Я* = (V, ^(0) + , ..., 9п(0) + ). Иначе
сеть будет функционировать бесконечно, останавливаясь в пределе.
Доказательство. На первом такте вершина г>1 отдает в собственную петлю полную ее пропускную способность г, в осталь-
q1(0)-r
ные вершины она отдает п_ 1 ресурса.
Пусть для вершины выполняется условие д2(0) < г, так как в противном случае, согласно утверждению 4, сеть функционирует бесконечно и останавливается в пределе. В таком случае вершина ь2 на первом такте имеет состояние д2 (1) = = 92(0) + .
Теперь если д2(1) ^ г, т.е. д2(0) ^ гпП~11(0), то, по определению, сеть остановится на такте Ь = 1, так как во всех вершинах ресурса не больше, чем пропускная способность петли. Заметим, что д2 (0) ^ гп11(0) < ^П—т = г, т.е. условие д2(0) < г заведомо
г—г п
выполнено.
Если же д2(0) > гп~—_ 11(0) (а это то же самое, что д2(1) > г), то, получив на первом такте 51(1) = г и перенумеровав вершины в порядке убывания ресурса (теперь д1(1) > г и д2(1) = г), имеем ситуацию, удовлетворяющую достаточному условию бесконечного функционирования сети и ее остановки в пределе (утверждение 4).
Таким образом, если д2(0) ^ тп111(0), то сеть остановится на первом такте с состоянием
О* = (г, .2(0) + ^М^.....,п(0) + ^^).
V п — 1 п — 1/
Если д2 (0) > 1пП-\(0), то сеть будет функционировать бесконечно, останавливаясь в пределе.
Приведем в пример, иллюстрирующий утверждение 6.
Пример 2. Рассмотрим полную однородную сеть с пятью «жадными» вершинами, г = 1.
1. Рассмотрим динамику распределения ресурса между вер-
шинами при следующем начальном состоянии:
Q(0) = (3; 0,5; 0,3; 0,2; 0).
Здесь суммарный ресурс Ш = 4, что удовлетворяет условию «недостаточного» ресурса: Ш ^ гп = 5, где п = 5, г = 1. При этом 91(0) = 3 > г и 92(0) = 0,5 < гп~-11(0) = = 0,5, а, значит, согласно утверждению 6, сеть остановится на первом же такте с состоянием
О* = (г. „2(0) + ^, 03(0) + А?,
«.(0) + , ,5(0) + ) =
п — 1 п — 1 /
/ 3 — 1 3 — 1 3 — 1 3 — 1N
= (^о^ + д—[ ;0,3 + 3—[ 1 ;0+5—г) =
= (1; 1; 0,8; 0,7; 0,5).
Результаты программной реализации процесса функционирования сети приведены в таблице 2.
Таблица 2. Протокол работы сети с начальным состоянием д(0) = (3; 0,5; 0,3; 0,2; 0)
и У2 Уз У4 У5
0 3,0000 0,5000 0,3000 0,2000 0,0000
1 1,0000 1,0000 0,8000 0,7000 0,5000
2 1,0000 1,0000 0,8000 0,7000 0,5000
2. Рассмотрим динамику распределения ресурса между вершинами при другом начальном состоянии:
д(0) = (3;1;0;0;0).
Здесь суммарный ресурс также Ш = 4. При этом д1 (0) = 3 > г и 92(0) = 1 > гпП!\(0) = 5-3 = 0,5, а, значит, согласно утверждению 6, сеть будет функционировать бесконечно, останавливаясь
19
в пределе. При этом, так как д3(0) = д4(0) = д5(0) = 0, выполняются условия утверждения 5, и предельное состояние вычисляется по формуле (2). Имеем:
№ — 2г 4 — 2
-=-= 0,6667.
п — 2 5 — 2
Тогда
д* = (1; 1; 0,6667; 0,6667; 0,6667).
Результаты программной реализации процесса функционирования сети представлены в таблице 3 и на рис. 4.^
Таблица 3. Протокол работы сети с начальным состоянием д(0) = (3;1;0;0;0)_
и У1 У2 Уз У4 У5
0 3,0000 1,0000 0,0000 0,0000 0,0000
1 1,0000 1,5000 0,5000 0,5000 0,5000
2 1,1250 1,0000 0,6250 0,6250 0,6250
3 1,0000 1,0313 0,6563 0,6563 0,6563
4 1,0078 1,0000 0,6641 0,6641 0,6641
5 1,0000 1,0020 0,6660 0,6660 0,6660
6 1,0005 1,0000 0,6665 0,6665 0,6665
7 1,0000 1,0001 0,6666 0,6666 0,6666
8 1,0000 1,0000 0,6667 0,6667 0,6667
9 1,0000 1,0000 0,6667 0,6667 0,6667
Теорема 1. Первое пороговое значение ресурса Т1 полной однородной сети с «жадными» вершинами, отделяющее недостаточный ресурс от достаточного, единственно и равно гп для любого начального состояния.
Доказательство. Утверждения 3 и 6 рассматривают все случаи распределения суммарного ресурса, не превосходящего гп. В каждом из этих случаев предельное состояние существует и сеть в нем останавливается.
Рис. 4. Динамика ресурса в сети с начальным состоянием д(0) = (3;1;0;0;0)
При суммарном ресурсе, строго большем гп, на каждом такте £ хотя бы одна вершина будет иметь ресурс, больший г, а это означает, что Т\ = гп.
3.5. ДОСТАТОЧНЫЙ РЕСУРС
Рассмотрим функционирование сети при Ш > Т\ = гп. В этом случае сеть не останавливается.
Утверждение 7. При суммарном ресурсе Ш > Т\ все вершины полной однородной сети с «жадными» вершинами насыщаются за конечное число тактов.
Доказательство см. в приложении.
Поскольку в сети с «жадными» вершинами при достаточном ресурсе через конечное число тактов все вершины насыщаются, после насыщения ресурс, равный гп, находящийся в петлях, не участвует в потоке между вершинами. Рассмотрим сеть без петель с суммарным ресурсом Ш — гп. Полученная сеть в таком случае будет функционировать как стандартная модель с «нежадными» вершинами.
Утверждение 8. Полная однородная ресурсная сеть c «жадными» вершинами при Ш > гп с момента насыщения всех вершин функционирует эквивалентно соответствующей
21
стандартной ресурсной сети без петель, у которой ресурс в каждой вершине уменьшен на величину пропускной способности петли г.
Замечание 2. Здесь и далее будем говорить о сетях с количеством вершин п > 2, так как в случае п = 2 при переходе к эквивалентной стандартной сети, «отбросив» петли, мы получим нерегулярную сеть - эта сеть будет циклической с периодом 2; при малом ресурсе она не имеет предельного состояния.
В силу утверждения 8, на полную однородную сеть с «жадными» вершинами при Ш > Т! переносятся свойства полной однородной стандартной ресурсной сети без петель.
Рассмотрим свойства полных однородных стандартных ресурсных сетей, описанные в [11], которые верны и для сетей без петель.
Свойство 1. Если для некоторого выполняется равенство Яг(¿0 = 1] (£'), то для всех также выполняется = ^ (Ъ).
Это следует из того, что с момента Ь обе вершины получают одинаковое количество ресурса.
Свойство 2. Если для некоторого вершина уг находится в зоне 2— (£'), то для всех £ > 1? вершина уг будет находиться в зоне 2 — (¿).
Это следует из того, что г>г в момент времени отдает весь свой ресурс, а получить может не более чем г(п — 1). То есть Яг(¿' + 1) ^ г(п — 1) = г™1, значит, вершина не может покинуть зону 2— (¿).
Свойство 3. Если все вершины находятся в зоне 2 + (£), то состояние устойчиво.
Это следует из того, что все вершины получают и отдают г(п — 1) ресурса.
Однако, ввиду отсутствия петель, следующие два свойства полной однородной сети не могут быть применимы к нашей модели.
Свойство 4. Для любых £ входные потоки всех вершин совпадают: ¡¡п(Ь) = ¡р(Ь), 1,] = 1,..., п.
Свойство 5. Если в момент £ вершины г^, ...,Угт (т ^ п) находятся в зоне Z-(1), то (£ + 1) = ... = (£ + 1).
Для проверки свойства 4 рассмотрим две любые вершины, например, и у2.
В силу однородности сети, ¡ц{Ъ) = ... = г = 1,п.
Расписав входные потоки вершин ы и г>2, получим:
ЯП^) = /21(4) + Ь (*) + ... + /ш(*),
№(!) = ¡12(1) + /32 (*) + ... + /п2^), так как У^ц= /22^) = 0 (петли отсутствуют).
лп(1) — ЯТ(1) = /21(4) — /12 Ю,
что в общем случае не равно нулю, а, значит, свойство 4 не выполняется для полных однородных сетей без петель.
Поскольку вершины зоны ^- (¿) отдают весь свой ресурс, то их ресурс в момент £ + 1 равен поступающему к ним входному потоку. Поэтому свойство 5 в общем случае также не выполняется.
Непосредственно из утверждения 8 следует центральный результат этого раздела.
Теорема 2. Второе пороговое значение полной однородной сети с «жадными» вершинами Т2 совпадает с пороговым значением соответствующей стандартной сети с петлями: Т2 = Т = гп2.
В самом деле, сеть с «жадными» вершинами при Ш > гп ведет себя так же, как стандартная сеть без петель с ресурсом Ш — гп. Ее пороговое значение Т = гп2 — гп (Теорема 2.1 [10]). Вернув назад петли вместе с циркулирующим в них ресурсом, равным гп, получим: Т2 = Т = гп2.
3.5.1. Малый ресурс
Рассмотрим функционирование сети при малом ресурсе: Т1 < Т2.
Случаи Ш < Т2 и Ш = Т2 рассмотрим отдельно. Для описания процесса распределения ресурса между вершинами удобно говорить о ресурсной сети как о потоковой модели.
Введем два вектора:
оиЩ) = и<Т\1), ..., п\1)), т(Ш) = (/1п(^+1), ..., С(Ш)).
Тогда процесс функционирования сети можно записать в следующем виде:
Таким образом, переход из состояния в состояние можно описать в терминах матрицы Р (¿).
Следующие утверждения доказываются в терминах потоков. Доказательства приведены в приложении.
Утверждение 9. При Ш < Т2 вне зависимости от начального состояния Q(0) все вершины полной однородной ресурсной сети с «жадными» вершинами переходят в зону 2 — (¿) за конечное число тактов.
Утверждение 10. Полная однородная ресурсная сеть с «жадными» вершинами при Ш < Т2 с момента насыщения всех вершин функционирует согласно формуле
Теорема 3. Полная однородная ресурсная сеть с «жадными» вершинами при Ш < Т2 вне зависимости от начального состояния Q(0) стабилизируется с предельным состоянием
Пример 3. Рассмотрим однородную ресурсную сеть с пятью «жадными» вершинами, г = 1 (см. рис. 1).
+1) = Я(г) — оиЩ) + т(г +1).
(3) <3(1 + 1) = д(1)в!,в!
0 1 ... 1
1 1 0 ... 1
п — 1
1 1 ... 0
Рассмотрим динамику распределения ресурса между вершинами при начальном состоянии, подобранном таким образом, чтобы все вершины были уже насыщенными:
д(0) = (3;2;2;1;1).
Суммарный ресурс Ш = 9 удовлетворяет условию гп <Ш ^ ^ гп2, а значит, согласно утверждению 12, сеть стабилизируется
/ Ш Ш \
с предельным состоянием О* = —, ..., — = (1,8; ...;1,8).
\ п п )
Протокол работы сети приведен в таблице 4.
Таблица 4. Протокол работы сети с начальным состоянием Q(0) = (3; 2; 2; 1; 1) (модель «жадных» вершин)
и У1 У2 Уз У 4 У5
0 3,0000 2,0000 2,0000 1,0000 1,0000
1 1,5000 1,7500 1,7500 2,0000 2,0000
2 1,8750 1,8125 1,8125 1,7500 1,7500
3 1,7813 1,7969 1,7969 1,8125 1,8125
4 1,8047 1,8008 1,8008 1,7969 1,7969
5 1,7988 1,7998 1,7998 1,8008 1,8008
6 1,8003 1,8000 1,8000 1,7998 1,7998
7 1,7999 1,8000 1,8000 1,8000 1,8000
8 1,8000 1,8000 1,8000 1,8000 1,8000
9 1,8000 1,8000 1,8000 1,8000 1,8000
Теперь рассмотрим функционирование соответствующей стандартной сети. Здесь начальное состояние не будет включать ресурс, задействованный петлями:
д(0) = (2;1;1;0;0).
Заметим, что для стандартной сети матрица пропускных способностей будет иной, так как мы не учитываем петли: г^ = 1, если г = ], и гу = 0, если г = ].
Таблица 5. Протокол работы сети с начальным состоянием
3(0
(2; 1; 1; 0; 0) (стандартная модель)
¿г VI У2 У3 У4 У5
0 2,0000 1,0000 1,0000 0,0000 0,0000
1 0,5000 0,7500 0,7500 1,0000 1,0000
2 0,8750 0,8125 0,8125 0,7500 0,7500
3 0,7813 0,7969 0,7969 0,8125 0,8125
4 0,8047 0,8008 0,8008 0,7969 0,7969
5 0,7988 0,7998 0,7998 0,8008 0,8008
6 0,8003 0,8000 0,8000 0,7998 0,7998
7 0,7999 0,8000 0,8000 0,8000 0,8000
8 0,8000 0,8000 0,8000 0,8000 0,8000
9 0,8000 0,8000 0,8000 0,8000 0,8000
Рисунки 5 и 6 позволяют сравнить функционирование сети с «жадными» вершинами и соответствующей ей стандартной сети без петель.»
щ
I
Рис. 5. Динамика ресурса в сети с начальным состоянием Q(0) = (3; 2; 2; 1; 1) (модель «жадных» вершин)
Рассмотрим случай Ш = гп2. Для него утверждение, что все вершины сети за конечное число тактов перейдут в зону Z- (¿) неверно: ресурс хотя бы одной из вершин в этой ситуации будет 26
Рис. 6. Динамика ресурса в сети с начальным состоянием Q(0) = (2; 1; 1; 0; 0) (стандартная модель)
стремиться к го сверху. Поэтому формула для вычисления следующего состояния сети по предыдущему также неприменима: хотя бы одна из вершин будет функционировать по правилу 1. Однако можно показать, что и в таком случае зоны «достаточного малого» ресурса сеть стабилизируется с тем же самым предельным состоянием:
/ W \
Q*(t) = —, ..., — = (го, ..., го).
\п п ;
Теорема 4. Полная однородная ресурсная сеть c «жадными» вершинами при W = го2 вне зависимости от начального состояния ф(0) стабилизируется с предельным состоянием
Q*(t) = (го, ..., го).
3.5.2. Большой ресурс
Рассмотрим функционирование сети при Ш > Т2 = го2. «Большой ресурс» означает, что при любом начальном состоянии Ф(0) в любой момент времени зона Z + (£) будет непуста.
Функционирование сети с «жадными вершинами» при большом ресурсе наиболее близко к функционированию соответству-
ющей стандартной сети без петель.
Свойства стандартных полных однородных ресурсных сетей без петель при большом ресурсе выполняются и для модели «жадных» вершин.
Утверждение 11. Если все вершины находятся в зоне 2+(£), то состояние Q(t) устойчиво.
Доказательство. Справедливость этого свойства вытекает из того факта, что входной и выходной потоки каждой вершины одинаковы и равны гп.
Утверждение 12. При большом суммарном ресурсе Ш >Т2 для произвольного начального состояния Q(0) предельное состояние существует.
Утверждение 13. Пусть при большом суммарном ресурсе Ш > Т2 для произвольного начального состояния Q(0) ресурс в вершинах упорядочен по убыванию и первые к вершин находятся в зоне 2 + (0), а остальные п — к вершин - в зоне 2-(0):
91(0) ^ ...дк(0) >гп > Як+1(0) ^ ЯП(0).
Тогда вектор предельного состояния будет иметь вид
= (гп + 5д1,..., гп + 5д1, гп, ...гп),
где I ^ к; Ьдг,% = 1, ...,1 - излишки ресурса вершин, оставшихся в зоне 2+*.
4. Заключение
Предложена и формально описана неклассическая модель ресурсной сети - сеть с «жадными» вершинами. «Жадность» вершин заключается в том, что на каждом такте они передают ресурс сначала в свою петлю, а после этого остаток распределяют в исходящие рёбра. Была исследована сеть, представленная полным однородным графом (полным графом, у которого веса ребер одинаковы).
Описано функционирование сети при любом количестве ре-
сурса, доказана сходимость процессов перераспределения ресурса, найдены предельные состояния. Показано, что при некоторых значениях суммарного ресурса сеть останавливается. Это особенно примечательно тем, что в стандартной регулярной сети остановка функционирования невозможна.
Найдены пороговые значения Т\ = гп и = гп2, разграничивающие зоны различного поведения сети:
1) «недостаточный» ресурс: Ш ^ Т\,
2) «достаточный» ресурс: Ш > Р\
а) «достаточный малый» ресурс: Т\ ^ Т2,
б) «достаточный большой» ресурс: Ш > Т2.
В дальнейшем планируется изучение сетей с «жадными» вершинами для других топологий; нахождение их основных характеристик, исследование устойчивости и нахождение векторов предельных состояний и потоков для разных величин суммарного ресурса и начальных состояний сети.
Литература
1. АБДУЛРАХМАН Х., ЕРУСАЛИМСКИЙ Я.М., СКОРОХОДОВ В.А. Ресурсные сети с вентильной достижимостью // Инженерный вестник Дона. - 2018. - №4. - 12 с.
2. АБДУЛРАХМАН Х., СКОРОХОДОВ В.А. Ресурсные сети с магнитной достижимостью // Известия ВУЗов. СевероКавказский регион. Естественные науки. - 2016. - №4. -С. 4-10.
3. АБДУЛРАХМАН Х., СКОРОХОДОВ В.А. Полные двух-ресурсные сети с петлями // Известия ВУЗов. СевероКавказский регион. Естественные науки. - 2016. - №2. -С. 10-16.
4. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Сходимость и устойчивость в задачах согласования характеристик (обзор базовых результатов) // Управление большими системами. -2010. - Вып. 30.1 - С. 470-505.
5. АГАЕВ Р.П., ЧЕБОТАРЕВ П.Ю. Метод проекции в задаче о консенсусе и регуляризованный предел степеней стохастической матрицы // Автоматика и телемеханика. - 2011. -№12. - С. 38-59.
6. АНТОНОВА В.М., ЗАХИР Б.М., КУЗНЕЦОВ Н.А. Моделирование графов с различными видами достижимости с помощью языка Python // Информационные процессы. -2019. - Том 19, №2. - С. 159-169.
7. ГУБАНОВ Д.А., НОВИКОВ Д.А., ЧХАРТИШВИЛИ А.Г. Социальные сети. Модели информационного влияния, управления и противоборства. - М.: Изд-во физ.-мат. литры, 2010.-228 с.
8. ЕРУСАЛИМСКИЙ Я.М., СКОРОХОДОВ В.А., КУЗЬМИ-НОВА М.В., ПЕТРОСЯН А.Г. Графы с нестандартной достижимостью: задачи, приложения. - Ростов-на-Дону: Южный федеральный университет, 2009. - 195 с.
9. ЖИЛЯКОВА Л.Ю., КУЗНЕЦОВ О.П. Теория ресурсных сетей: монография.- М.: РИОР: ИНФРА- М, 2017. - 283 с.
10. ЖИЛЯКОВА Л.Ю. Ресурсные сети с ограничениями на ёмкость вершин [Электронный ресурс]: монография. - М.: ИЦ РИОР, НИЦ ИНФРА-М, 2018. - 160 с.
11. КУЗНЕЦОВ О.П. Однородные ресурсные сети. I. Полные графы // Автоматика и телемеханика. - 2009. - №11. -С. 136-147.
12. СКОРОХОДОВ В.А., АБДУЛРАХМАН Х. Динамические ресурсные сети. Случай малого ресурса // Вестник Воронежского государственного университета. Серия: Физика. Математика. - 2018. - №4. - С.186-194.
13. BAK P., TANG C., WIESENFELD K. Self-organized criticality //Phys. Rev. A. -1988. -Vol. 38, No. 1.-P. 364-374.
14. BAK P., CHEN K. Self-organized criticality // Scientif. Am. -1991. - No. 264. - P. 46-53.
15. BJORNER A., LOVASZ L. Chip-firing game on directed graphs //J. Algebraic Combinatorics. - 1992.-No. 1.-P. 305-328
16. 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.
17. 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.
18. 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]
19. GLASS D., KAPLAN N. Chip-Firing Games and Critical Groups // Cornell University 12 Aug 2019 arXiv 1908.04395 [math.CO]
20. IVASHKEVICH E.V., PRIEZZHEV V.B. Introduction to the sandpile model // Physica. - 1998. - A 254. - P. 97-116.
21. MASUDA N., PORTER M.A., LAMBIOTTE R. Random walks and diffusion on networks / Physics Reports. - 2017. -Vol. 716-717. - P. 1-58.
22. PRISNER E. Parallel Chip Firing on Digraphs // Complex Systems. - 1994. - No. 8. - P. 367-383.
23. SMYTH W.D., NASH J.D., MOUM J.N. Self-organized criticality in geophysical turbulence // Scientific Reports. - 6 March 2019. - URL: https://www.nature.com/articles/s41598-019-39869-w.
24. SPEER E.R. Asymmetric abelian sandpile models // Journal of Statistical Physics. - April 1993. - Vol. 71, Iss. 1-2. - P. 61-74.
5. Приложение
Доказательство утверждения 4.
1. Докажем бесконечное функционирование сети.
Пусть 91(0) > Г и 92(0) ^ г.
На такте Ь = 1 вершина г>1 отдает каждой другой вершине ресурс Я11(>}-г, а вершина ь2 - ресурс . По утверждению 2,
суммарный ресурс, полученный вершинами г>1 и г>2 от остальных вершин vi, г = 3,п, одинаков. Обозначим его через -шГе.^(0). Тогда к концу первого такта имеем:
91(1) = г + 92(0) - Г + ыге8^0), тъ — 1
«2(1)= Г + ^—Г + ^(0). К концу второго такта аналогичным образом получим:
/о\ , ^1(0) — Г , ЫгезЛ0)
91(2) = Г + ^+ + ^(1),
. 92(0) — г У)гез^0)
92(2) =г + 72^—1)2 + + ^(1).
Соответственно на такте Ь имеем:
,91(0) — Г У)гез^0) 'ШгеЛ^ — 2) ^ = Г + (П — 1)Ь +(П — 1 у-1 + ." + —^-+ — 1),
/,\ I У2(0) — ' ^ге«4(0) . . шГевг(ь — 2) . /, ч
^ = Г + —-^ + ^ 1 + ... + —.-;-+ —1),
92(0) — Г + (0)
(п — 1)1 (п — 1) если £ - четное, и
т . 92(0) — Г И]гез1(0) , ^геы^ . ¡,
^ = Г + {п — 1)ь +{п _1 у-1 + —;-+ — 1),
. 91(0) — Г >Ыге.4(0) „гка^ ^ =г+...+ г. — 1 + ^—1), если £ - нечетное. 32
п — 1
2)
п — 1
^е^ — 2)
п — 1
2)
Видим, что > г при любом четном £ и > т при любом нечетном 1 Значит, как минимум одна вершина будет перенасыщенная, а это значит, что сеть не остановится и будет функционировать бесконечно.
2. Докажем асимптотическую остановку сети. Рассмотрим динамику суммарного свободного ресурса. Пусть кг - число перенасыщенных вершин на такте ¿, кг+1 - число перенасыщенных вершин на такте Ь + 1. При этом кг < п — 1 (так как Ш ^ гп, то все вершины сети не могут быть перенасыщенными). Суммарный свободный ресурс на так-
кг
те £ равен AQ(t) = £ Дqí(t). Суммарный свободный ресурс на
г=1
такте Ь + 1 тогда будет выглядеть следующим образом:
кг
дд(£ + 1) = ^ Дф + 1)+ ^ Дф + 1) =
г=1 г=кг + 1
кг кг « ил кг «
= Е £ Д—? + £ (ЕД—? — с — *(«)>)•
г=1 .7 = 1, ^'=г г=й(+1 .7 = 1
Упростим первое слагаемое:
кг д /
£ £ = Л £ (£ ^м - ^м) =
г=1 .7=1, г=1 .7=1
кь кг
1 ^г 1 -I
Дщ о) — Е Д^)) = д®(1)-
3=1 1=1
-1 \ Ь / -1 I ) / , -1 Л
п — ' ) п — 1
Второе слагаемое не только упростим, но и оценим сверху:
£(£ - (г—*<>)) = Е (Д—? — (г - ,(*))) <
г=кг+1 3=1 г=кг+1
< £ ДМ = ~ — 1 — к, Д
^ п — 1 п — 1
г=fct+1
Подставляя в исходную формулу результаты преобразований, имеем:
Д^ + 1) < 1 Д(2(1) + п — 1 — к Д$(*),
п — 1 п — 1
п — 2
Дяа + 1) < ~2 дQ(t).
Таким образом, последовательность значений суммарного свободного ресурса ограничена сверху убывающей геометрической прогрессией, а, значит, предел этой последовательности равен нулю. Это, в свою очередь, означает, что стремится к нулю и общий поток в сети (исключая поток в петлях), так как он создается только свободным ресурсом. Очевидно, что изменение ресурса в вершинах при этом также стремится к нулю, т.е. сеть стабилизируется и останавливается. □
Доказательство утверждения 7.
Пусть в сети в момент времени ¿о имеется к перенасыщенных вершин VI, г = 1, к, к < п. Если вершины , ] = к + 1, п, являются насыщенными, то утверждение уже доказано. Рассмотрим ситуацию, когда среди г^, ] = к + 1,п, есть хотя бы одна ненасыщенная, и докажем, что эта вершина насытится за конечное число тактов.
В силу условия Ш > гп верно следующее: Д<3(1) ^ Ш — гп (равенство наступает тогда, когда все вершины сети насыщены).
От каждой вершины ы, г = 1,к, ненасыщенной вершине Vj на такте ¿0 + 1 передается г ресурса, если ^ г, и
ресурса, если < г. Объединив эти условия, по-
лучим, что V] получает тгп{г, } ресурса от каждой пере-
насыщенной вершины. Тогда
к
Яз(¿о + 1) — Ъ(¿о) = ^ тт^г, |.
г=1
Если хотя бы для одной вершины ы, г = 1,к, min{г, } = г, то имеем
Qj ( t о + 1) - Qj (t о ) = ... + Г > Г.
Если же для всех вершин vî, г = 1,к, min{г, An_iio)} =
Aq_ (to) = v ', то имеем
п— 1 '
Clj (tо + 1) - Qj (tо) = £ A=T = È А*(*о) =
î=1 î=1
= AQ(tо) > W - m n - 1 n - 1
Получаем, что на каждом такте ресурс в ненасыщенной вершине увеличивается не менее чем на min{г, }. Тогда любая ненасыщенная на такте ¿о вершина за число шагов не более чем
г - qî(h) min{г, }
достигнет своего насыщения. Отметим, что если
W—г п п— 1
min{г, W—n } = г, то
г - qî(h) _ г - ^(¿о) = Qî(tо)
тгп{г, } г
и, значит, вершина насытится за 1 такт. □
Доказательство утверждения 9.
Перейдем в координаты стандартной ресурсной сети: имеем полную однородную сеть без петель, ресурс в каждой вершине уменьшен на г. Тогда ограничение на суммарный ресурс принимает вид Ш < гп2 — гп = гп(п — 1), и в сети всегда будет хотя бы одна вершина из зоны 2_ (¿). Пусть в начальном состоянии имеется к вершин находятся в зоне 2 +(0) и п—к вершин - в зоне 2 _(0):
VI,...,Ук е2+(0), ук+1,...,ип (0).
То есть ресурс распределен следующим образом:
| > гс[и1 = г(п — 1), г = 1, к, г г°иг = г(п — 1), г = к + 1,п.
1. Согласно свойству 2 полной однородной сети без петель, вершина, попавшая в зону 2-(£), уже из нее не выйдет.
2. Рассмотрим вершину, которая в начальном состоянии находится в зоне 2 + (0). Не нарушая общности, положим, что это вершина и\. Выпишем её выходной и входной поток:
/Г4(0) = г1иЬ = г(п — 1),
/Г(1) = ££ Пв^ — 1) + Л £ *(0).
г=2 г=к+1 г=к+1
Разность между выходным и входным потоком есть некоторая величина А/, показывающая, насколько уменьшился ресурс в вершине:
1 п
а / = Гки\(0) — /к+1(1)=Кп — 1) — г(к — 1) — — £ 9г(0) =
п - 1
г=к+1
1п
= г(п — к)--- V %(0).
п 1
г=к+1
Используя тот факт, что <?г(0) > г(п — 1), г = 1, к, получим:
пк
£ Яг(0) = Ш — £^(0) <Ш — гк(п — 1).
г=к+1 г=1
Тогда можно оценить А / снизу:
1 Ш А f > г(п — к)---(Ш — гк(п — 1)) = гп--- > 0.
Таким образом, на каждом такте ресурс вершины 1 уменьшается на некоторую ограниченную снизу величину А / > гп — . То-
гда любая вершина Уг, находящаяся в начальном состоянии в зоне 2+(0), за число шагов менее чем
^(0) - Г(П - 1)
ГП —
п— 1
перейдет в зону 2 (¿). □
Доказательство утверждения 10.
При ш < гп2 с определенного момента времени все вершины сети передают ресурс по правилу 2, тогда матрица потока выглядит следующим образом:
Р =
Я2 (*) — п— 1
Я! (*)-п— 1
. Яп (*)- Яп (*)-\ п—1 п—1
яФ)-г\
п-1 Я2^)-г п— 1
То есть полную пропускную способность г вершина Ьг отдает в петлю, а оставшийся ресурс Адг{Ъ) = (£ ) — г делит пропорционально между остальными вершинами (работает по правилу 2 стандартной ресурсной сети). Таким образом, вершины отдают весь свой ресурс, и каждое новое состояние состоит из ресурса, пришедшего на данном такте: <(£ + 1) = т(£ + 1). Тогда для вершины V] новое состояние qj(£ + 1) будет + 1), т.е. суммой элементов ^'-го столбца матрицы Р(¿).
^(¿ + 1) = г+ £
г=1, г=j
ЧгО — Г П 1
= Г +
£
=1, =
д^ ) г(п — 1) 1
П 1
п — 1
у
п — 1
=1, =
гп(г + 1) = (/Г(* + 1), ..., /пп(* + 1)) =
п— 1
П — 1
г=2
г=1, г=2
П 1
..., Яп(^)
/о 1
I о
II
Получили искомую формулу (3):
г=1
1 \
1 .о..
= Я(ь)К'
+1) = Ж^п', К =
о1 1о
П 1
1 1
. □
V1 1 ... о)
К доказательству этого утверждения можно подойти по-другому. Перейдем в координаты соответствующей стандартной ресурсной сети.
Для нахождения вектора новых состояний вершин в координатах стандартной ресурсной сети + 1) можно применить имеющуюся из теории стандартных ресурсных сетей формулу: + 1) = Здесь матрица К' будет иметь вид:
К =
0
1
П 1
1 о
1 1
V1 1 ... о)
Переходя обратно в координаты исходной модели, имеем:
ж$ + 1) = ж8& + 1) + ьг,
где Ьг = ( г, ..., г), получим:
+1) = (Я(г) - ьг)К + ьг = Ж^к' - ьгК + ьг = Ж^к' - ь.
где
Ь = ЪГВ! - Ъг = Ъг(К - Е) =
1
1
1
п 1
(г, ..., г)
1 — п 1 1 1 — П
1
1
1 1
1— п
= (0,..., 0).
Тогда зависимость вектора нового состояния сети от вектора предыдущего состояния сети выражается формулой (3). □
Доказательство теоремы 3.
Полная однородная ресурсная сеть с «жадными» вершинами при Ш < гп2 стабилизируется в случае, если стабилизируется соответствующая стандартная ресурсная сеть. В [9] доказано, что регулярная сеть при малом ресурсе всегда имеет единственное предельное состояние. Из этого факта и утверждения 10 следует стабилизация исходной сети с «жадными» вершинами при Ш < гп2.
Найдем соответствующее предельное состояние сети. Используя формулу из утверждения 11, запишем состояния сети на тактах £ + 1, £ + 2, £ + 3, ... £ + к.
<(*+1) = <(1)К' = —Ц-(д^), ..., яп(г))
п 1
01 10
1 1
п 1
11
.. 0
п 1
г=2
=1
1
п — 1 1
(V — ..., Ш — 9п(*))
п 1
Ш, ..., Ш -
п 1
< ( ),
<(* + 2) = <(* + 1)Д' =
1
1
( П - 1)2
о 1. .1
1 о. . 1
1 1. .. .о..
1 п п— 1
(((П- 1)Ш, ..., (п- 1)ш)-(£<&(*), ..., £дг(*))
(П - ^ ¿=2 ,= 1
= (П-1- (ПП^Ж- ш) +
Аналогично находим + 3):
Я(* + 3) = Я(* + 2)Я' = ... = / 1 1 1 „Л 1
\п - 1 (п - 1)2 (п - 1)3А ' •••' ) (п - 1)3 Для такта Ь + к таким образом получим:
г=1
к
Перейдем к пределу при к ^ то. Учитывая, что для первого слагаемого
Иш —--г = 0,
к^ж (1 - П)к
а предел суммы убывающей геометрической прогрессии равен
1
к 1 (£ (1 -ПУ )
1 А 1 — П 1
Иш . / . , = = ,
к^-ж \ (1 -П)г/ 1 1 П
=1 - 1 - П
получим искомое предельное состояние системы:
я* = (-,..., -У □
П П
1
Доказательство теоремы 4.
Перейдем в координаты стандартной ресурсной сети: имеем полную однородную сеть без петель, ресурс в каждой вершине уменьшен на г. Тогда суммарный ресурс равен Ш = гп2 — гп = = гп(п — 1). Пусть в начальном состоянии ¿ = 0 имеется к(0) =ко вершин из зоны 2+(0) и п — ко вершин из зоны 2-(0):
VI,..., уко е 2+(0), уко+1,..., уп е 2-(0). То есть ресурс распределен следующим образом:
\ > г%* = г(п — 1), г = 1, к0,
Чг(0) = \ ^ ы
^ гоиЪ = г(п — 1), I = к0 + 1, П.
Целочисленная функция к(Ь) является монотонной невозрастаю-щей: вершины из зоны 2-(0) не могут ее покинуть, согласно свойству 2 полной однородной сети без петель, а, значит, и к(£) с увеличением времени не возрастает.
Рассмотрим входной и выходной потоки вершины Уг, которая в начальном состоянии находится в зоне 2 +(0), т.е. г е {1, ..., ко}:
(0) = г°иЬ = г(п — 1),
ко п п
/Г(1) = Е- + Е П—1 = г(ко — 1) + ^—г Е *(0).
г=2 г=к0+1 г=к0+1
Состояние вершины ьг на следующем такте Ь = 1: %(!)= 9г(0) — ¡0и (0) + /Г(1) =
1 п
= &(0) — г( П — 1) + Т (ко — 1) + -- V Яг(0) =
п — 1
= ко+1
1 п
= 9г(0)+г(ко — П) + —у ^ ^(0).
г=ко+1
На такте £ = 2 получим аналогичную формулу; вместо ко здесь
будет уже к(1) = к1 - количество вершин из зоны 2 +(1) на такте
г = 1:
1 п
®(2) = ®(1)+г(ь п) + п~1 £ ^(1)-
г=к! + 1
Для такта + 1 имеем соответственно
1 п
&(* + 1) = + к^) -п) + —- V ф).
П 1
г=к(*)+1
Рассмотрим динамику избыточного ресурса (обозначим его 5Яг(£)) в вершине Уг. Под избыточным ресурсом будем понимать ту часть ресурса вершины гц, которая позволяет ей находиться в зоне 2 + (£) на такте т.е. 6Яг^) = Яг^) - г(п - 1). Распишем избыточный ресурс в вершине Vi на такте Ь + 1:
5 Яг (¿ + 1) = Яг (¿ + 1) - г(п - 1) =
1 п
= % (г ) + г( щ) -п) +-- V - г(п - 1).
П 1
г=к(*)+1
Выразим 5+ 1) через 6%(£):
1 п
6я^ + 1) = 6я^) + г(Щ) - п) + —— £
П 1 ,
г=к(*)+1
Последнее слагаемое в формуле распишем:
. п 1 к(4)
—П Е *<«> = —-Е*<«>) =
г=к( ¿)+1 «=1
( тп(п - 1) = ГП - —&(*).
— - 1 V ' / п - 1
=1 =1
Учитывая, что = г(п — 1) + 5дг^), г = 1, ..., к(1), имеем
п , Щ)
^ д№ = гп--гт Е (г(п — 1) + (^))
п — 1 ^ п — 1
г=к(г)+1 1=1
т
= гп — гк(ь)--- У $дг(1).
п — 1 —^
г=1
Тогда
к( )
бдг(* + 1) = 6дг(г) + г(к(г) — п) + гп — гк(Ъ)--Ц- V 5д^) =
П — 1 ' ^
г=\
1 т 1 2
= <ЗД)--- V 5ф) < 5дг(I)---5ф) =-25ф),
п — 1 — п — 1 п — 1
г=1
П — 2
5ф + 1) <-25ф).
п — 1
Последовательность избыточных ресурсов вершины г>г ограничена сверху убывающей геометрической прогрессией, а значит, предел этой последовательности равен нулю. Отсюда следует, что ресурсы в вершинах, принадлежащих зоне 2 +(0) в начальный момент времени, стремятся к состояниям д*(Ь) = г(п — 1), г = 1,..., к0. По закону сохранения ресурса в сети, ресурсы в вершинах, принадлежащих зоне 2-(0) в начальный момент времени, также будут стремиться к состояниям д*(Ь) = г(п — 1), г = к0 + 1, ..., п. Таким образом, сеть стабилизируется с предельным состоянием <3*(1) = (г(п — 1), ..., г(п — 1)), или, в переводе обратно в координаты «жадных» вершин, с предельным состоянием <3*(1) = (гп, ..., гп). □
1
RESEARCH OF COMPLETE HOMOGENEOUS «GREEDY-VERTICES» RESOURCE NETWORKS
Ludmila Zhilyakova, V.A. Trapeznikov Institute of Control
Sciences of RAS, Moscow , Moscow,
Doct.Sc.(zhilyakova.ludmila@gmail.com).
Nadezda Chaplinskaya, V.A. Trapeznikov Institute of Control
Sciences of RAS, Moscow , Moscow, student
(nadya1462@gmail.com).
Abstract: The modification of the graph dynamic model "resource network" -"resource network with greedy vertices" is described. In this model, the graph vertices at each discrete-time moment exchange resources through the edges with limited throughput, first passing the available resource to themselves via the loop and then distributing the remaining resource to adjacent vertices according to the "standard" resource network rules. These are two rules with threshold switching: if the vertex resource exceeds the total throughput of all its outgoing edges, it gives away the full throughput to each edge; otherwise, it gives the entire available resource, distributing it in proportion to the throughputs of the outgoing edges. The process of functioning of a complete homogeneous resource network with "greedy" vertices at different values of total resource and different initial states is investigated. Possible network states are described; a non-standard state - the shutdown of the network has been identified. Two total resource thresholds, separating zones of different network behavior, have been found: the first threshold divides the zones of "insufficient" and "sufficient" resources, the second divides the zones of "sufficient small" and "sufficient large" resources. For each zone the functioning of the network is described and the asymptotic states and flows are investigated. The examples demonstrating numerical experiments are given for all typical situations.
Keywords: resource net, «greedy-vertices» model, graph dynamic threshold model.
УДК 519.1 ББК 22.176
DOI: 10.25728/ubs.2021.89.1
Статья представлена к публикации членом редакционной коллегии Г.А. Угольницким.
Поступила в редакцию 17.11.2020. Дата опубликования 31.01.2021.