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

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

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

Из второго уравнения системы получаем х\ = ягЕсли сложить пер­ вое уравнение со вторым, то придем к уравнению

50zi -|- 2(4 + Зг)х3 = 0.

Отсюда находим

Пусть хз = —25. Тогда х\ = хъ = 4 + Зг, следовательно, х —

=(4+Зг; 4+Зг; —25)т Поэтому одним из ортогональных базисов про­ странства [/з, содержащим векторы ai, 02, является базис, состоящий из векторов ai, 02, х.

Квадратную матрицу Q, для которой сопряженная матрица Q* =

=Q” 1, называют унитарной. Иначе, квадратную матрицу Q назы­ вают унитарной, если выполнены условия

Q*Q = QQ* = Е.

В унитарной матрице столбцы (строки) составляют ортонормированную систему столбцов (строк). Произведение унитарных матриц является унитарной матрицей. Унитарные матрицы обладают и дру­ гими важными свойствами. Здесь существенным является то, что матрицей перехода от одного ортонормированного базиса унитар­ ного пространства к другому ортонормированному базису является унитарная матрица. Очевидно, действительная унитарная матрица является ортогональной.

6.7.Линейные операторы в евклидовом п ростран стве

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

Линейный оператор у?*, действующий в евклидовом пространстве ЕП) называют сопряж енным с оператором , если для любых век­ торов х и у из Еп выполняется равенство

0<рх,у) = (х,(р*у).

(6.27)

П - 1 3 0 7

Для любого линейного оператора у?, действующего в Еп, сопряжен­ ный оператор <р* существует и единственный. Операторы и <р* являются взаимно сопряж енными. Из определения сопряженного оператора * вытекают следующие его свойства:

1.(?•)* = <р,

2.(<р + фу = <р* + Г ,

3.(окр)* = &<р* при любом действительном числе а,

4.=

Матрицы А и Ai соответственно операторов у? и у?* в произвольном базисе е пространства Еп связаны соотношением

Ах = Г-1 А ТГ,

(6.28)

где Г — матрица Грама базиса е. В частности, если базис ортонормированный, то

Ах = Ат,

(6.28')

т.е. в ортонормированном базисе матрицей сопряженного операто­ ра <р* является транспонированная матрица к матрице оператора (р. Приведем пример на применение этих формул.

П ример 1. Линейный оператор в базисе

e/1 = ( l , l ) l ) T ,

e'2 = ( 0 , 1 , 1)Т .

4 = ( 0 , 0 , 1 ) т

имеет матрицу

 

1

0

0

А =

2

1

0

 

3

2

1

Найти матрицу А\ сопряженного оператора ip* в этом же базисе, если векторы е\, е2, е3 заданы координатами в ортонормированном базисе ei, е2, ез.

Реш ение. Векторы е^, е2, е3 заданы координатами в ортонорми­ рованном базисе ei, е2, ез. Поэтому мы можем вычислить скалярные произведения (е(-,е'), г, j = 1,2,3 по формуле (6.11) из п. 6.3 и соста­ вить матрицу Грама

(

(4 ,« i)

 

(ei,e'2)

(е[,е'3) \

 

3

2

1

 

Г =

( е 2

>

4

)

)

((^Л)е 2 > 4

) =

2

2

1

/

\

(

4

.

4

( е з ,

е 2 )

1 1 (

14 , 4 )

векторов базиса е'. Теперь по формуле (6.28) можно найти матрицу А\ оператора <р* в базисе е':

 

- 1

О

1

2

3 \

Ai = Г“ 1ЛТГ

2

-1

О

1

2

х

 

- 1

2

0

0

1 /

 

 

 

6

 

5

3

 

1

\ )

-3

 

-2

-1

 

1

1 /

- 2

 

- 2

- 1

Задачу можно решить иначе следующим образом. Выпишем маг

трицу перехода

 

 

1

0

0

Г = | 1

1

0

1

1

1

от базиса е к базису е' и найдем матрицу

А = ТАТ- 1 =

оператора ip в базисе е. Теперь по формуле (6.28') выпишем матрицу

1

2

3

А х = Ат- 0

1

2

0

0

1

оператора * в ортонормированном базисе е. По найденной матрице Ai в базисе е получим матрицу

 

1

 

0

0 \

1 1

2

3

 

Ах = Т -'А хТ

-1

 

1

0

0

1 2

 

 

0

 

-1

1 )

V 0

0

1

 

 

/ 1

0

 

0

 

6

5

3

х

1

 

1

0

 

-3

-2

- 1

 

V 1

 

1

1

 

-2

-2

-1

оператора <р* в базисе е'.

Пример 2. В базисе е скалярное произведение задано формулой

(х, у)

3xii/i + Х2У2 + 2хзуз + Я1У2 + ®2У1

 

-

2х 1уз - 2х3у! - х2уз - х3у2,

(6.29)

а линейный оператор — матрицей

Найти матрицу А сопряженного оператора <р* в том же базисе е.

Реш ение. По коэффициентам формулы (6.29) задания скалярного произведения в базисе е составляем матрицу Грама

базиса е. Теперь по формуле (6.28) находим матрицу

Для решения задачи другим способом можно построить ортонормированный базис пространства и поступить так, как сделано при решении вторым способом задачи предыдущего примера. Однако при заданном законе (6.29) скалярного умножения векторов вычисле­ ния будут более громозкими.

Областью значений сопряженного оператора <р* является подпро­ странство, ортогональное к ядру оператора <р. Это следует из того, что для любого х из ядра оператора и любого у выполняется соот­ ношение

(*,^*tf) = (р*.У) = (О,У) = 0, т.е. <pmyLx.

Основное свойство сопряженного оператора состоит в следующем:

если некоторое подпространство L инвариантно относительно оператора <р, то ортогональное дополнение LL этого подпро­ странства инвариантно относительно сопряженного опера­ тора (р*

Действительно, пусть х Е L и у G LL. Тогда из условия <рх G L следует (у?х,у) = 0 , а так как (ipx,y) = (х,<р*у), то и (х,(р*у) = 0 . Отсюда следует, что tp*y Е L1 . Это доказывает утверждение.

Характеристические многочлены, а следовательно, и собственные значения сопряженных операторов одинаковы. На самом деле, по­ скольку определитель не меняет значения при транспонировании, то характеристические многочлены —\Е\ и \АТ —ХЕ\ матриц А и Ат будут одинаковыми. Далее, так как характеристические многочлены подобных матриц равны, то получаем

|А - ХЕ\ = \АТ - \Е\ = |Г"1АТГ - ХЕ\ = \Аг - ХЕ\.

Каждый собственный вектор сопряженного оператора ip* ортого­ нален ко всем собственным векторам оператора у?, принадлежащим другим собственным значениям.

Действительно, пусть (рх = А,я, <р*у = Ау т/, А,- ф Aj. По определе­ нию сопряженного оператора выполняется условие (<рх, у) = (х, <р*у)> которое в данном случае принимает вид At(x,j/) = Aj(x,y). При А,- ф Aj это возможно лишь при (х, у) = 0.

Для матриц последнее утверждение перефразируется так:

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

Линейный оператор р ) действующий в евклидовом пространстве ЕПу называется сим м етрическим или самосопряж енным, если для любых векторов х и у из Еп выполняется равенство

(<рх,у) = (х,<ру).

(6.30)

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

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

означает, что симметрический оператор является оператором простой структуры, а среди матриц Т, приводящих симметрическую матрицу А к диагональному виду, т.е. удовлетворяющих соотноше­ нию Т~1АТ = = А, где А - диагональная матрица, есть ортогональ­ ная матрица или, что то же самое, симметрическая матрица А обладает каноническим разложением А = ТАТ~г, в котором транс­ формирующая матрица Т ортогональная. Правило построения та­ кой ортогональной матрицы Т остается таким же, как в случае любых операторов простой структуры (см. п. 3.8) с той лишь раз­ ницей, что базис из собственных векторов матрицы А здесь еще и ортонормируют. Поясним это правило на примере.

П ример 3. Построить каноническое разложение симметрической матрицы

0

1

1

- i

\

1

0

-1

1

 

1

-1

0

1

 

-1

1

1

0

)

с ортогональной трансформирующей матрицей.

Реш ение. Характеристический многочлен матрицы А

 

1

1

- 1

 

ХЕ\ =

1

- 1

1

= ( А - 1 ) 3(А + 3)

1

- 1

1

 

 

 

- 1

1

1

 

имеет корни Ai = А2 = A3 = 1, А4 = —3. Все они являются собствен­ ными значениями матрицы А.

При А = 1 система — АЕ)Х = 0 имеет ФСР, состоящую из трех решений, например, из решений ai = ( 1 , 1 , 0, 0)т , аг = (1 , 0, 1 ,0)т , а3 = (—1 ,0 ,0 ,1)т Так как мы строим каноническое разложение ма­ трицы А с ортогональной трансформирующей матрицей, то ортонормируем систему собственных векторов а\, аг, а3. В результате получим ортонормированную систему векторов

ei = -L(1,1,0,0)\ «2 = i ( l , - l , 2 , 0 ) T, 4 = ^ ( - 1 , 1 , U ) T

При A = - 3 система (A - AE)X = 0 имеет ФСР, состоящую из одного вектора, например, из вектора 04 = (1, —1, —1 ,1)т Нормируя его, получим в\ = ^(1, 1, 1, 1Собственные векторы е*, е1^, е3, е\

матрицы А составляют ортонормированный базис в £ 4. Из столбцов координат этих векторов составим матрицу Т и запишем искомое каноническое разложение

 

 

 

/ X

 

X X L

1 ч

 

 

 

 

/

v/2

 

ч/б2^5

_2

\

 

А = ТА Г - 1

= ТАТт=

\?2

 

t/б 2 ^ 3

~

 

 

о

 

4=

1

- 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ч/б

г/з

 

 

 

 

 

 

 

»

 

 

°

3 I

7/ з

 

/1

0

0

0 \

^ 72

*

»

 

 

 

 

-1

42т

 

о

0

1

0

Тб

 

0

у/б

у/б

 

 

0 0 1 0

_-1_

_ J L _

_ j _

_ з _

0

0

- ч

2 у/ 3

2 у/ 3

2 у/ 3

2 у/ 3

 

 

 

 

 

 

 

 

 

 

 

1 = 1 .

=1

1

/

 

 

 

 

 

2

 

2

2

2

Представление симметрических матриц в виде канонического раз­ ложения с ортогональной трансформирующей матрицей имеет самые широкие применения в теории и приложениях. Для построения та­ кого канонического разложения применимы и косвенные способы, на­ пример, метод вращений (см. п. 9.2).

Симметрический оператор р называют неотрицательны м (п о­ ложительно определенным), если для любого ненулевого вектора х выполняется неравенство (рх}х) > 0 ( (<рх,х) > 0 ). Неотрицатель­ ный и положительно определенный операторы обозначают соответ­ ственно через р > 0 и > 0.

Симметрический оператор является неотрицательным (положи­ тельно определенным) тогда и только тогда, когда все его собствен­ ные значения неотрицательные (положительные). Для любого нео­ трицательного симметрического оператора р существует такой нео­ трицательный оператор / , что f m = p при любом натуральном га. Симметрический оператор / , удовлетворяющий условию / т = р } на­ зывают ариф м етическим корнем га-и степени из неотрица­ тельного си м м етри ческого оператора р .

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

Операторы у / р * р и у/ р р * называют соответственно правым и левым модулями оператора р .

На матричном языке все сказанное здесь перефразируется следу­ ющим образом (сравни с результатами гл. 7).

Симметрическая матрица А называется неотрицательной (по­ лож ительно определенной), если для любого вектора х ф О выпол­ няется условие (Ахух) = хтАх > 0 ((Аж,х) = хтАх > 0). Неотрица­ тельная (положительно определенная) матрица А обозначается через А > 0 > 0). Симметрическая матрица является неотрицатель­ ной (положительно определенной) тогда и только тогда, когда все ее характеристические числа неотрицательные (положительные).

А риф м етическим корнем m-й степени из неотрицатель­ ной матрицы А называют такую матрицу J3, что Вт = А. При вычислении арифметического корня Гу/А следует построить канони­ ческое разложение А = ТА Г ” 1 и положить Тл/А = Т т\/АТ~ 1 , где

1

/ т

 

о

\

( К

 

 

г\/А =

 

 

. 1 / т

 

 

 

 

 

 

о

 

о

/

Поясним это правило на примере.

 

 

П ример 4. Для матрицы

/1 7

16

 

 

 

16

 

А =

16

41

32

 

V16

32

41

 

построить каноническое разложение с ортогональной трансформиру­ ющей матрицей и, пользуясь им, найти \[А.

Реш ение. Характеристический многочлен

 

|А - ХЕ\

17- А

16

16

16

41 - А

32

1632 41 - А

-А 3 + 99А2 - 1539А + 6561

имеет корни Ai = 81, Аг = A3 = 9. Поэтому матрицей Л является

 

81

0

0

Л =

0

9

0

 

0

0

9

Далее будем конструировать ортогональную матрицу Т. Для этого будем строить ФСР однородных систем AiE)X = 0 И ортонормировать их векторы-решения.

При Ai = 81 эта система имеет вид

 

 

—64a:i + 16*2 + 16*з

=

О,

16*х

- 40*2

+ 32*з

=

0,

16*i

+ 32*2

—40*з

=

0.

Ее общим решением является X = (|*з, *з, *з)т Здесь одно сво­ бодное неизвестное. Поэтому ФСР состоит из одного решения, наг пример, из решения bi = (1 , 2, 2)т Нормируя его, получим е[ =

= i ( l , 2, 2)т

При А2 = Аз = 9 рассматриваемая система имеет вид

{ 8*1 + 16х2 + 16ж3

=

0,

16xi + 32х2 + 32хз

=

0,

16xi

+ 32х2

+ 32хз = 0.

Ее общим решением является X = (—2х2 — 2хз, х2,хз)т Здесь два

свободных неизвестных.

Поэтому ФСР состоит из двух решений,

например, из решений

 

 

Ь2 = ( - 2 , - 1 , 2

и 63 = (2, - 2, 1

Эти векторы уже ортогональны. Поэтому остается их лишь нормиро­ вать. В результате получим ортонормированную систему векторов

4 = | (-2 , -1 ,

2)т ,

е'3 = 1(2, -2 , 1)т

 

Из столбцов координат векторов

 

, е^, е!3 строится ортогональная

И

матрица

 

21 - -21

- 22

 

 

2

2

1

А имеет вид

Поэтому искомое каноническое разложение матрицы

Теперь находим

у/А =

Примечания.

1.При вычислении rv^4 можно пользоваться каноническим разло­ жением матрицы А необязательно с ортогональной трансфор­ мирующей матрицей (см. п. 3.8).

2.При вычислении на ЭВМ корней из симметрических матриц ка­ ноническое разложение матрицы естественно получать методом вращений (см. п. 9.2).

3.О всех корнях из матрицы см.[7, гл. VIII].

Линейный оператор (р, действующий в евклидовом пространстве Еп, называется ортогональны м , если он сохраняет скалярное про­ изведение любых векторов х ъ у из Еп, т.е. если выполняется равен­ ство

(<рх, <ру) = (X, у).

(6.31)

Полагая в равенстве (6.31) х = у, получаем \фх\2 = |х|2. Это означает, что ортогональный оператор сохраняет длины векторов.

Поскольку угол между векторами определяется из соотношения

С 0 8 ¥>= и П и ’ 0 < < Р < 7 Г ,

1*1 •|У1

(см. формулу (6.8) из п.6.2) и поскольку числитель' и знаменатель в этом соотношении не меняются при действии ортогональных опера­ торов, то ортогональные операторы сохраняют углы между векто­ рами. Из сказанного вытекает, что ортогональные операторы любую ортонормированную систему векторов переводят в ортонормированную систему векторов, а ортонормированный базис - в ортонормированный базис.

Оператор р*, сопряженный с ортогональным оператором , со­ впадает с ~ 1

Корни характеристического многочлена ортогонального операто­ ра (в том числе и комплексные) по абсолютной величине равны еди­ нице; действительные корни равны ±1.

Собственные векторы ортогонального оператора, принадлежащие различным собственным значениям, ортогональны.

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