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

linalgebra

.pdf
Скачиваний:
37
Добавлен:
13.03.2015
Размер:
318.02 Кб
Скачать

ЛИНЕЙНАЯ АЛГЕБРА

41

В примерах 13.2(1) и 13.2(3) квадратичные формы в соответствующих базисах имеют канонический вид.

Возникает задача найти базис, в котором данная квадратичная форма имеет канони- ческий вид. Такую задачу называют задачей о приведении квадратичной формы к каноническому виду. Мы укажем алгоритм для решения этой задачи (в России за этим

алгоритмом закрепилось название метод Лагранжа).

Метод Лагранжа приведения квадðатичных форм к каноническому виду.

Пусть в пространстве V с базисом e1, . . . , en задана квадратичная форма q. Таким образом, в координатах x1, . . . , xn, соответствующих базису e1, . . . , en, форма q имеет вид

(13.1). При реализации метода Лагранжа работают не с базисами, а с координатами, которые им соответствуют. Приведем пример реализации метода Лагранжа.

Пример 13.3. Квадратичная форма

q = x21 + 4x1x2 + x3x4.

1-ый шаг. Выделим полный квадрат по x1

x21 + 4x1x2 + x3x4 = x21 + 4x1x2 + 4x22 4x22 + x3x4 = (x1 + 2x2)2 4x22 + x3x4

и сделаем соответствующую замену переменных

 

 

 

 

 

x1

+ 2x2 = y1

 

 

 

 

 

 

 

 

 

 

 

 

 

= y2

 

 

 

 

x2

 

 

 

 

x3 = y3

Получим

2

2

.

 

 

 

 

 

 

 

 

 

 

 

 

x4 = y4.

 

q = y1 4y2 + y3y4

 

 

 

 

2-ой шаг. Выделим полный квадрат по y2

y12 4y22 + y3y4 = y12 (2y2)2 + y3y4

и cделаем замену переменных

 

 

 

 

 

 

 

 

 

y1 = z1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2y2 = z2

 

 

 

 

 

 

 

 

y3 = z3

Получим

 

 

 

 

 

 

 

 

 

 

q = z

2

 

z

2

+ z3z4.

 

y4 = z4.

 

 

 

 

 

 

 

 

1

 

2

 

2

 

2

3-èé шаг. Коэффициенты при z3 è ïðè z4 равны нулю и выделить полный квадрат по

z3 è ïî z4 невозможно. В этом случае делаем замену

 

 

 

 

 

 

 

 

 

= u1

 

 

 

 

 

 

 

 

z1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= u2

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

z3 = u3 + u4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= u3 − u4

 

 

 

 

 

 

 

 

z4

èполучаем q = u21 − u22 + u23 − u24. Форма q приведена к каноническому виду.

Закон инерции. Пусть q - квадратичная форма в пространстве V и пусть

 

1, . . . ,

 

n

è

f1, . . . , fn

e

e

42

ЛИНЕЙНАЯ АЛГЕБРА

- базисы в котором форма q имеет канонический вид. Пусть

 

x1, . . . , xn

è

y1, . . . , yn

- соответствующие базисы. Переупорядочив (если надо) эти базисы получим q(x1, . . . , xn) = x21 + . . . + x2p − x2p+1 − . . . − xp+q,

q(y1, . . . , yn) = y12 + . . . + yr2 − yr2+1 − . . . − yr+s,

ãäå p, q, r, s - целые неотрицательные числа.

Теорема 13.4 (Закон инерции). В описанной выше ситуации имеем:

p = r è q = s.

Определение 13.5. Квадратичная форма q называется положительно определенной, если q(x1, . . . , xn)>0 для любых x1, . . . , xn причем q(x1, . . . , xn) = 0 тогда и только тогда,

когда x1 = . . . = xn = 0. Квадратичная форма q называется отрицательно определенной, если q(x1, . . . , xn)60 для любых x1, . . . , xn причем q(x1, . . . , xn) = 0 тогда и только тогда,

когда x1 = . . . = xn = 0.

ßñíî, ÷òî â n-мерном пространстве квадратичная форма

положительно определена тогда и только тогда, когда она приводится к канони-

ческому виду

x21 + . . . + x2n,

отрицательно определена тогда и только тогда, когда она приводится к канониче-

скому виду

−x21 − . . . − x2n.

Имеется критерий Сильвестра, который позволяет для данной квадратичной формы выяснить является ли она положительно или отрицательно определенной.

Теорема 13.6 (критерий Сильвестра). Пусть q = q(x1, . . . , xn) - квадратичная форма и

 

q11

q12

. . . q1n

 

 

.

.

... .

Q =

q21

q22

. . .

q2n

 

 

 

 

 

 

 

 

qn1

qn2

. . .

qnn

- матрица формы q в каком-нибудь базисе. Рассмотрим угловые миноры

2

=

 

 

 

q

 

 

,

 

,

 

 

∆ =

 

 

 

 

 

 

 

 

1

 

q11

11q12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q22

 

 

 

 

 

 

 

q21

 

 

 

 

 

 

q11

 

q12

 

q13

 

 

 

 

 

 

 

q22

 

 

 

 

 

 

3 = q21

 

 

q23 ,

 

 

 

 

 

 

q32

 

q33

 

 

q31

 

 

 

..........................

.

Тогда

(1)форма q положительно определена тогда и только тогда, когда все угловые миноры не равны нулю, причем

1 > 0, 2 > 0, 3 > 0, 4 > 0, . . .

ЛИНЕЙНАЯ АЛГЕБРА

43

(2)форма q отрицательно определена тогда и только тогда, когда все угловые миноры не равны нулю, причем

1 < 0, 2 > 0, 3 < 0, 4 > 0, . . .

Пример 13.7. Рассмотрим квадратичную форму q = −x21 4x1x2 5x22 + 2x2x3 2x23. Матрицей этой формы является

 

 

1

 

2

 

0

 

 

 

 

0

 

 

 

1

 

 

2

Q = 2

 

5

 

1

.

Вычисляем:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

1 = 1

=

 

 

1 < 0,

 

 

 

 

 

 

 

 

2 =

2

 

 

5

= 1 > 0,

 

1

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 = 2 5 1

= 1 < 0.

По критерию Сильвестра отсюда следует,

что форма

q отрицательно определена.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 13.8. Рассмотрим квадратичную форму q = −x21 4x1x2. Матрицей этой формы

является

 

 

 

(2

0

)

 

 

 

 

 

Q =

1

 

2 .

Вычисляем:

1 = 1

 

 

 

 

 

 

2

 

1 < 0,

 

 

2

1

=

 

 

 

 

0

 

 

 

2

=

 

 

 

= 4 < 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По критерию Сильвестра отсюда следует, что форма q не является положительно определенной и не является отрицательно определенной.

14. Квадратичные формы в евклидовом пространстве

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

горитм для нахождения канонического ортонормированного базиса для квадратичной формы в Rn.

Итак, пусть q - квадратичная форма в Rn,

 

q11

q12

. . . q1n

 

 

.

.

... .

Q =

q21

q22

. . .

q2n

 

 

 

 

 

 

 

 

qn1

qn2

. . .

qnn

- матрица квадратичной формы q в стандартном базисе пространства Rn.

Лемма 14.1.

44

ЛИНЕЙНАЯ АЛГЕБРА

Найдем характеристический многочлен

q11 − λ

q21

χQ(λ) := det(Q − λE) =

qn.1

матрицы Q.

 

..

 

 

 

q12

. . .

q1n

 

.

 

.

.

q22 − λ . . .

q2n

 

 

 

 

 

 

 

 

 

 

 

qn2

. . .

qnn

 

λ

Для второго шага важен следующий факт.

В описанной выше ситуации имеем: характеристический многочлен χQ(λ) имеет n вещественных корней (с учетом кратностей).

Вторым шагом мы находим корни характеристического многочлена χQ(λ). Допустим этими корнями являются λ1, . . . , λm, ãäå λi ≠ λj ïðè i ≠ j.

Для каждого корня λi находим ортонормированный базис пространства решений системы однородных линейных уравнений

 

q11 − λi

q12

 

. . .

q1n

x1

 

 

0

.

. ... .

.

 

.

 

q21

q22

λi

. . . q2n

x2

 

=

0 .

 

qn1

qn2

 

. . .

qnn

 

 

 

 

 

 

 

λi xn

 

 

0

Образуем систему векторов объединив базисные векторы полученные на 3-м шаге для всех корней λi. Эта система векторов и будет искомым базисом пространства Rn. Матрицей формы q в этом базисе будет диагональная матрица, у которой на диагонали

расположены корни λi, причем каждый корень встречается столько раз, какова кратность этого корня.

Пример 14.2. Пусть q = 3x12 4x1x2. Матрицей этой формы является

(2

0)

3

2

Найдем канонический ортонормированный базис для формы q. 1-й шаг. Находим характеристический многочлем матрицы Q:

 

 

3 − λ

 

 

 

 

 

 

χQ(λ) =

 

 

2

 

= λ2

 

3λ 4.

2-й шаг. Находим корни характеристического

многочлена:

 

 

2

 

0

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ1 = 1,

 

 

 

λ2 = 4.

 

3-й шаг. Для корня λ1 = 1 имеем систему

 

 

 

 

 

3

(

1)

 

 

2

 

 

 

 

 

x1

 

0

(

2

 

 

 

 

0 (1))(x2) = (0)

èëè

(2 1)(x2)

=

(0).

 

 

 

 

 

4

2

x1

 

 

 

0

 

 

Решая эту систему находим нормированное решение

 

 

 

 

 

 

 

1

 

5

(

2

)

 

 

 

 

 

 

 

 

 

1

 

 

 

1 .

 

 

 

 

 

 

f

 

=

 

 

 

 

 

 

 

 

ЛИНЕЙНАЯ АЛГЕБРА

45

Для корня λ1 = 4 имеем систему

 

 

0 4)(x2) (0)

 

(

2

 

 

3

4

2

 

x1

=

0

èëè

( 2 4)(x2) (0)

 

 

 

1

2

 

x1 =

0 .

Решая эту систему находим нормированное решение

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

f

2 =

 

(1).

 

 

 

 

5

 

4-é шаг. Базис f1, f2 - искомый. В этом базисе матрица формы q имеет вид

()

1

0 .

0

4

Дополнение 1. Формулы перехода

Пусть V векторное пространство размерности n и пусть e1, e2, . . . , en è f1, f2, . . . , fn

базисы пространства V . Разложим векторы второго базиса по векторам первого базиса:

 

 

l

 

 

n

(14.1)

fk = alk

e

l, k = 1, . . . , n

 

=1

 

 

èиз коэффициентов этих разложений образуем матрицу

 

a11

a12

. . .

a1n

 

 

a21

a22

. . .

a2n

 

 

 

 

 

 

 

A = a. .n.1. . a. .n.2. .

...... . .a.nn. .

(коэффициенты разложения вектора fk образуют k-ый столбец матрицы A). Матрицу A

называют матрицей перехода от базиса e1, . . . , en к базису f1, . . . , fn. Пусть

 

a11

a12

. . .

a1n

 

A1 =

a21

a22

. . .

a2n

 

 

 

 

 

 

 

a. .. . . a. .. . .

...... . .a.. .

 

n1

n2

 

nn

обратная к A матрица.

Лемма 14.3. Матрицей перехода от базиса f1, . . . , fn к базису e1, . . . , en является матрица A1, òî åñòü,

 

 

i

 

 

n

(14.2)

e

m = aim

f

i, m = 1, . . . , n.

 

=1

 

 

46 ЛИНЕЙНАЯ АЛГЕБРА

Лемма 14.4. Пусть x1, . . . , xn координаты вектора v V в базисе e1, . . . , en è y1, . . . , yn

координаты этого же вектора v

в базисе f1, . . . , fn. Тогда

 

i

 

 

 

n

 

(14.3)

yk =

akixi,

k = 1, . . . , n.

 

=1

 

 

i

 

 

 

n

 

(14.4)

xk =

akiyi,

k = 1, . . . , n.

 

 

=1

 

Эти формулы в матричном виде записываются следующим образом:

(14.5)

y = xA,

(14.6)

x = yA,

ãäå

 

 

 

x = (x1, . . . , xn),

y = (y1, . . . , yn).

Рассмотрим оператор

φ : V → V.

Лемма 14.5. Пусть Φ матрица оператора φ в базисе e1, . . . , en и Ψ матрица оператора φ в базисе f1, . . . , fn. Тогда

(14.7) Ψ = A1ΦA.

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

Φ =

φ11

φ12

. . .

φ1n

 

,

φ21

φ22

. . .

φ2n

 

 

 

 

 

 

 

 

 

φ. .n.2. .

 

 

 

 

 

φ. .n.1. .

.. .. .. . .φ.nn. .

 

Тогда

 

 

 

n

 

 

 

(14.8)

φ(

 

l) =

φml

 

m,

e

e

 

 

 

m=1

 

 

 

 

 

 

 

 

 

Ψ =

ψ11

ψ12

. . .

ψ1n

 

ψ21

ψ22

. . .

ψ2n

 

 

 

 

 

 

 

 

ψn.2

 

 

 

 

ψ. .n.1. .

.. .. .. .

.ψ.nn. .

l = 1, . . . , n,

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(14.9)

 

 

 

 

φ(fk) =

 

ψikfi,

k = 1, . . . , n.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя (14.1), (14.8) и (14.2) имеем:

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ(

 

k) = φ ( n

alk

 

l)

 

n

 

 

 

 

n

 

( n

φml

 

m) =

f

e

=

 

 

 

alk φ(

e

l) =

 

alk

e

 

 

l=1

 

 

 

l=1

 

 

 

 

=1

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

m ( n

aim

 

i) =

n

alkφmlaim

 

i.

 

 

 

alkφml

e

m =

 

 

 

 

 

alkφml

e

f

 

f

l;m=1

 

 

 

l;m=1

 

 

 

 

 

 

 

i=1

 

 

 

l;m;i=1

 

 

 

 

 

Сравнивая это с (14.9) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψik =

 

alkφmlaim,

i, k = 1, . . . , n.

 

 

 

 

 

 

 

 

 

 

 

 

l;m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нетрудно заметить, что в матричном виде эти равенства записываются равенством (14.7).

ЛИНЕЙНАЯ АЛГЕБРА

47

Матрицы Φ и Ψ называют подобными, если существует матрица A такая, что

(14.10)

Ψ = A1ΦA.

Лемма 14.6. Характеристические многочлены подобных матриц совпадают.

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

χ (λ) = det(Ψ − λ) = det(A1ΦA − λ) = det(A1− λ)A) = det(A1) det(Φ − λ) det(A) = det(Φ − λ) = χ (λ).

Характеристическим многочленом оператора φ называют

χ (λ) := det(Φ − λE),

где Φ матрица оператора φ в базисе e1, . . . , en. Из лемм 14.5 и 14.6 следует, что харак-

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

Рассмотрим квадратичную форму

q : V → R.

Лемма 14.7. Пусть Q матрица квадратичной формы q в базисе e1, . . . , en è R матрица квадратичной формы q в базисе f1, . . . , fn. Тогда

R = AQA .

Доказательство. Пусть x1, . . . , xn координаты вектора v V в базисе e1, . . . , en è

y1, . . . , yn координаты этого же вектора v в базисе f1, . . . , fn. Рассмотрим строки x = (x1, . . . , xn), y = (y1, . . . , yn)

С одной стороны,

q(v) = yRy ,

 

а с другой стороны, учитывая (14.6),

 

q(v) = xQx = (yA)Q(yA) = y(AQA )y

 

и, следовательно, R = AQA .

 

Дополнение 2. Ортогональные матрицы

Матрица A называется ортогональной, если ее строки образуют ортонормированную систему векторов. Нетрудно заметить, что матрица A ортогональна тогда и только тогда, когда она удовлетворяет матричному уравнению

(14.11) A · A = E.

Лемма 14.8. Определитель ортогональной матрицы равен ±1.

Доказательство. Возьмем определители правой и левой частей равенства (14.11) и получим

(14.12)

det(A · A ) = det(E).

Учитывая, что

det(E) = 1,

det(A · A ) = det(A) det(A ) (определитель произведения матриц равен произведению определителей),

48

ЛИНЕЙНАЯ АЛГЕБРА

det(A ) = det(A) (при транспонировании матрицы ее определитель не меняется) из (14.12) получаем

 

det(A)2 = 1

откуда и следует лемма.

 

Следствие 14.9. Ортогональная матрица невырождена.

Дополнение 3. Словарь некоторых терминов

Система линейных уравнений [system of linear equations].

Совместная (несовместная) система линейных уравнений [consistent (inconsistent) system of linear equations].

Однородная (неоднородная) система линейных уравнений [homogeneous (inhomogeneous) system of linear equations].

Однородная система B, соответствующая системе A [homogeneous system B associated to system A].

Частное решение [particular solution].

Элементарные преобразования системы по строкам [elementary row operations]. Ступенчатый вид [echelon form].

Приведенный ступенчатый вид [reduced echelon form]. Главное неизвестное [pivot unknown].

Свободное неизвестное [free unknown].

Общее решение системы линейных уравнений [не имеет аналога на английском]. Общее решение это формулы, выражающие главные неизвестные через свободные.

Векторное пространство [vector space]. Нередко используют неправильное название "линейное пространство".

Линейное подпространство [subspace].

n-мерное координатное пространство [ n-dimensional coordinate space]. Нередко используют неправильно название "n-мерное арифметическое пространство".

Система векторов [list of vectors].

Фундаментальная система решений [не имеет аналога на английском]. Фундаментальная система решений это базис пространства решений системы однородных линейных уравнений.

Скалярное произведение [inner product]

Евклидово пространство = пространство со скалярным произведением [inner product

space]

Стандартное скалярное произведение в Rn [scalar product, dot product]

Формулы перехода [transition rules] Ортогональная матрица [orthogonal matrix]

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