книги / Матричные и асимптотические методы в теории линейных систем
..pdfРавенства (6.7) можно представить в виде
А% = %(Х0ЕР + Нр),
где
g = (ег |
е 2 |
ер), |
|
~0 |
1 |
0 |
0“ |
0 |
0 |
1 |
О |
|
1 |
0 0 0 |
О |
Щр — матрица сдвига порядка р).
Таким образом, оператору А в / 0 в базисе (6.6) отвечает матрица Х0ЕР+ Нр. Линейно независимые векторы elt
...» ер, для которых имеют место равенства (6.7), образуют так называемую жорданову цепочку. Из жордановых цепо
чек, взятых в каждом |
из подпространств, можно составить |
|
базис (жорданов базис) |
в /?. В этом базисе |
матрица опера |
тора А имеет жорданову нормальную форму |
||
J = diag \Х1Ее1+ Я С1, ... » XsEl$ |
Hls). |
Матрицы Л и У, отвечающие одному и тому же линейно му оператору А в разных базисах, связаны друг с другом соотношением подобия:
А = TJT~\
Если в полном расщеплении пространства на цикличе ские подпространства минимальные многочлены всех этих подпространств линейны, то жорданова форма является диа гональной матрицей, и в этом случае имеем
А = Т diag (Л»1г Х,а> . . . , Xs) Т *.
Таким образом, линейный оператор А имеет простую структуру тогда и только тогда когда, пространство R рас щепляется на инвариантные подпространства с линейными минимальными многочленами.
З а м е ч а н и е . Из подобия матрицы Л, соответствую щей жордановой матрице J, следует, что
det А — Х{‘ . . . A.J*.
§ 7. Инвариантные многочлены. Единственность нормальных форм линейного оператора
Через Dp (X) обозначим наибольший общий делитель всех миноров р-го порядка характеристической матрицы
ХЕ — А (р = 1,2, ..., |
п). В ряду |
|
|
Dn(X), |
Z?n_j(A»), |
D,(X) |
(7.1) |
каждый многочлен делится на последующий без остатка. Действительно, минор /-го порядка можно разложить по элементам какой-либо строки. Каждое слагаемое этого раз ложения есть с некоторым множителем минор (/ — 1)-го порядка и потому делится без остатка на Dj-\ (X). Следова тельно, любой минор /-го порядка, а значит и D/ (Я), делится
без остатка |
на |
D /^i (X). |
Т е о р е |
м а |
7.1. Наибольший общий делитель Dk (Я,) |
миноров k-го порядка матрицы ХЕ — А, где А — матрица оператора в каком-нибудь базисе, не зависит от выбора базиса.
Д о к а з а т е л ь с т в о . Пусть А и А — две матрицы оператора А в разных базисах. Соответствующие характе ристические матрицы связаны друг с другом соотношением
ХЕ— А = Т~1(ХЕ — А) Т.
Покажем сначала, что общие наибольшие делители мат риц ХЕ — А и С (ХЕ — А)у где С — произвольная невы рожденная матрица, совпадают. Пусть С = (ctj)t А = (ац). Тогда 1-я строка матрицы С (ХЕ — А) имеет вид
т. е. является линейной комбинацией строк матрицы ХЕ — А с независящими от X коэффициентами cilt с[2, ..., с1п. Поэто му минор матрицы С (ХЕ — А) разлагается на сумму мино ров матрицы ХЕ — А с некоторыми не зависящими от X коэффициентами. Следовательно, всякий делитель миноров
k-ro порядка матрицы |
ХЕ — А будет делителем |
миноров |
k-го порядка матрицы |
С (ХЕ — А). Точно так же |
всякий |
делитель миноров k-ro порядка матрицы С {ХЕ — А) явля
ется делителем миноров k-ro порядка матрицы С-1 [С (КЕ —
— А)] = ХЕ— А. Это означает, что у матриц ХЕ — А и С {ХЕ — А) общие делители миноров k-ro порядка {k = = 1, 2, п) совпадают.
Аналогичным образом устанавливается, что у матриц ХЕ — А и {ХЕ — A) Ct где С — произвольная невырожден ная матрица, общие делители миноров k-ro порядка совпа дают.
Применяя полученный результат к матрицам ХЕ — А и
Т~х {ХЕ — А), а |
затем к |
матрицам Т~1{ХЕ — А) и |
Т~1 {ХЕ — А) Ту |
приходим |
к заключению, что у матриц |
ХЕ — Л и Т~х{ХЕ — А) Т |
общие делители (в том числе |
иобщие наибольшие делители) миноров k-ro порядка {k =
—1, 2, ..., п) совпадают. Теорема доказана.
Итак, Dk {X) {k = 1, 2, п) являются инвариантами линейного оператора и не зависят от выбора базиса в R. Разделив каждый член ряда (7.1) на последующий, получим
другую группу |
инвариантов оператора А: |
|
|||
h № |
Р , « |
|
|
№ = |
|
P„-iW |
’ |
о„_2М ’ |
|||
|
|
||||
|
|
|
O. (Ц |
(О0 (X) = 1). |
|
|
|
|
P . W |
Многочлены ip (X) (р — 1, 2, ..., п) называются инва риантными многочленами. Произведение всех инвариантных многочленов равно характеристическому многочлену:
Д(Х) = | Х £ _ Л | = Д ,(Х )= Ш ,(Ч -
|
|
|
р=1 |
Разложим инвариантный многочлен ip {X) на неприво |
|||
димые в поле Ж многочлены: |
|
||
|
h (Ч = [Ф1 (*)Г‘ [Ф. № • |
(Ф, № * |
|
Здесь |
cpi {X), |
фа (А-)» ..., ф* (Я) — различные неприводимые |
|
в поле |
ЗС многочлены. Степени этих |
многочленов 1фх (А,)]г‘, |
|
[фа (А.)Г', ...» |
[ф, {X)Ts, отличные от постоянной, называются |
элементарными делителями характеристической матрицы ХЕ — А или просто матрицы А .
Возьмем теперь в качестве матрицы оператора А матри цу, имеющую первую естественную нормальную форму:
А = d iag (Ли , Л22, . . . . Att),
где Ац (i = 1, 2,...» О — матрицы вида (6,4) с минималь ными многочленами
Ф* (X) = А**-J- ССл^ * ~f* |
~f- Ct/ т^—lA»H“ |
|
|
(i = 1, |
2, |
Определитель |
матрицы |
|
|
|
/А , |
0 |
|
Р 1 ь. |
II |
1 |
>» |
|
\ |
0 |
0 |
• **> 0*
0 |
GCi mi |
0 |
ОС/ т^—1 |
— 1 |
К ОС/1 |
как легко видеть, совпадает с ее минимальным многочленом:
IХЕт ~ А „ \= (Ц |
( 1 = 1 . 2 ......... |
t). |
(7.2) |
Определитель квазидиагональной матрицы равен про изведению определителей диагональных блоков. Поэтому, учитывая (7.2), будем иметь
Д , (А,) » А (А,) = |
Ъ (Я) ф2 А) |
(X). |
(7.3) |
Вычислим далее D„_i (X). |
|
|
|
Минор (п — 1)-го |
порядка |
элемента |
матрицы |
diag (ХЕт, — Лп , .... %Emt — Л„), расположенного вне диа гональных клеток, равен нулю. Действительно, для получе ния этого минора нужно вычеркнуть из определителя \КЕ—А\ строку и столбец, на пересечении которых расположен рассматриваемый элемент. Эти линии пересекают два раз ных диагональных блока. Пусть, например, вычеркнута строка, которая пересекает блок ХЕт{ — Ац. Тогда верти кальная полоса, содержащая остаток блока XEmf — Ац
и состоящая из т / столбцов, будет иметь только ту — 1 ненулевую строку. Разлагая рассматриваемый определи тель (п — 1)-го порядка по теореме Лапласа на миноры ту-го порядка указанной вертикальной полосы, убеждаемся, что он равен нулю (каждое слагаемое в этом разложении содер жит определитель, у которого одна строка сплошь заполне на нулями).
Рассмотрим теперь минор элемента, принадлежащего од
ному из диагональных |
блоков, |
например, блоку ХЁт /— Ац. |
||||
Такой минор равен |
|
|
|
|
|
|
Фх (Я) |
Ф/-1(Я) Фж |
(Я) |
ф* (Я) х (Я), |
|
||
где х (Я) — определитель блока, полученного |
из блока |
|||||
ХЕт. — А// после исключения |
строки и столбца, |
на пересе |
||||
чении которых расположен данный элемент. |
|
|
||||
В частности, минор элемента aimравен |
|
|
||||
(— l)w/-1 t i |
(Я) |
Ф/-1 (Я) фж |
(Я) |
ф, (Я). |
||
Поэтому многочлен |
|
|
|
|
||
'Фх (Я) |
|
ф/_1(Я) ф/-|_! (Я) |
М Я ) |
|
(7.4) |
является общим наибольшим делителем миноров (п — 1)-го порядка, отвечающих элементам матрицы ЯЕт. — Ац.
Так |
как ф< (X) делится |
без остатка на |
ф/_1 |
(Я), |
то |
все |
|||||
произведения вида |
(7.4) (/ |
= 1, 2, |
..., |
t) делятся |
без остат |
||||||
ка на произведение ф2 (Я) Фз (Я)... ф* (Я). Значит, |
|
|
|
||||||||
|
|
|
1 (Я) — ф2 (Я) |
|
|
ф/ (Я). |
|
|
(7-5) |
||
Аналогичным образом |
можно |
получить |
|
|
|
|
|||||
|
|
Dn- 2 (Я) = |
фз (Я) |
|
Ф/ (Я), |
|
|
|
|
||
|
|
Dn—t+i (Я) = |
фi (Я), |
|
|
|
|
|
^ |
^ |
|
|
|
Dn-t (Я) = |
= |
|
(Я) ■= 1. 1 |
|
|
|
|||
Используя (7.3), (7.5) и (7.6), находим |
|
|
|
|
|||||||
‘i W |
= |
~ ’I’1 W* |
‘2W |
= |
D"_J (*,) |
“ ♦» W .......... |
|||||
tt w |
= |
°в ~ ^т = |
’h w . |
'<+■ w |
= |
= |
^ w |
= |
»• |
|
|
Таким образом, |
многочлены ф< (Я) |
(/ = 1, 2, |
..., |
/) сов |
падают с отличными от единицы инвариантными многочле нами оператора А, а отличные от единицы многочлены
[фА(Я)Г*, [<р* (Я) ] \ |
.... [ф*(Я))'й (k = |
1,2...... s) в разложении |
(5.13) совпадают с |
элементарными |
делителями оператора |
А (и соответствующей матрицы А). |
|
Отсюда следует единственность (с точностью до порядка расположения диагональных блоков) как первой и второй естественных нормальных форм, так и нормальной жор-
дановой формы |
матрицы |
оператора А , |
так как минималь |
|
ные многочлены |
% (Я) (/ |
= 1, 2, |
t) |
с точностью до по |
рядка расположения диагональных блоков однозначно опре деляют все эти нормальные формы.
З а м е ч а н и е . Выше мы видели, что две подобные мат рицы имеют одни и те же инвариантные многочлены. Но тог да эти матрицы подобны одной и той же нормальной (на пример, жордановой) матрице и потому подобны друг другу. Таким образом, справедлива следующая
Т е о р е м а 7.2. Для того чтобы две матрицы с эле ментами из поля di были подобны, необходимо и достаточно, чтобы они обладали одними и теми оке инвариантными многочленами.
Г л а в а V
ПРЕОБРАЗОВАНИЕ МАТРИЦЫ
ККВАЗИДИАГОНАЛЬНОМУ ВИДУ
ИРАЗЛОЖЕНИЕ ЕЕ НА СОСТАВЛЯЮЩИЕ
§ 1. Дефект матричного многочлена
Если |
в многочлене |
|
|
f (^) “ |
«(Лт 4~ сс1Хт~^ -|- |
-}- ат—j?i -Г сст |
(а/ £ |
скалярный аргумент заменить квадратной матрицей А, то получится матричный многочлен
f (^) — а0^ 4" |
4~ |
4- a m —i A -f ат£л. |
Вычислим дефект этого матричного многочлена. Будем считать, что Di — поле комплексных чисел. Пусть
|
J = diag |
J2 (^2)* • • •» |
(К)} |
— жорданова матрица, подобная матрице А, где |
|||
|
|
O^i) = ^ i ^ r j 4~ H r h |
|
г/ — степень элементарного делителя, |
отвечающего жор- |
||
дановой клетке Jj (X/) (порядок блока Jj (kj)). |
|||
Тогда вследствие того, что А = TJT~\ |
|||
|
f(A)=Tf(J)T~'. |
(1.1) |
|
Здесь |
f(J) = diag I/(Л), / W , |
/ W ) - |
|
Согласно (1.1) дефект матрицы f (А) |
равен дефекту мат |
||
рицы / (J). Дефект же квазидиагональной матрицы / (/) ра |
|||
вен сумме |
дефектов |
диагональных блоков. Поэтому, обо |
значив через dj дефект блока / (У/), будем |
иметь следующее |
выражение для дефекта матрицы f (Л): |
|
d = 2 4 |
( 1. 2) |
/-=1
Вычислим дефект dj матрицы f (J/). Имеем
А=
л= i ; W r W ,
д=о
(1.3)
А =
|*=0
т
J? =
В-0
Здесь
с* = * (* - 0 • |
H + l) |
(ft= 1,2, . . . . т) |
—биномиальные коэффициенты.
Сложив равенства (1.3), умноженные соответственно на
числа |
а1П, а т _ !........а 0, получим, |
учитывая, |
что |
|
= |
0 |
|
||||||
при \i |
> |
г,* (г/ — порядок матрицы Я Г/), |
|
|
|
|
|
||||||
/№ ) = П Ч |
Ег, + |
г/ |
1 |
|
|
|
|
|
|
|
|||
Е < “ ^ Г |
Ц + |
|
+ |
|
|
|
|
||||||
|
|
|
|
|
Ц=1 |
|
+ а т _ |1с{;ХЙЯ’;.. |
(1.4) |
|
||||
|
|
|
|
|
|
|
|
|
|||||
С другой стороны, |
|
|
|
|
|
|
|
||||||
|
= |
/ (fl) (^/) = |
а0т (т — 1 ) |
(т — р. - |
j - К?“1 |
) |
м |
4 |
- |
||||
|
+ а i (m — 1) ( т |
— 2) |
(т — р) |
+ |
|
|
|
|
|||||
или |
|
|
|
|
|
|
|
4* ат—др (р — 1) |
|
П °/, |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с(Ю |
= |
K |
d |
^ |
r 4*-1 |
|
4 - |
4 - < |
w |
! M |
^ |
||
Г (Ь/) |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
(1.5) |
|
|
С учетом |
(1.5) |
равенство (1.4) |
принимает вид |
|
|
|
|
||||||
|
|
|
|
|
|
|
Г/-1 |
|
|
|
|
|
|
|
|
|
т |
- |
т |
в . ' + |
ъ |
|
|
|
|
|
|
ц—J
§ 21 |
ТЕОРЕМА |
ГАМИЛЬТОНА — КЭЛИ |
|
99 |
||||
ИЛИ |
|
|
г(ад |
|
|
|
|
|
|
|
/(а д |
|
/ < ' / - » » / ) - |
|
|||
|
|
II |
|
(/7-1)1 |
|
|
||
/ |
Ш = |
0 |
/(ад |
|
/‘Ч -21(V) |
• |
(1-6) |
|
|
(>7-2)1 |
|||||||
|
|
|||||||
|
_ |
0 |
0 |
|
|
/Я ,) |
_ |
|
Обозначим через V/ кратность |
Я/ как корня многочлена |
|||||||
/ (Я) (в частности, |
v;- может быть |
равна и нулю). Тогда |
||||||
/ (ад = |
/' (ад - |
|
= |
(ад = о, |
/^/» (ад ф о. |
|||
Поэтому дефект матрицы / (J/) (см. 1.6)) равен v/} если |
||||||||
V/ < г/, и равен г/, если |
v;- > |
/7 . Итак, |
|
|
|
|||
|
|
d/ = min(v/, |
гу), |
|
|
|
и, принимая во внимание (1 .2), приходим к следующей тео реме.
Т е о р е м а |
1.1. Дефект матричного многочлена / |
(Л), |
||
где А — матрица с элементарными делителями |
|
|||
(Я |
^i)rS (Л> |
Я2)г*, |
(^ W 5> |
|
определяется формулой |
|
|
|
|
|
d = |
£ m in (v /f |
о). |
(1.7) |
/- 1
Здесь у,-— кратность Я/ как корня / (X).
§ 2. Теорема Гамильтона — Кэли
Пусть Д (Я) — характеристический многочлен матрицы А порядка п. Кратность v;- корня Я/ характеристического многочлена не меньше, чем степень любого элементарного
делителя матрицы А вида (Я — Я/)г/. Учитывая это, соглас но (1.7) будем иметь
= 2 г,. /=i
Но сумма степеней элементарных делителей квадратной матрицы порядка п равна п. Поэтому d = п, а это означает, что А (Л) — нулевая матрица, т. е.
Д (Л) = 0.
Итак, всякая квадратная матрица удовлетворяет своему характеристическому уравнению (теорема Гамильтона — Кэли).
Эта теорема непосредственно следует и из соотношения (4.7.3). В самом деле, характеристический многочлен А (Я) оператора А в /?, которому в некотором базисе отвечает рассматриваемая матрица А , содержит в качестве множи теля минимальный многочлен всего пространства ifo (Я). По этому из равенства
'К (А) = 0
немедленно вытекает равенство
А (Л) = 0
и, значит,— равенство (2 .1 ).
§ 3. Построение матрицы, преобразующей квадратную матрицу к квазидиагональному виду
Рассмотрим квадратную матрицу U порядка п с собст венными значениями Я!, Я2, ..., Я„, среди которых могут быть и равные *). В дальнейшем эти собственные значения будем разбивать на группы, отмечая принадлежность к той
или иной группе верхним |
индексом. Так, |
Я^0) обозначает |
|||||
/-е собственное значение группы о. |
|
|
|
||||
Пусть собственные значения матрицы U разбиты на р |
|||||||
групп Я1°\ ...» ЯЦ |
(а = |
1, |
...» |
р; 2 |
= |
п)> так> что |
|
|
|
|
|
<7=1 |
|
|
|
—x f|> 0 |
( а ф т |
i |
- l |
....... А„; / = |
1, |
k s). |
|
|
|
|
|
|
|
|
(3.1) |
Каждой группе |
о поставим |
в соответствие |
многочлен |
||||
Д с ( Ч = П П |
(X — |
A f ) |
(<Т = 1 , 2 , |
|
|
р ) . ( 3 .2 ) |
|
/-=1 |
|
|
|
|
|
|
Вычислим ранг матрицы Aff (U). Для этого воспользу емся теоремой о дефекте матричного многочлена, согласно
*) Собственное значение кратности т здесь рассматривается как т равных собственных значений, и каждому из них приписывается свой индекс.