Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Funkts_analiz.pdf
Скачиваний:
90
Добавлен:
09.05.2015
Размер:
760.8 Кб
Скачать

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′=T1 A T X.

Сравнивая последнее равенство с (4.12) и используя определение равенства матриц, сможем написать выражение для матрицы рассматриваемого оператора в новом базисе

B =T1 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 )

= yx,e

i )

=0 ,

i=1,2,...,k

)

.

(5.3)

(

(

 

(

 

 

48

Подставив выражение для вектора x в (5.3), получим k равенств

(yx,ei )=(yx1e1 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 равенств

zb, a

=0;

zb, a

2 )

=0

;

zb, 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 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отличен от нуля и, следовательно, по теореме Крамера получаем выражения

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