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

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

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

$ 7]

определению некоторого числа квадратных матриц меньшего порядка.

Исходя из указанного в начале данного параграфа раз­ биения собственных значений матрицы U на группы, мож­

но построить матрицу К = (Кц К2

КР), преобразующую

матрицу U к квазидиагональ ному виду (6.5). При этом груп­

пе s (s = 1 , 2 , ..., р), состоящей из

ks равных собственных

значений Jis, будет отвечать диагональный блок As матрицы Л с ks кратным собственным значением A.s.

Итак, наряду с (6.2) имеет место легко реализуемое пред­

ставление

U =

КШ .

(6 .6)

 

При этом блоки матриц Л и / с одинаковыми индексами

подобны между собой, так что

 

Л. = NJ, (\)

( 5 = 1 . 2 .......... р).

(6.7)

Сравнивая (6.2) с (6 .6), находим

 

 

Т = KN,

(6 .8)

где N = diag (Nlt

W2, . . . .

Np).

 

Таким образом, использование преобразования, приво­

дящего данную

матрицу U к квазидиагональному

виду,

позволяет свести задачу по построению преобразующей мат­

рицы Т гс-го порядка к более простой задаче: к

построению

удовлетворяющих равенствам (6.7) матриц Ns (s = 1,2,

..., р) меньшего порядка.

 

Если кратность собственного значения

невелика

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

AsNt = N J t (K).

§ 7. Случай матрицы простой структуры

Из теорем предыдущего параграфа непосредственно сле­ дуют два результата.

Если /га-кратному собственному значению матрицы от­ вечают в форме Жордана ka клеток Жордана первого по­ рядка, то

I*.® = 4 1».

и обратно.

Если U — матрица простой структуры порядка п, то

2 й!" = 2 4 ° = п,

S—1 S=1

и обратно.

Приведем еще некоторые предложения, касающиеся матриц простой структуры.

Л е м м а 7.1. Пусть собственные значения матрицы U разбиты на р групп при условии (3.1) и tip какое-нибудь

из собственных значений,

включенных в группу о.

Если

(U -

Xf'Ej До (V) = 0,

(7.1)

то 1) все собственные значения группы а равны между собой;

2)матрица Л а является диагональной матрицей.

До к а з а т е л ь с т в о . 1) Из (7.1) следует, что все ненулевые столбцы матрицы Да (U) являются собственными

векторами, соответствующими собственному значению tip. Так как ранг матрицы Дс (U) равен ka, то собственному значению tip соответствуют ka линейно независимых

собственных векторов. Это значит, что tip является по крайней мере 6а-кратным собственным значением матрицы.

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

tip в точности равна ka.

2) Представим равенство (7.1) так:

(U - X f E n) КаМм = 0.

Отсюда, так как среди миноров

матрицы Мао ранга ka

имеется минор ранга ka, то

 

 

(U - Я}°>£„) Ко =

0.

(7.2)

Используя последнее соотношение, получим

Л„ = MaUKa = MoXf'Ko =

Х } ^ .

Лемма доказана.

Т е о р е м а 7.1. Пусть собственные значения матрицы U разбиты на р групп при условии (3.1), причем в каждую группу включены только равные между собой собственные^ значения. Для того чтобы V была матрицей простой

структуры, необходимо и достаточно, чтобы

(U -

1?Е„) Д5 (У) =

0

(s =

1,2..........р),

(7.3)

где Kj собственное значение группы s.

 

Д о с т а т о ч н о с т ь .

Согласно

лемме 7.1

все мат­

рицы As (s =

1, 2, ..., р) имеют диагональную структуру.

Значит, матрица U приводится к диагональному виду (при

этом, кстати,

преобразующая матрица К в силу равенства

(7.2) состоит из собственных

векторов матрицы U). Следо­

вательно, U как матрица, подобная диагональной, явля­

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

 

 

Н е о б х о д и м о с т ь .

 

Пусть U — матрица

простой

структуры и

она подобна

диагональной матрице

Л. Эту

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

Л = diag(Л1э Ла, . . . , Лр),

где As— диагональная матрица, на главной диагонали ко­ торой расположены собственные значения матрицы U, вклю­ ченные в группу s. Пусть К — матрица, преобразующая V к диагональному виду Л, и М К~“1:

 

 

U = КАМ.

 

 

(7.4)

Рассмотрим

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

 

 

 

 

 

 

Л(£/) =

П ( * / - У ? Я),

 

 

 

 

s=l

 

 

 

где Xlt

...» Хр— все различные

собственные значения

матрицы U (Я, =

X\s) =

=

 

 

 

С учетом (7.4) имеем

 

 

 

 

 

l ( U ) = K U ( A - X sEn)M.

 

 

 

S = 1

 

 

 

Здесь

 

diag (Ai

\ sEklt

 

 

Ар Х3Е ^ ,

Л

=

 

. . . ,

Так как

As =

XsEka (s = 1, ..., /?), то

s-й диагональный

блок матрицы Л — XsEn — нулевая

матрица. В силу этого

 

 

П (Л —

= 0,

(7.5)

и потому

А (£/) = 0.

Матрица А (U) является множителем левых частей всех равенств (7.3), что и подтверждает справедливость этих ра­ венств. Теорема доказана.

Равенство (7.5), как мы видели, является необходимым условием того, чтобы U имела простую структуру. Можно показать, что (7.5) является в то же время и достаточным условием.

Действительно, из (7.5) следуют равенства

 

П (Ло

XsEka) =0

(о = 1,2,

р).

 

S=1

 

 

 

 

Но при а Ф s матрицы Ла — \Е ь а являются

невырожден­

ными,

значит,

 

Ла ^aEka о.

 

 

 

 

 

Следовательно, диагональные блоки матрицы Л, по­

добной U, являются

диагональными матрицами, а это озна­

чает,

что U — матрица простой структуры. Таким образом,

имеет

место и

такая.

Х2, ..., Хр все различные

Т е о р е м а

7.2.

Пусть

собственные значения матрицы V . Для того чтобы матрица U имела простую структуру, необходимо и достаточно, чтобы

Д (У) = п (U — >,£„) = 0 .

5=1

С л е д с т в и е Если U матрица простой струк­ туры с одинаковыми собственными значениями, равными X, то непременно

U = ХЕ.

З а м е ч а н и е

1.

Если U — матрица

простой струк­

туры, то

 

 

 

 

 

А(Х) = П ( Ь - К )

 

 

 

 

s—\

 

 

является ее минимальным многочленом.

 

 

З а м е ч а н и е

2.

Равенство A (V)

=

0 является необ­

ходимым и достаточным

условием линейности всех элемен­

тарных делителей

характеристической

матрицы ХЕ V.

§ 8. Разложение квадратной матрицы на составляющие

Равенство (3.9) после перемножения блочных матриц /С, Л и М можно записать в виде

и =

(8 .1)

Ов1

 

где

 

Uа ==/СоЛдЛ^о’.

 

Квадратные матрицы Ua (а = 1, 2,

р) будем называть

составляющими матрицы U. Каждая составляющая Ua

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

составляющей

UQ\

остальные п ka собственных

значе­

ний матрицы U

(ka — число собственных значений

матри­

цы U, включенных

в группу а) равны нулю. Это

следует

из того, что матрица MUaK, подобная матрице Ua, явля­ ется квазидиагональной матрицей, в которой диагональным

блоком с

номером а,

как и в матрице А, служит ka X ka-

матрица

Ла, а все

остальные

диагональные блоки — ну­

левые.

 

 

 

Составляющие матрицы U, как нетрудно проверить, ис­

пользуя

(3.5), взаимно ортогональны:

 

UaUs = 0

(о Ф s),

и, кроме того, удовлетворяют еще следующим соотношениям:

U'S = КаКМ„ (т = 1 , 2 , ...),

UU„ = UJJ = U*.

Как уже отмечалось, в выборе матриц /Со (сг = 1 , 2 ,...»р) при данном разбиении собственных значений матрицы U на группы имеется определенный произвол (см. § 5). Со­ ставляющие же при данном разбиении собственных чи­ сел на группы определяются однозначно; они инвариантны относительно произвола, имеющегося в выборе /(с- Дейст­ вительно, принимая вместо Ка матрицу KaNa, будем иметь

й„ = КАаМо = KnNaN7'KNaN^M a = К Л Ж = V,.

Пусть f (А,) — многочлен скалярного аргумента Л. Если две матрицы А и В подобны и матрица Т преобразует А в В,

то матрицы / (А) и f (В) также подобны, причем та ж е мат­

рица Т преобразует f (Л) в f (5).

К преобразует U

 

Матрицы U и Л

подобны, и матрица

в

Л. Поэтому f (U)

и f (Л) подобны и К

преобразует f (U)

в

/ (Л), т. е.

 

 

f(U) =

Отсюда непосредственно вытекает следующее разложение матричного многочлена f (£/) на взаимно ортогональные со­ ставляющие:

/((/) - % K a f ( A a ) M a.

(8.2)

ff=l

 

Разложение (8.2) справедливо при произвольном разбие­ нии собственных значений матрицы U на группы согласно условию (3.1).

Посмотрим, какой вид приобретает формула (8.2) в том частном случае, когда в каждую группу объединены только равные между собой собственные значения.

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

Д а(U) = п ( £ / —

Я

Д0 = . )1 , 2 , . . . , ( / » )

S—1

 

 

— матрицы ранга 1 и потому каждая из этих матриц раз­ лагается на произведение п х 1-матрицы Ка на 1 х X л-матрицу Мое:

Дa(U) ** КаМоа.

Преобразует (8.2). Учитывая, что в данном случае А а = = Ха и, следовательно, f (Ла) = f (Яо) — скаляр, и исполь­ зуя (3.4), будем иметь

/ ( ( / ) = ^К оН Х а)Ма =

а=1

= у f (К) К, ( М о Л Г 1Мы —

j \ f (Ха)

,

а=1

 

o d

0

или

 

 

 

/ д о * 2 / < м

МоаКа ■

 

(8.3)

<т=1

 

 

 

Вычислим М о а К а .

Матрицу, присоединенную к характеристической мат­

рице ХЕп U, обозначим через F (Я). Имеем

 

F (X) (ХЕ„ U) = (ХЕ„ -=V) F (X) = Д (I)

(8.4)

где

А(Х) = П(Х — V

S=1

характеристический многочлен матрицы U. Прежде всего покажем, что

 

Да (U) = F (^а)-

 

(8.5)

 

Положим

 

 

 

 

 

 

Л (Я)

Д ( х ) _

/ч,

х).

(8.6)

 

Х ^с

=

ф

 

 

 

В этом равенстве заменим X на С/, а х на

ХЯ^чогда,

так

как Д (V) = О,

 

 

 

 

 

 

( X E „ - U ) < p ( U , Я) = Д ( Х ) Е „ .

 

 

Отсюда видно, что ф (U, X) есть матрица,

присоединен­

ная к матрице ХЕп U.

 

 

 

 

 

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

из

равенства

(8.6), если положить

х — Ха, получаем

 

 

 

 

 

 

Ф(Я, х а) =

 

= д„(X),

 

и,

значит,

 

Ха) = ДaiJJ),

 

 

q>(U,

 

что и доказывает справедливость равенства (8.5). Продифференцируем равенство (8.4) по X и в полученном

соотношении положим X = Ха- Будем иметь

(ХаЕп U)

d F (Я)

Н“ F (^а) — Да (Ха) Вп,

 

dk

 

Умножим последнее равенство слева на F (А,а) и,

учиты­

вая, что F (Ачт) (ХаЕп— U) =

0 (это следует из равенства

(8.4)), получим

F2(Xa) = ba(Xo)F(Xa),

 

 

 

или

Д?(У) =4<,(11)Да(У).

(8.7)

 

С другой сторойы, так как в данном случае МооКа

скаляр,

 

 

 

Да (U) = КаМоаКаМоа = MQaKa^o(U).

(8.8)

Сравнивая равенства (8.7) и (8.8), получаем До (Я.о) — МоаКа>

С учетом полученного соотношения (8.3) приобретает вид

пА (И\

/ ( t O - S / W A T O -

<8'9>

<т=1

 

что представляет собой формулу Сильвестра для случая, когда все собственные значения матрицы U простые. Из соотношения (8.9) путем предельного перехода может быть получена общая формула Сильвестра, справедливая и при наличии у матрицы U равных друг другу собственных зна­ чений:

 

1

 

о = 1(*а -1 )!

Здесь

обозначает общее значение всех ka равных собст­

венных значений матрицы U. включенных в группу а (ta =

=

*

1Г= 1 .,

. .2 .

. kc). .

к

квазидиагональному виду и

 

П р

и м е р.

Приведем

разложим на

составляющие матрицу

 

 

 

 

 

— 8

2

5

— 4

\

 

 

 

U =

— 15

5

6

— 3

1

 

 

 

 

— 12

4

9

— 10 Г

 

 

 

 

— 6

2

5

— 6

 

 

а) Собственные значения этой матрицы А,х = 1 , Х2= —1,

А,3 = 2, Я4 = —2 разобьем, например, на следующие две группы:

группа

1:

J40 - 1 .

ЧР.—

1,

 

группа

2:

Я(,а -■= 2,

Я(2а =

— 2;

 

 

 

 

 

— 6

6 — 3 0 \

Ai (У) =

(U

Я1а £ 4) (У -

=

— 9

9 - 6 3 1

— 12

12 — 9 6 Г

 

 

 

 

— 6 6 - 3 О /

§ 8] РАЗЛОЖЕНИЕ КВАДРАТНОЙ МАТРИЦЫ НА СОСТАВЛЯЮЩИЕ!19

Примем

д2(U) = (U - W E ,) (U - W E,) =

3

б

9

12

12

12

 

 

6

б

Примем

Преобразующая матрица

К= (К,Кг) =

м= к ~ ' =

Отсюда

— 2

1

 

 

1

2

 

1

2

 

М* =

— 1

1 0 - 1

 

- 1 О

 

О 3

 

 

 

и, далее,

 

 

 

 

 

 

 

8>

A j= м р к \ = ( _^

_ ; 2;

 

л„ =

м ,и/с,

 

 

 

=

( 2

2)

 

 

 

 

 

 

 

\0

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

матрица

К преобразует

матрицу U к

виду

 

7

12

0

0

 

 

 

 

 

 

 

 

 

Л =

— 4 — 7

0

0

 

 

 

 

0

0

2

8

 

 

 

 

 

 

 

 

 

 

0

0

0 — 2,

 

 

 

Составляющие матрицы U:

 

 

 

 

 

'2

 

7

 

Ux =

К1А1М1=

 

10

 

 

13

 

 

 

 

 

 

 

 

 

 

7

 

 

 

'— 10

4

— 2

8'

£/а =

К Л Ж

— 18

8

- 4

14

=

8

— 4

12

 

 

— 16

 

 

- 8

4

- 2

6,

б)

Приведем теперь матрицу

U к диагональному виду.

Для этого

разобьем

собственные

значения

на

четыре

группы:

 

 

 

 

 

 

 

 

Х‘" =

1,

Х<2>=

— 1, Я,3| =

2,

Л(4>=

— 2.

Д4 (У) =

(У — Я<2)

(У — Х,3,£^) (У —

 

=

 

 

 

 

 

 

— 12

12

— 24

36

 

 

 

 

— 18

18

— 36

54

 

 

 

 

— 24

24

- 4 8

72

Примем

 

 

 

— 12

12

— 24

36

 

 

 

 

 

 

 

 

* , =

Да (У) = (У - Х '% ) (У -

Примем

/ С , -

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