linalgebra
.pdfЛИНЕЙНАЯ АЛГЕБРА |
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.
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′ |
. . . |
a1′ n |
|
A−1 = |
a21′ |
a22′ |
. . . |
a2′ n |
|
|
|
|
|
|
|
|
a. .′ . . . a. .′ . . . |
...... . .a.′. . |
|||
|
n1 |
n2 |
|
nn |
обратная к A матрица.
Лемма 14.3. Матрицей перехода от базиса f1, . . . , fn к базису e1, . . . , en является матрица A−1, òî åñòü,
|
|
∑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 = |
aki′ xi, |
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) Ψ = A−1Φ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) |
Ψ = A−1ΦA. |
Лемма 14.6. Характеристические многочлены подобных матриц совпадают.
Доказательство. Имеем:
χ (λ) = det(Ψ − λ) = det(A−1ΦA − λ) = det(A−1(Φ − λ)A) = det(A−1) 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]