Управление предельными состояниями в поглощающих ресурсных сетях Текст научной статьи по специальности «Математика»

УДК 519.857

УПРАВЛЕНИЕ ПРЕДЕЛЬНЫМИ СОСТОЯНИЯМИ П ПОГЛОЩАЮЩИХ РЕСУРСНЫХ СЕТЯХ1

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

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

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

ВВЕДЕНИЕ

Ресурсная сеть — динамическая параллельная пороговая модель, предложенная в работе [1] и получившая развитие в работах [2—5]. Она представляет собой ориентированный взвешенный граф, в котором вершины обмениваются ресурсами в дискретном времени. В отличие от классической потоковой модели Форда—Фалкерсона [6] и ее динамических модификаций [7], в ресурсной сети не существует направления движения ресурса от источников к стокам; кроме того, распространение ресурса происходит по всем путям, а не только по некоторым выделенным и удовлетворяющим заданным условиям. В этом смысле модель ближе к процессам рассеяния на графах [8, 9], названных «интегральными по путям» [9]. В работе [9] дан большой обзор этих моделей. Еще один класс моделей, с которыми ресурсная сеть имеет сходство, — это пороговые модели «игра выстреливания фишек» [10, 11] и «абелева куча песка» [12, 13], находящие применение в различных областях, в том числе, при моделировании «самоорганизующейся критичности» [14—16].

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

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

Ресурс является бесконечно делимым, благодаря этой особенности разнообразные переходные процессы сходятся к предельным состояниям, которые в целочисленных пороговых моделях существуют не всегда. Кроме локальных пороговых значений для вершин, в сети имеется генеральная характеристика: пороговое значение Т суммарного ресурса, такое, что при ресурсе, не превосходящем это значение, все вершины сети за конечное число шагов оказываются функционирующими по одному и тому же правилу. В этом случае переходные процессы в сети описываются однородной цепью Маркова. Если суммарный ресурс сети больше Т, в ней появится хотя бы одна вершина, суммарный ресурс которой превзойдет ее выходную пропускную способность. В этом случае вершина переходит на другое правило функционирования, отличное от того, по которому отдают ресурс остальные вершины. Цепь Маркова, описывающая такую сеть, уже не является однородной. В эйлеровых сетях [5], в которых каждая вершина имеет одинаковые суммарные входную и выходную способности, значение Т совпадает с суммой пропускных способностей всех ребер сети. В регулярных несимметричных сетях [3, 4] значение Т строго меньше этой суммы.

Все сети, исследованные в работах [1—5], имеют регулярные стохастические матрицы (термин «регулярные» соответствует здесь определению, данному в работе [17], и отличается от его использования в работах [18, 19]), начиная с некоторого натурального к, все степени таких матриц строго положительны. Сети, описанные в настоящей работе, в этом смысле регулярными не являются. Поглощающие сети — это сети со стоковыми вер-

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

1. ОСНОВНЫЕ ОПРЕДЕЛЕНИЯ

Ресурсная сеть представляет собой ориентированный взвешенный граф О = (V, Е), \ У| = п, веса

ребер которого г.. обозначают их пропускные спои

собности. Матрица пропускных способностей — неотрицательная матрица Я = (гЦ)пХп.

Вершинам приписаны неотрицательные числа называемые ресурсами и изменяющиеся в дискретном времени t.

Состоянием 0(1) сети в момент t называется вектор = (д^), ..., дп(ф.

Состояние 0 * = (д*1, ..., д*п) будем называть предельным, если для любого е > 0 существует такое, что для всех t > \ д* — д.(1)\ < е, I = 1, 2, ..., п.

Суммарной пропускной способностью сети называется сумма пропускных способностей всех ее ребер:

шина V. в ребро, соединяющее ее с вершиной V

к

г,,, =

I I гц.

I = 1 ] = 1

Суммарную пропускную способность входных ребер вершины будем называть ее входной пропус-

п

кной способностью и обозначать через г¡п = I г.;

И = 1

суммарную пропускную способность выходных ребер, соответственно, назовем выходной пропуск-

п

ной способностью и обозначим через г= I г..

И = 1

Пропускная способность петли (если она существует) входит в обе суммы. Эти величины являются обобщением полустепеней захода и исхода вершин в не взвешенном ориентированном графе.

Суммарный ресурс, находящийся во всех вершинах, обозначим через W. В сети выполняется закон сохранения: при ее функционировании ресурс остается постоянным — не поступает извне и не расходуется.

Распределение ресурса в сети происходит по одному из двух правил, выбор которых зависит от значения ресурса в вершинах. В момент t + 1 вер-

отдаст:

г., единиц ресурса, если #.(*) > г"ш (правило 1);

^ д(Р), если д^) 1 г{ (правило 2).

Стохастическая матрица, соответствующая матрице пропускной способности Я, задается формулой:

г11 г1п

Я' =

вШ г1

!п1

оШ

ои1 г1

оЫ

Если все вершины сети функционируют по правилу 2, для вектора состояния выполняется:

ш + 1) = от'.

Множество вершин, для которых д(1) < г"ш, называется зоной Z~(1). Вершины из Z~(1) функционируют по правилу 2. Z— множество вершин, ресурс которых больше их выходной пропускной способности, они функционируют по правилу 1.

Т — пороговое значение ресурса, такое, что при № < Т все вершины, начиная с некоторого переходят в зону Z (/); при № > Тзона Zне пуста, начиная с некоторого Г.

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

2. ПОГЛОЩАЮЩИЕ РЕСУРСНЫЕ СЕТИ

2.1. Свойства поглощающих сетей

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

Я =

Б 01

Я1 Я2

(1)

где Б — диагональная матрица размера I х I с произвольными неотрицательными диагональными элементами, равными пропускным способностям петель в стоках, О1 — нулевая матрица размера I х (п — I), Я1 — матрица размера (п — I) х I и Я2 — квадратная матрица (п — I) х (п — I). Матрица Я1 состоит из пропускных способностей ре-

п

п

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

Стохастическая матрица R', соответствующая такой матрице R, будет иметь вид

R' =

Ei Oi

Ri R2

(2)

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

Перечислим свойства поглощающих сетей, вытекающие непосредственно из их топологии.

Свойство 1. Предельное состояние в поглощающих сетях существует. Координаты вектора предельного состояния, соответствующие переходной компоненте, равны нулю. Предельный поток состоит из потока в петлях стоковых вершин. ♦

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

Свойство 2. Наличие или отсутствие петель в стоках не влияет на количество ресурса в стоках в предельном состоянии. ♦

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

Свойство 4. Если в сети имеется один сток, он в предельном состоянии соберет весь суммарный ресурс W.

Свойство 5. В поглощающей сети с несколькими стоками и произвольными пропускными способностями ребер в общем случае предельное состояние зависит от начального при любом суммарном ресурсе. ♦

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

2.2. Поглощающие сети с единственным предельным состоянием

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

R =

D Oi

H R2 '

(3)

Ее отличие от матрицы (1) заключается в блоке H — он имеет ранг, равный единице: rankH = 1. Это означает, что пропускные способности ребер, ведущих из переходной компоненты в стоки, пропорциональны.

Теорема 1. В поглощающей сети с матрицей (3), в которой rank H = 1, для любого начального состояния Q(0) = (qx(0), ..., q1(0), ql + Д0), ..., qn(0)) предельное состояние определяется формулой:

(

Q * =

h

qi(0) + И- W , ..., qi(0) +

hsum

1 4 h,

+ -h- W-, 0, ..., 0 h

1 * о/ж

(4)

п п -1 п - II

гОе ЦТ = I qг'(0), У =1 V = 11 V

I = I + 1 к = 1 к = 1 у = 1

Доказательство. 1. Вклады в предельный ресурс #1(0), ..., #¿(0) очевидны, так как ресурс из стока выйти не может.

2. Если ранг матрицы равен 1, то все ее строки и столбцы пропорциональны. Матрицу Н можно записать в виде:

H =

a2 h1 a2 h2

hn -1 a2hn -l - alhn - l

afti alh2

\

Если некоторая вершина V. из переходной компоненты функционирует по правилу 1, она будет отдавать в стоки ресурс, равный пропускным способностям ребер: к, а2к., ..., а,к.; т. е. в каждый сток он попадет в пропорции 1 : а2 : ... : а,. Если вершина функциони-

рует по правилу 2, она отдает в стоки ресурс, равный

( h¡ . . ^а, ^

I -±, q<t)' -Ьа q<t)' •••' -Ь q/t) I, и стоки получат ресурс

^ ^ ^ )

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

1

W

а2

W-.....

1 + а.2 + ... + а, '1 + а.2 + ... + а, ' '

а,

1 + а2 + ... + а

W-.

;

Входные пропускные способности стоков равны соответственно:

h, = h + h + ... + h ,, 1 1 2 n — P

К = + h2 + ... + hn _ ,), k = 2, ..., l.

Суммарная пропускная способность матрицы Н: h = (h, + К +...+ h ,)(1 + а + ... + а,),

sum v 1 2 n — //v 2

h

i _

1

h

k —

1 + а2 + ... + а,

1 + а2 + ... + а,

,

hs

1 + а2 + ... + а,'

откуда и вытекает формула (4). ♦

Следствие 1. Если в начальном состоянии весь ресурс находится в переходной компоненте, то предельное состояние единственно и не зависит от начального распределения ресурса по вершинам.

Следствие 2. Из формулы предельного состояния видно, что матрица Я2 в сети с гапкН = 1 может быть любой. Перетоки ресурса внутри переходной компоненты никак не влияют на предельное состояние.

2.3. Поглощающие сети общего вида. Пороговое значение Т

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

Я2 строго положительна.

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

Таким образом, в поглощающей сети с rankR > 1 предельное состояние всегда зависит от начального.

Из сказанного вытекает еще одно основное отличие поглощающих ресурсных сетей от регулярных и эргодических — они не имеют порогового значения Т.

Теорема 2. В поглощающей сети порогового значения Тресурса не существует.

Доказательство. По определению пороговое значение Т — это такое количество ресурса, что при W > T для некоторого момента времени t' и некоторой вершины v. выполняется: V t > t' v; е Z+(t). В поглощающих сетях такими вершинами могут быть только стоки. Но принадлежность вершины зоне Z+(t) означает, что ее ресурс

Г out

превосходит ее выходную пропускную способность ri .

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

2.3. Предельные состояния в поглощающих сетях

Найдем предельную матрицу произвольной пог-

ч г» f Ю ч

лощающей сети R и опишем ее свойства.

Лемма 1. Пусть поглощающая сеть с l стоками

представлена матрицей (1). Тогда R'ю — предел степеней стохастической матрицы R' — определяется по формуле:

R'ю =

( E2 - R2) 1R1

Oi

O

(5)

2

где Е2 — единичная матрица размера (п — I) х (п — I).

Доказательство. Степени матрицы R' имеют вид:

Rk =

E1 Oi

R, + R2R, + (R2 )2R1 + . . + (R2 )k - 1R1 (R2)k

к

Для цепи Маркова (R2) ^ O2 при k ^ да, где O2 —

да

к

нулевая матрица размера (п — l) х (п — l); ^ (R2) =

к = 0

= (^2 — R2) , причем, матрица — R2) существует [17, теоремы 3.1.1 и 3.2.1)].

Тогда R; + R2 R^ + (R2 )2 R; + ... + (R2 )к Rí + ... = = — R2 )-1 R^, и матрица R'ш имеет вид:

R, Ю =

( e2 - R2) 1r1

o,

O

2 У

а

k

h

h

sum

sum

а

Теорема 3. В поглощающей сети с I стоками матрица R'" остается неизменной при любых изменениях диагональных элементов.

Доказательство. Матрице Л, заданной формулой (3), соответствует стохастическая матрица

R ' =

E 0;

R1 R2

1. Для изменения диагональных элементов матрицы Б доказательство очевидно.

2. Докажем это утверждение для диагональных элементов матрицы Л2.

Пусть В2пек = - Б(Дг,..), где Б(Дг,..) — диагональная матрица размера (и — /) х (и — /), элементы которой характеризуют изменение пропускных способностей петель. Будем предполагать, что 0 < Дг. < и хотя бы для одной вершины выполнится строгое неравенство Дгкк > 0. Если для некоторой вершины Ук выполняется равенство Дгкк = гкк, петля этой вершины исчезает вообще.

Тогда новая матрица пропускных способностей сети

R =

new

D °i| - 03 °1 | = D Oi

R; R2 у V 04 D(Arii) у v Ri R2 new

Стохастическая матрица, соответствующая матрице Лп™, будет отлична от Л.

Зависимость между матрицами Л1 и Л2 и матрицами

Л1пеы и Л2пек :

Выразим (E2 - R2„êw ) 1R [К

lnew через исходные матрицы:

(Е — )-1 Л 1пем = (Е — (Б (— Бд))-1Б Л1 .(7) Все диагональные элементы матрицы Б ' отличны от 0. Это означает, что матрица Б ' 1 существует. Тогда в первом множителе правой части равенства (7) Б ' можно внести в скобку:

(Б ■( Л2 — Бд))- 1Б =

Е2 ^д/

-Д// "1 V- 'Л — Л2)

(E2

= (D ■ -1£2 - D ■ -1D ( R2 - Da))- 1 = = (D '-1 - (R2 - Da))- 1R1 = (D '-1 + Dh - R2 )-1 R; .

1

+ da-

Вычислим сумму D ' Это две диагональные матрицы. Элементы матрицы

Т\> out и out . ч

D равны г, /( ri - Ar;;), тогда соответствующие эле-

out

—1 r- — Ar--менты D ' равны —-- . Элементы матрицы Da рав-

out

ны Ar;;/ r,- .

Для каждого i выполнится:

out out out

( ri - Ar,,)/ri + Ar,/ ri = 1

Отсюда

(d ■ 1 + Da - R2 )) 1 R; = (E2 - R2 ) 1 R; ,

а это означает, что

(E2 - R2new

и соответственно

)—1R

= (E2 - R2 ) 1R1 '

lnew

Л 1™ = Б( гГ /(С — Дг,,)) Л1,

^ = Б( /Г /(г°ш — Дг,.,.))( Л2 — Б(Дг,,/ г,0")), (6)

ДОШ и ОШ .

Г1 /(Г1 — Дг,,)) — диагональная матрица размера (и — /) х (и — /) с элементами (г°ш/(— Дг,,)). Именно в такой пропорции пересчитываются элементы двух нижних блоков матрицы Л '. Обозначим ее через Б ';

Б(Дгй/ ггОмГ) — диагональная матрица, соответствующая изменению матрицы Л2 после изменения пропускных способностей петель. Обозначим ее через БД.

Тогда равенства (6) можно переписать в более компактном виде:

R l new D R1 ,

R.

2new

= D ( R2 - Da).

Новая матрица R 1new =

E Oi

R R 1 new R'2new

по-прежне-

му остается стохастической, и поэтому предел ее степеней существует и равен

R =

R new =

(E2 - R2new) lR 1

O1

O

2 У

R "» = Rnew =

( E2 - R2 ) 1r 1

O,

O

= R,ю.

2 У

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

Теорема 4. Пусть поглощающая сеть с I стоками представлена матрицей (1). Тогда для любого ресурса W и любого начального состояния Q(0) = = ^(0), ..., qn(0)) предельное состояние рассчитывается по формуле:

Q * = Q(0)R'

(8)

где Я'" — предельная матрица, определяемая по формуле (5).

Доказательство 1. Если начальное состояние таково, что все вершины из переходной компоненты функционируют по правилу 2, доказательство теоремы очевидно. Для каждого I выполнится: 0(?) = 0(0)Л''. Совершив предельный переход, получим искомое равенство (8).

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

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

г,,(0 =

/ ,ч , вит г№ если t) < ^ ,

t) - X гф если t) > ГГ

] * I

Нормируя матрицу R(t), получим стохастическую матрицу R'(t), которая задает неоднородную цепь Маркова. Для этих матриц выполняются условия:

• Q(t + 1) = Q(t)R'(t);

• вектор Q(t) совпадает с вектором состояния сети, заданной матрицей R, с начальным состоянием 0(0) = (91(0), qд(0))•

Тогда функционирование сети описывается правилом:

+ 1) = Q(t)R (0 = Q(0)(R '(0)R '(1)•••R (/))• За конечное число шагов т все вершины из переходной компоненты перейдут на правило 2. Тогда, начиная с шага т + 1, сеть начинает описываться однородной цепью Маркова: R(t) = R и R '(t) = R •

Вектор 0(т + к) можно записать в виде

0(т + к) = 9(0) I П R '(t)| R '

1 т = о

При к ^ да получим:

т = 0

0* = 0(0)I П R (о) 11ш R 'К = 0(0)I П R '(/)1 R

т = о

По теореме 3 матрица Rне изменяется при любых изменениях диагональных элементов.

Тогда для всех t = 0, ..., т выполнится R '(t)R= R

Отсюда 0 * = 0(0^♦

Следствие. Вектор предельного состояния О * не зависит от наличия или отсутствия петель в вершинах.

Теорема 5. Пусть поглощающая сеть с I стоками представлена матрицей (1). Тогда элементы 1-й строки матрицы Я'" (/ > I) соответственно равны компонентам вектора предельного состояния при начальном состоянии Ог(0) = (0, ..., 0, 1, 0, ..., 0), где единице равна 1-я компонента:

(

Я'" =

л

О

1 +1

О*

^п

(е1, ..., е1 — первые I вектор-столбцов единичной матрицы (Е)пХп).

Доказательство вытекает непосредственно из формулы (8).

Ранее было сформулировано свойство поглощающих сетей: предельное состояние в них не зависит от наличия/отсутствия петель и, соответственно, диагональных элементов в матрице Я. Это

7"» ЮТ ч

означает, что матрица Я' не зависит от значений диагональных элементов г.. (I > I) матрицы Я, а

следовательно выражение (Е2 — Я'2) 1Я1 остается неизменным, хотя при изменении диагональных элементов Я все элементы матриц Я1 и Я2 также изменяются.

Замечание. Предельная матрица Я'™ является собственным проектором лапласовской матрицы (Е — Я') [20, 21]. Также этот проектор совпадает с нормированной матрицей максимальных исходящих лесов орграфа, полученного из исходного орграфа сети изменением направлений ребер [22]. Для вычисления собственного проектора существуют конечные алгоритмы. В работе [21] доказана теорема о том, что собственный проектор любой квадратной матрицы А может быть вычислен с помощью аннулирующего многочлена для матрицы

Лк, где к 1 V (V = тёЛ — размер наибольшего блока нулевого собственного значения в жордановой форме для Л). В нашем случае Л = (Е — Я') — лап-ласовская матрица, для которой всегда V = 1. Таким образом, предельная матрица Я'™ может быть найдена с помощью любого аннулирующего многочлена для лапласовской матрицы (Е — Я').

3. УПРАВЛЕНИЕ В ПОГЛОЩАЮЩИХ СЕТЯХ

Предел степеней стохастической матрицы Я' поглощающей сети с I стоками (формула (2)) имеет вид:

Я'т =

( Е2 - ^ ) ^ 1

О,

О

2 I

Доказано (см. теорему 3), что предельное состояние зависит от начального при любом количестве ресурса в сети, если ранг матрицы (Е2 — Я2 )-1 Я1 больше единицы. Помещая ресурс в разных пропорциях в вершины из переходной компоненты, можно получить разное распределение его по стокам.

Будем рассматривать задачу управления в поглощающих сетях, в которой переходные вершины являются управляющими, стоки — целевыми вершинами. Управление заключается в распределении ресурса по вершинам переходной компоненты. Управляющих воздействий в процессе функционирования сети не происходит: закон сохранения ре-

сурса в сети не нарушается. Рассмотрим две различные задачи управления, которые назовем прямой и обратной задачей.

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

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

3.1. Прямая задача управления. Распределение фиксированного суммарного ресурса между стоками

Такая постановка задачи управления не всегда

71 г 00

имеет точное решение, однако, зная матрицу к и используя свойство линейности оператора перехода в предельное состояние, можно задать начальное состояние, при котором решение будет иметь минимальное расстояние до целевого в евклидовом пространстве. Если целевое предельное состояние О = (#1, ..., #1, 0, ..., 0), то задача управления формулируется в виде задачи оптимизации:

Тогда справедливо равенство:

Q - G'll2 = (qg

qg )2 + ... + (qg

q* )2 ^ min.

Матрица (Е2 — К2) 1 Я[, обладает тем свойством, что сумма элементов любой ее строки равна единице. Таким образом, ее можно назвать квазистохастической («квази» — поскольку в общем случае она не является квадратной).

Введем обозначения:

-1

(E - R2) 1R = R

yass '

подвектор начального состояния с компонентами, соответствующими переходным вершинам, обозначим через Qpass(0), его длина равна (n — l);

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

Целевой вектор Qg тоже редуцируем до l ненулевых компонент.

По условию W = const, следовательно, векторы Qpass(0), Q* и Q можно записать в виде:

Qpass(O) = ..., _ W X

x, l 0, i = 1, ..., n — l,

Q* = (P1W, ..., P7W), p, > 0, i = 1, ..., l,

Qg = (YiW ..., Yi W), Yi > 0, i = 1, ..., l.

Q = Qpass(0) R

lpass

или (PiW, ..., P7W) = (XiW, ..., Xn- WR

'pass

отсюда

l

(Pl, ..., р/) = (хр ..., хп - ); р, = £ х(^ )г

) = 1

Задача состоит в поиске такого вектора

II II2

(х1, ..., хп - 7), что значение ||01 - 0*||2 минимально. Отсюда для фиксированного ресурса

Yi, ..., Yi) - (xi .», xn-1)RP

pass 2

^ min,

n - l

I X, = 1

i = 1

I x, > 0.

Здесь у. > 0, I = 1, ..., I; £ у. = 1, у1 : у2 : ... : у1 —

I = 1

целевая пропорция распределения ресурса W в стоковых компонентах; х, : х.,: ... : х , — искомая

' 1 2 п — 1

пропорция распределения ресурса W в переходных

вершинах, RZss = (Ег — R2 )-1 R1.

Поскольку целевая функция квадратична, а ограничения линейны, эта задача относится к задачам квадратичного программирования [23, 24]. Ее

нужно решать, только если матрица Rpass имеет

ранг, больший единицы. Если rank Rp^ss = 1, предельное состояние не зависит от начального, и управлять ресурсом нельзя.

Функция fx) = || (Y1,..., Yi) - (Xi, ..., Xn - i)R'p2s 112 -выпукла.

Доказано, что для любых векторов x и y выполняется неравенство Йенсена:

f (x + У) m f( x ) + f( У).

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

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

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

цхр ..., %п - р X) =

I Г п-1 \2 Гп-1 \

= I у/- I г;?х\ + 1 ^х(-1|. / = 14 ]=1 ] ч = 1 ]

Необходимые условия локального экстремума:

dL

l ( П - l

д- = -2 X rkr Yi- I xj I + a. = 0

k i = 1 v j = 1 k = 1,..., n -1;

dL

da

(9)

n-l

= I x-1 = 0;

i = 1

хк 1

Это система п — I + 1 линейных уравнений с п — I + 1 неизвестными.

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

В общем случае существуют разнообразные численные методы для решения задачи наименьших квадратов с ограничениями-равенствами и неравенствами. Алгоритмы решения задачи наименьших квадратов с линейными ограничениями приведены в работе [25].

3.2. Обратная задача управления. Достижение заданного значения ресурса

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

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

г да г да

полниться условие г н = тах тк1 .

' к

Пусть д* — количество ресурса, которое должен получить сток с номером I в предельном состоянии. Тогда общее количество ресурса, которое в начальном состоянии должно находиться в переходной вершине V, будет равно д* /'да. Пос-

кольку весь ресурс в начальном состоянии должен быть сосредоточен в вершине v., выполняется равенство W = q* / j.

Это значение суммарного ресурса минимально для заданного значения q*.

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

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

Если q* , ..., q* — ресурсы, которые должны

'1 1k

получить стоки с номерами i1, ..., ik (i1, ..., ik m l), то суммарный необходимый ресурс обозначим через Wmin: Wmin = q* + ...+ q*k. Вектор предельного состояния Qg тогда можно записать как

Qg = (Y1 Wmin, ..., YlWmin), где Y' = q* /Wmin. Если для некоторого стока с номером i условие на ресурс отсутствует, считаем, что qi* = 0, и соответствующая

компонента Qg также равна нулю.

Поиск начального состояния для получения наилучшего приближения вектора Qg соответствует прямой задаче из п. 3.1. Решив ее и найдя вектор начального состояния Qpass(0), для которого вектор предельного состояния Q* будет иметь минимальное расстояние от вектора Qg, оценим разность AQ = Qg — Q*.

Нас будут интересовать только те компоненты вектора AQ, которые соответствуют индексам i1, ..., ik. Все они неотрицательны. Пусть Aqm =

= maxAq. > 0, j e {/,, ..., ik}. Тогда новое значение

j j 1 k

суммарного ресурса зададим как W = Wmin + kAqm, где k — число стоков, в которых ресурс в предельном состоянии должен быть отличен от нуля. Значение W гарантирует, что одна компонента нового вектора предельного состояния будет равна соответствующей компоненте целевого вектора Qg = (Y1Wmin, ..., YlWmin), а остальные будут не меньше. По построению значение W минимально.

ЗАКЛЮЧЕНИЕ

В работе рассмотрены ресурсные сети, содержащие стоковые вершины, т. е. вершины без исходящих ребер. Класс таких сетей назван погло-

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

свойств матрицы К^ и соотношения п и I, т. е. от числа стоков и нестоковых вершин в сети.

Задачи управления могут быть поставлены сходным образом и для регулярных несимметричных сетей с несколькими аттракторами [4, 5]. Но их решение более сложно. Распределение ресурса между аттракторами происходит по закону индуцированной поглощающей сети (в которой элиминированы все ребра, исходящие из аттракторов). Ресурс распределяется в пропорции, заданной индуцированной матрицей перехода К^. Однако точно эта пропорция соблюдается лишь для ряда начальных состояний. При произвольном начальном состоянии в несимметричной сети существует так называемая «поправка на регулярность», которая позволяет находить предельное значение ресурса в аттракторах аналитически лишь с некоторой погрешностью, которая может быть оценена сверху.

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

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

ЛИТЕРАТУРА

1. Кузнецов О.П. Однородные ресурсные сети. I. Полные графы // Автоматика и телемеханика. — 2009. — № 11. — С. 136—147.

2. Кузнецов О.П, Жилякова Л.Ю. Двусторонние ресурсные сети — новая потоковая модель // Доклады Академии наук, 2010. — Т. 433, № 5. — С. 1—4.

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

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

5. Жилякова Л.Ю. Исследование эйлеровых ресурсных сетей // Управление большими системами. — 2012. — Вып. 40. — С. 28—50.

6. Форд Л-P, Фалкерсон Д.Р. Потоки в сетях. — М.: Мир, 1996. — 334 с.

7. Ahuja R.K., Magnanti T.L., Orlin J.B. Network Flows: Theory, Algorithms, and Applications. — New Jersey: Prentice Hall, 1993. — 864 p.

8. Lovasz L, Winkler P. Mixing of Random Walks and Other Diffusions on a Graph // Surveys in Combinatorics / Ed. P. Row-linson. — London Math. Soc. Lecture Notes Series 218, Cambridge Univ. Press, 1995. — P. 119—154.

9. Blanchard Ph., Volchenkov D. Random Walks and Diffusions on Graphs and Databases: An Introduction (Springer Series in Synergetics). — Berlin—Heidelberg: Springer-Verlag, 2011.

10. Bjorner A., Lovasz, L. Chip-firing game on directed graphs // J. Algebraic Combinatorics. — 1992. — N 1. — P. 305—328.

11. Prisner E. Parallel Chip Firing on Digraphs // Complex Systems. — 1994. — N 8. — P. 367—383.

12. Ivashkevich E.V., Priezzhev V.B. Introduction to the sandpile model // Physica A 254. — 1998. — P. 97—116.

13. Speer, E.R. Asymmetric abelian sandpile models // J. of Statistical Physics. — April, 1993. — Vol. 71, iss. 1-2. — P. 61—74.

14. Bak P. How Nature Works: The Science of Self-Organized Criticality. — N.-Y.: Copernicus, 1996.

15. Bak P., Tang C, and Wiesenfeld K. Self-organized criticality // Physical Review A 38. — 1988. — P. 364—374.

16. Dhar D. Self-organized critical state of sandpile automaton models // Physical Review Letters 64. — 1990. — P. 1613—1616.

17. Кемени Дж., Снелл Дж. Конечные цепи Маркова. — М.: Наука, 1970. — 272 с.

18. Гантмахер Ф.Р. Теория матриц. — М.: Физматлит, 2004. — 560 с.

19. Seneta E. Non-negative Matrices and Markov Chains. — N.-Y.: Springer-Verlag, 2006. — 279 p.

20. Rothblum U.G. Computation of the eigenprojection of a nonnegative matrix at its spectral radius // Stochastic Systems: Modeling, Identification and Optimization II, ser. Mathematical Programming Study / R.J.-B. Wets, ed. — 1976. — Vol. 6. — P. 188—201.

21. Агаев Р.П., Чeботарев П.Ю. О нахождении собственного проектора и компонент матрицы // Автоматика и телемеханика. — 2002. — № 10. — С. 3—12.

22. Чеботарев П.Ю, Агаев Р.П. Об асимптотике в моделях консенсуса. — URL: http://ubs.mtas.ru/bitrix/components/ bitrix/forum.interface/show_file.php?fid=7640 (дата обращения 4.04.2013).

23. Кузнецов А.В., Сакович В.А., Холод Н.И. Высшая математика. Математическое программирование. — Минск: Вы-шэйшая школа, 1994. — 286 с.

24. Boyd, S., Vandenberghe, L. Convex Optimization. — Cambridge: Cambridge University Press, 2004. — 727 p.

25. Лоусон Ч., Хенсон Р. Численное решение задач метода наименьших квадратов. — М.: Наука, 1986. — 232 с.

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

А.С. Манделем.

Людмила Юрьевна Жилякова — канд. физ.-мат. наук,

ст. науч. сотрудник, Институт проблем управления

им. В.А. Трапезникова РАН, г. Москва,

® (495) 334-76-39, И zhilyakova.ludmila@gmail.com.