Теорема о возможности приведения квадратичной формы к каноническому виду. Приведение квадратичной формы к каноническому виду Приведение квадратичной формы к каноническому виду пример
Данный метод состоит в последовательном выделении в квадратичной форме полных квадратов.
Пусть дана квадратичная форма
Напомним, что, ввиду симметричности матрицы
,
Возможны два случая:
1. Хотя бы один из коэффициентовпри квадратах отличен от нуля. Не нарушая общности, будем считать(этого всегда можно добиться соответствующей перенумерацией переменных);
2. Все коэффициенты,
но есть коэффициент , отличный от нуля (для определённости пусть будет).
В первом случае преобразуем квадратичную форму следующим образом:
,
а через обозначены все остальные слагаемые.
представляет собой квадратичную форму от (n-1) переменных .
С ней поступают аналогичным образом и так далее.
Заметим, что
Второй случай заменой переменных
сводится к первому.
Пример 1:Квадратичную форму привести к каноническому виду посредством невырожденного линейного преобразования.
Решение. Соберём все слагаемые, содержащие неизвестное , и дополним их до полного квадрата
.
(Так как .)
или
(3)
или
(4)
и
от неизвестных
формапримет вид.
Далее полагаем
или
и
от неизвестных
формапримет уже канонический вид
Разрешим
равенства (3) относительно
:
или
Последовательное
выполнение линейных преобразований
и
,
где
,
имеет матрицей
Линейное
преобразование неизвестных
приводит
квадратичную форму
к каноническому виду (4). Переменные
связаны с новыми переменными
соотношениями
С LU - разложением мы познакомились в практикуме 2_1
Вспомним утверждения из практикума 2_1
Утверждения (см.Л.5, стр. 176)
Данный скрипт призван понять роль LU в методе Лагранжа, с ним нужно работать в блокноте EDITOR с помощью кнопки F9.
А в прилагаемых ниже заданиях лучше создать свои М-функции, помогающие вычислению и осознанию задач линейной алгебры (в рамках данной работы)
Ax=X."*A*X % получаем квадратичную форму
Ax=simple(Ax) % упрощаем ее
4*x1^2 - 4*x1*x2 + 4*x1*x3 + x2^2 - 3*x2*x3 + x3^2
% найдем LU разложение без перестановки строк матрицы A
% При преобразовании матрицы к ступенчатому виду
%без перестановок строк, мы получим матрицу M1 и U3
% U получается из A U3=M1*A,
% вот такой матрицей элементарных преобразований
0.5000 1.0000 0
0.5000 0 1.0000
%мы получим U3=M1*A, где
4.0000 -2.0000 2.0000
% из M1 легко получить L1, поменяв знаки
% в первом столбце во всех строках кроме первой.
0.5000 1.0000 0
0.5000 0 1.0000
% L1 такое, что
A_=L1*U % вот это и есть нужное нам LU разложение
% Элементы, стоящие на главной диагонали U -
% это коэффициенты при квадратах y i ^2
% в преобразованной квадратичной форме
% в нашем случае, есть один только коэффициент
% значит, в новых координатах будет только 4y 1 2 в квадрате,
% при остальных 0y 2 2 и 0y 3 2 коэффициенты равны нулю
% столбцы матрицы L1 - это разложение Y по X
% по первому столбцу видим y1=x1-0.5x2+0.5x3
% по второму видим y2=x2; по третьему y3=x3.
% если транспонировать L1,
% то есть T=L1."
% T - матрица перехода от {X} к {Y}: Y=TX
0.5000 1.0000 0
1.0000 -0.5000 0.5000
% A2 – матрица преобразованной квадратичной формы
% Заметим U=A2*L1." и A=L1* A2*L1."
4.0000 -2.0000 2.0000
1.0000 -0.5000 0.5000 |
% Итак, мы получили разложение A_=L1* A2*L1." или A_=T."* A2*T
% показывающее замену переменных
% y1=x1-0.5x2+0.5x3
% и представление квадратичной формы в новых координатах
A_=T."*A2*T % T=L1." матрица перехода от {X} к {Y}: Y=TX
isequal(A,A_) % должно совпасть с исходной A
4.0000 -2.0000 2.0000
2.0000 1.0000 -1.5000
2.0000 -1.5000 1.0000
Q1=inv(T) % находим матрицу перехода от {Y} к {X}
% Найдем преобразование,
% приводящее квадратичную форму Ax=X."*A*X
% к новому виду Ay=(Q1Y)."*A*Q1Y=Y." (Q1."*A*Q1)*Y=Y." (U)*Y
Ay =4*y1^2 - y2*y3
x1 - x2/2 + x3/2
% матрица второго преобразования,
% которая составляется значительно проще.
4*z1^2 - z2^2 + z3^2
% R=Q1*Q2, X=R*Z
R=Q1*Q2 % невырожденное линейное преобразование
% приводящее матрицу оператора к каноническому виду.
det(R) % определитель не равен нулю - преобразование невырожденное
4*z1^2 - z2^2 + z3^2 ok
4*z1^2 - z2^2 + z3^2
Сформулируем алгоритм приведения квад ратичной формы к каноническому виду ортогональным преобразованием:
определяет на плоскости кривую. Группа членов называется квадратичной формой, – линейной формой. Если в квадратичной форме содержатся только квадраты переменных, то такой ее вид называется каноническим, а векторы ортонормированного базиса, в котором квадратичная форма имеет канонический вид, называются главными осями квадратичной формы.
Матрица называется матрицей квадратичной формы. Здесь a 1 2 =a 2 1 . Чтобы матрицу B привести к диагональному виду, необходимо за базис взять собственные векторы этой матрицы, тогда , где λ 1 и λ 2 – собственные числа матрицы B.
В базисе из собственных векторов матрицы B квадратичная форма будет иметь канонический вид: λ 1 x 2 1 +λ 2 y 2 1 .
Эта операция соответствует повороту осей координат. Затем производится сдвиг начала координат, избавляясь тем самым от линейной формы.
Канонический вид кривой второго порядка: λ 1 x 2 2 +λ 2 y 2 2 =a , причем:
а) если λ 1 >0; λ 2 >0 – эллипс, в частности, при λ 1 =λ 2 это окружность;
б) если λ 1 >0, λ 2 <0 (λ 1 <0, λ 2 >0) имеем гиперболу;
в) если λ 1 =0 либо λ 2 =0, то кривая является параболой и после поворота осей координат имеет вид λ 1 x 2 1 =ax 1 +by 1 +c (здесь λ 2 =0). Дополняя до полного квадрата, будем иметь: λ 1 x 2 2 =b 1 y 2 .
Пример
. Дано уравнение кривой 3x 2 +10xy+3y 2 -2x-14y-13=0 в системе координат (0,i,j), где i
=(1,0) и j
=(0,1).
1. Определить тип кривой.
2. Привести уравнение к каноническому виду и построить кривую в исходной системе координат.
3. Найти соответствующие преобразования координат.
Решение
. Приводим квадратичную форму B=3x 2 +10xy+3y 2 к главным осям, то есть к каноническому виду. Матрица этой квадратичной формы . Находим собственные числа и собственные векторы этой матрицы:
Характеристическое уравнение:
; λ 1 =-2, λ 2 =8. Вид квадратичной формы: .
Исходное уравнение определяет гиперболу.
Заметим, что вид квадратичной формы неоднозначен. Можно записать 8x 1 2 -2y 1 2 , однако тип кривой остался тот же – гипербола.
Находим главные оси квадратичной формы, то есть собственные векторы матрицы B. .
Собственный вектор, отвечающий числу λ=-2 при x 1 =1: x
1 =(1,-1).
В качестве единичного собственного вектора принимаем вектор , где – длина вектора x
1 .
Координаты второго собственного вектора, соответствующего второму собственному числу λ=8, находим из системы
.
1 ,j
1).
По формулам (5) пункта 4.3.3. переходим к новому базису:
или
; . (*)
Вносим выражения x и y в исходное уравнение и, после преобразований, получаем: .
Выделяем полные квадраты : .
Проводим параллельный перенос осей координат в новое начало: , .
Если внести эти соотношения в (*) и разрешить эти равенства относительно x 2 и y 2 , то получим: , . В системе координат (0*, i 1 , j 1) данное уравнение имеет вид: .
Для построения кривой строим в старой системе координат новую: ось x 2 =0 задается в старой системе координат уравнением x-y-3=0, а ось y 2 =0 уравнением x+y-1=0. Начало новой системы координат 0 * (2,-1) является точкой пересечения этих прямых.
Для упрощения восприятия разобьем процесс построения графика на 2 этапа:
1. Переход к системе координат с осями x 2 =0, y 2 =0, заданными в старой системе координат уравнениями x-y-3=0 и x+y-1=0 соответственно.
2. Построение в полученной системе координат графика функции.
Окончательный вариант графика выглядит следующим образом (см. Решение :Скачать решение
Задание
. Установить, что каждое из следующих уравнений определяет эллипс, и найти координаты его центра С, полуоси, эксцентриситет, уравнения директрис. Изобразить эллипс на чертеже, указав оси симметрии, фокусы и директрисы.
Решение
.
При рассмотрении евклидового пространства мы вводили определение квадратичной формы. С помощью некоторой матрицы
строится многочлен второго порядка вида
который называется квадратичной формой, порождаемой квадратной матрицей А.
Квадратичные формы тесно связаны с поверхностями второго порядка в n - мерном евклидовом пространстве. Общее уравнение таких поверхностей в нашем трехмерном евклидовом пространстве в декартовой системе координат имеет вид:
Верхняя строка - это не что иное, как квадратичная форма, если положить x 1 =x, x 2 =y, x 3 =z:
- симметричная матрица (a ij = a ji)
положим для общности, что многочлен
есть линейная форма. Тогда общее уравнение поверхности есть сумма квадратичной формы, линейной формы и некоторой постоянной.
Основной задачей теории квадратичных форм является приведение квадратичной формы к максимально простому виду с помощью невырожденного линейного преобразования переменных или, другими словами, замены базиса.
Вспомним, что при изучении поверхностей второго порядка мы приходили к выводу о том, что путем поворота осей координат можно избавиться от слагаемых, содержащих произведение xy, xz, yz или x i x j (ij). Далее, путем параллельного переноса осей координат можно избавиться от линейных слагаемых и в конечном итоге свести общее уравнение поверхности к виду:
В случае квадратичной формы приведение ее к виду
называется приведением квадратичной формы к каноническому виду.
Поворот осей координат есть не что иное, как замена одного базиса другим, или, другими словами, линейное преобразование.
Запишем квадратичную форму в матричном виде. Для этого представим ее следующим образом:
L(x,y,z) = x(a 11 x+a 12 y+a 13 z)+
Y(a 12 x+a 22 y+a 23 z)+
Z(a 13 x+a 23 y+a 33 z)
Введем матрицу - столбец
Тогда
- гдеX
T
=(x,y,z)
Матричная форма записи квадратичной формы. Эта формула, очевидно, справедлива и в общем случае:
Канонический вид квадратичной формы означает, очевидно, что матрица А имеет диагональный вид:
Рассмотрим некоторое линейное преобразование X = SY, где S - квадратная матрица порядка n, а матрицы - столбцы Х и У есть:
Матрица S называется матрицей линейного преобразования. Отметим попутно, что всякой матрице n-ного порядка при заданном базисе соответствует некоторый линейный оператор.
Линейное преобразование X = SY заменяет переменные x 1 , x 2 , x 3 новыми переменными y 1 , y 2 , y 3 . Тогда:
где B = S T A S
Задача приведения к каноническому виду сводится к отысканию такой матрицы перехода S, чтобы матрица В приобрела диагональный вид:
Итак, квадратичная форма с матрицей А после линейного преобразования переменных переходит в квадратичную форму от новых переменных с матрицей В .
Обратимся к линейным операторам. Каждой матрице А при заданном базисе соответствует некоторый линейный оператор А . Этот оператор имеет, очевидно, некоторую систему собственных чисел и собственных векторов. Причем, отметим, что в евклидовом пространстве система собственных векторов будет ортогональна. Мы доказывали на предыдущей лекции, что в базисе собственных векторов матрица линейного оператора имеет диагональный вид. Формула (*), как мы помним, это формула преобразования матрицы линейного оператора при смене базиса. Положим, что собственные вектора линейного оператора А с матрицей А - это вектора у 1 , y 2 , ..., y n .
А это означает, что если собственные вектора у 1 , y 2 , ..., y n взять за базис, то матрица линейного оператора в этом базисе будет диагональной
или В = S -1 А S, где S – матрица перехода от первоначального базиса {e } к базису {y }. Причем в ортонормированном базисе матрица S будет ортогональной.
Т. о. для приведения квадратичной формы к каноническому виду необходимо найти собственные числа и собственные векторы линейного оператора А, имеющего в первоначальном базисе матрицу А, которая порождает квадратичную форму, перейти к базису собственных векторов и в новой системе координат построить квадратичную форму.
Обратимся к конкретным примерам. Рассмотрим линии второго порядка.
или
С помощью поворота осей координат и последующего параллельного переноса осей это уравнение можно привести к виду (переменные и коэффициенты переобозначены х 1 = х, х 2 = у):
1)
если линия центральная, 1
0, 2
0
2)
если линия нецентральная, т. е.
один из i
= 0.
Напомним виды линий второго порядка. Центральные линии:
Нецентральные линии:
5) х 2 = а 2 две параллельные линии;
6) х 2 = 0 две сливающиеся прямые;
7) у 2 = 2рх парабола.
Для нас представляют интерес случаи 1), 2), 7).
Рассмотрим конкретный пример.
Привести к каноническому виду уравнение линии и построить ее:
5х 2 + 4ху + 8у 2 - 32х - 56у + 80 = 0.
Матрица квадратичной
формы есть
.
Характеристическое уравнение:
Его корни:
Найдем собственные векторы:
При
1
= 4:
u 1
= -2u 2 ;
u 1
= 2c, u 2
= -c или
g 1
= c 1 (2i
– j).
При
2
= 9:
2u 1
= u 2 ;
u 1
= c, u 2
= 2c или
g 2
= c 2 (i
+2j).
Нормируем эти векторы:
Составим матрицу линейного преобразования или матрицу перехода к базису g 1 , g 2:
- ортогональная матрица!
Формулы преобразования координат имеют вид:
или
Подставим в наше уравнение линии и получим:
Сделаем параллельный перенос осей координат. Для этого выделим полные квадраты по х 1 и у 1:
Обозначим
.
Тогда уравнение приобретет вид: 4х 2 2
+ 9у 2 2
= 36 или
Это эллипс с полуосями 3 и 2. Определим угол поворота осей координат и их сдвиг для того, чтобы построить эллипс в старой системе.
Построим:
Проверка: при х = 0: 8у 2 - 56у + 80 = 0 у 2 – 7у + 10 = 0. Отсюда у 1,2 = 5; 2
При у =0: 5х 2 – 32х + 80 = 0 Здесь нет корней, т. е. нет точек пересечения с осью х !