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

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

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

ГЛАВА 3

ЛИНЕЙНЫЕ ОПЕРАТОРЫ В /2-МЕРНОМ ВЕКТОРНОМ ПРОСТРАНСТВЕ

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

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

Рассмотрим совокупность всех линейных операторов в R. В этой совокупности операции сложения и умножения естественно ввести так. Пусть А и В — операторы, соотносящие вектору х векторы У1 = Ах, у2 = Вх. Тогда суммой операторов А и В назовем опера­

тор С = А + В такой, что

У1 + у2= (А + В)х = Сх.

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

z = Ау,

у = Вх.

Тогда произведением АВ операторов назовем оператор С = АВ, соотносящий вектору х вектор z:

z = АВх = Сх.

Эти операции сложения и умножения, как легко проверить, об­ ладают следующими свойствами:

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

 

А + 0 = А,

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

<ЗЛЛ)

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

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

 

В рассматриваемой совокупности, наряду с оператором А, име­ ется также обратный (противоположный) оператор (—А).

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

Пусть у = Ах (х, у G R). Через х и у обозначим столбцовые матрицы, элементами которых служат координаты векторов х и у в базисе ер е2, ...» еп. Тогда у = А х, где А — квадратная матрица

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

АГ = ГА.

(3.1.2)

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

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

Рассмотрим в пространстве R два базиса

Г = (et ё2 ... е„) и

= (е/ е2 ... е^), связанные друг с другом соотношением

Г ,= Г 7 \

(3.2.1)

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

ра А в базисах Г и Г 2 соответственно. Тогда согласно (3.1.2)

АГ = ГА, АГ, = ГАР

(3.2.2)

Умножая второе равенство (3.2.2) справа наГ”1, получим с учетом (3.2.1)

АГ = %>ТАуТ~1.

Сравнивая последнее соотношение с первым равенством (3.2.2), на­ ходим

А - Т А у Т ' К

(3.2.3)

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

А х= Т~1АТ.

(3.2.4)

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

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

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

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

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

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

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

Оператор А называется особенным (неособенным), если |А | = 0 (соответственно | А | ^ 0 ). Если оператор неособенный, то

1)из Ах=0 следует х = 0;

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

В самом деле, если Ах = 0 , то в некотором базисе %

Ах = 0,

откуда, так как |А | ^ 0 , х = 0. Далее, пусть у — произвольный вектор пространства R, у — столбцовая матрица, составленная из координат вектора у в базисе а А — матрица линейного опера­ тора А в этом базисе. Так как А — невырожденная матрица, то для любой столбцовой матрицы у существует столбцовая матрица х, определяемая равенством

х — А~^у.

(3.3.1)

Отсюда у = Ах.

Полученному матричному соотношению соответствует вектор­ ное равенство

У — Ах ( x . y e R ) ,

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

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

данном базисе пространства R. Оператор А-1 также является ли­ нейным в R и

АА"1 = А_1А = Е,

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

у = Ах, х = А-1у.

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

Вектор х G R называется собственным вектором линейного оператора А, а число X е X — его собственным значением, если

Ах = Хх.

(3.4.1)

Выберем в R некоторый базис & = (е, е2 ... е„). Пусть А

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

а х — столбцовая

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

Ах = Alfjc = ЖАх, Хх = = 2>\х.

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

&Ах = &\х,

и, значит,

Ах = Хх.

(3.4.2)

Матричное равенство (3.4.2) в свою очередь эквивалентно системе алгебраических уравнений

(«п

4" &12 % 2 "Ь •••

 

а1пх п = О,

 

а21х 1 +

2 2

Ь )х2 + . . .

+

а2пх п =

0,

(3.4.3)

amx i + <*„2*2 +

••• +

*■)*„ =

0 •

 

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

 

11 х

«12

«1п

 

| А —Х£|

«21

«22 “ Х

«2п

= 0. (3.4.4)

 

а п1

а п2

«„„

Х

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

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

му корню X уравнения (3.4.4) соответствует ненулевое решение xlt х2, ..., х п системы (3.4.3), и, значит, числу X отвечает собствен­

ный вектор х = £ x fr = оператора А. Столбцовую матрицу х, составленную из чисел x v xv ..., хп — решения системы (3.4.3),

будем называть собственным вектором матрицы.

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

Пусть А1— матрица, отвечающая тому же оператору А, но при другом базисе в R. Матрицы А{ и А подобны: А} = Т~1АТ. Отсюда

A, - ХЕ = Т~1АТ — ХТ~1Т = Т~1( А - \ Е ) Т ,

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

| л х- Х £ | = | Л - Х £ | .

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

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

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

Отсюда видно, что сх ф 0 тоже является собственным вектором, от­ вечающим собственному значению X.

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

А(ах + Ру + ... + 6и) = аАх + РАу 4*... + 6Аи = = Х(ах + Ру + ... + 6и) (а, р, ..., д Е 5£).

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

Ле м м а 3 . 4 . 1 . Собственные векторы линейного оператора

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

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

АХ| = Х1х/

(/ = 1, 2,...» к; X, = Xj при i ф j). (3.4.5)

Допустим противное, а именно что в условиях леммы собственные

векторы хр х2, ...»xfc линейно

зависимы,

т.е. имеются

числа

а р а2, .... аАЕ

не все равны нулю и такие, что

 

 

а ,х 1 + а2х2 +

... + а*хл =

0.

(3 .4 .6)

Пусть, например, ак Ф 0.

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

(3.4.5),

 

к

 

2 п{Х{х. = 0,

(3.4.7)

i- i

 

Равенство (3.4.6) умножим на X, и вычтем затем из (3.4.7). Будем

иметь

 

к

 

Б аДХ,- - Х()х,. - 0.

(3.4.8)

12

 

Теперь равенство (3,4.8) умножим слева на А. Придем к равенству

к

 

Б аДХ, - X,)X,x, = 0 .

(3.4.9)

1= 2

 

Из равенства (3.4.9) вычтем равенство (3.4.8), умноженное на Х2. Получим

БаА ”

“ *г)х1= °-

i = 3

 

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

’ ^тН^/к ~

~ —i)x* =

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

множители отличны от нуля. Полученное противоречие доказывает лемму. ■

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

Линейный оператор А в л-мерном пространстве R может иметь не более чем л линейно независимых собственных векторов. Если характеристическое уравнение имеет п различных корней, то one-

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

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

Пусть А — оператор простой структуры и gp g2, ..., g„ — ли­

нейно независимые собственные векторы оператора А:

 

 

 

Agk = Xkgk

(k — 1, 2,

n),

(3.5.1)

Примем

эти

векторы

в

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

х = З’х, где

3’ = (g,,g2,

g

;i), а х — столбцовая матрица, эле­

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

то

У - Ах = AS?* =

(Ag, Ag2 ... AgRT) а = (Xjgj Xjg2 ... Х,£л)х = S?y,

где

( V i

У =

Хих„

п п

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

векторами gp g2, ..., gHс коэффициентами

Х2, ...,

Соотношения (3.5.1) эквивалентны одному матричному равен­ ству

A& = &At

где

О

V

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

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

А = КЛК~1.

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

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

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

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

х = х1 + х2 (Xj G R,, х2 е R2),

(3.6.1)

то говорят, что пространство R расщепляется на два подпростран­ ства R: и R2 или что пространство R разлагается в прямую сумму

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

R = R, + R2.

(3.6.2)

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

х = Xj + х2,

(3.6.3)

после вычитания (3.6.3) из (3.6.1) придем к равенству двух векто­ ров: Xj — iij е Rx и х2 х2 е R2, что невозможно, ибо у Rx и R2 нет

общего ненулевого вектора.

Т е о р е м а 3 . 6 . 1 . Если п-мерное пространство R расщеп­ ляется на два подпространства Rx и R2, т.е.

R = Rx + Rj»

то сумма размерностей Rj и R2 равна п.

Д о к а з а т е л ь с т в о . Выберем

некоторый базис flf f2, ..., f k в

подпространстве Rj и базис gp g2,

g, в подпространстве R2. Век­

торы fp f2, ..., fk, gA, g2, .... g/ линейно независимы. Действитель­ но, пусть

а 1*1 +

а2*2 +

+

a *f* + Pi8i + Р282 + •••+ Р/8/ = ®>

тогда

 

 

 

а Л +

a 2f2 +

••• +

aJk — ~ O 1S1 + Р282 + ••• + fyg/)*

Левая часть последнего равенства есть вектор из Rlt а правая часть — из Rj. Поскольку у непересекающихся подпространств R, и R2 общим является только нулевой вектор, то

a,fj + a2f2+

+ а АГА= 0,

ftgi + Р282 + — + Р/В/=

Отсюда,

в

силу линейной

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

f2,

f*,

a i = а 2 ==

= ал 0» а

из

линейной

независимости

векторов

g p g2i ...

» g i следует, что Pi = Р2 = ••• — Р/ = 0. Следовательно, век­

торы f p

f2, ...,

fk, g,, g2, ...,

g, линейно независимы. Так как они

векторы л-мерного пространства R, то

 

 

 

 

 

 

 

к + / «S л.

 

(3.6.4)

Рассмотрим

теперь

л

линейно

независимых

векторов

е,, е2, ...» еп пространства R. Каждый из этих векторов может быть представлен как сумма двух векторов из Rj и R2 и, значит, как ли­

нейная комбинация векторов f p ...» f*,

g , , ..., g,:

 

e i =

 

a i2^2 “b ••• “b

M

i

P1282

Pi/8/>

e2 =

**21^1

a 22^2 “I" •••

a2kfk

P2181

P2282 "b ••• "Ь P2/8/>

 

anlfi+

*

 

 

»

«

e « =

a n2f2 + ■•• +

a . tff* +

M

 

l + М 2 +

” •+ M l '

Применяя к системе векторов ер е2, ...» еп и fp ..., g/ лемму 2.2.2, будем иметь

л к + /.

(3.6.5)

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

к + 1 = п.

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

Пусть дано произвольное расщепление линейного пространства R на два подпространства S и Т: R = S + Т. Тогда любой вектор х G R раскладывается, и притом единственным образом (см. § 3.6), на сумму двух векторов из S и Т:

x = x5 + xr (x5 GS, xr GT).

(3.7.1)

Вектор x s называется проекцией вектора х на подпространст­ во S, параллельной подпространству Т. Аналогично, хг называет­

ся проекцией вектора х на подпространство Т, параллельной под­ пространству S.

Пусть Р — оператор, осуществляющий проектирование про­ странства R параллельно подпространству Т. Этот оператор опреде­ ляется равенством Рх = х5, где х — произвольный вектор из R, а

xs — его проекция на S. Очевидно, что Р является линейным опе­ ратором. Равенство (3,7.1) можно записать так:

X = Рх + хг.

(3.7.2)

Если х е S, то разложение (3.7.1), в силу единственности разло­ жения вектора из R на сумму (3.7.1), принимает вид x = xs (хт = 0), и из (3.7.2) в этом случае получаем xs = Pxs, т.е. оператор Р, примененный к вектору из S, действует как единичный оператор.

Если х G Т, то разложение (3.7.2) принимает вид

 

хт= Рхт+ х7,

 

и, значит,

(3.7.3)

Рхг = 0.

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

тор Р к обеим частям равенства (3.7.2). Будем иметь

 

Рх = Р ^ + Рхг.

 

Отсюда, учитывая (3.7.3), получим Рх = Р2х. Следовательно,

Р* = Р.

(3.7.4)

Оператор Р в R, удовлетворяющий равенству (3.7.4), называется

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

Произвольный проекционный оператор Р в R осуществляет про­ ектирование R на подпространство S = PR параллельно подпрост­

ранству

Т = R — S.

Действительно, множество

векторов x s — Рх

(х е R)

образует

некоторое подпространство S

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

Все пространство R можно рассматривать как прямую сумму двух подпространств S и Т = R — S. Произвольный вектор х G R раскла­ дывается на сумму

х = xs + х Т (xs G S, хтG Т).

(3.7.5)

Применим к (3.7.5) оператор Р: Рх = Pxs + Рхг. Но

 

Рх = Р2х = Р(Рх) = Рх5,

 

поэтому Рхг = 0 и, кроме того, Pxs = Рх = х5.

 

Квадратная матрица Р называется проекционной, если

 

Рг = Р .

(3.7.6)

Проекционному оператору в произвольном базисе отвечает про­ екционная матрица. Действительно, если Р — матрица, отвечаю­ щая оператору Р в базисе % — (е, е2 ... еп), то Р% = &Р. Но, с дру­ гой стороны,

Pg = p2g> - pg’p = g p \

Значит, %Р = &Р2, откуда и следует равенство (3.7.6). Докажем некоторые свойства проекционных матриц.

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