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

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

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

где

Яe-f- ос|аЯ°

4~

4-

4~ <*ft0o.

фо (А)=

Положим

£lo == К.аО’Ч-

 

 

 

 

 

Здесь аа —

некоторый

й0-мерный

вектор (столбцовая

матрица типа ka X 1).

Тогда первое равенство (2.11) можно преобразовать к

виду

 

/Сафа (Ао) fl(j — 0.

(2.12)

Так как инвариантное подпространство /?а, соответствую­ щее группе а собственных значений матрицы U, цикличе­ ское, минимальный многочлен этого подпространства яв­ ляется многочленом степени ka, коэффициенты которого определяются по формулам Виета:

фа (^) =

4~ (х\а)Хка 1+

+ 0*0—iA + ос}£,

а ‘,<” = М<” +

+ 4 J .

а ^ - М 0)4 ° '+

+ * з 5 - ,« г .

< > = ( - 1 ) Ч Г

 

■4 ? •

Примем

 

 

а/о = сс}а).

Тогда

 

 

 

фа(Л.) в

Фо (Я),

 

 

и, значит, согласно лемме 1.1

 

 

 

фо(Лд) —фа(Ло) = 0.

В этих условиях

равенство (2.12) выполняется. Из осталь­

ных же равенств (2.11), зная

£1ст, последовательно можно

определить

| 2о,

Ь 0а'-

 

%2о=

иКайо + aioKoCtoКо (Ло -f aiaEka)аа>

|зa =

UКо (Ла 4" ttia^fca) ^а 4~ &2оКойа=

 

 

 

^

/Со (Ла 4" ®1оЛо 4“ ®2а^Ад)^а»

Если

в качестве а0 взять вектор, для

которого

ф0 (Я)

является

минимальным многочленом, то

векторы

|ю,

|* ао

будут линейно независимы. Тогда

преобразование

(2.8) , приводящее систему (2.1) к расщепленной системе

(2.9) , будет невырожденным. Теорема доказана.

2.3.Случай матрицы простой структуры. Если U име­

ет простую структуру, то равенства

(1.2)

остаются в силе

и в том случае, когда р = п. При

этом

К есть матрица,

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

сл :дующая

Те о р е м а 2.3. Если U матрица простой структу­

ры, то решение уравнения (2.1)

может быть представлено

в виде

 

 

п

d q a

х== 2 %\oQa>

—^

b a lo<7a = 0,

о=1

0 — собственное значение, а

где | i а — Ко, акт = —Я0

Ко — ему отвечающий собственный вектор матрицы U).

Эта теорема следует из теоремы 2.2, если (при условии, что U — матрица простой структуры) принять р = п.

§ 3. Преобразование однородной нестационарной системы дифференциальных уравнений к расщепленной системе

Рассмотрим однородную нестационарную систему вида

Л(-1 ,е ) - £ - = В (т ,г )*

(т = е*),

(3.1)

где А (т, е), В (т, е) — квадратные матрицы порядка п, допускающие разложение по степеням параметра 8 (сходя­ щееся или по крайней мере асимптотическое):

А (т, е) = 2 &% (т). В (т, е) =

2 вкВ„(т) (т € [0, £,]),

А=0

А=0

причем det Л„ (т) фО.

(3.2)

 

3.1.Преобразование к расщепленной системе уравнений

второго порядка. Предположим,

что

матрицы

А

(т,

е),

В (т, е) четного порядка

п и собственные значения матри­

цы U (т) = ЛГ 1 (т) В0 (т)

могут

быть

разбиты

на р

=

п/2

групп, по два собственных значения в каждой группе, так, что

|

(т) — Я/8) (т) | > 0

(£, / =

1,2; s Ф а; х £ [О, L)). (3.2а)

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

 

 

Т е о р е м а З Л . Если на сегменте [О, L] а) матрицы

Ak (т), Вк (т) (k = 0,

1 ,2 ,...)

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

порядков, б) инвариантные подпространства R±, R2.......Rp п-мерного пространства R, соответствующие указанным выше р группам собственных значений матрицы U, явля­ ются циклическими, то формальное решение уравнения ( 3 . 1 ) может быть представлено в виде

ИГ

(3.3)

где скалярные функции qa являются решениями уравнений

d?qn

 

-

 

do„

-

 

(3.4)

—£1 +

а ю (т>8)

----Ь а 2а fa, 8) Qo= 0

 

 

(ст — 11 2, .». ,

р),

 

£/а fa, е) (/ =

1,

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

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

 

 

 

 

 

Ъо (х. е) =

2 a*g? (т)

(1 = 1,2);

(3.5)

 

 

 

k—0

 

 

 

 

ajtj fa, е) (j =

1,2) —скалярные функции, представляемые

рядами

 

 

 

 

 

 

 

Що(х, в) =

2 Л

# 1 (х)

(1=1,2).

(3.6)

 

 

 

h=Q

 

 

 

 

Д о к а з а т е л ь с т в о .

Подставим выражение

(3.3) в

 

 

 

d2qa

помощью (3.4) и в

полу­

(3.1), исключая при этом

с

ченном таким образом тождестве

приравняем коэффициенты

при q0

dаа

(а — 1.......р). Будем

иметь

и -jj—

i4V8 ' ^ ---- Ьо«|11+ &о) = ВЁ|«

(<Х=1.............. ),

 

А

— ||о012оj-

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

или, компактнее,

 

 

 

^

£io°&/o +

S/+10) = ^ £/0

(3*7)

( / = 1 , 2 ;

а = 1 , 2 ,

, р \ | з о = 0 ) .

 

Подставим ряды (3.2), (3.5) и (3.6) в (3.7) и приравняем коэффициенты при одинаковых степенях е. Получим

 

Ж * =

v$ 2 +

«jffeW.

 

 

r if t .

-

 

+

«

+ « $ $ ?

“уот,/fl0)

 

s f t о = и й ? + «

+ «

- d'-у.

 

 

 

 

 

 

 

(3-8)

 

д а . = у д а + д а д а + в

 

где

( /= 1 ,2 ;

о г= 1 ,2 ,

р ; Ь о 1■ = 0),

 

 

 

 

 

 

 

д а - - v

Iе д а +

« д ад а - д а . и + i f - .

д а —

v

£ l£ /o

+ ^ 2 ^ /с

+ ^ 2 (« /а Й о

— |/+ ]1а) +

,

 

>1

[Д1]ьЮ] .

„[°ы и t£i]

 

I а /0 5ю т* a /ci Si® — 5

/+I0

и,

вообще,

 

 

да-11= — а°'

BiiSc-11 +

 

V.

« № > ] 4 -

 

 

^ ~ ) \

~ >а1х° + dT

+ в , д а +

+ а , ( « д а - д а j + л _ / д а д а +

-

д а . - §dx) I+

 

 

tft-2]

Первое из соотношений (3.8) в развернутом виде запи­

шется так:

 

 

Й Р -Ы Й Р + аВИР. )

(3.9)

о = щ \У+

« Ж 1.

 

Умножим первое равенство

(3.9) слева на U и сложим

со вторым. Получим

 

 

+ <4Ж ]= о,

ИЛИ

(U2+ а.\0,}и + <4о £„) l!”1 = 0.

Поэтому система (3.9) эквивалентна системе

(U * + «S№/+<4!?£J ^ ] - O,

® = ЦЦ + аВ Д .

Тем же путем преобразуем (k -f- 1)-е соотношение (3.8). Имеем

Е» -

У ®

+

« й ®

+ а Й ® -

4 J - 11,

(3.10)

 

0 = V &

+ « й ®

+ «Й1?!? -

d i r 1.

Отсюда

 

 

 

 

 

 

 

 

 

( V ® +

«

» +

»

)

+

 

 

ИЛИ

 

+ ай1! ®

+ аЙ ® — U d \r]— 4 ‘_1] = о.

 

 

 

 

 

 

 

 

(U*+ а\Уи + с4“1£„)5Й1=

_

(ай>1/ + „« £ „) ®

+

где

 

dlak~u =

 

4- 4*-1].

 

 

 

 

 

 

 

Таким

образом,

система

 

(3.10) эквивалентна

системе

(и г + <АУи + <$}Е„) ЕЙ1 = -

(«Й%/ + с$ Е „ ) ЕЙ1 + 4*~ч .

® - I ® + « й ® 1 + а Й ® -

Вместо уравнений (3.8) далее будем рассматривать эк­ вивалентные им соотношения

Фа (U)

=

О,

|f f

=

U |ff + a f f |f f ,

фо (О) iff =

- (ocffu + aff£„) |ff +

&к~'\

(3.11)

р Ш

__ п р Ш

J_ „МЛО]

, „[0]tLfej

A h - 1 ]

 

Ь2а

и Ь№

I

a la g ia

T r i a b l e

— ^l<x

 

 

 

(* =

1,2, ...),

 

 

Фа( 1 . ) ^ у + а т я + а ВД

Используя ( 1 . 2 ) (при соответствующем разбиении соб­ ственных значений матрицы U на непересекающиеся груп­

пы),

равенства ( 3 . 1 1 )

приводим

к виду

2

Ks4,„(As) M s|f f =

О,

 

 

Sff =

i C ] + < «

],

 

 

 

р

Ksфс (AJ Ms|f f

=

-

г

I (3.12)

2

(a ff{/ +

а$Е„) iff + < ft'\ '

S=1

 

 

 

 

 

 

iff =

Oil!1+ affsff +

aff£ff -

d i r 1

(* = 1,2, ...) .

По условию теоремы инвариантное подпространство J R a , соответствующее группе о собственных значений матрицы £/, циклическое. Поэтому минимальный многочлен этого подпространства есть многочлен второй степени с коэффи­ циентами, определенными по формулам Виета:

фо (Я) =

W +

+ а?>,

a ff =

— (iff +

 

Примем

GC/o = ССу°'

(у = 1 , 2 ) .

Тогда

фа (^ )= ф 0(Х),

т.е. фа (А,) — минимальный многочлен подпространства Ra* следовательно,

фа(Лд) = О,

Положим

 

 

 

 

& -

Кль

 

 

 

(З-'З)

где аа — некоторый двумерный вектор-столбец.

 

 

Легко видеть,

что

при таким

образом

определенных

<р0 (к) и gfS1 первое равенство первой

группы соотношений

(3.12) выполняется тождественно. Действительно,

 

 

2

(Л,) Ms$ ?

= 2

/С,<рст (Л8) MsKoa<s=Ка<р<>(Ли) а„= 0 .

 

 

 

S = 1

 

 

 

 

 

 

Второе же равенство определяет £га-

 

 

2; i =

Допустим,

что

уже

построены

& +

«58 (/*=!»

= 0,

1.......Л — 1). Покажем,

что

1 группа

равенств

(3.12) вполне определяет

и «{S1.

 

 

 

 

Имеем

 

 

 

 

 

 

 

 

 

2

/<5ф„ (AJ Щ ,1? — — (о$1/ +

£») й°

+

d'"

"•

s«=l

 

 

 

 

 

 

 

 

 

Здесь d^k~l]

является

известной величиной.

получим

Принимая

во

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

(3.13),

2

К&а (AJ

------\9КоЪ + о&'К*) «о +

 

 

S—1

 

 

 

 

 

 

 

 

 

 

или,

обозначив

 

 

 

 

 

 

 

 

 

 

£бо=

(Л0ЯаД<т),

 

 

 

 

 

 

2 к ,ф, (л5)

При умножении обеих частей последнего равенства сле­ ва на М это равенство распадается на р независимых матричных соотношений

фа (Л5) Qs? = — MsKa££aa>^ +

^

(3*14)

(s — 1, • • • , р),

где

<2* = М$У.

Обозначив

будем иметь

Q S,*1= ■

так что если удастся из (3.14) определить Q?], то Легко по­ строить и искомый вектор |[{? по формуле

При s ф о MsKa = 0, а фо (As) в силу условия (3.2а) — невырожденная матрица. Поэтому

ОЙ1 = <р7' (Л,)м /с * - '1

(S фсг).

При s == а Ма Ка = Е<&, фа(Аа ) =

0. Поэтому из (3.14)

следует, что

 

£ос4ч = Л М ? -11.

(3.15)

Так как Ra — циклическое подпространство, то соглас­ но лемме 1.2 существует такая столбцовая матрица аа , ко­ торой отвечает линейно независимая система векторов аа, Астя0 . Пусть аа выбрана из этого условия. Тогда векторы Ааво и аа линейно независимы и, значит, &£а — невырож­ денная матрица. При этом матричное равенство (3.15) раз­

решимо относительно

Решая его, находим

Что касается субматриц QES» то здесь имеется известный произвол. В качестве Qo3 может быть взята произвольная квадратная матрица второго порядка, имеющая достаточное

число производных. В частности, можно принять Qo*d = = 0.

Зная |SCT и aSJ1, легко определить и g^no формулам (3.12). Таким образом, изложенный метод позволяет последо­ вательно построить члены формальных рядов (3.5) и (3.6), с помощью которых представляется формальное решение уравнения (3.1) в форме (3.3), (3.4). Тем самым теорема до­

казана.

 

 

 

 

3.2. Общий

случай.

 

L] а) матри­

Т е о р е м а

 

3.2. Пусть на сегменте [0,

цы Ak (т), Bk (т)

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

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

всех порядков;

б) собственные значения матрицы U (т) =

= А~1 (ir) В0 (т)

разбиты на

р групп

А,<°> , ..., Ag

\

и

ka

= п)

так,

что

 

(а = Ц

...>р, ^

 

 

°

1|^ e,(t) —Я,?*(х)|>0

(3.16)

 

( я ф о ;

i =

l ,

, ko,

j = 1,

, fcs),

причем соответствующие этим группам инвариантныеподпространства /?х, /?2, Rp являются циклическими под­ пространствами п-мерного пространства R.

Тогда формальное решение уравнения (3.1) может быть представлено в виде

где <7о (а = 1, 2 ,...,/?) —скалярные функции, удовлетво­ ряющие уравнениям

+ Otia(T, е) k

+

-f-

 

 

dta

 

Ч ^0

 

I

/

 

“Ь

a k g —\[О fa , £) — ^

Ь

 

+ а*ао(т, 8)^ = 0

(<т = 1

/7) (3.18)

а £/а (т, е), а/0 (т, е) — соответственно векторные и ска­ лярные функции, представляемые формальными рядами

lie (Т, е) = 2 8*1}? (г),

а„ (т, е) =

£

(т). (3.19)

Д о к а з а т е л ь с т в о .

Применяя

метод,

использо­

ванный в § 2, подставим выражения (3.17) и (3.18) в урав­ нение (3.1) и приравняем коэффициенты при qa, d{[a -,

Аналогично тому, как были получены формулы (3.8), отсюда находим

 

(3.20)

№ ' = и & + « 5 ® + «

- < # -«,

( / = 1 ...........ka; * = 0,1,

а™ 0).

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

достаточно показать, что

при соответствующем выборе членов рядов (3.19) матрич­ ные соотношения (3.20) выполняются тождественно.

Первое соотношение (3.20) в развернутом виде представ­ ляется так:

М _ 7/tW] I

S2a — u Zio - r a iaSia»

$ ! = у $ ч д а ,

(3.21)

о = у ^ + « д а .

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

u ha~l, Uka-‘2, ...,

£/, Еп

И СЛОЖИМ. ПолучИМ

(У‘о + ai°o]C/V-4-

+

с$1ю и + a$,E„) = 0. (3.22)

Заменив последнее равенство системы (3.21) равенством (3.22), далее вместо системы (3.21) будем рассматривать эквивалентную ей систему

фо (U)йа = о,

1

$ , „ = £ / $ ? + « д а

о = 1 ,

, *0 - 1 ),

где

фа(А-) =

1-f-

+ abq—I “Ь

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