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

Е. Е. Тыртышников

127

 

 

поменять их местами с помощью поворота на угол π/2). Åñëè ïðè ýòîì c > 0, то изуча- емое множество пусто. Если c = 0, в нем только одна точка (0, 0). Предположим, что

pp

c < 0. Тогда, положив a = −c/λ1, b = −c/λ2, уравнение (1) можно записать в виде

x2

 

y2

(10)

 

+

 

 

= 1, a ≥ b > 0.

a2

b2

Определение 1. Множество точек (x, y), удовлетворяющих уравнению (10), называ-

åòñÿ эллипсом с полуосями a è b.

Точки F= (−c, 0) è F+ = (c, 0), ãäå c = a2 − b2 ≥ 0, называются, соответственно, отрицательным и положительным фокусами эллипса. Число e = c/a называется ýêñ-

центриситетом эллипса. Заметим, что 0 ≤ e < 1. Прямые l: x = −a/e è l+ : x = a/e называются, сответственно, отрицательной и положительной директрисами эллипса.

Пусть точка M = (x, y) удовлетворяет уравнению (10). Найдем сумму расстояний от нее до фокусов:

|MF| + |MF+| =

p

 

p

 

 

 

 

 

 

 

(x + c)2 + y2 + (x − c)2 + y2

 

p

 

+ p

 

.

=

x2(1 − b2/a2) + 2xc + b2 + c2

x2(1 − b2/a2) − 2xc + b2 + c2

Заметим, что 1 − b2/a2 = e2 è b2 + c2 = (c/e)2 = a2. Кроме того, |ex| ≤ a. Поэтому

|MF| + |MF+| =

p

 

 

 

 

 

 

 

 

= (

+

)

p

(ex)2 − 2(ex)(c/e) + (c/e)2

|

(ex)2

+ 2(ex)(c/e) + (c/e)2

+

 

=

a + ex

|

+

|

a

ex

|

 

a ex + (a

ex) = 2a.

 

 

 

 

 

 

 

Таким образом, сумма расстояний от любой точки эллипса (10) до его фокусов постоянна и равна 2a.

Определение 2. Множество тех и только тех точек плоскости, для которых сумма расстояний до заданных точек постоянна, называется эллипсом.

Мы уже выяснили, что все точки эллипса как множества из определения 1 принадлежат множеству из определения 2. Рассмотрим теперь эллипс как множество, данное

определением 2. Выберем декартову систему, в которой заданные точки Fè F+ ïîëó- чают координаты (−c, 0) è (c, 0). Постоянную сумму расстояний будем считать равной

2a. Тогда

p

 

2

= 2a − p

 

2

a p

 

= a2 + xc.

(x − c)2 + y2

(x + c)2 + y2

(x + c)2 + y2

Еще одно возведение в квадрат дает определения 1 и 2 эквивалентны.

b2x2 + a2y2 = a2b2

 

(10). Следовательно,

 

 

В случае e = 0 эллипс есть окружность радиуса a = b. Пусть e > 0. Выше мы получили равенства

|

MF

|

=

a + ex

= e

x + (a/e

 

|MF|

= e,

 

|

 

 

|

|

 

|

|x + (a/e)|

 

|

MF

 

=

a

ex

= e x

(a/e)

 

|MF+|

= e.

 

+|

|

 

 

|

| −

|

|x − (a/e)|

 

128

Лекция 19

 

 

Возводя каждое из последних равенств в квадрат, получаем (10). Таким образом, доказано следующее

Утверждение. Множество тех и только тех точек плоскости, для которых отношение расстояний до заданной точки и заданной прямой постоянно и равно 0 < e < 1,

является эллипсом.

Ясно, что для выбора точки (фокуса) и соответствующей прямой (директрисы), определяющих один и тот же эллипс, имеются в точности две возможности.

Задача. Написать общее уравнение прямой, проходящей через точку M(x0, y0) эллипса x2/a2 + y2/b2 = 1 и имеющей с ним единственную общую точку (такая прямая называется касательной ê

эллипсу в точке M). Доказать, что прямая, ортогональная к данной прямой и проходящая через точку

M, является биссектрисой угла AMB, ãäå A è B фокусы эллипса.

Задача. Докажите, что преобразование z 7→12 (z + z−1) комплексной плоскости 3 переводит точки окружности радиуса r > 1 с центром в начале координат в точки некоторого эллипса.

Задача. Докажите, что точки эллипса x2/a2 + y2/b2 = 1 допускают параметрическое представление x = a cos φ, y = b sin φ, 0 ≤ φ < 2π.

Задача. Докажите, что геометрическое место точек пересечения взаимно ортогональных касательных к эллипсу x2/a2 + y2/b2 = 1 есть окружность x2 + y2 = a2 + b2.

19.6Гипербола

По-прежнему, пусть в декартовой системе координат уравнение f(x, y) = 0 получает âèä (1), íî λ1 è λ2 имеют разные знаки. Если при этом свободный член оказался равен

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

x2 y2 = 1, (100) a2 b2

ãäå a, b некоторые положительные числа (возможно, для того потребуется дополнительно повернуть систему координат на угол π/2).

Определение 1. Множество точек (x, y), удовлетворяющих уравнению (100), называется гиперболой с полуосями a è b.

Легко видеть, что точки (x, y) гиперболы (100) находятся в объединении двух непересекающихся областей плоскости (как говорят, распадаются на две ветви)

D+ = {(x, y) : x ≥ a, |y| ≤ (b/a)|x|}, D= {(x, y) : x ≤ −a, |y| ≤ (b/a)|x|}.

Прямые h+ : y = (b/a)x è h: y = −(b/a)x называются асимптотами гиперболы. Пусть x > 0 è y = y(x) единственное значение для y такое, что y > 0 и точка (x, y)

удовлетворяет уравнению (100). Очевидно, расстояние от точки (x, y((x)) до асимптоты h+ не превышает

b2

|(b/a)x − y(x)| = |(b/a)x + y(x)| → 0 ïðè x → +∞.

3Функция z 7→12 (z + z−1) называется функцией Жуковского и широко применяется при решении задач гидро- и аэродинамики.

M(x0, y0) гиперболы x2/a2

Е. Е. Тыртышников

129

 

 

Ïðè x > 0 è y < 0 соответствующие точки (x, y(x)) гиперболы приближаются к асимп-

òîòå h. Аналогичные наблюдения справедливы также для точек гиперболы при x < 0.

Точки F= (−c, 0) è F+ = (c, 0), ãäå c = a2 + b2 > 0, называются, соответственно, отрицательным и положительным фокусами гиперболы. Число e = c/a называется эксцентриситетом гиперболы. Заметим, что в случае гиперболы e > 1. Прямые

l−1 : x = −a/e è l+ : x = a/e называются, соответственно, отрицательной и положительной директрисами гиперболы.

Найдем расстояния от произвольной точки M = (x, y), удовлетворяющей (100), до фокусов (выкладки проводятся в полной аналогии со случаем эллипса):

 

 

 

 

|MF| = p

 

 

 

(A)

 

 

 

 

(x + c)2

+ y2 = |a + ex|,

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

|MF+| = p

 

= |a − ex|.

 

 

 

 

 

(x − c) + y

(B)

Поскольку |x| ≥ a и e > 1, получаем

 

 

 

|

 

| − |

 

|

 

−(ex + a) + (ex − a) = −2a, x < 0.

 

 

a + ex

 

a

 

ex

=

(ex + a)

− (ex − a) = 2a, x > 0,

 

Таким образом, абсолютная величина разности расстояний от любой точки гиперболы (100) до ее фокусов постоянна и равна 2a.

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

Пусть точка (x, y) принадлежит множеству из определения 2. Введем декартовы координаты таким образом, что заданные точки получают координаты (−c, 0) и (c, 0). Постоянную абсолютную величину разности расстояний обозначим через 2a. Тогда

 

 

 

 

 

 

 

 

= 2a

 

(x + c)2 + y2 − (x − c)2

+ y2

p

 

 

 

 

p

 

 

 

 

2

 

2

+ y

2

=

 

 

 

 

 

 

 

(100).

(x + c)

 

 

 

2a − p(x − c)2 + y2

Таким образом, определения 1 и 2 эквивалентны. Формулы (A) и (B) делают очевидным также следующее

Утверждение. Множество тех и только тех точек плоскости, для которых отношение расстояний до заданной точки и заданной прямой постоянно и равно e > 1,

является гиперболой.

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

Задача. Написать общее уравнение прямой, проходящей через точку

y2/b2 = 1 и имеющей с ней единственную общую точку в области точек (x, y) с координатой x того же знака, что и x0 (такая прямая называется касательной к гиперболе). Доказать, что данная прямая является биссектрисой угла AMB, где A и B фокусы гиперболы.

Задача. Докажите, что никакая прямая не может иметь ровно одну общую точку с каждой ветвью гиперболы.

Задача. Докажите, что точки ветви гиперболы x2/a2 − y2/b2 = 1 при x > 0 допускают параметрическое представление x = a ch (φ), y = b sh (φ), −∞ < φ < +∞. По определению,

ch (φ) =

1

(eφ + e−φ), sh (φ) =

1

(eφ − e−φ).

2

 

2

130

Лекция 19

 

 

19.7Парабола

Пусть уравнение f(x, y) имеет вид (2). Можно считать, что λ2 > 0 è b < 0 (этого всегда можно добиться умножением уравнения на (−1) и дополнительным поворотом

системы координат на угол π). Уберем штрихи, рассматривая новую систему в качестве исходной. Положив p = −b/λ2, получаем уравнение

y2 = 2px, p > 0.

(20)

Определение 1. Множество точек (x, y), удовлетворяющих уравнению (20), называ-

åòñÿ параболой с фокальным параметром p.

Точка F = (p/2, 0) называется фокусом параболы (20). Прямая l : x = −p/2 называется директрисой параболы (20).

Пусть M = (x, y) произвольная точка параболы. Расстояние от нее до фокуса имеет вид

|MF | = p

(x − p/2)2 + y2

=

p

x2 − px + (p/2)2 + 2px

p

=x2 + 2x(p/2) + (p/2)2 = |x + p/2|.

Èòàê, расстояние от любой точки параболы до фокуса |F | равно расстоянию от нее до директрисы l.

Определение 2. Множество тех и только тех точек, для которых расстояние до заданной точки F равно расстоянию до заданной прямой l, называется параболой.

Путь расстояние от заданной точки до заданной прямой равно p. Выберем систему координат таким образом, что F = (p/2, 0) è l : x = −p/2. Åñëè |MF | = |x + p/2|, то, возводя это равенство в квадрат, получаем (20). Таким образом, определения 1 и 2 действительно эквивалентны.

Задача. Написать общее уравнение прямой, проходящей через точку M(x0, y0) параболы y2 = 2px и имеющей с ней единственную общую точку (такая прямая называется касательной к параболе). Доказать, что прямая, ортогональная данной прямой и проходящая через точку M, делит пополам

угол между прямой F M, где точка F фокус параболы, и прямой, параллельной оси x и проходящей через точку M.

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

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

Задача. Кривая S это эллипс, одна из ветвей гиперболы или парабола, E произвольная фиксированная точка на S. Для произвольных точек A è B íà S прямая l(A, B) определяется как прямая AB в случае A 6= B и как касательная в точке A в случае A = B. Докажите, что прямая, проведенная через точку E параллельно l(A, B), имеет с S не более одной общей точки C, помимо

E. Пусть C = E, если общая точка только одна. Таким образом, любой паре точек A, B ставится в

соответствие точка C назовем ее суммой точек A è B. Докажите, что множество S относительно этой операции является абелевой группой.

Лекция 20

20.1 Квадратичные многочлены от трех переменных

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

f(x, y, z) = a11x2 + 2a12xy + 2a13xz + a22y2 + 2a23yz + a33z2 + 2a14x + 2a24y + 2a34z + a44

от декартовых координат x, y, z в геометрическом пространстве и исследуем множест-

во решений уравнения f(x, y, z) = 0 его принято называть поверхностью второго

порядка.

 

 

 

 

 

 

 

 

 

 

Легко проверить, что

 

a12

 

 

 

y

 

 

f(x, y, z) = x y z 1

a22

a23

a24

,

 

 

 

a11

a12

a13

a14

 

x

 

a14 a24 a34

a44

1

 

 

 

 

 

a13

a23

a33

a34

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

а квадратичная часть многочлена f(x, y, z) имеет вид

y .

 

 

f2(x, y, z) =

x y z

a12

a22

a23

 

 

 

 

 

a11

a12

a13

 

x

 

 

 

 

a13 a23 a33

z

 

 

Как и в случае двух переменных, попробуем перейти к более удобной декартовой системе координат.

20.2Декартовы системы и ортогональные матрицы

Пусть e1, e2, e3 è e1, e2, e3 базисные векторы двух декартовых систем коорди-

нат с общим началом. Выразим векторы второй системы в виде линейных комбинаций векторов первой системы

e1

= p11e1 + p21e2

e2

= p12e1 + p22e2

e3

= p13e1 + p23e2

+p31e3,

+p32e3,

+p33e3

и заметим, что

p12

p22

p32

p21

p22

p23

=

(e2

, e1) (e2

, e2) (e2

, e3)

=

0

1

0 .

p11

p21

p31

p11

p12

p13

 

(e1

, e1)

(e1

, e2)

(e1

, e3)

 

1

0

0

p13

p23

p33

p31

p32

p33

(e3, e1)

(e3, e2)

(e3, e3) 0

0

1

131

132

Лекция 20

 

 

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

P >P = I.

( )

Ясно и то, что если матрица перехода обладает свойством ( ), то декартова система переходит в декартову.

Определение. Квадратная вещественная матрица P , удовлетворяющая равенству ( ),

называется ортогональной.

Ортогональные матрицы порядка 2 осуществляют переход между базисами декар-

товых систем на плоскости. Таковы, в частности, матрицы перехода, реализующие поворот (см. лекцию 18).

Данное нами определение применимо и для матриц, порядок которых больше 3.

Согласно ( ), для ортогональных матриц произвольного порядка обращение сводится

к транспонированию:

P −1 = P >.

Кроме того, произведение двух ортогональных матриц остается ортогональной матри-

öåé: åñëè P >P = Q>Q = I, òî (P Q)>(P Q) = Q>(P >P )Q = Q> I Q = Q>Q = I.

Очевидно, единичная матрица I является ортогональной матрицей.

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

Утверждение. Пусть y = P x, где P произвольная ортогональная матрица порядка n и x Rn×1. Тогда сумма квадратов элементов матрицы-столбца y равна сумме квадратов элементов матрицы-столбца x.

Доказательство.

y12 + . . . + yn2 = y>y = (P x)>(P x) = x>(P >P ) x = x>x = x21 + . . . + x2n. 2

Следствие. Пусть B = P AQ, где P, Q произвольные ортогональные матрицы порядка n и A произвольная вещественная матрица порядка n. Тогда сумма квадратов элементов матрицы B равна сумме квадратов элементов матрицы A.

20.3Метод вращений

Попробум упростить квадратичную часть f2(x, y), используя ту же идею поворота сис-

темы координат, как и в случае плоскости. Теперь, однако, у нас есть три координатных плоскости, порождаемые тремя парами координатных осей. Цель вращения получить

нуль вместо какой-нибудь одной пары элементов aij

= aji

ïðè i 6= j. Рассмотрим три

возможности:

 

 

 

a12

 

 

 

 

 

 

sin φ

cos φ

 

 

=

 

 

 

 

 

 

 

,

 

− sin φ

cos φ

 

0

a22

a23

 

0

0

a22

a23

(1)

 

cos φ

sin φ

 

0

a11

a12

a13

 

cos φ

 

sin φ

0

 

 

a11

0

a13

 

 

 

 

0

 

 

0

 

1

a13

 

23

 

33

 

0

 

 

0

 

 

 

1

 

 

e13

a23

e33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

e

 

 

 

 

 

 

 

 

 

 

a12

a

 

 

a

 

 

 

 

 

 

 

 

 

 

 

=

 

a

 

 

 

a

 

,

 

 

0

φ

1

 

0

φ

a22

a23

0

φ

1

 

0

φ

a12

a22

a23

(2)

 

cos

0

sin

a

11

a

12

a

13

 

cos

0

 

 

sin

 

 

a

11

a

12

e0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e e e

 

 

 

sin φ

0

cos φ

a13

a23

a33

 

sin φ

0

 

cos φ

 

 

e0

a

 

a

33

 

 

 

 

 

 

 

 

e23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

e

 

 

Е. Е. Тыртышников

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

133

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

cos φ

 

sin φ

a12

a22

a23

 

0

cos φ

− sin φ

=

a12

a22

0

 

. (3)

 

1

0

 

 

0

 

 

a11

a12

a13

 

1

0

 

0

 

a11

a12

a13

 

 

0

− sin

 

 

cos

 

 

13

23

 

33

0

sin

 

cos φ

e13

e0

e33

 

Обозначим через

 

 

a

è

a

a

 

 

 

 

φ

 

 

e

внедиагональныхe

 

 

 

φ

 

φ

 

 

 

 

 

 

 

 

 

a

 

a

 

 

 

 

 

 

 

 

 

d0, h0

 

d1, h1 суммы квадратов диагональныхeè

 

e

 

 

элементов исходной и новой матриц в каждом из трех случаев:

 

 

 

 

 

 

 

 

 

 

 

 

d0 = a112 + a222 + a332 , h0 = 2a122 + 2a132 + 2a232 ,

 

 

 

 

 

 

 

 

 

 

 

d1 = a112

+ a222

+ a332 , h1

= 2a122

+ 2a132

+ 2a232 .

 

 

 

 

Согласно отмеченным

 

e

e

 

e

 

 

 

e

e

e

 

 

 

 

 

 

 

 

 

 

 

 

 

выше свойствам ортогональных матриц,

 

 

 

 

 

 

 

 

 

 

 

 

d1 + h1 = d0 + h0

h1 = h0 − (d1 − d0).

 

 

 

 

По той же причине в случае (1) имеем a112 + a222 + 2a122 = a112 + a222

и, поскольку a33

= a33,

d1 − d0 = 2a122

 

 

i, j определяют координатную

 

e

e

 

 

 

 

e

. В случае (2) d1 − d0 = 2a132 ,

а в случае (3)

d1 − d0 = 2a232 .

 

 

 

Пусть индексы

 

 

 

 

 

 

 

 

 

 

 

плоскость, в которой проводится вра-

щение (и указывают на то, какое из соотношений (1), (2) èëè (3) имеет место). Выберем

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

d1 − d0 = 2aij2 ≥ 2(h0/6) = h0/3

h1

2

h0.

 

3

Пусть A0 = [aij] è A1 = [eaij]. Рассматривая A1 в качестве новой исходной матрицы, выберем в ней максимальный по модулю внедиагональный элемент и, занулив его с

помощью вращения, получим матрицу A2. Продолжая действовать таким же образом и далее, построим последовательность матриц Ak = [a(ijk)], k = 0, 1, . . . .

Пусть hk обозначает сумму квадратов внедиагональных элементов матрицы Ak. Òîã-

äà

 

 

 

hk

2

 

k

 

h0 → 0 ïðè k → ∞.

3

Следовательно, при любых фиксированных i 6= j последовательность внедиагональных элементов a(ijk) сходится к нулю при k → ∞.

20.4Вложенные подпоследовательности

Лемма об ограниченных последовательностях. Пусть имеется конечное число

ограниченных последовательностей {s(1k)}, . . . , {s(mk)}, k = 1, 2, . . . . Тогда можно

выбрать последовательность номеров k1

< k2

< . . .

таким образом, что каждая

из подпоследовательностей {s1(kl)}, . . . ,

{sm(kl)},

l = 1,

2, . . . , будет сходящейся.

Доказательство. Из ограниченной последовательности {sk} выбираем сходящуюся

подпоследовательность skl и вместо исходных последовательностей рассматриваем под-

последовательности {s(1kl)}, . . . , {s(mkl)}, l = 1, 2, . . . . Они остаются, конечно, огра-

ниченными и при этом первая из них будет сходящейся. Теперь уже из ограниченной последовательности {sk2l } выберем сходящуюся подпоследовательность (подпоследова-

тельность подпоследовательности по отношению к исходной последовательности) и переходим к подпоследовательностям {s(1kli )}, . . . , {s(mkli )}, i = 1, 2, . . . . Полученные вложенные подпоследовательности будут, по-прежнему, ограниченными, но сходящи-

мися являются уже первые две. И так далее. 2

будет сходящейся. Для упрощения
Ak и есть та самая подпос-

134

Лекция 20

 

 

20.5Диагонализация в пределе

Вернемся к методу вращений. Будут ли сходиться к конечным пределам последовательности диагональных элементов a(iik) для нашей ближайшей цели не очень важно.

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

подпоследовательность матриц Ak, в которой каждая из последовательностей диагональных элементов сходится к какому-то пределу.

Чтобы не загромождать обозначения, будем считать, что

ледовательность, для которой все последовательности a(ijk) являются сходящимися (как мы знаем, при i 6= j к нулю). Пусть

lim a(iik) = λi, i = 1, 2, 3.

k→∞

Понятно, что

 

Ak = Pk>A0Pk, k = 1, 2, . . . ,

(#)

где матрицы Pk = [p(ijk)] являются произведениями использованных матриц вращения

(из соотношений (1), (2) èëè (3)). Поэтому при любом k матрица Pk является ортого- нальной (как произведение ортогональных матриц). Следовательно, сумма квадратов всех элементов матрицы Pk при любом k одинакова (и равна 3). Значит, каждая последовательность p(ijk) является ограниченной при k → ∞ и поэтому обладает сходящейся

подпоследовательностью.

По лемме об ограниченных последовательностях, существует подпоследовательность матриц Pk, в которой каждая последовательность p(ijk)

обозначений будем считать, что Pk и есть именно такая подпоследовательность. Пусть

lim p(ijk) = pij, i, j = 1, 2, 3.

k→∞

При каждом k выполняется равенство (#). Переходя к пределу в соответствующих поэлементных равенствах, получаем

 

 

λ1

0

0

 

 

P >A

 

Λ ≡

0

λ

0

=

P.

0

02

λ3

0

 

Кроме того, для каждого k имеем P >P

k

= I

 

P >P = I. В итоге доказана следующая

k

 

 

важная

 

 

 

 

Теорема. Для любой вещественной симметричной матрицы A порядка 3 существуют ортогональная матрица P и диагональная матрица Λ такие, что

Λ = P >AP.

Следствие. Существует декартова система координат, в которой уравнение поверхности второго порядка имеет вид

2

2

2

f(x,e y,e ze) = λ1xe + λ2ye

+ λ3ze + 2b1xe + 2b2ye+ 2b3ze+ a = 0.

Замечание. В силу ортогональности матрицы P , равенство Λ = P >AP выполняется

в том и только том случае, когда

AP = P Λ. Последнее означает, что j-й столбец pj

путем умножения слева и справа на
Веществен-

Е. Е. Тыртышников

135

 

 

матрицы P удовлетворяет уравнению Apj = λjpj (A − λjI)pj

= 0. Учитывая,

÷òî pj 6= 0, находим det(A − λjI) = 0. Таким образом, числа λ1, λ2

, λ3 это корни

кубического многочлена f(λ) = det(A − λI). Если они уже найдены, то столбец pj можно получить как решение однородной системы линейных алгебраических уравнений

(A − λjI)pj = 0.

20.6Диагонализация вещественных симметричных матриц

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

Теорема о диагонализации вещественных симметричных матриц.

ная симметричная матрица A произвольного порядка n приводится к диагональной матрице Λ с помощью некоторой ортогональной матрицы P :

Λ = P >AP.

Доказательство. Начиная с A0 = A, построим последовательность матриц Ak k = 0, 1, . . . , в которой Ak получается из Ak−1

матрицы вращения:

Ak = Rk>Ak−1Rk,

ãäå Rk отличается от единичной матрицы I лишь четырьмя элементами подматрицы, раположенной на пересечении строк и столбцов с номерами i

равной

sin φ

cos φ .

 

 

cos φ

sin φ

= [a(ijk)],

( )

2 × 2-

< j è

Матрица Rk осуществляет поворот на угол φ в координатной плоскости, определяемой номерами i 6= j. Любая матрица вращения такого вида является, очевидно, ортогональ-

íîé.

Ясно, что симметричность матрицы A0 = A наследуется всеми матрицами Ak. Èç

предыдущих исследований мы уже знаем, что φ можно выбрать таким образом, что a(ijk) = a(jik) = 0. Обозначим через dk è hk суммы квадратов диагональных и внедиаго- нальных элементов матрицы Ak. Тогда

2 2

a(iik) + a(jjk)

2

=a(iik−1)

2

+a(jjk−1)

2

+2 a(ijk−1)

2

dk − dk−1 = 2 a(ijk−1) .

Для каждого k будем выбирать плоскость вращения (номера i < j) такми образом, чтобы исключаемый элемент a(ijk−1) был максимальным по модулю среди всех внедиа- гональных элементов матрицы Ak−1. Общее число внедиагональных элементов равно n2 − n. Поэтому

aij

1)

 

2

n2 − n.

(k

 

 

 

hk−1

Отсюда, учитывая равенство dk + hk = dk−1 + hk−1, получаем

hk

hk−1

 

2

 

hk−1

1 −

 

2

 

k

h0 → 0 ïðè k → ∞.

n2

n

n2

n

 

 

 

 

 

 

 

 

 

 

 

 

Из соотношений ( ) вытекает, что

Ak = Pk>APk, k = 1, 2, . . . ,

Для любых фиксированных

136 Лекция 20

ãäå äëÿ âñåõ k матрицы Pk = [p(ijk)] являются ортогональными (как произведения ортогональных матриц).

i, j последовательности a(ijk), p(ijk) являются ограничен-

ными. По лемме об ограниченных последовательностях, существует последовательность

номеров k1 < k2 < . . . такая, что каждая из подпоследовательностей aij(kl)

, pij(kl)

будет

сходящейся. Заметим, что aij(kl) → 0 ïðè i 6= j. Пусть

 

 

lim aii(kl) = λi,

i = 1,

. . . , n,

lim pij(kl) = pij, i, j = 1, . . . , n.

 

l→∞

 

 

l→∞

 

 

 

Введем матрицы

Λ = λ1 ...

,

P = [pij].

 

 

 

 

λn

 

 

 

 

 

 

 

 

 

 

 

Äëÿ âñåõ l = 1, 2, . . . имеем Akl = Pk>l APkl . Переходя к пределу в поэлементных равенствах, получаем

Λ = P >AP.

Из условий ортогональности Pk>l Pkl = I вытекает, что в пределе P >P = I. Значит, матрица P является ортогональной. 2

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

приближенного вычисления матриц Λ è P . Это один из ранних практических методов вычислительной алгебры, предложенный К. Якоби в 1846 году. 1

Задача. Дана симметричная матрица A Rn×n с ненулевой суммой элементов главной диагонали. Доказать существование ортогональной матрицы Q Rn×n такой, что в матрице Q>AQ все элементы главной диагонали одинаковы.

1Последние результаты по изучению метода вращений принадлежат совсем недавнему прошлому: в 1990-х годах были обнаружены его особые возможности, связанные с высокоточным вычислением

малых по модулю элементов матрицы Λ.

Соседние файлы в папке лИТЕРАТУРА