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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
97
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

330

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ . X I I I

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

II

ВкУ= ( S ВіУві, ek) = (By, ek).

(=i

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

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

П

 

By = H i(y,A ek)ek,

(3.3)

fc=l

 

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

П

 

П

 

 

 

( X, By) = ( X ,

(у, А ек) ek) = (х, Е

(А*у, ek) ek) =

 

А= I

 

к—1

 

 

откуда В = А*.

 

 

 

=

(лг, А*у),

 

 

 

 

 

Равенство (3.3) определяет оператор, сопряженный опе­

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

оператора А*

 

Найдем матрицу сопряженного

в орто-

нормированном

базисе

g.

А

в базисе g. Через В

Пусть А — матрица

оператора

обозначим матрицу сопряженного оператора А* в том же базисе. Имеем

 

 

 

A g = gA,

А*g = gß.

 

Отсюда, так как g — ортонормированный базис и потому

g x g =

Еп, находим

 

 

Далее,

 

А = g x,4g,

ß = g xyl*g.

 

 

 

 

 

 

 

еГ

 

f е?А*е1 ...

é\A*en

В =

I

••• I

(A*exA * e ,... A*en) = .............................

 

 

 

 

. . .

e*A*enj

Так

как А и

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

 

е?А*е] =

(A*ej, et) = (е/} A et) = (Aeit е/) =

efA e^

$ 3] л и н е П н ы е о п е р а т о р ы

в у н и т а р н о м п р о с т р а н с т в е

331

Поэтому

 

 

 

е*Ае 1 .

. впАе 1

еі А ег ... è\ A еп

 

В =

 

=

 

_ е ? А е п .

еп А е п_

_ е%Аег . . . впАеп__

 

 

 

=

А*.

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

Аи А*.

3.2.Собственные векторы и собственные значения эрми­

това оператора. Линейный оператор Н называется эрмито­ вым, если он совпадает со своим сопряженным: Н = Н*.

Для изучения свойств собственных значений и собствен­ ных векторов эрмитова оператора Н нам понадобится

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

Д о к а з а т е л ь с т в о .

 

Пусть

размерность

подпро­

странства /

равна

k

и

& =

{g1

g 2 --gk) — какой-нибудь

базис этого

подпространства.

 

представляется

в

виде

Произвольный

вектор

х

из /

 

* = X\g\ +

4g2 +

• • •

+ xkg k.

 

 

(3-4)

Так как / — инвариантное

подпространство, то A g& I,

и потому

 

 

 

 

 

 

 

(i= 1, 2,

 

 

 

A g t = с1(§Т + c2ig 2+

• • •

+

ckig k

. . . ,

k).

Учитывая

это,

получаем

 

 

 

 

 

 

 

k

 

к

 

 

 

к

к

2Cjigf =

 

 

 

А X = А 2xtgt =

2xtA gl =

2

 

 

 

 

i=i

 

i=i

 

 

 

;=i

і=і

 

 

 

 

 

 

 

 

 

 

=

è ( è x , - c p ) g p

(3.5)

 

 

 

 

 

 

 

 

;=1 \t'=l

/

 

 

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

А х = Хх,

3 3 2 К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы [ГЛ . X II I

в силу (3.4) и (3.5) можно записать в виде

к / к \ п

 

S

2

СцхЛ «7 = Ь 2

Xjgj.

 

 

 

/=1v=i

/

 

J=1

 

 

 

Так как векторы gj {j — 1,2,..., k) линейно независимы,

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

 

к

срх{

hxj

U ~

 

 

п),

 

2

П 2, . . .

 

і=і

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

( С 1 1 —

^ ) Х 1 +

С12Х 2 “ Ь

• • ■

+

С1кХк

~

 

С21Х1

(с22

х 2 +

' - '

+

С2кХк —

^ 0^

с к1Х 1 +

С к2Х 2 +

-

+ {С/ік — ^) X k

О-

 

Для доказательства леммы достаточно показать, что су­

ществует число Я £ Ж.и числа хѵ лг2, ..., хк, не

все равные

нулю и удовлетворяющие системе (3.6).

 

 

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

£ ц

X

'

*

C l k

 

 

 

 

 

 

 

 

(3.7)

 

Ck l

CkZ

’ ‘

C k k - Ь

Но (3.7) представляет собой уравнение степени k отно­

сительно Я с коэффициентами

из поля комплексных чисел

и потому имеет по крайней мере один (вообще говоря, комп­

лексный) корень

Я0.

 

 

 

Значит,

существует такое

число Я0, что

при Я =

Я0 си­

стема (3.6)

имеет

ненулевое

решение х°и

х°, ..., х\.

Чис­

ло Я0 является собственным значением, а вектор

 

 

=

x°igi + 4 g 2 + • • ■ + xlg k

 

собственным вектором оператора Д.так как А х 0 = Я0л:0. Лемма доказана.

Ле м м а 3.2. Собственные значения эрмитова операто­ ра вещественны.

Д о к а з а т е л ь с т в о . Пусть х — собственный век­ тор, а Я — соответствующее собственное значение эрмитова оператора Н, так что

Н х = Хх (х Ф 0).

$ 3] Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы В У Н И Т А Р Н О М П Р О С Т Р А Н С Т В Е 3 3 3

Так как Н* = Н, то

 

 

 

 

 

Отсюда

(Нх, X) = (х, Н*х) =

(л:, Нх).

 

%(х, х) =

%{X, X),

 

 

 

и, поскольку (х, х) ф О, К =

X, что и требовалось доказать.

Л е м м а

3.3. Пусть Н эрмитов оператор в п-мер-

ном унитарном пространстве

R,

а е — его собственный

вектор. Тогда совокупность R 1

векторов х,

ортогональных

к е, есть (п \)-мерное

подпространство,

инвариантное

относительно

оператора

Н.

 

 

проверить, что R 1есть

Д о к а з а т е л ь с т в о . Легко

подпространство пространства

R.

 

 

Пусть X — произвольный вектор из R. Представим его

в виде

X = X! 4

(X

Xj),

 

где

 

Вектор X} принадлежит одномерному подпространству /, порожденному вектором е. Вектор же х х х принадлежит подпространству R lt так как

(X — х ѵ е) = {X, е) (хѵ е) = (х, е) — (х, е) = 0.

Таким образом, произвольный вектор х

из R представ­

ляется в виде суммы двух векторов: х х £ /

и х

Xj £ R v

Значит, R есть прямая сумма подпространств I

и R v По­

скольку R

я-мерно, / одномерно,

то, значит, размерность

R x равна

п — 1.

 

 

 

Покажем, что R ±инвариантно относительно Н.

Пусть X £ R lt так что (х, е) =

0. Тогда

 

 

(Нх, е) = (х, Н*е) = (х, Не) =

(х,Хе) = X (х, е) = 0.

Это значит, что Н х £ R v т. е. оператор Н переводит век­ торы из R r в векторы того же подпространства R v что до­ казывает инвариантность подпространства R 1 относитель­ но эрмитова оператора Н.

Т е о р е м а 3.1. В п-мерном унитарном пространстве R существует п попарно ортогональных собственных векто­ ров эрмитова оператора Н. Соответствующие им собствен­ ные значения вещественны.

334

К В А Д Р А Т И Ч Н Ы Е

И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ. X III

 

 

Д о к а з а т е л ь с т в о .

По лемме 3.1 в R существует

хотя бы один собственный вектор ех оператора Н. Совокуп­ ность векторов из R, ортогональных вектору еѵ согласно лемме 3.3 образует (п — 1)-мерное подпространство R v инвариантное Н. В этом подпространстве оператор Н имеет хотя бы один собственный вектор е2. Далее, совокупность векторов из R lt ортогональных е2, образует (п — 2)-мерное инвариантное подпространство R z, в котором оператор Н имеет по крайней мере один собственный вектор е3, и т. д. Следуя этим путем, мы получим искомую систему п попарно ортогональных собственных векторов е1г е2....... еп. Соглас­ но лемме 3.2 соответствующие собственные значения вещест­

венны. Теорема

доказана.

 

 

 

ва

Ортогональную систему собственных векторов эрмито­

оператора Н

примем

в

качестве базиса в R: g =

=

(ßi, е2...еп). В данном случае

 

 

Н ес=

( і— 1,2,

. . . , п),

или (в компактной записи)

 

 

 

 

 

//g =

gA,

(3.8

где

 

 

 

 

 

 

 

Л = diag (X.J, Я2, . . . ,

Х„).

Соотношение (3.8) свидетельствует о том, что в выбран­ ном базисе матрица эрмитова оператора Н имеет диагональ­ ную структуру, причем диагональные элементы матрицы вещественны.

Допустим теперь, что матрица некоторого оператора Н в ортогональном базисе имеет вид

где "к-, — вещественные числа. Матрица

этого оператора не

изменится,

если предположить, что

векторы базиса

еъ

ег....... еп

пронормированы. Матрица

сопряженного

опе­

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

§ 3] Л И Н Е Й Н Ы Е

О П Е Р А Т О Р Ы В У Н И Т А Р Н О М П Р О С Т Р А Н С Т В Е

335

отвечает

одна

и та же матрица. Значит, Н = Н*. Таким

образом,

имеет место

 

Т е о р е м а

3.2. Для того чтобы оператор Н е унитар­

ном пространстве R был эрмитовым, необходимо и доста­

точно, чтобы в R существовал ортогональный базис, в

ко­

тором матрица оператора диагональна и вещественна.

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

Тогда

Н е1 = Х1е,

Не2 = Х2е2

(к1ф Х2).

 

 

 

 

 

 

 

в2)=

{ві,

 

(&1*Н в2)~

^ 2 (^1> ß2)*

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

(е1г Не2), то

 

 

 

(^і — Ю (ßv е2) =

О,

 

и, значит,

(еѵ е2) = 0.

 

 

 

 

 

3.3.

Унитарный оператор. Линейный оператор U называ­

ется унитарным, если

 

 

 

 

 

 

 

U U* =

U* U = Е

 

 

— единичный оператор).

унитарного

пространства R.

Пусть

х и у

— векторы

Тогда

( U x , Uy) =

(X, U *U y) =

(л г,з і).

 

Значит, унитарный оператор сохраняет скалярное про­

изведение

векторов.

 

 

 

 

 

И обратно, оператор, сохраняющий скалярное произве­

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

векторов X и у

(Ux, Uy) =

(*,}>).

 

Тогда

 

 

 

(U*Ux,y) =

(Ex, у),

 

 

 

 

и, следовательно,

((U*U— E )x ,y ) = 0.

Но это равенство выполняется для любых х и у только тог­ да, когда U*U = Е, т. е. когда U унитарен.

При X = у имеем

(Ux, Ux) = (х, х).

336

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы ФОРМЫ

ІГЛ. X II I

Значит, унитарный оператор сохраняет длину векторов. Выберем в R ортонормированный базис g = (eh e2... e„), и пусть U — матрица унитарного оператора U в этом ба­

зисе. Тогда

t/'g = W -

(3.9)

Сопряженному оператору U* соответствует в ортонормированном базисе сопряженная матрица ІІ*\

U*g = g<7*.

(3.10)

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

 

{ /* ig = t/t/*g =

g£,

где E — единичная матрица. С другой стороны, принимая во внимание (3.9) и (3.10), имеем

и и *g =

и%и* =

%ии*,

 

и*и% =

и*%и =

g и*и.

 

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

 

 

 

UU* = U*U = E.

(3.11)

Матрица U, обладающая свойством (3.11), называется унитарной матрицей. Таким образом, унитарному опера­ тору в ортонормированном базисе отвечает унитарная мат­ рица.

Пусть U — унитарный оператор в R, а g = (ег, е2...е„) — ортонормированный базис в R. Оператор U переводит си­ стему векторов ех, е2, ..., еп в новую систему g lt g 2.......g n, так что если & = (g1 g 2 ... g„), то

# = g£/.

(3.12)

Принимая во внимание (3.9), из (3.12) находим

= £/*gxg£/.

(3.13)

В силу ортонормированности базиса g g xg =

Е. Учиты­

вая также и соотношение U*U = Е, получаем

 

= Е.

 

Значит,

унитарный оператор переводит ортонормирован-

пый базис

снова в ортонормированный базис. Из (3.13)

видно, что справедливо и обратное утверждение:

оператор,

переводящий ортонормированный базис g (gxg =

Е) в орто­

нормированный базис $ ( $ х$ = Е), унитарен: U*U = Е.

§ 3]

Л И НЕЙНЫ Е ОПЕРАТОРЫ В УНИТАРНОМ ПРОСТРАНСТВЕ

3 3 7

Итак, для того чтобы оператор

U был унитарным,

необ­

ходимо

и достаточно, чтобы он

переводил какой-нибудь

ортонормированный базис в ортонормированный базис. Так как матрица унитарного оператора в ортонормиро-

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

З а м е ч а н и е . Пусть g и $ — заданные ортонорми­ рованные базисы. Тогда матрица перехода от базиса g к базису $ в соответствии с равенством

& = &

определяется

формулой

 

 

 

 

f ( g v ß i ) • • •

( g n , e i ) \

 

 

U

= g x# = .......................................

(3.14)

 

 

W o e„) ...

(gn, en)J

 

3.4.

Преобразование эрмитовой матрицы

к диагональ­

ному виду

с

помощью унитарной

матрицы.

Рассмотрим

эрмитову матрицу Н порядка п. Будем рассматривать Н как матрицу эрмитова оператора Н в ортонормированном ба­

зисе $

= (g1 g 2 ■.. g n) «-мерного унитарного пространства

/?, так

что

 

(3.15)

В пространстве R существует ортонормированный базис g = (ег, е2...еп), в котором матрица оператора Н диагональна и вещественна (см. п.3.2). Обозначим эту матрицу через Л. Тогда

Н g = gA.

(3.16)

Далее, существует унитарная матрица U, которая пре­ образует ортонормированный базис gß ортонормированный базис $■:

$ = W -

(3.17)

Подставляя (3.17) в (3.15), получаем

Щ = g и н и ~ \

Отсюда, сравнивая с (3.16), находим

Л = и н і г \

3 3 8

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ . X I I I

ИЛИ

л = и ни*,

так как в силу унитарности матрицы U U* = U~l. Разрешая полученные соотношения относительно мат­

рицы Н, имеем

Н= LT'AU = U*AU.

§4. Линейные операторы в евклидовом пространстве

4.1.Транспонированный оператор. Симметрический оп ратор. Рассмотрим /г-мерное евклидово пространство R и линейный оператор А в нем.

Оператор А' называется транспонированным по отно­ шению к А, если для любых векторов х и у пз R

{Ах, у) = (х .А 'у).

Аналогично тому, как это было сделано в п. 3.1 для со­ пряженного оператора, устанавливается существование и единственность транспонированного оператора А'. Если elt е2, •••> еп — ортогональный базис, то транспонированный оператор можно определить формулой (см. (3.3))

П

А 'у = l] ( y ,A e k)ek, fc=i

где у — произвольный вектор из R.

Далее, в ортонормированном базисе транспонированным операторам А я А' отвечают транспонированные матрицы А и А'. Заметим, что матрица линейного оператора в евк­ лидовом пространстве вещественна.

Линейный оператор А называется симметрическим, если А' = А. Исследуем свойства собственных векторов и соб­ ственных значений симметрического оператора А.

При исследовании свойств собственных векторов и соб­ ственных значений эрмитова оператора в п. 3.2 была ис­ пользована лемма 3.1. В случае векторного пространства над полем вещественных чисел, вообще говоря, нельзя утвер­ ждать, что любой линейный оператор в соответствующем инва­ риантном подпространстве имеет хотя бы один собственный вектор. Однако в отношении симметрического оператора

§ 4 ] Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы В Е В К Л И Д О В О М П Р О С Т Р А Н С Т В Е 339

в евклидовом пространстве такое утверждение справед­ ливо. Докажем это.

Л е м м а

4.1. Пусть А симметрический оператор

в евклидовом пространстве R,

а I инвариантное подпро­

странство

пространства R.

Тогда

оператор А имеет в

подпространстве I хотя бы один собственный вектор.

Д о к а з а т е л ь с т в о .

Пусть

размерность подпро­

странства / равна k, а $ = (gx g"2 ... g k) — ортонормирован­

ный базис этого подпространства.

 

 

Произвольный вектор х

из /

представляется в виде

X = *1gl + X2g t +

■• •

+ Xkg k.

 

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

A g t I,

и поэтому

 

 

 

 

A gi = clig 1 + c2!g 2 + ■■■ + ckig k

(i = 1 , 2 , . . . ,

k). (4.1)

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

 

 

 

 

A x = 2

( 2

 

gp

 

/=I

w=i

 

/

 

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

спхі +

с\%хг +

■•

'

+ cikxk —

^хѵ

С 2 1 Х 1 + С 2 2 Х 2 +

' • •+ О 2 k X k =

 

4 2)

CklXl +

Ck2X2 +

• • •+

CkkXk =

^xk- ■

 

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

сп — Х

С22 Af •

Clk

 

 

С21

C2k

= 0.

(4.3)

 

 

 

Ck1

Ск2

. . ckk К

 

 

Для доказательства леммы нужно показать, что урав­ нение (4.3) (/г-й степени относительно X) имеет вещественный корень А,0 и что этому корню отвечает вещественное решение

х°, х\, ..., х\ системы (4.2).

Соседние файлы в папке книги из ГПНТБ