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

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

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

Произвольный проекционный оператор Р в 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),

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

Многочлен скалярного аргумента А, со старшим коэффициентом, равным единице,

Ф (^) = ^ + 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)

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