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

книги / Линейная алгебра

..pdf
Скачиваний:
32
Добавлен:
12.11.2023
Размер:
13.06 Mб
Скачать

в характеристическом многочлене) называют сп ектром линейного оператора. Если матрицу отождествляют с оператором, то множе­ ство всех ее собственных значений называют сп ектром матрицы .

Из сказанного следует, что длл отыскания всех собственных зна­ чений оператора с матрицей А нужно найти все характеристиче­ ские числа матрицы А и из них выбрать лишь те, которые принад­ лежат основному полю, а длл отыскания всех собственных векто­ ров оператора с матрицей А нужно найти все ненулевые решения системы (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)т ;