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

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

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

Равенства (6.7) можно представить в виде

А% = %(Х0ЕР + Нр),

где

g = (ег

е 2

ер),

~0

1

0

0“

0

0

1

О

 

1

0 0 0

О

Щр матрица сдвига порядка р).

Таким образом, оператору А в / 0 в базисе (6.6) отвечает матрица Х0ЕР+ Нр. Линейно независимые векторы elt

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

чек, взятых в каждом

из подпространств, можно составить

базис (жорданов базис)

в /?. В этом базисе

матрица опера­

тора А имеет жорданову нормальную форму

J = diag \Х1Ее1+ Я С1, ... » XsEl$

Hls).

Матрицы Л и У, отвечающие одному и тому же линейно­ му оператору А в разных базисах, связаны друг с другом соотношением подобия:

А = TJT~\

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

А = Т diag (Л»1г Х,а> . . . , Xs) Т *.

Таким образом, линейный оператор А имеет простую структуру тогда и только тогда когда, пространство R рас­ щепляется на инвариантные подпространства с линейными минимальными многочленами.

З а м е ч а н и е . Из подобия матрицы Л, соответствую­ щей жордановой матрице J, следует, что

det А — Х{‘ . . . A.J*.

§ 7. Инвариантные многочлены. Единственность нормальных форм линейного оператора

Через Dp (X) обозначим наибольший общий делитель всех миноров р-го порядка характеристической матрицы

ХЕ А (р = 1,2, ...,

п). В ряду

 

 

Dn(X),

Z?n_j(A»),

D,(X)

(7.1)

каждый многочлен делится на последующий без остатка. Действительно, минор /-го порядка можно разложить по элементам какой-либо строки. Каждое слагаемое этого раз­ ложения есть с некоторым множителем минор (/ — 1)-го порядка и потому делится без остатка на Dj-\ (X). Следова­ тельно, любой минор /-го порядка, а значит и D/ (Я), делится

без остатка

на

D /^i (X).

Т е о р е

м а

7.1. Наибольший общий делитель Dk (Я,)

миноров k-го порядка матрицы ХЕ А, где А матрица оператора в каком-нибудь базисе, не зависит от выбора базиса.

Д о к а з а т е л ь с т в о . Пусть А и А — две матрицы оператора А в разных базисах. Соответствующие характе­ ристические матрицы связаны друг с другом соотношением

ХЕ— А = Т~1(ХЕ А) Т.

Покажем сначала, что общие наибольшие делители мат­ риц ХЕ А и С (ХЕ А)у где С — произвольная невы­ рожденная матрица, совпадают. Пусть С = (ctj)t А = (ац). Тогда 1-я строка матрицы С (ХЕ А) имеет вид

т. е. является линейной комбинацией строк матрицы ХЕ А с независящими от X коэффициентами cilt с[2, ..., с1п. Поэто­ му минор матрицы С (ХЕ А) разлагается на сумму мино­ ров матрицы ХЕ А с некоторыми не зависящими от X коэффициентами. Следовательно, всякий делитель миноров

k-ro порядка матрицы

ХЕ А будет делителем

миноров

k-го порядка матрицы

С (ХЕ А). Точно так же

всякий

делитель миноров k-ro порядка матрицы С {ХЕ — А) явля­

ется делителем миноров k-ro порядка матрицы С-1 [С (КЕ

А)] = ХЕА. Это означает, что у матриц ХЕ — А и С {ХЕ А) общие делители миноров k-ro порядка {k = = 1, 2, п) совпадают.

Аналогичным образом устанавливается, что у матриц ХЕ А и {ХЕ A) Ct где С — произвольная невырожден­ ная матрица, общие делители миноров k-ro порядка совпа­ дают.

Применяя полученный результат к матрицам ХЕ А и

Т~х {ХЕ — А), а

затем к

матрицам Т~1{ХЕ А) и

Т~1 {ХЕ А) Ту

приходим

к заключению, что у матриц

ХЕ — Л и Т~х{ХЕ А) Т

общие делители (в том числе

иобщие наибольшие делители) миноров k-ro порядка {k =

1, 2, ..., п) совпадают. Теорема доказана.

Итак, Dk {X) {k = 1, 2, п) являются инвариантами линейного оператора и не зависят от выбора базиса в R. Разделив каждый член ряда (7.1) на последующий, получим

другую группу

инвариантов оператора А:

 

h

Р , «

 

 

№ =

P„-iW

о„_2М ’

 

 

 

 

 

O.

(О0 (X) = 1).

 

 

 

P . W

Многочлены ip (X) (р — 1, 2, ..., п) называются инва­ риантными многочленами. Произведение всех инвариантных многочленов равно характеристическому многочлену:

Д(Х) = | Х £ _ Л | = Д ,(Х )= Ш ,(Ч -

 

 

 

р=1

Разложим инвариантный многочлен ip {X) на неприво­

димые в поле Ж многочлены:

 

 

h (Ч = [Ф1 (*)Г‘ [Ф. № •

(Ф, № *

Здесь

cpi {X),

фа (А-)» ..., ф* (Я) — различные неприводимые

в поле

ЗС многочлены. Степени этих

многочленов 1фх (А,)]г‘,

[фа (А.)Г', ...»

[ф, {X)Ts, отличные от постоянной, называются

элементарными делителями характеристической матрицы ХЕ А или просто матрицы А .

Возьмем теперь в качестве матрицы оператора А матри­ цу, имеющую первую естественную нормальную форму:

А = d iag (Ли , Л22, . . . . Att),

где Ац (i = 1, 2,...» О — матрицы вида (6,4) с минималь­ ными многочленами

Ф* (X) = А**-J- ССл^ * ~f*

~f- Ct/ т^—lA»H“

 

 

(i = 1,

2,

Определитель

матрицы

 

 

/А ,

0

Р 1 ь.

II

1

 

\

0

0

• **> 0*

0

GCi mi

0

ОС/ т^—1

— 1

К ОС/1

как легко видеть, совпадает с ее минимальным многочленом:

IХЕт ~ А „ \= (Ц

( 1 = 1 . 2 .........

t).

(7.2)

Определитель квазидиагональной матрицы равен про­ изведению определителей диагональных блоков. Поэтому, учитывая (7.2), будем иметь

Д , (А,) » А (А,) =

Ъ (Я) ф2 А)

(X).

(7.3)

Вычислим далее D„_i (X).

 

 

Минор (п — 1)-го

порядка

элемента

матрицы

diag (ХЕт, — Лп , .... %Emt — Л„), расположенного вне диа­ гональных клеток, равен нулю. Действительно, для получе­ ния этого минора нужно вычеркнуть из определителя \КЕА\ строку и столбец, на пересечении которых расположен рассматриваемый элемент. Эти линии пересекают два раз­ ных диагональных блока. Пусть, например, вычеркнута строка, которая пересекает блок ХЕт{ Ац. Тогда верти­ кальная полоса, содержащая остаток блока XEmf Ац

и состоящая из т / столбцов, будет иметь только ту — 1 ненулевую строку. Разлагая рассматриваемый определи­ тель (п — 1)-го порядка по теореме Лапласа на миноры ту-го порядка указанной вертикальной полосы, убеждаемся, что он равен нулю (каждое слагаемое в этом разложении содер­ жит определитель, у которого одна строка сплошь заполне­ на нулями).

Рассмотрим теперь минор элемента, принадлежащего од­

ному из диагональных

блоков,

например, блоку ХЁт /— Ац.

Такой минор равен

 

 

 

 

 

 

Фх (Я)

Ф/-1(Я) Фж

(Я)

ф* (Я) х (Я),

 

где х (Я) — определитель блока, полученного

из блока

ХЕт. А// после исключения

строки и столбца,

на пересе­

чении которых расположен данный элемент.

 

 

В частности, минор элемента aimравен

 

 

(— l)w/-1 t i

(Я)

Ф/-1 (Я) фж

(Я)

ф, (Я).

Поэтому многочлен

 

 

 

 

'Фх (Я)

 

ф/_1(Я) ф/-|_! (Я)

М Я )

 

(7.4)

является общим наибольшим делителем миноров (п — 1)-го порядка, отвечающих элементам матрицы ЯЕт. Ац.

Так

как ф< (X) делится

без остатка на

ф/_1

(Я),

то

все

произведения вида

(7.4) (/

= 1, 2,

...,

t) делятся

без остат­

ка на произведение ф2 (Я) Фз (Я)... ф* (Я). Значит,

 

 

 

 

 

 

1 (Я) — ф2 (Я)

 

 

ф/ (Я).

 

 

(7-5)

Аналогичным образом

можно

получить

 

 

 

 

 

 

Dn- 2 (Я) =

фз (Я)

 

Ф/ (Я),

 

 

 

 

 

 

Dn—t+i (Я) =

фi (Я),

 

 

 

 

 

^

^

 

 

Dn-t (Я) =

=

 

(Я) ■= 1. 1

 

 

 

Используя (7.3), (7.5) и (7.6), находим

 

 

 

 

‘i W

=

~ ’I’1 W*

‘2W

=

D"_J (*,)

“ ♦» W ..........

tt w

=

°в ~ ^т =

’h w .

'<+■ w

=

=

^ w

=

»•

 

Таким образом,

многочлены ф< (Я)

(/ = 1, 2,

...,

/) сов­

падают с отличными от единицы инвариантными многочле­ нами оператора А, а отличные от единицы многочлены

[фА(Я)Г*, [<р* (Я) ] \

.... [ф*(Я))'й (k =

1,2...... s) в разложении

(5.13) совпадают с

элементарными

делителями оператора

А (и соответствующей матрицы А).

 

Отсюда следует единственность (с точностью до порядка расположения диагональных блоков) как первой и второй естественных нормальных форм, так и нормальной жор-

дановой формы

матрицы

оператора А ,

так как минималь­

ные многочлены

% (Я) (/

= 1, 2,

t)

с точностью до по­

рядка расположения диагональных блоков однозначно опре­ деляют все эти нормальные формы.

З а м е ч а н и е . Выше мы видели, что две подобные мат­ рицы имеют одни и те же инвариантные многочлены. Но тог­ да эти матрицы подобны одной и той же нормальной (на­ пример, жордановой) матрице и потому подобны друг другу. Таким образом, справедлива следующая

Т е о р е м а 7.2. Для того чтобы две матрицы с эле­ ментами из поля di были подобны, необходимо и достаточно, чтобы они обладали одними и теми оке инвариантными многочленами.

Г л а в а V

ПРЕОБРАЗОВАНИЕ МАТРИЦЫ

ККВАЗИДИАГОНАЛЬНОМУ ВИДУ

ИРАЗЛОЖЕНИЕ ЕЕ НА СОСТАВЛЯЮЩИЕ

§ 1. Дефект матричного многочлена

Если

в многочлене

 

 

f (^) “

«(Лт 4~ сс1Хт~^ -|-

-}- ат—j?i -Г сст

(а/ £

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

f (^) — а0^ 4"

4~

4- a m —i A -f ат£л.

Вычислим дефект этого матричного многочлена. Будем считать, что Di — поле комплексных чисел. Пусть

 

J = diag

J2 (^2)* • • •»

(К)}

— жорданова матрица, подобная матрице А, где

 

 

O^i) = ^ i ^ r j 4~ H r h

 

г/ — степень элементарного делителя,

отвечающего жор-

дановой клетке Jj (X/) (порядок блока Jj (kj)).

Тогда вследствие того, что А = TJT~\

 

f(A)=Tf(J)T~'.

(1.1)

Здесь

f(J) = diag I/(Л), / W ,

/ W ) -

Согласно (1.1) дефект матрицы f (А)

равен дефекту мат­

рицы / (J). Дефект же квазидиагональной матрицы / (/) ра­

вен сумме

дефектов

диагональных блоков. Поэтому, обо­

значив через dj дефект блока / (У/), будем

иметь следующее

выражение для дефекта матрицы f (Л):

 

d = 2 4

( 1. 2)

/-=1

Вычислим дефект dj матрицы f (J/). Имеем

А=

л= i ; W r W ,

д=о

(1.3)

А =

|*=0

т

J? =

В-0

Здесь

с* = * (* - 0 •

H + l)

(ft= 1,2, . . . . т)

биномиальные коэффициенты.

Сложив равенства (1.3), умноженные соответственно на

числа

а1П, а т _ !........а 0, получим,

учитывая,

что

 

=

0

 

при \i

>

г,* (г/ — порядок матрицы Я Г/),

 

 

 

 

 

/№ ) = П Ч

Ег, +

г/

1

 

 

 

 

 

 

 

Е < “ ^ Г

Ц +

 

+

 

 

 

 

 

 

 

 

 

Ц=1

 

+ а т _ |1с{;ХЙЯ’;..

(1.4)

 

 

 

 

 

 

 

 

 

 

С другой стороны,

 

 

 

 

 

 

 

 

=

/ (fl) (^/) =

а0т (т 1 )

р. -

j - К?“1

)

м

4

-

 

+ а i (m — 1) ( т

— 2)

(т — р)

+

 

 

 

 

или

 

 

 

 

 

 

 

4* ат—др (р — 1)

 

П °/,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с(Ю

=

K

d

^

r 4*-1

 

4 -

4 - <

w

! M

^

Г (Ь/)

 

 

 

 

 

 

 

 

 

 

 

 

(1.5)

 

С учетом

(1.5)

равенство (1.4)

принимает вид

 

 

 

 

 

 

 

 

 

 

 

Г/-1

 

 

 

 

 

 

 

 

т

-

т

в . ' +

ъ

 

 

 

 

 

 

ц—J

(2. 1)

§ 21

ТЕОРЕМА

ГАМИЛЬТОНА — КЭЛИ

 

99

ИЛИ

 

 

г(ад

 

 

 

 

 

 

 

/(а д

 

/ < ' / - » » / ) -

 

 

 

II

 

(/7-1)1

 

 

/

Ш =

0

/(ад

 

/‘Ч -21(V)

(1-6)

 

(>7-2)1

 

 

 

_

0

0

 

 

/Я ,)

_

 

Обозначим через V/ кратность

Я/ как корня многочлена

/ (Я) (в частности,

v;- может быть

равна и нулю). Тогда

/ (ад =

/' (ад -

 

=

(ад = о,

/^/» (ад ф о.

Поэтому дефект матрицы / (J/) (см. 1.6)) равен v/} если

V/ < г/, и равен г/, если

v;- >

/7 . Итак,

 

 

 

 

 

d/ = min(v/,

гу),

 

 

 

и, принимая во внимание (1 .2), приходим к следующей тео­ реме.

Т е о р е м а

1.1. Дефект матричного многочлена /

(Л),

где А матрица с элементарными делителями

 

^i)rS (Л>

Я2)г*,

(^ W 5>

 

определяется формулой

 

 

 

 

d =

£ m in (v /f

о).

(1.7)

/- 1

Здесь у,-— кратность Я/ как корня / (X).

§ 2. Теорема Гамильтона — Кэли

Пусть Д (Я) — характеристический многочлен матрицы А порядка п. Кратность v;- корня Я/ характеристического многочлена не меньше, чем степень любого элементарного

делителя матрицы А вида (Я — Я/)г/. Учитывая это, соглас­ но (1.7) будем иметь

= 2 г,. /=i

Но сумма степеней элементарных делителей квадратной матрицы порядка п равна п. Поэтому d = п, а это означает, что А (Л) — нулевая матрица, т. е.

Д (Л) = 0.

Итак, всякая квадратная матрица удовлетворяет своему характеристическому уравнению (теорема Гамильтона — Кэли).

Эта теорема непосредственно следует и из соотношения (4.7.3). В самом деле, характеристический многочлен А (Я) оператора А в /?, которому в некотором базисе отвечает рассматриваемая матрица А , содержит в качестве множи­ теля минимальный многочлен всего пространства ifo (Я). По­ этому из равенства

(А) = 0

немедленно вытекает равенство

А (Л) = 0

и, значит,— равенство (2 .1 ).

§ 3. Построение матрицы, преобразующей квадратную матрицу к квазидиагональному виду

Рассмотрим квадратную матрицу U порядка п с собст­ венными значениями Я!, Я2, ..., Я„, среди которых могут быть и равные *). В дальнейшем эти собственные значения будем разбивать на группы, отмечая принадлежность к той

или иной группе верхним

индексом. Так,

Я^0) обозначает

/-е собственное значение группы о.

 

 

 

Пусть собственные значения матрицы U разбиты на р

групп Я1°\ ...» ЯЦ

(а =

1,

...»

р; 2

=

п)> так> что

 

 

 

 

<7=1

 

 

 

—x f|> 0

( а ф т

i

- l

....... А„; / =

1,

k s).

 

 

 

 

 

 

 

(3.1)

Каждой группе

о поставим

в соответствие

многочлен

Д с ( Ч = П П

(X —

A f )

(<Т = 1 , 2 ,

 

 

р ) . ( 3 .2 )

/-=1

 

 

 

 

 

 

Вычислим ранг матрицы Aff (U). Для этого воспользу­ емся теоремой о дефекте матричного многочлена, согласно

*) Собственное значение кратности т здесь рассматривается как т равных собственных значений, и каждому из них приписывается свой индекс.

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