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

книги / Матричное исчисление с приложениями в теории динамических систем

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
36.1 Mб
Скачать
(6.1.2)

лексное число, модуль которого равен 1 / 1х | . Два вектора х и у на­ зываются ортогональными (x-Ly), если (х, у) = 0 .

В случае численного пространства, элементами которого служат столбцовые матрицы, скалярное произведение определяется так: (а, Ь) = Ь*а. По аналогии с этим, наряду с общепринятым обозна­ чением скалярного произведения векторов пространства R (х, у) мы будем пользоваться и другим условным обозначением ухх, по­ лагая (х, у) = ухх. Здесь ух выступает в роли вектора, «сопряжен­ ного» вектору у.

Свойства (1)—(5) при этом представляются в следующем виде:

1)ухх — (х, у) = (у, х) = хху,

2)ух(ах) = (ах, у) = а(х, у) = аухх,

2')

(ау)хх = (х, ау) = а(х, у) = аухх,

3)

zx(x + у) =

(х + у, z) =

(х, z) +

(у, z) = zxx + zxy,

3')

(у + z)xx =

(х, у + z) =

(х, у) +

(х, z) = ухх 4- zxx,

4)

X х х =

(х, х) £5 0,

 

 

5)

ххх =

(х, х) > 0 (х?ь 0).

 

Длина вектора х равна | х | = Vxxx, а условие ортогональности век­ торов х и у принимает вид

ухх = у = 0.

Обозначая звездочкой переход к эрмитово-сопряженной величи­ не, отметим еще такое свойство:

Пусть пространство R — n-мерное унитарное пространство и & — (е, е2 ... ett) — какой-нибудь базис в нем. Обозначим через х

столбцовую матрицу, составленную из координат хр х2, хп век­

тора х, а через у — столбцовую матрицу, составленную из коорди­ нат уу, у2, ..., уп вектора у в этом базисе. Тогда

X = Гх,

у = £у.

(6.1.1)

В силу свойств (1)—(4)

 

 

/ и

\ Ж / и

\

или

ухх = ух2’хх,

где у* — у2 ... уп) — строчная матрица, эрмитово-сопряженная столбцовой матрице у, а

вх ci

ех Г хх = с 2

• ♦ а х

Учитывая еще первое соотношение (6.1.1), из (6.1.2) находим

 

 

 

ухх = у*Н'х.

(6.1.3)

Здесь

 

 

 

 

 

ef ei

ехе2

ехе,

 

Н' =

=

е2*2

 

 

 

 

 

 

 

 

 

 

e?ei

*;«2

 

 

(е„

(е2, Cj)

(е„. e i )

(е„ е2)

(®2’ ®2)

(е„>« г )

( е Р е„)

(е2’ е„)

(е„.е„)

так что скалярное произведение (6.1.3) можно записать и в форме

У** = X

(A « = (e/. e <) = (e J.e,)).

(6.1.4)

i\ £=1

 

 

В частности, имеем

 

 

ххх =

X V *!**'

(6.1.5)

 

 

i, к — 1

 

На основании свойства (1) эрмитовой метрики

 

^ik ~~eAei — е?Ск ~ h/cr

(6.1.6)

Форма

 

 

п

 

 

Е

hikxi*k’

 

i, к — 1

где hik = hki (iy к = 1, 2 ,..., п), называется эрмитовой формой;

выражение, стоящее в правой части равенства (6.1.4), называется

билинейной эрмитовой формой. Таким образом, квадрат длины вектора унитарного пространства представляется в виде эрмитовой

формы его координат. Матрица Я эрмитовой (билинейной эрмито­ вой) формы в силу (6.1.6) является эрмитово-сопряженной (эрми­ товой), т.е. Н* — Я.

Если базис % в R составлен из ортонормированных векторов

ер е2, е „ ,

т.е. таких,

что

 

е \ — х _

J0 при

i Ф к,

Кег

ск>

при

i= k (/, к = 1 ,2 ,...,« ) ,

то в этом случае Я — диагональная (единичная) матрица и формы (6.1.4), (6.1.5) принимают вид

ухх = (х,у) - 2

xtYlt

ххх = (х, х) = 2 | х(|2.

i

=1

i=1

В следующем параграфе будет показано, что ортонормированный базис существует в каждом унитарном пространстве.

Если в качестве числового поля Ж принято поле вещественных чисел, то метрика, удовлетворяющая условиям (1)—(5), называет­ ся евклидовой. Векторное пространство R над полем вещественных чисел с положительно-определенной евклидовой метрикой называ­ ется евклидовым пространством.

Вевклидовом пространстве скалярное произведение векторов х

иу с координатами соответственно xt n y t (i = 1, 2..., п) представ­

ляется равенством

(х, у) = ухх = у'&х£ х = J

sikxkyk.

i.

k= l

Здесь у' — строчная матрица, полученная из у транспонированием, sik = su (i, k = 1, 2 ,...» n) — вещественные числа. В частности,

П

( х , х ) = X X X = x ' & ^ l f x = 2 S i k X i X k . i, А=1

n

Выражение ^ sutxixk называется квадратичной формой относи-

*,А=1

п

тельно х,, х2, ...» хп. Квадратичной форме J sikxixk отвечает били-

нейная форма sikxtyk. Квадрат длины вектора евклидова про-

<.*«1 странства представляется в виде квадратичной формы его коорди­

нат.

З а м е ч а н и е 6.1.1. Из положительной определенности метрик, введенных в унитарном и евклидовом пространствах, вытекает по­ ложительная определенность соответствующих эрмитовых и квад­ ратичных форм, т.е.

(знаки равенств имеют место только при х = 0).

§ 6.2. Ортонормированные базисы в унитарном и евклидовом пространствах

Здесь мы покажем, что как в унитарном, так и в евклидовом пространстве имеется ортогональный базис, т, е. такая система век­

торов е1? е2, е п (л — размерность пространства), что

 

при

i Ф к,

(6.2.1)

при

i = к.

 

Заметим, что векторы, обладающие свойствами (6.2.1), линейно независимы. В самом деле, равенство

а1е1 + а2е2

+

+ алеп = 0

возможно тогда, когда cij = а2

=

= ап = 0, так как по умноже­

нии этого равенства скалярно, например, на ек, получаем в силу (6.2.1) ак(ек, е*) = 0 , откуда ак = 0 .

Доказательством существования ортогонального базиса в рас­ сматриваемых пространствах может служить приводимый ниже процесс ортогонализации *, который позволяет из данных л линей­ но независимых векторов построить л попарно ортогональных век­ торов.

Пусть 2? = (gx, g2, ...,*gre) — какой-нибудь базис в простран­

стве (унитарном или евклидовом). Построим систему ненулевых векторов е1} е2, ..., еп, удовлетворяющих условию (6.2.1). Поло­

жим

Si*

1• • М

 

82 У12^11

• • •»

(6.2.2)

en = 8« + 7iae1+ ... + 7n-i«eu-i-

Условие (6.2.1) приводит к следующим выражениям для чисел у.

>• • • I л - 1 ; j = 2, 3, ..., л).

(6.2.3)

Этот процесс в иных обозначениях приведен и в § 8.3

Соотношения (6.2.2), (6.2.3) позволяют последовательно стро­ ить взаимно ортогональные векторы ер е2, ...,еп, что и подтверж­

дает существование в R ортогонального базиса. Равенства (6.2.2) можно записать в виде матричного соотношения Е = & + <§Т, где

0

712

Vis

Уы

 

0

0

7гз

Угп

 

— (е2е2 ... еп), Г —

 

 

 

 

0

0

0

Уп- 1

п

0

 

 

0

0

0

 

Отсюда 1? = &(Еп — Г )-1, где Еп — единичная матрица порядка п. От ортогонального базиса можно перейти к ортонормирован-

ному базису путем умножения & справа на диагональную матрицу

diag

’ V(e„, e j

Таким образом, в унитарном (евклидовом) пространстве существу­ ет и ортогональный, и ортонормированный базис.

§6.3. Линейные операторы в унитарном пространстве

6.3.1.Сопряженный оператор. Если для любых векторов х и у из унитарного пространства между линейными операторами А В

В в R выполняется соотношение (Ах, у) — (х, By), то оператор

называется сопряженным по отношению к оператору А. Оператор, сопряженный оператору А, обозначается через А*.

Каждому линейному оператору А отвечает единственный сопря­

женный А*. Покажем это. Выберем в R ортогональный

базис

& = (в! е2 ... еп). Если В — линейный оператор в R, то для произ­

вольного вектора у из R справедливо соотношение

 

П

 

By = 2 (By, е*)е*.

(6.3.1)

jt=i

 

В самом деле, имея в виду, что В — матрица оператора В в ба­ зисе Bv В2, ..., Вп — строки матрицы В, а у — столбцовая мат­

рица координат ур у2, ..., уп вектора у в том же базисе, получаем

П

П

 

By = В&у = &Ву = 2

екВку = 2 В*у ек

(6.3.2)

k=l к=1

(так как В-у — скаляр). Но, с другой стороны, так как & — ортонормированный базис, то

Вку = 2 в 1Уеп е* = (5у, е*).

Учитывая последнее соотношение, из (6.3.2) непосредственно при­ ходим к (6.3.1).

Теперь легко установить, что оператор В, определенный соотно­

шением

 

By = 2 (у, Аек)ек,

(6.3.3)

*=1

 

является сопряженным по отношению к оператору А. Действитель­ но, если А* — оператор, сопряженный оператору А, то для произ­ вольных х и у

(х, By) = (х, 2 (у, Ае*)е*) =* (х, £ (А*у, е*)еА) =

(х, А*у),

*=1

к = 1

 

откуда В = А*. Равенство (6.3.3)

определяет оператор,

сопряжен­

ный оператору А, единственным образом.

Найдем матрицу сопряженного оператора А* в ортонормированном базисе I?. Пусть А — матрица оператора А в базисе £?. Через В обозначим матрицу сопряженного оператора А* в том же базисе. Имеем А& = &A, A*S’ = &В. Отсюда, так как & — ортонормиро-

ванный

базис

и

потому

= Еп,

находим

А = £‘ХА£’,

В = #’ХА*^’. Далее,

 

 

 

 

 

 

В =

(ер

 

 

 

fef А*е,

еГА'еп'

» » • (A Cj

А е2

А*еп) =

* •

• •

 

ех

 

 

 

efA*e,

е?А*еп

 

п)

 

 

 

 

 

П

ft j

 

 

 

 

 

 

 

 

Так как А и А* — взаимно сопряженные операторы, то

 

ех-А*е/. =

(А*е; , ef) =

(еу., Ае,.) = (Ае,, еу) = ex/Aef.

Поэтому

 

 

 

 

еГАе„" *

 

 

ef Ае,

e fA e;

ef Ае,

 

 

 

 

 

 

 

*

= (g*xAё ) ' = А*

efA e,

 

 

 

 

 

е^Аея

eJAe,

ехАе„

 

 

 

 

 

 

 

п

п

 

 

Таким образом, в ортонормированном базисе сопряженным опера­ торам А и А* отвечают сопряженные матрицы А и А*.

6 .3 .2 . Собственные векторы и собственные значения эрмито­ ва оператора. Линейный оператор Н называется эрмитовым, если он совпадает со своим сопряженным: Н = Н*. Для изучения свойств собственных значений и собственных векторов эрмитова оператора

Ннам понадобится следующая

Ле м м а 6.3.1. Пусть А — линейный оператор в векторном пространстве R над полем комплексных чисел Ж, а I — инвари­ антное подпространство пространства R. Тогда оператор. А имеет в подпространстве I хотя бы один собственный вектор.

До к а з а т е л ь с т в о . Пусть размерность подпространства I рав­ на к и & = (g, g2 ... g*) — какой-нибудь базис этого подпростран­

ства. Произвольный вектор х из I представляется в виде

 

х = *i8 i + *282 +

••• +

 

 

 

(6.3.4)

Так как I — инвариантное подпространство, то Agi e l ,

и поэтому

Agj =

^i/Si 4 * ^2/82 4" •••

4*

 

(*

1»

•••»

^)*

Учитывая это, получаем

 

 

 

 

 

 

 

к

к

к

к

 

к

I к

 

\

ах = а 2*,.g. ~ 2^Ag, =

2

2

c j $ j

= 2

2

x t c )i

Sy

i= 1

i=l

f = 1

J = 1

y =ll/ =l

 

(6.3.5)

Условие того, что x является собственным вектором оператора А, отвечающим собственному значению X, т. е. равенство Ах = Хх, в си­

лу (6.3.4) и (6.3.5), можно записать в виде

 

 

к

 

 

 

 

 

2

2

 

cH*i 8; = >-2 x fir

 

j =1 i=l

J

j =1

 

Так как векторы %. ( j 1, 2,..., Л) линейно независимы,

то по­

следнее равенство возможно,, если только

 

к

 

 

 

 

 

 

S CjiXi =

^xj

( / =

1, 2, ..., и),

 

I=1

 

 

 

 

 

 

или

 

 

 

 

 

 

(^ 1 1

 

4 “ ^ 1 2 х 2

’ **

^1 кх к

(6.3.6)

^21Х 1

( ^ 2 2

 

“ “ . ^ ,) ^ с2

**•

С2кХ к ~

C*l*l + C*2*2+‘“ + ( С**

Для доказательства леммы достаточно показать, что существует число X Е X и числа х2> не все равные нулю и удовлет­ воряющие системе (6.3.6). Условием существования ненулевого ре-

шения однородной системы (6.3.6) является равенство нулю его определителя:

С11

^

С12

 

 

(6.3.7)

Ckl

 

Ск2

Но (6.3.7) представляет собой уравнение степени к относительно X с коэффициентами из поля комплексных чисел и потому имеет по крайней мере один (вообще говоря комплексный) корень Х0.

Значит, существует такое число Х0> что при X= Х0 система (6.3.6)

имеет ненулевое решение аг°, х®»

Число Х0 является собст­

венным значением, а вектор

 

*о ^?8 i ^282 “Ь *

"Ь **8 *

— собственным вектором оператора А, так как Ах0 = Хх0. Лемма

доказана. ■ Л е м м а 6.3.2. Собственные значения эрмитова оператора ве­

щественны.

Д о к а з а т е л ь с т в о . Пусть х — собственный вектор, а X — со­ ответствующее собственное значение эрмитова оператора Н, так что Нх = Хх (х ^ 0). Так как Н* = Н, то

(Нх, х) = (х, Н*х) = (х, Нх).

Отсюда Х(х, х) = Х(х, х), и, поскольку (х, х) ф 0, то X = X, что и требовалось доказать. ■

Л е м м а 6.3.3. Пусть Н — эрмитов оператор в п-мерном уни­ тарном пространстве R, а е — его собственный вектор. Тогда совокупность R( векторов х, ортогональных к е, есть ( п — 1)-

мерное подпространство, инвариантное относительно операто­ ра Н.

Д о к а з а т е л ь с т в о . Легко проверить,что Rt есть подпростран­

ство пространства R. Пусть х — произвольный вектор из R. Йредставим его в виде х = Xj + (х — х,), где

(х, е) е. 1 (е, е)

Вектор х, принадлежит одномерному подпространству I, порож­ денному вектором е. Вектор же х — х1 принадлежит подпрост­ ранству Rx, так как

(х — X,, е) = (х, е) - (х,, е) = (х, е) — (х, е) = 0 .

где kj — вещественные числа. Матрица этого оператора не изме­ нится, если предположить, что векторы базиса е,, е2, е и пронор­

мированы. Матрица сопряженного оператора Н* в ортонормирован­ ием базисе получается из матрицы оператора Н транспонированием и заменою каждого элемента комплексно-сопряженным. Проделав это, убеждаемся, что Л = Л*, т. е. операторам Н и Н* отвечает одна и та же матрица. Значит, Н = Н*. Таким образом имеет место

Т е о р е м а 6.3.2. Для того чтобы оператор Н в унитарном пространстве R был эрмитовым, необходимо и достаточно, что­ бы в R существовал ортогональный базис, в котором матрица оператора диагональная и вещественная.

Наконец, отметим еще одно свойство собственных векторов эр­ митова оператора: собственные векторы эрмитова оператора, отве­ чающие различным собственным значениям, взаимно ортогональ­

ны. В самом деле, пусть Het =

XLelv Не2= Х2е2 (X, ^ Х2). Тогда

(Hej,e2) = Х,(е,, е2),

(Cj, Не2) = Х2(в|, е2).

Отсюда, так как (Нер е2) = (ер Не2), то

(Xj —Х2)(ер е2) = 0 и, значит, (е^ е2) = 0.

6.3.3. Унитарный оператор. Линейный оператор U называется унитарным, если UU* = U*U = Е (Е — единичный оператор).

Пустьх и у — векторы унитарного пространства R. Тогда

(Ux, Uy) = (х, ITUy) = (х, у).

Значит, унитарный оператор сохраняет скалярное произведение векторов. И обратно, оператор, сохраняющий скалярное произведе­ ние векторов, унитарен. В самом деле, пусть (Ux, Uy) = (х, у) для

любых векторов х и у. Тогда (U*Ux, у) = (Ех, у), и, следователь­ но, ((U*U — Е)х, у) = 0. Но это равенство выполняется для любых х и у только тогда, когда U*U = Е, т. е. когда U унитарен. При х = у имеем (Ux, Ux) = (х, х). Значит, унитарный оператор сохра­ няет длину векторов.

Выберем в R ортонормированный базис

& = (ej е2 ... еп), и

пусть U — матрица унитарного оператора U в этом базисе. Тогда

U& = 9U.

(6.3.9)

Сопряженному оператору U* соответствует в ортонормированном

базисе сопряженная матрица U*:

 

b*g = £ U \

(6.3.10)

В силу унитарности U

 

U*Uif = UU =

 

Соседние файлы в папке книги