Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

itmo114

.pdf
Скачиваний:
20
Добавлен:
21.03.2016
Размер:
5.44 Mб
Скачать

Итак, линейному оператору A в данном базисе соответствует квадратная матрица

a11

a12

...

a1n

 

 

 

 

 

 

 

(7)

A = a 21

a 22

...

a 21

 

,

 

.

...

.

 

 

.

 

 

an 1

an 2

...

ann

 

 

которая называется матрицей линейного оператора A , i -й столбец которой состоит из координат вектора A ei (i =1,2,...,n ) относительно данного базиса.

Отметим, что матрица A оператора A зависит от выбора базиса e1,e2 ,...,en

пространства Λn.

Введём теперь в рассмотрение одностолбцовые матрицы

X = (x1 x 2

... xn )T иY = (y1 y 2

... y n )T

 

 

,

соответствующие векторам x Εn и y Εn . Тогда соотношения (6) в матричном виде можно записать так:

Y = A X

(8)

Итак, мы показали, что всякому линейному оператору A

в евклидовом

пространстве Εn соответствует матрица A ; можно доказать и обратное утверждение: всякую квадратную матрицу A можно рассматривать как матрицу некоторого линейного оператора A в данном базисе e1 ,e2 ,...,en .

Представляют интерес невырожденные линейные операторы, т.е. такие

операторы, матрицы которых имеют обратную A1 , т.е. также являются невырожденными. В этом случае каждому вектору y (образу), определённому

соотношением (1), отвечает единственный вектор x (прообраз) и при этом имеет место матричное равенство:

X= A1 Y .

2.Примеры линейных операторов

e2

r2

M 2 (x 2 ,y 2 )

 

 

 

 

 

 

M 1

(x1

,y 1 )

 

α

 

r1

 

 

 

 

 

 

ϕ

O e1

Рис.5.3.1

120

1. Возьмём в пространстве Ε2 ортонормированный базис e1 ,e2 и рассмот-

рим в этом базисе вектор r1 =x1e1 +y1e2 (или точку M 1(x1,y 1 ) ) (рис.5.4.1). Повернём вектор r1 вокруг начала координат на угол α против часовой

стрелки. Он займёт положение r2 =x 2e1 +y 2e2 , а точка M 1 перейдёт в точку M 2 (x 2 ,y 2 ) , т.е. r2 = A r1 , где A - оператор поворота против часовой стрел-

ки на угол α относительно точки 0. Очевидны равенства

x 2 =r cos(ϕ +α) =x1 cosα y 1 sinα ,

y 2 =r sin(ϕ +α) =x1 sinα +y 1 cosα .

Или в матричном виде

x

2

 

cosα

sinα

x

1

 

 

 

=

cosα

 

 

 

y 2

 

sinα

 

y1

 

Здесь матрица

cosα

sinα

A =

cosα

 

sinα

 

является матрицей поворота.

Тождественным называется преобразование, определяемое соотношением E x = x, x Εn. В частности E ei = ei (i =1,2,...,n ).

Матрица тождественного линейного оператора в любом базисе имеет

вид

1

0

0

...

0

 

 

 

 

 

0

1

0

...

0

 

 

 

 

E =

 

0

0

1

...

0

 

 

 

 

.

.

.

...

.

 

 

 

0

0

0

...

1

 

 

 

 

3.Действия над операторами

1.Сложение линейных операторов. Пусть x Εn, A и B - два линей-

ных оператора в этом пространстве.

Определение 1. Суммой линейных операторов A и B в Εn называ-

ется оператор C, определяемый равенством Cx = A x + Bx , где x - любой вектор из Εn .

Очевидно, что сумма линейных операторов является линейным оператором, причём его матрица C = A + B , где A и B - матрицы линейных операторов A и B .

121

2. Умножение линейного оператора на число. Пусть x Εn , линейный оператор A определён в Εn , α - некоторое число.

Определение 2. Произведением линейного оператора A на число α

называется оператор αA , определяемый равенством

(αA ) x =α(A x), x Εn.

Очевидно, что αA является линейным оператором, а матрица этого линейного оператора получается из матрицы A умножением её на число α , т.е. она равна α A .

3. Умножение линейных операторов. Пусть x Εn, y Εn, z Εn и кроме

того в Εn определены линейные операторы A и B таким образом, что y = Bx, z = A y .

Определение 3. Произведением A B линейных операторов A и B называется оператор C, определяемый соотношением Cx = A (Bx) . Та-

ким образом, перемножение линейных операторов состоит в последова-

тельном их применении по отношению к вектору x .

 

Рассмотрим матрицы – столбцы:

 

 

 

x1

 

y1

 

z1

 

 

 

 

 

 

 

X = x2

;

Y = y2

;

Z = z2

 

 

 

 

 

 

 

...

 

...

 

...

 

xn

yn

zn

и обозначим через A , B и C - соответственно матрицы линейных операторов A , B и C. Тогда очевидно, что Z =A (B X ) = (A B ) X =C X , та-

ким образом, C = A B , т.е. матрица произведения линейных операторов также является линейным оператором.

Действительно,

a) (A B )(x + y) = A (B (x + y)) = A (Bx +By) = A (Bx) + A (By) = = (A B ) x + (A B ) y

б) (A B )(αx) = A (B (αx)) = A (αBx) =αA (Bx) =α(A B )x

Свойства умножения линейных операторов вытекают из свойств умножения матриц.

§4. Замена базиса.

1. Матрица преобразования координат.

Возьмём в пространстве Εn два различных базиса e1 ,e2 ,...,en и

E1, E2 ,..., En .

122

Рассуждение проведём для случая n = 3 . Очевидно, что один и тот же вектор x относительно различных базисов имеет различные координаты.

Действительно, ограничиваясь случаем n = 3 , можем написать:

x =x1e1 +x 2e2 +x 3e3

 

(1)

x =x E +x

E +x E

(2)

1

1

2

2

3

3

 

Любой вектор второго базиса можем разложить по первому базису, т.е.

E1 =τ11e1 +τ21e2 +τ31e3

E2 =τ12e1 +τ22e2 +τ32e3

(3)

E3 =τ13e1 +τ23e2 +τ33e3

Подставим (3) в (2):

x =x (τ e +τ e +τ e ) +x (τ e +τ e +τ e ) +x (τ e +τ e +τ e ) =

=(τ11x1+τ12x2+τ13x3)e1 +(τ21x1+τ22x2+τ23x3)e2 +(τ31x1+τ32x2+τ33x3)e3 (4)

Всилу единственности разложения по данному базису мы должны при-21 2 23 222 2 33 31 11 1 32 3 3 13 131 3 2 12 1

равнять коэффициенты при векторах e1, e2 , e3 в правых частях формул (1) и

(4). Тогда получим:

 

 

x

1

=τ x

 

+τ x

+

τ x

 

 

 

 

 

 

11

1

12

2

13

3

 

(5)

 

 

x

2

=τ x

 

+τ x +τ

x

 

 

 

 

21

1

22

2

 

23

3

 

 

 

 

x

3

=τ x

 

+τ x

+τ x

 

 

 

 

 

31

1

32

2

 

33

3

 

 

Введём в рассмотрение матрицы

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

τ11

 

τ12

...

τ1n

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X = x2 ;

 

 

 

 

 

 

 

 

τ22

...

τ2n

X ′= x2

; T = τ21

 

...

 

 

 

 

...

 

 

.

 

 

.

...

. .

 

 

 

 

 

 

 

 

 

 

 

 

τn 2

...

 

xn

 

 

 

 

 

 

 

τn1

 

τnn

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

Тогда соотношения (5) можно записать в матричном виде

 

 

 

 

 

X =T X .

 

 

 

 

 

Матрица T называется матрицей преобразование координат при пе-

реходе от старого базиса

к

 

 

новому,

т.е.

от

базиса

e1,e2 ,...,en к базису

E1,E2 ,...,En . Причём, столбцами матрицы преобразования координат являются координаты вектора нового базиса E1,E2 ,...,En относительно старого базиса e1,e2 ,...,en .

Если преобразования координат состоит в повороте координатных осей, то матрицаT называется матрицей поворота.

Пример. При повороте координатных осей xOy на угол α мы имели

(рис.5.4.1)

123

x=x1 cosα y1 sinα

y=x1 sinα +y1 cosα

Здесь

 

 

 

 

 

 

 

 

 

cosα

 

 

sinα

- матрица поворота.

 

 

 

 

 

 

 

T =

 

 

 

 

cosα

 

 

 

 

 

 

 

 

 

 

 

sinα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

M (x ,y )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.5.4.1

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Ортогональный оператор и замена базиса

 

 

 

 

Определение 1. Оператор A , матрица которого относительно дан-

ного ортонормированного базиса

e ,e

2

,...,e

n

 

евклидова пространства Εn

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

ортогональна, называется ортогональным оператором.

 

 

 

 

 

Предположим, что в пространстве Εn переход от ортонормированного

базиса

e 0

,e

0 ,...,e 0 к другому ортонормированному базису

E 0 ,E

0 ,...,E 0

 

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

n

осуществляется с помощью преобразования координат X =T X .

Так как

базис E 0 ,E 0 ,...,E 0

- ортонормированный, то

 

 

 

 

 

 

 

 

 

 

 

 

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ei 0 ,Ej 0 )

= 1, i

 

= j

(i , j

=1,2,...,n )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,i j

 

 

 

 

 

 

 

 

 

 

 

 

Выражая скалярные произведения

(E

0 ,E 0 )

(i , j =1,2,...,n )

через коорди-

наты этих векторов, получим:

 

 

 

 

 

 

i

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ 2

+τ

2

+ ... +τ

2

=1

 

 

 

τ

 

τ

 

+τ τ

22

+ ... +τ

τ

n 2

= 0

 

 

11

 

21

 

n 1

 

 

 

 

 

 

11

12

21

 

 

 

n 1

 

 

 

 

2

+τ

2

 

2

 

 

 

 

τ11τ13

+τ21τ23

+ ... +τn 1τn 3 =

 

(14)

τ12

22

+ ... +τn 2

=1

 

 

 

0

 

.

 

.

.

 

.

 

 

 

 

 

 

.

 

 

 

.

 

 

.

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

2

+τ

2

 

2

=

 

τ

 

τ

+

τ

 

 

τ

 

+ ... +τ

 

τ

 

=

0

 

 

τ1n

2n

+ ... +τnn

1

 

1,n 1 1n

 

 

 

2,n 1 2n

 

 

 

n ,n 1 n ,n

 

 

 

 

Из равенств (14) следует, что T T T =E или T T

=T 1 , т.е. матрица T ,

осуществляющая переход от одного ортонормированного базиса к другому, ортогональная.

§5. Изменение матрицы линейного оператора при переходе к новому базису.

124

Пусть в пространстве Εn определён линейный оператор A , т.е.

 

y = A x

 

(1)

или

 

 

 

Y = A X ,

 

(2)

где Y = (y1

 

y 2 ... y n )T , X = (x1

x 2 ...

xn )T - матрицы-столбцы, со-

ставленные из координат векторов

x и y

относительно данного базиса

e1 ,e2 ,...,en , A - матрица линейного оператора A .

 

Выберем в том же пространстве Εn другой базис E ,E ,...,E . Относи-

 

 

 

1 2

n

тельно нового базиса матрица линейного оператора A будет иной. Обозна-

чим через T

матрицу преобразования координат, а через X

и Y - одно-

столбцовые матрицы, составленные из координат векторов x и y относи-

тельно нового базиса, т.е.

(3)

X =T X ,

Y =T Y

(4)

Подставим (3) и (4) в (2), тогда получим:

(5)

T Y ′ = A T X

Умножая левую и правую части равенства (5) слева наT 1 , получим:

(6)

Y ′=T 1 A T X .

Если к тому жеT - ортогональная матрица, т.е. осуществляет переход от

одного ортонормированного базиса к другому, то

(7)

Y ′=T T A T X .

Итак, если в Εn перейти к новому базису, то матрица линейного оператора также изменится и в самом общем случае будет равна

T 1 A T .

§6. Сопряженный и самосопряженный оператор

Пусть в вещественном евклидовом пространстве Εn определён линейный оператор A .

Определение 1. Оператор A в вещественном евклидовом пространстве Εn называется сопряженным по отношению к линейному оператору A в том же пространстве, если его матрица в любом ортонормированном базисе этого пространства является транспонированной по отношению к матрице оператора A .

1.Свойства сопряженного оператора

1.E =E , где E - тождественный оператор, т.е. оператор, матрица которого E единичная в Εn ;

2.(A +B ) = A +B

125

3. (A B ) = B A

5. если A 1 существует, то (A 1 ) = (A )1 .

Определение 2. Линейный оператор A , определённый в вещественном евклидовом пространстве Εn , называется самосопряженным, или симметрическим, если он совпадает со своим сопряженным оператором A , т.е. если A = A .

Очевидно, что матрица самосопряженного оператора совпадает с транспонированной в любом ортонормированном базисе, т.е. является симметричной относительно главной диагонали.

 

2. Свойства самосопряженного оператора

1.

если A = A , B = B , то (A +B ) = A +B = A +B ;

2.

если A - невырожденный самосопряженный оператор, то

(A 1 ) = (A )1 = A 1 .

Доказательство. Действительно, если существует A 1 и кроме того A = A , то в силу свойства 4 сопряженного оператора, получим

(A 1 ) = (A )1 = A 1 ;

3.Если A -самосопряженный оператор в вещественном пространстве Εn , то имеет место равенство:

(A x,z) = (x,A z), x,z Εn

(1)

Действительно, вводя в рассмотрение одностолбцовые матрицы

X и Z ,

и учитывая, что (A B)T

= BT

AT , для скалярного произведения (A x,z) по-

лучим: (A X )T Z =X T

AT

Z .

 

В свою очередь

для

скалярного произведения (x,A z)

имеем

X T AT Z .

 

 

 

Следовательно, равенство (1) верно.

§7.Собственные векторы и собственные значения линейного оператора.

Пусть A - линейный оператор. Пусть x Ε1, где Ε1 некоторое подпространство пространства Εn. Вектор y = A x может принадлежать подпространству Ε1, а может и не принадлежать.

Определение 1. Подпространство Ε1 называется инвариантным по отношению к оператору A , если A x Ε1 , x Ε1.

Определение 2. Ненулевой вектор x называется собственным вектором линейного оператора A , если найдётся такое число λ , что будет выполняться равенство

126

A x = λx .

При этом число λ называют собственным значением (собственным числом) оператора A , соответствующим вектору x . Множество всех собственных значений оператора A называется его спектром.

Остановимся на отыскании собственных значений и собственных векторов линейного оператора A .

Рассмотрение проведём для случая n = 3 .

Итак, пусть в некотором базисе оператор A имеет матрицу

a11

A = a 21a 31

и пусть одностолбцовая матрица X вектору x . Тогда в силу определения

a12

a13

 

a 22

 

 

a 23

a 32

 

 

a 33

= (x1

x 2

... xn )T соответствует

A X = λX => A X λEX = 0 =>

(1)

(A λE) X = 0

Итак, дело свелось к решению системы линейных однородных уравнений, записанной в матричном виде. Очевидно, что эта система имеет ненулевое решение, если det(A λE) = 0 . Уравнение det(A λE) = 0 называется

характеристическим, или вековым уравнением оператора A ; многочлен det(A λE) называется соответственно характеристическим многочле-

ном оператора A . В координатной форме характеристическое уравнение выглядит так:

 

a11 λ

a12

a13

 

 

 

 

 

a 21

a 22 λ

a 23

 

= 0

 

a 31

a 32

a 33 λ

 

 

Решив его, найдём λ1,λ2 ,λ3 - собственные значения линейного оператора. Можно показать, что собственные значения оператора A не зависят от

выбора базиса, т.е. матрицы A и T 1 A T имеют одинаковый набор собственных значений. Далее, для суммы диагональных элементов матрицы A ,

которую называют следом этой

матрицы trA или следом оператор

A (trA ) , справедлива формула

λ1 + λ2 + λ3 =a11 +a 22 +a 33 . Кроме того,

det A = λ1λ2λ3 .

 

После того как найдены собственные значения линейного оператора A , остаётся подставить их по очереди в уравнение (1) и найти соответствую-

щие собственные векторы x(1) ,x(2) ,x(3) .

127

Пример 1. Найти собственные значения и собственные числа линейного оператора, матрица которого

1

2

A =

2

1

 

 

 

Решение. В силу определения собственного вектора можем написать

 

A X λ X = 0 => (A λ E) X = 0 , где

X = (x1 x 2 )T - матрица

столбец, соответствующая искомому вектору x линейного оператора A ;

 

В матричной форме получим:

 

x

 

 

 

 

1λ

2

 

 

0

(2)

 

 

 

1

=

 

2

1λ

x 2

 

0

 

Система однородная, следовательно, она имеет бесчисленное множество решений, если определитель системы равен нулю, т.е. имеем характеристическое уравнение:

 

1λ

2

λ

 

= 0

 

 

 

2

1

 

 

Решая его, получим такие собственные значения λ1 = −1; λ2 = 3 .

Найдём соответствующие собственные векторы. 1) λ1 = −1 подставим в (2), получим

2

2

x (1)

0

(1)

(2)

(1)

= −x 2

(1)

=t (1) ,

 

 

1

 

=

Ù x1

+x 2

= 0 =>x1

 

2

2

 

(1)

0

 

 

 

 

 

 

x 2

 

 

 

 

 

 

 

где t (1) - некоторый параметр. Таким образом, имеем множество коллинеарных векторов, соответствующих первому собственному числу λ1 = −1:

X (1) = t (1)t (1)

Этот вектор нетрудно пронормировать, тогда мы получим единичный собственный вектор, соответсвующий первому собственному числу λ1 = −1

т.е.

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

X (10) =

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

2) λ2 =3 подставим в (10), получим

128

2 2

 

x

(1)

0

(2)

 

(2)

= 0

(2)

 

(2)

=t (2) , т.е.

 

 

1

 

=

=> x1

x 2

 

=>x1

=x 2

2 2

 

 

(2)

0

 

 

 

 

 

 

 

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

; X

(20) =

 

 

 

 

 

 

 

 

X (2) = t

 

 

 

 

 

 

 

 

 

 

t (2)

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В заключение заметим, что множество всех векторов y =A x , где x Εn,

называется областью значений линейного оператора A в Εn, а множество всех векторов x Ε1 Εn, таких, что A x = 0 , называется ядром линейного оператора.

1.Свойства собственных чисел и собственных векторов самосопряженного оператора

Рассмотрим самосопряженный оператор A , определённый в вещественном евклидовом пространстве Εn. В силу определения матрица его A - симметрическая.

Теорема 1. Собственные числа самосопряженного оператора A есть вещественные числа. (Без доказательства)

Теорема 2. Собственные векторы, отвечающие двум различным собственным значениям самосопряженного оператора, ортогональны.

Доказательство. Пусть λ1 и λ2 - различные собственные значения самосопряженного оператора A , а x1 и x2 - соответствующие им собствен-

ные значения. Очевидны равенства:

A x1 = λ1x1

 

 

(1)

A x2 = λ2x2

 

 

(2)

Следовательно,

 

 

 

 

 

(A x ,x

 

) = λ (x ,x

 

)

(3)

1

2

1 1

2

 

(x1,A x2 ) = λ2 (x1,x2 )

 

Но (A x1,x2 ) = (x1,A x2 ) т.е. левые части равенств (3)

равны, следова-

тельно, вычитая их почленно, получим:

(λ1 λ2 ) (x1,x2 ) = 0 => (x1,x2 ) = 0 ,

129

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