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

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

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

Докажем теперь достаточность. Пусть имеет место формула (18.5.12). Отсюда

 

K(t) = X(t)Y~l(t).

(18.5.13)

В уравнении (18.5.10) производим замену переменных

 

x = X ( t) Y ~ 1(t)y.

(18.5.14)

Получим

 

 

 

% y - ‘(t)y + X(l) * £ - у +

%

_ P(t)X(,t)Y-4l)y.

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

 

 

 

§

= Р(1)Х,

= Y -\t)B(t),

получаем

 

 

 

P ( t ) X ( t ) Y - \ t ) y - Y ~ l(t)B(t)y+ X(t)Y~\t) & = PxY~y

и, значит, dy/dt =

B{t)y. Таким

образом,

системы (18.5.10) и

(18.5.11) эквивалентны. Теорема доказана. ■ В частном случае, когда одна из систем, например, (18.5.11) яв­

ляется системой с постоянными коэффициентами dy/dt = By (В = const), мы имеем Y — eBt, и соотношение (18.5.12) принимает

вид X(t) — K(t)eB t при этом, как нетрудно видеть, теорема 18.5.3 переходит в известную теорему Н. П. Еругина о приводимости ли­ нейной системы [75].

Продолжим наши исследования по проблеме приводимости ли­ нейных систем. Рассмотрим систему с постоянными коэффициентами

4* = Dx.

(18.5.15)

где D — матрица простой структуры с чисто мнимыми собственны­

ми значениями. Фундаментальная матрица этой системы eDt, как нетрудно показать, является матрицей Ляпунова. На основе выше­ изложенного все системы, фундаментальные матрицы которых яв­ ляются матрицами Ляпунова, приводимы к системе вида (18.5.15), т.е. приводимы к системе с постоянными коэффициентами. Этот результат сформулируем в виде следующей теоремы.

Т е о р е м а 18.5.4. Все системы, фундаментальные матрицы которых суть матрицы Ляпунова, являются приводимыми сис­ темами.

18.5.3.Об устойчивости линейных систем.

Оп р ед ел ен и е . Линейную систему вида (18.5.1), фундамен­ тальная матрица которой является матрицей Ляпунова, назовем

ляпуновской системой. U

384

ПЛ. 18. ^-УСТОЙЧИВОСТЬ НА НЕОГРАНИЧЕННОМ ИНТЕРВАЛЕ

Остается рассмотреть вопрос об устойчивости системы (18.5.16) в случае, когда эта система не является ляпуновской. В этом случае вопрос об устойчивости нулевого решения линейной однородной си­ стемы решается в зависимости от того, выполняется или не выпол­ няется условие (18.5.5). Рассмотрим раздельно эти случаи.

1. Пусть система (18.5.18) не является ляпуновской, а условие (18.5.5) выполняется. Тогда система (18.5.5) неустойчива. В самом деле, если допустить устойчивость, то это в сочетании с выполне­ нием условия (18.5.5), согласно лемме 18.5.1, определяло бы, что система (18.5.16) является ляпуновской, что противоречит исход­ ной предпосылке.

2. Наконец, рассмотрим случай, когда система (18.5.16) не явля­ ется-ляпуновской и в то же время не выполняется условие (18.5.5). Учитывая, что матрица P(t) является непрерывной матрицей, невы­ полнение условия (18.5.5) означает, что имеет место следующее:

t

 

J Sp P(x)dx-+—оо при <-*•«>.

(18.5.19)

*0

 

В этом случае тривиальное решение линейной системы устойчиво*. Таким образом, может быть сформулирована следующая итоговая те­ орема об устойчивости тривиального решения системы (18.5.16).

Т е о р е м а 18.5.7. Пусть матрица P(t) системы 2 вида (18.5.16) непрерывна на [<0, »). Тогда тривиальное решение этой

системы устойчиво по Ляпунову, если Г. л е KL

или

2". 2 ё Ки но имеет место (18.5.19).

Оно неустойчиво, если 2 ё KL и выполняется условие (18.5.5).

§18.6. Устойчивость неавтономных систем

Вэтом параграфе, как и в предыдущем, в качестве заданной по­ ложительной функции a>(t) принимаем постоянную величину, т.е.

ш= const.

18.6.1.Линейная однородная система. Рассмотрим систему

§ = P ( t ) x ,

(18.6.1)

где P(t) — непрерывная и ограниченная на [а, ») л х п-матрица та­ кая, что каждое решение системы (18.6.1) ограничено на [а, »). Рас-

* Mi n - Ye n W и. On Stability of linear Time — varying Sistems//Inter. J. Syst. Sci — 1984.—N 9 2.

IIA'Cjll ’

смотрим матрицу 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 Л(т)с/т] **(т).

Матрицу B(i) в соответствии с теоремой о скелетном разложе­ нии эрмитовой матрицы представим в виде

B(t) =

(18.6.8)

где M(t) — квадратная матрица, столбцы которой h{{t) , ..., hn(t) имеют на [*0, оо) одну и ту же эрмитову норму.

Полагая ЦАу-(ОН ==%(0 (/ = 1* 2,..., л), получим

(18.6.9)

Здесь

Легко видеть, что H(t0) = K(t0) = HQ. Таким образом, рассмат­

риваемый пучок решений уравнения (18.6.1) представляется соот­ ношением

(H~l(t)x(t), H~l(t)x(t)) р2 (t0 ^ t < оо).

(18.6.10)

Используя приведенные выше построения, можно установить относительно простые условия устойчивости и асимптотической ус­ тойчивости на [f0, оо).

Т е о р е м а 18.6.1.1. Если при всех t> t0

i 2 e x p M ' ) ( r - < 0) ] < l .

(18.6.11)

0=1

 

то невозмущенное движение (тривиальное решение системы

(18.6.1)) устойчиво на [*0, оо).

Д о к а з а т е л ь с т в о . Пусть соотношение (18.6.11) выполняется при всех i > tQ. Рассмотрим рш-трубку

(G~l(i)xt G~l(t)x) < р2,

(18.6.12)

где G(t) = (оо/co2(f))#(r), a H(t) — матрица, определенная разло­ жением (18.6.8). Очевидно, G(t) G К%. Пусть х° Ф 0 — какое-ни­

будь нетривиальное решение уравнения (18.6.1), принадлежащее пучку (18.6.10). На промежутке [*0, оо) имеем

« г ‘ ( 0 * Ч 0 . е г Ч 0 * ' ( 0 ) « ^ р г-

Отсюда, поскольку со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х\ т{> 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, оо).

 

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