книги / Матричное исчисление с приложениями в теории динамических систем
..pdfлексное число, модуль которого равен 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, а условие ортогональности век торов х и у принимает вид
ухх = Xху = 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 = |
|