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

Lineynye_operatsii

.pdf
Скачиваний:
166
Добавлен:
05.06.2015
Размер:
5.31 Mб
Скачать
. Для это-

2.5. Приведение матрицы к диагональному виду

 

 

71

Записывая координаты

элементов e1 , e2 , …, en

 

образов базисных элементов

e1 , e2 ,

 

, en

по столбцам, получаем

матрицу A оператора в базисе e1, e2 ,

, en :

 

 

 

0

 

0

 

 

 

 

1

 

 

0

 

 

 

A

 

2

 

 

.

 

 

 

 

 

 

 

 

 

.

.

.

.

 

 

 

 

0

0

 

n

 

 

 

 

 

 

 

 

Теорема 2.15 доказана.

Следствие (достаточное условие приведения матрицы линейного оператора к диагональному виду). Если линейный оператор

, действующий в n -мерном линейном пространстве V , имеет n

различных собственных значений, то существует базис, в котором матрица этого оператора имеет диагональный вид.

Доказательство. Пусть 1 , 2 ,

, n

различные собственные

значения оператора , а e1 , e2 , , en

 

отвечающие им собствен-

ные векторы. По свойству 3 предыдущего параграфа они линейно независимы. Следовательно, они образуют базис. По теореме 2.15 в базисе из собственных векторов матрица оператора имеет диагональный вид.

0

1

0

 

 

 

1

0

0

 

матрица оператора в

Пример 2.10. Пусть A

 

 

0

0

0

 

 

 

 

 

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

Решение. Найдем собственные значения оператора го составим и решим характеристическое уравнение

72

 

 

 

 

Глава 2. ЛИНЕЙНЫЕ ОПЕРАТОРЫ

 

-

1

0

 

0 2 1 0 1 0, 2 1, 3 1.

 

 

 

1

-

0

 

 

0

0

-

 

 

Так как мы получили три различных собственных значения, то по следствию к теореме 2.15 существует базис из собственных векторов, в котором матрица оператора имеет диагональный вид.

Найдем их.

Обозначим базис, в котором задана матрица оператора , через

e1 , e2 , e3 . Найдем собственный вектор

f1 , отвечающий соб-

ственному значению 1

0 .

Он является ненулевым решением

однородной системы линейных алгебраических уравнений

 

 

 

0

 

1

0 1

 

 

0

 

0

 

A 0E 0

 

 

 

 

 

 

2

 

 

 

0

 

1

0

0

 

1

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

0

 

0

0

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Полагая 3 1 , получим

f1 0, 0,1 e3.

 

 

 

 

 

 

 

 

Аналогично найдем собственные векторы

 

f2

и f3 , отвечающие

собственным значениям 2

1 и 3

1, соответственно.

 

 

 

 

 

 

 

-1 1

 

 

0 1

 

0

 

 

A E 0

 

1 -1

 

 

0

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

-1

3

 

 

 

 

1 2 0

 

1 2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая 2

1,

получим

f2 1,1, 0 e1

e2.

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

 

0 1

 

0

 

 

A E 0

1 1

 

0

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

3

 

 

0

 

 

 

 

 

 

 

 

 

 

 

-1

 

 

 

 

2.5. Приведение матрицы к диагональному виду

73

 

 

 

0

 

 

 

2 .

 

1

 

2

 

1

0

 

 

0

 

 

3

 

 

3

 

 

 

 

 

 

Полагая 2

1, получим

f3

1,1, 0

e1

e2 .

Таким образом, мы получили базис из собственных векторов

оператора :

f1 e3 , f2

e1 e2 ,

f3

e1 e2 . В этом базисе

матрица оператора имеет вид

 

 

 

 

 

0

0

0

 

 

f

 

 

 

 

 

 

A

 

0

1

0

.

 

 

 

0

0

 

 

 

 

 

-1

 

1

1

1

1

 

 

 

 

-1

 

 

Пример 2.11. Пусть A

1 1

-1

матрица линейного

1

-1

1

-1

 

 

 

 

 

 

 

1

-1

-1

1

 

оператора из примера 2.9. Требуется выяснить, можно ли при-

вести матрицу этого оператора к диагональному виду путем перехода к новому базису. В случае положительного ответа найти этот базис и соответствующую ему матрицу.

Решение. В примере 2.9 мы нашли собственные значения оператора и отвечающие им линейно независимые собственные век-

торы. Собственному

значению

1 2 отвечают векторы

f1 1,1, 0, 0 ,

f2

1, 0,1, 0 и

f3 1, 0, 0,1 ,

а

собственному

значению 2

2

 

f4 1,1,1,1 . Покажем,

что

векторы f1 ,

f2 , f3 , f4 образуют базис. Составим матрицу, строками которой

являются координаты этих векторов, вычислим ее ранг, приведя матрицу к ступенчатому виду элементарными преобразованиями:

74

 

 

 

 

 

 

 

Глава 2. ЛИНЕЙНЫЕ ОПЕРАТОРЫ

1

1

0

0

1

0

0

0

1

1

0

0

 

 

1

0

-1

 

 

0

-1

-1

 

 

0

-1

 

 

 

 

-1

 

-1

 

-1 -1

.

 

1

0

0

-1

 

0

-1

0

-1

 

0

0

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-1 1

1

1

 

0

2

1

1

 

0

0

0

-1

 

Мы получили, что ранг матрицы равен 4. Это означает, что векторы f1 , f2 , f3 , f4 линейно независимы и, следовательно, образуют базис. В этом базисе матрица оператора имеет вид

2

0

0

0

 

 

0

2

0

0

 

A f

.

 

0

0

2

0

 

 

 

 

 

 

 

0

0

0

-2

 

Заметим, что полученный результат (диагонализуемость оператора ) не является случайностью. Имеет место

Теорема 2.16 (критерий диагонализуемости линейного опера-

тора). Оператор L(V ,V ) диагонализуем тогда и только тогда, когда каждый корень характеристического многочлена принадлежит F и его геометрическая кратность совпадает с алгебраической

кратностью.

 

 

 

 

 

Пусть L(V ,V )

 

Доказательство. Необходимость.

диаго-

нализуемый оператор. Выберем в V

базис

e1 , e2 ,

, en , в

котором матрица A оператора имеет диагональный вид

1

0

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

A

1

 

 

diag( ,..., , ,..., ,..., ,..., ) ,

 

 

 

.

.

.

.

 

1

1 2

2

k

k

 

 

 

 

 

0

0

...

k

 

 

 

 

 

здесь число

i

F, i 1, 2,..., k встречается в матрице

A ni раз,

n1 n2

nk n .

 

 

 

 

 

 

 

2.5. Приведение матрицы к диагональному виду

75

Характеристический многочлен матрицы A (оператора ) равен

P ( ) det A E ( )n1 ( )n2

( )nk .

n

1

2

k

Отсюда следует, что число

i

является собственным значением

оператора алгебраической кратности ni , i 1, 2,..., k. Рассмотрим матрицу A i E . Ее ранг равен n ni . Тогда по теореме 2.13

собственному значению i отвечает n r n n ni ni линейно независимых собственных векторов, следовательно, геометрическая кратность i также равна ni .

Достаточность. Пусть любой корень характеристического уравнения оператора принадлежит множеству F (если F множество комплексных чисел, то это условие выполняется автомати-

чески, а если F множество вещественных чисел, то оно означает, что все корни характеристического уравнения являются вещественными числами). Тогда характеристический многочлен оператора имеет вид

P ( ) det A E

( )n1 ( )n2

 

( )nk ,

n

1

2

 

k

n1 n2

nk n .

 

 

Из теоремы 2.13 следует, что числа 1 , 2 , …, k

являются собст-

венными значениями оператора алгебраической кратности

n1 ,

n2 ,..., nk ,

соответственно. По условию алгебраические кратности

собственных

значений

совпадают

с геометрическими:

ni li ,

i 1, 2,..., k.

Пусть

ei , ei , …,

ei

линейно независимые собст-

 

 

 

 

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

венные

векторы,

отвечающие

собственному

значению

i ,

i 1, 2,..., k. Покажем, что векторы

e1 ,

e1

, …, e1

, …, ek ,

ek , …,

 

 

 

 

 

 

 

 

1

2

n

1

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

ek

линейно независимы. Рассмотрим их линейную комбинацию и

n

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

приравняем ее нулевому вектору:

 

 

 

 

 

 

 

 

 

1e1

1

e1

k ek

 

k ek .

(2.19)

 

 

1

1

n

n

 

1

1

 

n

n

 

 

 

 

 

 

1

1

 

 

 

 

k

k

 

 

76

 

 

 

Глава 2. ЛИНЕЙНЫЕ ОПЕРАТОРЫ

Обозначим xi i ei

i

ei , i 1, 2,..., k и перепишем (2.19)

 

1 1

 

n

n

 

 

 

 

 

 

i

i

 

 

 

в виде

 

 

 

 

 

 

 

 

 

x1 x2

xk

.

(2.20)

Покажем, что

x1 x2

xk

.

Предположим противное.

Пусть среди векторов

x1 , x2 , …,

xk есть ненулевые векторы (из

(2.20) очевидно,

что k 1), например, x1 , x2 , …, xm

( 2 m k ).

Тогда, с одной стороны, они являются собственными векторами, отвечающими различным собственным значениям и, следовательно, по свойству 3 линейной независимости собственных векторов, отвечающих различным собственным значениям, они линейно не-

зависимы. С другой стороны, имеем x1 x2 xm , что означает их линейную зависимость. Получили противоречие. Следо-

вательно, x1 x2

 

xk , т.е.

i ei

 

 

i

ei

, i 1, 2,..., k .

1 1

 

 

 

 

n

n

 

 

 

 

 

 

i

i

 

Так как по условию

ei

,

ei

 

, …, ei

линейно независимы, то

 

1

 

2

 

n

 

 

 

 

 

 

 

i

 

i i

 

 

i

0 , i 1, 2,..., k .

1

 

2

 

 

 

n

 

 

 

 

 

 

 

i

 

Итак, мы получили, что в равенстве (2.19) все коэффициенты равны нулю. Следовательно, собственные векторы e11 , e12 , …, e1n1 ,

…, e1k , e2k , …, enkk линейно независимы. А так как их количество равно размерности пространства V , то они образуют базис. Со-

гласно теореме 2.15 в базисе из собственных векторов матрица оператора имеет диагональный вид. Теорема 2.16 доказана.

Замечание 2.10. Линейная независимость собственных векторов f1 , f2 , f3 , f4 в примере 2.11 автоматически следует из доказанной выше теоремы.

77

Глава 3. ЖОРДАНОВА НОРМАЛЬНАЯ ФОРМА МАТРИЦЫ

Пусть V

 

линейное пространство над множеством F,

dimV n,

 

линейный оператор, отображающий V в себя

( L(V ,V )).

 

В этой главе

будем предполагать, что корни характеристиче-

ского многочлена Pn ( ) оператора лежат в F (в комплексном

линейном пространстве в силу основной теоремы алгебры это требование выполнено автоматически).

Цель данной главы доказать, что в этом случае в пространстве V существует базис, в котором матрица оператора имеет спе-

циальный вид, называемый жордановой нормальной формой матрицы оператора :

J1

 

0

 

0

 

 

 

 

 

 

0

J2

 

0

 

diag J1, J2 ,

, Jd ,

J

 

 

 

 

.

 

. . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

Jd

 

 

 

 

 

где матрицы Ji ,

называемые жордановыми клетками, имеют вид

 

i

1

0

 

 

0

0

 

 

 

 

 

0

i

1

 

 

0

0

 

 

 

 

 

 

 

 

, i

F, −

Ji . . .

 

. .

.

 

 

 

0

0

0

 

...

i

1

 

 

 

 

 

 

 

 

 

 

 

0

0

0

 

...

0

 

 

 

 

 

 

 

i

 

 

матрица порядка mi , i 1, 2,

 

, d , m1 m2

 

md n.

78Глава 3. ЖОРДАНОВА НОРМАЛЬНАЯ ФОРМА МАТРИЦЫ

3.1.Корневые векторы и корневые подпространства

Определение 3.1. Ненулевой элемент (вектор) x0 V называет-

ся корневым вектором оператора , отвечающим числу 0 , если при некотором натуральном s выполняется равенство

0 I s (x0 ) .

Определение 3.2. Натуральное число m называется высотой корневого вектора x , если

0 I m (x0 ) , а 0 I m 1 (x0 ) .

Утверждение 3.1. Число 0 из определения 3.1 корневого век-

тора, является собственным значением оператора .

 

Доказательство. Пусть m − высота вектора

x0 . Обозначим

y0 0 I m 1 (x0 ). В силу определения 3.2 y0

. С учетом

того, что

0 I m (x0 ) 0 I 0 I m 1 (x0 ) 0 I ( y0 ),

имеем

0 I ( y0 )y0

( y0 ) 0 y0 .y

0

Это и означает, что 0 является собственным значением оператора

, отвечающим собственному вектору y0 .

Утверждение 3.2. Собственный вектор является корневым вектором высоты 1.

Доказательство. Пусть x0 собственный вектор оператора ,

отвечающий собственному значению 0 :

(x0 ) 0 x0x .0

Перепишем это определение в виде

3.1. Корневые векторы и корневые подпространства

79

 

I

(x )

 

 

 

0

 

0

,

 

 

 

0

(x0 )

 

0 I

 

 

 

 

 

 

 

 

 

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

x0 корневой вектор высоты 1. Утверждение 3.2

доказано.

 

 

Определение 3.3. Подпространство L линейного пространства

V называется инвариантным относительно оператора L(V ,V ) ,

если x L

(x) L.

Пример

3.1.

Пусть 0 собственное значение оператора

L(V ,V ) , а x0 отвечающий ему собственный вектор. Покажем,

что L x0 {cx0},

c F,

является одномерным подпространст-

вом, инвариантным относительно оператора .

 

 

Действительно,

для

любого

элемента

x L x0

(x) (cx0 ) c (x0 ) c 0 x0 L(x0 ).

 

 

 

 

Теорема 3.1. Пусть 0

собственное

значение

оператора

L(V ,V ) геометрической

кратности l,

а

x1 , x2 ,

, xl

линей-

но независимые собственные векторы, отвечающие этому собственному значению. Тогда подпространство

L x1, x2 ,

, xl {c1x1 c2 x2

cl xl },

 

ci F,

инвариантно относительно оператора ,

его размерность равна l.

Доказательство. Возьмем произвольный элемент

 

x c1x1 c2 x2

cl xl L x1, x2 ,

, xl ,

и покажем, что (x) L x1, x2 ,

, xl

:

 

 

 

(x) (c1x1 c2 x2

cl xl )

 

 

c1 (x1 ) c2 (x2 )

cl (xl )

 

 

c1 0 x1 c2 0 x2

cl 0 xl

 

 

0 (c1x1 c2 x2

cl xl ) 0 x L x1, x2 ,

, xl .

V 0 .

80

Глава 3. ЖОРДАНОВА НОРМАЛЬНАЯ ФОРМА МАТРИЦЫ

 

Поскольку

x1, , xl − линейно независимые порождающие эле-

менты подпространства L x1, x2 ,

, xl , то в силу теоремы 1.12

его размерность равна l. Теорема 3.1 доказана.

 

 

Определение 3.4. Подпространство L x1 , x2 ,

, xl называется

собственным подпространством оператора , отвечающим собственному значению 0 , и обозначается

Из определения собственного подпространства V 0 следует, что его размерность равна геометрической кратности собственного значения 0 .

Теорема 3.2. Множество U , состоящее из всех корневых векторов оператора , отвечающих собственному значению 0 и ну-

левого элемента , является линейным подпространством, инвариантным относительно оператора .

Доказательство. Пусть x и y корневые векторы, отвечающие собственному значению 0 , высоты r и s, соответственно. Возьмем m max r, s . Вычислим 0 I m x y :

0 I m x y 0 I m (x) 0 I m ( y)

(здесь мы воспользовались линейностью оператора 0 I m ). Следовательно, x y либо нулевой элемент, либо корне-

вой вектор, отвечающий собственному значению 0 , высоты m .

Это означает, что множество U является подпространством линейного пространства V . Докажем, что U инвариантно относительно оператора . Возьмем произвольный элемент x U и по-

кажем, что (x) U. Если x , то в силу линейности оператора и, следовательно, ( ) U. Пусть теперь x . То-

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]