Файл: Абгарян К.А. Матричные и асимптотические методы в теории линейных систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 178

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

174

 

С И СТ Е М Ы Л И Н Е Й Н Ы Х

У Р А В Н Е Н И Й

 

[ГЛ . V II

Пусть

 

 

 

 

 

 

 

 

 

 

JA =

d i a g ^

(Я,), 4 Л) ( U

. . . .

J{PA) ß p)),

 

 

JB =

diag (j\B) а д , 4 B) (m), . . . .

4 B> Ы) .

 

Здесь

Я2,

 

Xp — собственные

значения

матрицы А, а

Pi, р2>

Р?— собственные

значения матрицы

В. Через

kt обозначим порядок клетки Жордана jjA){kt),

через

1{

порядок клетки Жордана 4 В)(мД-

Матрицу

Z

и матрицу

R = Т - 1 CS представим в виде блочных матриц:

 

 

г11

гХ2

• ■■

г,Д

 

/''и

г12

. • •

/V

 

Z =

Z21

Z22

• '■•

г* Н ,

я =

Г2Х

Г22

• ..

r2q

 

 

Zp\

 

• ■••

ZpJ

 

\\гр\

Гр2

..

ГРЧ

 

 

— матрицы типа k( X

/у.

 

 

 

 

 

Ги ■

 

 

 

 

 

 

 

 

При этом уравнение (7.5) расщепляется на тп матрич­

ных уравнений:

 

 

 

 

 

 

 

 

 

 

 

4 Л)(4 )* / =

ztjJ}B) fr,) + Гц

 

 

(7.7)

 

( і =

1, 2, . . . , р \

/ =

1, 2 .............q).

 

 

Каждое из матричных уравнений системы (7.7) представ­

ляет собой уравнение

типа

 

 

 

 

 

 

 

+ Я*) и =

и (ц£, + Я,) + г

(Я, =зь ц),

(7.8)

где Et — единичные матрицы порядков k и I соответ­ ственно, Hk, Hj — матрицы сдвига порядков k и I соответ­ ственно, и и г k X /-матрицы.

Для доказательства леммы достаточно показать, что уравнение (7.8) имеет единственное решение. Пусть

и = (uij),

г = {гц)

(і = 1,2,

1,2, . . . , / ) .

Тогда матричное уравнение (7.8) эквивалентно следующей системе скалярных уравнений:

( к — ( А )

U i j +

« 1 + 1

/ —

Щ і -

1 =

г и

 

(7.9)

(/= 1,2, . .. ,

А;

/ =

1,2,

. . . ,

/),

 

причем

 

 

 

 

 

 

 

(7 .1 0 )

Ир+і/ =

0

( / =

1,2,

. .. ,

/),

I

«ю = 0

( / = 1 , 2 ,

 

 

k). J


1 г]

Т Е О Р И Я

В О З М У Щ Е Н И Й

175

Разрешая

(7.9) относительно ис!, имеем

 

 

иц =

г ц +

и і / - 1 ~ и і+ 1 І

(7.11)

 

 

X — (i

 

 

 

 

Рекуррентное соотношение (7.11) однозначно определя­ ет все U[j. В самом деле, учитывая (7.10), находим

Up\

X ц

Зная Upi, далее последовательно при

і = р — 1, р — 2, ...

..., 1 определяем значения всех элементов первого столбца

матрицы и, пользуясь соотношением

 

«п

гі\ ~ ці+ і 1

X — ji

 

Затем строим элементы второго столбца, начиная с эле­ мента Upг, посредством соотношения

 

ГІ2+ ип ~ и1+ 12

=

Р, р — 1,

р — 2

1)

W/2 =

А, — р.

 

 

 

 

 

И т. д.

Этот процесс подтверждает существование и единствен­ ность решения матричного уравнения (7.8), что в свою оче­ редь доказывает лемму.

7.3.Асимптотический метод для однородной системы.

Пусть К = 1 Кг КР) — матрица, преобразующая мат­ рицу А к квазидиагональному виду Л = diag (Alt Л2, ..., Ар) при условии, что у матриц As и Л0 (s ф а) нет общего соб­ ственного значения (в частности, Л может быть матрицей

Жордана).

решение

однородного уравнения

(7.1)

Будем строить

(h (t) =э 0) в виде

 

 

 

 

* =

2,КоУо,

(7.12)

 

 

0= 1

 

предполагая, что

уа — решения уравнений

 

~ ~

= АоУо

(о = 1.......... р),

(7.13)

а постоянные матрицы Ко, Аа представлены формальными рядами

Ко = 2 &kK[o \

~Ао = 2 е*МЧ

(7.14)

4=0

4=0

 


176

 

С И СТ Е М Ы Л И Н Е Й Н Ы Х

у р а в н е н и й

[ГЛ. ѵ п

Подставим (7.12) в однородное уравнение (7.1) и исключим

ЛУа

с помощью равенства (7.13).

Получим

 

di

 

 

 

 

 

 

 

 

 

 

У. КаАоУа -f- eß)

2

КаУо-

 

 

<т=1

 

 

а=I

 

 

Это соотношение будет выполняться тождественно, если

 

-f eß) Ко =

КаАа

=

1..........р).

(7.15)

Подставляя сюда ряды (7.14) и приравнивая члены, со­

держащие е в одинаковых степенях, имеем

 

 

АК™ =

 

 

 

 

 

 

 

А ф =

/СУ]АГа] +

Ф А[о1] -

ß/C[o0],

 

(7.16)

 

А ф ] = K[a]Alo]+

К1а]ф + Ф

- В К [о'\

 

Положим

 

 

 

 

 

 

 

 

Ф =

Ко, Аса0] = Ао

(а =

1,2,

. . . , р).

(7.17)

При этом первое равенство (7.16) выполняется тождест­ венно. Из остальных равенств последовательно могут быть

определены Ка \ А д іф ,

A ^;

... В самом деле,

пусть уже

найдены /(ст], Ag] — 0,

1,2,

.... k

— 1).

Определим I\lak],

используя (k + 1)-е равенство

(7.16),

которое с уче­

том (7.17) можно представить так:

 

 

 

АК1а] = К[о ]Ао +

КоА[а 1+ D[ok-U

(7.18)

Здесь D[k- ' ] =

А'о^Ло'1-11+ • • • +

 

-

ВК[о~и -

уже известная

матрица.

 

 

 

 

 

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

К~' = М =

Получим

AQlok] =С ф А а +М К оА Ікі + M D [k- n

(7 .1 9 )

где

 

QW = MKlo].

 


§ 7 ]

ТЕО РИ Я ВОЗМУЩ ЕНИЙ

177

Ввиду квазидиагональной структуры матрицы Л ра­ венство (7.19) распадается на р независимых равенств

ЛД о =

Qlo Ла + М'КаЛ!ьк] + МДак- ' ]

(s = 1,

. . . , р).

При s =

а

МаКа — Eka, и мы имеем

 

(7.20)

 

 

Отсюда

 

Л„<2аа =

QaaAa +

+ ЛСД

 

 

 

 

 

 

 

 

ЛС„*] = -

М Д * “ *1+ AoQao1-

QroAa-

(7.21)

При s Ф a

MsKo =

0, и потому

 

 

 

 

AsQs*] = Q^a’Aa + М Д Д " .

(7.22)

Так как Л5 и Лст не имеют общих собственных значений, то матричное уравнение (7.22) по лемме 7.1 имеет единствен­ ное решение.

Неопределенной осталась лишь матрица Qoa1Так как не осталось никаких невыполненных условий, то в качест­

ве Qaa можно взять произвольную постоянную матрицу порядка ka (ka — порядок блока Аст) и, в частности, можно

принять QIJC =

0. Зная

легко

определить

До^'-

 

К\?] =

AQ[-ft] .

 

(7.23)

Полученные

рекуррентные

соотношения (7.21)

и. (7.23)

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

ные значения

..., Кп. При

этом согласно

(7.22)

 

 

 

[*і

л Д * ~ 1]

 

 

 

 

Qso

 

 

и (7.23) принимает

вид

 

 

K W

=

у

K sM s

D [fe-n + K o Q m

_

—J As — Ka

$=1 sv»(J

Упрощается и выражение (7.21). Так как теперь Л0, Qoo — скалярные величины, то

ЛЕ,*1 = Х Р = — M CTD [aft- 1 ] .


178

СИСТЕМЫ ЛИНЕЙНЫ Х у р а в н е н и й

[ГЛ . VII

Приближенное решение однородного уравнения (7.1) можно получить из формул (7.12) и (7.13), удерживая в формальных рядах (7.14) некоторое число первых членов:

* „ = а-2=1 к ? Ѵ .

^ h - A S ” Ѵ Г ,

где

 

К ст)( = 2 BW

,

Л'Г = 2 е ^ 1.

= 0

 

 

 

/2=0

 

 

 

 

Поскольку Аат) — постоянная матрица,

Уіт) =

,(т).

 

 

е

С(7

 

 

*,„ =

2

r(m)

A<m>f

 

в

C-

 

0=1

 

 

 

7.4. Асимптотический метод для неоднородной системы. Решение неоднородной системы (7.1) будем строить в виде

 

 

* = j^KoiJo,

 

(7.24)

где t/a — решение

СТ=1

 

 

 

уравнения

 

 

 

 

= А ауа + Mah (t)

(а = 1 , 2 .......... р).

 

(7.25)

После подстановки (7.24) и (7.25) в (7.1) имеем

 

2

КЛаУа + 2

KoMoh (0 =

(А + гВ) 2

+

Л (Q.

0=1

0=1

 

0=1

 

 

Отсюда, приравнивая коэффициенты при уа и свободные члены, получим соотношение (7.15) и равенство

2koMah(i) = h(t).

0=1

Чтобы выражения (7.24) и (7.25) представляли собой

формальное решение уравнения (7.1), нужно Ко и Л„ опре­

делить, как это указано выше, а Ма выбрать так, чтобы выполнялось равенство2

2КоМо = Е.

(7 .2 6 )

0=1