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

A_G_2014

.pdf
Скачиваний:
43
Добавлен:
10.02.2015
Размер:
2.75 Mб
Скачать

§ 5. Приведенная форма квадратичной функции

271

где

 

(5.25)

a˜i = ai + xˆi0aii,

i = 1, 2, . . . , r,

a˜i = ai, i = r + 1, . . . , n,

(5.26)

 

n

 

a˜0 = a0 +

aixˆi0.

(5.27)

 

=1

 

 

i

 

Вследствие (5.19), (5.26) справедливо равенство Ir+2(BU) = Ir+2(B). Точно так же проверяется, что умножение матрицы вида BU слева на матрицу UT не меняет инварианта Ir+2(BU). Отсюда, очевидно, вытекает равенство (5.24). Второе утверждение леммы с использованием соотношений (5.25)–(5.27) доказывается аналогично. При этом надо учесть, что определитель (5.23) не меняется при переходе от

e

матрицы B к матрице B.

3.4. Теорема. Пусть ранг матрицы A квадратичной функции F равен r, ранг матрицы B не превосходит r + 1. Пусть при помощи замены переменных x = x0 + T y с ортогональной матрицей T квадратичная функция F приведена к виду

ˆ

2

2

2

+ aˆ0.

(5.28)

F (y) = α1y1

+ α2y2

+ · · · + αmym

Тогда: 1) m = r, αi = λi, где λi, i = 1, 2, . . . , r, — все ненулевые собственные числа матрицы A; 2) справедливо равенство

aˆ0 = Ir+1(B)/Ir(A).

(5.29)

Доказательство. Используем формулы (4.4)–(4.8), связывающие исходную и приведенную формы квадратичной функции. В рас-

T ˆ

сматриваемом случае T AT = A = diag(α1, α2, . . . , αm). Поскольку

ˆ

матрица T ортогональна, то матрицы A, A подобны, откуда вытекает справедливость утверждения 1). Используем теперь тот факт, что, в рассматриваемом случае aˆ = 0. Поэтому по лемме 3.2 полу-

 

 

ˆ

· · ·

λraˆ0. Матрица A имеет ровно r нену-

чаем, что Ir+1(B) = λ1λ2

левых собственных чисел, следовательно, Ir(A) = λ1λ2 · · · λr

 

(см.

п. 2, с. 192). Таким образом,

ˆ

T

BQ

aˆ0 = Ir+1(B)/Ir(A). Матрица Q

 

подобна матрице B. Поэтому все их инварианты совпадают. С дру-

гой стороны матрица QT BQ удовлетворяет условиям леммы 3.3, зна-

чит, Ir+1(Q

T

 

ˆ

ˆ

 

 

 

BQ) = Ir+1(B), т. е. Ir+1(B) = Ir+1(B), и утвержде-

ние 2) также доказано.

 

 

 

 

272

Глава 14. Квадратичные формы и квадратичные функции

3.5. Теорема. Пусть ранг матрицы A квадратичной функции F равен r, r 6 n − 1, ранг матрицы B равен r + 2. Пусть при помощи замены переменных x = x0 + T y с ортогональной матрицей T квадратичная функция F приведена к виду

ˆ

2

2

2

+ 2bym+1.

(5.30)

F (y) = α1y1

+ α2y2

+ · · · + αmym

Тогда: 1) m = r, αi = λi, где λi, i = 1, 2, . . . , r, — все ненулевые собственные числа матрицы A; 2) выполнено равенство

b2 = −Ir+2(B)/Ir(A).

(5.31)

Доказательство. Справедливость утверждения 1) обосновывается точно так же, как и при доказательстве теоремы 3.4. В рассматриваемом случае лишь одна компонента вектора aˆ отлична от

нуля. Поэтому по лемме

3.1 получаем, что b

2

 

 

 

ˆ

 

= −Ir+2(B)/Ir(A).

Равенство Ir+2

ˆ

Ir+2

(B)

получаем аналогично

доказательству

(B) =

 

 

 

T

BQ, а так-

предыдущей теоремы, используя подобие матриц B и Q

 

же лемму 3.3.

Глава 15

Кривые второго порядка

§ 1. Приведение к простейшему виду уравнения кривой второго порядка

1. Как и в § 6, гл. 4, будем рассматривать плоскость, отнесенную к декартовой системе координат x1, x2.

Множество всех точек x = (x1, x2) плоскости, удовлетворяющих уравнению

a11x12 + 2a12x1x2 + a22x22 + 2a1x1 + 2a2x2 + a0 = 0,

(1.1)

называют кривой второго порядка. Здесь aij, i, j = 1, 2, ai, i

= 0,

1, 2, — вещественные числа, называемые коэффициентами уравнения. Для сокращения записей, как и в предыдущей главе, введем в

рассмотрение симметричную ненулевую матрицу

()

A =

a11

a12

(1.2)

a12

a22

и вектор a = (a1, a2). Тогда уравнение (1.1) запишется в виде1)

(Ax, x) + 2(a, x) + a0 = 0.

(1.3)

2.В соответствии с общей теорией квадратичных функций

(см. § 4, § 5, гл. 14) упрощение этого уравнения мы будем выполнять с помощью замены переменных

x = x0 + T y,

(1.4)

где T — ортогональная матрица вида

()

T =

cos φ

sin φ .

(1.5)

 

sin φ

cos φ

 

Геометрически эта замена переменных может быть интерпретирована как поворот координатных осей против часовой стрелки на угол φ, если считать при этом, что исходная декартова система координат

1)В этой главе под скалярным произведением всюду понимается стандартное скалярное произведение в пространстве R2.

274

Глава 15. Кривые второго порядка

правая, т. е. поворот от оси x1 к оси x2 — поворот против часовой стрелки (см. п. 2, с. 253), и последующий перенос начала системы координат в точку x0.

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

x = T z

(1.6)

в уравнении (1.3). Проводя элементарные выкладки, получим

 

(T T AT z, z) + 2(ˆa, z) + a0 = 0,

(1.7)

где aˆ = T T a.

Построим ортогональное преобразование T так, чтобы матрица T T AT приняла диагональный вид (см. упражнение 2 на с. 255). С этой целью решим характеристическое уравнение

|A − λI| = 0.

Корни его легко выписываются в явном виде:

 

a11 + a22

 

(a11

a22)2 + 4a2

 

 

λ1,2 =

 

±

2

12

.

(1.8)

 

 

Если a12 = 0, то λ1 = a11, λ2 = a22. Положим в этом случае T = I.

Уравнение (1.7) принимает вид

 

 

 

 

 

λ1z2

+ λ2z2 + 2ˆa1z1 + 2ˆa2z2 + a0 = 0

(1.9)

1

2

 

 

 

 

 

 

(здесь aˆ1 = a1, aˆ2 = a2).

Если a12 = 0, то, очевидно, λ1 = λ2. Найдя λ1, λ2, определим

̸

 

 

 

̸

1

2

:

соответствующие им единичные собственные векторы e

, e

ek = (cos φk, sin φk),

k = 1, 2,

 

 

как решения уравнений

 

 

 

 

 

 

 

 

Aek = λkek, k = 1, 2,

 

 

или, более подробно,

 

 

 

 

 

 

 

 

(a11 − λk) cos φk + a12 sin φk = 0,

 

 

a12 cos φk + (a22 − λk) sin φk = 0,

 

 

откуда получаем уравнения для определения углов φ1, φ2:

 

tg φ

k

=

a11 − λk

, k = 1, 2.

 

(1.10)

 

 

a12

 

 

§ 1. Приведение к простейшему виду уравнения кривой второго порядка

275

Будем считать, что φ1, φ2 (−π/2, π/2). Причем, поскольку собственные векторы симметричной матрицы, отвечающие различным собственным числам ортогональны (см. теорему 8, с. 227), то обязательно

φ1 − φ2 = ±π/2.

Элементарные вычисления дают

 

 

 

 

 

 

 

 

 

 

 

 

 

tg φ

1

tg φ

2

=

λ1 − λ2

=

(a11 − a22)2 + 4a122

.

(1.11)

 

 

 

a12

a12

В соответствии со знаком a12 занумеруем собственные числа (и соответствующие им углы) так, чтобы tg φ1 6 tg φ2, т. е.

φ2 = φ1 + π/2.

Используя общие построения (см. упражнение 2 на с. 255), мат-

рицу T составим из собственных векторов e1, e2:

 

cos φ1

cos φ2

cos φ1

sin φ1

 

T = (sin φ1

sin φ2 ) = (sin φ1

cos φ1).

(1.12)

При указанном выборе матрицы T получаем

 

 

 

λ

0

 

 

 

T T AT = ( 01

λ2),

 

(1.13)

иуравнение (1.7) вновь принимает вид (1.9).

3.Дальнейшие упрощения уравнения (1.3) используют перенос начала системы координат. Будем различать два случая: det(A) ≠ 0

иdet(A) = 0.

3.1.Предположим сначала, что det(A) ≠ 0. Это условие эквива-

лентно тому, что λ1, λ2 ≠ 0. При выполнении этого условия, выделяя полные квадраты, т. е. полагая y = xˆ0 + z, где xˆ01 = aˆ11, xˆ02 = aˆ22,

уравнение (1.9) можно записать в виде

λ1y12 + λ2y22 + aˆ0 = 0.

(1.14)

Здесь aˆ0 = a0 − aˆ211 − aˆ222.

3.2. Пусть теперь det(A) = 0. Напомним, что мы считаем, что A ≠ 0. Понятно, что в этом случае либо λ1 = 0, либо λ2 = 0. Одновременно λ1 и λ2 не могут равняться нулю (почему?). Предположим для определенности, что λ1 ≠ 0. Тогда уравнение (1.9) можно представить в виде

λ1(z1 + aˆ11)2 + 2ˆa2z2 + aˆ0 = 0.

(1.15)

276

Глава 15. Кривые второго порядка

где aˆ0 = a0 − aˆ121.

 

Здесь опять надо различать два случая.

1) Если aˆ2 = 0, положим y = xˆ0 + z, где xˆ01 = aˆ11, xˆ02 = 0. Такая замена переменных приведет уравнение (1.15) к виду

 

λ1y12 + aˆ0 = 0.

(1.16)

2) Если aˆ2 ̸=

0, представим уравнение (1.15) в форме

 

λ1

(z1 + aˆ11)2 + 2ˆa2(z2 + aˆ0/(2ˆa2)) = 0

 

и положим y = xˆ0 +z, xˆ01 = aˆ11, xˆ02 = aˆ0/a2. Тогда уравнение (1.15) примет вид

λ1y12 + 2ˆa2y2 = 0.

(1.17)

Предполагая, что λ1 = 0, а λ2 ̸= 0, можно точно так же преобра-

зовать уравнение (1.9) либо к уравнению

 

λ2y22 + aˆ0 = 0,

(1.18)

либо к уравнению

(1.19)

λ2y22 + 2ˆa1y1 = 0.

4. Таким образом, полагая, что z в (1.6) равно −xˆ0 + y, где в зависимости от свойств собственных чисел матрицы A вектор xˆ0 выбирается по одной из выше приведенных формул, а T — матрица вида (1.5) такова, что выполнено (1.13), получим, что общее уравнение (1.1) кривой второго порядка при помощи замены переменных x = x0 + T y, где x0 = −T xˆ0, можно преобразовать к одной из следующих форм:

λ1y12 + λ2y22 + aˆ0 = 0,

λ1, λ2 ̸= 0,

(1.20)

λ2y22 + 2ˆa1y1

= 0,

λ2

̸= 0, λ1

= 0,

(1.21)

λ2y22 + aˆ01

= 0,

λ2

̸= 0,

λ1

= 0,

(1.22)

λ1y12 + 2ˆa2y2

= 0,

λ1

̸= 0, λ2

= 0,

(1.23)

λ1y12 + aˆ01

= 0,

λ1

̸= 0,

λ2

= 0.

(1.24)

Важно подчеркнуть, что угол φ и вектор x0 определяются по коэффициентам уравнения (1.1) при помощи простых явных формул. В соответствии с общей теорией квадратичных функций каждое уравнение (1.1) при помощи преобразования вида (1.4) может быть однозначно приведено лишь к одной из форм (1.20)–(1.24).

§ 2. Геометрические свойства кривых второго порядка

277

5. Из общих результатов, полученных при изучении квадратичных функций, следует, что коэффициенты уравнений (1.20)–(1.24) однозначно определяются при помощи простых формул по коэфффициентам уравнения (1.1), так что построение преобразования (1.4), фактически, требуется лишь для того, чтобы установить, как располагается исследуемая кривая по отношению к исходной декартовой системе координат.

Наряду с матрицей A, определяемой соотношением (1.2), введем в рассмотрение матрицу

 

a11

a12

a1

 

 

B =

a12

a22

a2

,

 

a1

a2

a0

 

 

соответствующую квадратичной функции, определяемой левой частью уравнения (1.1), и выпишем выражения для коэффициентов

уравнений (1.20)–(1.24) (см. теоремы 3.4, 3.5, с. 271):

aˆ0 = I3(B)/I2(A), aˆ1 = aˆ2 = −I3(B)/I1(A), aˆ01 = I2(B)/I1(A).

Здесь (см. общие формулы для инвариантов оператора, с. 191)

 

 

I

 

a1

a2

 

a0

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11

a12

a1

 

 

 

 

 

 

 

 

 

 

 

 

a11

a12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3(B) = det(B) =

a12

a22

a2

 

,

2(A) = det(A) =

 

a12

a22

 

,

 

 

 

 

 

 

 

 

 

 

 

 

a

 

a

 

 

 

a

 

a

 

 

 

 

 

(A) = tr(A) = a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a

 

,

 

(B) =

 

 

11

1

 

+

 

 

22

 

2

 

.

 

 

I1

 

 

 

11

 

22

 

 

I2

 

a1 a0

 

 

a2 a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По сравнению с (7.5), с. 192, в выражении I2(B) на одно слагаемое меньше, поскольку в рассматриваемом случае det A = 0.

§ 2. Геометрические свойства кривых второго порядка

Опираясь на уравнения (1.20) – (1.24), исследуем геометрические свойства кривых второго порядка.

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

λ1x2 + λ2y2 + d = 0, λ1, λ2 ̸= 0,

(2.1)

y2 = 2px,

(2.2)

y2 + d = 0.

(2.3)

278 Глава 15. Кривые второго порядка

1. Начнем с уравнения (2.3). Возможны три случая: d = 0, кри-

вая совпадает с осью y; d < 0, кривая распадается на две параллель-

ные прямые y =

−d, y = − −d; d > 0, множество точек плоскости,

удовлетворяющих уравнению (2.3), пусто, говорят, что в этом случае кривая распадается на две мнимые параллельные прямые.

2. Исследуем уравнение (2.1). Здесь нужно различать такие случаи:

1)знаки собственных чисел λ1, λ2 матрицы A совпадают, при этом, не ограничивая общности, можно считать, что λ1, λ2 > 0;

2)знаки собственных чисел λ1, λ2 различны.

Кривые, соответствующие первому случаю, называют эллипсами. Здесь опять нужно различать три случая: d = 0, кривая вырождается в точку, совпадающую с началом координат; d > 0, уравнение определят так называемый мнимый эллипс; d < 0, в этом случае уравнение (2.1) запишем в виде

 

 

x2

+

y2

= 1,

 

(2.4)

где

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ1

 

 

λ2

 

a =

−d

,

 

b =

−d

.

 

 

 

 

Кривую, описываемую уравнением (2.4), называют эллипсом. Кривые, соответствующие случаю, когда λ1, λ2 имеют различные

знаки, называют гиперболами. Будем для определенности считать, что λ1 > 0, λ2 < 0 и рассмотрим три случая. Если d = 0, то уравнение (2.1) можно записать в виде

λ1x = ± −λ2y,

т. е. в данном случае кривая распадается на две прямые, пересекающиеся в начале координат. Случаи d < 0, d > 0, фактически, можно не различать, так как они сводятся один к другому за счет переименования осей координат.

Будем для определенности считать, что d < 0. Тогда уравнение (2.1) можно записать в виде

 

 

x2

y2

(2.5)

 

 

 

 

 

= 1,

где

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

a =

−d

,

b =

−d

.

 

λ1

 

−λ2

Кривую, описываемую уравнением (2.5), называют гиперболой.

§ 2. Геометрические свойства кривых второго порядка

279

3. Опишем геометрические свойства эллипса (см. рис. 1). Непосредственно из уравнения (2.4) вытекает, что для всех точек эллипса справедливы неравенства: |x| 6 a, |y| 6 b, т. е. эллипс — ограниченная кривая, расположенная в соответствующем прямоугольнике.

Рис. 1. К описанию геометрических свойств эллипса

Точками пересечения этой кривой с осями координат являются точки (±a, 0), (0, ±b). Они называются вершинами эллипса. Оси координат — оси симметрии эллипса, так как если точка (x, y) принадлежит эллипсу, то точки (−x, y), (x, −y) также лежат на эллипсе. Начало координат — центр симметрии эллипса, так как, если точка (x, y) принадлежит эллипсу, то и точка (−x, −y) лежит на эллипсе.

Числа a, b называют длинами полуосей эллипса. Будем для опре-

деленности считать, что a > b. Понятно, что при a = b эллипс пре-

вращается в окружность (радиуса a). Положим c = a2 − b2. Величи-

на e = c/a = 1 − b2/a2 [0, 1) характеризует степень вытянутости эллипса вдоль большой полуоси и называется эксцентриситетом эллипса.

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

 

 

 

 

(2.6)

(x + c)2 + y2 +

(x − c)2 + y2 = 2a.

 

 

Равенство (2.6) означает, что сумма расстояний от точки эллипса до фокусов одна и та же для всех точек эллипса (см. рис. 2). Это свойство эллипса можно принять за его определение, так как, исходя из (2.6), очевидно, можно получить уравнение эллипса.

Докажем справедливость равенства (2.6) для точек, принадлежащих эллипсу. Используя равенства c2 = a2 − b2, y2 = b2 − b2x2/a2, можем написать

(x + c)2 + y2 = x2 + 2cx + a2 − b2 + b2 − b2x2/a2 =

280

Глава 15. Кривые второго порядка

Рис. 2. К определению эллипса и гиперболы

=x2(1 − b2/a2) + 2cx + a2 = x2c2/a2 + 2cx + a2 =

=(xc/a + a)2. (2.7)

Точно так же

(x − c)2 + y2 = (−xc/a + a)2.

Заметим, что c < a. Учтем также, что |x| 6 a для любой точки эллипса. Поэтому справедливы неравенства

 

xc/a + a > 0,

−xc/a + a > 0,

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

 

 

 

 

 

 

 

 

(x + c)2 + y2 = xc/a + a,

 

(x − c)2 + y2 = −xc/a + a,

откуда непосредственно вытекает (2.6).

4. Опишем геометрические свойства гиперболы (см. рис. 3). Из уравнения (2.5) непосредственно вытекает, что если точка (x, y) лежит на гиперболе, то x2 > a2, y2 6 b2x2/x2, т. е. кривая, описываемая уравнением (2.5), лежит вне полосы |x| < a и внутри соответствующих (вертикальных) углов, образованных прямыми y = ±(b/a)x.

Как и в случае эллипса, проверяется, что кривая симметрична относительно осей координат. Начало координат — центр симметрии кривой. Точки (−a, 0), (a, 0) пересечения с осью x называются вершинами гиперболы.

Прямые y = ±(b/a)x асимптоты соответствующих ветвей гиперболы (рис. 3). Покажем это применительно к ветви, определяемой

уравнением

x2 − a2, x > a,

 

b

(2.8)

y = a

и прямой y = (b/a)x. Для остальных ветвей выкладки полностью аналогичны. В соответствии с определением асимптоты (см. курс математического анализа) достаточно проверить справедливость следующих равенств:

 

ab

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 − a2

 

b

 

 

 

b

 

 

 

 

 

b

x

 

lim

 

,

lim

 

 

x2

 

a2

= 0.

 

x

= a

(a

 

x→∞

 

 

x→∞

 

a

)

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