книги / Матричное исчисление с приложениями в теории динамических систем
..pdf§ 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 |
|