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

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

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

Если L0 — невырожденная матрица, то А — также невырожденная матрица. В самом деле, умножая матрицу А справа на матрицу

О О

и»

/ =

 

 

 

Ят 0

0

 

получаем

 

0

0

0

 

А>

AJ =

0

А)

0

0

 

 

= L,

 

0

0

EQ 0

 

Ек-У ^ -2

Еу

Е0

Отсюда А LJ, так как

J~l = J. В соответствии с этим det А =

= det L det J. Но det L =

(det L0)k>

 

det/= (-l)**-*)"1•

«

( _

^(fc-A+O^wdetEm=

 

 

__ ^_j^m2|l+2 + ...+(ifc-1)1 _ (_k(k- 1)

Поэтому

 

 

 

 

det A — (—I

_ l)- (det L0)k>

откуда следует, что матрица А одновременно с матрицей LQ вы­

рождена или невырождена.

Предполагая, что L0 — невырожденная матрица, построим об­

ратную А~1. Пусть

II т

( * 1 1

Nl2

Nlk)

N21

*12

*2*

*41 * 4 2

Имеем

Отсюда

 

 

 

 

 

II

JI

SS

н»

1

1 II

*2k"

**■О

1

II

А>*»

II

 

 

 

 

 

=

... = Nk2 = 0

 

Итак,

 

 

 

 

 

н 4° ^2к-2 ~••

f III

II

и

•• II©

- О

1 <*1

1

h

И ^к-2^0

J-,Q ^j-^o

 

 

 

_/ _1 [ г - 1

L~lL /Г 1

А~1=

0

 

0

V

 

 

 

 

 

0

 

V

0

 

^о1

 

0

0

Зная А 1, легко привести (7.2.9) к форме Коши:

< £ = U x + A-if,

где

*31о 0

0

0

(7.2.10)

 

А)1 ^ 2

Е т

0

0

0

0

0

< 1*1l>1

0

0

~VA 0

0

0

З а м е ч а н и е 7.2.1. Характеристический многочлен системы (7.2.8):

 

Д(Х) =

| L0Xk + L{kk~l +

... + Lt \.

Характеристический многочлен системы (7.2.10):

1

otI 1— 1

- Е Л

_/ -1 7

1

 

ТП

•Н-

д а > =

Е т

 

L - E n k

0

 

 

 

 

1 сП ь___1 0

0

0

~ ЕтХ

Оба эти многочлена одной и той же степени тк и, более того

точностью до множителя ± | L01 совпадают.

б) В торой способ, Положим

 

 

0

0

0 '

0

£

0

0

0

0

Е

0

0

0

0

L0

'

0

Е

0

0

*

 

а

 

0

0

Е

0

 

 

dg_

в =

 

 

 

»

Х =

dt

0

0

♦*•

 

0

Е

 

 

dk~ lg

 

 

 

2

 

 

 

 

 

~ L k - 1

~ L i

 

 

d t * '1

В этих обозначениях уравнение (7.2.8) можно представить так:

A < £ = B x + f ,

или в нормальной форме:

 

 

 

 

fjr =

их + А-Ч,

 

 

где

 

 

 

 

0

E

0

0

0

0

0

E

0

0

0

0

0

0

E

_l ~ 1/

—L~l T

""A) *Afc-2

_г —1 1

_f ~l /

 

J-‘0 t-’k - l

^0 ij2

J~'Q

§ 7.3. Норма матрицы

Нормой прямоугольной матрицы А (в частности, и столбцовой, и строчной матриц) называется неотрицательное действительное число ||Л||, удовлетворяющее условиям:

1)

|| А\\ > 0,

если А ^

0 и || Л|| = 0, если А = 0;

2)

||Л + * ||< |Н || +

И ||;

3)

||ХЛ|| =

| Я11| ^4||

(X — число из поля Ж);

4) ||Л £ И |И |||/? ||.

Условиям (1)—(4) можно удовлетворить многими способами. Например, можно положить

Il^ll = max 2

/к

И Л И

|М|| = max 2 \ ajkI*

 

*

И Л И

1

 

Mil

2

2 K * i2

 

j.t

Норма, определенная последним равенством, называется эрмито­ вой (в случае вещественных а|у- — евклидовой).

Отметим два свойства нормы матрицы:

 

К /И М 1 1

(7.3.1)

|Х ,.|«М Н ,

(7.3.2)

где Xj — собственное значение матрицы А.

Свойство (7.3.1) очевидно, а (7.3.2) можно установить так. Пусть Xj — собственный вектор, отвечающий собственному значе­

нию Ху. kjXj = Axj. Переходя к нормам, получаем

1*/И1*у11 - М*;11 MIIMyll-

Отсюда, так как х } Ф 0, следует (7.3.2).

§ 7 .4 . Матричные ряды

Рассмотрим последовательность матриц Ср С2, ..., Сп, ...

(Ср = (c$J>)) одного и того же типа т х п . Пределом этой последо­ вательности называется матрица

С — lim Ср = (lim с\р)) (i — 1, 2, ..., т; j = 1,2, ..., л),

р -»ао р -*« о

если, конечно, она имеет смысл, т.е. существуют пределы числовых рядов последовательностей

с<у, с ? } ,..., с\р) , ...

(i'= 1 , 2 ,..., т ; j =

1 , 2, ..., п ).

Пусть Uv U2, .... Un, ...

матрицы одного

и того же типа.

Матричный ряд

 

-f1и р + ...

(7.4.1)

^ + 6Г2 +

называется сходящимся, если существует предел последовательно­ сти его частичных сумм S[t S2, ..., Sn, ... (£р = + ... + Up). Предел этой последовательности называется суммой ряда (7.4,1).

Наряду с матричным рядом (7.4.1) введем в рассмотрение ряд,

составленный из норм матриц Up = (uW)

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

 

+ \W2\\ +*•• +

+••*

(7.4.2)

Если ряд (7.4.2) — сходящийся, то матричный рад (7.4.1) — также сходящийся. Действительно, так как при всех i и / | uW \ <

(см. §

3), то согласно признаку сравнения скалярных рядов все ря-

со

 

 

 

ды ^

и\р) — абсолютно сходящиеся. Следовательно, ряд

 

p=i

/

 

 

 

00

 

 

Рт1

2«Й>

 

 

р = 1

 

также сходится.

 

 

§ 7.5. Теорема существования и единственности решения

 

однородного векторно-матричного уравнения

 

•Т еор ем а 7.5.1. Если U(t) непрерывна на t0^ t ^ Т ,

mo суще­

ствует единственное решение уравнения

 

 

% =U (t)x,

(7.5.1)

 

 

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

(7.5.2)

 

x(t0) — с.

Это решение непрерывно и дифференцируемо на t0*Zt**T.

7.5.1. Существование решения. От дифференциального уравне­ ния (7.5.1) перейдем к соответствующему интегральному уравнению

t

 

 

x(t) = с + J U(s)x(s)ds,

/0 < f ^ 71.

(7.5.3)

*0 Уравнение (7.5.3) будем решать методом последовательных при­

ближений. Пусть х(ЛГ_1) и х(лг> — последовательные приближения

уравнения

(7.5.3). Тогда

 

 

 

t

 

 

x ^ ( i )

= с + J U (s)x^~ l\s )d s

(N = 1, 2, 3, ...).

(7.5.4)

 

 

 

Отсюда, полагая х ^ — с, при N = 1 будем иметь

 

 

t

 

(7.5.5)

 

хW — *№> = J f/(s)x^f/5,

а при N > 1

* « _ * < « - ■ ) = $

(7.5.6)

Из равенств (7.5.5) и (7.5.6) вытекают следующие неравенства:

||x(l)_ * (0)|| ^ j ||t/(s)||||x<0>||</s,

(7.5.7)

||л(лг)_х(л(-1)|| ^ $ ||f/(5)|[||x(Ar- ,) - ^ JV- 2>||£/s.

(7.5.8)

Нормы векторов и матриц определим так:

М = 2

l*il.

Il^ll = тах 2 IaikI-

i=1

 

/

k

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

матрицы U на

замкнутом промежутке

[/0, Т] следует ограниченность ее нормы. Пусть m — max ||i/(/) ||; ге[/0, г\

тогда из (7.5.7) и (7.5.8) вытекает

t

||*(1) - х(о)ц ^ J m\\c\\ds = mc(t — Z0),

цх(2) _ х(1)ц ^ J m2||c ||( s - t0)ds = m2c2 U ~ tQ)2

и далее по индукции

(mTt)N

 

N\ *

(7.5.9)

{Тх = Т - 1 0).

Ряд mNT*/ N\ — сходящийся. Действительно, используя признак ЛГ«1

сходимости Даламбера, будем иметь

mN+lTK+lN, mT.

О при N-+ оо.

Ш+ l)\mNT^ N + l

Из сходимости рада (7.5.9) следует равномерная сходимость на [г0,Г] рядов

tf=l

ЛГ=1

 

и, наконец, ряда

 

 

со

 

 

*№) + 2) (*<*) -

(7.5.10)

N = 1

 

 

Поскольку ряд (7.5.10) сходится равномерно на [tQ,T], то суще­

ствует предел последовательности

хО),

lim *W (0 = x(t)

N-*«о

 

 

(XW есть сумма первых N + 1

членов ряда (7.5.10)). В силу рав­

номерной сходимости последовательности х^°\ *0), ...}x^N\ ... в левой и правой частях равенства (7.5.4) можно перейти к пределу при JV-* оо. В результате получим

x(t) = с +

t

 

J U(s)x(s)ds.

(7.5.11)

 

*0

 

Функция x(t) как предел

равномерно сходящейся

последова­

тельности непрерывных функций сама является непрерывной фун­ кцией, и, как это видно из (7.5.11), она дифференцируема и удов­ летворяет уравнению (7.5.1) и начальному условию (7.5.2).

7.5.2.Единственность решения.

Лем м а Г р о н у о л л а — В еллм ан а. Если

с{> 0, и(1)^ 0,

v(t)& 0

(u(t), v(t) Е С(а,Ь)),

 

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

t

 

 

 

 

 

u(t) < Cj +

J u(s)v(s)ds

(a < t0>t < B)

(7.5.12)

следует неравенство

t

u(t) ^ c^xp J v(s)ds.

Д о к а з а т е л ь с т в о . Из (7.5.12) имеем

u(f)

t

Cj +

Умножим обе части последнего неравенства на v(t) и проитегрируем от t0 до /:

 

u(t')v{t’)dt'

t

v(s)ds.

 

S

^ ^

 

I

 

 

 

+ ^ u ( s ) v ( s ) d s

 

 

 

Отсюда

 

 

 

 

 

t

 

t

 

t

 

In [c, + 5 M(s)u(s)f/(s)] | {0 ^

j v(s)ds,

In u(t) — In cY< $ v(s)ds

и, следовательно,

 

 

t

 

 

 

 

 

 

 

 

u(t) sS Cjexp J v(s)ds.

 

Лемма доказана. ■

 

 

*o

 

 

 

 

 

 

 

Пусть х и у — два решения уравнения (7.5.1), удовлетворяю­

щие одному и тому же начальному условию (7.5.2). Тогда

 

t

 

 

 

i

 

х — с + $

U(s)x(s)ds, у — с + J U(s)y(s)ds.

 

^0

 

 

 

*0

 

Вычитая эти выражения друг из друга, получаем

 

 

*

 

 

 

 

 

х у = $ U(s)(x — y)ds.

 

Отсюда

*0

 

 

 

t

 

 

 

 

 

 

 

 

 

II* - У\\ * 5

l|tf(s)llll* - y||rfs

(7-5-13>

и, тем более,

*0

 

 

 

 

 

 

 

 

t

II* - у\\ < с, + $ ||£/(s)||||* - y||</s *0

для любого положительного числа с,. Используя лемму, получаем

t

||х - у || < с,ехр $ ||t/(s)||tfs. 'о

Это соотношение справедливо для любого сколь угодно малого по­ ложительного с,. Поэтому ||* — у|| = 0 и, значит, х = у.

Единственность можно доказать и другим путем, не прибегая к лемме Гронуолла—Веллмана. Так как х и у — два непрерывных и дифференцируемых решения уравнения, то их нормы на [*0, Г] ог­

раничены. Пусть т1=ш ах||х — у||. Из (7.5.13) находим

II*- У\\ ** mmx(t - t0).

Найденную оценку для \\х — у|| снова подставим в (7,5.13). Полу­ чим

II-» - УН«

т2тit Г1(<лT

_

-

(0)2,

Повторяя этот процесс, будем иметь

 

 

 

т . т "

 

Полагая «-»<» (при

t< оо),

получаем ||х —у||<0. Значит,

|| х — у|| = 0 и х — у.

Лемма доказана. ■

§7.6. Фундаментальная матрица системы

7.6.1.Решение матричного дифференциального уравнения.

Построим теперь решение матричного уравнения

^ - = U ( t ) X

(7.6.1)

при начальном условии

X(t0) = С,

(7.6.2)

где С — постоянная невырожденная матрица порядка п. Эта задача эквивалентна построению п решений векторно-матричного уравне­ ния

 

= U(t)x)

(7.6.3)

соответствующих п начальным условиям

 

*(*о) — ci

(i = 1, 2, ..., л),

(7.6.4)

где ct (i = 1, 2, ..., л) — столбцы матрицы С.

По теореме существования и единственности решения уравне­ ния (7.6.3) каждому сг соответствует единственное решение xt(t).

Ясно, что матрица, составленная из этих решений, а именно мат­ рица X = (xt х2 ... *„), представляет собой решение матричного

уравнения (7.6.1) при начальном условии (7.6.2).

7.6.2. Формула Остроградского—Лиувилля. Продифферен­ цируем определитель матрицы X(t) = (х, х2 ... хп), представляю­

щей решение уравнения (7.6.1) при условии (7.6.2):

 

* n (0

*i*(0

*i„(0

<ша = у

dx^tt)

dxik(<)

d*jnW

dt

dt

dt

dt

j = 1

 

 

 

Ho

*„i(0

*И*(0

*n„(0

 

 

 

d x ik

"

(;', &=

1, ..., л),

dt

 

 

 

 

где ы -5 (s = 1, 2 ,..., л) — элементы /-й строки матрицы U. Поэтому

ft.

II

 

* 11

 

* l / f c

* 1 *

n

n

n

 

n

= 2

2

2

5fc

2 **jsXsn =

y = i

5 = 1

5 = 1

 

5 = 1

 

* « 1

 

 

*nn

 

* 1 1

* 1 *

* 1 »

 

n n

 

n n

= 2 2 и Уж * , i

 

=2 2 > y , ( o a y , m .

У = 1 5 = 1

 

7 = 1 5 = 1

* « l

X

n k

где Ь!г — символ Кронекера, и, далее,

Js

 

 

d- ^

= i u

ij{t ) \ x \ = S v u \ x \ .

 

1

 

Интегрируя последнее соотношение, получим формулу Остроград­ ского—Лиувилля:

t

| x ( t ) | = | *(/<,) | exp $ Sp U(t)dL

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