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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
22
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

6 0

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

[ГЛ. III

Тогда любой вектор х £ R разлагается, и притом един­ ственным образом (см. § 6), на сумму двух векторов из 5 и Т:

* = .* :s + .vr

(ATS £ S, Хт £ T).

(7.1)

Вектор Xs называется проекцией вектора х на подпрост­ ранство S параллельно подпространству Т. Аналогично, Хт называется проекцией вектора х на подпространство Т параллельно подпространству S.

Пусть Р — оператор, осуществляющий проектирование пространства R параллельно подпространству Т. Этот опе­ ратор определяется равенством

Р х = xs,

где X — произвольный вектор из R, a x s его проекция на S. Очевидно, что Р является линейным оператором.

Равенство (7.1) можно записать и так:

X = Р х -+• Хт-

(7.2)

Если л: £ S, то разложение (7.1), в силу единственности разложения вектора из R на сумму (7.1), принимает вид

.V = Xs {Хт = 0),

и из (7.2) в этом случае получаем

x s = Pxs,

т. е. оператор Р, примененный к вектору из S, действует как единичный оператор.

Если X £ Т, то разложение (7.2) принимает вид

Хт = Р хг + Хт,

и, значит,

Рхт = 0.

(7.3)

 

Пусть теперь х — произвольный вектор из R. Применим

оператор Р

к обеим частям равенства

(7.2). Будем иметь

 

Р х = Р2Х + Рхт.

 

Отсюда,

учитывая (7.3), получим

 

 

Р х = Р'1Х.

 

Следовательно,

 

 

Р2 = Р.

(7.4)

Оператор Р в R, удовлетворяющий равенству (7.4),

называется

проекционным оператором.

 

§ 7]

П Р О Е К Ц И О Н Н Ы Е О П Е Р А Т О Р Ы И М А Т Р И Ц Ы

61

Произвольный проекционный оператор Р в R осуществ­ ляет проектирование R на подпространство S’ = PR парал­ лельно подпространству Т = R S.

Действительно, множество векторов

Xs = Р х

(х£ R)

образует некоторое подпространство 5 пространства R. Все пространство R можно рассматривать как прямую сумму двух подпространств S и Т = R — 5. Произвольный вектор X £ R разлагается на сумму

 

X = Xs + Хт

(Xs £S, Xт£ Т).

(7.5)

Применим к (7.5) оператор Р:

 

 

Но

 

Р х = P xs + Рхт.

 

 

 

Р х = Р2х — Р (Рх) = Pxs,

 

поэтому

 

Р х,

= 0

 

 

и, кроме того,

 

 

P xs = Р х = xs-

 

 

 

 

 

 

Квадратная

матрица Р называется

проекционной, если

 

 

Р 2 = Р.

 

(7.6)

Проекционному оператору в произвольном базисе отве­

чает проекционная матрица.

 

отвечающая

опера­

Действительно, если Р — матрица,

тору

Р в базисе g = (ехе2... еп), то

 

 

 

 

Р% =

gP.

 

 

Но,

с другой

стороны,

 

 

 

 

 

P g = p 2 g = / > g P = g P 2 .

 

Значит, gP = g P 2, откуда и следует равенство (7.6). Докажем некоторые свойства проекционных матриц. Л е м м а 7.1. Пусть Р1 и Р2 — две проекционные ма­

трицы. Для того чтобы матрица

Р = Р, + Р2

также была проекционной, необходимо и достаточно, чтобы

Р 1Р 2 = Р 2Р 1 = 0.

(7.7)

6 2 Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы [ГЛ . Ш

Д о с т а т о ч н о с т ь .

Пусть справедливы

равенства

(7.7). Тогда

 

 

 

 

 

 

 

( P i +

P J * = Р 2І +

P S * +

P S i +

P i =

P l + P i = P l + P 2l

и, значит, Px

P%— проекционная

матрица.

 

Н е о б х о д и м о с т ь .

Если

Рг + Р2 — проекцион­

ная матрица, то необходимо

 

 

 

 

 

 

 

Р ^

+

Р .Р і^ О .

 

(7.8)

Умножая равенство (7.8) справа на Р2, а слева на Рх,

получим

Р 1Р 2(Р +

Р 1Р 2) =

0.

(7.9)

 

 

Далее, умножая равенство (7.8) справа на Р1( а слева

на Р 2, получим

Р 2Р 1(£ +

Р 2Р 1) =

0.

(7.10)

 

 

Складывая (7.9) и (7.10) и учитывая (7.8), будем иметь

 

 

х^2)2 +

( ^

і)2 =

0.

(7.11)

Из двух равенств (7.8) и (7.11) находим

 

 

л.Рг)2 + (РіРг)2 =

2 (РХР2)2 — 0.

 

Отсюда

 

(РіР*)2=о,

 

 

 

а из (7.11) тогда и

 

 

 

(Р2Рі)2 =

0.

 

 

 

 

 

 

 

 

 

На основании последних двух соотношений из (7.9) и

(7.10)

получаем

 

 

 

 

 

 

 

 

 

Р1Р 2 =

Р 2Р1 = 0.

 

 

Лемма доказана.

Проекционную матрицу Р порядка п ранга г, как и вся­ кую прямоугольную или квадратную матрицу, можно раз­ ложить на множители:

Р = КМ,

где К и М — матрицы ранга г и с размерами п X г и г х п соответственно.

Л е м м а 7.2. Пусть Рх и Р2 — проекционные матрицы порядка п и рангов гх и г2 соответственно и

P i =

P s = K aMt,

5 7]

 

 

П Р О Е К Ц И О Н Н Ы Е О П Е Р А Т О Р Ы И М А Т Р И Ц Ы

6 3

где

Ки

Мъ

К2, М2 матрицы размеров п X гх,

гх х п,

п X

г,,

г„

X

соответственно. Тогда,

если

 

то

 

 

 

Р іР 2 =

1=

О

 

 

 

 

 

[ЕГ

І =

і,

 

 

 

 

 

 

 

(W

= 1,2).

(7.12)

 

 

 

M‘Kl = { o t

j¥si

 

Д о к а з а т е л ь с т в о .

Так

как Рх — проекционная

матрица,

то

К1М1К1М1 = К1М1.

 

(7.13)

 

 

 

 

 

Отсюда следует, что ранг квадратной матрицы МгКі порядка гх равен рангу матрицы К1М1, т. е. гх. Умножим (7.13) справа на /Сх и, учитывая, что М^Кх — невырожденная матрица, получим

К х М х К ^ К х ,

или

К Л М х К х - Е г ^ 0.

Ранг матрицы МгКх ЕГі, который мы обозначим че­ рез г', равен нулю. В самом деле, используя неравенства Сильвестра, получим

гі + г' — Г і < 0 .

Отсюда

г' = 0,

и, значит,

М х К х ^ Е г , .

Тем же путем можно показать, что

М2К 2 = Ег%.

Докажем теперь второе соотношение (7.12). По условию леммы имеем, например,

КхМхК2М2 = 0.

Ранг матрицы МуКгМ2 размера rx X п обозначим через г". Согласно неравенствам Сильвестра

гх + г" — гх < 0.

Значит,

г* = 0,

и, следовательно,

М х К щ М ^ О .

6 4

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

[ГЛ . I l l

Обозначив через r'" ранг матрицы МХК*Убудем иметь

г ' + г2 — г2-^0. Отсюда г " ' 0, и потому

МХК 2 = 0.

Точно так же

М2КХ= 0.

Лемма доказана.

Л е м м а 7.3. Пусть Рх и Р2 проекционные матрицы порядка п и рангов r x , г 2 соответственно и

^

=

^ 1 = 0.

Тогда матрица Р =

Рх -f

Р2 является проекционной и ее

ранг равен г = rx -f г2.

 

То, что Р есть проекционная

Д о к а з а т е л ь с т в о .

матрица, непосредственно следует из леммы 7.1, так что оста­ ется показать, что ранг матрицы Р равен r x -f r 2 .

Матрицы Рх и Р2 разложим на множители:

 

 

Рх= КХМХ,

Р2 - К2М2,

 

 

 

где К.1 , Мх — матрицы ранга г х

и размеров п X

гх и г х X п

соответственно, К2>М2 — матрицы ранга

г2

и

размеров

п X г 2 и r 2 X п соответственно.

 

 

 

Матрица

 

 

 

 

Р = Рі + Р 2 = К ,м х + К М 2 = к м ,

(7. ] 4)

где

 

 

 

 

 

/ М л

 

 

 

 

м = ( м ;) .

 

 

 

Матрицы К и М размеров п

X г' и г' X

п (r' =

гх + г2)

соответственно имеют ранг, равный г'. Действительно, учи­ тывая (7.12), имеем

(МЛ

(Er, 0 \

М К = [ м Т К ^ =

 

Отсюда видно, что ранг

произведения МК равен г’.

Но ранг сомножителей не меньше, чем ранг произведения, и, поскольку ранг матриц М и К не может быть и больше г', их ранг в точности равен г'.

Остается к произведению КМ применить неравенства Сильвестра:

г < г < г ', г '.

ST)

ПРОЕКЦИОННЫ Е ОПЕРАТОРЫ И МАТРИЦЫ

65

Отсюда

 

 

 

 

 

 

 

 

 

г = г' =

гх 4- г%.

 

 

Лемма

доказана.

Пусть Plt Р2, ..., Рр проекционные

С л е д с т в и е .

матрицы рангов rlt г2, ..., гр соответственно и

 

 

Р(Р1=

0

(і ф /;

£ ',/= 1 ,2 ,

. . . , р).

 

Тогда P

Р

является проекционной

матрицей

ранга

S P ,

/=і

р

г =

/=1

Л е м м а 7.4. Пусть Ри Р2, ..., Лр — проекционные матрицы порядка п, рангов rL, г2, ..., гр соответственно и такие, что

а)

=

0

(£ =Л /;

£ ',/= 1 ,2 .......... р),

б)

 

0 +

^ 2+

•••

+Гь = п.

Тогда

 

 

 

 

-

 

=

-^і +

^2 +

•••

+

= Дг-

Д о к а з а т е л ь с т в о . Матрица Л является проек­ ционной матрицей ранга п. Разложим Р на множители:

Р = КМ.

(7.15)

Здесь К, М — невырожденные квадратные матрицы поряд­ ка п. Умножая равенство (7.15) справа на К и учитывая, что

МК — Еп

(см. (7.12)), получим

( Р - Е п) К = 0.

Отсюда, так как К — невырожденная матрица,

Р = Еа.

Лемма доказана.

Л е м м а 7-5. Пусть Ри Л2, ..., Рр проекционные матрицы, удовлетворяющие соотношениям

а)

^ 1 + ^ 2 + ••• + -Р р = Дп,

и А некоторая квадратная матрица такая, что

б)

Л/Л = 0

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

р).

Тогда А =

0.

 

 

8 К, А. Абгіряі

6 6

 

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

 

[ГЛ III

Д о к а з а т е л ь с т в о

следует из цепочки

равенств:

 

А =

ЕА =

(Pj +

Р 2 +

•••

+ -Р р)^ = 0-

 

Л е м м а

7.6.

Пусть

Ри

Р2,

...,

Рр квадратные

матрицы порядка п, удовлетворяющие равенствам

а)

Р,-Р/ = 0

(t, /

=

1,2..........р;

іфі),

 

б)

 

-Р1 —}- Z32 —[— - - - —

= Еп.

 

Тогда Рх, Р2, ...,

Рр проекционные матрицы.

Д о к а з а т е л ь с т в о .

Умножая равенство

б) на Р,,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

Р* =

Pt,

 

 

 

что доказывает лемму.

 

 

 

 

Р р — операторы,

Л е м м а

7.7.

Пусть

Р х,

Р 2, ....

удовлетворяющие условиям

 

 

 

 

 

 

а)

 

 

P {P j =

0

(/ ф /),

 

 

б)

 

Рі + Р г+

 

+ Р р = Е ,

 

где 0 и Е соответственно нулевой и единичный операторы. Тогда Р г, Р 2, Ррпроекционные операторы, расщепля­ ющие пространство R на подпространства R lt R 2, ..., R p (Ri = P (R ), m. e.

R = R , + R 2+ ... + R p.

Д о к а з а т е л ь с т в о . Так же как и в предыдущей лемме, по умножении равенства б) на Р,- убеждаемся в том, что операторы Р с = 1, 2, ..., р) — проекционные.

Далее, из условия б) получаем

х = х г + х 2 +

■■■ + х р

(х, = РіХ).

 

Если X £ Ri, то или х £ R, (і ф

/), или х — 0. В

са­

мом деле, допуская, что х

£ R t и х £ Rj, будем иметь P tx

=

PjX. Но тогда Р\х = Р іХ = P tPjX = 0 и, значит, л: =

0.

Лемма доказана.

 

 

 

Г л а в а IV

РАСЩЕПЛЕНИЕ ПРОСТРАНСТВА НА ИНВАРИАНТНЫЕ ПОДПРОСТРАНСТВА.

НОРМАЛЬНЫЕ ФОРМЫ МАТРИЦЫ

§ 1. Минимальные многочлены вектора, векторного пространства, матрицы

Рассмотрим п-мерное векторное пространство R над чис­

ловым полем Ж и линейные операторы в нем. Если

А

линейный оператор в/?, то А А = А 2, ААА = А 3,

А ...

...А = А т — тоже линейные операторы в /?.

Будем счи­

тать, что нулевая степень любого линейного

оператора в

R равна единичному

оператору Е: А° = Е.

Многочлен

f W = а іДт +

а 1Я,т ~1-)- • • • + сст —1к

ат,

где ^ Ж, после замены в нем скалярного аргумента к линейным оператором А тоже представляет собой некоторый линейный оператор, а именно:

/ (Л) = ct0Am +

о^А'"-1 -}-

• • •

+ ат_і А -j- &тЕ.

З а м е ч а н и е .

Пусть f (к)

и g

(к) — два многочлена

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

A pA q= Ap+q = A qA p,

имеем

f (А) g (А) = g (А) f (А),

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

»*

6 8

Р А С Щ Е П Л Е Н И Е П Р О С Т Р А Н С Т В А

[ГЛ . IV

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

ф (X) = + a fo -' + • • • + aq_ xX + aq, alt а 2, . .. , aq£ ді,

называется аннулирующим многочленом вектора х , если (р(А)х = 0.

Аннулирующий многочлен вектора х наименьшей степе­ ни называется минимальным аннулирующим многочленом векторах или просто минимальным многочленом вектора х .

Допустим, что степень минимального многочлена век­ тора X равна р, т. е. ф (А,) = Хр + а^Кр~1-f- ... + арХХ

+ ар. Тогда

А рх-{- <х1А р~'1х +

■■■ + ар_хА х

арх = 0.

Отсюда вытекает, что векторы

 

 

X ,

Ах, . ..

, А рх

 

линейно зависимы. В то же время векторы

 

X , Ах, . .. ,

А р~'х

 

линейно независимы, так как нет многочлена степени р — 1,

который был бы для вектора х аннулирующим.

р, а

Если степень минимального

многочлена равна

ф (X) — многочлен степени меньше, чем р, то из

 

ф (Л) X — 0

(1.1)

следует

0.

 

ф (X) =

 

Действительно, если допустить, что ф (X) =j=0, то соглас­ но (1.1) ф (X) — аннулирующий многочлен для х, что невоз­ можно, так как степень минимального аннулирующего мно­ гочлена вектора больше, чем степень многочлена г|э (Ä,).

Каждому вектору х отвечает только один минимальный многочлен. В самом деле, пусть

ф (X) = Хр -)- а]Хр~1-f- • • • +

и

ф (?*,) = Хр

■ + ßp

— два минимальных многочлена вектора х, т. е.

ф (А )**=0, ф (А)а; = 0.

$ П

Тогда

М И Н И М А Л Ь Н Ы Е М Н О Г О Ч Л Е Н Ы ВЕКТОРА

69

 

[ф (і4 ) — < р (Л )]л : = 0.

(1 .2 )

Но степень многочлена ф (к) — ф (к)

во всяком случае

меньше, чем р,

поэтому в силу равенства (1.2)

и, значит,

Ф (к) — ф (к) = О,

 

a. = ß<

(* = 1,2, . . . , р)•

 

Пусть ф (к) — произвольный аннулирующий многочлен

вектора х, а

ф ( \) — его

минимальный

многочлен. Тогда

Ф (к) нацело делится на ф (к). Действительно, степень ф (к)

не меньше, чем степень ф (к). Разделив ф (к) на ф (?і), полу­ чим

Ф ( к ) = ф W х ( к ) + г (к),

(1.3)

где т(Ji.) — остаток от деления. Из (1.3) находим

г (А ) х = 0.

Отсюда

г ( к ) = 0,

так как степень многочлена г (к) меньше, чем степень ми­ нимального многочлена ф (к).

Многочлен ф (к), который является аннулирующим для любого вектора х из R , называется аннулирующим много­ членом пространства R.

Пусть g = (е1 е%... еп) базис пространства R , аф (к) — аннулирующий многочлен пространства /?:

Ф (А)х = 0

(Y x £ R ) .

(1.4)

В силу (1.4)

 

 

ф(Л)бі = 0

(t= 1,2, . . . , п).

(1.5)

Пусть фг (X), ..., ф„ (к) — минимальные многочлены век­

торов ег...... еп. Из (1.5) следует, что ф (к) делится без остат­ ка на наименьшее общее кратное многочленов ф!, ..., ф„. Это наименьшее общее кратное обозначим через ф (к). Тогда

Ф (Л)е( = 0 ( £ = 1 ,2 .......... я) (1.6)

Соседние файлы в папке книги из ГПНТБ