книги / Матричные и асимптотические методы в теории линейных систем
..pdfВышеизложенное позволяет сформулировать следующие условия устойчивости.
Т е о р е м а 4.1. Если |
|
J - 2 е2“а « «-'•> а« (0 < (a® (t) |
(t g р0, Г)), |
то невозмущенный процесс {решение уравнения (2.2)) устой чив на промежутке U0, Т).
Т е о р е м а 4.2. Пусть на промежутке [t0, Т)
и |
а (0 < to(t) |
|
|
р (0 < 0 |
(р (0 = шахора (О)- |
Тогда невозмущенный процесс {решение уравнения (2.2)) ус тойчив на [tQ Т).
Те о р е м а 4.3. Пусть на промежутке U0, Т)
а(0 <со(0
и
Ро (0 ^ ^ |
(Ро (0 = niaXtf Re {£)). |
|
Тогда невозмущенный процесс {решение уравнения (2.2)) |
||
устойчив на U0, |
с»). |
|
Т е о р е м а |
4.4. Пусть на промежутке [ у оо) |
|
|
|
а (/)< (о (0 |
“ |
|
( i( 0 < — |
где Ъ— положительная постоянная. Тогда невозмущенный процесс {решение уравнения (2.2)) асимптотически устойчив на U0,oo )•
Наконец, призедем еще одну теорему, определяющую условие существования конечного промежутка времени, на котором процесс устойчив.
Те о р е м а 4.5. Если на промежутке [/0, у
а(0 < со (*)
“ |
(*»(<.)<о, |
(4.2) |
то существует |
конечный промежуток [/0, |
Т) cz [t0, У , |
на котором невозмущенный процесс {решение уравнения (2.2)) устойчив.
Д о к а з а т е л ь с т в о . Условие (4.2) эквивалентно неравенству
Р (У <0,
а последнее соотношение влечет за собой, в силу непрерыв ности функции р, (0. выполнение неравенства
р( / ) < 0
впределах некоторого конечного промежутка U0. Т), и по тому, согласно теореме 4.2, невозмущенный процесс (реше
ние уравнения (2.2) устойчив.
§ 5. Случай стационарной системы
Рассмотрим процесс, представленный уравнением
- Т Г - и * . |
(5-1) |
где U— постоянная квадратная матрица порядка п. Для простоты ограничимся случаем, когда U — матрица прос той структуры. В этом случае фундаментальную матрицу системы (5.1) можно представить в виде
|
X = Ге>('-« Г " 1, |
|
|
|
|
где Г — квадратная |
матрица, |
составленная |
из нормиро |
||
ванных собственных |
векторов |
Tlf Г2, ..., |
Гл |
матрицы U, |
|
J — жорданова форма матрицы U, в данном случае — диаго |
|||||
нальная матрица, по диагонали которой расположены |
соб |
||||
ственные значения v1, v2, ..., v„ матрицы U. |
|
|
|
||
Пусть со (t) — заданная положительная |
функция, |
опре |
деляющая класс п х л-матриц Ад, a a (t) — некоторая непрерывно дифференцируемая функция, удовлетворяющая на заданном промежутке U0, Т) условиям *)
оь (t) |
® (0> а (to) ” ® (^о)' |
В формуле (2.4), представляющей матрицу преобразова ния уравнения (5.1) к уравнению с диагональной матрицей А, положим С = Г. Тогда в соответствии с (2.10)
*) Если 0(/) — сама непрерывно дифференцируемая функция, то в качестве а(/), которая определяет норму столбцов матрицы преобразо вания уравнения (5.1) к уравнению с диагональной матрицей, можно принять заданную функцию ш(/).
устойчивости на промежутке U0, оо):
maxff Rev0 < 0 .
Условие асимптотической устойчивости на U0, оо):
max* Re v„ < — b |
(b> 0). |
§ 6. Об устойчивости на конечном промежутке нелинейного процесса по линейному приближению
Ниже устанавливаются некоторые условия устойчивости на конечном промежутке [/„, Т) (Т <; оо) процесса, пред ставленного тривиальным решением (я == 0) векторно-мат ричного уравнения
-lL |
= U(t)x + h(t,x), |
(6.1) |
где U — квадратная |
матрица порядка л, |
непрерывная на |
[/0, Т), h — столбцовая матрица, |
элементы которой — не |
||
линейные функции отклонений xs — таковы, что |
равномер |
||
но по t на промежутке U0, Т) |
|
|
|
Hm h{t,x) |
= |
0. |
(6.2) |
*-*■0 || х || |
6.1.Теоремы об устойчивости по линейному приближе
нию. Пусть (о (/) — заданная функция, порождающая класс
л х л-матриц К ! а К (t) = (/(i /С2 ... Кп) — невырожден ная и дифференцируемая на [10, Т ] матрица преобразования уравнения линейного приближения
-ТВ— " ® *
к уравнению
с |
непрерывной диагональной матрицей |
А — diag (ЯА, |
|
Я2, |
...» А.,,) при условии |
|
|
|
II Я /(/)[ = «(/) |
( / - 1 , 2 ........... |
л), |
где a (t) — положительная функция, непрерывно дифферен цируемая на U0t Т), причем a (t0) = со (t0).
Те о р е м а 6.1. Пусть на промежутке [/<>> Т)
а(*)<а>(0
и
*1 ( 0 = шахоТ=Г7 |
(6.3) |
где b — положительное число. Тогда невозмущенный про цесс (тривиальное решение уравнения (6.1)) устойчив на про межутке U0, Т).
Д о к а з а т е л ь с т в о . Введем в рассмотрение мат рицу
0(t) = К (t) <а(0a(t)
Ясно, что G (/) есть матрица класса /(д. Для доказатель ства теоремы достаточно показать, что все решения урав нения (6.1), удовлетворяющие условию
( / Г 1(У *0, К~' (tc)х0)< |
ра, |
(6.4) |
при V* £ Uo. Т) удовлетворяют условию |
|
|
(CT'(t)x, G-'(t)x) = ^ r (fC\t)x, |
К~‘(t)x)<p*. |
(6.5) |
В уравнении (6.1) произведем замену переменных
х = Ку.
Получим
-%- = Л(t)!/ + M(t)h(i,Ky),
где М = /С-1. Функция
V(t, х) = (K~l (t)х, /С-1 (t)х) =Цу ||а
является положительно определенной. Ее производная по t, вычисленная в силу уравнения возмущенного процесса (6.1), представляется в виде
~ |
= |
2 |
2 Re К | УсIs + 2 Re (у*МН), |
(6.6) |
|
ас |
|
0=1 |
|
|
|
где уа (ст = 1,2, |
.... п) — элементы столбцовой матрицы у. |
||||
Интегрируя (6.6) вдоль решения уравнения возмущен |
|||||
ного процесса, |
получим |
t |
|
||
|
|
|
t |
|
|
V(t, X) = V(t„, x0) + |
f 2 2 Re Ла 1 [ 2 Л + |
f 2 Re (y'Mh) dt. |
|||
|
|
|
U c |
t. |
|
Преобразуем интеграл
(
/ = f 2 2 R e * B|0,|*<#. t O’
Так как A (t) — диагональная матрица, векторно-мат ричное уравнение относительно у можно представить в виде следующей системы уравнений первого порядка:
^hr = К9а + M0h ( о = 1 , 2 , . . . , п).
где Ма — строка а матрицы М .
Отсюда, переходя к дифференциальному уравнению от носительно модуля | £/а | и интегрируя это уравнение, полу чим
г
|
J 2Re*,0rf* |
|
J |
2Re\edx |
||
lyoj2 = ef* |
|
|
ta |
Re {уаМаК)dt' f |
||
где Уао = |
Уа (g. |
|
|
|
|
|
Тогда |
|
V |
|
|
|
|
|
|
|
|
|
|
|
Г |
|
J 2Re\0rfT |
|
|
|
|
^ — j S 2 |
R |
e е*° |
\у0о|2dt’ -f- |
|
|
|
U ° |
|
V |
|
г |
|
|
|
|
|
|
|
||
<»' |
|
^ 2Re7tadx /<' |
-{ 2 R e M * |
\ |
||
Н~ J 2 2 ReX0ef" |
|
2е *• |
|
Re (уаМ0Ь) dt"J dt' = |
||
|
|
|
f 2*eКdx |
|
|
|
|
|
|
|
\yсо |
+ |
|
|
|
V |
|
Г |
|
J*e |
|
|
|
|
|
||
|
|
J 2Re?.0dx |
— J 2ReK0dx |
Re (yaM„h) d f |
||
+ |
2 |
e>’ |
j 2e |
|
||
|
‘ |
J IReKodz |
- J 2Rel(,i/t |
|
||
— J S « '- |
2e |
'• |
Re (yaMaK)dt' |
|
|
J 2ReXadx |
|
|
|
= |
2 ( е ‘- |
-1)|» «Р + |
|
|
|
G |
|
|
t |
|
|
|
|
J 2ReX0dx |
* |
- $ 2RеЛ,0Л |
R£(gcMah)dt'— |
|
+ 2 е*° |
] 2е |
'» |
||
9 |
и |
|
|
I |
|
|
|
|
|
|
|
|
— 2 |
f 2 Re (ffoMoh) dt\ |
|
|
|
a |
i |
Принимая во внимание полученное выражение для ин теграла /, будем иметь
V (It, *) =
|
|
|
С 2ReX0dx |
I „ |
.2 |
|
1 |
|
= v (#0, *«) |
I + |
2 |
(е*° |
- 1)-{j^jjr + |
(* - |
/о)Ч> (/. ?)], |
||
|
|
|
|
|
|
|
|
(6.7) |
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
|
|
|
|
|
9 |
л |
I 2ReAdT |
|
|
|
♦ Й = |
« - ц ц Д |
i Re |
|
|
Л4Л1 dV- |
|||
|
|
|
|
tn |
|
|
|
|
Покажем, |
что |
на |
промежутке |
t0 < |
t < |
Т |
равномерно |
|
по i |
|
|
lrnnj>(rf, у) —0. |
|
|
(6.8) |
||
|
|
|
|
|
у-* 0
Действительно,
С другой |
стороны, |
|
|
|
||
|1уНа |
_ |
|
|
|
|
|
1|УоП3 |
|
|
|
|
|
|
|
t |
|
|
|
t |
|
= ехр |
5? |
2Re w |
dt' + 1 im p Re v m # |
< |
||
|
|
I |
|
Г |
|
|
< e x p |
j |
£ 2 Re M f ' + J |
Re (R*MA) df' |
< |
||
|
|
Лл |
^ |
tn |
|
|
|
|
|
< e x p |
|
|
JIM |
|
|
|
|
|
|
Ik II |
Так как К — матрица, невырожденная на замкнутом промежутке U0, 74, то М на этом промежутке — ограничен ная матрица с ограниченной нормой. Кроме того, в силу
(6.2) |
lim |
= 0. |
1 |
„,0 |
II» II |
В этих условиях JIIу0Ml — ограниченная величина, и
так как ограничены и все другие множители подынтеграль ного выражения, то
II У II*
ИУо II3
где с — некоторая положительная |
постоянная. |
|
|||||
Итак, |
|тр(/, у)\ <!gsup ■ |
, |
откуда |
и следует (6.8). |
|||
При условии |
ца,т] |
ПУП |
|
|
0, |
что |
|
(6.3) существует такое 6 > |
|||||||
2 |
(ехр [ 2 Re Xadi - |
1 j |
|
< - |
26 (t - („). |
||
С другой стороны, учитывая |
(6.8), можно указать такое |
||||||
ро > 0, |
что для |
всех у, |
удовлетворяющих |
неравенству |
М < Р о » будем иметь |
|я|> (/, у) | < 26, и тогда |
(см. (6.7)) |
||||
V (/> j c ) < y |
(t0, х0), а |
это означает, |
что |
любое |
решение |
|
уравнения (6.1), удовлетворяющее условию |
(6.4), |
где |
р — |
|||
произвольное положительное число из промежутка |
0 < |
|||||
<С р < ро» в |
пределах |
промежутка [/0, Т) |
удовлетворяет |
|||
неравенству |
(i(_1 (/) х, |
г-1 |
Это гарантирует вы |
|||
K~l (t)х) *< ра. |
полнение и неравенства (6.5), так как на промежутке [/0, Т)
по условию теоремы а (?) с |
ш (t). Теорема |
доказана. |
||
Т е о р е м а |
6.2. Если на промежутке Н0*^i) |
|||
и |
a (t) |
(a (f) |
|
(6*9) |
, |
Л |
г |
||
то существует |
f4*o)<0, |
г |
||
конечный промежуток |
U0, |
Т) с |/0, у , |
на котором невозмущенный процесс {тривиальное решение уравнения (6.1)) устойчив.
Д о к а з а т е л ь с т в о . При условии (6.9) по непре рывности в пределах некоторого замкнутого промежутка Но, Т] ( t o < T < у
IA( 0 < O. |
(6.10) |
Согласно неравенству (6.10) существует такое положи тельное число b, что p ,(f)< — Ь {t £ [/0, Г]).
Таким образом, на промежутке [t0, Т) условия теоремы 6.1 выполняются, и, значит, на этом промежутке невозму щенный процесс устойчив.
6.2.Обобщение теорем об устойчивости по линейному
приближению. Условия устойчивости, установленные в п. 6.1, основаны на теореме 2.2 о диагонализации линейной системы. Эта теорема определяет общий вид матрицы преоб разования линейной системы к диагональному виду. Одна ко, чтобы воспользоваться представлением (2.4), нужно располагать фундаментальной матрицей X линейной си стемы. В некоторых случаях, например в случае линейной стационарной системы, определение X, а значит и матрицы преобразования линейной системы к диагональному виду, не представляет труда. Но все же случаи, когда могут быть найдены точные выражения для X в конечном виде, исклю чительны. В то же время имеется возможность построения матрицы преобразования линейной дифференциальной си стемы к системе, «близкой» к диагональной. В связи с этим представляется целесообразным построение достаточных ус ловий устойчивости, основанных на преобразованиях тако го рода.
Допустим, что |
К (0 = (/Ci Кг |
Кп) — невырожденная |
|||
и дифференцируемая на [/0, |
Т 1 матрица, столбцы которой |
||||
имеют одинаковую норму, а именно: |
|
|
|||
и |
\\K,(t)\\ = *(t)>0 |
(/ = |
1,2, |
п) |
|
|
|
|
|
|
|
|
|
*!L = U K - K A + N > |
|
||
где А = |
diag |
Х2, ..., Ял), |
а N — некоторая |
квадратная |
|
матрица |
порядка |
п. |
|
|
|
Замена переменных
х — К у
приводит уравнение (6.1) к виду
=A(t)y — M(t)N(t)y + M(t)h(t, Ky).
Полная производная от положительно определенной функции
V(t,x) = (K~'(t)x, K~'(f)x) = l y f
по t в силу уравнения возмущенного процесса в данном случае представляется в виде
~ = 2\y{l)fv/(t, у (0) + 2 Re(y*Mh), |
(6.11) |
где
ф(*. </(0) = S ReX|,li7 iF + |
’ |
Р = —-i-(AW + /V*M*)-
Интегрируя (6.11) в пределах от /„ до 1, получаем
|
|
Ь |
|
l/(/.*) = K(!0,* 0)!l + |
|
exp j 2<p(f, g(f))dt' — 1 |
+ |
|
|
+ — *о)Ф(^, I/) |
, (6. 12) |
где |
|
|
|
|
|
f 2 § Ф (Т,у (т)) dx |
|
ф (i■9)= « - « и |
|
Re(y*/W/i)<#' |
|
л |
IP *Jо |
|