Полные гиперграфы. Спектры лапласианов. Мультиагентные системы Текст научной статьи по специальности «Математика»

УДК 519.179.1 ББК 22.176

ПОЛНЫЕ ГИПЕРГРАФЫ. СПЕКТРЫ ЛАПЛАСИАНОВ. МУЛЬТИАГЕНТНЫЕ СИСТЕМЫ 1

Блюмин С.Л.2

(Липецкий государственный технический университет,

Липецк)

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

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

Введение

Лапласианы орграфов, а также их спектры играют важную роль в теории и приложениях (см., например, [1, 5, 6, 8, 10]). В [2] предложено формирование матриц инцидентности, валентности, смежности и лапласианов оргиперграфов с использованием комплексных корней из единицы подходящих степеней. В данной работе рассмотрено вычисление спектров так сформи-

1 Работа поддержана РФФИ, проект №09-07-00220-а.

2 Семен Львович Блюмин, доктор физико-математических наук, профессор (slb@stu.lipetsk.ru).

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

1. Циркулянты и их спектры

Циркулянты являются частными случаями теплицевых матриц и занимают особое место в областях математики, связанных с разработкой эффективных алгоритмов [3]. Стандартными общими методами вычисления спектров циркулянтов являются [3] применение дискретного преобразования Фурье или использование представления циркулянта многочленом от матрицы циклического сдвига; в ряде частных случаев применимы и более простые способы. Все они могут быть использованы для вычисления спектров лапласианов гиперграфов в тех случаях, когда лапласианы оказываются циркулянтами, что имеет место для полных гиперграфов.

Для целей данной работы удобно определение циркулянта [3], основанное именно на применении дискретного преобразования Фурье, задаваемого матрицей

0-0 Єт єм т є0<т-1) т

^ = т єьо т Є'1 т є1-(т-1) т ( 2Р ,єт = ехрі ^ т

1 3? 1 0 є(т-і)-і т 1 7 7 £ є є

Обратное дискретное преобразование Фурье задается матрицей

у—> — 1 —1 у—г* —1 1 ' Т Т

Fm = m • Fm = m • F , где Fm - сопряженная матрица.

Для квадратной (вообще говоря, комплексной) матрицы А порядка т через Л(А) обозначается диагональная матрица, на диагонали которой расположен спектр матрицы А, т. е. ее собственные значения:

D( A) =

5( A) 0

0 d2( A)

0

0

0 0 ... Sm (A)_

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

Д(С) = Fm • C • F- = m - • Fm • C • F*.

Такое определение удобно как для выяснения циркулянт-ности матрицы (если Fm-AFml не диагональна, то А не цирку-лянтна), так и, в случае ее циркулянтности, для вычисления ее спектра. Впрочем, циркулянтность матрицы может быть выяснена и непосредственно, чисто визуально [3]: каждые ее столбец (строка), кроме первых, получаются из предыдущих путем циклического сдвига вниз (вправо) на одну позицию; в соответствии с этим циркулянт порядка m может быть записан в виде

С =

с

т-1 С0

"т-1 ^т-2 0 _

Наряду с применением дискретного преобразования Фурье спектр циркулянта может быть вычислен следующим способом [3]: каждый циркулянт может быть представлен в виде

С = Е(=0 ср , где Р - матрица циклического сдвига - основная

циркулянтная матрица перестановок с первой строкой (0, ..., 1),

собственные значения которой равны є]т ; тогда

55 (С) = ЕГО С (Єт У > J = 1>-> m •

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

с

с

0

с

с

2

член от циркулянта являются циркулянтами. В то же время эти операции над нециркулянтами могут дать циркулянты. Для диагональной матрицы О матрица Ет-Б-Ет1 является циркулянтом. Все это используется при вычислении спектров лапласианов полных и полных однородных неориентированных и ориентированных гиперграфов, лапласианы которых являются циркулянтами.

2. Полные и полные однородные гиперграфы: матрицы инцидентности, лапласианы и их спектры

Пусть V - конечное множество, \У\ = т. Пусть элементы V е V этого множества каким-либо (произвольным, но фиксированным на протяжении всего рассмотрения) образом помечены, упорядочены, пронумерованы, например, V = (V!, ..., vm}.

Пусть В - булеан этого множества, т. е. полное множество всех его различных подмножеств, состоящих из различных элементов, \В\ = 2т. Пусть В[к] с В - ^-однородный булеан -полное множество всех различных ^-элементных подмножеств, состоящих из различных элементов, множества V, так что В[к] оВ[^] = 0, к ФI, \В[к]\ = Скт, к = 0,1,...,т , и имеет место разложение

В =^т=0В[к], 2т =хт=0С‘„ .

Гиперграф, или система множеств [7], определяется как пара Н = (V, НЕ), где НЕ с В - некоторый набор гиперребер -различных подмножеств, состоящих из различных элементов, множества V, \НЕ\ = п < 2т Гиперребра также предполагаются каким-либо (произвольным, но фиксированным на протяжении всего рассмотрения) образом помеченными, упорядоченными, пронумерованными, например, НЕ = (Не1, ..., Неп}.

Матрица инцидентности 1(Н) = I(V, НЕ) определяется как (0, 1)-матрица размера т х п, строки которой помечены вершинами, столбцы - гиперребрами гиперграфа, а на месте (V, Не) находится 1 тогда и только тогда, когда V е Не. Иначе говоря, 8

столбец матрицы инцидентности, помеченный некоторым гиперребром, является характеристическим вектором этого гиперребра как подмножества множества вершин гиперграфа.

Лапласиан Ь(Н) = Ь( V, НЕ) определяется как квадратная матрица порядка т, строки и столбцы которой помечены вершинами гиперграфа, по формуле Ь(Н) = 1(Н)-(1(Н))Т.

Применение термина «лапласиан» здесь и далее может быть оправдано тем, что он используется, в соответствии с [2], как единое название для произведения матрицы инцидентности на ее транспонированную (в следующем разделе - на сопряженную). В частном случае графа О применение этого термина не соответствует традиционному, так как приводит к соотношению Ь(О) = О(О) + А (О) (где О(О) — матрица валентности, А(О) — матрица смежности), тогда как классический лапласиан равен О(О) - А(О) (определение лапласиана, использованное в следующем разделе, в частном случае графа приводит к классическому определению лапласиана). Матрица типа О(О) + А(О) встречается в теории графов (например, она допускает интерпретацию в связи с реберными графами), но специального названия не имеет. В предлагаемом контексте она может быть названа «лапласианом без учета ориентации», тогда как классический лапласиан - «лапласианом с учетом ориентации».

Полный гиперграф определяется как пара Н = (V, В) и обозначается Н(т), а полный однородный (к-однородный) гиперграф - как пара Н = (V, В [к]) и обозначается Н(т; [к]). Их матрицы инцидентности и лапласианы обозначаются:

1(Н(т)) = 1(т), Ь(Н(т)) = Ь(т),

1(Н(т; [к]))=1(т; [к]), Ь(Н(т; [к]) = Ь(т; [к]).

При этом матрица 1(т) состоит из матриц 1(т; [к]) как из блоков, а лапласианы связаны соотношением

Ь(т) = ^т=0 Ь(т;[к ])

и имеют следующую специальную структуру:

Ь(ш) =

2т-1 2т~2 2т~2

2т~2 2т-1 2т-2

¿(т;[к ]) =

С- Ск-22

т-1 т-2

Стк--22 Стк--11

т-2 т-1

Ск-2 Ск-2

т-2 т-2

Стк--22

т-2

Стк--22

т-2

Стк--11

т-1

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

1(4)=

¿(4)=

"о" "1 0 0 0" "1 1 1 0 0 0" "1 1 1 0

0 0 1 0 0 1 0 0 1 1 0 1 1 0 1

0 0 0 1 0 0 1 0 1 0 1 1 0 1 1

0 0 0 0 1 0 0 1 0 1 1 0 1 1 1

]), к = 0, 1, 2, 3, 4:

[0] [1] [2] [3]

"8 44 4" "0 0 0 0" "1 0 0 0"

4 84 4 0 0 0 0 0 1 0 0

= + +

4 48 4 0 0 0 0 0 0 1 0

4 44 8 0 0 0 0 0 0 0 1

[4]

-1 -1

¿(4;[к]), к = 0, 1, 2, 3 , 4: [0] [1]

"3 1 1 1" "3 2 2 2" "1 1 1 1"

1 3 1 1 2 3 2 2 1 1 1 1

+ + +

1 1 3 1 2 2 3 2 1 1 1 1

1 1 1 3 2 2 2 3 1 1 1 1

2

т-2

2

т-2

2

т-1

[2] [3] [4]

Собственные значения этих лапласианов обозначаются далее через 1j(m; [к]), j = 1, m, к = 0, 1, m. Очевидно,

что

1j(m; [0]) = 0, 1(m; [1]) = 1, j = 1, ..., m.

Спектры лапласианов L(m), L(m; [к]) как циркулянтов могут быть вычислены описанными выше стандартными методами. Однако их частный вид А = (a - b)I + bJ, где I -единичная матрица, J -матрица из единиц, имеющая собственные значения m кратности 1 и 0 кратности m - 1, с учетом того, что A(A) вычисляются двучленом (a - b) + bAj(J), позволяет заключить, что А имеет собственные значения a - b + mb кратности 1 и a - b кратности m - 1.

Это приводит к следующим наборам собственных значений лапласианов:

1(m) = (m + 1)2m-2, Aj(m) = 2m-2, j = 2,...,m ;

1(m; [к]) = Cm-1 + (m -i)^,

Aj (m;[k]) = Ckm- - Ckm-2, j = 2,...,m .

Они связаны соотношениями

Aj (m) = Z" 0 Aj (m; [к]), j = 1,..., m .

Представляют интерес суммы

sm) = zm=1 Aj (m) =m2m-1, s(m; [к]) = A (m; [к]) =mCkm\ ,

связанные соотношением s(m) = Zm=0s(m; [к]).

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

3. Полные и полные однородные оргиперграфы: матрицы инцидентности, лапласианы и их спектры

В случае полного к-однородного гиперграфа H(m; [к]), в соответствии с предложенным в [2] способом перехода от него к полному к-однородному оргиперграфу ОН^; [к]), переход от

гиперребер he к гипердугам ha отражается в матрице инцидентности I(OН(m; [£])) = OI(m; [^]) следующим образом: в столбце, отвечающем гипердуге ha, единицы заменяются, вообще говоря, комплексными корнями степени k из единицы

( 2т"

к-1

имеющими нулевую сумму. Лапласиан ОL(m; [k]) = OI(m; [k])•(OI(m; И))* оказывается при этом эрмитовой циркулянтной матрицей, элементы которой формируются из указанных корней из единицы. Лапласиан OL(m), как сумма таких циркулянтов, также оказывается циркулянтом. Пример, допускающий сравнение с предыдущим и используемый в дальнейшем:

OI(4) =

" 1 1 1 0 0 0 '

-10 0 110

0 -10 -10 1

0 0 -10 -1 -1

OI(4; И), k = 0, 1, 2, 3, 4:

[2] [3]

"0" 1 0 0 0

0 0 0 0

0 0 0 0

1 О і 0 0 0 1

"1110" "1"

ее 0 1 і

е 0 ее -1

0 е е е - і

[0] [1]

[4]

OL(4) =

8 -1 + 2е -і -3 -1 + 2е + і"

-1 + 2е + і 8 1 + 2е - і -3

-3 -1 + 2е + і 8 -1 + 2е - і

-1 + 2е - і -3 -1 + 2е + і 8

"0 0 0 0" "1 0 0 0" " 3 -1 -1 -1"

0 0 0 0 0 1 0 0 -1 3 -1 -1

+ + +

0 0 0 0 0 0 1 0 -1 -1 3 -1

0 0 0 0 0 0 0 1 -1 -1 -1 3

OL(4;[k]), ¿=04,2,3,4: [0]

[1]

[2]

+

где

" 3 2е S + S 2е 1 - i -1 i

2е 3 2s S + S + i 1 - i -1

е + е 2S 3 2S -1 i 1 - i

2е S + S [3] 2S 3 - i -1 i 1 [4]

S = S3 = S3 = expl

(¥:

S =S

1 = S3 = S3 •

Собственные значения этих лапласианов обозначаются далее через m(m), Mj(m; И), j = 1, m, k = 0, 1, m. Очевидно,

что

m(m; [0]) = 0, m(m; [1]) = 1, j = 1, ..., m.

Собственные значения связаны соотношениями M; (m) = Zm=0 M; (m; [k]), j = 1,...,m •

Представляют интерес суммы

tm) =2'm=1 Mj (m), t(m;[k]) = (m;[k]),

связанные соотношением z(m) = Zm=0r(m;[k]).

Кроме того, справедливы соотношения tm) = Sm), tm; [k]) = o(m; [k]).

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

Для вычисления спектров так сформированных лапласианов оргиперграфов могут быть применены описанные выше стандартные методы вычисления спектров циркулянтов.

В работах [8, 10] приведены обширные таблицы спектров лапласианов графов различной структуры. В следующем разделе приведены вычисленные описанными методами спектры лапласианов полных и полных однородных неориентированных и ориентированных гиперграфов для 1 < m < 5. Проиллюстрированы указанные выше соотношения между собственными значениями и их суммами.

4. Таблица спектров лапласианов полных и полных однородных гиперграфов и оргиперграфов: т=1,...,5

т = 1

1(1; [0]) = 0 о(1; [0]) = 0

1(1; [1]) = 1 о(1; [1]) = 1

1(1) = 1 0(1) = 1

М1; [0]) = о т(1; [0]) = 0

т(1; [1]) = 1 <1; [1]) = 1

м(1) = 1 7(1) = 1

т = 2

1(2; [0]) = 1(2; [0]) = 0 о(2; [0]) = 0

1(2; [1]) = 1(2; [1]) = 1 о(2; [1]) = 2

1(2; [2]) = 2, 12(2; [2]) = 0 о(2; [2]) = 2

1(2) = 1, 1(2) = 3 0(2) = 4

М2; [0]) = т(2; [0]) = 0 <2; [0]) = 0

М(2;[1]) = т(2; [1]) = 1 <2; [1]) = 2

М2; [2]) = 0, т(2; [2]) = 2 <2; [2]) = 2

м(2) = 1, т(2) = 3 4 = 2)

т = 3

11(3; [0]) = ... = 1(3; [0]) = 0 о(3; [0]) = 0

11(3; [1]) = ... = 1(3; [1]) = 1 о(3; [1]) = 3

11(3; [2]) = 4, 12(3; [2]) = 13(3; [2]) = 1 о(3; [2]) = 6

11(3; [3]) = 3, 12(3; [3]) = 13(3; [3]) = 0 о(3; [3]) = 3

1(3) = 1(3) = 2, 1(3) = 8 2 1 = 3)

т(3; [0]) = ... = т(3; [0]) = 0 <3; [0]) = 0

т(3; [1]) = ... = т(3; [1]) = 1 <3; [1]) = 3

т(3; [2]) = 0, т(3; [2]) = т(3; [2]) = 3 <3; [2]) = 6

т(3; [3]) = т(3; [3]) = 0, т(3; [3]) = 3 <3; [3]) = 3

т(3) = 1, т(3) = 4, т(3) = 7 2 1 = 3)

т = 4

11(4; [0]) = ... = 14(4; [0]) = 0 о(4; [0]) = 0

11(4; [1]) = ... = 14(4; [1]) = 1 о(4; [1]) = 4

11(4; [2]) = 6, 12(4; [2]) = ... = 14(4; [2]) = 2 о(4; [2]) = 12

11(4; [3]) = 9, 12(4; [3]) = ... = 14(4; [3]) = 1 о(4; [3]) = 12

11(4; [4]) = 4, 12(4; [4]) = ... = 14(4; [4]) = 0 о(4; [4]) = 4

1(4) = ... = 1(4) = 4, 14(4) = 20 о(4) = 32

т(4; [0]) = ... = т(4; [0]) = 0 т(4; [0]) = 0

т(4; [1]) = ... = т(4; [1]) = 1 т(4; [1]) = 4

М(4; [2]) = 0, т(4; [2]) = ... = т(4; [2]) = 4 т(4; [2]) = 12

т(4; [3]) = 0, т(4; [3]) = 4-2*13, тз(4; [3]) = 4, т(4; [3]) = 4 + 2л/3 т(4; [3]) = 12

т(4; [4]) = ... = т(4; [4]) = 0, т(4; [4]) = 4 т(4; [4]) = 4

т(4) = 1, т(4) = 9 - 2л/3, т(4) = 9, т(4) = 13 + 2л/3 2 = 4)

т = 5

11(5; [0]) = ... = 15(5; [0]) = 0 о(5; [0]) = 0

11(5; [1]) = ... = 15(5; [1]) = 1 о(5; [1]) = 5

11(5; [2]) = 8, 12(5; [2]) = ... = 15(5; [2]) = 3 о(5; [2]) = 20

11(5; [3]) = 18, 12(5; [3]) = ... = 15(5; [3]) = 3 о(5; [3]) = 30

11(5; [4]) = 16, 12(5; [4]) = ... = 15(5; [4]) = 1 о(5; [4]) = 20

11(5; [5]) = 5, 12(5; [5]) = ... = 15(5; [5]) = 0 о(5; [5]) = 5

11(5) = ... = 14(5) = 8, 15(5) = 48 о(5) = 80

т(5; [0]) = ... = т(5; [0]) = 0 т(5; [0]) = 0

т(5; [1]) = ... = т(5; [1]) = 1 т(5; [1]) = 5

т(5; [2]) = 0, т(5; [2]) = ... = т(5; [2]) = 5 т(5; [2]) = 20

15 л/З М5; [З]) = 0, т(5; [З])= ——а, 2 4 15 л/З 15 л/З тз(5; [З]) = - + Ї-/, т(5; [З]) = /, 2 4 2 4 15 л/З т(5; З)= — + — а 2 4 т(5; [З]) = З0

т(5; [4]) = 0, т(5; [4]) = 5 ^ л/5 -1 а, 2 т(5; [4]) = 5-л/5+1 /, т(5; [4]) =5 ^л/5 -1 /, 2 2 т(5; [4]) = 5 + л/5 +1 а 2 т(5; [4]) = 20

М(5; [5]) = ... = т(5; [5]) = о , т(5; [5]) = 5 т(5; [5]) = 5

З7 Г л/З + 2 ^х(5) = 1, т(5) = + л/5 а , 24 тз(5) =З7 - л/5 ^+2/, т(5)= З7 -75 -^ +2/, 2 4 2 4 ,,/г\ 47 г- л/З + 2 т(5) = —+ л/ 5 + а 2 4 т(5) = 80

где а = З^Ю + 2л/5 + л/їо-2/5 , / = V10 + 2л/5 -З^10-2/5.

5. Мультиагентные системы: проблемы согласия

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

Решение проблем согласия в мультиагентных системах с учетом парных взаимодействий агентов использует спектр лапласиана взвешенного орграфа, ассоциированного с системой [5]. Пусть мультиагентная система состоит из m агентов, xi(t) -

состояние /-го агента, i = 1, ..., m. В модели достижения консенсуса (согласия) с учетом парных взаимодействий агентов [5]

Х (О = "Хт=1 а/ (1)'[х/(1) - х/ (*)]’/ =1’ -’т ,

число а^) > 0 означает вес, с которым /-ый агент учитывает парное бинарное расхождение значения своего состояния со значением состояния /-го агента, определяемое как разность этих значений, т. е. как их сумма с коэффициентами ±1 или

= ехр^ к ^ = ехррк), к = 1,2,

- корнями степени 2 из 1, имеющими нулевую сумму:

Д( X/ 0); х}. О)) = X/ 0) - х}. О) = (+1)' X/ Ц) + (-1)' х}. О) =

= е22 ' х/ ^) + е2 ' х3 ^).

Лапласиан Ь(?) упомянутого взвешенного орграфа - взвешенный лапласиан указанной модели [5], с использованием которого она записывается в матричной форме

х(^) = - Ь(1) х(^),

по определению, имеет нулевые строчные суммы, а потому сингулярен и его собственному значению 0 отвечает собственный вектор из единиц 1 = [1, ..., 1] . Если 0 - простое собственное значение, то соответствующее собственное подпространство одномерно и вектор 1 является единственным с точностью до некоторой константы с, что и означает полное согласование состояний всех агентов с учетом их парных взаимодействий.

Учет групповых взаимодействий агентов группами, содержащими по 2 < к < т из общего числа т агентов - так что каждый агент непосредственно взаимодействует с к - 1 агентами -осуществляется следующим образом.

Групповое к-арное расхождение значения состояния /-го агента со значениями состояний всех остальных агентов, входящих в группу, определяется как сумма этих значений с коэффициентами е/, / = 1, ., к, - корнями степени к из 1, имеющими нулевую сумму. Иллюстрирующим это примером аналога

приведенной выше модели достижения консенсуса для случая т = 4, к = 3 является следующая модель, в записи которой время ^ опущено, а для весов, с которыми каждый агент учитывает групповое расхождение значения своего состояния со значениями состояний остальных агентов группы, использованы упрощенные обозначения:

х1 =-{а1(х1 + ех2 + ех3) + а2(х1 + ех2 + £х4) + а3(х1 + ех3 + ех4)}, ,х2 = -{¿1( х2 + ех3 + £х4) + Ь2( х2 + ех3 + £х1) + Ъ3( х2 + ех4 + £х1)}, х3 = -{с1( х3 + ех4 +£х1) + с2( х3 + ех4 +ех2) + с3( х3 + ех1 + ех2)}, х4 = -{й?1(х4 + ех1 + ех2) + й2( х4 + ех1 + £х3) + ё3( х4 + ех2 +£х3)}. В матричной форме эта модель принимает вид х = - Ь(4;[3]) х, где использован взвешенный лапласиан данной модели

а1 + а2 + а 3 е(а1 + а2) еа1 + еа3 е(а2 + а3)

Ь(4[3]) = е(Ь2 + Ь3) Ъ + Ь2 +Ъ 3 е(Ь1 + Ъ2) еЪ1 + ^3

ес1 + ес3 е(с2 + с3) с1 + с2 +с 3 е(с1 + с2)

ё(^1 + с12) ейх + е(й2 + ^3) + ^2 +^ 3

Как и взвешенный лапласиан Ь(0 выше, взвешенный лапласиан Ь(4;[3]), очевидно, имеет нулевые строчные суммы, а потому сингулярен и его собственному значению 0 отвечает собственный вектор из единиц 1 = [1, ..., 1] . Если 0 - простое собственное значение, то соответствующее собственное подпространство одномерно и вектор 1 является единственным с точностью до некоторой константы с, что и означает полное согласование состояний всех агентов с учетом их групповых взаимодействий.

Сказанное выполняется в простейшем случае, когда все веса равны единице: в этом случае взвешенный лапласиан Ь (4;[3]) совпадает с лапласианом ОЬ(4; [3]) соответствующего оргиперграфа, приведенным в разделе 3, спектр которого, приведенный в разделе 4, имеет простое собственное число 0 и,

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

Частным случаем групповых взаимодействий агентов (наряду с парными, когда к = 2) являются их совокупные взаимодействия, когда к = т. В этом случае рассмотрение, подобное проведенному выше, приводит к лапласианам ОЬ(т; [т]) соответствующих оргиперграфов. Как показывают примеры их спектров, приведенные в разделе 4, собственное число 0 в этих случаях не является простым, соответствующие собственные подпространства не являются одномерными и обоснование достижения консенсуса, вообще говоря, непосредственно выполнено быть не может.

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

X = -а • (х1 + е • х2 + е • х3), х2 = -Ъ • (х2 + е • х3 + е • х1) ,

Х3 = -с • (х3 + е • х1 + е • х2) ,

которая в матричной форме использует соответствующий взвешенный лапласиан:

х = - Ь (3;[3]) х = -

а е • а е• а

е • Ъ Ъ е • Ъ

е • с е• с с

х.

В простейшем случае, когда все веса равны единице, этот взвешенный лапласиан Ь (3;[3]) совпадает с лапласианом ОЬ(3; [3]) соответствующего оргиперграфа, спектр которого, приведенный в разделе 4, имеет двукратное собственное значение 0, так что соответствующее собственное подпространство является двумерным. Оно описывается соотношением х1 + е • х2 + е • х3 = 0,

( 1 ^л/3) ( 1 .л/5 0

х1 + (-2-^ х2 + (-2-^"2) ^хз = 0,

или

х1 - 2 (х2 + х3) + ¡^2- (х3 - х2) = 0,

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

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

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

В частности, представляет интерес установление взаимосвязи как между лапласианами невзвешенного и взвешенного оргиперграфов, так и между их спектрами. Примером может служить случай к = т, когда

Ь (т; [т]) = А • ОЬ(т; [т]), где А = diag {а, Ъ, с, ..., ё} - диагональная матрица весов модели. Так как матрица ОЬ(т; [т]) имеет, по определению, единичный ранг и спектр 0, ., 0, т, а при умножении диагональной матрицы А слева на матрицу единичного ранга ранг последней не увеличивается, а ее .-ая строка умножается на 7-ый диагональный элемент матрицы А, то матрица Ь (т;[т]) имеет собственные значения 0 кратности т - 1 и а + Ъ + с+ ... +ё кратности 1.

В исследовании мультиагентных систем важную роль играют обобщения [5] матричной теоремы Кирхгофа из теории

графов [4]. В ее основе лежит лемма о равенстве между собой алгебраических дополнений всех элементов матрицы с нулевыми суммами строк и столбцов [4]. Использование комплексных корней из единицы при формировании матриц инцидентности оргиперграфов и их лапласианов обеспечивает последним это ключевое свойство. Представляет интерес получение и интерпретация подобной теоремы в контексте гиперграфов.

Заключение

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

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

Литература

1. БЛЮМИН С.Л. Мультиагентные системы: проблемы и протоколы согласия, псевдообращение лапласианов графов // Системы управления и информационные технологии.

- 2007. - №2(28). - С. 4-9.

2. БЛЮМИН С.Л. Оргиперграфы: матричные представления

[Электронный ресурс]. - Режим доступа:

http://ubs.mtas.ru/forum. - Сетевые модели в управлении. -/2009/10/17/- 19 с.

3. ВОЕВОДИН В.В., ТЫРТЫШНИКОВ Е.Е. Вычислительные процессы с теплицевыми. - М.: Наука, 1987. - 230 с.

4. ЕМЕЛИЧЕВ В А., МЕЛЬНИКОВ О.И., САРВАНОВ В.И., ТЫШКЕВИЧ Р.И. Лекции по теории графов.- М.: Наука, 1990. - 384 с.

5. ЧЕБОТАРЕВ П.Ю., АГАЕВ Р.П. Согласование характеристик в многоагентных системах и спектры лапласовских матриц орграфов // Автоматика и телемеханика. - 2009. -№3. - С. 136-151.

6. BIYIKOGLU T. LEYDOLD J., STADLER P. Laplacian Eigenvectors of Graphs. Perron-Frobenius and Faber-Krahn Type Theorems. - Berlin: Springer, 2007. - 115 p.

7. BOLLOBAS B. Combinatorics: Set Systems, Hypergraphs, Families of Vectors and Combinatorial Probability. - Cambridge: Cambridge University Press, 1986. - 178 p.

8. CVETKOVIC D., DOOB M., SACHS H.. Spectra of Graphs: Theory and Applications. 3rd Ed. - Leipzig: Johan Ambrosius Bart Verlag, 1995. - 368 p.

9. LIPATOV L., RAUSCH DE TRAUBENBERG M., VOLKOV G. On the Ternary Complex Analysis and its Applications [Электронный ресурс]. - Режим доступа:

http://arxiv.org/0711.0809v1. - 28 p.

10. NEWMAN M. The Laplacian Spectrum of Graphs: Master Thesis. - Winnipeg: University of Manitoba, 2000. - 126 p.

COMPLETE HYPERGRAPHS. LAPLACIAN SPECTRA. MULTI-AGENT SYSTEMS

Sam L. Blyumin, Lipetsk State Technical University, Lipetsk, Dr. Sc. in Math. & Phys., professor (slb@stu.lipetsk.ru).

Abstract: Computation of Laplacian spectra is considered for complete and complete uniform undirected and directed hypergraphs which Laplacians are circulants. Standard methods are used for computation of circulant spectra. Applications are discussed to

consensus problems in multi-agent systems taking account of agents’ group interactions.

Keywords: complete and complete uniform undirected and directed hypergraphs, Laplacians, circulants, spectra, multi-agent systems, consensus problems, agents’ group interactions.

Статья представлена к публикации членом редакционной коллегии А. К. Погодаевым