книги / Линейная алгебра
..pdfв характеристическом многочлене) называют сп ектром линейного оператора. Если матрицу отождествляют с оператором, то множе ство всех ее собственных значений называют сп ектром матрицы .
Из сказанного следует, что длл отыскания всех собственных зна чений оператора с матрицей А нужно найти все характеристиче ские числа матрицы А и из них выбрать лишь те, которые принад лежат основному полю, а длл отыскания всех собственных векто ров оператора с матрицей А нужно найти все ненулевые решения системы (3.21) при каждом собственном значении А. Это правило поясним на примере.
П ример. Для оператора с матрицей |
|
|
|||
|
/ з |
1 |
- 3 |
\ |
|
А = |
3 |
1 |
- 1 |
, |
0 / |
|
V |
2 |
|
- 2 |
действующего в действительном пространстве, найти собственные значения и собственные векторы.
Реш ение. Характеристический многочлен
3 |
- А |
1 |
- 3 |
= —А(3 - |
А)(1 - А) + |
|А - ХЕ\ = |
3 |
1 - А |
- 1 |
||
|
2 |
- 2 |
-А |
|
|
+18 + 6(1 - |
А) - |
2(3 - |
А) + ЗА = |
(А2 + 4)(4 - А) |
матрицы А имеет корни Ai = 4, А2,з = ± 2 i. Так как рассматриваемый оператор действует в действительном линейном пространстве, то его собственным значением будет лишь А = 4. При этом значении А система (А — АЕ)Х = 0 имеет вид
“ *1 |
+ |
Х2 |
-- |
Зхз |
= |
0, |
3xi |
— |
3x2 |
- |
*3 |
= |
0, |
2xi |
— |
2X2 |
-- |
4хз |
= |
0. |
Ее общим решением является X = (xi, xi, 0) с произвольным посто янным x i. При xi, пробегающем все действительные значения, оно дает общий вид собственных векторов оператора с матрицей А, при надлежащих собственному значению А = 4. Других действительных собственных векторов оператор с матрицей А не имеет, так как у него нет других собственных значений.
Собственные векторы линейного оператора ip с матрицей А, при надлежащие одному и тому же собственному значению А, вместе с ну левым вектором образуют подпространство, которое называют со б ственны м п одп ростран ством оп ератора <р по А. Это подпро странство является ядром оператора р — Хе с матрицей А — АЕ.
Размерность собственного подпространства оператора (р по соб ственному значению А равна п—г\, где п - размерность пространства, в котором действует оператор у?, т.е. порядок матрицы А , г\ - ранг оператора р — Хе или, что то же самое, ранг матрицы А —ХЕ. Эту размерность называют геом етри ческой к р а тн остью собствен н ого значения А. Другими словами, геометрической кратностью собственного значения А называют максимальное число линейно не зависимых собственных векторов оператора р, принадлежащих соб ственному значению А. Геометрическая кратность собственного зна чения не превосходит его алгебраической кратности, т.е. кратности, с которой А входит корнем в характеристический многочлен \А —\Е\.
Собственные векторы линейного оператора р , принадлежащие различным собственным значениям, линейно независимы. Их линей ные комбинации, вообще говоря, не являются собственными векто рами линейного оператора (р.
Если отождествлять оператор с его матрицей, то естественно го ворить о собственных значениях и собственных векторах матрицы. На практике так обычно и поступают. Квадратную матрицу назы вают п ростой , если для каждого собственного значения матрицы его геометрическая кратность совпадает с алгебраической кратностью. В противном случае матрицу называют деф ектной .
Собственный вектор х, определенный условием Ах = Ах, еще на зывают правым собственны м век тор ом м атрицы А , принадле жащим собственному значению А. Рассматривают и левые собствен ные векторы матрицы. Ненулевой вектор у называют левым со б ственны м вектор ом м атрицы А , принадлежащим собственному
значению А, если
у'А = Ау*
Напомним, что здесь звездочка означает транспонирование и ком плексное сопряжение. Если это равенство подвергнуть транспони рованию, то придем к равенству Ату = Ау. Если провести еще и комплексное сопряжение, то получим А*у = Ху. Это означает, что левый собственный вектор у матрицы А, принадлежащий собствен ному значению А, является правым собственным вектором матрицы А* по А, а вектор у - правым собственным вектором матрицы Атпо
А. В дальнейшем, говоря о собственных векторах, без указаний ” пра вый” , ” левый” , будем понимать, что речь идет о правых собственных векторах.
3.8.Линейные операторы п р остои стр у к ту р ы
Линейный оператор, действующий в линейном пространстве Х п, называют оп ератором п р остои структуры , если в пространстве Х п существует базис, состоящий из собственных векторов этого опе ратора. В силу п. 3.2 в базисе из собственных векторов матрица оператора простой структуры имеет вид
А1 |
0 ^ |
|
(3.22) |
Ап /
где Ai, А2, ..., Ап - собственные значения оператора.
Если в исходном базисе е оператор простой структуры имеет ма
трицу А ) а в базисе е' из собственных векторов - |
матрицу Л, то в |
силу соотношения (3.13) из п.3.4 имеем |
|
А = Т~1АТ1 |
(3.23) |
где Т - матрица перехода от базиса е к базису е' Она состоит из столбцов координат векторов базиса е! в базисе е , Л - матрица вида (3.22). На матричном языке соотношение (3.23) означает, что ма трица А приводится матрицей Т к диагональному виду.
Разрешив соотношение (3.23) относительно матрицы А, придем к соотношению
А = ТАТ~\ |
(3.24) |
которое называют каноническим разложением матрицы А.
П р и п о с т р о е н и и м а т р и ц ы Т д л я с о о т н о ш е н и й (3 .2 3 ) и (3 .2 4 ) н у ж н о н а й т и в се с о б с т в е н н ы е з н а ч е н и я м а т р и ц ы А и п р и к а ж д о м с о б с т в е н н о м з н а ч е н и и А* п о с т р о и т ь Ф С Р о д н о р о д н о й с и с т е м ы у р а в н е н и й ( А - Х [Е )Х = 0; и з р е ш е н и й в с е х п о с т р о е н н ы х Ф С Р , к а к и з с т о л б ц о в , с о с т а в и т ь м а т р и ц у Т . П р и ч е м в м а т р и ц е Т ст о л б ц а м и з а п и с ы в а ю т с я р е ш е н и я п о к а ж д о м у А* в п о р я д к е н у м е р а ц и и с о б с т в е н н ы х з н а ч е н и й А], Аг, ., Ап ( о д и н а к о в ы е Ai с ч и т а ю т с я с т о л ь к о
раз, каковы их крайности; все А,- можно занумеровать так, что будет Xi, Х2, ..., Аг ф 0, \r+i = = Лп = 0). Если матрица Т окажется квадратной, то она будет удовлетворять соотноше ниям (3.28) и (3.24)- Если же матрица Т окажется неквадратной, то соотношения (8.28) и (8.24) <?лл матрицы А будут невозможны, т.е. матрица А не приводится к диагональному виду и, следова тельно, не имеет канонического разложения. Этот способ равноси лен нахождению невырожденной матрицы Г из матричного уравне ния АТ = ТА. Из правила построения матрицы Т видно, что она будет квадратной лишь в случае, когда каждое характеристическое число А,- матрицы А является ее собственным значением и для каждого А,- совпадает его алгебраическая кратность с геометрической кратностью, т.е. с максимальным числом линейно независимых соб ственных векторов матрицы А по А,-, равным /,• = п — г,-, где г,- - ранг матрицы А — А{Е. Лишь в таком случае оператор с матрицей А является оператором простой структуры, а матрица А приводится к диагональному виду.
Операторами простой структуры, в частности, являются опе раторы простого спектра, т.е. операторы, имеющие п различных собственных значений, и симметрические операторы (см. п. 6.7), а также нормальные операторы (см. п. 6.8).
П ример 1 . Привести, если возможно, действительные матрицы
/ - 1 |
3 - 1 \ |
/ 0 |
3 |
3 \ |
||
Аг = —3 5 - 1 |
, |
А2 = |
- 1 |
8 |
6 |
|
V - 3 |
3 1 |
/ |
V |
2 |
—14 |
—10 / |
к диагональному виду и построить для них канонические разложения.
Реш ение. а) Корнями характеристического многочлена
1 — А |
3 |
- 1 |
|
\А\ — ХЕ\ = - 3 |
5 - А |
- 1 |
= _ ( А - 2 ) 2( А - 1) |
- 3 |
3 |
1 - |
А |
матрицы А\ являются числа Ai = 2, Aj = 1 соответственно крат ности ki = 2, &2 = 1. Все они действительные и потому являются собственными значениями матрицы А\. Причем при А* = 2 матрица
- 3 |
3 |
- 1 |
\ |
Аг - AiE = —3 |
3 |
—1 ) |
|
- 3 |
3 |
- 1 |
/ |
имеет ранг г\ = 1 и потому /i = п — г\ = 3 — 1 = 2 = Jbi.
При А2 = 1 матрица |
|
|
|
|
( |
- 2 |
3 |
- 1 |
\ |
Ai - \ 2Е = I |
- 3 |
4 |
- 1 |
|
\ |
- 3 |
3 |
о |
/ |
имеет ранг т*2 = 2 и потому /2 = п — гг = 3 — 2 = 1 = &2*
Таким образом, у матриць! Ai геометрическая кратность каждого А,- совпадает с его алгебраической кратностью. Поэтому матрица А\ приводится к диагональному виду
2 0 0
А = 0 2 0
0 0 1
В этом можно убедиться и непосредственным конструированием ма^ трицы Т, удовлетворяющей соотношению Т~1А\Т = А. Действи тельно, при А = 2 система (А\ —АЕ)Х = 0, т.е. система
{ —3xi |
+ 3X2 |
— |
Х3 |
= |
0, |
—3xi |
+ 3x2 |
— |
хз |
= |
0, |
-3 x i |
+ 3x2 |
- |
х3 = |
0 |
|
имеет общее решение X |
= (xi, |
Х2, |
—3xi + Зхг)7 , в котором два |
(/1 = &i) свободных неизвестных.
Поэтому возьмем, например, определитель
Полагая в общем решении сначала xi = 1 , хг = 0, затем xi = 0,
=2, соответственно получим частные решения Х\ = (1 , 0, 3)т ,
Х2 = (0, 1, 3)т , которые составляют ФСР рассматриваемой системы уравнений.
При А = 1 система (А\ — АЕ)Х = 0, т.е. система
—2xi |
+ |
3x2 — хз |
= 0, |
||
—3xi |
+ |
4x2 — |
хз |
= |
0, |
—3xi |
+ |
3x2 |
|
= |
0 |
имеет общее решение X = |
|
(xi, xi, |
x i)T, в котором одно (/2 = £2) |
свободное неизвестное. Поэтому ФСР этой системы состоит из од ного решения, например, из решения Х 3 = (1, 1, 1)т Из решений Х\%
Х 2, Хг, как из столбцов, составляется невырожденная матрица
Поэтому матрица А\ приводится к диагональному виду
А= Т~1А1Т =
иимеет каноническое разложение
б) Корнями характеристического многочлена
- А |
З |
3 |
|
\А2 - ХЕ\ = - 1 |
8 - А |
6 |
= - А ( А + 1)2 |
2 |
- 1 4 |
- 1 0 - А |
|
являются Aif2 = —1 , A3 = 0 соответственно кратности к\ = 2, &2 = 1 . Все они действительные и поэтому являются собственными значени ями матрицы А2. Причем при Ai = —1 матрица
имеет ранг г\ = 2 и потому |
= гг — п = 3 — 2 = 1 ^ fci, т.е. |
геоме |
трическая кратность /i = l |
характеристического числа Ai = |
— 1 не |
равна его алгебраической кратности ki = 2. Поэтому матрица А2 не приводится к диагональному виду. Если бы строили для матрицы А2 матрицу Т, удовлетворяющую соотношениям (3.23) и (3.24), то она получилась бы неквадратной.
Действительно, при А = —1 система (А2 — АЕ)Х = 0, т.е. система
+ |
3X2 |
+ |
Зхз |
= |
0, |
—х\ + |
9x2 |
+ |
6x3 |
= |
0, |
2xi — 14x2 - 9хз = 0
имеет ФСР, состоящую из одного решения (/i ф &i), например, из решения Х\ —(3,3, —4)т
При А = 0 система (А2 — AJE7)J\T = 0, т.е. |
система |
|||||
{ |
+ |
3x2 |
+ |
Зхз |
= |
О, |
—®i |
8х2 |
+ |
бхз |
= |
О, |
|
2xi |
- |
14х2 |
- |
Юхз |
= |
О, |
также имеет ФСР, состоящую из одного решения (/2 = Агг), например, из решения Х 2 = (2, 1, — 1)т Решений Х\ и Х 2 недостаточно для конструирования квадратной невырожденной матрицы Т третьего порядка. Поэтому матрица А2 не приводится к диагональному виду и не имеет канонического разложения.
Приведение матриц к диагональному виду и каноническое разло жение матриц широко используется в теории и вычислительной прак тике. Например, если известно каноническое разложение А = ГА Т”"1 матрицы А, то ее m-я степень при натуральном числе m легко нахо дится по формуле
|
XT |
0 ) |
|
|
1 = Т |
(3.25) |
|
|
Г " 1, |
||
|
\ 0 |
А? ) |
|
так как А = 1 |
ТКТ- 1 |
ТА Г’ 1 = |
|
(3.25) сохраняется при m целом отрицательном для невырожденной матрицы А. В частности,
АГ1 |
0 ^ |
|
|
|
А~ 1 = Т U |
Г -1 . |
(3.26) |
||
Ай1 ) |
Один из корней т-й степени из матрицы А определяется формулой
V A = T UVXT |
0 |
|
Г -1 . |
(3.27) |
Действительно, возведя правую часть равенства (3.27) по формуле (3.25) в т-ю степень, получим А. Бели в формуле (3.27) все А; > О, то, беря арифметические значения корней m-й степени из каждого
7 - 1 3 0 7
А», получим единственный корень m-й степени из матрицы А, у которого все характеристические числа положительные. О всех кор нях из матрицы см.[7, гл. VIII].
Решение системы АХ = Ьлинейных уравнений также значительно упрощается, если известно каноническое разложение А = ТА Г -*1. В этом случае от системы ТАТ~гХ = Ъпереходят к системе АТ~гХ = Т _16. Затем вводят обозначение Z = Т~1Х и решают систему AZ = Т~1Ъ. Причем, неизвестным zr+i, ..., zn, при которых множителями стоят Ar+i, ..., А„, придают соответственно произвольные значения Си Cbi •••>Сп-г- В результате получают
Z — (z j, Z2y •••, Z f , Си С2j •••j Cfi—r)
По найденному Z находят
X = r Z = (X1,X2,...,X r ,X r+1, ... ,X n ) x |
(3.28) |
Х(^1, -2^2, . . . , Zry Си С2, •••, Сц—г) |
= |
=(z\Xi + + zrXr) + (C iX r+i + ... + Сп-гХп) = х° + Х 0<?н
(сравни с формулами (2.16) и (2.18) из п.2.6).
П ример 2. Решить систему АХ = (12, 12, —8)т , если известно каноническое разложение:
А = ТА Г” 1 = |
/ 3 2 |
|
1 \ / 2 0 0 \ / |
1 —2 —1 \ |
|||||
0 1 |
|
1 0 2 0 |
|
- 2 5 |
3 |
||||
|
\ 2 0 —1 / \ 0 |
0 0 / \ 2 - 4 - 3 / |
|||||||
Реш ение. |
От системы АХ = |
(12, 12, |
—8)т перейдем к системе |
||||||
АТ~гХ = Т ” 1(12, 12, —8)т ,т.е. к системе |
|
|
|
||||||
|
/ 2 |
0 |
0 \ |
/ |
- 4 |
\ |
|
||
|
V |
0 |
2 |
0 |
Т~1Х = 1 |
2 |
о / |
|
|
|
о |
о |
о |
У |
|
V |
|
||
Полагая здесь Т~гХ = Z, получим систему |
|
|
|
||||||
|
2 |
0 |
0 |
|
*1 |
|
- 4 |
|
|
|
0 |
2 |
0 |
|
*2 |
|
12 |
|
|
|
0 |
0 |
0 |
|
|
|
0 |
|
|
или, в подробной записи,
|
|
|
2*i |
= |
-4, |
|
|
|
2Z2 |
= |
12, |
|
|
|
0*з |
= |
0. |
Отсюда находим Z = (—2, б, С)т. Поэтому |
|||||
|
/ 3 |
2 |
1 |
-2 |
б |
X = TZ = I 0 1 |
1 |
6 |
6 + С |
||
|
V 2 |
0 |
-1 |
С |
-4 |
П ример 3. |
Найти квадратный корень из матрицы А при известном |
||||
ее каноническом разложении |
|
|
|||
/ |
17 |
16 |
16 \ |
|
|
А = ( |
16 |
41 |
32 = |
|
|
\ |
16 |
32 |
41 ) |
|
|
= |
/ 1- 2 |
- 2 \ |
/ |
81 0 0 \ . / |
1 |
2 |
2 \ |
||||
1 2 - 1 |
0 |
о ) |
0 |
9 0 |
|
•i |
- 2 5 - 4 |
||||
' \ 2 |
|
1 / \ |
0 90 \9 /—2 |
—4 5 / |
|||||||
Реш ение. |
По формуле (3.27) получаем |
|
|
|
|
||||||
|
/ 1 |
- 2 |
- 2 |
\ |
/ л/81 |
0 |
0 \ |
|
|
||
|
V A = 1 |
2 |
—1 |
0 |
|
|
0 |
л/9 |
О |
х |
|
|
\ |
2 |
0 |
1 |
/ |
\ |
0 |
0 л/9 / |
|
|
|
|
|
- / |
1 |
|
2 2 \ |
|
. / 33 12 12 \ |
|
|||
|
ХН |
" 2 |
|
5 " 4 |
= 9 |
12 |
51 |
24 |
|
||
|
у V |
—2 |
|
—4 |
5 / |
у V 12 |
24 |
51 / |
|
3.9.Упражнения
1.Оператор <р переводит вектор х = (xi, хг, хз) в вектор <р(х). В каких
случаях оператор (р является линейным ?
1 ) |
<р(х) = (x i + 2x2 + Зхз; |
Х2 — 2хз; х\ — 2хг — Зхз)т , |
2) |
<р(х) = (xi +1; Х2 + 2; |
х з+ 3 )т , |
3) |
^(х) = (2 xi + Зха + ®з; xi —Х2 —хз\ xi + 2 x2 + Зхз)т , |
|
4) |
<р(х) = (xi + X2 j Х2 +Х3; Xl + Х2 + Хз)Т i |
|
5) |
<р(х) = (2xi + х2; х \ + х2; хз)т , |
|
6) |
¥>(х) = (xi + 2 X2 + Хз; |
Xl + Х2\ Х\ + Х2 + 2хз)Т |
Вслучае линейности оператора указать его матрицу в том же базисе,
вкаком заданы координаты векторов х и <р(х).
2. Линейный оператор (р переводит векторы ai, a2, аз соответственно в векторы bi, Ь2, Ьз. Найти матрицу оператора в том же базисе, в каком заданы координатами все векторы:
1 ) |
®i |
= |
( 1 , |
2, - 3 ) T, |
02 |
= (0, 1 , 2) T , |
03 = |
( 1 , 0, |
4 )T, |
|||||
|
bi = ( 1 , 1 , 1 ) T , |
62 = ( 1 , 2, 1 ) T , |
Ьз = ( 0, 1 , 1 ) T , |
|||||||||||
2) |
ai |
= |
( 1 , |
2, 1 ) T, |
02 = |
(4, 3, |
—2 )T , |
ез = |
( - 5 , |
- 4 , - 1 ) |
||||
|
bi = (1 , 1 , 1 ) T , |
b* = ( 1 , 0, 1 ) T , |
b3 = (0, - 1 , 1 ) T , |
|||||||||||
з) |
oi |
= |
( 1 , |
1 , 1 ) T , |
02 |
= |
(2, - 3 , 1)T, |
03 = |
(4, |
1, |
—5 )T , |
|||
|
h |
= (о, |
1 , |
o )T , |
62 |
= ( 0, 1 , 1 ) T , |
Ьз = |
(1, |
1, |
0) T , |
||||
4) |
ei |
= |
(4, |
- |
8, —5 )T , |
02 = |
( - 4 , |
7, - 1 ) T, |
03 = |
( - 3 , |
5, |
1)T , |
||
|
bi = ( 1 , 1 , o )T , |
h = (0, 2, 1 ) T , |
b3 = ( 0 , 1, 3 )T , |
|||||||||||
5) |
Ol = |
( 1 , |
2, |
4 )T, |
02 = |
(1, - 3 , 1)T, |
03 = |
( 1 , 1 , |
—5 )T, |
|||||
|
bi = ( 1 , 1 , 1 ) T , |
b2 = (0, 1 , 2) T , |
Ьз = (0, 1 , 3 )T , |
|||||||||||
6) |
oi |
= |
( 1 , |
2, |
4 )T, |
02 = |
( 1 , 1 , |
—5 )T, |
ез = |
(1, |
- 3 , |
1)T, |
||
|
bi = |
( 1 , |
2, |
0) T , |
62 |
= |
(1, 3, |
0 )T , |
Ьз = |
( 0, |
0, |
1 ) T |
3. Линейный оператор (р в базисе ei, е2, ез имеет матрицу А. Найти матрицу этого оператора в базисе еi, e£, е'3, если:
|
|
1 |
|
- 1 |
|
|
1 |
|
е\ |
(1,о,1)т , |
e'i = |
(1, “ 1,1)Т> |
1) |
А = |
2 |
|
1 |
|
- 1 |
> |
е2 |
(2,1,0)т , |
е2 = |
(0,1, —1)Т, |
|
|
|
3 |
|
- 1 |
|
|
1 |
|
ез |
(-3,2,4)т , |
= |
(0,1,1)т , |
|
А = |
( |
1 |
|
2 |
|
3 |
|
ei |
(1>4, —5)т , |
езei = (2,1,2)т , |
|
2) |
|
0 |
|
1 |
|
2 |
> |
е2 |
(2,3, —4)т , |
е2 = (1, —1,2)т , |
||
|
|
- 1 |
|
2 |
|
1 |
|
ез |
(1 ,-2 ,-1 )т , |
е'3 |
=(0,1,—1)т , |
|
|
А = |
( |
1 |
0 |
|
1 |
> |
|
ei |
(1>2,4)т , |
ei =(1,3,2)т , |
|
3) |
|
1 |
1 |
|
1 |
|
|
е2 |
(1,-3,1)т , |
*2 = (1, —2, —1)т , |
||
|
|
- 1 |
2 |
|
1 |
) |
|
ез |
(1,1, -5)т , |
ез = (0,1,1)т , |
||
|
А = |
(2 |
|
3 |
|
- 1 |
|
ei |
(4, —4, —3)т , |
е( =(3,1,2)т , |
||
4) |
1 |
|
- 1 |
|
|
2 |
> |
е2 |
(—8,7,5)т , |
е2 = (1,2,1)т , |
||
|
|
1 |
|
0 |
|
|
1 |
|
ез |
(—5, —1,1)Т , |
ез = (0,0,1)т . |
|
|
А = |
(1 |
|
2 |
|
3 |
> |
|
ei |
(1,2,4)т , |
e i = ( l , l , l ) T, |
|
5) |
3 |
|
1 |
|
2 |
|
|
е2 |
(1,-3,1)Т , |
ei = (0,l,2)T, |
||
|
|
1 |
|
3 |
|
2 |
) |
|
ез |
(1,1,- б ) т , |
ез = (0,1,3)т , |
|
|
А = |
(3 |
|
1 |
1 |
\ |
|
ei |
(1,0,1)т , |
ei = (2,1,1)т , |
||
|
2 |
|
3 |
1 |
|
|
|
е2 |
(2,1,0)т , |
ei = (3,2,l)T, |
||
6) |
|
(1 |
|
1 |
2 |
) |
|
|
ез |
(—3,2,4)т , |
ез = (0,0,1)т |
4. Линейный оператор (рпереводит вектор х = (zi, х2, хз)т в вектор (рх. Найти образ вектора а и прообраз вектора у, если
1 ) |
<рх = (xi + |
2х2 + |
Хз, х2 + 2Z 2 + 2z3, xi + х3) т , |
2) |
a = (1, 1, |
1)т , |
у = (1, 2, 3)т ; |
<fix = (2xi, |
3xi + х 2, Х1 + 2х2 + хз)т, |
||
|
a = (2, 1, |
1)т , |
» = (1, 1, 0)т ; |