книги / Матричное исчисление с приложениями в теории динамических систем
..pdfГЛАВА 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) |
1—2 |
|
Теперь равенство (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). Докажем некоторые свойства проекционных матриц.