книги / Матричное исчисление с приложениями в теории динамических систем
..pdfЕсли 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