книги / Матричные и асимптотические методы в теории линейных систем
..pdfПроизвольный проекционный оператор Р в R осуществ ляет проектирование R на подпространство S — PR парал лельно подпространству Т = R — S.
Действительно, множество векторов
X s = P x |
( X $ R ) |
образует некоторое подпространство 5 пространства R. Все пространство R можно рассматривать как прямую сумму двух подпространств 5 и Т = R — S. Произвольный вектор х £ R разлагается на сумму
х = xs Н- хт |
(x s £$, |
х г € Т). |
(7.5) |
Применим к (7.5) оператор Р: |
|
|
|
Р х = P xs 4- Рхт. |
|
|
|
Но |
Р (Рх) = |
Pxs, |
|
Р х = Р2х = |
|
||
поэтому |
|
|
|
Рхт = 0 |
|
|
|
и, кроме того, |
|
|
|
Pxs = Р х = Xs. |
|
|
|
Квадратная матрица Р называется проекционной, если |
|||
Р2 = Р. |
|
(7.6) |
Проекционному оператору в произвольном базисе отве чает проекционная матрица.
Действительно, если Р — матрица, отвечающая опера
тору Р в базисе g = |
— О » то |
|
|
Р% = gP. |
|
Но, с другой |
стороны, |
|
|
/ > g = p 2g = P g P = g P 2. |
|
Значит, g P = |
g P 2, откуда и следует равенство (7.6). |
|
Докажем некоторые свойства проекционных матриц. |
||
Л е м м а |
7.1. Пусть |
Рх и Р а — две проекционные ма |
трицы. Для того чтобы матрица
Р = Р г + Р ,
также была проекционной, необходимо и достаточно, чтобы
Р1Р1= РгР, = 0. |
(7.7) |
Д о с т а т о ч н о с т ь . Пусть справедливы равенства (7.7). Тогда
(P i+ p j * = р? + р , р г + р гр , + р \ = P ? + P I = P ,+ р„
и, значит, Рг -f- Р2— проекционная матрица. |
|
|||||
Н е о б х о д и м о с т ь . |
Если |
Рг -f Рг — проекцион |
||||
ная матрица, то необходимо |
|
|
|
|
|
|
V * |
+ |
V |
i |
= |
0 - |
(7.8) |
Умножая равенство (7.8) справа на Р2, а слева на Р19 |
||||||
получим |
|
|
|
|
|
|
7 у > 2 ( Е + / у > а) = 0 . |
|
|||||
Далее, умножая равенство |
(7.8) |
справа на |
Plt а слева |
|||
на Р2, получим |
|
|
|
|
|
|
Р2Р1(£ + РйР1) = 0 . |
(7.10) |
|||||
Складывая (7.9) и (7.10) и учитывая (7.8), будем иметь |
||||||
(/V >2)2 + |
(РгРУ = 0- |
17.11) |
||||
Из двух равенств (7.8) и (7.11) находим |
|
|||||
СPiPа)2 + (РгР2)2 = 2 (PSJ* - - 0 . |
|
|||||
Отсюда |
|
|
|
|
|
|
( |
/ W |
- |
0 , |
|
|
|
а из (7.11) тогда и |
|
|
|
|
|
|
( 7 W |
= |
0 . |
|
|
||
На основании последних двух соотношений из (7.9) и |
||||||
(7.10) получаем |
|
|
|
|
|
|
^ Л |
= |
Л |
Л |
= 0 . |
|
Лемма доказана.
Проекционную матрицу Р порядка п ранга г, как и вся кую прямоугольную или квадратную матрицу, можно раз ложить на множители:
Р = КМ,
где К и М — матрицы ранга г и с размерами п X г я г х п соответственно.
Л е м м а 7.2. Пусть Р1и Р2— проекционные матрицы порядка п и рангов гг и г2 соответственно и
Pl = KlMlt Р ,= К ,М „
где Къ |
Mlt |
М2— матрицы размеров п X rlf |
rx х л, |
|||
п X г2, |
г2 X п соответственно. Тогда, |
если |
|
|||
то |
|
РХР2= РгРг= О, |
|
|
||
|
f£ ,(, |
/ = |
( , |
|
|
|
|
|
|
(7.12) |
|||
|
Af,/С/ = |
|
( / . / - 1 . 2 ). |
|||
|
|
О, |
/V=i |
|
|
|
Д о к а з а т е л ь с т в о . |
Так как |
Рг — проекционная |
||||
матрица, то |
KiM1K1M1= KXMV |
|
(7.13) |
|||
|
|
|
Отсюда следует, что ранг квадратной матрицы MxKi порядка гх равен рангу матрицы K iM lt т. е. гг. Умножим (7.13) справа на Кхи, учитывая, что МХКХ— невырожденная матрица, получим
КгМхКх= Къ
или
Ранг матрицы М1К1— Яг,» который мы обозначим че рез г', равен нулю. В самом деле, используя неравенства Сильвестра, получим
ri + r’ — гх< 0 .
Отсюда
г' = 0 ,
и, значит,
Тем же путем можно показать, что
М2Кг— ЕГ9.
Докажем теперь второе соотношение (7.12). По условию леммы имеем, например,
KXMXK*M2= 0 .
Ранг матрицы МхКгМ2 размера гг X п обозначим через г". Согласно неравенствам Сильвестра
Гж+ Г* — Га< 0 .
Значит,
г* = 0,
и, следовательно,
Обозначив |
через г ранг матрицы МХК2, будем иметь |
|||
|
г"' -1- г2 —г2< 0. |
|||
Отсюда г"' = |
0 , и потому |
|
|
|
Точно так же |
МХК2= |
0 . |
||
М2Кг = |
0 . |
|||
|
||||
Лемма доказана. |
|
|
||
Л е м м а |
7.3. Пусть Рх иР 2— проекционные матрицы |
|||
порядка п и рангов rlt г2 соответственно и |
||||
|
Р \Р ъ — Р %Р1= 0. |
|||
Тогда матрица Р = Рх -f |
Р2 является проекционной и ее |
|||
ранг равен г = гх + г2. |
То, что Р есть проекционная |
|||
Д о к а з а т е л ь с т в о . |
матрица, непосредственно следует из леммы 7.1, так что оста ется показать, что ранг матрицы Р равен rx -f г2.
Матрицы Рх и Р2 разложим на множители:
Рг = KxMlf Р2 = К2М2,
где /Ci, Мх — матрицы ранга гх и размеров п X гх и rx X п соответственно, /С2, М2— матрицы ранга г2 и размеров п X г2 и г2 X п соответственно.
Матрица
Р = Л + Р* = КХМХ+ КгМ2 = КМ, |
(7.14) |
где |
|
(МЛ |
|
K = (Kl f K2), |
|
Матрицы К и М размеров п X г' и г' х П (Г' = |
Г, + Гг) |
соответственно имеют ранг, равный Я. Действительно, учи тывая (7.12), имеем
Отсюда видно, что ранг произведения МК равен г'. Но ранг сомножителей не меньше, чем ранг произведения, и, поскольку ранг матриц М и К не может быть и больше г', их ранг в точности равен г'.
Остается к произведению КМ применить неравенства Сильвестра:
К -f г' — г <^г
Отсюда
г = г' = г1 + г2.
Лемма доказана. |
|
|
|
Рр — проекционные |
||||
С л е д с т в и е . |
Пусть Ръ Р2> ..., |
|||||||
матрицы рангов гъ |
г2, ..., гр |
соответственно и |
|
|||||
|
Р(Р}= 0 |
|
|
i, j — 1,2, |
. . . , p). |
|
||
|
p |
|
проекционной |
матрицей |
ранга |
|||
Тогда P = 2 Л является |
||||||||
|
/=i |
|
|
|
|
|
|
|
г = % . |
|
|
|
|
|
|
|
|
/= 1 |
|
Пусть |
Р1У Р2, |
|
Рр — проекционные |
|||
Л е м м а 7.4. |
..., |
|||||||
матрицы порядка nt рангов гъ г2, |
..., гр |
соответственно и |
||||||
такие, что |
|
|
|
|
|
|
|
|
а) |
PtP/ = 0 |
|
|
* , / = 1 , 2 , |
. . . . /?), |
|
||
б) |
г\-\~ Г2 |
|
4 " г0= п. |
|
||||
Тогда |
Р = РХ+ Р*+ |
|
+Рр = Еп. |
|
||||
|
|
|
||||||
Д о к а з а т е л ь с т в о . |
Матрица |
Р является |
проек |
|||||
ционной |
матрицей ранга п. Разложим |
Р на множители: |
||||||
|
|
Р = КМ. |
|
|
|
(7.15) |
Здесь К, М — невырожденные квадратные матрицы поряд ка п. Умножая равенство (7.15) справа на К и учитывая, что
МК = Еа
(см. (7.12)), |
получим |
|
{ Р - Е п)К = 0 . |
Отсюда, так как К — невырожденная матрица, |
|
|
Р = Еп. |
Лемма доказана. |
|
Л е м м а |
7.5, Пусть Р^, Р2* •••> Рр — проекционные |
матрицы, удовлетворяющие соотношениям |
в) 4~ Дi 4" 4~ Р р ~ Я».
и А — некоторая квадратная матрица такая, что
б) |
PLA=*0 |
( / = 1 , 2 ...........р). |
Тогда А = |
0. |
|
Д о к а з а т е л ь с т в о |
следует из цепочки равенств: |
|||||||
|
А = ЕА — (Рi -{- Р2+ |
|
+ Р р) Л = 0. |
|||||
Л е м м а |
7.6. |
Пусть Р1У Я2, |
..., |
Рр — квадратные |
||||
матрицы порядка п, удовлетворяющие равенствам |
||||||||
а) |
P i P , = о |
( i , / = |
1 |
, 2 , |
р; |
гф /), |
||
б) |
|
Pi + Pt + |
|
+Р„ = £«• |
||||
Тогда Plt Я2, |
.... Рр — проекционные матрицы. |
|||||||
Д о к а з а т е л ь с т в о . |
Умножая равенство б) на р» |
|||||||
получаем |
|
|
|
|
|
|
|
|
что доказывает лемму. |
Р\ = |
Р{, |
|
|
||||
|
|
Рг* |
|
Рр — операторы, |
||||
Л е м м а |
7.7. |
Пусть Рц |
—I |
|||||
удовлетворяющие условиям |
|
|
|
|
||||
а) |
|
|
PtPi = 0 |
(i¥=j), |
|
|||
б) |
|
|
P i + |
P * + |
|
+ Рр= Е, |
где 0 и Е — соответственно нулевой и единичный операторы. Тогда P lt Р 2, ..., Рр— проекционные операторы, расщепля ющие пространство R на подпространства R ly /?2, ..., R p (Rt = PiR)t m. e.
Д о к а з а т е л ь с т в о . Так же как и в предыдущей лемме, по умножении равенства б) на P t убеждаемся в том, что операторы Р { (i = 1, 2, ..., р) — проекционные.
Далее, из условия б) получаем
х = |
* i -Ь х г + |
+ х р |
(xt = Ptx). |
|
Если х £ |
R {, то или х |
£ /?, (i Ф /), |
или х = 0 . В са |
|
мом деле, допуская, что.* £ R t и х £ Rj, |
будем иметь Я,.* — |
|||
Я/.*. Но тогда Я?* = Р (х |
= Я*Я/л; = |
0 |
и, значит, х = 0 . |
|
Лемма доказана. |
|
|
|
Г л а в а IV
РАСЩЕПЛЕНИЕ ПРОСТРАНСТВА НА ИНВАРИАНТНЫЕ ПОДПРОСТРАНСТВА.
НОРМАЛЬНЫЕ ФОРМЫ МАТРИЦЫ
§ 1. Минимальные многочлены вектора, векторного пространства, матрицы
Рассмотрим я-мерное векторное пространство R над чис ловым полем ЗС и линейные операторы в нем. Если А —
линейный оператор в/?, то АЛ |
= Л а, ААА = Л 3, А |
...А — Ат — тоже линейные |
операторы в R. Будем счи |
тать, что нулевая степень любого линейного оператора в
R равна единичному оператору |
Е\ А0 = Е. Многочлен |
|
f (^) = а(Лт 4“ |
4“ |
4- ат—1^4 “ &т, |
где at £ ЗС, после замены |
в нем скалярного аргумента Я |
линейным оператором А тоже представляет собой некоторый линейный оператор, а именно:
f(A) = а0Ат4- «1Л'П- 1 4 - 4 - a,„_iА 4- а „ Д
З а м е ч а н и е . Пусть / (Я) и g (Я) — два многочлена относительно скалярного аргумента с коэффициентами из поля 3i. Тогда в силу того, что для любых целых неотрица тельных р и q (а в случае невырожденного оператора для любых целых р и q)
Л М * = Л Р+* = А М Р,
имеем
f(A)g(A)=g(A)f(A),
т. е. любые два многочлена с коэффициентами из поля 3£ относительно одного и того же линейного оператора ком мутативны.
Многочлен скалярного аргумента А, со старшим коэффициентом, равным единице,
Ф (^) = ^ + ctjA?- 1 -f- -f- «21 • ■• *
называется аннулирующим многочленом вектора х , если
ср (А)х = 0 .
Аннулирующий многочлен вектора х наименьшей степе ни называется минимальным аннулирующим многочленом векторах или просто минимальным многочленом вектора х .
Допустим, что степень минимального многочлена век
тора х равна р, т. е. ф (А) = |
Ар -f с^Ар- 1 4- |
4- ap_iA 4- |
4 - ар. Тогда |
|
|
А рх 4 - « 1А р~Ах 4- |
4- a P- i А х 4- %Х — 0 . |
|
Отсюда вытекает, что векторы |
|
|
х , Ах, |
... , А рх |
|
линейно зависимы. В то же время векторы |
|
|
х , А х, |
А р~'х |
|
линейно независимы, так как нет многочлена степени р — 1 ,
который был бы для вектора х аннулирующим. |
р, |
|
|
Если степень минимального |
многочлена равна |
а |
|
ф (А,) — многочлен степени меньше, чем р, то из |
|
|
|
ф(Л)д? = |
0 |
(1 . |
1 ) |
следует |
|
|
|
ф(А)=гО.
Действительно, если допустить, что ф (А,) Ф 0 , то соглас но (1.1) ф (А,) — аннулирующий многочлен для х , что невоз можно, так как степень минимального аннулирующего мно гочлена вектора больше, чем степень многочлена ф (А,).
Каждому вектору х отвечает только один минимальный многочлен. В самом деле, пусть
Ф(А») = АР 4~ о^Ар- 14~ |
4~ otp |
|
и |
|
|
ф (^) = А.р 4" |
1 И- |
4~ Рр |
— два минимальных многочлена вектора х , т. е.
у ( А ) х — 0 , |
— 0 . |
S 1] |
МИНИМАЛЬНЫЕ МНОГОЧЛЕНЫ ВЕКТОРА |
69 |
Тогда |
|
|
|
[<р (Л) — ф (Л)] л: = 0 . |
(1.2) |
Но степень многочлена ф (А) — <р (А,) во всяком случае меньше, чем р, поэтому в силу равенства (1 .2 )
ф W ~ <РW = 0,
и, значит,
=(/ = 1 , 2 , . . . , р).
Пусть ф (А,) — произвольный аннулирующий многочлен вектора х , а ф (А)— его минимальный многочлен. Тогда
Ф (А,) нацело делится на ф (А,).. Действительно, степень ф (А,)
г*
не меньше, чем степень ф (А,). Разделив ф (А.) на ф (А,), полу чим
ф (А.) = Ф (А,) и (К) + г (А,), |
(1.3) |
где г (А) — остаток от деления. |
|
Из (1.3) находим |
|
г(А )х = 0. |
|
Отсюда |
|
г (I) а 0 , |
|
так как степень многочлена г (А,) меньше, |
чем степень ми |
нимального многочлена ф (А,).
Многочлен ф (А,), который является, аннулирующим для любого, вектора х из R, называется аннулирующим много членом пространства R.
Пусть g = (e\ ez ...en) — базис пространства R , а ф (А,) — аннулирующий многочлен пространства R:
ф(Л)лг = |
0 |
(V *€#?). |
(1.4) |
В силу (1.4) |
|
|
|
ф(Л)е, = 0 |
( / = |
1,2, . . . , п). |
(1.5) |
Пусть ф! (А,), ...» Ф„ (Я) — минимальные многочлены век- |
|
торов el t ...»еп. Из |
(1.5) следует, что ф (А,) делится без остат |
ка на наименьшее |
общее кратное многочленов q^, ..., ф„. |
Это наименьшее общее кратное обозначим через ф (А.). Тогда |
ф (Л ) * - 0 (<— 1 , 2 ............я) (1 .6)
и, следовательно, |
|
\ р ( А ) х = ур(А)%х = (ф(Л)<?! |
ф ( Л ) е л)х = 0. |
Таким образом, ф (Я) является аннулирующим многочле ном пространства R . Из всех многочленов с равными едини це старшими коэффициентами, удовлетворяющих соотноше ниям (1.6), ф (Я) является многочленом наименьшей степени. Этот многочлен называется минимальным многочленом про странства R.
Заданием линейного оператора минимальный многочлен пространства определяется единственным образом. Из един ственности минимального многочлена всего пространства следует, что он не зависит от выбора базиса.
Минимальный многочлен пространства /?, являясь ан нулирующим многочленом для любого вектора из /?, де лится на минимальный многочлен любого вектора х £ R без
остатка. |
(ех е 2... е п) — базис, |
|
Пусть g = |
а |
|
ф (Я) |
Яр ~|—с^Яр ^ “J- |
“4” оср_|Я -j- otл |
— аннулирующий многочлен пространства /?, т. е. для лю |
||
бого х £ R |
ф ( Л ) л ; = 0 . |
(1.7) |
|
||
Имеем (см. (3 .1 .2 )) |
|
|
|
Л 8 = & 4 , |
(1 .8) |
где А — квадратная |
матрица, |
отвечающая оператору А в |
||||
базисе g. |
|
|
|
|
|
|
|
Из (1 .8) находим |
|
|
|
|
|
и, |
Л *8 = Л g Л = g Л a, |
|
||||
вообще, |
Л |
|
%АК |
|
||
|
Вследствие этого |
|
|
|||
|
*(А)8 = 8 |
$ (А), |
|
|||
|
|
|
||||
и из равенства (1.7) |
получаем |
|
|
|
||
|
|
gi|) (А) * = |
0. |
|
||
Отсюда, так как векторы |
е г..... е п линейно независимы, |
|||||
|
|
ф (Л) X = |
0 , |
|
||
и, |
поскольку х — произвольный |
вектор |
из R , |
|||
|
|
ф(Л) -= 0. |
(1.9) |