- •Предисловие
- •Введение
- •Глава 1. Линейные пространства
- •§ 1. Введение
- •§ 2. Определение линейного пространства
- •§ 3. Свойства линейного пространства
- •§ 4. Линейная зависимость
- •§ 5. Базис и координаты
- •§ 6. Размерность
- •§ 7. Подпространства
- •Глава 2. Евклидовы пространства
- •§ 1. Введение
- •§ 2. Определение евклидова пространства
- •§ 3. Длина вектора
- •§ 4. Неравенство Коши-Буняковского
- •§ 5. Неравенство треугольника
- •§ 6. Угол между векторами
- •§ 7. Ортонормированный базис
- •Глава 3. Линейные операторы
- •§ 1. Определение линейного оператора
- •§ 2. Примеры линейных операторов
- •§ 3. Действия над линейными операторами
- •Глава 4. Преобразование координат
- •§ 1. Замена базиса
- •§ 2. Ортогональные преобразования
- •§ 3. Матрица оператора при замене базиса
- •§ 1. Задача о проекции вектора и перпендикуляре к нему
- •§ 2. Несовместные системы линейных уравнений
- •§ 3. Метод наименьших квадратов
- •Глава 6. Собственные векторы и собственные числа
- •§ 1. Определение собственных векторов и собственных чисел
- •§ 3. Собственные векторы симметричных операторов
- •§ 1. Приведение квадратичной формы к каноническому виду
- •§ 2. Приведение двух квадратичных форм к каноническому виду
- •§ 3. Малые колебания механических систем
- •Глава 8. Элементы теории метрических пространств
- •§ 1. Определение метрического пространства
- •§ 2. Сходимость. Полные метрические пространства
- •§ 3. Принцип сжимающих отображений
- •Библиографический список
- •Предметный указатель
46
Таким образом, определитель ортогональной матрицы равен либо +1, либо −1.
§ 3. Матрица оператора при замене базиса
Легко видеть, что один и тот же линейный оператор в разных базисах имеет различные матрицы. Выясним, как изменяется матрица линейного оператора при изменении базиса. Обозначим через A матрицу некоторого линейного оператора в старом базисе, а через B – матрицу того же самого оператора в новом базисе. Если обозначить через X и Y – одностолбцовые матрицы, элементами которых являются координаты векторов прообраза и образа в старом базисе, а через X′ и Y′ – одностолбцовые матрицы, элементами которых являются координаты векторов – прообраза и образа в новом базисе, то сможем написать
Y =A X; |
(4.11) |
|
′ |
′ |
(4.12) |
Y =B X . |
||
Обозначив через T – матрицу поворота координатной |
системы, |
|
будем иметь |
|
|
|
′ |
(4.13) |
X =T X ; |
||
|
′ |
(4.14) |
Y =T Y . |
||
Подставив (4.11) и (4.13) в (4.14), получим |
|
|
′ |
′ |
|
TY |
=A T X , |
|
откуда следует
Y′=T−1 A T X′.
Сравнивая последнее равенство с (4.12) и используя определение равенства матриц, сможем написать выражение для матрицы рассматриваемого оператора в новом базисе
B =T−1 A T. |
(4.15) |
47
Глава 5. Несовместные системы линейных уравнений и метод наименьших квадратов
§ 1. Задача о проекции вектора и перпендикуляре к нему
Введем |
в рассмотрение евклидово |
пространство R и |
произвольный |
||||||
|
|
|
|
|
′ |
|
|
|
|
вектор y этого пространства. Обозначим через R |
некоторое подпространство |
||||||||
R. Требуется представить вектор y в виде суммы |
|
|
|
|
|||||
|
|
|
y =x +h , |
|
|
(5.1) |
|||
где вектор x |
|
|
|
|
′ |
h ортогонален к |
|||
принадлежит подпространству R |
, а вектор |
||||||||
этому подпространству. |
Вектор x |
называется |
проекцией |
вектора y |
на |
||||
подпространство |
′ |
а вектор |
h – |
перпендикуляром к |
вектору |
x |
|||
R , |
|||||||||
|
|
|
|
|
|
|
′ |
|
|
(перпендикуляром к проекции вектора y на подпространство R ). |
|
||||||||
Примем |
для |
определенности, |
что |
подпространство |
′ |
k -мерное |
и |
||
R |
выберем в нем ортонормированный базис e1, e2,..., ek . Тогда некоторый вектор x можно представить в виде
|
x = x1e1 +x2e2 +...+xkek , |
|
(5.2) |
где числа x1, x2,..., xk |
подлежат определению. Согласно |
условию |
вектор |
h =y −x должен быть |
ортогональным к подпространству |
R′, то |
есть он |
должен быть ортогонален к каждому из векторов базиса e1, e2,..., ek .
Необходимым и достаточным условием этой ортогональности должно быть выполнение k равенств
h,e |
i ) |
= y−x,e |
i ) |
=0 , |
i=1,2,...,k |
) |
. |
(5.3) |
( |
( |
|
( |
|
|
48
Подставив выражение для вектора x в (5.3), получим k равенств
(y−x,ei )=(y−x1e1 −x2e2 −...−xkek ,ei )=
(y,ei )−x1(e1,ei )−x2 (e2,ei )−...−xi (ei ,ei )−...−xk (ek ,ei )=0, (i=1,2,...,k).
|
|
|
|
|
|
(5.4) |
Так как векторы e1, e2,..., ek |
ортогональны и нормированы, то |
|||||
равенства (5.4) могут быть записаны в виде |
|
|
||||
(y,e |
i |
)−x =0 |
, |
i=1,2,...,k |
), |
|
|
i |
( |
|
|||
откуда следует k выражений для искомых чисел xi . |
|
|||||
|
|
x = y,e |
i ). |
(5.5) |
||
|
|
i |
( |
|
Отсюда следует существование и единственность разложения (5.2).
§ 2. Несовместные системы линейных уравнений
Рассмотрим несовместную систему линейных уравнений
a11x1+a12 x2 +...+a1mxm =b1;
a x +a x +..+a x =b ;
21 1 22 2 2m m 2 (5.6)
an1x1+an2 x2 +...+anmxm =bn,
относительно неизвестных x1, x2,..., xm . Так как система (5.6) несовместна, то это значит, что не существует такого набора чисел c1, c2,..., cm , которые при подстановке в систему (5.6) вместо неизвестных x1, x2,..., xm обращали бы каждое уравнение системы в тождество.
49
Подставляя различные наборы чисел c1, c2,..., cm вместо неизвестных x1, x2,..., xm в левые части уравнений (5.6), мы будем получать наборы чисел
z1, z2,..., zm .
Требуется найти такой набор чисел c1, c2,..., cm , чтобы среднее квадратное отклонение соответствующих чисел z1, z2,..., zm от данных величин b1, b2,..., bn , т.е. значение выражения
δ2 = |
n |
(z j −bj )2 |
(5.7) |
∑ |
|||
|
j=1 |
|
|
было наименьшим по сравнению с другими возможными значениями. Заметим, что требование отыскать наименьшее значение величины (5.7) означает требование найти такие значения коэффициентов x1, x2,..., xm , для которых
абсолютные величины ошибок z j −bj были в каком-то смысле малыми в
совокупности.
Для решения поставленной задачи введем в рассмотрение m векторов a1, a2,..., am , компонентами которых являются столбцы коэффициентов при
x1, x2,..., xm соответственно, то есть
a1 ={a11, a21,..., an1},a2 ={a12, a22,..., an2},…,am ={a1m, a2m,..., anm},
b = b , b ,..., b |
. |
(5.8) |
|
{1 2 |
n} |
|
|
Обозначим через z линейную комбинацию векторов (5.8), так что
z =α1a1 +α2a2 +...+αmam , |
(5.9) |
где числа α1, α2,..., αm принимают любые значения.
50
Совокупность всех линейных комбинаций (5.9) образуют подпространство L = L(a1, a2,..., am ). С геометрической точки зрения ясно,
что выражение (5.7) принимает наименьшее значение тогда, когда вектор z −b
совпадает с перпендикуляром к проекции вектора b на подпространство L, а
это значит, что вектор z −b должен быть ортогонален каждому из векторов a1, a2,..., am , т.е. должны выполняться m равенств
z−b, a |
=0; |
z−b, a |
2 ) |
=0 |
; |
z−b, a |
m ) |
=0. |
(5.10) |
|
( |
1) |
|
( |
|
|
( |
|
|
Заменив вектор z во всех m уравнениях (5.10) на соответствующие выражения из (5.9) и произведя очевидные операции, получим систему m линейных неоднородных уравнений относительно неизвестных α1, α2,..., αm .
|
1( 1 1) |
+α |
2 ( 2 |
|
1) |
+ |
...+α |
m |
( |
m |
|
1) |
|
( |
|
1) |
|
|
|
|
||||||||||
α |
a |
,a |
|
|
|
a |
,a |
|
|
|
a |
|
,a |
|
= b,a |
, |
|
|
|
|
||||||||||
α |
a |
,a |
2 ) |
+α |
2 |
a |
|
,a |
2 ) |
+...+α |
|
a |
|
,a |
2 ) |
= b,a |
|
, |
(5.11) |
|||||||||||
|
1( 1 |
|
|
|
( 2 |
|
|
|
|
m ( |
m |
|
|
( |
|
2 ) |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m ( |
m m ) ( |
|
|
m ) |
|
||||||||
1( 1 m ) 2 ( 2 m ) |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
+α |
a ,a |
|
|
+...+α |
|
a |
|
,a |
|
= b,a . |
|
||||||||||||||
α a ,a |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
Так как поставленная задача имеет единственное решение, то |
||||||||||||||||||||||||||||||
определитель системы (5.11) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
(a1,a1) |
|
(a2,a1) |
|
|
|
|
|
(am,a1) |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
D = |
(a1,a2 ) |
|
(a2,a2 ) |
|
|
|
|
|
(am,a2 ) |
|
|
(5.12) |
|||||||||||||
|
|
|
|
|
|
|
|
(a1,am ) (a2,am ) |
|
|
|
|
|
(am,am ) |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
отличен от нуля и, следовательно, по теореме Крамера получаем выражения