![](/user_photo/_userpic.png)
книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем
.pdf6 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) = 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)