CC BY
37
УДК 512
О НИЖНЕЙ ГРАНИЦЕ ДОМИНИРУЮЩЕГО СОБСТВЕННОГО ЗНАЧЕНИЯ НЕРАЗЛОЖИМОЙ НЕОТРИЦАТЕЛЬНОЙ МАТРИЦЫ
© В.И. Фомин
Ключевые слова: неразложимая неотрицательная матрица; доминирующее собственное значение; нижняя граница; положительный собственный вектор; спектральный радиус.
Предложена некоторая оценка снизу для доминирующего собственного значения неразложимой неотрицательной матрицы.
ВВЕДЕНИЕ
Как известно, аппарат неотрицательных матриц находит широкое применение в различных областях знания, например, при изучении случайных процессов, описываемых цепями Маркова, используются матрицы переходных вероятностей [1]; при построении межотраслевой модели Леонтьева в математической экономике рассматриваются матрицы расходных коэффициентов [2]; при принятии решений в теории игр используются платежные матрицы с неотрицательными элементами [3]; при исследовании малых колебаний упругих систем применяются осцилляционные матрицы [4]. В связи с этим актуальна задача получения оценок спектрального радиуса таких матриц. Верхние оценки можно получать с помощью известного неравенства p(A) < || A || [5] за счет удачного выбора матричной нормы. При получении нижних оценок для p(A) можно использовать теорему Фробениуса-Перрона [4]: неразложимая неотрицательная матрица имеет положительное собственное значение X , которое мажорирует модуль любого другого собственного значения этой матрицы, т. е. X = p(A) (по этой причине X называется доминирующим собственным значением).
В данной работе предлагается одна оценка снизу для доминирующего собственного значения неразложимой неотрицательной матрицы.
ОСНОВНАЯ ЧАСТЬ
Теорема 1. Пусть A = (a¡j ) - неразложимая неотрицательная матрица размера n х n, X - доминирующее собственное значение матрицы A . Тогда верна оценка
Ґ
X > max
1< р < n
max
1< Г ,•••, < n
п«
m = 1
(1)
где Гр+і = ri для любого 1<p<n ; p, r^...,Гр - произвольные целые числа, изменяющиеся в указанных пределах.
Доказательство. Пусть р - некоторое фиксированное целое число, 1<р<п ; г гр - некоторые
фиксированные числа, 1 <г 1,...,гр < п .
Покажем, что
Xр >
г
П <
Гт+л Гт ■
(2)
Возможны два случая.
1) г 1,...,гр таковы, что найдется хотя бы одно такое то = / , 1 < то < р , что
<з_ _ = 0 .
Ч+1 ч
Тогда
У
П
П а
rm + 1 rm rt + 1 rt 11 rm + 1 rm
m = 1 m Ф t
Следовательно, в силу положительности X, p
Xp > II ar r .
1 1 rm + 1 rm
m = 1
2) г 1,...,rp таковы, что
а г ^ 0
rm+1 rm
для любого 1<т <р (напомним, что гр+: = гх).
По определению собственного значения, выполняются равенства
У ' atj Xj = X x¡, 1 < i < n ,
(3)
j=1
m = 1
m=1
1
m+1 'm
n
где X = (Ху) - положительный собственный вектор, отвечающий собственному значению X: Ху > 0,
1< у < п .
Запишем (3) для г = г2 в виде
X
аг гхг + 8, = аг г-------------хг
г2 г1 г1 1 г2 г1 „ г2
аг2 г
где
(4)
Из (9) следует, что
X- <- X,
Гп г
р аг г г1 гр
(10)
Тогда из (8) и (10) следует, что
XX XX
хг <---------------...--------------хг
1 аг г аг г аг г аг г 1
г2 г1 г3 г2 гр гр-1 г1 гр
81 = 2 0-2 уХу > 0 . у =1 у *г
Из (4) следует оценка
Из последнего неравенства следует, в силу положительности х , что
Xр >Пат+1 г
хг <-----------------------Хг
Г\ Т~)
1 аг г
г2 г1
Запишем (3) для г = г3 в виде
(5)
где гр+1 = г1.
Неравенство (2) доказано.
Из (2) следует, в силу произвольности выбора р и г-1,...,гр , что
агз г Хг2 + 82 аг г2 Хг3
агз г
где
(6)
X > шах
1< р < п
1<г г„ < п
82 = 2 Оз уХу > 0 у =1 у * г2
Из (6) следует оценка
X- <------------- X,
г2 гЗ
2 аг г 3 г3 г2
(7)
Тогда из (5) и (7) следует, что
где гр+х = ^ для любого 1 < р < п2 .
Теорема доказана.
Сделаем некоторые замечания по поводу доказанной оценки (1).
1. Вид пределов изменения для р обусловлен тем, что количество элементов произвольного непустого подмножества множества всех элементов матрицы
размера п х п заключено между 1 и п2 .
2. Оценка (1) равносильна оценке вида
X X
Х <--------------Хг3
1 а а
иг г1 аГз г2
Продолжая этот процесс, получим на предпоследнем шаге оценку
Хг < г1
XX X
°г2 г1 °г3 г2 °гргр-1
Запишем теперь (3) для г = г в виде
а _ хг + 8 „ = а _ -хг
г1 гр гр р г1 гр п г1
где
8р = 2 аг1 уху > 0 ■ у =1 у * гр
(8)
(9)
X > шах
1< р<п
1<г ,..., г„ < п
у
П«
(11)
где гр+1 = гх для любого 1<р<п ; р, г1,...,гр - произвольные целые числа, изменяющиеся в указанных пределах, и гг * г у для любых 1< г, у < р , г * у .
Это следует из того, что при отыскании максимума в правой части (1) можно не учитывать любой набор индексов, в котором имеется хотя бы одна пара совпадающих индексов.
Действительно, пусть д - некоторое фиксированное целое число, 1<д<п2 ; Iь...,- некоторый набор индексов, 1 < / 1,..., < п, причем существуют такие
1< г, у < д , г * у , что Гг- = 1у . Пусть, для определенности, г < у .
В силу соотношения (2) имеем
т = 1
X
гт +1 гт
п
Г_ , 1 г.
т + 1 'т
т = 1
Х
Г
X
п г
1р
п
Xq >
Л + ... Л + ... Л+ + Л+ + ... .
г2 Ч TiTi-1 Ti+1Ti tjtj-1 tj+1tj t1tq
В силу равенства tj = tj получаем
Xq > л t... Л , Л t ... Л , a t ... Л / ■ (12)
t2 *1 tjti-1 ti+1tj tjtj-1 tj+1 V t1 tq V '
Заметим, что оценка (12) является следствием следующих оценок:
Xq+! j > л t... л t л t ... л t
l2 Ч IjIi-1 lj+1lj 4 lq
X J-1 >
f ... f ,
li+1 lj lj lj-1 '
которые получаются из (2), соответственно, при Р = q + *'- j , r1 = t1 , r2 = t2 , ..., ri-1 = ti-1, ri = ^ , = tq и ПРи Р = j - '■ , r1 = tj ,
ri+1 = tj+1 ,
'q+i-j
г2 = ^+1 , ■■■, гу-г “ (/-1 .
Следовательно, оценку (12) можно не учитывать при отыскании максимума в правой части (1). Значит, в силу произвольности выбора д и /:,...,I , при отыскании максимума в правой части (1) достаточно учитывать лишь те наборы индексов г1,...,гр , в каждом из
которых все индексы различны между собой. Тем самым доказано, что оценка (1) равносильна оценке (11).
3. Если А - положительная матрица, то в оценке (1) имеет место знак строгого неравенства.
Это следует из того, что для положительной матрицы
где а43 = а13 .
5. Ослабляя оценку (1), можно получить следующие наглядные неравенства:
Х>
Х>
(aji aij)
max (a ji }
1 < i, j < n
max (aji ay ац
1< i, j, I < n
(13)
(14)
Действительно, для p = 2 1
X > max (arr a__ ) ,
- „ r2r1 r1r2 > ’
или, полагая гг = г, г2 = у , получим (13). Аналогично
получается (14).
6. Из (13) следует неравенство
Х>
1
Л-2
П
т. е. доминирующее собственное значение неразложимой неотрицательной матрицы ограничено снизу средним геометрическим ее матричных элементов.
5m = 2аГщ+1 jXj > 0 , Кm<Р ■
j=1
j*rm
4. Существуют неразложимые неотрицательные матрицы, для которых оценка (1) является неулучшаемой. Рассмотрим, например, следующую матрицу:
ЗАКЛЮЧЕНИЕ
Двусторонние оценки доминирующего собственного значения неразложимой неотрицательной матрицы дают представление о расположении ее спектра в некотором кольце на комплексной плоскости.
ЛИТЕРАТУРА
' 0 0 a13 ^ 1
A = a21 0 0 , 2
I 0 a32 0 , 3
где а1з, а21, 032 > 0 . А является неразложимой неотрицательной матрицей. Максимумы, указанные в (1), достигаются для нее при р = 3, г1 = 1, г2 = 2 , гз = 3 и
Волков И.К., Зуев С.М., Цветкова Г.М. Случайные процессы. М., 2000. С. 167.
Никайдо Х. Выпуклые структуры и математическая экономика. М., 1972. С. 124.
Волков И.К., Загоруйко Е.А. Исследование операций. М., 2004. С. 324.
Гантмахер Ф.Р. Теория матриц. М., 2010. С. 380, 339.
Хорн Р., Джонсон Ч. Матричный анализ. М., 1989. С. 359.
Поступила в редакцию 16 мая 2014 г.
Пa
m = 1
П ai+1 i i = 1
Fomin V.I. ON LOWER BOUND OF DOMINANT OWN VALUE OF INDECOMPOSABLE NONNEGATIVE MATRIX Offer some lower bound for the dominant eigenvalue indecomposable nonnegative matrix.
Key words: indecomposable nonnegative matrix; dominant own value; lower limit; positive own vector; spectral radius.
)
2
j
3
3
X
rm + 1 rm
Фомин Василий Ильич, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, кандидат физико-математических наук, доцент, доцент кафедры прикладной математики и механики, e-mail: [email protected]
Fomin Vasiliy Ilyich, Tambov State Technical University, Tambov, Russian Federation, Candidate of Physics and Mathematics, Associate Professor, Associate Professor of Applied Mathematics and Mechanics Department, e-mail: [email protected]
