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

УДК 519.71; 62.50

В. М. Осипов, В. В. Осипов

УСТОЙЧИВОСТЬ ТЕПЛИЦЕВЫХ МАТРИЦ И ИХ ПОЛОЖИТЕЛЬНАЯ ОПРЕДЕЛЕННОСТЬ

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

Введем в рассмотрение одностороннюю (треугольную) и-диагональную теплицеву матрицу N х N N > и) вида

1 а о 0 0 0 0 1 0

2а1 а 0 0 0 0 0

2ак 2а1 а0 0 0 0

2а п-1 2ак 2а1 а 0 0 0

о ••• 2а «- 2а к 2а1 а0 0

1 0 0 2а «-1 2ак 2а1 а0 _

«-і .

= а0Е: + 22акгк (: х N)

к=1

(1)

с элементным ^-вектором

а^) = Со1оп[а0, 2а1, ••• 2ак, ••• 2аи-1, 0, ••• 0],(2)

являющимся первым столбцом матрицы, поэтому мы можем записать:

а(м) = ти (а((); г)()), е1(N) = со1оп[1,0, • 0]. (3)

Символом гк (к = 0,1,... (и -1)) обозначим канонические матрицы сдвига размерности N х *), введенные в [1; 2]. Это теплицевы матрицы вида (1) с элементными *-векторами в виде единичных векторов

0, 2. 0, 1, 0, ••• 0

к

(к = 0,п - ГЪ (4)

т. е.

2к = Т«

"0

0 ••• 0

1 о 0

0 1 0 • 0

0 ••• 0 1 0 ••• 0

(: х :) (к = 0, п -1). (5)

При к = 0 очевидным образом получаем единичную матрицу Е: = 20.

Канонические матрицы (5) образуют линейно независимую систему теплицевых матриц, играющую роль базисных элементов в пространстве теплицевых матриц (: х N вида (1), аналогично тому, как система единичных :-векторов (4) выполняет такую же роль в :-мерном векторном пространстве. Очевидно, полиномиальное представление в (1) есть разложение теплицевой матрицы по базису.

Введем Д-норму элементного :-вектора (2)

IIа (: 1 = к|+2«2 Iа к|

к=1

ІІ1

(6)

и максимальную столбцевую норму матрицы (1)

Тп (а

(а(:); 2 ) 1 = |а0| Е:| І1 + 2§| а І2к\I =11 1 к=1 1

(7)

которая подчинена и согласована с векторной Д-нормой [3], а в нашем случае эти нормы будут просто совпадать, так как для столбцевой нормы канонических матриц имеют место равенства

||гк|| = 1 (У к = 0,1,...) (г0 = Е*).

Заметим, что нормы (6) и (7) при и < N не зависят от N - порядка матриц.

Множество теплицевых матриц (Ы х N вида (1) при всяких и < N образуют линейное нормированное пространство. Но в нем определена и коммутативная операция умножения матриц, не выводящая за пределы указанного множества в силу свойств нильпотентности канонической матрицы г N х _^):

гк = 1 У к > N, (8)

причем при всех т + и < N — 1 для введенной матричной нормы выполняется кольцевое свойство:

т (а ^); г Ут 0b(N); 21 < ||у (); г )) (); г) =

= а(:)

Таким образом, указанное множество тепли-цевых матриц вида (1) образует коммутативную нормированную алгебру с единицей (нормированное кольцо) [1; 2]. Роль единицы, естественно, выполняет единичная матрица Е:. Обозначим эту алгебру символом АТ:(2).

Т

к

ні іі1

Всякой теплицевой матрице вида (1) может быть поставлен в однозначное соответствие порождающий полином

Р« (а(:); г )= а 0 + 2«£а кгк

(9)

комплексной переменной г с теми же вещественными коэффициентами, которые имеет матрица (1) в своем полиномиальным представлении по степеням матрицы 2 (: х:), т. е. в своем разложении по каноническому базису (5). Отметим, что именно это представление дало основание в [1; 2] назвать теплицеву матрицу (1) полиномиальной или просто Р-матрицей. Эти термины, однако, мы использовать не будем.

Если полиномы вида (9) как функции комплексного переменного при всяком п < N задать в единичной круговой области |г| < 1 комплексной

плоскости г и снабдить переменную г свойством нильпотентности порядка N, т. е. считать, что

гк = 1 Ук > N (10)

(как уже отмечалось, это характерно для канонической матрицы 2 (: х :)), то максимальное значение модуля комплексного полинома (9), которого он достигнет на границе единичной круговой области |г| = 1, причем и |гк| = 1 (к = 0,п -1), будет равно

Рп а

(а(:); г )

= \\Т„ ((N); 2 ) =1 |а

і п-1і | + 2 2 Iа к

к=1

= 1 Іа (: )|| =

(11)

и совпадает с ^-нормой элементного :-вектора (6) матрицы (1) и такой же ее нормой (7), поэтому может быть названо также нормой полинома (9).

Во множестве Р: всевозможных полиномов вида (9), определенных в круге |г| < 1 и имеющих

степени, не превышающие (: - 1) (п < :), содержится, очевидно, единица как полином нулевой степени и определены операции сложения и умножения на число, а также коммутативная операция умножения полиномов, которая, в силу условия нильпотентности (10) комплексной переменной г, также не выходит за пределы множества Р:, причем для нормы произведения при всех т + п < N — 1 выполняется условие

р ('; Щ ). (■; Щ )І1 <1 р (аиг )| 1 р ('; г ) =

—(N)

= а

аналогичное такому же условию для норм соответствующих теплицевых матриц.

Указанные выше операции обладают всеми свойствами, необходимыми для определения множества PN как нормированной и коммутативной алгебры с единицей (нормированного кольца)

полиномов вида (9), заданных в круговой области |г| < 1 и имеющих степени, не превышающие величины ^ — 1). Обозначим эту алгебру символом АР^щ).

В силу однозначного соответствия

Тп а1

(а ^); 2)

п-1 , у п-1 ,

2)= а0EN + 22 ак2 <------------а0 + 22 акг

к=1

= Ри (а ^); У) (12)

и равенства (11) для норм, алгебра ЛТм(г) тепли-цевых матриц вида (1) становится изометрически изоморфной алгебре ЛPN(.z) соответствующих порождающих полиномов вида (9) [1; 2]. В частности, это означает, что всякие операции над матрицами - элементами алгебры ЛТм(г) - одновременно осуществляются и над соответствующими порождающими полиномами - элементами алгебры АР^У), причем результаты этих операций как элементы своих алгебр сохраняют между собой прежние взаимно однозначные соответствия.

Вернемся к матрице (1). При любой ее конечной размерности и любом наборе { ак } ее конечных элементов, причем а0 ф 0, эта и-диагональная теплицева матрица окажется невырожденной, так как

Бе1;Ти (а ^); г )= а N ф 0 УМ > и.

При а0 = 1 будем иметь матрицу вида (1), но с единичной главной диагональю

Ти(1)(а^);г) = EN + 22акгк (#X N) (N > и), (13)

к=1

определитель которой равен единице при любых N включая и предельный случай N ^ да . Матрицу (1) можно выразить через матрицу (13), рассматривая а0 > 0 в качестве общего множителя всех других элементов:

N

Тп (а

(а(N); 2)

EN + 2 2 ак2к

(а(N); 2) ак = — (к = 1, п -1). (14)

Очевидно, при любых а0 > 0 свойства матриц (13) и (14) совпадают, но матрица (13) с единичной главной диагональю более удобна при рассмотрении этих свойств, поэтому далее вместо матрицы (1) будем использовать матрицу (13), опуская в обозначении все введенные верхние индексы, т. е. считая в (1) а0 = 1 при прежнем обозначении:

Тп (а(); 2 )

п-1

= EN + 2 2 ак2

к =1

(Ы х N) (N > п). (15)

Этой теплицевой матрице будет соответствовать порождающий полином

Р

(а ^); г)

п-1

г )=1 + 2 2 а кЩ

к=1

к

(16)

к=1

1 N 111

п-1

а

0

к=1

0

с единичным свободным членом, заданным в единичной круговой области |г| < 1. При необходимости элемент а 0 Ф1 может быть учтен путем умножения на него всех элементов матрицы (15) и коэффициентов соответствующего полинома (16). Очевидно, теплицевы матрицы вида (15) остаются элементами алгебры ATд(Z), а соответствующие им полиномы (16) - элементами алгебры АР^г), изометрически изоморфной алгебре ATд(Z).

Итак, будем работать с теплицевыми матрицами вида (15), т. е. с теми же матрицами вида (1), но при а0 = 1. Для таких матриц при любых N > п существуют обратные матрицы ^х^ того же вида, но с N диагоналями:

Тп-(С;2)= Ен + 22^ (N х N), (17)

к=1

и полновесными элементными ^векторами Р^) = Со1оп[1,Р1,..., Рк, ..., р N-1 ]= ] (( 2 )) (18) с ^-нормами

||p(N)|| = \т-();2) = 1 + 2|Рк| = "ї-ІРк|, Р0 = 1. (19)

1 1 к=1 к=0

По выражению (18) следует, что элементный ^вектор р(N) обратной матрицы (17) при любых N > п, есть единственное решение уравнения

Тп (її^); 2)(N) = ё1(N). (20)

Это, в частности, означает, что (п — 1) конечных элементов ак (к = 1, п -1) исходной

п-диагональной матрицы (15) определяют все N — 1) элементов рк (к = 1, п -1), р0 = 1 обратной матрицы (17). При всяком конечном N и любом наборе конечных элементов ак (к = 1, п -1) нормы (19) обратной матрицы (17) также являются конечной величиной. Однако при N ^ да сумма в (19) становится числовым рядом с положительными членами:

Шп \Т-(а^);2) = Нш ||р^)|| = ]Т|Рк|, Р0 = 1, (21)

N ^даІІ 111 N ^даІІ 111 к=0

и для конечности /1-нормы матрицы (17) в этом случае уже требуется сходимость ряда в (21), что возможно лишь при определенных ограничениях на члены ряда, т. е. необходимы некоторые ограничения на элементы ак (к = 1, п -1) исходной матрицы (15), которые определяют все свойства элементов р к обратной матрицы (17) при любом N.

Таким образом, может быть выделен класс те-плицевых матриц вида (15) (или вида (1) при а0 = 1), элементы которых ак (к = 1, п -1), кроме конечности обладают еще дополнительными

свойствами, обеспечивающими конечность норм соответствующих обратных матриц при N ^ да .

Введем следующее определение. Одностороннюю n-диагональную теплицеву матрицу (15) назовем устойчивой, если Д-норма обратной ей матрицы (17) при N ^да окажется конечной величиной, т. е. если будет сходиться ряд (21).

Введенное понятие устойчивости односторонней теплицевой матрицы связано с использованием точечных моделей в теории линейных динамических систем [1; 2]. Дело в том, что векторноматричные равенства вида

Tn (a(N);ZX(N) = U(N) о X(N) = т- (((N);Z)(N) (22)

могут рассматриваться как точечные модели линейных динамических систем, описываемых линейными дифференциальными уравнениями с постоянными коэффициентами и связывающих точечные изображающие N-векторы входного

сигнала U(N) и выходного X(N) [1; 2].

Для /j-нормы точечного изображающего N-вектора выхода имеем оценку

||X(N^ < Цт,- (а(N);Z|J|U(N^ (23)

Норма окажется конечной величиной при любых N, включая и предельный случай N ^ да ,

если такие же нормы входного N-вектора U(N) и матрицы T— (a(N); Z) окажутся ограниченными при N ^ да, т. е. если теплицева матрица Tn (a(N); Z) окажется устойчивой в смысле данного определения. Как установлено в [1], это означает устойчивость моделируемой динамической системы в обычном смысле, точнее в смысле «ограниченный вход - ограниченный выход».

Не вдаваясь в детали, отметим, что последнее связано с принципиально важным фактом: гомоморфизм алгебраических структур временных сигналов и их точечных представлений при N ^ да переходит в изометрический изоморфизм. При этом точечные модели динамических систем типа (22) становятся точными [1].

Введем еще одно понятие. Порождающий полином (16) односторонней n-диагональной теплицевой матрицы (15) будем называть также характеристическим полиномом этой матрицы, а алгебраическое уравнение степени (и — 1)

Pn (a(N); z ) = 1 + 22 a kzk = 1 + 2a1z + 2a 2 z 2 +... +

k=1

+ 2a kzk +... + 2a n-1zn-1 = 0 (24)

- ее характеристическим уравнением.

Уравнение (24) не зависит от размерности (N х N) матрицы и определяется ее элементами ak как коэффициентами, причем an-1 ф 0. Число

диагоналей матрицы (15), а это и < N, определяет степень (и — 1) уравнения (24). А свойства характеристического полинома, естественно, определяют свойства соответствующей теплицевой матрицы.

Теорема 1. Для устойчивости односторонней и-диагональной теплицевой матрицы вида (15) необходимо и достаточно расположения всех корней ее характеристического уравнения (24) за пределами единичной круговой области, т. е. в области |и| > 1.

Доказательство. Оно может быть наиболее эффективно выполнено на основе фундаментального положения об изометрическом изоморфизме алгебры ЛТ^г) теплицевых матриц вида (1) (или (15)) и алгебры ЛPN(.z) соответствующих порождающих (характеристических) полиномов. Как уже отмечалось, это означает однозначное соответствие всех алгебраических операций, осуществляемых над соответствующими элементами этих алгебр. Для операции обращения можем записать:

N-1

Т-ф(N);2) = | EN + 22'ак2к

< 2 ^г >| 1 + 2 2 акгк к=1

= рп

1 (а(*); г)

п-1

(25)

к=1

Если характеристический (порождающий) полином

Рп

(а(д); г)

=1 + 2 2 а кг

(26)

заданный в единичной круговой области И <1,

не имеет в ней нулей, т. е. уравнение (24) не имеет корней в этой области, то обратная функция Р— (к(д); г) комплексного переменного и окажется

аналитической в круговой области |и| < 1 и, следовательно, может быть представлена в ней степенным рядом Тейлора при г = 0:

Р- (а ('; г )-■

1

п-1

2 2 с

к=1

1 + 2 2 а кщк

=1+2 в г

к=1

(27)

Ряд сходится абсолютно для всех И из единичного круга И| < 1, включая и его границу |и| = 1,

поэтому для максимума модуля функции (27), принятого за его норму, при любых конечных N будем иметь

||р;' (а(к >; г), = НРи"’ (к ’; г) = 1 + еК I (28)

II III |г|=^ I к=1

В силу изометрии алгебр, вместо соответствия (25) получим равенство для норм указанных элементов

,2 ]) =|р- (а(д)-

= 2 |Рк |; в0 = 1.

(29)

к=0

Очевидно, что в силу абсолютной сходимости ряда в (27) при |и| = 1 в предельном случае при N ^ да окажется

1іш ІІТ71 (а ^); 2 ) = 1іш| \р- (а(д); г)

Nп ' 1І1 Nп ' 1І1

N-1

= ііш 2 Iе к| = 2 Iе к| <да.

(30)

к=0

к=0

По определению, это означает устойчивость матрицы Ти (а(д); г). Теорема доказана.

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

Вернемся к характеристическому полиному (26) теплицевой матрицы (15). Найдем его представление на единичной окружности |и| = 1 - граничной кривой круговой области |и| < 1, в которой определен полином. Положим

г = е'п0 (-1 < 0 < 1) и отделим вещественную и мнимую части:

Рп а

(а(д); вм )=

п-1

-л + 2 2 а ке“

к=1

= Рп(0) + О, (9), (31)

где

п-1

Рп (9) = Ке| 1 + 2 2 а кЄи к=1

Оп (9) = 1ш| 1 + ^ке‘

п-1

= 1 + 2 2 а к сов кж9; к=1

п-1

= 22ак віп кп9.

к=1

(32а)

(32б)

Вещественная часть Ри (0) есть четный тригонометрический полином порядка п с такими же коэффициентами, какие имеет порождающий полином (26). Это означает, что тригонометрический полином (32а) способен характеризовать свойства как порождающего полинома (26), так и свойства соответствующей матрицы вида (15), в частности ее устойчивость в смысле данного ранее определения.

Теорема 2. Если четный тригонометрический полином (32а), как 2-периодическая функция окажется положительным, то теплицева матрица (15)будет устойчивой.

Доказательство. Согласно принципу аргумента теории функций комплексного переменного [4], при полном обороте вектора г = в”10 вокруг начала координат плоскости г в положительном

1

п

к=1

к=1

направлении (против часовой стрелки) вектор (31), расположенный в комплексной плоскости (Pn, iQn), сделает столько полных оборотов вокруг начала координат этой плоскости, сколько нулей имеет полином (26) внутри единичного круга |z| < 1. Предположим, что при 0 = 0 (z = 1) тригонометрический полином (32а) как вещественная часть вектора (31) оказался положительным. Тогда мнимая часть (32б) как нечетная функция будет равна нулю и вектор (31) займет горизонтальное положение, совпадающее с положительным направлением оси Pn комплексной плоскости (Pn, iQn):

Pn (a(N); 1)= Pn(0) = 1 + 2]-1ak > 0. (33)

k=1

Если при изменении переменной 0 от 0 до 2 тригонометрический полином (32а) сохранит свой положительный знак (вектор z = е'п0 при этом повернется на угол 2п, делая один оборот вокруг начала координат и перемещаясь своим концом по единичной окружности в положительном направлении), то это будет означать, что годограф вектора (31) останется в пределах правой полуплоскости комплексной плоскости (Pn, iQn) и не пересечет ее мнимой оси iQn, т. е. никакого вращения этого вектора вокруг начала координат происходить не будет. Согласно принципу аргумента, это означает отсутствие нулей у характеристического полинома (26) в единичной круговой области |z| < 1 и, следовательно, устойчивость теплицевой матрицы Tn (a(N);Z), согласно теореме 1. Теорема доказана.

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

(Tn (a(N)()(N),U(N))= (U(N),T +(a(N))(N)) (34)

здесь для сокращения записи использовано обозначение Tn (a(N) ) = Tn (a(N); Z)

Будем предполагать, что произвольный

N-вектор

U(а) = Colon[U1, ••• Uv, ••• /:

(ЗЗ)

при любых N > п имеет конечную Д-норму. Равенство (34) возникает в силу свойства скалярного произведения, которое порождает транспонированную верхнетреугольную теплицеву матрицу

т;(т (а); Z )= T+Z(N); Z)

Z = E,

+ 22a k (Z

(N x N) (N У n).

(36)

Равенство (34) позволяет представить квадратичную форму в виде

0т. Z( " ))[/ (: ),u (а) )=1

Z " та:

+((а(а ("/и(N ),U(N)

(N) /(N)

= | 1T (a(N)+T.+(z)a>))(N),U

:(sTn (a(N) U(N ),u (N)), (ЗУ)

где

ST. (a(N))=1 [t. Z(N))+c"(N))]=

n-1

= En + Z a k

k=1

Zk +

0z 'Ї

(ЗВ)

есть симметричная (2и - 1)-диагональная тепли-цева матрица ^ х Щ ^ > и) с единичной главной диагональю.

Квадратичные формы в (37), построенные на указанных выше теплицевых матрицах, будут теплицевыми формами четного тригонометрического полинома (32а) [5], причем при всех N > и имеет место представление

(рТи (а^))у ^),й^))= (Ти (а(д))^),й^))=

NN

= ZZ а k-mU kU m.

k=1m=1

(39)

Если квадратичные формы в (39) положительны при любых ^мерных векторах (35) с конечными /1-нормами, то они называются положительно определенными. Это понятие распространяется и на теплицевы матрицы, на которых строятся указанные квадратичные формы.

Известен следующий результат [5; 6].

Утверждение 1. Для положительного четного тригонометрического полинома

.-1

P. (0) = 1 + 2 Z ak cos k%0

k=1

(4О)

необходимо и достаточно положительной определенности теплицевых форм (39), т. е. выполнения условия

(т. (a(N))U(N),U(N)) = (sr. (а(а))(N),U(а))> о (41)

для всех N У n .

Обратное утверждение также верно: положительность полинома (4О) с необходимостью означает выполнения условия (41).

Это утверждение и теорема 2 позволяют сформулировать следующий результат.

Теорема З. Для устойчивости односторонней n-диагональной теплицевой матрицы (1З) необходимо и достаточно выполнения условия (41), т. е. положительной определенности при всех N У n самой матрицы (1З) и, что равносильно, симметричной теплицевой (2. — 1)-диагональной матрицы (ЗВ).

19В

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

Следствие. Если односторонняя п-диагональная теплицева матрица (15) устойчива, то при всех N > п будет устойчивой и обратная ей матрица

T«-1(а(д); 2)= Eд + 2Г2Рк2к (д х д). (42)

к =1

Рассмотрим векторно-матричное уравнение с невырожденной теплицевой матрицей (15) и произвольным Д-вектором и(N) в правой части:

Т« (а(N);2)(N) = и(N). (43)

Умножим скалярно обе стороны этого уравнения на Д-вектор X (д):

(т« (а(д); 2 )х(д), X(д) )= (и(д), X(д))

Если выполняется условие

(т« (а(д);2)х(д),X(д))=

= (и(N),X(N))> 0 (УД > «), (44)

то, согласно теореме 3, матрица Т« (а(д); 2) будет устойчивой. Единственное решение уравнения (43)

X(д) = Т- (a(д); 2 )и(д) (45)

умножим скалярно на Д-вектор и(N):

(X(д), и (д)) = (((д), X(д)) = (т- (а(д); 2 )(д), и (д))

В силу (44)окажется

(и(д), X(д)) = (г-1 (а(д); 2 )• и(д ),и(д)) > 0

(Уд > «), (46)

что означает положительную определенность матрицы (42) и, следовательно, ее устойчивость, согласно теореме 3.

Замечание. Утверждение следствия непосредственно вытекает и по факту однозначного соответствия алгебраических операций, осуществляемых над элементами изоморфных алгебр АТд(2) и АРд(г), т. е. над теплицевыми матрицами вида (15) и их порождающими (характеристическими) полиномами. Это фундаментальное положение было использовано при доказательстве теоремы 1.

Для операции обращения односторонней теп-

лицевой матрицы Т« (х(д); 2) будем иметь следующее представление для характеристических полиномов:

Р- (а(д); г) =--«I----= 1 + д2 ркгк, (Уд > «)..

1 + 2 2 акгк к=1

к=1

Оно возможно, если полином

Pn (a(N);z)= 1 + 22akzk не будет иметь нулей внут-

k=1

ри единичной круговой области |z| < 1, т. е. если

соответствующая матрица Tn (a(N); Z) устойчива.

Формально имеет место и равенство обратных сумм

1 n-1 k

——--------=1 + 2 2 akz , (VN > n),

1 +2 P kzk k=1

k=1

которое также возможно лишь при условии отсутствия нулей внутри единичного круга |z| < 1 у

полинома P-1 (a(N); z) при всех N > n, что означает устойчивость обратной матрицы T- (a(N); Z).

Итак, свойство устойчивости односторонних теплицевых матриц вида (1) (или (15)), в смысле данного ранее определения, эквивалентно свойству положительной определенности как самих этих матриц, так и их симметричных двойников (38), что сводится к требованию положительности соответствующих квадратичных (теплицевых) форм при всех N > n, т. е. к условию

(t (a(N) )u (N ),U(N) )=(sTn (a(N) )u (N ),U(N) )> 0. (47)

Непосредственная проверка этого условия означала бы либо определение всех (вещественных) собственных значений симметричной матрицы STn (a(N)) при различных N > n, которые все должны быть положительными, либо проверку по критерию Сильвестра. И то, и другое в общем случае оказывается задачей аналитически неразрешимой. Это означает крайнюю неконструктивность условия (47). Значительно более конструктивным является утверждение теоремы 2, согласно которому устойчивость матрицы (15) и, следовательно, ее положительная определенность гарантируется положительностью тригонометрического полинома (32а), т. е. выполнением условия

n-1

Pn(0) = 1 + 22 ak cos Ы0 > 06 (1 < 0 < 1)= (48)

k=1

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

Теорема 4. Если Д-норма односторонней n-диагональной теплицевой матрицы (15) удовлетворяет условию

T (a(N); Z) 1 = 1+22|a^ < 2, (49)

1 k=1

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

Доказательство. Попытаемся найти те ограничения, которые нужно наложить на элементы ак симметричной теплицевой матрицы (38), чтобы все ее вещественные собственные значения оказались положительными, что означало бы положительную определенность матрицы, положительность тригонометрического полинома (32а) и, следовательно, устойчивость матрицы (15).

Все круги Гершгорина для симметричной теплицевой матрицы (38) сливаются в один круг

|г -1 < 22|ак| = Я (50)

к=1

с центром на вещественной оси комплексной плоскости г в точке г = 1 и радиусом

«-1

2 2 Ы = Я. (51)

к=1

Круг (50) отсекает на вещественной оси отрезок [1 - Я, 1 + Я], в котором и располагаются все вещественные собственные значения симметричной теплицевой матрицы 8Т« (а(д)). Других собственных значений матрица не имеет. Очевидно, указанный отрезок будет целиком располагаться на положительной половине вещественной оси и все собственные значения матрицы окажутся положительными, если будет выполняться условие

n-1

1 -R > 0 о R = 22|ak| < 1, (52)

k=1

что эквивалентно условию (49) для Д-нормы матрицы (15). Теорема доказана.

Отметим, что условие (52) определяет множество положительно определенных матриц вида (38) со строгим диагональным преобразованием. Условие (49) или (52) легко проверяемо для любой теплицевой матрицы вида (15), но является довольно жестким ограничением, поскольку оно выделяет узкий класс устойчивых теплицевых матриц.

Библиографический список

1. Осипов, В. М. Моделирование линейных динамических систем методом точечных представлений / В. М. Осипов, В. В. Осипов. M. : МАКС-Пресс, 2005. 296 с.

2. Osipov, V. M. Pointwise representation method / V. M. Osipov, V. V. Osipov // Electronic J. of Differential Equations. 2005. № 12. P. 103-116.

3. Полиа, Г. Задачи и теоремы анализа / Г. Полиа, Г. Сеге. М. : Наука, 1978. Ч. 2. 431 с.

4. Ахиезер, Н. И. Классическая проблема моментов / Н. И. Ахиезер. М. : Физматгиз, 1961. 310 с.

5. Хорн, Р. Матричный анализ / Р. Хорн, Ч. Джонсон. М. : Мир, 1989. 655 с.

6. Маркушевич, А. И. Краткий курс теории аналитических функций / А. И. Маркушевич. М. : Гостехтеоретиздат, 1957. 335 с.

V. M. Osipov, V. V. Osipov STABILITY OF TOEPLITZ TYPE MATRIXES AND THEIR POSITIVE DETERMINATION

For the first time the term and the meaning of Toeplitz type matrixes stability are given in. It’s connected with a stability of linear dynamic systems phenomenon, which models worked out on the point method have Toeplitz type matrixes as a transmission. It is also set the connection of the matrixes with the known notion ofpositive matrixes determination.