книги / Матричное исчисление с приложениями в теории динамических систем
..pdfсмотрим матрицу K(t) = X{t)CZ(t)i где X(t) — единственное реше ние матричного уравнения dX/dt = P(t)X (X(t0) = Е), В — еди ничная матрица; С = (с,, с2, си) — постоянная п х «-матрица;
Z(t) = codiag
’ W X c J
— непрерывно-дифференцируемая на [a, °°) диагональная матри ца; 0,(0. •••. 0„(О — непрерывно-дифференцируемые на [а, <») ве щественные скалярные функции. Легко видеть, что столбцы
Х,(0» .... Kn(t) |
матрицы |
K(t) |
удовлетворяют |
условию |
|| А:.(ОН — на [а, <») ( i — 1, 2, |
«). Замена переменных |
|||
|
x=K(t)y(f) |
|
(18.6.2) |
приводит уравнение (18.6.1) к уравнению с диагональной матрицей
^ = Д (0 у, |
(18.6.3) |
где |
|
М О = z-'(t) 422- = diag (Х,(0. М О ...... МО)- |
|
В соответствии с (18.6.2) и (18.6.3) имеем |
|
/ |
|
дс(0 = /ВДехр J Л(г)с/ту0 (у0 = у(*0)). |
(18.6.4) |
*0 |
|
Построим пучок решений уравнения (18.6.1), берущих начало |
|
внутри и на поверхности эллипсоида |
|
Н^'х0) « р2, |
(18.6.5) |
где HQ— постоянная невырожденная п X «-матрица, столбцы кото рой имеют эрмитову норму, равную со.
Совокупность вектор-функций (18.6.4), ограниченная условием
(Уо> Уо) < Р2. |
(18.6.6) |
определяет искомый пучок решений уравнения (18.6.1). |
|
Подставляя (18.6.4) в (18.6.6), получим |
|
х*Л“,(/)х(/) < р2, |
(18.6.7) |
где |
|
t
B(t) = А(/)ехр [ 2$ Re Л(т)с/т] **(т).
Отсюда, поскольку со0(*)<а>, как это 'следует из (18.6.9) и (18:6.11) при всех t> tQ, то на [*0, ») выполняется условие
(С -Чфс"(0. G-'0)x’( t ) ) * p \
а это означает, что решение дс°(*) при всех t > t 0 не покидает пределы рш-трубки (18.6.12); последнее свидетельствует об ус
тойчивости невозмущенного движения (тривиального решения уравнения (18.6.1)). Теорема доказана. ■
Т е о р е м а 18.6.1.2. Если в какой-нибудь точке |
G [£0, ®) вы |
полняется неравенство |
|
1 2 <=*Р [1 M ' I K 'I - 'о)1 > 1. |
(18.6.13) |
О- 1 |
|
то невозмущенное движение (тривиальное решение системы (18.6.1)) неустойчиво.
Д о к а з а т е л ь с т в о . Пусть в некоторой точке е [*0, «>) вы-.
полняется соотношение (18.6.13). Тогда в соответствии с (18.6.9) со0(/,) > о» и повторяя рассуждения при доказательстве теоремы
18.6.1.1, можно показать, что какова бы ни была ря-трубка (см. (18.6.12)), вне этой трубки в момент окажутся некоторые из тех решений JC°(0 системы (18.6.1), которые в момент t0 находились внутри или на поверхности эллипсоида
(G -4.t,)x'(to)'G -'(t„)*•('„))« Р*.
что и доказывает теорему. ■ Без труда может быть доказана и следующая
Т е о р е м а 18.6.1.3. Если при t> tQ
Ц0(0 = max Re \ a(t),
О
то невозмущенное движение устойчиво.
Т е о р е м а 18.6.1.4. Если на [f0, °°) |
|
ц(0«*—Ъ ( |i(0 — maxp0), |
(18.6.14) |
О |
|
где Ъ — положительная постдянная, то невозмущенное движение (тривиальное решение системы (18.6.1)) асимптотически устой чиво на [Г0, £»).
Д о к а з а т е л ь с т в о . Так как при условии (18.6.14) выполня ется и неравенство (18.6.11), то по теореме 18.6.1.1 невозмущен ный процесс устойчив.
Пусть теперь xe(t) — произвольное нетривиальное решение уравнения (18.6.1), удовлетворяющее неравенству
(С -‘(1„)*•(<„). O~'(t0)x'(t0)) S р*
где G(t) — некоторая матрица из класса К%. Докажем, что
lim||x°(OII=0. |
(18.6.15) |
/-♦00 |
|
Согласно (18.6.4) имеем |
|
11*4011 « ШОНехр | Л(т)Л||у0||. |
(18.6.16) |
Здесь ||ЛТ(/)|| =со\ЛГ. В силу этого из (18.6.16) при всех t > t0 по лучаем
|
II АОII < |
exp [-b(t - <0)]||у0||, |
|
|
|
откуда и следует (18.6.15). Теорема доказана. ■ |
|
|
|||
Т е о р е м а |
18.6.1.5. Если на [<0, «) выполняются неравенства |
||||
|
t |
|
|
|
|
|
|ехр J Sp А{х)<1х\ 5» т{> 0, |
|
(18.6.17) |
||
|
t |
|
|
|
|
— f |
$ ^(х)йх^ —b< 0 |
(* = 1, 2,..., |
п), |
(18.6.18) |
|
где p.j(/),..., р„(0 — все собственные значения |
эрмитово-сим- |
||||
метризованной матрицы AH{t) = |
^ (A(t) + A*(t)), |
то невозму |
щенное движение (тривиальное решение системы (18.6.1)) асим птотически устойчиво.
Д оказательство. В преобразовании (18.6.2), где K(t) =
— X(t)CZ(t), примем
z (0 = О) diag j p c j ’ —»||*cj|J
Тогда элемент Хо(0 диагональной матрицы Л(/) в (18.6.3) опреде ляется формулой
K V ) = £ ln ||ХС„|| (а = 1, 2 ,.... п). |
(18.6.19) |
В соответствии с (18.6.2) и (18.6.3) имеем
x(i) = К(Oexp J A(x)dxy(t0). |
(18.6.20) |
Матрицы X l(t) и Z ‘(f) существует на [/0, «>). Существование
первой матрицы следует из формулы Остроградского—Лиувилля в силу (18.6.17), а второй — из ограниченности каждого решения си
стемы (18.6.1). Следовательно, на [/0, °°) матрица А-1(г) также
существует, а функция V(t, х ) = (A-1(/)JC, K~l(l)x) = ||у||2 явля ется положительно определенной. Полная производная по / от этой функции в силу уравнения (18.6.1) представляется в виде
= 2 2 X J y J \ |
(18.6.21) |
0—1 |
|
где уа — элемент столбцовой матрицы у(1).
Интегрируя (18.6.21) вдоль решения уравнения возмущенного движения (18.6.1), получим
t П
V(t, х) = K(l0, x(t0)) + \ X 2 Х»Ы *Л . '. о - 1
Из уравнения (18.6.3) имеем d y j d t = \ ьуа (о = 1, 2,..., л), откуда
ы |
\ 2 |
|
|
|
В силу этого |
|
|
|
|
|
п |
t |
I W I 2 |
|
х) = V(tQ,xVо» |
1 + 2 |
(expJ 2 Xadx — 1) |
(18.6.22) |
Имеют место неравенства ^ ^а(0 ^ И'тах(^) (С = 1, 2, ..., п), (18.6.23)
где цт|П(0 и цтах(/) — соответственно минимальное и максимальное
собственные значения матрицы AH(t). Из соотношений (18.6.18) и (18.6.23) следует, что на [*0, ») выполняется неравенство
— |
С\ a(x)dx ^ —b < 0 |
(о = |
1, 2, ..., л), |
(18.6.24) |
t |
tQJ |
|
|
|
|
<Q |
|
|
|
и в соответствии с этим из (18.6.22) получаем |
|
|||
|
7(/, x(t)) «а v(t0, x(t0)) |
при |
Vf Е [/0, оо). |
|