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

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

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

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

Квадратную матрицу U (т), собственные значения кото­

рой

разбиты

на р групп Я{а),

 

 

^<т =

1, 2,

р;

2

ka = п\ так, что на промежутке О-Ст <; L выполняют­

с я

 

}

 

 

 

 

 

 

 

 

 

ся

условия

 

| № -

 

 

| >

0

 

 

(1.1)

 

 

 

 

 

 

 

 

 

( а ф у , 1 =

1 , 2 , . . .

,

k Q\

/ =

1 , 2 , . . .

,

fcs),

 

при каждом фиксированном т £ [О, L1 можно

представить

так

(см. гл. V):

 

 

 

 

 

 

 

 

 

 

 

 

a = = 2 K c A ffMff|

 

 

(1.2)

 

Ко,

 

 

0=1

 

 

 

 

 

п X

где

Ла,

М0 — матрицы

типа

соответственно

ka X k0t

ka X п,

удовлетворяющие

равенствам

 

 

 

 

 

MaKs =

Ekat

s — о*

 

 

(1.3)

 

 

 

 

 

 

О,

Бфи

 

 

 

(Ег — единичная

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

 

 

 

Собственные значения матрицы Л0 суть собственные зна­

чения матрицы U, включенные в группу а.

 

 

 

Введя блочные матрицы

 

 

 

 

 

 

 

 

 

 

 

л

 

 

 

 

 

 

'ЛГ,

 

 

 

 

л =

 

 

 

 

 

 

K = ( K t

 

Ю ,

 

О

°

м -

ж

 

 

 

 

 

 

 

 

 

 

 

будем также иметь

 

 

 

р—

 

 

 

 

 

 

 

 

 

(1.4)

 

 

 

 

U = KAM,

 

 

 

 

 

 

 

МК = КМ = Еп.

 

 

(1.5)

В качестве матрицы /Со может быть взята любая матри­ ца, составленная из ka линейно независимых линейных ком­ бинаций столбцов матрицы

Р k s

(1.6)

Да(У)= П П (У—

5=1 } = \

S + 0

ранг которой равен kc, и, в частности, из kQлинейно неза­ висимых столбцов этой матрицы. Зная матрицу /С, легко определить М и А, используя соотношения (1.4) и (1.5).

Пусть /Со, Ла, Мо (<т = 1,

р) — построенные таким

путем матрицы. Тогда выражения

КоМа, М^1Ма, М71Л0Ма,

где Me — произвольная невырожденная матрица порядка ka, представляют общий вид матриц, осуществляющих раз­ ложение (1.2). Соответственно, если N — квазидиагональная матрица, составленная из матриц Mat то КМ является

общим выражением матриц, преобразующих

матрицу U

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

 

Для дальнейшего представляют интерес те матрицы Ко,

которые нужное число раз дифференцируемы на

промежут­

ке [О, L1. Вопрос о дифференцируемости матрицы Ко тесно связан с дифференцируемостью матрицы Ла (U). Покажем,

что матрица

Д* (U)

дифференцируема на [О, L] по

край­

ней мере

столько же

раз, сколько

раз дифференцируема

матрица

U (т).

 

 

 

 

Рассмотрим многочлен с коэффициентами, зависящими

от параметра т:

 

 

 

 

Д (А) = А"

аг(т) А.”-1 -f

+

ап-\ (т) A -f ап(т).

(1.7)

Допустим, что корни этого многочлена, не обязательно все разные, разбиты каким-нибудь образом на две группы:

..., А*1) и AI(2), ..., А^. В соответствии с этим много­ член Д (А) можно представить в виде

Д (А) •—Дх (А) Д2 (А),

(1.8)

где A i(A )— многочлен степени k с корнями А!^, ..., А*\

а Д2 (А) — многочлен степени

тс корнями

Ар, ...,

AJ?.

Пусть

 

 

 

 

 

Дх М =

Afe + а л.(т) Aft_1 +

-f a ft_, (т) A +

a ft(т),

|

^

Дг(^) =

^т Н" Pi (т) AW 1-f-

~h Pm-l (Т) A +

Pm (T).

,

 

Имеет место следующая

Ле м м а 1.1. Если на промежутке [О, L]

1)коэффициенты aj (/ = 1, ..., п) многочлена Д (А) яв­

ляются I раз дифференцируемыми функциями от t (/ > 1) и

2)

I х[1)(т) - xf>(т) | > О

( 1. 10)

1 , . . . , J — 1

то коэффициенты а/ (/ = 1, .... k) и Р/ (/ = 1, ...,

т) мно­

гочленов Д2 (X) и Л2 (^) являются по крайней мере I раз диф­

ференцируемыми на [О, L] функциями от т.

много­

Для д о к а з а т е л ь с т в а леммы подставим

члены (1.9) в равенство (1.8) и приравняем коэффициенты при одинаковых степенях X. Получим

Фх =

«1 +

Pi — Gj =

О,

 

Фг =

0С2 -j- OCJPJ -J- Р2

#2 — 0,

 

Фз =

аз +

а 2Р1 + а 1Р2 + Рз — о3 = 0,

( 1. 11)

фи—1=

«fcPm—1-f- &k—lPm— #n—1— o,

 

Фл =

«ftPm

an = 0.

 

Функциональный определитель системы (1.11) имеет вид

 

 

1

0

D (q>lt .. • »Фл)

Pi

1

Р2

Pi

D (alf ..

»Pm)

 

*

 

 

0

0

 

 

0

0

0

1

0

0

«1

1

о

«2

«1

 

 

 

Рет-1 0 0 Р,„ 0 0

0 0

0 0

0 0

«А

0

( 1. 12)

Определитель в правой части равенства (1.12) в точности равен результанту R (Аь Д2) многочленов (1.9), который в силу условия 2) леммы на [0, L] нигде не обращается в нуль. Следовательно,

D(фг, . . . , фд)

, л

(0 < т < L).

(1.13)

D (av

pffI) ^

U

 

 

Пусть т0 — произвольная точка отрезка [0, L]. В силу (1.13) в некоторой окрестности точки (а±(т0), ..., ап (т0)) существуют непрерывные функции

«1 = flfav • • ** ^л)»

«А

f k (flj, • • . ,

(1.14)

 

 

Pi = fk-\-1(^i> • •« » ^л)»

Pm = fn ( ^ 1» • • • * ®л)»

удовлетворяющие системе (1.11) и условиям

а 1 (То) = /l (^1 (То)> * ап(уо))*

P m W = / « ( а 1 <то)........... ап (to)).

причем эти функции имеют непрерывные частные произ­

водные по alt а2г

ап в

некоторой окрестности точки

(fli (TQ), ... | Qn (TQ)).

= 1, ..., n) — дифференцируемые функ­

Поскольку CLJ (j

ции от т, то из

равенств

(1.14) следует существование

производных по т от а 1? ...,

Рт в окрестности точки т0. Учи­

тывая это, продифференцируем равенства (1.11) по т. По­

лучим

 

 

 

 

 

dat

 

dPi

 

dax

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

dx

 

dx

 

 

 

 

da,

,

dp,

Pi

dax

+

dP„

 

da2

 

 

 

 

~аГ + “ i

dx

dx

dx

 

dx

 

 

 

 

 

 

 

~лГ~ Pm +

a A dPm

 

dan

 

 

 

 

 

 

 

dx

 

 

 

dx

 

dx

 

 

 

Решая эту систему, будем иметь

 

 

 

 

 

 

dax

 

i

t l

(®1.

• • •

» Pm>

dar

’ ** * 5

dan \

}

dx

R( Al t A2)

dx

dx

) '

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

dat

 

 

*

 

dpm

_____________J

+Л (“ и

• ••

» Pm»

 

 

dan \ .

J

 

 

 

dx

* ‘

»

dx

 

 

 

 

 

 

 

 

 

 

 

 

(1.15)

здесь ф/ — функции,

дифференцируемые по своим аргумен­

там любое число раз.

 

 

 

 

 

 

 

 

 

 

Если существуют

вторые

производные

 

 

 

 

то из

(1.15)

можно

непосредственно

получить

выражения

для

 

 

 

Для этого достаточно

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

цировать равенства (1.15) пот и в правых

частях

исключить

 

 

с помощью тех же равенств (1.15). Этим

путем можно

получить все производные

по

т от a l}

|3т

 

dla y

 

dl$

 

 

 

в произвольной точке

ВПЛОТЬ ДО

 

 

■>т . Итак,

dx1 dx1

т0 £ [О, L] коэффициенты многочленов ДАи Д2 дифференци­ руемы по крайней мере столько раз, сколько раз дифферен­ цируемы коэффициенты многочлена Д. Лемма доказана.

Пусть теперь Д — характеристический многочлен матрицы U (т). Коэффициенты этого многочлена дифференци­ руемы по х по крайней мере столько же раз, сколько раз дифференцируема матрица U (т). Коэффициенты матрично­ го многочлена (1.6) равны коэффициентам многочлена Да (Я), корнями которого являются собственные значения матрицы 0, не включенные в группу а. Многочлен Да (Я) является множителем многочлена Д (Я), причем имеют ме­

сто неравенства (1.1). Условия леммы, таким

образом, вы­

полняются. Согласно лемме

коэффициенты

многочлена

ДС(Я) дифференцируемы по

крайней мере

столько раз,

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

Из вышеизложенного ясно, что если на рассматриваемом промежутке [О, L) линейная независимость каких-нибудь фиксированных kaстолбцов матрицы Дс (U) не нарушает­ ся, то в качестве искомой матрицы Ко (дифференцируемой столько же раз, сколько раз дифференцируема матрица U) можно взять матрицу, составленную из этих ka столбцов. Если же таких кд столбцов не окажется, то матрицу Кд следует набрать из ka линейных комбинаций столбцов матрицы Дg(U) соответствующее число раз с дифференци­ руемыми коэффициентами, выбор которых должен быть подчинен условию линейной независимости этих ka линей­ ных комбинаций.

§ 2. Построение формального процесса для расщепления системы дифференциальных уравнений

на независимые подсистемы меньшего порядка

Полагая, что

элементы матриц A (/) и В (t) в уравне­

нии

(0.1) как функции от t являются медленно меняющими­

ся

функциями,

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

х =

е/ и рассмотрим уравнение несколько более общего

вида:

А (т, е)

= В (т, е) х + / (t, т, е)

(т = е*). (2.1)

Здесь А (т, е), В (т, е) — квадратные матрицы порядка п, допускающие разложение по степеням параметра е:

 

А (т, е) = f

еЧ

W,

В (т, 8) =

f

 

(т)

(2.2)

 

 

 

ft= 0

 

 

 

 

 

fe=0

 

 

 

 

 

 

 

(Т € №, L1).

 

 

 

 

 

При

ЛА(т) =з Bk(т) =

0 (k = 1, 2,...)

и

е = 1

уравне­

ние (2.1) принимает вид (0.1).

на

сегменте

 

 

т •< L

Т е о р е м а 2.1.

Пусть

0

•<

a) f (t,

т,

е) — непрерывная

вектор-функция,

регулярная

относительно е в окрестности точки

е — 0,

б) матрицы

Ak (т),

Bk (т)

(k = 0, 1 ,2 ,...)

имеют производные по т всех

порядков, a AQ(т),

кроме того,

является

невырожденной

матрицей.

Тогда.,

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

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

матрицы V (т) = Ао1(т) В0(т) разбиты на р групп Я[а>, . . .

. . . ,

(а =

1 , . . . ,

р; У ka = п\при условии,

что

 

°

\

 

 

а==1

/

 

 

 

 

 

 

 

 

| Х,<сг) (т) -

Я}4 (т) | >

0

 

 

 

(2.3)

(s Ф о;

 

i = 1,

,

ka; } =

1,

,

ks; т £

[0, L]),

формальное решение системы (2.1) можно представить так:

 

 

 

* = 2 ^ а ( т ,

е) уа,

(2.4)

 

 

 

сг=*1

 

 

где уа (о =

1,2,

...,/?) есть решения уравнений

 

= Аа (т, б) Уа -J- Ма(т, 8) R (т,

в) f (t, т, е) (2.5)

 

 

 

( а = 1,

р),

 

а Ко, Аа,

Ма,

R

матрицы

типа

соответственно

п X ka,

ka X kQ,

ka X n,

n X п,

представляемые

формальными рядами

Ко (Т, е) =

2

| W (*). л „ (т, е) =

2

е‘ лЕ,61 (т),

 

Л=0

£=0

 

Л , (т, е) =

f

Л <т- «) =

2

(т).

 

6=0

А=»0

 

Д о к а з а т е л ь с т в о . Подставим выражение (2.4) в (2.1) и используем при этом равенства (2.5):

А(х, Б) S

8)-j/g + Ко(т, е) Ла (т, е) уа +

+

Ка (т, е) Ма (“Г, е) R (т,

е) / (t, т, в)

 

 

= В (т, е) ^

Ка(т, е) f/0 + / (*, т, е).

 

 

а=1

В полученном равенстве приравняем нулю отдельно коэффициенты при уа (о = 1, 2, ..., р) и свободный член. Будем иметь

А(х, в)

dK0 (т, в)

 

В (т, в) Ко(х, 8) (2.7)

+ Ко (X, е) Ла(т, 8)1=

 

( 0 = 1 , . . .

, р),

 

Л (т, е) 2

Ко (Г, 8) Ма (т, в) Я (т, в) - Ея] f (t, х, в) = 0. (2.8)

Для доказательства теоремы

нужно

показать возмож­

ность построения членов рядов (2.6), удовлетворяющих ра­ венствам (2.7) и (2.8).

Займемся сначала равенствами (2.7). Подставим в эти равенства ряды (2.2) и (2.6) и отделим коэффициенты при

одинаковых степенях

в. Получим

А—1

wrrlA—V—1], .

k

k—v

2

л . м

dt - 1+

2 Л, (Т ) £ Л Р*Я (1)А 5Р(1) =

v=o

 

v=o

/—о

=

У В,(X)Kla~n (Т)

(о = 1,

, р\ k = 0 , 1 , 2 , ...).(2.9)

На

Равенства (2.9) умножим слева на Ао1 и,

учитывая, что

AolB0 =

U, представим их в виде

 

 

 

V (х) Ко”1(Т) =

/СУ1(х) Ag4 (т) +

у

/СУ- '1(т) ЛУ1 (х) +

 

 

 

 

 

т а

 

 

 

 

 

л

(т)

d K [ak~ vl (Т)

R /чл

*-[*-*] ли ,

 

 

i4v—1

----- --------------(т) Д а

 

(Т) +

 

V— 1

 

 

 

 

 

 

+ у

A, (t) /С У "^ '1(Т)луз (т)

 

 

 

'“ °

обозначив

(0 = 1,

 

р ; * = 0 , 1 , 2 ,

Наконец,

 

 

 

 

 

 

0 у -Ч (х) = у

/СУ“ /] (Т) лУ 1 (г) +

 

 

 

 

т а

 

 

<o<y-v>(т)

 

 

 

+ Ао1(т) ^

/4v_i

 

- Bv (Т) /С У '* 1 (т) +

( т )

.

 

ат

 

 

 

 

 

V=1

 

 

 

 

 

 

 

 

[fe-v-/] /яЛ А(/]

 

= 1,

р), (2.10)

+ *£ А (т)/С&

(т) ЛУ1(х) |

/to

 

 

 

 

 

 

 

 

будем иметь

 

 

 

 

 

 

 

U (х) Я У1 <т) =

/СУ1 W ЛУ1 (х),

 

 

 

 

и (х) /сУ1 (т) =

/СУ1(X)лУ 1 (X) +

/СУ1 № лУ 1 (х) + о у -11 (X)

 

 

 

 

(а =

1,

Pt

&

^»2f .. •). J

 

 

 

 

 

 

 

 

(2. 1-1)

Пусть /С = (/Ci (х)

/СР

(х)) — матрица,

преобразующая

при условии (2.3) квадратную матрицу U к квазидиагональ-

ному виду

 

 

 

 

 

 

 

 

 

ЛН

Л

‘*

 

о

 

 

 

 

 

о

 

 

 

 

 

 

 

4

 

 

Л р(тГ

 

 

и дифференцируемая столько же раз, сколько раз дифферен­

цируема матрица

£/.

----------------- --—

*

*) При k < тгнужно принять J s O .

А—г.

Первое из равенств (2.11) будет выполняться тождест­ венно, если принять

К 1С ] ( х ) 3=К о (х),

 

Ag4 (х) =Л о (X)

(о —

1, р ).(2,.12)

В самом деле,

используя равенства (1.2) и (1.3), имеем

w d ?1 = и Ко =

2 K A M SKO =

КоАо =

л р д р .

 

 

 

S = 1

 

 

 

Пусть /Сст], Л[,0];

; Кок~1\ Ад*-1-1уже найдены. При

этом Daft_11 будет

известной матрицей. Определим

и

№ .

1)-е равенство (2.11) слева на М. Полу­

Умножим (k +

чим

 

 

 

 

 

Л ф ' =

(ф ]Ао + М К о №

+ M DLck~ l\

(2.13)

Здесь

 

 

 

 

 

 

 

QP] - A fK ? ].

 

 

(2.14)

Матрицу Qo] , состоящую из п строк и kg столбцов, коагулируем (разобьем на блоки) следующим образом:

где

<?й] = м д 1.4

— субматрица типа ks х ka.

Так как Л имеет квазидиагональную структуру, то ра­ венство (2.13) распадается на р независимых матричных равенств

л 5< $ ]

- Qg]A0 +

MsKoЛ»1 +

МД * -Ч

(s= 1...........р).

Отсюда, принимая во внимание равенства (1.3),

полу­

чаем

 

 

 

 

 

 

ЛоСЙ1 - <3$л +

Ай1 +

 

(2.15)

 

ASQ $

= «ЗЙ’Л,, +

МД " -11

(s Ф а).

(2.16)

Из

(2.15) непосредственно

находим

 

 

 

-

A.QW _ ф

Ао _

 

(2.17)

Здесь Qo*3 — произвольная, достаточное число раз диффе­ ренцируемая матрица порядка £0. В частности, можно

принять QoS = 0.

 

 

 

 

 

 

 

Остальные субматрицы матрицы

однозначно опре­

деляются

алгебраическими

соотношениями (2.16),

ибо

в

этих матричных равенствах Ла и Ла (s Ф а) не

имеют

об­

щих собственных значений (согласно лемме 7.7.1).

 

 

Определив из

(2.16)

 

(s ф о)

и задавшись

матри­

цей Qaa,

будем

иметь Qak\

после

чего легко

построить

матрицу К1ок] по формуле

 

 

 

 

 

 

 

 

=

К ф \

 

 

(2.18)

которая получается из равенства (2.14) умножением обеих его частей слева на К и использованием соотношения (1.5).

Таким образом, с помощью рекуррентных соотношений (2.16), (2.17) и (2.18) можно последовательно определить

члены формальных рядов Л£Р, /С0 3; Л0], Лс?1;

Остается определить еще матрицы Ма и R, фигурирую­ щие в равенстве (2.8), которое можно записать в виде

(т, е) К (т, е) М (т, е) R (т, е) — Е„] f (t, т, е) = 0, (2.19)

где

 

/М Л

 

К =

Кр\ М = \

.

 

'Л Г /

 

Поскольку f {t, т, е) — ограниченная функция, то для выполнения равенства (2.19) достаточно, чтобы имели мес­ то следующие соотношения:

 

К (*, е) М (t, е) =

Ет

 

(2.20)

Обозначим

Л(Т, е)Я (т, е) = Еп.

 

(2.21)

 

 

(м \к\

 

 

 

 

 

 

 

а:'*1 =

(к!*1

м™- 1

 

 

 

 

 

1

V

 

Тогда

 

 

\м 1кУ

 

 

 

 

 

 

К =

V е‘к “ ' W ,

УЙ=

2 e*Mw (t).

(2.22)

 

ь&>

 

*-о

 

 

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