Линейная алгебра (методические указания), часть 2

advertisement
ОДЕССКИЙ НАЦИОНАЛЬНЫЙ УНИВЕРСИТЕТ ИМ. И. И. МЕЧНИКОВА
ИНСТИТУТ МАТЕМАТИКИ, ЭКОНОМИКИ И МЕХАНИКИ
ЛИНЕЙНАЯ АЛГЕБРА
( РЕШЕНИЕ ТИПОВЫХ ЗАДАЧ)
ЧАСТЬ 2
Методические указания для студентов 1 курса
ОДЕССА – 2008
Составители: д-р ф-м н., проф. Варбанец П.Д.,
к-т ф-м н., доц. Савастру О.В.
Рецензенты: д-р ф-м н., проф. Евтухов В.М.,
к-т ф-м н., доц. Белозеров Г.С.
Рекомендовано к печати
Ученым советом ИМЭМ Одесского национального
университета им. И. И. Мечникова
протокол № 1 от 5 февраля 2008 г.
.
СОДЕРЖАНИЕ
Обозначения…………………………………………………4
1. Линейные пространства …………………………………...5
1.1. Линейные пространства и подпространства………….5
1.2. Базис пространства, его размерность…………………6
1.3. Координаты вектора в данном базисе…………….…11
1.4. Сумма и пересечение подпространств………………12
2. Евклидовы и унитарные пространства ………….…........17
2.1. Процесс ортогонализации Шмидта………………….17
2.2.Ортогональные дополнения…………………………..19
2.3. Ортогональная проекция и перпендикуляр на подпространство……………………………………………………..20
3. Операторы в линейных пространствах…………….........23
3.1. Образ, ядро линейного оператора……………………28
3.2. Матрица линейного оператора в данных базисах…..29
3.3. Собственные векторы и собственные значения..…...31
3.4. Канонический корневой базис и жорданова нормальная форма…………………………………………………….34
4. Операторы в евклидовых и унитарных пространствах..40
5. Приведение двух квадратичных форм к каноническому
виду…………………………………………………………...45
Список литературы………………………………………….51
Линейные пространства и линейные операторы представляют собой начало абстрактной части математики, с
которой студенту в дальнейшем неоднократно придется
иметь дело.
Эти методические указания по самостоятельной работе студентов предполагают использование следующего задачника:
И.В.Проскуряков. Сборник задач по линейной алгебре. М., Наука, 1974.
ОБОЗНАЧЕНИЯ
В дальнейшем мы будем придерживаться следующих
обозначений (если в тексте нет специальной оговорки):
 V , X , Y - произвольные пространства над некоторым
полем P ;
 P n - пространство n - мерных строк (столбцов) с
элементами из поля P над полем P (арифметическое пространство).
В частности
 n - действительное n - мерное арифметическое
пространство;
 n - комплексное n - мерное арифметическое пространство;
 Wk  k  1, 2,3 - пространства геометрических векторов (прямой, плоскости, пространства);
 En , E - евклидовы пространства (с указанием размерности или без него);
 L, Lk - подпространства данного пространства ( k индекс, не связанный с размерностью);
 a, b, , x, , ai , b j , xk ,  векторы рассматриваемого
пространства; 0 - нулевой вектор;
4
 , ,
, ,
,i ,  j , k ,
,  ik ,
 скаляры из данного
поля, 0 - нуль этого поля;
 A, B, C ,  линейные операторы, в отдельных случаях – матрицы;
 AÁ , BÁ1 , CÁ2 ,  матрицы линейных операторов в базисах соответственно Á , Á1 , Á2 , ;
 dim V , dim L,  размерности пространств V , L,
 rA , rB ,  ранги операторов (матриц) A, B, ;
;
 .,.  скалярное произведение в данном пространстве;
 .,.  векторное произведение в данном пространстве W3 .
1. ЛИНЕЙНЫЕ ПРОСТРАНСТВА.
Основными типами задач этого параграфа являются
следующие:
А) выяснение вопроса, будет ли данное множество с указанными операциями линейным пространством, подпространством;
В) выделение базиса пространства, определение его размерности;
С) вычисление координат вектора в данном базисе;
D) нахождение суммы, пересечения подпространств, их
размерностей и базисов.
1.1. Линейные пространства и подпространства.
Для решения задач первой группы необходимо знание
аксиом линейного пространства (вообще, не следует приниматься за решение задач любого раздела, не ознакомившись предварительно с основными понятиями и теоремами
5
данного раздела). Заметим, что в группе аксиом линейного
пространства содержатся требования неограниченной применимости, однозначности и замкнутости линейных операций, которые не выделены под отдельными номерами. Распространенная ошибка: забывают проверить выполнение
этих условий.
В тех условиях, когда данное множество состоит из векторов некоторого известного пространства, полезной является следующая теорема (критерий подпространства):
Теорема. Подмножество L векторов пространства
V над полем P является подпространством тогда и только
тогда, когда
1. L замкнуто относительно сложения, т.е.
a, b  L a  b  L ,
2. L замкнуто относительно умножения векторов
на любые скаляры из основного поля P :
a  L,   P  a  L .
Некоторые из задач требуют хорошего знания других
разделов курса (элементарной теории матриц, квадратичных форм, систем линейных уравнений). Ниже мы подробнее остановимся на одной из этих задач.
1.2. Базис пространства, его размерность.
Построение базиса пространства, подпространства несколько упрощается, если мы располагаем некоторыми
представлениями о размерности пространства, подпространства. Одним из наводящих соображений здесь может
быть следующее. Подмножество L векторов пространства
V выделяется из V с помощью дополнительных условий,
накладываемых на векторы. При этом, чем больше таких
условий, тем меньшей, вообще говоря, будет размерность
подпространства L . Если dimV  n , а L выделено с помощью r условий специального вида, то есть основания
ожидать, что dim L  n  r .
6
Задача 1.1. (№1297[4]) Доказать, что множество L пмерных векторов, у которых первая и последняя координаты равны между собой, образует линейное подпространство
пространства n .
Решение. Множество L образует линейное подпространство пространства n , так как удовлетворяет критерию подпространства. Действительно, L выделяется из n
с помощью одного условия 1  n , поэтому
1. a   a1 , a2 ,
, an1 , a1  , b  b1, b2 ,
a  b   a1  b1 , a2  b2 ,
2. a   a1 , a2 ,
, bn1, b1   L
, an1  bn1 , a1  b1   L ,
, an1 , a1   L,  
 a   a1 , a2 , , an1, a1   L .
Кроме того, нетрудно показать, что dim L  n 1 . Для этого
рассмотрим
векторы
стандартного
базиса
ek   0, ,0,1,0, ,0 k  1, n . Векторы e1 , en не принадлежат L . Но построение базиса подпространства в ряде случаев удобно выполнить, исходя из стандартного базиса самого пространства, изменяя его векторы так, чтобы они
«попали» в подпространство. Поэтому преобразуем векторы e1 , en так, чтобы у них первая и последняя координаты
были
равны.
Например,
пусть
e1  1,0,...,0,1 ,
en  1,0,...,0,1 . Рассмотрим систему векторов e1, e2 ,..., en1 .
Она образует базис L , так как нетрудно проверить, что она
является линейно независимой и каждый вектор из подпространства линейно выражается через вектора этой системы. А так как количество векторов системы равно n 1 ,
то и dim L  n 1 . Итак, наше предположение оказалось
верным.
Линейные подпространства, размерности которых на 1
меньше размерности самого пространства называются гиперплоскостями.
7
В следующей задаче условий больше.
Задача 1.2. (№1298[4]) Доказать, что множество L пмерных векторов, у которых координаты с четными номерами равны нулю, образует линейное подпространство пространства n .
Решение. Для доказательства того, что L является
подпространством, нужно также воспользоваться критерием подпространства. Так как 2  4   2k  0  2k  n  ,
dim L  n  k , где k 1
наибольшее четное число, не превышающее n ( k  n , ес2
1
ли n - четное, и k   n  1 , если n - нечетное). Базисом L
2
является подсистема стандартного базиса пространства n ,
содержащая векторы только с нечетными номерами.
поэтому следует ожидать, что
Задача 1.3. Проверить, является ли множество L многочленов степени 3 с вещественными коэффициентами
подпространством пространства многочленов степени  3
( 3 t  ).
Решение. Воспользуемся критерием подпространства.
Проверим условие f  t  , g  t   L
?
f  t   g  t  L .
Пусть f  t   t 3  t , g t   t 3  1 L , тогда
f  t   g  t   t  1 L ,
так как степень суммы этих двух многочленов равна двум.
Итак, множество L не является подпространством.
Задача 1.4. (№№1291, 1308[4]) Найти какой-нибудь
базис и размерность линейного подпространства L пространства n , если L составляют все векторы из n , у которых сумма координат 1  2   n  0 .
Решение. Очевидно векторы стандартного базиса
8
ek   0, ,0,1,0, ,0 (1 на k - ой позиции ) множеству L
не принадлежат ни при каком k . Однако, замена на векторах e1 , e2 ,..., en1 последнего нуля числом (-1) дает нам векторы из L . Таким образом мы получаем систему n 1 векторов
e1  en , e2  en ,..., en1  en
из L , которая линейно независима (почему?) и обязана
быть базисом L , ибо из условия задачи явно следует, что из
L  n и, следовательно, dim L  dim n  n .
Попутно решен вопрос (и подтвердилась гипотеза) о разn
мерности L ( dim L  n  1, L выделено из
одним условием).
Задача 1.4. (№1306[4]) Пусть f - неотрицательная
квадратичная форма от n неизвестных ранга r . Доказать,
что все решения уравнения f =0 образуют  n  1  мерное
линейное подпространство пространства n .
Поиск решения. Вспоминаем основные понятия теории квадратичных форм (матрица формы, ранг формы,
определение формы). Очевидно, что более подробные запиn
n
си данного уравнения в виде f  x   xT Ax    ik i k  0,
 x   ,  ,...,    , никак
T
1
2
n
k 1 i 1
не указывают на способ решения
задачи.
В процессе дальнейших размышлений начинаем понимать, что мы должны исходить из неотрицательной определенности формы f . Нормальный вид такой формы
f  12  22  ...  r2  0r21  ...  0n2 ,
(1)
а множество решений уравнения f =0 в этом случае состоит из векторов вида
(2)
 0,..., 0,r 1 ,...,n  ,
9
Где  j  r  1  j  n  - произвольные числа из
. Имеющийся опыт (задача 1.2) подсказывает, что множество векторов такого вида есть ( n  r )-мерное подпространство
пространства n . Но данная нам форма не обязательно
нормальная. И здесь мы вспоминаем, что каждая неотрицательно определенная форма ранга r невырожденным линейным преобразованием приводится к виду (1). Создается
план решения: преобразовать форму f к виду (1) , найти
решения (2) уравнения f =0 для преобразованной формы, а
затем с помощью обратного преобразования построить решения уравнения f =0 для данной формы f .
Решение. По теореме о приведении квадратичной
формы к нормальному виду существует невырожденное
линейное преобразование
 x   , ,...,  , y   , ,...,  , R  r  ,
x  Ry,
T
1
2
T
n
1
2
n
ik
приводящее форму f к виду
f  x   xT Ax  yT  RT AR  y 
 12  22  ...  r2  0r21  ...  0n2  g  y  .
Множество решений уравнения f  x   0 состоит из
векторов x  Ry , где g  y   0 , то есть из векторов
 0 


x  R 0  ,  j  , j  r  1, n .
 r 1 
 
 n 
Обозначим ek   0, ,0,1,0, ,0 (1 на k - ой позиции) и докажем, что множество L решений уравнения
f  x  =0
T
r 1
есть линейная оболочка системы векторов
T
r 2
Re , Re
,..., RenT
10
L  L  RerT1 , RerT 2 ,..., RenT

.
Пусть f  x   0 . Тогда
 0 


x  R 0   R r 1erT1  ...  n enT   r 1RerT1  ...  n RenT  L.
 r 1 
 
 n 
Очевидно и другое:
f  R ekT   ek RT ARek T  0, k  r  1,..., n.
Кроме того, система RerT1 , RerT 2 ,..., RenT линейно независима (проверяется непосредственно). Составляем линейную комбинацию  r 1 RerT1  ...   n RenT  0 . Получаем
 0  0

  
 
0
R       . Мы пришли к матричному уравнению, ко r 1   
   0
 n 
торое имеет единственное решение, так как матрица R является невырожденной.
 0 
0 0


 
1  
0
   R
  .


 r 1 
0 0
 
   
 n 
Отсюда  r 1  ...   n  0 . Тем самым мы показали, что
система RerT1 , RerT 2 ,..., RenT является линейно независимой.
Следовательно, L - линейное пространство (по построению)
и его размерность dim L  n  r.
1.3. Координаты вектора в данном базисе.
Решение вопроса о ранге системы векторов, заданных координатами в некотором базисе, выделение из си11
стемы ее максимальной линейно независимой подсистемы,
выражение остальных векторов в виде линейных комбинаций векторов этой подсистемы сводится к решению этих же
задач для системы строк (столбцов) координатной матрицы,
которые подробно обсуждались в соответствующем параграфе первой части.
1.4.Сумма и пересечение подпространств.
Пусть L1 , L2 - данные подпространства пространства.
Обычно их задают в виде линейных оболочек систем векторов или как множества решений некоторых однородных
систем линейных уравнений, а сами векторы- координатными строками в некотором базисе. Вычисление
dim  L1  L2  не составляет особого труда: это ранг объединения базисов или порождающих систем подпространств
L1 и L2 . dim  L1 L2  находится по формуле
(3)
dim  L1  L2   dim L1  dim L2  dim  L1 L2  .
Несколько сложнее обстоит дело с поиском базиса пересечения L1 L2 . В общем виде этот вопрос рассматривается в задаче №1319 [4]. Здесь же мы укажем, как найти
решения конкретных задач (№№ 1320-1322 [4]). Задачу 1.6
мы решим двумя способами, второй - с помощью схемы
Штифеля (предполагаем, что №1319 вы уже разобрали).
Задача 1.6. Найти базис суммы и пересечения подпространств, натянутых на системы векторов
a1   2,1,0 , a2  1, 2,3 , a3   5, 2,1 и
b1  1,1, 2 , b2   1,3,0 , b3   2,0,3 .
Решение. Обозначим L1  L  a1 , a2 , a3  ,
L2  L  b1 , b2 , b3  . Будем считать, что координаты векторов
12
заданы в единичном базисе
e1  1,0,0 , e2   0,1,0 , e3   0,0,1 .
1 способ. Как известно, базисом суммы служит любая база
системы векторов a1 , a2 , a3 , b1 , b2 , b3 . Его построение сводится к вычислению ранга матрицы, строками которой являются координаты векторов последней системы. Кроме
того, базис суммы можно получить, добавляя к базису первого подпространства некоторые из векторов базиса второго подпространства.
 2 1 0
2 1 0 




Итак, dim L1    1 2 3     0 3 6   2 . Базис
2 c2  c1
 5 2 1  c3 5c2  0 8 16 




L1 составляют a1 , a2 .
 1 1 2
 2 1 0




dim L2    1 3 0     0 4 2   2 .
c2  c1
 2 0 3  c3  2c2  0 6 3 




ставляют b1 , b2 .
Базис L2 со-
2 1 0 


1 2 3 

dim  L1  L2   {a1 , a2 , b1 , b2 }  
 3.
1 1 2 


 1 3 0 
Базис L1  L2 составляют a1 , a2 , b1 . По формуле (3) получаем
dim  L1
L2   2  2  3  1. Базис пересечения будем искать
из условия x  L1 L2  x  L1  x  L2 . Значит, x представим в виде x  1a1   2 a2 и x  1b1   2b2 . Приравниваем правые части 1a1   2 a2  1b1  2b2 . Это равенство эквивалентно системе трех линейных однородных уравнений
с четырьмя неизвестными. Нужно решить эту систему и по-
13
строить ФСР. Тогда z  1a1   2 a2  1b1  2b2 будет образовывать базис пересечения.
21   2  1   2  0,

1  2 2  1  3 2  0,
3  2  0.
1
 2
Решив систему, строим ФСР.
1  2 1  2
2 14

7 1
3
3
Вектор z  7b1  b2   6,10,14 образует базис L1 L2 .
2 способ. 1) Составим таблицу Штифеля для объединенной
системы векторов a1 , a2 , a3 , b1 , b2 , b3 и перебрасываем
наверх сначала векторы ai , пока это возможно (квадратиками выделены разрешающие элементы). Векторы ek , переходящие налево, не пишем и их координаты не вычисляем.
а)
e1
e2
e3
б)
e1
a1
e3
a1
2
1
0
a2
-3
2
3
a2
1
2
3
a3
-1
-2
1
a3
-5
-2
1
b1
-1
1
2
b1
1
1
2
b1
-7
3
0
b2
-1
3
0
b1
2
0
3
b3
2
0
3

в)
e1
a1
a3
a2
0
8
3
14

5
2
b2
1
-7
3
0
b3
5
6
3
b1
Перебросить a2 наверх вместо e1 невозможно. Следовательно, dim L1 =2, а базис L1 составляют a1 , a3 . Исключаем
из таблицы строку a2 и перебрасываем наверх b j вместо
оставшихся ek .
Из
таблицы
г)
г)
b1
a1
a3
b2
-7
38
14
b3
5
-19
-7
получаем:
dim  L1  L2   3 ,
то
есть
L1  L2  3 и базис суммы образуют векторы b1 , a1 , a3 .
2) Продолжаем работу с таблицей г), перебрасывая
наверх b j вместо находящихся наверху ai , пока это воз-
можно. Как и выше, векторы, уходящие налево, опускаем.
д)
b1
a1
b2
b3
119
0
-7
Вектор b3 перебросить наверх вместо a1 невозможно. Приходим к выводу, что dim L2  2 , базис L2 составляют b1 , b2 .
По (3) dim  L1
L2   2  2  3  1.
3) Возвращаемся к таблице г). Вектор b2 , вошедший в
базис L2 , представим через базис суммы L1  L2 в виде:
b2  38a1  14a3  7b1.
Отсюда находим 38a1  14a3  7b1  b2 : z .
15
Вектор z  0 и z  L1
L2 , а так как dim  L1
то z образует базис пересечения L1
L2   1 ,
L2 . Оба представле-
ния вектора z дают один результат z   6,10,14 , что подтверждает правильность вычислений. Задача решена.
Для более полного усвоения понятия суммы, прямой
суммы подпространств полезно решить задачи №№13231329 [4].
Задача 1.7. Для подпространства L1 , натянутого на
векторы a1   2,3,0,1 , a2   4,3, 2,1 , a3  8,9, 2,3 , найти
дополнительное подпространство.
Решение. Для любого подпространства L1 линейного
пространства V всегда найдется дополнительное подпространство L2 , то есть такое подпространство, что
V  L1  L2 . Причем, оно определяется неоднозначно.
Найдем одно из таких подпространств. Для этого мы должны найти базис a1 ,..., ak подпространства L1 и дополнить
его до базиса всего пространства V .
Пусть
a1 ,..., ak , ak 1 ,..., an - базис V . Тогда L2  L  ak 1 ,..., an  .
Найдем базис и размерность L2 .
2 3 0 1
2 3 0 1




dim L1    4 3 2 1     0 3 2 1  2 .
c2  2 c1
 8 9 2 3  c3  4c2  0 3 2 1




4
Базис L2 - a1 , a2 . Так как L1  L2 
- сумма прямая, то
dim L2  2 . Чтобы найти базис L2 дополним базис L1 до базиса всего пространства
4
e4   0,0,0,1 .
16
векторами
e3   0,0,1,0  ,
2

4
dim  L1  L2    
0

0
3
3
0
0
0
2
1
0
1

1
 4 . Итак, L2  L  e3 , e4  .
0

1
2. ЕВКЛИДОВЫ И УНИТАРНЫЕ ПРОСТРАНСТВА.
Основные типы задач этого параграфа:
 проверка выполнения аксиом скалярного произведения и доказательство его различных
свойств (№№1351-1354, 1384);
 ортогонализация данной системы векторов, построение ортогональных и ортонормированных
базисов (№№1355-1363);
 построение ортогональных дополнений данных
подпространств (№№1364-1368);
 нахождение ортогональных проекций и перпендикуляров на подпространство (№№13691372);
 вычисление длин, расстояний, углов (№№13731406).
2.1. Процесс ортогонализации Шмидта.
Обычно метод ортогонализации Шмидта рассматривают и обосновывают в лекциях. Тем не менее, подчеркнем,
что данная система векторов x1 , x2 ,..., xm и ортогональная,
т.е. полученная из данной методом Шмидта y1 , y2 ,..., ym ,
являются эквивалентными системами - их линейные оболочки совпадают. Поэтому ортогонализация системы векторов, порождающей подпространство L , приводит к построению ортогонального базиса L . Обратим внимание на
некоторые частные случаи, встречающиеся в задачах:
17
1. если подлежащая ортогонализации система x1 , x2 ,..., xm
распадается на две взаимно ортогональные подсистемы
x1 , x2 ,..., xk и xk 1 , ..., xm , то для решения задачи достаточно ортогонализировать каждую из этих подсистем
независимо от другой;
2. если выяснилось, что подсистема x1 , x2 ,..., xk уже ортогональна, то ортогонализацию начинаем с вектора xk 1 ,
полагая
yk 1  xk 1  k xk  k 1 xk 1  ...  1 x1
и дальше по стандартной схеме;
3. если в процессе ортогонализации, полученная система
векторов y1 , y2 ,..., ym содержит нулевой вектор, то можно сразу сказать, что исходная система является линейно зависимой.
Задача 2.1. Применить процесс ортогонализации к
следующей системе векторов из 4 :
x1   2,3,0,1 , x2  1, 2, 1,0 , x3   3, 2,1,0 , x4   2,1,0,1 .
Решение. Можно сразу заметить, что система распадается на две взаимно ортогональные подсистемы x1 , x2 и
x3 , x4 . Поэтому ортогонализируем каждую из подсистем
независимо друг от друга.
f1  x1 , f 2  x2   f1 ,
f3  x3 , f 4  x4   f3 .

 x2 , f1    4     x4 , f3    4
,
.
 f1 , f1  7
 f3 , f3  7
4
1 2
4
 2,3, 0,1    , , 1,   ,
7
7
 7 7
4
 2 1 4 
f 4   2,1, 0,1   3, 2,1, 0     ,  ,  ,1 .
7
 7 7 7 
f 2  1, 2, 1, 0  
18
2.2.Ортогональные дополнения.
Задачи этого раздела не вызовут трудностей, если
разобраться в свойствах решений линейной однородной системы как векторов евклидова (унитарного) пространства.
Рассмотрим пространство n и систему линейных однородных уравнений над :
111  ...  1n n  0,

(4)
...
   ...     0.
mn n
 m1 1
Обозначив ai   i1 ,...,  in 
T
i  1, m и
перепишем систему (4) в виде
 a1 , x   0,

...
 a , x  0.
 m 
x  1 ,...,  n  ,
T
(5)
Пусть L  L  a1 , ..., am  . Тогда уравнения (5) означают, что


x  ai i  1, m и, следовательно, x  L , а каждый вектор
из L является решением системы (4). Итак, множество
решений системы (4) и линейная оболочка ее строк коэффициентов являются ортогональными дополнениями друг
для друга в пространстве n . (Какие изменения надо внести в рассуждения в случае пространства n ?)
Задача 2.2. Найти базис ортогонального дополнения

L подпространства L , натянутого на векторы:
a1   2,1,3, 1 , a2   3, 2,0, 2 , a3  3,1,9, 1 .
Найти уравнения, задающие подпространство L .
Решение. Так как L  L  a1 , a2 , a3  , то L состоит из
множества решений системы уравнений
19
21   2  33   4  0,

 2 4  0,
31  2 2
3    9    0.
3
4
 1 2
Находим фундаментальную систему ее решений (ранг
системы 2)
x1   6, 9, 1,0 , x2   0,1,0,1 .
Следовательно, L  L  x1 , x2  , а система уравнений со
строками коэффициентов x1 и x2
61  92  3  0,

2
 4  0,

задает подпространство L , как множество решений этой
системы (убедитесь: системы векторов a1 , a2 , a3 и x1 , x2
взаимно ортогональны, а объединение их базисов есть базис 4 ).
Аналогичные соображения используются при дополнении ортогональной системы до ортогонального базиса.
2.3. Ортогональная проекция и перпендикуляр на
подпространство.
Известно, что
L  L  E ,
и потому каждый вектор x  E единственным способом
представим в виде суммы
x  y  z,

где y  L, z  L . Вектор y называют (ортогональной) проекцией вектора x на подпространство L и обозначают
prL x , а z - перпендикуляром (ортогональной составляю-
щей) из вектора x на подпространство L : z  ortL x . Очевидно, что
20
prL x  ort L x , ort L x  prL x .
(6)
y  1a1   2 a2  ...   k ak
(7)
Если L  L  a1 , a2 ,..., ak  , то
и тогда
x  1a1   2 a2  ...   k ak  z .
Умножаем последнее равенство скалярно на a j , 1  j  k , с
учетом  a j , z   0 , получаем
1  a1 , a1    2  a2 , a1   ...   k  ak , a1    x, a1  ,

1  a1 , a2    2  a2 , a2   ...   k  ak , a2    x, a2  ,
(8)

...

  a , a     a , a   ...    a , a    x, a  .
2
2
k
k
k
k
k
 1 1 k
Эта система в силу существования представления
(7) совместна. Определитель матрицы этой системы есть
определитель Грама G  a1 , a2 ,..., ak  . Если a1 , a2 ,..., ak - линейно независима, то G  a1 , a2 ,..., ak   0 и система (8) имеет единственное решение. В противном случае у системы
(8) решений бесконечно много. Но нам достаточно найти
одно (все другие решения дадут нам те же векторы y и z ).
Если система (8) получилась несовместной, ищите ошибку.
Вычисления сокращаются, если известны базисы L и L и
если, опираясь на соотношения (6), выбрать то подпространство, размерность которого меньше.
Умение находить prL x и ortL x позволит успешно
справиться с большинством задач на вычисление длин, расстояний и углов.
Задача 2.3. Вычислить расстояние от вектора
x   4, 2, 5,1 до плоскости H , заданной системой уравнений
21
2 x1  2 x2  x3  2 x4  9,
.

2 x1  4 x2  2 x3  3x4  12.
Решение. Напомним: плоскостью (линейным многообразием) называется множество векторов вида
H  x0  L  x0  y | y  L ,
где x0 - вектор сдвига, L - данное (направляющее) подпространство. Любая плоскость может быть задана как множество решений некоторой системы линейных уравнений.
Частное решение этой системы дает координаты вектора
сдвига. Общее решение приведенной системы - направляющее подпространство. Учитывая это, находим, что
H  x0  L , где
x0   3,0,3,0  , L  L  a1 , a2  ,
a1  1,0, 1, 4 , a2   0,1, 1,6 .
Расстояние от вектора x до плоскости H определяется как
inf   x, z  или как inf   x, x0  y  .
zH
yL
Так как
 2  x, x0  y   x  x0  y  prL   x  x0   y   ort L  x  x0  
2
2
 prL   x  x0   y   ort L  x  x0 
2
2
(применена теорема Пифагора), то
inf  2  x, x0  y   ort L  x  x0  .
2
yL
Мы получили формулу   x, H   ort L  x  x0  .
Остается вычислить ortL  x  x0  и его длину.
x  x0  1, 2, 8,1 .
prL  x  x0   y  1a1  2 a2
и тогда
x  x0  1a1   2 a2  z .
22
Умножаем последнее равенство скалярно на a j , 1  j  2 , с
учетом  a j , z   0 , получаем
181  25 2  13,

251  38 2  16.
94
37
Решая систему, находим 1  ,  2   .
59
59
Тогда
94
37
prL  x  x0   1, 0, 1, 4    0,1, 1, 6  
59
59
 94 37 57 154 
  , , ,
.
 59 59 59 59 
 35 155 415 95 
ort L  x  x0     ,
,
,  .
59
59 
 59 59
35
  x, H   ortL  x  x0   10
.
59
3. ОПЕРАТОРЫ В ЛИНЕЙНЫХ ПРОСТРАНСТВАХ.
Пусть X и Y - два произвольных линейных пространства. Как известно, оператором, действующим из X в
Y называется отображение пространства X в пространство
Y . Если отображение обозначить символом A , то это записывают так:
A: X Y .
Образ вектора x обозначают Ax или A  x  и называют
значением оператора A на векторе x . По определению
Ax Y .
Оператор A : X  Y называют линейным оператором, если X и Y - пространства над одним и тем же полем
P и при этом
1. A  x1  x2   A  x1   A  x2  x1 , x2  X
23
(аддитивность оператора);
2. A   x    A  x  x  X ,  P
(однородность оператора).
Понятие линейного оператора является одним из
важнейших в математике. Это подтверждается хотя бы тем,
что основные операторы, изучаемые в математическом анализе и алгебре (предельный переход, дифференцирование,
интегрирование, проектирование на подпространство,
умножение на матрицу и др.) являются линейными.
Оператор A : X  X ( X  Y ) называют также преобразованием пространства X .
Основные типы задач по этой теме:
a) проверка линейности заданного оператора;
b) нахождение образа, ядра, ранга и дефекта линейного оператора;
c) построение матрицы линейного оператора в
данных базисах (в данном базисе);
d) нахождение собственных значений и собственных векторов линейного оператора (№№14651484);
e) построение канонического базиса и жордановой нормальной формы линейного оператора
(№№1529-1536).
Основная трудность задач первой группы состоит в
том, что примеры операторов могут быть взяты из различных разделов математики и требуют от студента эрудиции
и определенной математической культуры. Приведем несколько примеров.
Задача 3.1. Проверьте линейность следующих операторов:
24
 1 
    
1. A : 4  2 ; A  2    1 2  .

 
 3   3 4 
 4
2. B : M 2  t   M 2  t  ; ( M 2  t  -пространство многочленов степени n  2 над некоторым полем).
Bf  t   f   t   t 2 .
3. C : W3 
3
 1  1 
; C 1i   2 j   3k    2  1 .
  
3 





; Dx   x, k  , j .
5. F : M 2  t   M 2  t  ;
4. D : W3 
1
f 1  u   f 1  u 
du .
u
1
6. X  L1  L2 . Определим оператор P так: если
t
Ff  t   
x  X и x  x1  x2
 xi  Li  , то
Px  x1 (опе-
ратор проектирования на L1 параллельно L2 ).
7. Q : W3  W3 ; Qx  x a ( a  0 - фиксированный вектор).
Решение.
1. A : 4  2 является отображением. Проверим аддитивность и однородность.
T
T
x   x1 ,..., xn  , y   y1 ,..., yn   4 ,  
 x1  y1 
 x  y2   x1  y1   x2  y2  
A x  y   A 2


x  y   x  y   x4  y4  
 x3  y3   3 3
 4
4
x x  y y
x x
y y
  1 2 1 2    1 2    1 2   A  x   A  y  .
x

x

y

y
x

x
3
4
 3 4
 3 4   y3  y4 
25
  x1 
  x   x   x2 
 x1  x2 
A   x   A  2    1
 x3   x3   x4     x3  x4    A  x  .
x 
 4
Все условия выполнены, значит, A является линейным
оператором.
5. f  t   M 2  t   f  t   0  1t   2t 2 .
f 1  t   f 1  t    0  1 1  t    2 1  t  
2


 0  1 1  t    2 1  t    21  4 2  t ,
t

1
2
f 1  t   f 1  t 
 21  4 2 .
t
t
f 1  u   f 1  u 
t
du    21  4 2  du   21  4 2  u 1 
u
1
  21  42  t 1 .
Ff  t   M 2  t  . Теперь проверяем аддитивность и однородность. Напомним: если   t   f1  t   f 2  t  , то
    f1    f 2   и  1  u   f1 1  u   f 2 1  u  .
Находим
t
F  f1  f 2   
f1 1  u   f 2 1  u    f1 1  u   f 2 1  u  
1
t

1
t
f1 1  u   f1 1  u 
du  
u
1
u
f 2 1  u   f 2 1  u 
du 
u
 Ff1  t   Ff 2  t  .
Точно так же
t
F   f t   
 f 1  u    f 1  u 
u
1
26
du 
du 
f 1  u   f 1  u 
du   F  f  t   .
u
1
Все условия определения линейного оператора выполнены.
F - линейный оператор.
t
 
Линейный оператор нулевой вектор отображает в нулевой ( A0  A  0 x   0 A  x   0 ). Поэтому, если A0  0 , то
A нелинейный. Рекомендуем в подозрительных случаях
прежде, чем начинать проверку аддитивности и однородности, вычислить значение оператора на нулевом векторе. Так
в упражнении 2 B 0  t 2  0  B - нелинейный. В упражне-
 1   0 
нии 3 C 0   1    0   C - нелинейный.
 0  0
   
Для доказательства нелинейности достаточно привести пример двух векторов, для которых нарушена аддитивность, или пример вектора или скаляра, для которых не выполнена однородность (равенство A0  0 может иметь место и для нелинейных операторов). Например: в упражнении 7 настораживает то, что текущий вектор x находится
под знаком . . Поэтому проверку ведем на конкретных
векторах.
Qi  j   i  j  a  2  a,
Q  i   Q  j   i  a  j  a  2  a.
Очевидное неравенство 2  2 доказывает неаддитивность Q и его нелинейность.
В этом же примере можно поступить и так:
Q  i   i  a  a , 

 a  a .
 Q  i    i  a  a , 

Поэтому оператор Q неоднороден, следовательно, и
нелинеен.
27
Проверку линейности операторов из упражнений 4 и 6
предоставляем читателю.
Чтобы глубже понять определение линейного оператора, придумайте примеры:
1. оператора аддитивного, но не однородного;
2. оператора однородного, но не аддитивного.
3.1. Образ, ядро линейного оператора.
Образом линейного оператора A называется множество всех векторов вида Ax, x  X . Если A : X  Y , то
образ A есть подмножество из Y . Его обозначают AX или
Im A .
AЕсли
линейный
оператор,
то
Im A  L  Ae1 , Ae2 ,..., Aem  , где e1 , e2 ,..., em - какой-либо базис
пространства X .
Ядро линейного оператора A : X  Y - это множество
тех x  X , для которых Ax  0 . Ядро линейного оператора
(обозначается Ker A ) – подпространство пространства X .
Полезно уметь находить ядра и образы линейных операторов, их размерности (дефект и ранг).
Задача 3.2. Найти образ, ядро, ранг и дефект оператора A : W3  W3 , Ax    x , j  , k  (оператор двойного векторного умножения).
Решение. Будем считать, что мы уже убедились в
линейности оператора A .
Вычисление образа. Возьмем стандартный базис пространства W3 : i , j , k . Находим
28



  Im A  L  0, 0, j   L  j 

Ak    k , j  , k     i , k   j , 

(подпространство одномерное).
rA  dim Im A  dim L  j   1 .
Ai    i , j  , k    k , k   0,
Aj    j , j  , k    0, k   0,
Вычисление ядра. Пусть x  1 i   2 j   3k , x  Ker A .
Это означает, что Ax    x , j  , k   0 или
Ax   1i   2 j  3k , j  , k   1 0   2 0  3 j .
Отсюда x  Ker A  3  0  x  1i   2 j , где 1 ,  2 
.
Другими словами Ker A  L  i , j  , а дефект nA  2 .
(В нашем примере Im A  Ker A , но это не общее правило).
Можно было воспользоваться формулой для двойного векторного произведения. Но решение вряд ли упростилось бы
от этого.
Как правило, нахождение ядра в конце концов сводится к решению системы линейных однородных уравнений относительно координат произвольного вектора ядра. В
рассмотренном нами примере эта система оказалась очень
простой
0 1  0 2  13  0,
что позволило нам сразу записать общее решение
1,2 ,0 .
3.2. Матрица линейного оператора в данных базисах.
Обязательно нужно научиться строить матрицу линейного оператора в данных базисах. Но кроме этого, еще
раз обратим наше внимание на следующую теорему: каждый линейный оператор из X в Y однозначно определяется своими значениями на каком-либо базисе пространства
29
X . Эта теорема позволяет строить примеры различных
операторов, удовлетворяющих наперед заданным свойствам.
Задача 3.3. Для каждого из нижеперечисленных условий постройте пример линейного оператора A : 4  4 :
1. Ker A  Im A .
2. Ker A  Im A, Ker A  Im A .
3. Ker A  Im A  4 .
4. AL  L , где L  L  e1 , e2  .
5. На L  L  e1 , e2 , e3  A действует как тождественный,
но A  E .
6. Каждое L  e1  , L  e2  , L  e3  , L  e4  переводит в себя,
но A  E .
Решение. 1. Возьмем какой-либо базис в 4 , например, стандартный
e1  1,0,0,0 , e2   0,1,0,0  , e3   0,0,1,0  , e4   0,0,0,1 .
Так как rA  nA  4 , то из условия Ker A  Im A следует rA  nA  2 . Для определенности возьмем
Ker A  Im A  L  e1 , e2  . Определим A на базисе так:
Ae1  0, 

Ae2  0, 

Ae3  e1 , 
Ae4  e2 .
Этими условиями линейный оператор A полностью
определен.
Если x  1 , 2 , 3 , 4   1e1  2e2  3e3  4e4 , то по
нашему определению
A 1 , 2 , 3 , 4   3e1  4e2  3 , 4 ,0,0 .
30
Легко убеждаемся, что Ker A  Im A  L  e1 , e2  .
Действительно, Im A  L  Ae1 , Ae2 , Ae3 , Ae4   L  e1 , e2 
Ker A - это множество тех x  X , для которых
Ax  3 , 4 ,0,0  0  3  0, 4  0 , то есть
x  1 , 2 ,0,0  L  e1 , e2  .
6. Так как необходимо построить такой линейный
оператор A , который каждое L  e1  , L  e2  , L  e3  , L  e4  переводит в себя, но A  E , то будем считать, что система
e1 , e2 , e3 , e4 является линейно независимой, а значит, является базисом
. Определим A на базисе так:
 Ae1  2e1 ,
 Ae  e ,
 2
2

 Ae3  4e3 ,
 Ae4  2e4 .
Можно проверить, что таким образом введенный операторм является линейным и удовлетворяет всем необходимым условиям.
4
3.3. Собственные векторы и собственные значения.
Процедура вычисления собственных значений и собственных векторов (собственных подпространств) линейного оператора A : X  Y вытекает из соответствующего теоретическоо материала. Продемонстрируем ее на конкретном примере.
Задача 3.4. Найдите собственные значения и собственные подпространства оператора A : 3  3 (необходимо самостоятельно проверить линейность)
 1   41  2  3 
A   2    21  4 2  3  .
  
2  43 
 3 
31
Решение. 1) Строим матрицу оператора A в стандартном
базисе Á пространства 3 (предполагаем, что линейность
оператора проверена):

1   4 
Ae1  A  0    2   4e1  2e2  0e3 , 
0 0

   

4 1 1
 0  1 



Ae2  A 1    4   1e1  4e2  1e3 ,  , AÁ   2 4 1  .
 0  1 

 0 1 4
   



 0  1 

Ae3  A  0   1   1e1  1e2  4e3 . 
1   4 

   
2) Составляем характеристическую матрицу AÁ   E , вычисляем ее определитель и находим корни характеристического многочлена.
1
1 
4


AÁ   E   2
4
1 ;
 0
1
4   

4
1
1
6 6 6
det  AÁ   E  
2
0
1
1
 6    2 4  
0
1
 6   
2
1
4
1 
1
4
2
0
4
1 
1
4
1
1
0
0
1  6    2 2  
1 
4
0
1
4
1
3
 6   
4
1
3

4
1
1
2
  6    3    ;
1 4
1  6; 2  3.
  6    3   
32
Оба корня принадлежат полю
и являются собственными значениями оператора; 1  6 - кратности 1;
2  3 - кратности 2.
3) Составляем систему линейных однородных уравнений с
матрицей AÁ  6 E и находим ее фундаментальную систему
решений:
21   2  3  0;

21  2 2  3  0;
  2  0.
3
 2
Система ранга 2. Множество ее решений - одномерное пространство, линейно независимых решений. Легко находим
ее фундаментальную систему решений - x1   3, 4, 2 . Собственное подпространство, относящееся к 1  6
X 6  L  x1      3, 4, 2   
.
4) Составляем систему линейных однородных уравнений с
матрицей AÁ  3E и находим ее фундаментальную систему
решений:
1   2  3  0;

21   2  3  0;
    0.
 2 3
Система ранга 2. Множество ее решений также является
одномерным пространством. Легко находим
ФСР x2   0,1, 1 . Собственное подпространство, относящееся к
2  3
X 3  L  x2      0,1, 1  
.
Задача решена.
Замечание 1. Если оператор A задан своей матрицей, то пункт 1) не нужен.
33
Замечание 2. Если в процессе решения получилась
несовместная система, то допущена ошибка (неверно вычислен характеристический многочлен, найденное  не является собственным значением или допущена ошибка в вычислении коэффициентов системы линейных уравнений
или в процессе решения системы).
Замечание 3. Полезно помнить, что размерность
собственного подпространства, относящегося к собственному значению  у любого линейного оператора не превышает (меньше или равна) кратности этого собственного
значения, как корня характеристического многочлена (геометрическая кратность собственного значения его алгебраической кратности). Например, в задачах №№1465, 1466,
1481 есть собственные значения с алгебраической кратностью 3, а их геометрическая кратность соответственно равна 1,2 и 3.
Замечание 4. В учебных примерах, как правило,
корни характеристического многочлена вычисляются точно. На практике часто приходится довольствоваться их
приближениями. Возникающие при этом проблемы достаточно сложны и здесь не обсуждаются.
3.4. Канонический корневой базис и жорданова нормальная форма.
Задачи построения канонического корневого базиса
(ККБ) и жордановой нормальной формы (ЖНФ) матрицы
линейного оператора наиболее сложные, так как 1) большой объем вычислений и 2) приходится опираться на обширный и серьезный теоретический материал.
Рассмотрим сначала задачу о построении ККБ. Так
как ККБ линейного пространства есть объединение ККБ
корневых подпространств данного оператора, то можно
предположить, что мы уже выделили корневые подпространства и имеем дело с одним из них - L , относящимся к
34
собственному значению  , кратность которого совпадает с
dim L  k .
В лекционном курсе ККБ корневого подпространства
был построен в виде системы башен убывающей этажности,
нижний этаж которых состоял из линейно независимых
собственных векторов, а в каждом столбце нижестоящий
вектор получался из непосредственно вышестоящего в результате применения оператора A   E .
Процедура практического построения ККБ такова:
1.
опираясь на какой-либо базис L , строим систему
башен, вообще говоря, линейно зависимую;
2.
элементарными преобразованиями, сохраняющими башенную структуру, преобразуем систему так, чтобы ее нижний этаж состоял из линейно независимых векторов (этим будет обеспечиваться ЛНЗ всей системы).
Для этого элементарные преобразования будем выполнять сразу над целым столбцом: его перемещение, прибавление к другому столбцу (с меньшей этажностью) и т.п.
Если x1 , x2 ,..., xk - некоторый базис пространства, то
исходная система башен, упорядоченных по высоте, имеет
вид:
x1
xq1
 A   E  x1
2
 A   E  x1
 A   E  xq
2
 A   E  xq
 A  E
 A  E
t 1
x1
xq1 1
1
1
t 1
xq1
xq2
 A   E  xq 1
1
 A  E
t 2
xq1 1
 A   E  xq
2
 A  E 
t 2
xq2
xqt 1 1 ... xqt
Преобразования выполняем слева направо, выбирая их по
векторам нижней строки.
35
Рассмотрим конкретный пример.
Задача 3.5. Постройте ККБ оператора A : 4  4

 1    2  3

 2  

2
.
A   
3
1  23   4 
  

4
 4 

Решение. 1) Строим матрицу оператора A в стандартном базисе e1 , e2 , e3 , e4 пространства 4 . Находим образы базисных векторов:
1   0 
 0  1 




0
0
Ae1  A      , Ae2  A 1   1  ,
0
1
0
0
0 0 
0  0
   
   
 0  1 
0  0




0
0
Ae3  A      Ae4  A  0    0  .
1
2
0
1
0 0
1  1 
   
   
Тогда
 0 1 1 0


0 1 0 0
AÁ  
.
 1 0 2 1 


 0 0 0 1
2) Вычисляем характеристический многочлен оператора и
собственные значения:
1
 

0 1 
det  AÁ   E   
 1
0

0
 0
1
0 

0
0 
4
 1    ;
2
1 

0
1  
  1 кратности 4.
36
Пространство 4 является корневым пространством, относящимся к собственному значению   1 .
3) Берем произвольный базис (например, стандартный) и к
каждому его вектору применяем оператор A  E до получения нуля:
 A E
 A E
 A E
 A E
e1   1, 0, 1, 0  ;
 A  E  e1   0, 0, 0, 0  ;
2
 A  E  e2   1, 0, 1, 0  ;
2
 A  E  e3   0, 0, 0, 0  ;
2
 A  E  e4  1, 0,1, 0  ;
2
e2  1, 0, 0, 0  ;
e3  1, 0,1, 0  ;
e4   0, 0,1, 0  ;
 A  E  e2   0, 0, 0, 0  ;
3
 A  E  e4   0, 0, 0, 0  .
3
Можно сделать некоторые предварительные выводы: максимальная высота корневого вектора равна 3, поэтому в
ЖНФ матрицы оператора будет клетка порядка 3 и, следовательно (так как 4-3=1), одна клетка порядка 1. Хотя этого
задача не требует, но мы можем написать ЖНФ матрицы
оператора:
1 1 0 0 


0 1 1 0

AÁ 
.
0 0 1 0


0
0
0
1


4) Полагаем  : A  E . Выписываем полученные в 3) векторы (кроме нулей, разумеется) в башню:
e2   0, 1, 0, 0 
e4   0, 0, 0,1
 e2  1, 0, 0, 0 
 e4   0, 0, 1, 0 
e1  1, 0, 0, 0 
e3   0, 0, 1, 0 
 e2   1, 0, 1, 0 
 e4  1, 0, 1, 0 
 e1   1, 0, 1, 0 
 e3  1, 0, 1, 0 
2
2
37
и рассматриваем нижнюю строку.
Замечаем:
а) третий столбец векторов пропорционален (равен) соответствующей части первого, а четвертый - соответствующей части второго. Выбрасываем без сожаления третий и
четвертый столбцы.
в) нижний вектор второго столбца пропорционален своему
соседу слева. Прибавляем ко второму столбцу первый,
умноженный на 1, и получившийся нуль во втором столбце
отбрасываем. Таблица примет вид:
 0, 1, 0, 0 
1, 0, 0, 0 
 0, 1, 0, 1 .
 1, 0, 1, 0  1, 0, 1, 0 
Нижняя строка все еще линейно зависима. Ко второму
столбцу прибавляем первый, умноженный на 1 (верхний
вектор первого столбца, у которого нет соседа справа, в
этой операции не участвует), отбрасываем получившийся
нуль и приходим к таблице:
 0, 1, 0, 0 
1, 0, 0, 0 
.
 1, 0, 1, 0  1, 1, 1, 0 
Нижняя строка этой таблицы линейно независима, поэтому
и вся система векторов линейно независима. Она, занумерованная снизу вверх и слева направо, образует ккб:
f1  e1  e3 , f 2  e1 , f3  e2 , f 4  e1  e2  e4 .
Упражнение. Постройте матрицу оператора A в базисе f .
Сравните ее с ЖНФ из 3).
38
ЖНФ матрицы оператора можно отыскать, не выполняя построения ККБ. Число жордановых клеток каждого порядка и максимальный порядок жордановых клеток
для каждого собственного значения  могут быть вычисk
лены, если известны ранги матриц  Ae   E   k  1, 2,..., t  ,
где t - показатель степени такой, что
r1  r2  ...  rt 1  rt  rt 1  rt  2  ...
( t - «момент» стабилизации ранга).
В нашем примере
 1

0
AÁ  E  
 1

0
 AÁ  E 
2
 AÁ  E 
3
1
0
0
0
0 1

0 0

 0 1

0 0
0

0

0

0
0
0
0
0
0

0
, r1  2;
1

0
1
0
1
0
1
0
0
0
0
0
0
0
0

0
, r2  1;
1

0
0

0
, r3  0.
0

0
«Момент» стабилизации ранга t  3 , так как r4  0 невозможно по определению ранга.
Есть клетка третьего порядка (и только одна) и, очевидно, еще одна клетка первого порядка (результат, полученный нами в 3)).
39
Напомним, что в общем случае ( dim X  n , t - «моk
мент» стабилизации рангов матриц  AÁ  E  ) число клеток, относящихся к собственному значению  определяется формулами:
клеток порядка t
: rt 1  rt ,
клеток порядка t 1 : rt 2  2rt 1  rt ,
клеток порядка t  2 : rt 3  2rt 2  rt 1 ,
…………………………………………
клеток порядка 2
: r1  2r2  r3 ,
клеток порядка 1
: n  2r1  r2 .
Следует помнить, что хотя жнф матрицы определена
однозначно с точностью до порядка клеток вдоль главной
диагонали, ККБ существует бесконечно много. Поэтому не
удивительно, если найденный вами ККБ не совпадает с ответом в сборнике задач (но проверить свое решение полезно).
Замечание. Если f1 , f 2 , f3 , f 4 - ККБ оператора A и
T - матрица перехода от базиса e1 , e2 , e3 , e4 к ККБ, то имеет
место равенство:
1 1 0 0 


0 1 1 0
1

T AÁT 
.
0 0 1 0


0 0 0 1 
Таким образом нами «попутно» найдена преобразующая матрица Е, приводящая данную матрицу AÁ к ЖНФ.
4. ОПЕРАТОРЫ В ЕВКЛИДОВЫХ И УНИТАРНЫХ
ПРОСТРАНСТВАХ
Пусть X и Y - два произвольных пространства, оба
евклидовых или оба унитарных. Рассмотрим линейный
40
оператор A : X  Y . Напомним, что оператор A : Y  X
называется сопряженным по отношению к оператору A ,
если x  X и y  Y выполняется равенство
 Ax, y    x, A y  .
Для всякого оператора A сопряженный оператор
A существует и единственный. Причем

n
A y    y, Aek  ek ,
(9)
k 1
где e1 , e2 ,..., en - какой-либо ортонормированный базис
(ОНБ) X .
Равенство (9) можно применять за определение сопряженного оператора.
Матрица A размерности n  m с элементами aij
называется сопряженной по отношению к матрице A размерности
с
элементами
если
aij ,
m n
aij  a ji , i  1, n, j  1, m .
Заметим, что в любых ОНБ унитарных пространств
X и Y сопряженному оператору соответствует соряженная
матрица, справедливо и обратное.
Если мы рассматриваем евклидовы пространства X
и Y , то таким же образом устанавливается соответствие
между сопряженными операторами и транспонированными
матрицами.
Задача 4.1. Найти сопряженный оператор для оператора A : 3  3

 x1   x1  ix2
A  x2    2 x1  1  i  x3  .
x  

x2
 3 

3
В
введено естественное скалярное произведение
3
 x, y    xi y j .
i 1
41
Решение. В заданном унитарном пространстве 3
стандартный базис e1 , e2 , e3 является ОНБ. Для построения
сопряженного оператора воспользуемся равенством (9).
1  1 
Ae1  A  0    2  ,  y, Ae1   y1  1  y2  2  y3  0  y1  2 y2 ,
 0  0
   
0 i 
Ae2  A 1    0  ,  y, Ae2   y1  i  y2  0  y3  1  iy1  y3 ,
 0  1 
   
0  0 
Ae3  A  0   1  i  ,  y, Ae3   1  i  y2 ,
1   0 
  

3
A y    y, Aek  ek 
k 1
  y1  2 y2  e1   iy1  y3  e2  1  i  y2  e3 .
Итак, искомый сопряженный оператор имеет вид
 y1   y1  2 y2 

A  y2    iy1  y3  .
 y   1 i y 
 2
 3  
Задача решена.
Задача 4.2. В пространстве
2
[t ] введено скалярное
1
произведение
 f  t  , g  t     f t  g t dt
1
1 0 0


и задана матрица AT   2 1 0  линейного оператора 
0 0 1


2
в базисе T  {1, t , t } .


Построить матрицу AT сопряженного оператора   в
базисе T  {1, t , t 2} .
42
Решение. Проверим, является ли базис T  {1, t , t 2}
ортонормированным в заданном евклидовом пространстве.
1
1
1
2
2
2
2
,
1,
t

1

tdt

0,
1,
t

t
dt

t
,
t

t 3dt  0 ,
  
  



3
1
1
1
1
1
2
 t , t    t 2 dt  ,
3
1
1
Ортогонализируем систему T .
f1  1 ,
1,1   1dt  2 ,
1
 t 2 , t 2    t 4dt  52 .
1
f2  t ,
f3  t 2    f1 ,   
 t 2 , f1 
 f1 , f1 
2
1
3  ,
2
3
1
f 3  t 2    f1    t 2 .
3
Осталось нормировать полученную систему
f1
1
,
f1 

 f1 , f1  2
f 2 
f3 
f2

3
t,
2
f3

5
1  3t 2  .

8
 f2 , f2 
 f3 , f3 
Мы построили ОНБ F  { f1, f 2, f3} . Теперь мы
должны, используя матрицу перехода от одного базиса к
другому, перейти к матрице оператора в ОНБ. Зная, как
связаны матрицы сопряженных операторов в ОНБ, можно
построить матрицу сопряженного оператора   в базисе F ,
а затем вернуться опять к исходному базису.
Матрица перехода от базиса T к базису F имеет вид
43
P  PT  F










P 1  PT1 F  PF T
1
0
5


8
2
0
0

3
2
0

 2


 0


 0

0 ,

5

3
8 
0
2
3
0
2 

3 

0 .

8 

3 5
0 0
 1



2
5
AF  P 1 AT P  
1 
.
 3
3


0 1
 0

Мы знаем, как связаны матрицы сопряженных операторов в
ОНБ
AF   AF  ,

T
поэтому

2
1
3



AF  0 1

5

0 
3

44

0

0 .


1

Используя матрицу перехода PT  F , возвращаемся к
исходному базису T
3


1 4 0




AT  PAF P 1   0 1 0  .


5
1
0 
2


Задача решена.
5. ПРИВЕДЕНИЕ ДВУХ КВАДРАТИЧНЫХ ФОРМ К
КАНОНИЧЕСКОМУ ВИДУ
Рассмотрим две вещественные квадратичные формы
f  x1 , x2 ,..., xn  и g  x1 , x2 ,..., xn  . Можно ли заданные формы
единым преобразованием привести к каноническому виду?
Эту задачу помогают решить результаты, относящиеся к
линейным операторам. Мы рассмотрим случай, когда одна
из этих квадратичных форм, например f , является положительно определенной. Тогда выполняем сначала преобразование X  QY , которое приводит форму f к нормальному виду (сумме квадратов переменных). При этом форма
g перейдет в новую форму от переменных y1 , y2 ,..., yn . На
следующем шаге выполняется ортогональное преобразование Y  PZ , которое приводит форму g к каноническому
виду. Квадратичная форма f при этом не изменится, так
как ее матрица является единичной, а PT EP  P 1EP  E .
Итак, результирующим преобразованием, которое приведет
обе квадратичные формы к каноническому виду, причем
положительно определенную представит в виде суммы
квадратов, будет X  QPZ .
Задача 5.1. Для заданной пары квадратичных форм
найти невырожденное линейное преобразование, которое
приводит эти формы к каноническому виду.
45
f  x1 , x2 , x3   2 x12  2 x22  x32  2 x1 x2  2 x2 x3 ,
g  x1 , x2 , x3   2 x12 
1 2 4
x2 
2 x32  2 x1 x2  2 2 x2 x3 .
2
3
Решение.
Перепишем формы
f
и g в виде
f  X T AX и


2
1
0 
 2 1 0


1


T

 2  - матg  X BX , где A   1 2 1  , B  1


2
0 1 1




4
0  2 
2


3
рицы соответствующих квадратичных форм.
2 1 0
2 1
 3  0, A  1 2 1  1  0 , то согласТак как 2  0,
1 2
0 1 1
но критерию Сильвестра, форма f является положительно
определенной. Поэтому по ней можно восстановить соответствующую билинейную форму и ввести в 3 скалярное
произведение
 a, b  2a1b1  2a2b2  a3b3  a1b2  a2b1  a2b3  a3b2 .
Оно удовлетворяет всем аксиомам скалярного произведения (положительная определенность формы необходима
для выполнения аксиомы 4 , а именно a  3  a, a   0;
 a, a   0  a  0 ).
3
Рассмотрим
стандартный
базис
в
:
e1  1, 0, 0 , e2   0,1, 0 , e3   0, 0,1 .
Используя введенное скалярное произведение, ортогонализируем его:
46
e1  e1  1, 0, 0  ,
 e  , e  e   e  1 e     1 ,1, 0  ,


2
 2

e , e  
 e  , e  e    e  , e  e   e  2 e   0e    1 ,  2 ,1 .
e e 


3
3 3 
e , e   e , e  
e2  e2 
1
2
1
1
1
2
3
2
3
3
1
1
3
1
1
2
2
2
1
3
2
1
3
Нормируем вектора e1 , e2 , e3 и получаем ОНБ в
, в котором билинейная форма (следовательно, и квадратичная
форма f ) будет иметь единичную матрицу.
1   1

e1 
e1  
, 0, 0  ,
 2

e1
e2 
e3 
1
e2
e2 
2   1
2 
e2   
,
,0,
3
6 3 

2
 1

e3  3e3  
,
, 3 .
3
 3

e3
1
Матрица перехода от старого базиса к новому задает матрицу Q невырожденного преобразования переменных
X  QY квадратичных форм f и g .
1 
 1  1
 2
6
3 


2
2 

Q 0

.
6
3


 0
0
3 




T
T T
2
2
2
f  X AX  Y Q AQY  Y T AY
1  y1  y2  y3 .
Действительно,
47
 1
0
 2

2
 1
T
A1  Q AQ   
6
6

2
 1


3
 3
1

 1
0  
 2
2

  2
3
2  

 0
 0
6
6

 
1   0

0
 0
 
3 


0 
  2 1 0
 

0  1 2 1 Q 
 

  0 1 1 

3


1
6
2
6
0


 1 0 0
2  

  0 1 0 .

3
  0 0 1 
3 


1
3

g  X T BX  Y T QT BQY  Y T B1Y  y12  4 y1 y2 .
Аналогично,
 1


0
0  
 2
2
1
0 

 

2
1
 1
 
T
B1  Q BQ   
0  1
 2 Q 
 

2
6
6

 

2
 1
 0  2  4 2 

3

 

3
3
 3


 1
1
1 
1

0  
 2
6
3 
2

  2
 1 0 0 

2  
2
2  

 0
0

  0 0 2  .
 0

3 
6
3


  0 2 0 


0
3 
2 2  0

 0  6 
 


3 

Далее используем метод приведения квадратичной формы к
главным осям.
48
Характеристический многочлен
  1 0 0 

FA1      0
 2      1   2    2  имеет три
 0
2  

корня 1  1, 2  2, 3  2 , которым соответствуют следующие собственные вектора: 1  1,0,0 ,  2   0, 1,1 ,
3   0,1,1 . Они являются попарно ортогональными, так
как соответствуют разным собственным значениям, и образуют собственный ортогональный базис. Осталось его пронормировать:
1
1
1 1 

1 
1  1, 0, 0  ,  2 
 2   0, 
,
,
1
2
2 2

1 1 
,
.
3
2 2

Теперь составляем ортогональную матрицу, столбцами которой являются векторы 1 ,  2 ,  3 ,
 3 
1

 3   0,

1
0

1

P  0 
2

1
0


2

0 

1 
.
2

1 

2

1
0

1

Тогда матрица T  QP   0 
2

1

0
2

49

0 

1 
2 
1 

2

 



1
0
0
2
6
3

 
2
2  
1
1 
0

 0 

6
3 
2
2

 
1
1 
0
0
3  0
 

 
2
2
 1
1 
1 
1 
1 

1 

1 

6
2
6
2 
 2


1
1
 0 
2 2 
2 2 
6
6




3
3
 0

6
6


и будет искомой матрицей невырожденного линейного преобразования переменных



T  QP  




1
1


1



1
1 
1 
1 
1 
z1 
1 
 z2 
1 
 z3 ,
2
6
2
6
2
1
1
x2  
2  2 z2 
2  2 z3 ,
6
6
3
3
x3 
z2 
z3 ,
6
6
приводящего формы f и g к каноническому виду
f  z12  z22  z32 ,
x1 




g  z12  2 z22  2 z32 .
Задача решена.
50
СПИСОК ЛИТЕРАТУРЫ
1. Завало С.Т., Костарчук В.Н., Хацет Б.И. Алгебра и теория чисел, ч.1. – К.: Вища школа, 1980.
2. Кострикин А. И. Введение в алгебру. - М.: Наука, 1977 .
3. Курош А.Г. Курс высшей алгебры. – М.: Наука, 1975.
4. Проскуряков. Сборник задач по линейной алгебре. М.,
Наука, 1974.
Download