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

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

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

Г л а в а III

ЛИНЕЙНЫЕ ОПЕРАТОРЫ

Вл-МЕРНОМ ВЕКТОРНОМ ПРОСТРАНСТВЕ

§1. Кольцо линейных операторов

Линейный оператор, который каждому вектору х из п - мерного векторного пространства R относит некоторый век­ тор у из того же пространства R , будем называть

о п е р а т о р о м в R .

Рассмотрим совокупность всех линейных операторов в R . В этой совокупности операции сложения и умножения естественно ввести так. Пусть А и В — операторы, относя­ щие вектору X векторы

 

 

У і = Ах,

у 2 = Вх.

Тогда суммой операторов А и В назовем оператор С =

= А + В такой,

что

 

 

 

Уі

У2 = (A -f- В) X = Сх.

Далее,

пусть А — оператор,

относящий вектору у век­

тор z, а

В — оператор, относящий вектору х вектор у:

 

 

z = Ay,

 

у = Вх.

Тогда

произведением AB операторов назовем оператор

С = AB,

относящий вектору х

вектор г:

z = А В х = Сх.

Эти операции сложения и умножения, как легко проверить,

обладают следующими

свойствами:

 

А + В = В + А, А + ( В + С) = ( А + В ) + С'

 

А + 0 = А,

А {ВС) = {AB) С,

( 1 . 1)

(А + В)С = АС + ВС,

А( В + С ) = АВ + АС.

 

$ 2]

МАТРИ ЦЫ Л И Н Е Й Н О Г О О П Е Р А Т О Р А В Р А ЗН Ы Х БАЗИСАХ 51

А ,

В рассматриваемой совокупности, наряду с оператором

имеется также обратный (противоположный) оператор

(-А ).

В силу свойств (1.1) и последнего замечания совокуп­ ность всех линейных операторов в «-мерном пространстве R образует некоммутативное кольцо. В этом кольце имеется единичный оператор Е, который каждому вектору х £ R относит тот же самый вектор:

E x = X .

Пусть

 

 

 

У = А х

( х ,у £

R).

 

Через х и у обозначим

столбцовые

матрицы,

элемента­

ми которых служат координаты векторов х и у

в базисе

^li ^2> • ■■) ßfl’

 

 

 

Тогда

у = Ах,

 

 

 

 

 

где А — квадратная матрица порядка «, отвечающая в дан­ ном базисе оператору А.

Линейный оператор А, матрица базисных векторов g и матрица А связаны друг с другом равенством (см. (2.6.9))

Л £ = &4.

(1.2)

Выбором базиса устанавливается изоморфное соответ­ ствие между кольцом линейных операторов и кольцом квад­ ратных матриц «-го порядка. В самом деле, сумме и произ­ ведению двух операторов А и В соответствуют, как это сле­ дует из (1.2), сумма и произведение квадратных матриц А и В, а произведению числа а из Зъ на линейный оператор А отвечает произведение того же числа на матрицу А. Наконец, единичному оператору Е отвечает единичная матрица Е =

=(ö//).§

§2. Матрицы линейного оператора в разных базисах

Рассмотрим в R два базиса g = (ег е2 ... еп) и ^ = = (ві в2 ... еп), связанные друг с другом соотношением

& = 87’,

(2.1)

где Т — неособенная квадратная матрица порядка «, и ли­ нейный оператор А, который произвольному вектору х £ R

5 2

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

[ГЛ. Ill

относит некоторый вектор у £ R. Пусть А и Ах — матрицы линейного оператора А в базисах g и gj соответственно. Тогда согласно (1.2)

i 4 g = g i 4 ,

. d g ^ g A -

(2.2)

Умножая второе равенство (2.2) справа на Т ~ \

получим

с учетом (2.1)

 

 

А 8 = $ТА]Г - \

 

Сравнивая последнее соотношение с первым равенством

(2.2), находим

 

 

А =

ТАгТ~1.

(2.3)

Разрешая (2.3) относительно Ах, получим

 

А, = Т~'АТ.

(2.4)

Две матрицы А и В, связанные друг с другом соотно­ шением вида (2.3) (или (2.4)), называются подобными.

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

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

Определители подобных матриц равны друг другу. Это следует из свойства определителя произведения матриц.

§ 3. Обратный оператор

Принимая во внимание, что определитель матрицы ли­ нейного оператора не зависит от выбора базиса в R , можно ввести понятие определителя линейного оператора, подра­ зумевая под этим определитель матрицы линейного опе­ ратора в каком-нибудь базисе. Определитель линейного опе­ ратора, как и определитель матрицы, обозначается симво­ лами IЛ I и det А.

Оператор А называется особенным (неособенным), если

I = 0 (соответственно 1=7^0). Если оператор неособенный, то 1) из А х = 0 следует х = 0;

§ 4]

С О Б С Т В Е Н Н Ы Е В Е К Т О Р Ы И С О Б С Т В Е Н Н Ы Е З Н А Ч Е Н И Я

5 3

2) AR = R, т. е. векторы А х (Ѵ х R) заполняют

все

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

А х = 0, то в некотором базисе g

В

самом деле, если

 

 

Ах — О,

 

откуда, так как |Л| Ф 0,

х = 0.

 

Далее, пусть у — произвольный вектор пространства R,

у — столбцовая матрица, составленная из координат

век­

тора

у в базисе g, а А — матрица линейного оператора А

в этом базисе. Так как А — невырожденная матрица, то для любой столбцовой матрицы у существует столбцовая матрица

X, определяемая

равенством

 

 

Отсюда

х = А ~ ]у.

(3.1)

У = Ах.

 

 

 

Полученному

матричному

соотношению

соответствует

векторное равенство

 

 

 

у = А х

(x ,y £ R ),

 

т. е. рассматриваемый (произвольный) вектор у £ R есть вектор вида А х (х £ R). Значит, действительно, векторы А х (Y x £ R ) заполняют все пространство R.

Матрицу А~] линейного преобразования (3.1) можно рассматривать как матрицу, соответствующую обратному

оператору Д“ 1в данном базисе пространства R. Оператор Л~! также является линейным в К и

АА~' = А ~ХА = Е ,

что немедленно следует из двух равенств:

у = А х и х — А ~ 1'у.

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

Вектор X £ R называется собственным вектором ли­ нейного оператора А , а число Я £ ді — его собственным зна­ чением, если

А х ^ 'К х .

 

(4.1)

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

(ех е2... еп).

Пусть

А — матрица, отвечающая оператору А

в базисе g,

х

54

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

[ГЛ . Il l

столбцовая матрица, элементами которой служат коор­ динаты вектора х в этом базисе. Имеем, учитывая (1.2),

А X = А %х = %Ах, Хх = XQx = %Хх.

Отсюда в силу (4.1)

&Ах = gta,

и, значит,

Ах =

Хх.

 

(4.2)

 

 

Матричное равенство (4.2) в свою очередь эквивалентно

системе алгебраических уравнений

 

 

(an — X) х1+ аи х„ +

• • •

+

а1пхп = О,

а2і*і+

(а22— ^)*2+

••• +

«2пХ„= О,

 

 

 

 

(4.3)

ciniX-L+

ап2Хо. + • • •

+ (апп Х)хп = 0.

Для того чтобы система линейных однородных уравнений (4.3) имела ненулевое решение, необходимо и достаточно, чтобы определитель системы был равен нулю:

IА ХЕ I =

Уравнение (4.4) представляет собой алгебраическое урав­ нение п-й степени относительно X и называется характери­ стическим уравнением. Многочлен | А ХЕ | называется характеристическим многочленом.

Каждое собственное значение X линейного оператора А является корнем характеристического уравнения (4.4). И на­ оборот, каждому корню X уравнения (4.4) соответствует ненулевое решение хъ х2, ..., хп системы (4.3), и, значит,

числу X отвечает собственный вектор х = ^ х сес = &х one-

ратора А. Столбцовая матрица х,

І

составленная из чисел

xlt х2, ..., хп — решения системы

(4.3),— называется соб­

ственным вектором матрицы.

Уравнение (4.4) имеет не более чем п корней, поэтому линейный оператор А в R имеет не более чем п собственных значений.

S 4] С О Б С Т В Е Н Н Ы Е В Е К Т О Р Ы И С О Б С Т В Е Н Н Ы Е З Н А Ч Е Н И Я 5 5

Пусть 4 — матрица, отвечающая тому же оператору А, но при другом базисе в /?. Матрицы Ах и А подобны:

4 = т ~ 1а т .

Отсюда

4 — ХЕ = Т -'А Т — ХТ~хТ = Т~х (А ХЕ) Т,

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

| 4 — Х£| = \А — ХЕ\.

Таким образом, подобные матрицы имеют один и тот же хар актеристический многочлен.

Собственный вектор оператора (матрицы) определяется с точностью до произвольного множителя, не равного нулю. Действительно, пусть х — собственный вектор оператора А ,

отвечающий собственному значению

X, а сф О . Тогда

 

А (сх) = сАх = сХх =

X (сх).

Отсюда

видно, что сх ф 0 тоже является собственным

вектором,

отвечающим собственному

значению X.

Данному собственному значению X могут соответствовать и несколько линейно независимых собственных векторов. Если собственному значению X отвечают собственные век­ торы X , у, ...., и оператора А, то любая линейная комби­ нация этих векторов либо сама является собственным век­

тором, либо равна нулю. Действительно,

 

А (ax + ßy + • ■■ + ба) = аА х +

$Ау + ■• • +

бА и =

= X(ах + ßy> + • • • + би)

(а, ß, . ..

, б £ Ж).

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

Л е м м а 4.1. Собственные векторы линейного опера­ тора А (матрицы А), отвечающие попарно различным собственным значениям, линейно независимы.

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

А х і ^ Х іХі

(і = 1,2, . . . , k\ X ^ X j при іф і). (4.5)

Допустим противное, а именно, что в условиях леммы собственные векторыл^, х 2, ..., x k линейно зависимы, т. е. имеются числа аи а2, ..., ak£ Ж, не все равные нулю и

5 6

Л И Н Е П Н Ы Е

о п е р а т о р ы

[ГЛ. Ill

такие, что

 

 

 

 

 

 

 

 

агх г + а 2х 2+

• • •

+

akx k = 0.

(4.6)

Пусть, например,

ak ф

0.

 

 

 

Равенство (4.6) умножим слева на Л. Получим, учитывая

(4.5),

 

k

 

 

 

 

 

 

 

 

о-іКх і =

 

 

 

 

 

2

0.

(4.7)

 

 

«=1

 

 

 

 

 

Равенство (4.6) умножим на ^

 

и вычтем затем из (4.7).

Будем иметь

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

і Л

- ^ і)->С(=0.

(4.8)

 

t=2

 

 

 

 

 

 

 

Теперь равенство (4.8)

умножим слева на Л. Придем к

равенству

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

( ^

- / 4 ) ^

= 0.

(4.9)

 

1=2

 

 

 

 

 

 

 

Из равенства

(4.9)

вычтем равенство (4.8),

умноженное

на А#2 * Получим

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ai

 

 

^i) (^т

 

^2 ) Xi — 0.

 

(= 3

 

 

 

 

 

 

 

Повторяя этот прием, в конце концов придем к равенству

ak(^A

^l) (^к

^2 )

• • •

 

l) Хк = 0.

Однако последнее соотношение не может выполняться, так как в левой части, согласно условиям леммы и сделан­ ному предположению о существовании не равного нулю коэффициента ak, все сомножители отличны от нуля. Полу­ ченное противоречие доказывает лемму.§

§ 5. Линейные операторы и матрицы простой структуры

Линейный оператор А в «-мерном пространстве R может иметь не более чем п линейно независимых собственных век­ торов. Если характеристическое уравнение имеет п различ­ ных корней, то оператор А имеет точно п линейно незави­ симых собственных векторов. Если характеристическое уравнение имеет k различных корней (/г •< п), то число ли­ нейно независимых собственных векторов может быть и больше, чем /г, и даже равно п.

§ 5] Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы И М А Т Р И Ц Ы П РОСТО Й С Т Р У К Т У Р Ы 57

Линейный оператор А в п-мерном векторном пространст­ ве называется оператором простой структуры, если А имеет п линейно независимых собственных векторов.

Пусть А — оператор простой структуры и g lt g 2, g n— линейно независимые собственные векторы оператора А:

Agk = hgk (Â = 1,2, . .. , л). (5.1)

Примем эти векторы в качестве базисных векторов. Если

X = # х ,

где $ = (gigz-.-gn), а X — столбцовая матрица, элементами которой служат координаты вектора л: в базисе $ , то

у = А х = А & х= {Ag1A g i

...

A g n)x =

где

= (K g lK g 2 ••• hgn)X = & y,

 

\

/

V i

У=

столбцовая матрица, элементами которой служат коор­ динаты вектора А х в базисе

Таким образом воздействие оператора простой структу­ ры А на вектор х сводится к «растяжению» составляющих этого вектора по собственным направлениям, порожденным

векторами gi, g e, ■■■, g n>с коэффициентами ^ХД 2, ..., Ъ,п. Соотношения (5.1) эквивалентны одному матричному ра­

венству

где

К

Значит, оператору простой структуры А в «собственном» базисе g lt g 2, ..., ^„соответствует диагональная матрица.

В произвольном базисе оператору простой структуры А соответствует матрица А, подобная диагональной матрице:

А = КАК~1.

5 8

Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы

[ГЛ. I l l

Матрица, подобная диагональной матрице, называется

матрицей простой структуры.

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

§ б. Расщепление «-мерного пространства

Пусть R x и /?2 — подпространства «-мерного простран­ ства R. Если и R 2не имеют общих векторов, кроме нуля, и любой вектор X из R представляется в виде

х = х х + х 2

(хх £ Я*, х 2 £ R 2),

(6.1)

то говорят, что пространство R расщепляется на два под­ пространства R 1 и Ro или что пространство R разлагается в прямую сумму подпространств R x и R 2.

Это разложение записывают так:

R = Яі + /?2-

(6-2)

Представление вектора л: в форме (6.1) единственно. Действительно, допуская, что возможно еще другое пред­ ставление

х = х х + х 2,

(6.3)

после вычитания (6.3) из (6.1) придем к равенству двух век­

торов: х х х х £ R x и х 2 — лг2 £ R 2, что невозможно,

ибо

у R i и R 2 нет общего ненулевого вектора.

 

рас­

Т е о р е м а

6.1.

Если п-мерное пространство R

щепляется на два подпространства R x и R 2, т. е.

 

 

 

 

 

R = R x+ R 2,

 

 

 

то сумма размерностей R x и R 2 равна п.

 

 

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

Выберем некоторый базис f x,

f 2,

в подпространствеR xи базис g x, g 2,..., g Lв подпро­

странстве R 2.

 

 

g x,

g 2,

 

 

 

Векторы f x, f 2, ....Л ,

линейно независи­

мы. Действительно, пусть

 

 

 

 

 

а і / і +

а г Л +

+

«*Л +

Рі£і +

 

 

 

тогда

 

 

 

 

+ РгІГг + • ‘ ‘ +

ßigT = 6;

 

 

 

 

 

 

 

 

« і / і +

а 2 Л +

••• +

a k f k

~

 

 

 

 

 

 

 

=

(ßiiTi +

РгЙ*2 +

+ ß l S r/)-

5 7]

П Р О Е К Ц И О Н Н Ы Е О П Е Р А Т О Р Ы И М А Т Р И Ц Ы

59

Левая

часть последнего

равенства

есть вектор из

R x,

а правая

часть — из /?2. Поскольку у подпространств

/?г

и R 2 общим является только нулевой вектор, то

 

 

 

“ і Л + “ г / г +

 

• •• ++

а

/

й

— О,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ßigi +

ß2Sr2 +

 

 

ß/g-/ = 0.

 

 

Отсюда, в

силу

линейной

независимости векторов f lt

/ 2,

=

а 2 =

... =

afe=

0, а из

линейной

независи­

мости векторов gi, g 2,

gi следует,

 

что ßx =

ß2 = ... =

=ß/ = 0.

Следовательно, векторы/^ . . . , / ft, g ^ , ..., g-, линейно не­

зависимы. Так

как все они—векторы /г-мерного

простран­

ства R , то

 

 

 

 

6

+ /С/г.

 

 

 

(6-4)

 

 

 

 

 

 

 

 

Рассмотрим

теперь

п

линейнонезависимых

векторов

е1л е2, ..., еп пространства R. Каждый из этих векторов мо­

жет быть представлен как сумма двух векторов

из

R x и

Rz и, значит, как линейная комбинация векторов/!,

 

g i.......gi-

 

 

 

 

 

 

 

 

 

 

 

ех = a u f! -j- а 12/

a +

• • •

+

aikfk +

 

 

 

 

 

 

 

 

 

 

+

ßllffl +

ßl2g"2 +

+

ßl/S'n

= “ 2 1 / 1+

2

2

/

2+

 

+

a2kfk +

ß2 iS'i +

 

 

 

 

 

 

 

 

 

+ ß22fir2 +

+ß2/g"/.

ßn — “nl/l +

“ /12/

2

+

• • •

+

ankfk +

ßnllTl +

 

 

 

 

 

 

 

 

 

+ ßn2g-2+ •••

+

ßn/gV

Применяя

к системе векторов еъ

ег,

еп и Д ,

gi

лемму 2.2.2,

будем

иметь

 

 

 

 

 

 

 

 

 

 

 

я ■<£ + /•

 

 

 

(6.5)

Сравнивая (6.4) и (6.5), получаем

 

 

 

 

 

 

 

 

 

k +

/ = п.

 

 

 

 

§ 7. Проекционные операторы и матрицы

Пусть дано произвольное расщепление линейного про­ странства R на два подпространства S и Т:

R = S + T .

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