книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем
.pdf330 |
К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы |
[ГЛ . 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» |
^ 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).
