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

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

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

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

рые нужное число раз дифференцируемы на промежутке [О, L). Вопрос о дифференцируемости матрицы Ка тесно связан с диффе­

ренцируемостью матрицы Дa(U). Покажем, что матрица Aa(U)

дифференцируемая на [О, L] по крайней мере столько же раз, сколько раз дифференцируема матрица U(т).

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

Д(Х) s l " + а,(х)Х»-‘ + ... + <J„_1(T)X+ e„(t). (9-1-7)

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

Я^\ ..., Яр. В соответствии с этим многочлен Д(Я) можно предста­ вить в виде

Д(Я) - Д1(Я)Д2(Я),

(9.1.8)

где А!(Я) — многочлен степени k с корнями Я р ,..., Яр, а А2(Я) —

многочлен степени т с корнями Я р ,..., Яр. Пусть

А^Я) э Я Ч а ^ х )^ "1+ ... + аА_1(т)Я +а4(т),

^

Д2(Я) = Я» + р1(т)Я*-1+ ... +рт _1(т)Я + Рт(^

 

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

 

 

 

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

являют­

1) коэффициенты ау (у = 1, ..., п) многочлена А(Я)

ся I раз дифференцируемыми функциями от %(12* 1) и

 

2) |Я р (т )-Я р (т )|

> 0

(/= 1 ,

у = 1,..., т),

(9.1.10)

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

(у =

1,..., к)

и |Зу- (у = 1,..., т) многочле­

нов ДА(Я) и Д2(Я) являются по крайней мере I раз дифференци­ руемыми на [0, Ь\ функциями от т.

Для доказательства леммы подставим многочлены (9.1.9) в ра­ венство (9.1.8) и приравняем коэффициенты при одинаковых сте­ пенях Я. Получим

<Pi =

° i "Ь Pi — ai

О»

 

=

^2 ”Н

■+* Р2

a>i —0,

 

Ф3 а3 + а2^1 + а102 + Рз ~ а3 —

(9.1.11)

Фп —1—

-1

- l^m ^n-1

 

fn S aAn -

e« = °*

 

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

(9.1.11) имеет вид

 

1

0

0

1

0

0

0

 

ft

1

0

a,

1

0

0

D <Ф1..... уя>

f t

f t

0

«2

a i

0

0

° S ...... Р,„>

 

 

 

 

 

 

(9.1.12)

0

0

 

0

0

 

 

 

ftn —i

a k

« * -

 

0

0

0

0

0

 

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

 

^

..... ч>„> Ф

О

( 0

^ т

L ) .

 

(9.! .13)

 

а1.....Рт^

 

 

 

 

 

 

Пусть т0 — произвольная точка

отрезка

[О, L],

В силу

(9.1.13)

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

(а,(т0) , д „ ( т 0))

существуют не­

прерывные функции

 

 

 

 

 

 

 

а 1=

/ j(^P •••> ац)> •••’

ак

f

■>•"> a>i)>

(9 |

ft

+

**• ’ “и)'

 

ftn =

 

 

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

и условиям

 

 

a,(t0) = A (ai(To)...... •••• ft„(To) = /„ O iK ) ’ •••» а,.(хо))>

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

а {, a v .... ап в некоторой окрестности точки (ц,(т0), ..., ап( т0)).

Поскольку о.j

(у = 1, — , п) — дифференцируемые функции от

т, то из равенств

(9.1.14) следует существование производных по

т от dj, ..., {Зт в окрестности точки т0. Учитывая это, продифферен­

цируем равенства (9.1.11) по т. Получим

rfcij

d pj

da j

i a 2 .

,

о

d a l j _ dh

da’2

dx

dx

dx

7 7 + а Ч 7 + ^ ~ 7 7 + 7 7 =

dx 1

 

 

 

dak n ,

 

_

dan

 

 

 

 

 

d т

^ dx

 

dx

 

 

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

 

 

 

 

 

d cij

_

i

 

 

da 1

d <

 

 

7 7

ЖД,. A2) 4>J

f a p ’ - * . f t „ ;

7 7 ’

dx

 

 

 

_

1

\

 

dax

<K

(9.1.15)

 

 

 

 

 

dr

f®P *•* *

 

 

 

«СДр Д2) ^11

ftn’ dx ’ • "

dx

 

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

бое число раз.

§ 9.2. Построение формального процесса для расщепления векторно-матричного уравнения

Полагая, что элементы матриц A(t) и B(t) в уравнении (9.0.1) как функции от t являются медленно меняющимися функциями, введем так называемое медленное время х = zt и рассмотрим урав­ нение несколько более общего вида:

Л(т, е ) ^ = £(т, е ) х + /(* , т, е ) (т =£*)•

(9-2.1)

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

со

со

 

 

 

А{х, с) = 2

B(s <0 = 2

(* s

1°. £])•

(9.2.2)

Jfc= 0

А = 0

 

 

 

При Ак(х) = Вк(х) = 0

(/г = 1 ,2 ,...)

и е= 1

уравнение

(9.2.1)

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

 

 

 

 

Т ео р ем а 9.2.1. Пусть на сегменте 0 < т ^ L

а) f ( t , т, е) непрерывная вектор-функция, регулярная от но­ сительно е в окрестности точки е = 0;

б) матрицы Ak(x), Вк(г) (А: = 0, 1, 2, ...) имеют производные по х всех порядков, а ,/40(т), кроме того, является невырожденной матрицей.

Тогда, предполагая, что собственные значения матрицы

т

= А -'(х)В 0(х)

 

 

разбит ы на р групп

 

(ог=

1,

р;

р

к0 — п) при уело-

 

£

6uuf что

 

о

 

 

0 = 1

 

 

 

 

 

 

 

| Х(")(х) -

Х<р(х) | г

с > 0

 

(9 2 3)

(s 9й a, i 1,

koi

j 1,

 

ks,

х £

[0, /<]),

Ф ормальное решение системы

(9.2.1)

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

(9.2.4)

0=1

где уа (а = 1,2, ..., р) есть решение уравнений

^ 7 = А0(т, е)уа + Ма(х, е ) Л ( т , z)f(t, г, е) (сг= 1, ..., р), (9.2.5)

а Ка, Аа, м а, Л — матрицы типа соответственно п х к0, ка X ка, к0Х п, и х п, представляемые формальными рядами

« во

е) = 2

(f), Аа(т, е) = 2 e*Aj*] (х),

к = 0

к=*0

 

(9.2.6)

00

во

М„(т, е) = 2

й(т, t) = £

А=0

* =0

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

М'* Е) = 2 [е *"7 ^' Е" Уо +

E)Aa(t, t)ya +

о = 1 L

+ *„(т, Е)М„(т, е)Л(х, е)/(1, т, е)] =

= В(*>E)S *о(Т>Е)Уа + /(<. «)• а=1

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

+ Ка(т, е)Л0(х, е)] =

^(х, ес(х, е)

 

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

 

(9.2.7)

 

 

А(х, е)2 Ка(г, е)А/а(х, е)Д(х, б) - £ „

/(Г, х, б)= 0 .

(9.2.8)

о- 1

 

 

Для доказательства теоремы можно показать возможность по­ строения членов рядов (9.2.6), удовлетворяющих равенствам (9.2.7) и (9.2.8). Займемся сначала равенствами (9.2.7), Подставим в эти равенства ряды (9.2.2) и (9.2.6) и отделим коэффициенты при одинаковых степенях б. Получим

к - 1 к к - у

2 л м

dx

/ = 0

 

v = 0

 

к

 

=

2 Ду.(т )4 * -'Ч 0

(сг = 1,.... р; /t = 0 ,1 .2 ,...)- (9.2.9)

 

У-о

 

Равенства (9.2.9) умножим слева на А0-1 и, учитывая, что

А ^ 1В0 = U, представим их в виде

U (r)K W (x ) =

КМ (т)АМ (т) + 2 А1*-Л(х)Л*Л(т) +

 

 

 

 

 

У-

1

 

 

 

+ V

W 2

[A - .W

 

------- +

 

 

J t - v

 

 

 

 

 

(сг= 1 ,...,р ; * = 0, 1, 2, ...)"■.

+ 2

А ,(т)Л « -’ ->| (т)Л'Л(х)

/"О

 

 

 

 

 

 

 

 

Наконец, обозначив

 

 

 

 

 

 

4 * '" ( * )

= 2

4 * " Л(*)Л1Л(*) +

 

 

 

 

;=i

 

 

 

 

 

 

 

 

Jk-vr

rfK!*_v,(t)

 

 

+

Л )"

*(х)2

 

Л -1 (*)

rft

- Bv(t)/C'‘ - y|(x) +

 

 

 

y=0L

 

 

 

 

 

Jfc -

V

 

 

 

 

 

 

 

+ 2

А ( х ) ^ _ ,_ л (х)Л“ (х)

( o = l .........P ) .

(9.2.10)

 

У-О

 

 

 

 

 

 

будем иметь

 

 

 

 

 

 

 

 

 

1/(т)Х«>'(т) =

4

0|(т)Л101(х),

 

1/(т)4«(т) =

4 “ (x )4 0l(t) + 4»'(х)Л ^(х) + 4 * '" ( х )

(9.2.11)

 

 

 

(аг = 1, ..., р;

к = 1 ,2 ,...).

 

Пусть К = (К 1(т)

 

УСр(т)) —

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

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

виду

 

 

 

 

 

 

 

 

 

 

 

 

/Л,(т)

0

\

 

Л(т) =

Ол »

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

к

* При к < г, нужно принять

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

/Ci01(x) гг Ка(т), Ai0,(x) = А0(х) (cr= 1 , р). (9.2Л2) В самом деле, используя равенства (9.1.2) и (9.1.3), имеем

 

V

 

 

t/K '0' = UKa = 2

KSASMSK. = К„А0 = х«ял*я.

 

5 = 1

Пусть

Ai°J,

К1ак~1], Aj*-11 уже найдены. При этом

jrjU-i] будет известной матрицей. Определим К.]® и Л ^1. Умно­ жим равенство (9.2.11) слева на М. Получим

AQj,*1= Q[ak] Аа + MKaA ak]l + MD[k~1].

(9.2.13)

Здесь

 

Q W = MK[k].

(9.2.14)

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

Q i«=col(Q |JI е ^ 1).

где Q\kJ = М '.К ^ — субматрица типа ks х ка. Так как Л имеет

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

A.QlS1= е ‘*'Л0 + М, КаА

( s = 1, .... л).

Отсюда, принимая во внимание равенства (9.1.3), получаем

A0Q'*I = <Ж'Л, + ЛШ + М0л1‘ -Ч,

(9.2.15)

AjQj*1=

+ MsDlJ- ~11

(s * a ) .

(9.2.16)

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

 

 

л'*1= к< $ $ -

e i‘'A0-

 

(9-2Л7>

Здесь Q[ak] — произвольная, достаточное число раз дифференциру­

емая матрица порядка ка. В частности, можно принять

= 0.

Остальные субматрицы матрицы Qlak] однозначно определяются

алгебраическими соотношениями (9.2.16), ибо в этих матричных равенствах A s и Aa (s ет) не имеют общих собственных значений

(согласно лемме 8.7.1).

Определив из (9.2.16) Q\® (s ^ о) и задавшись матрицей Q ^},

будем иметь Q ^1, после чего легко построить матрицу К\к] по фор­ муле

К\р = KQ[ak],

(9.2.18)

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

Таким образом с помощью рекуррентных соотношений (9.2.16), (9.2.17) и (9.2.18) можно последовательно определить члены фор­

мальных рядов Л*1], а£ ], ;Aj21, К ™ ,...

Остается определить еще матрицы М0 и Я, фигурирующие в ра­

венстве (9.2.8), которое можно записать в виде

 

[Л(т, е) К (т, е) М (т, е) Л(т, е) — Еп] f(t, т, е) = 0,

(9.2.19)

где

 

К = { К 1 Кр)у М = со1 (М, ... £р).

 

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

 

К{т, е)М(т, е) =

Еп,

(9.2.20)

 

А(т, е)Л(т, е) =

Еп.

(9.2.21)

Обозначим

 

 

 

 

к ш = ( к ш

Kl.k]),

— col (Mj*1

м '41)-

Тогда

 

 

 

 

Jt =0

 

 

Jt=0

(9.2.22)

 

 

 

Подставим разложения (9.2.22) в равенство (9.2.20) и прирав­ няем коэффициенты при одинаковых степенях г. Получим

к

К'°1мМ = Еп, 2 К ММ 1к~Я = 0 (* = 1 ,2 ,...) . I * 1

Отсюда, учитывая (9.2.12) и (9.1.5), находим

Лг

А/<01 = Л/,

= - М 2

( Л = 1, 2, ...).

 

Д 2(С/) = W -

\ ^ Е Ъ)( U - # > £ 3) =

о

г(< -£)2

(>-^)2

 

П о э то м у м о ж н о п р и н я ть 1 1

О

-2

Г

- Н ) :

 

 

 

 

 

 

 

 

 

 

 

 

* 1=

0

 

,

к 2=

' 01

 

 

 

 

 

 

 

 

0 -2И )1

 

-1

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х l2

 

 

 

 

 

 

 

 

 

 

 

-1

-1

к-

0 -

I

f - -7

I

1

JT-1.

 

1

 

 

 

/I

 

 

 

 

о

 

 

 

 

 

 

 

м = к

 

 

 

■ И) t —

 

О 2 к - - т | -1

 

f

t2

о

 

 

 

 

 

 

 

 

 

 

 

t

t

- f ,

О

 

 

 

 

 

 

 

 

 

A = M U K = о

 

t

о

 

 

 

 

 

 

 

 

 

0

 

0 - 7

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

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

1

* 1 -

' " 7

t - jti

- 1 -1

t t - j

 

, м 2= (0 2 1), л 1 =

^2—j2*

0

1 1

0 t

О п р ед ел и м т е п е р ь Л^Ч» ЛЦЧ» В д ан н о м сл у ч а е

 

 

О

О

 

 

 

 

О - м

. л р - ^ Ь - о .

 

 

2(1+?).

 

 

а

и Л /2Л ^ , к а к л егк о п р о вер и ть, — н у л ев ы е

м а тр и ц ы . У ч и ты в ая это , и з р а ­

в ен ств а

(9 .2 .1 6 )

п о л у ч и м

Q$—0 . П о л а га я

р ав н ы м н у лю и п р о и зв о л ь н ы е

м а т р и ц ы

и

б й * бУДем и м еть

б й “ °* То1Да

 

 

A14=A:QW-(A:1A:2)

=о.

Д а л е е ,

 

 

ей /

 

 

 

 

 

Лй----

О f* +2

, ЛЙ*-М 2^°1=*0.

 

Г* + 2

 

 

 

< « * - ! )

 

 

1(^-1)

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