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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

Г л а в а VIII

АСИМПТОТИЧЕСКОЕ РАСЩЕПЛЕНИЕ НЕСТАЦИОНАРНОЙ СИСТЕМЫ УРАВНЕНИЙ

Решение и исследование нестационарной системы диф­

ференциальных уравнений

 

A ( t ) ^ - = B(t)x + f(t)

(0.1)

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

В настоящей и следующей главах приводятся два раз­ личных метода асимптотического расщепления и интегри­ рования линейной дифференциальной системы, содержа­ щей параметр е и совпадающей, в частности, при е = 1 с системой (0.1) *). Применимость этих методов для при­ ближенного расщепления и интегрирования дифференциаль­ ных систем вида (0.1) можно ограничить классом довольно распространенных в приложениях систем с медленно меняю­ щимися коэффициентами.

*) Асимптотическому расщеплению линейных дифференциальных систем, содержащих параметр, на независимые подсистемы уравнений пер­ вого порядка и асимптотическому интегрированию таких систем посвя­ щено большое количество работ. Некоторые из них указаны в нашей краткой библиографии. Значительно более подробный перечень работ в этом направлении содержится в библиографии, приведенной в работах [7, 8, 49, 52].

$ 1]

Д И Ф Ф Е Р Е Н Ц И Р У Е М О С Т Ь М А Т Р И Ц Ы

181

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

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

рой

разбиты на р групп А.1а),

. . . , ^

( 0 = 1

, ............2

р;

2

ka =п ] так, что на промежутке 0 -<т ■< L выполняют­

с я

 

J

 

 

 

 

 

 

 

\

ся

условия

 

|Ä,}*4 -

Я,)4 |> О

 

(1.1)

 

 

 

 

 

 

 

(o^=s; г

=

1,2..........ka\

/ = 1 , 2 , . . . .

kJ,

 

при

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

можно

представить

так

(см. гл. V):

 

 

 

 

 

 

 

 

 

 

 

 

U = 2

КоЛаМо,

 

 

(1.2)

 

 

 

 

а=І

 

 

 

 

 

где

Ко,

Ко, М а — матрицы

типа соответственно п X ka,

ko X ka,

ko X n,

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

 

 

 

 

 

MoKs =

 

 

s =

ст,

 

(1.3)

 

 

 

 

0,

 

s =

а

 

 

 

 

 

 

 

 

 

(Er — единичная

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

 

 

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

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

 

 

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

 

 

ол„

 

 

 

 

 

 

‘Л,

 

т г\

 

 

 

 

 

 

 

K = ( K 1 . . . K J ,

 

Л =

О

м =

 

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

м ..

 

 

 

 

и = КАМ,

Еп.

 

 

(1.4)

 

 

 

 

М К =

КМ =

 

(1.5)

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

бинаций столбцов матрицы

 

До (£/)= П гі { U - % fE n),

( 1.6)

s= \ /=1 БфО

182 А С И М П Т О Т И Ч Е С К О Е Р А С Щ Е П Л Е Н И Е С И С Т Е М Ы [ГЛ. V III

ранг которой равен ka, и, в частности, из kQлинейно неза­ висимых столбцов этой матрицы. Зная матрицу К, легко

определить М и Л, используя соотношения (1.4) и (1.5).

 

Пусть Ко, Ад, Мо (ст = 1, ..., р) — построенные таким

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

N 7xN0Na,

где

Na — произвольная невырожденная

матрица порядка

ko,

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

ложение (1.2). Соответственно, если N — квазидиагональ­

ная матрица, составленная из матриц Na, то KN является

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

матрицу

U

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

 

Ко,

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

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

промежут­

ке [О, L). Вопрос о дифференцируемости матрицы Ко тесно

связан с дифференцируемостью матрицы Да (U). Покажем, что матрица Д0 (U) дифференцируема на [О, L] по край­ ней мере столько же раз, сколько раз дифференцируема матрица U (т).

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

Д ( Я ) ^ Г + аг ( т ) Г - ‘ + . . . + ал_,(т)Я, + ал(т). (1-7)

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

ХІ'\ ..., Xi'1 и

Хр1, ..., Хщ*.

В

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

этим

много­

член А (X) можно представить в виде

 

 

 

 

 

 

Д (X) =

Дх(Х) Д2 (X),

 

 

 

(1.8)

где Aj_ (X) — многочлен степени k

с

корнями

Х[1),

...,

Х*,1’,

а Дг (X) — многочлен

степени

т с

корнями XS2),

...,

Хт.

Пусть

 

 

 

 

+ ak—\ (т) X 4- ak(т),

 

 

 

Ді(Х) = Х

-f-a1(r)Xft

+

 

 

1

(19)

Аг (X) — X

+

ßj (т) Xm

1+ •••

-)- ßm_i (т) X + ßm (т).

I

 

Имеет

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

 

 

 

 

 

 

 

 

Л е м м а

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

 

 

(X) яв­

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

а/ (j =

1 / ... ,

п)

многочлена

Д

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

1) и

2)

 

|Я/(1)(х) — Я,)8>( т ) | > 0

 

 

(1.10)

 

 

(/ = 1,

 

 

1,

. . . . т),

 

 

 

 

§ П

Д И Ф Ф Е Р Е Н Ц И Р У Е М О С Т Ь М А Т Р И Ц Ы

183

то коэффициенты a,j (j — 1,

k) и ß/ (/ = 1,

m) мно­

гочленов

(К) и A2 (1) являются по крайней мере I раз диф­

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

много­

Для

д о к а з а т е л ь с т в а

леммы подставим

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

Фі =аі +ßi — аі =о,

ф2 =

а 2 +

a ißi +

ß2 — а г =

0>

 

 

Фз =

аз +

«2ßi +

«iß2 +

ß3 — °з =

 

( 1. 11)

фа—1=

ttfcßm

-f- Ctfe—lß,,i — Cln—I — О,

 

 

Фа =

“*ßm

~1

«а =

°-

 

 

 

 

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

 

1

0 . .

0

1

0 . . .

0

0

Р (ф і...........фп )

ßl

1 . .

0

«i

1 . . .

0

0

ß2

ßl

.

0

a 2

“ i . . .

0

0

D (CÜJ, . . . . ßm)

 

 

 

 

 

 

 

 

 

0

0

.

ßm—I

0

0 . . .

a k

OCft-l

 

0

0

.

ßm

0

0 . ..

0

a k

 

 

 

 

 

 

 

 

( 1. 12)

Определитель в правой части равенства (1.12)

в точности

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

D (фі., - ■■■Фа)

/ п

( 0 < т < L).

(1.13)

D ( а , ...........ßm) ^

U

 

 

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

ai ~

fi (ßi>

• • ■

Oa).

а* =

fk (a-v

■■■

( 1 .1 4 )

ßi =

/*+1 (ai> • •

. a«)

ßm f n i f l x ' ■• ■I ^а)>

184 А С И М П Т О Т И Ч Е С К О Е Р А С Щ Е П Л Е Н И Е С И С Т Е М Ы [ГЛ. V III

удовлетворяющие системе (1.11) и условиям а і Ы = h К (т0)...........а„(т0)),

р„,(т0) = /ЛМто). ■• ■. ал(т0)),

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

произ­

водные по аъ а2, ...,

ап

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

точки

(at (т0),

а„ (т0)).

1,

п) — дифференцируемые функ­

Поскольку а/ (/ =

ции от т, то из равенств (1.14) следует существование производных по т от а,, ..., ßm в окрестности точки т0. Учи­ тывая это, продифференцируем равенства (1.11) по т. По­

лучим

 

 

 

da.

 

4ßi _

da.

 

 

 

 

I

da2

 

 

 

dx

1

dx

dx

P a

dx

-1- В

da.

I

rfß«

da2

dx

г

1

dx

 

dx

dx

 

 

 

dak

ßm +

ak

dßmdx

dan

 

 

 

dx

dx

Решая эту

систему,

будем иметь

 

 

da,

 

1

 

 

dx

' R ( Д „ Д 2)

• ,

ß m ;

dßm

 

I

 

 

dx

R

' ' h , ( 0 4 ,

. .

ß m ;

( А , , Д 2)

 

 

da,

 

 

dan

dx

*

dx

da.

 

 

dan

dx

r

*

dx

(1.15)

здесь ф/ — функции, дифференцируемые по своим аргумен­ там любое число раз.

 

i-ч

существуют вторые

производные

(і2ал

d2 ап

 

Ьсли

 

» ■• • »

то из (1.15) можно непосредственно

получить

выражения

для

d2a,

 

 

d~ßm

TT

 

 

 

,

 

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

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

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

исключить

с/ах

dß.

с помощью

тех же равенств (1.15). Этим

dx

 

 

 

путем можно

получить

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

от а,, ..., ßm

вплоть д

d1а ,

 

d!ßm

, 7

 

 

 

о /

 

..........- , ,

. Итак, в произвольной точке

 

 

dr

 

 

dx

 

 

 

 

§ 2] П О С Т Р О Е Н И Е Ф О Р М А Л Ь Н О Г О П Р О Ц Е С С А 185

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

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

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

образом, вы­

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

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

многочлена

Аа(к) дифференцируемы по

крайней мере

столько раз,

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

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

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

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

Полагая, что элементы матриц A (t) и В (/) в уравне­ нии (0.1) как функции от t являются медленно меняющими­ ся функциями, введем так называемое медленное время т = ві и рассмотрим уравнение несколько более общего

186

А С И М П Т О Т И Ч Е С К О Е

Р А С Щ Е П Л Е Н И Е

С И С Т Е М Ы

ІГЛ . V III

вида:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А(%, е ) - ^ - = В(т, г)х + f(t, г, е)

 

= et).

(2.1)

Здесь А (т,

е),

В (т,

е) — квадратные матрицы порядка л,

допускающие разложение по степеням параметра е:

 

 

А (Г, 8) =

2

е Ч

(*),

В (х, е) =

I

е Ч

(х)

(2.2)

 

 

 

 

k=0

 

 

 

 

 

 

А=0

 

 

 

 

 

 

 

 

 

 

 

(г б [О, Ц).

 

 

 

 

 

 

 

При Ak (x) =

Bk (т) s

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

и

е =

1

уравне­

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

на

сегменте

0

с

т

■< L

Т е о р е м а

2.1.

 

Пусть

а) f (t,

г,

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

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

 

регулярная

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

г — 0,

б) матрицы

Ak (х),

Bk (т)

(k =

0,

1,2,...)

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

порядков, а

Л0 (х),

кроме того,

является

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

матрицей.

Тогда,

 

предполагая, что собственные

значения

матрицы U (х) =

ЛсГ1(х) В0(х) разбиты на р групп ?ча), . . .

. . . , XjfJ |ог =

1......... p - ^ k a = n j при условии,

что

 

 

 

 

 

 

 

\Kila)( x ) - X f )(T)\>0

 

 

 

 

 

(2.3)

(s ф а;

 

і = 1 ,

. . . , А : а; / = 1, . ..

,

ks; х g

[0,

Ц),

 

формальное решение системы (2.1)

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

 

 

 

 

 

 

X =

2

 

г) уо,

 

 

 

 

 

 

(2.4) .

 

 

 

 

 

 

 

а=і

 

 

 

 

 

 

 

 

 

где уо (а =

1,2,

...,

 

р) есть решения уравнений

 

 

 

 

L

 

= Ла (X, е) У о

+ Mo ( X ,

е) R (х, е) / (t,

 

х, е)

(2.5)

 

 

 

 

 

 

 

(ст= 1,

.. . ,

р),

 

 

 

 

 

 

 

а Ко,

Ла,

Mo,

R

матрицы

типа

 

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

п X ko,

ka X ka,

 

ko X n,

n

X n,

 

 

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

5 2]

П О С Т Р О Е Н И Е Ф О Р М А Л Ь Н О Г О

П Р О Ц Е С С А

187

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

 

 

 

 

 

 

Ко (т, в) =

2 гкК1ок] (т),

л„(т, е) =

2

е‘Л ?] (т),

 

 

 

 

fc=0

 

 

*■—"

(2.6)

 

 

 

 

 

 

 

УЙСТ(т, е) =

2

(т),

/? (т, е) =

J

М-

 

 

 

 

А=0

 

 

А=0

 

 

Д о к а з а т е л ь с т в о . Подставим

выражение

(2.4)

в

(2.1) и используем при этом равенства (2.5):

 

 

 

р

d K G (т,

е)

 

 

 

 

 

 

 

 

 

 

 

 

Л ( т ,

б ) 2

dx

Уа Ко (т, е) Ла (т, е) Уо +

 

 

 

СГ= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Ко (t, е) Mo (Т, е) R (х, е) / (t, т, е)

 

 

 

 

 

=

В (х, е) ^ Ко (г, е) уа+ f (t, г,

е).

 

 

 

 

а=1

 

 

 

 

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

А (т,

dK„ (х, е)

~

 

В(Х, е) Коіу, е ) (2.7)

е) -------- ^

-------------- Ь К о ( т , е ) А о

( т , е )

 

 

 

( 0 = 1 ,

р),

 

л (т,

е) 2

К

е) Кіо (т, е) /? (т, е) —

f(t, т, е) = 0. (2.8)

 

ст=І

 

 

 

 

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

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

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

е . Получим

А — 1

 

 

j f t ' t * — ѵ — Ч /

*

А — ѵ

£

А

„т - +

£

А, (X)2 Й* - ѵ-д (X)Л И (X) =

ѵ=0

 

 

ѵ=0

/= О

=

p

/

( T ) i ( ' M ] (T)

( 0

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

 

/ - 9

 

 

 

 

188

 

А С И М П Т О Т И Ч Е С К О Е

Р А С Щ Е П Л Е Н И Е С И С Т Е М Ы

[ГЛ V III

Равенства (2.9) умножим слева на Ло1 и, учитывая, что

AöxB0 =

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

 

 

 

 

и (Т) K la ](Т) =

Kik] (т) Л™ (т) +

у

А

- ' 1 (г) Ад/] (г) +

+

Ао1(т) ^

Аѵ- 1(t)

ак[к- ѵЦх)

■ßv (T) K[cA' Vl(T) +

 

âx

 

 

k—V

Ѵ = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

V

Лѵ(т)А А ѵ~ /](тА ' А )

 

 

 

 

Наконец,

обозначив

(0 =

1, ...J ,

p\

£ = 0 , 1 , 2 ,

 

 

 

 

 

 

 

 

 

Й Ы 1(т) =

^ Ѵ а М

( т ) 4 Л (T) +

 

 

 

 

 

 

 

/=І

 

 

 

 

 

 

 

 

 

+

Л^1(т)

*

Лѵ-І (Т) -

 

/r)

-

в ѵ (т) /йА_ѵ] (Т) +

V

- dx

 

 

k — V

 

ѵ = I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

V Лѵ (т) К1а Ѵ

(т) А^1 (т)

 

( 0 = 1 .......... р), (2.10)

 

/= О

 

 

 

 

 

 

 

 

 

 

 

будем иметь

 

 

 

 

 

 

 

 

 

 

U (т)

А ° ] (т) =

А ° ] (т) А о 1 (X).

 

 

 

 

 

У (т) К1о ](т) =

(т )А ^

(т) + /^0] (т)

(т) + А *“ 11 (г)

 

 

 

 

 

 

 

( 0 = 1 .......... р\

k = 1,2,

.. .).

 

 

 

 

 

 

 

 

 

 

 

 

(2.11)

Пусть А =

(/<і (т) ... Ар

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

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

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

ному виду

 

 

А

(т)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

А(Ѵ =

І

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

л

р (т)'

 

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

*) При А< Т у нужно принять У, g 0.

S 2]

П О С Т Р О Е Н И Е Ф О Р М А Л Ь Н О Г О П Р О Ц Е С С А

189

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

^Сс] (т) = Ко (т),

А^0] (т) = А0 (т)

(а = 1.........р). (2.12)

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

UKla] = и Ко = 2

К А М Д а = КаАа =

 

 

 

S=1

 

 

 

 

Пусть /Са°], Ла°]; . . . ; Как

Лп*

|] уже найдены. При

этом Da*-1-1 будет

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

и

Умножим

(k +

1)-е равенство

(2.11) слева на М. Полу­

чим

 

Qla ]Ac + M K cA ^ +

MDlak~l\

 

AQaft] =

(2.13)

Здесь

 

 

Q W = M K la \

 

(2.14)

 

 

 

 

Матрицу

, состоящую из n строк и ka столбцов,

коагулируем

(разобъем

на блоки)

следующим образом:

[ft]\

'Qxo

Ф ^

\ Q pl a J

где

=MsKok]

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

Так как А имеет квазидиагональную структуру, то ра­

венство (2.13) распадается на р независимых матричных равенств

Лф а1= Qlc]A„ +

+ м Л к~']

(S =

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

Отсюда, принимая

во внимание равенства

(1.3),

полу­

чаем

 

 

 

 

A(;Q!£] =

Q&]A + Л ^1+

 

 

(2.15)

ASQSC] = (ф А 0 + М Л ~ 1]

(s Фа).

(2.16)

Из (2.15) непосредственно находим

 

 

 

Л™ = AoQao — Qao]Aa —

 

 

(2.17)

Соседние файлы в папке книги из ГПНТБ