Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Матричное исчисление с приложениями в теории динамических систем

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
36.1 Mб
Скачать

§ 4.6. Нормальные формы матрицы

Пусть дано расщепление «-мерного пространства R на два под­ пространства: R = Ij + 12, где I, и 12 — инвариантные относительно

А подпространства. Пусть, далее, = (е ^

е ^ ) — базис в 1р

а Шг — (е^2)

е ^ ) — базис в 12. Имеем

 

 

А^ = %А%

(4.6.1)

где %— (<?, # 2), а А — матрица, отвечающая в выбранном базисе & оператору А. Представим эту матрицу в виде блочной матрицы:

где Аи , А12у А2[, А22 — матрицы типа соответственно т х т, т х д , д х m, q х q. Тогда (4.6.1) можно представить в виде:

Отсюда

А * , = аг, Аи + 1Г2Л21,

Ag*2 = ШхА12 + %2А 22.

Так как Ij и 12 — инвариантные пространства, то из получен­ ных равенств ясно, что А21 = 0 и А12= 0, т.е. матрица А оператора А в выбранном базисе имеет квазидиагональную структуру:

Тем же путем легко устанавливается, что если пространство R расщеплено на инвариантные подпространства 1Р 12, ..., 1Р то, на­

брав базис в R из базисов этих подпространств, мы будем иметь матрицу оператора А в таком базисе в виде квазидиагональной матрицы А = diag (Ап >А22, ..., Ап).

4.6Л. Естественные нормальные формы. Согласно теореме 4.5.1 пространство R может быть расщеплено на циклические под­ пространства 1Р 12, ..., I, с минимальными многочленами соответст­

венно

^ J(A.) = + ct11Ami“ 1 4-... + a lmj,

ЧЪ(*) = Хтг + a21Xm2~1+ ... + a2m ,

^(X ) = Xmr + a<1X"*<- 4 - . . . + a

Здесь т1> т2 > > mt, причем каждый многочлен я|к(Я) есть де­ литель предыдущего. Пусть е,, е2, ..., ef — порождающие векторы подпространств I,, 12, ...»1Г Составим базис всего пространства R из

базисов этих циклических пространств: Ж

(I?,, ^ 2, г д

е

= (е, Ае, ... A'"i“ le,),

&2 = (е2 Ае2 ... Ата-2е2), ....

(е, Ае/ ... Ам<"1е,).

 

 

Равенство A2? = 8?А приводит к следующим соотношениям:

 

( i - 1 , 2 ......

t).

А|у-

0

( i* y ) .

(4.6.2)

Из (4.6.2) находим

 

 

 

 

 

А(АИ ’е,) = e(a5j,> + Ае,в^ +

+А '".-‘е,о«,

(4 6 3)

 

fflj) <

 

 

 

Здесь afy — элементы р-го столбца матрицы Аи.

При Ц =я: 1 9

, ftlf -

1 из L6.3) получаем

 

а<Дн- = :1,

ajv- = 0

 

+ 1).

 

При р = w ■, учитывая, что

 

 

 

 

А

=

с'„А”1' “ V

...

aimЛ*

 

будем иметь

 

 

 

и 11 *2

 

 

аСО = - ,2/т(-/+1

• • •>

 

Итак, диагональные блоки Аи квазидиагональной матрицы име-

ЮТ вид

 

'0

0

0

aimf

 

 

 

 

 

 

1

0

0

 

- 1

 

 

Ап =

0

1

0

^im|- 2 •

(4.6.4)

 

 

• • •

 

 

 

 

 

 

 

0

0

1

 

 

 

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

Если принять расщепление пространства согласно теореме 4.5.6, то в соответствующем базисе матрица оператора также будет квазидиагональной матрицей с диагональными блоками вида (4.6.4), только в данном случае характеристический многочлен каждого ди­ агонального блока будет степенью неприводимого в поле Э? много­

члена. Про квазидиагональную матрицу такого типа говорят, что она имеет- вторую естественную нормальную форму.

4.6.2. Жорданова нормальная форма. Пусть пространство R

расщеплено согласно

теореме 4.5.6

на подпространства

...

..., 1^, ... , 1^,..., 1^,

минимальные

многочлены которых

пред­

ставляют собой неприводимые в поле Ж многочлены (см. (4.5.13)). Пусть X — поле комплексных чисел. Тогда эти минимальные мно­ гочлены будут степенями линейных двучленов:

( Х - Х ,) с., ( Х - Х ^ , . . . , (X -X ,)S ,

(Х-Х.УЧ (X —X2) V . . , ( х - х у „

(4.6.5)

(X -X ^ 'i,

(X — Х2У*.......

(X — XJ) /i,

(ск > dk & ... 2*

ск > 0,

к = 1, 2,..., 5).

Возьмем один из многочленов (4.6.5), например,

а -Х о )р ,

где XQ — одно из чисел X,, ..., Х5, а р — один из отличных от нуля показателей ск, ..., 1к (к = 1,.... s). Этот многочлен является ми­ нимальным многочленом определенного циклического подпростран­

ства 1 0

(одного из подпространств 1 ^ , ..., 1 ^ ) .

Пусть е

— порож­

дающий вектор этого подпространства. Тогда векторы

 

=

(А — Х0Е)#>_1е, е2= (А — Х0Е)/>_2е ,...,

= е,

(4.6.6)

где р — размерность подпространства 10, линейно независимы.

Примем систему векторов (4.6.6) в качестве базиса 1 0.

Воздействуя

на векторы (4.6.6) оператором А — ХдЕ, будем иметь

 

(А — Х0Е)е1— О, (А

Х^Е)е2 = е^,..., (A

 

XQE)C^ = е^_ р

Отсюда

 

 

 

 

 

 

Ае^ Хдвр Ае2

Хде2 ■+■

,..., Ае^

XQер -Ь

(4.6.7)

Равенства (4.6.7) можно представить в виде

 

 

 

АГ = ^(Х0£ /, + Я/,),

 

 

 

где

 

 

'0 1 0

0'

 

 

 

 

 

% — (с, е2

ер).

Н =

0 0

1

0

 

• *

 

0 0 0

1

0 0 0

0

( Нр — матрица сдвига*). Таким образом, оператору А в 10 в базисе (4.6.6) отвечает матрица Х0Ер -f Нр. Линейно независимые векто­ ры ер е2, ...» ер, для которых имеют место равенства (4.6.7), обра­

зуют так называемую жорданову цепочку. Из жордановых цепочек, взятых в каждом из подпространств, можно составить базис (жор- данов базис) в R. В этом базисе матрица оператора А имеет жор­ данову нормальную форму

/ = diag {Х.Е,

+ Н . , ..., ХЕ. + ЯЛ.

1 *\

J

J J

Матрицы Ли У, отвечающие одному и тому же линейному опе­ ратору А в разных базисах, связаны друг с другом соотношением подобия:

Л = T J T ~ '.

Если в полном расщеплении пространства на циклические под­ пространства минимальные многочлены всех этих подпространств линейны, то жорданова форма является диагональной матрицей, и

вэтом случае имеем

Л= 7'diag (X,, Х2, ..., Х ^Г "1.

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

З а м е ч а н и е 4.6.1. Из подобия матрицы Л, соответствующей жордановой матрице /, следует, что

det Л = Х£i...Xy*.

§ 4.7. Инвариантные многочлены. Единственность нормальных форм линейного оператора

Через D p (X) обозначим наибольший общий делитель всех мино­ ров р-го порядка характеристической матрицы ХЕ — А (р — 1, 2, ..., п). В ряду

D n(X), ЯП_](Х ),..., i>i(X)

(4.7.1)

каждый многочлен делится на последующий без остатка. Действи­ тельно, минор /-го порядка можно разложить по элементам какойлибо строки. Каждое слагаемое этого разложения есть с некоторым

* Название связано с тем., что произведением некоторой матрицы (Аа... Ап), со­ ставленной из столбцов Al t ..., Ап на матрицу сдвига Нп справа является матрица (О A v Ап^) со «сдвинутыми вправо» столбцами.

множителем минор (/ — 1)-го порядка и потому делится без остат­ ка на D j _ [(X). Следовательно, любой минор /-го порядка, а значит

и

делится без остатка на D j _ l(X).

 

Т е о р е м а 4.7.1. Наибольший общий делитель Dk(X) миноров

к-го порядка матрицы ХЕ А, где А — матрица оператора А в каком-нибудь базисе, не зависит от выбора базиса.

Д о к а з а т е л ь с т в о . Пусть А и А — две матрицы оператора А в разных базисах. Соответствующие характеристические матри­ цы связаны друг с другом соотношением

Х Е - А = Т ~ 1(ХЕ — А)Т.

Покажем сначала, что общие наибольшие делители матрицы ХЕ — А и С ( \ Е — А), где С — произвольная невырожденная мат­ рица, совпадают. Пусть С ~ А = (aiy). Тогда /-я строка мат­

рицы С(ХЕ— А) имеет вид

т.е. является линейной комбинацией строк матрицы ХЕ А с не зависящими от X коэффициентами с.,, сп, ..., cin. Поэтому минор

матрицы

С(ХЕ— А)

разлагается на. сумму

миноров

матрицы

ХЕ — А с некоторыми не зависящими' от X коэффициентами. Сле­

довательно,

всякий

делитель миноров к-то

порядка

матрицы

ХЕ — А

будет делителем миноров к-го

порядка

матрицы

С ( Х Е - А).

 

 

 

 

Точно так же всякий делитель миноров к-то порядка матрицы

С(ХЕ— А)

является

делителем миноров Х-го

порядка

матрицы

С ~1[С(ХЯ— А)] = ХЕ — А. Это означает, что у матриц ХЕ — А и С(ХЕ — А) общие делители миноров к-го порядка = 1, 2, ..., п) совпадают.

Аналогичным образом устанавливается, что у матриц ХЕ — А и С{ХЕ — А), где С — произвольная невырожденная матрица, общие делители миноров к-то порядка совпадают.

Применяя полученный результат к матрицам ХЕ— А и Т ~ 1(ХЕ— А), а затем к матрицам Т ^ ^ Х Е — А) и Т ~ 1(ХЕ— А)Т,

приходим к заключению, что у матриц ХЕ — А и Т~'(ХЕ — А)Т об­ щие делители (в том числе и общие наибольшие делители) миноров Х-го порядка = 1, 2 , п) совпадают. Теорема доказана. ■

Итак, Dk{\) (к = 1, 2,..., п) являются инвариантами линейно­

го оператора и не зависят от выбора базиса в R. Разделив каждый член ряда (4.7.1) на последующий, получим другую группу инва­ риантов оператора А:

*i(X)

Dn™

г2(Х) = Dn-x™

D^X)

(Х ) =

Ш Х ) = 1).

Многочлены

гр(Х) (р —1, 2,...» л)

называются

инвариантными

многочленами. Произведение всех инвариантных многочленов рав­ но характеристическому многочлену:

Д(Х) = |Х Я - А | = D„(X) = П ip(X).

p - i

Разложим инвариантный многочлен i (X) на неприводимые в поле Ж многочлены:

i,tt)

= (f,(X)]r'[92(X)]'>...

 

Здесь (X), ф2(X),

... — различные неприводимые в поле

много­

члены. Степени этих многочленов [<p|(X)]ri, [ф2(Х)]Гг»

[фДХ)]*'»,

отличные от постоянной, называются элементарными делителями характеристической матрицы А или просто матрицы А.

Возьмем теперь в качестве матрицы оператора А матрицу, име­ ющую первую естественную нормальную форму:

А = diag (Ajj, А22>

где Ац (i — 1, 2, ..., t) — матрицы вида (4.6.4) с минимальными многочленами

ф.(Х) = кт14- алХт/-1 + ...

+

 

 

(* =

Определитель матрицы

 

 

 

 

' X

0

0

 

 

-1

X

0

а,™ -1

Х ^ , — Ац =

 

 

1

1

 

 

 

 

0

0

-1

X + а(1

как легко видеть, совпадает с ее минимальным многочленом:

|Х£,П- Аи \ = ф,(Х)

(/= 1, 2......

t).

(4.7.2)

Определитель квазидиагональной матрицы равен произведению определителей диагональных блоков. Поэтому, учитывая (4.7.2), будем иметь

Dn(X) = Д(Х) = ч1>ДХ)1|>2(Х)...4>ДХ).

(4.7.3)

Вычислим далее Dn_ ДХ). Минор (л — 1)-го порядка элемента матрицы diag (ХЕт — Аи , Х Е т — At(), расположенного вне ди­

агональных клеток, равен нулю. Действительно, для получения этого минора нужно вычеркнуть из определителя | ХЕ — А | строку и столбец, на пересечении которых расположен рассматриваемый элемент. Эти линии пересекают два разных диагональных блока. Пусть, например, вычеркнута строка, которая пересекает блок ХЕ — А--. Тогда вертикальная полоса, содержащая остаток блока

/У/

ХЕт AJJ и состоящая

из

столбцов, будет иметь

только

j — 1

ненулевую строку.

Раскладывая

рассматриваемый

опреде­

литель

(л — 1 )-го порядка

на

миноры

/л.-го порядка указанной

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

Рассмотрим теперь минор элемента, принадлежащего одному из диагональных блоков, например, блоку ХЕт — Л/у.. Такой минор

равен

 

...

+](*■) ••• 4\(>.)х(Л),

где х(Х) — определитель блока, полученного из блока ХЕт AJJ

после исключения строки и столбца, на пересечении которых рас­ положен данный элемент. В частности, минор элемента aim равен

( - i ) » / - 4 i W ... ^ _ ,( Х ) ^ +1(Х )... ц>,(Х).

 

Поэтому многочлен

 

^(Х ) ... ^ - . Д Х ^ +ДХ) — 4>,(*)

(4.7.4)

является общим наибольщим делителем миноров (л — 1)-го порядка, отвечающих элементам матрицы ХЕт — АГ}. Так как ч|)ДХ) делится

без

остатка на ^ _ ДХ), то все произведения вида (4.7.4)

( / =

1,2, —, г) делятся без остатка на произведение 'ФгСЧ*з(^-) •••

4>ДХ). Значит,

(4.7.5)

Аналогичным образом можно получить

 

 

0 „ - 2(А) =

Ц<зШ ... N>,(4.....

(4.7.6)

д » - <+|(*)= Ч>((Х),

D„ - iW =

= O l(X) = l.

Используя (4.7.3), (4.7.5) и (4.7.6), находим

 

 

DJX)

 

 

Г>

,(Х)

 

'/>„_,<Х) “

^ 1 ^ ’

 

1г№ = £)„_2(Х) =

^2(^)> •••»

£> / + 1(Х)

= ^ (X ),

г/ +1(Х) =

=

in(X) = 1.

^(>-) = --д - ^ ( Х) -

Таким образом, многочлены грДА.) ( /= 1, 2,...,

t) совпадают с

отличными от единицы инвариантными многочленами оператора

А, а отличные от единицы

многочлены

...

(TAWH* (А = 1, 2 ,. . . . S ) в

разложении

(4.5.13) совпадают с

элементарными делителями оператора А (и соответствующей мат­ рицы А). Отсюда следует единственность (с точностью до порядка расположения диагональных блоков) как первой и второй естест­ венных нормальных форм, так и нормальной жордановой формы матрицы оператора А, так как минимальные многочлены 4>Г-(Х) (i = 1, 2, ..., t) с точностью до порядка расположения диаго­

нальных блоков однозначно определяют все эти нормальные фор­ мы.

З а м е ч а н и е . Выше мы видели, что две подобные матрицы имеют одни и те же инвариантные многочлены. Но тогда эти мат­ рицы подобны одной и той же нормальной (например, жордановой) матрице и потому подобны друг другу. Таким образом, справедлива следующая

Т е о р е м а 4.7.2. Для того чтобы две матрицы с элементами из поля Ж были подобны, необходимо и достаточно, чтобы они обладали одними и теми же инвариантными многочленами.

ПРЕОБРАЗОВАНИЕ МАТРИЦЫ К КВАЗИДИАГОНАЛЬНОМУ ВИДУ И РАЗЛОЖЕНИЕ ЕЕ НА СОСТАВЛЯЮЩИЕ

§ 5.1. Дефект, матричного многочлена

Если в многочлене

/(X ) = а0Хт + а,Х™- 1 + ... + а т _,Х + ат ( а , € Ж )

скалярный аргумент заменить квадратной матрицей А, то получит­ ся матричный многочлен

/( А ) — а0Ат + а {А т~ 1+ ... + а т _ И + а тЕ п.

Вычислим дефект этого матричного многочлена. Будем считать, что 5Z — поле комплексных чисел. Пусть

/— diag {/^(Xj), / 2(^2)» •••»

жорданова матрица, подобная матрице А, где

W = К Е,, +

г. — степень элементарного делителя, отвечающего жордановой

клетке / у-(Яу.) (порядок блока Уу(Ау-)). Тогда вследствие того, что

А = TJT~l,

 

/(А ) = Tf(J)T~l.

(5.1.1)

Здесь / ( / ) = diag {/(/,), / ( / , ) ......./(У,)).

Согласно (5.1.1) дефект матрицы /(А ) равен дефекту матрицы / ( / ) . Дефект же квазидиагональной матрицы / ( / ) равен сумме де­ фектов диагональных блоков. Поэтому, обозначив через dj дефект

блока /(/у), будем иметь следующее выражение для дефекта мат­ рицы /(А):

d = i d j -

(5.1.2)

/ 1

 

Соседние файлы в папке книги