Достаточные условия устойчивости решений систем обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени. Часть. Iii. Нелинейные уравнения Текст научной статьи по специальности «Математика»

МАТЕМАТИКА

УДК 517.929

DOI 10.21685/2072-3040-2019-2-1

И. В. Бойков

ДОСТАТОЧНЫЕ УСЛОВИЯ УСТОЙЧИВОСТИ РЕШЕНИЙ

СИСТЕМ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ

УРАВНЕНИЙ С ЗАПАЗДЫВАНИЯМИ, ЗАВИСЯЩИМИ ОТ ВРЕМЕНИ. ЧАСТЬ. III. НЕЛИНЕЙНЫЕ УРАВНЕНИЯ

Аннотация.

Актуальность и цели. Работа посвящена анализу устойчивости в смысле Ляпунова установившихся решений систем нелинейных дифференциальных уравнений с коэффициентами и с запаздываниям, зависящими от времени. Рассматриваются случаи непрерывного и импульсного возмущения.

Материалы и методы. Исследование основано на использовании связи между устойчивостью исходных систем нелинейных дифференциальных уравнений и устойчивостью специальным образом построенных систем линейных дифференциальных уравнений. При анализе построенных таким образом систем линейных дифференциальных уравнений используются свойства логарифмических норм.

Результаты. Предложены алгоритмы, позволяющие получать достаточные условия устойчивости решений конечных систем нелинейных дифференциальных уравнений с коэффициентами и с запаздываниями, зависящими от времени. Достаточные условия представлены в виде неравенств, связывающих коэффициенты линеаризованных систем уравнений. Алгоритмы эффективны как в случае непрерывных, так и в случае импульсных возмущений.

Выводы. Предложенный метод может быть использован при исследовании нестационарных динамических систем, описываемых системами обыкновенных линейных дифференциальных уравнений с запаздываниями, зависящими от времени.

Ключевые слова: устойчивость, нелинейные системы, обыкновенные дифференциальные уравнения, запаздывания, зависящие от времени.

I. V. Boykov

SUFFICIENT CONDITIONS FOR SUSTAINABILITY OF SOLUTIONS OF SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS WITH TIME DELAY.

PART III. NONLINEAR EQUATIONS

© Бойков И. В., 2019. Данная статья доступна по условиям всемирной лицензии Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), которая дает разрешение на неограниченное использование, копирование на любые носители при условии указания авторства, источника и ссылки на лицензию Creative Commons, а также изменений, если таковые имеют место.

Abstract.

Background. The paper is devoted to the analysis of stability in the sense of Lyapunov steady-state solutions of systems of nonlinear differential equations with coefficients and with time delays. The cases of continuous and impulsive perturbations are considered.

Materials and methods. The study is based on the use of the relationship between the stability of the initial systems of nonlinear differential equations and the stability of specially constructed systems of linear differential equations. When analyzing systems of linear differential equations constructed this way, the properties of logarithmic norms are used.

Results. Algorithms are proposed that allow one to obtain sufficient conditions for the stability of solutions of finite systems of nonlinear differential equations with coefficients and with time delays. Sufficient conditions are presented in the form of inequalities connecting the coefficients of linearized systems of equations. The algorithms are effective both in the case of continuous and in the case of impulsive perturbations.

Conclusions. The proposed method can be used in the study of nonstationary dynamic systems described by systems of ordinary linear differential equations with time delays.

Keywords: stability, nonlinear systems, ordinary differential equations, time delays.

Введение

Данная статья является продолжением статей [1-3].

В работе [1] исследовалась устойчивость решений систем нелинейных неавтономных обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени

^ = fk (t, X (t ), • • •, xn (t ), xi(t - hki (t )),..., xn (t - hkn (t ))), (1)

dt

k = 1,2,...,n, t > to, с начальным множеством 9(t) = ^(t),...,фп(t)), to - H < t < to, где x(t ) = 9(t ) при to - H < t < to - непрерывная вектор-функция; x(t) = (X1(t),...,xn(t)) - искомое решение. Подобное решение обозначим через хф (t ).

В работе [1] получены достаточные условия устойчивости решения системы уравнений (1) в метрике пространства Rn векторов х = (x1,...,xn)

с нормой || x ||з= max | xk |. Достаточные условия устойчивости были

1<k <n

выражены через логарифмические нормы семейств матриц, построенных на частных производных первого порядка функций fk.

Соответствующие логарифмические нормы были ассоциированы

з

с метрикой пространства Rn. Логарифмическая норма матрицы A = {aj},

з

i, j = 1,2,.,n, в пространстве Rn с нормой || x ||з= maX | xk | определяется

1<k <n

формулой [4, 5]:

f

Л( A) = max

n

aii + Z | aj

j=1, j

При получении достаточных условий устойчивости решений систем уравнений вида (1) в работе [1] существенно использовались свойства

упомянутой выше метрики пространства Я^.

Кроме того, на функции Иуу (^) были наложены условия:

*

1) уравнения t — Иу (^) = 0, 7 , у = 1,2,..., п, имеют только одно решение t^у;

2) t > Иу (0 при t > ^.

Представляет значительный интерес получение достаточных условий устойчивости решений систем неавтономных дифференциальных уравнений с запаздываниями, зависящими от времени, справедливых в любых п -мерных пространствах, на основе методологии, свободной от использования свойств метрики конкретного пространства и от выполнения вышеупомянутых условий.

В первую очередь представляет интерес построение достаточных условий устойчивости решений систем уравнений в п -мерных пространствах

ЯП, у = 1,2,3, векторов х = (х1,...,хп) соответственно с нормами

= V I к=1

ira = ^ I хк I. 11x12 =

■ n -|1/2

tlxkl2 к=1

, X 3 = max 1 xk

1<k<n

В работе [2] получены достаточные условия устойчивости решений систем уравнений вида

= ^Оу ^)ху (t) + ^Ьу Ц)ху ^ — И«)), 7 = 1,2,..., п, (2)

У"1 у=1

в п-мерных пространствах Яу,у = 1,2,3, векторов х = (х1,...,хп). Здесь

ау ^), Ьу ^), 7,у = 1,2,...,п, - непрерывные функции; И^) - непрерывная

*

функция, удовлетворяющая следующим условиям: 0 < И^) < Н +1 при 0 < t < го. Здесь и ниже полагаем ^ =0.

В работе [3] получены достаточные условия устойчивости решений систем уравнений вида

- (,) п п

= ^ау ^)ху(0 + ^Ьу(0ху ^ — Ъц(0), 7 = 1,2,...,п, (3)

у=1 у=1

2

в п -мерном пространстве Яп векторов х = (х1,...,хп).

В данной работе результаты, полученные в статьях [1-3], распространяются на нелинейные уравнения.

Приведем определения, используемые в статье.

Через Вк g (^ М1,.., ип) обозначена частная производная Dkg ^, М1,..., ип) = = дg ^, и1,..., ип)/ дик, к = 1,2...,п.

Пусть X — банахово пространство; К — оператор, действующий из X в X; В (а, г) = {х, а е Х:| |х — а||< г}; £(а, г) = {х, а е X :| |х — а| |= г}; Л(К) — логарифмическая норма линейного оператора К, определяемая [4] выражением Л(К) = Нтн^о(| 11 + НК 11 —1) / Л, где символ Л ^ 0 означает, что Л стремится к нулю, убывая.

Для матриц в часто используемых пространствах логарифмические нормы известны.

Пусть дана вещественная матрица А = {ау}, /, у = 1,2,..., п, в п -мерном пространстве Я^, у = 1,2,3, векторов х = (х1,...,хп) с нормой

-,1/2

п

llxl1 = ZI xk I' к=1

Ixll2 =

Z|xk к=1

х 3 = max 1 хк

1<к<п

Логарифмическая норма матрицы A равна [5]:

( п \

Л1(A) = max a

j

JJ

+

Z |

ч

, Л2(A) = ^

(a AT \

A + A

Л3( A) = max

aii + Z

j=1, j *i

Здесь Xтах ((А + А ) / 2) — наибольшее собственное значение матрицы

Т

(А + А )/2; символ Лу(А(0), у = 1,2,3, означает логарифмическую норму

в прстранстве ЯП соответственно,

Предварительно приведем, следуя [1, 6], достаточные условия устойчивости систем нелинейных неавтономных дифференциальных

уравнений в метрике пространства Я^.

Рассмотрим систему нелинейных неавтономных дифференциальных уравнений

dxl (t) , , s , ss , .

= ai (t, X1(t),..., Xn (t)), i = 1,2,..., n,

dt

(4)

где функции а{(^,М1,...,ип) непрерывны по первой переменной и имеют частные производные по остальным переменным, удовлетворяющие условию Липшица с коэффициентом А :

Djai (t, х *,..., хП) - Djai (t, y*,..., y*)

<

< А(| х1 — у* | +...+ | х„ — Уп |), У = 1,..., п. (5)

Пусть система (4) имеет установившееся решение

* * *

х (^) = (х1^),..., хп (^)) при t > 0.

*

Исследуем устойчивость установившегося решения x (V) при возмущении

X (0) = X0,1 = 1,2,..., п. (6)

Справедливо следующее утверждение.

Теорема 1. Пусть система уравнений (4) имеет установившееся

* * *

решение х (t) = (Xl(t),...,хп(V)). Пусть функции щ(V,XI,...,хп), 1 = 1,2,...,п, непрерывны по первой переменной, непрерывно-дифференцируемы по остальным переменным и частные производные удовлетворяют условию Липшица (5). Пусть при всех V > 0 выполнены условия

Da(t,x*(o,...,x*(t)) + ^ °iai(t,x*(t),...,x*(t))

j=1, i &

i = 1,2,..,n, x = const.

- -x <0,

Тогда установившееся решение x (t) системы уравнений (4)

асимптотически устойчиво в метрике пространства Rn.

Для дальнейшего исследования понадобятся аналоги этой теоремы в других метриках.

По аналогии с доказательством теоремы 1 доказывается следующее утверждение.

Теорема 2. Пусть система уравнений (4) имеет установившееся

* * *

решение x (t) = (xi(t),...,xn(t)). Пусть функции at(t,xi,...,xn), i = 1,2,...,n, непрерывны по первой переменной, непрерывно-дифференцируемы по остальным переменным и их частные производные удовлетворяют условию Липшица (5). Пусть при всех t > 0 выполнены условия

Л j (A(t)) <-x <0, j = 1,2,3, x = const.

*

Тогда установившееся решение x (t) системы уравнений (4) асимптотически устойчиво в метрике пространства Rn, j = 1,2,3.

Здесь A(t) - матрица с элементами aj(t) = DJ ai(t,x* (t),...,xn(t)),

i, j = !,•••,n .

Замечание. Теорема 2 доказывается по аналогии с доказательством теоремы 1 (см. разд. 3 работы [1]. Отличие заключается в оценке

|| G(s,T,x (s),u(s)) || в различных метриках.

Приведем эти оценки.

В работе [1] показано, что устойчивость установившегося решения

*

x (t) при t > T связана с устойчивостью системы уравнений du(t)

dt

- = A(T)u(t) + G (t,T,x (t),u(t)), (7)

Известия высших учебных заведений. Поволжский регион

где А(Т) = {ау (Т)}, ау (Т) = Буа{ (Т, х*(Т),..., х**(Т)), /, у = 1,2,., п;

и ^ ) = (ul(t),..., ип ^));

* * ** \Т

G (t,T,х (t),u(t)) = (,Т,x (t),u(t)),...,gn(t,T,x (t),u(t))) . Здесь

gi(t,T, x (t),u(t)) = gi(t,T, x (t),u(t)) +

*, *

n

+Z (Djai (t, x* (T),., x*-1 (T), x* (T) + Qjuj (t)), j=1

x*+1 (T) + u j+1 (t),..., x* (T) + un (t)) -- Dai (t , x* (T),., x*-1 (T), x* (T), x*+1 (T),..., x* (T))) uj (t), i = 1,2,...,n, 0< Qj <1, j = 1,2,...,n;

4$ / 4$ 4$

gi (t,T, x (t),u(t)) = (a (t,x* (t) + u1 (t),., xn (t) + un (t)) -

4$ 4$ \ / 4c 4c

-ai (T, x* (T) + u1 (t),..., x*(T) + u (t))) - (ai (t, x (t),..., x* (t)) -

-a (T , X*(T ),..., x*(T ))), i = 1,..., n.

Требуется показать, что при t > Т | |u (t) | j < | |u (T)| j, j = 1,2,3.

В метрике пространства Rn это неравенство доказано в работе [1]. Требуется

1 2

доказать его справедливость в пространствах Rn, Rn . Уравнение (7) при t > Т имеет решение

t

u (t) = eA(T )(t -Т )u (t) + feA(T )(t-s)G*(s, Т, x*(s), u( s))ds. (8)

Т

Переходя к нормам в пространствах Rn, j = 1,2, имеем

Л (A(T))(t-Т)

| | u(t) ||j< e J | | u(T)||j +

t Л ; ( A(T))(t-S) + e '

Л

^fe j |G(s,T,x (s),u(s)) |jds, j = 1,2. (9) Т

Нужно показать, что для как угодно малого £ такого, что Л j (A(T)) + £<-%, х > 0, найдется промежуток времени t е [Т, Т + ДТ ], в течение которого

* *

|| О &Т,х (V),ы(0) ||у<е|| ы^) ||у, у = 1,2.

Докажем справедливость этого неравенства при у = 1,2. Пусть у = 1, тогда

п

|| G*(t,Т,x*(t),ы(t)) |1= ^ | й (V,Т,x*(t),ы(t)) + 1=1

+^ (Djai (T, x (T), •, x*-1 (T), x* (T) +

j=1

+Qju j (t), x*+1 (T) + uj+1 (t), •, x*(T) + un (t)) --Dat (T, x* (T), •, x*-1 (T), x* (T), x*+1 (T), •, x* (T))) uj (t) |<

< £|g,. (t, T, x*(t), u(t ))| +

i=1

££(D at(T,x (T),x*-1 (T),x*(T) + i=1j=1

+Qjuj (t), x*+1 (T) + u j+1 (t), •, x*(T) + un (t)) --D] ai (T, x* (T),•.., x*-1 (T), x* (T), x*+1 (T),•.., xl (T))) uj (t)

= I1 +12.

Из построения функций gj (V,Т, х (V), ы(t)) и предположения, что || ы(Т) |1= 5, следует, что для £1 = £ /2 найдется интервал [Т,Т + ДТ] такой, что

п

(V,Т,х*(0,ы(V)) |<£1 || ы(V) ||1.

1=1

Оценим выражение /2. Используя неравенство (5), имеем

72 = ЕЕ I Djai (T, x* (T), •, x*-1 (T), x* (T + Qj,uj (t)), x*+1 (T + uj+1 (t)), •, j=1 i=1

x* (T + un (t)) - DJai (T, x*(T),•.., x*-1 (T), x* (T), x*i+1 (T),•.., x* (T)) || uj (t)|<

< EAnduj(t)|+- +|un(t) |) |uj(t)|< j=1

< AnEE^i(t) |)| uj(t)|< An ||u(t)| j=1 i=1

Отсюда следует, что существует промежуток времени [Т ,Т + Л2Т ], Л2Т < Л{Г, такой, что | | и(0 | 1< (1 + £2) 11 и(Т) 11= (1 + £2)5. Следовательно, при t е [Т,Т + Л2Т]

12 < Ап(1 + £2)5 11 и^)| | 1. Полагая 5 таким, что Ап(1 + £2)5< £1 = х /2, имеем

12 < £1 11 и^)| 1 = х 11 иЦ)| | 1/2. Подставляя эти значения в неравенство (9), имеем

,, Л1( А(Т ))(t -Т К ._,.„ 1 \ Л1( А(Т ))(t-5)„ . ... ,

| | и(0 | 1 < е ^ '| | и(Т)| 1 + -х^ у | и(5)| 1 Ж.

Т

Отсюда стандартными рассуждениями убеждаемся в справедливости теоремы 2 в метрике пространства Я.

Пусть у = 2. Применяя неравенство Минковского, имеем

Z

i=1

I2 <|| G(t,Т,x (t),u(t))| Ь +

+Qjuj (t), x*+1 (Т) + uj+1 (t),..., xn (Т) + un (t)) -

Z (Da, (Т, x* (Т),., x*-1 (Т), xj (Т) + j=1

- D] a (Т, x* (Т),., x*-1 (Т), x* (Т), x*+1 (Т),..., xn (Т))) uj (t)

1/2

=121 +I22.

Нетрудно видеть, что для любого как угодно малого £(£ > 0) найдется такой промежуток времени [Т, Т + ЛТ ], в течение которого

121 <£ 11 и ^ )| | 2 /2. Приступим к оценке /22. Очевидно,

п п

2

122 < ZI Z (DJai (Т, x* (Т),., x*-1 (Т), xj (Т) + Q.u. (t), xj+1 (Т) + u.+1 (t),., i=1 j=1

xj (Т) + un (t)) - DJai (Т, xj (Т),., x*-1 (Т), x* (Т), x*+1 (Т),., xj (Т))) X n n n

xuj(t)|2 Zlu;(t)|2 = ||u(t)| 2 ZZ (DJai(Т,x*(T),..„x*-1(T),x*(T) + j=1 j=1 i=1

+Qjuj (t), x*+1 (Т) + uj+1 (t),..., x*(T ) + un (t)) -

- DJ ai (Т, x* (Т),..., xj-1 (Т), xj (Т), xj+1 (Т),..., x* (Т))

<

< | \и ^) | 2 ^ )||2 < а2 | ) 12 п31 \и (о| 2.

/=1 ¿=1

Так как по предложению | |и (Т)| 2= 6, то существует такой интервал [Т, Т + ЛТ], А1Т < АТ , в котором

| |и^ )| 2 < (1 + £1) | |и (Т )| | 2= (1 + £1)6.

Так как 6 определяет норму начального возмущения, то можно

выбрать 6 таким образом, чтобы Ал/п3(1 + £1)6< х/2. Тогда

722 < -2 х | |и(()| 2,

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

h ^2X I\u(t)1 2-

В результате имеем при t е [Т, Т + ЛТ] неравенство

II /чи.- Л2(А(Т))(1-Т) ,, 1 Г Л2(А(Т))(/-5) ,, /чи ,

| |и(0 | 2 < е 2 у | и (Т) | 2 +-х\в 2 '| |и(у) | 2 Л, (10)

2

Т

из которого стандартным способом завершается доказательство теоремы.

2. Устойчивость решений систем нелинейных дифференциальных уравнений с постоянными запаздываниями, зависящими от времени

Рассмотрим систему нелинейных неавтономных дифференциальных уравнений с запаздываниями

= а1 С, *1('),..., хп ( t)) + Ь (^ Х ((- Н(()),..., Хп (í - Н(())), / = 1,2,...,п, (11)

где функции а(¿,Х1,...,хп),Ь{Х1,...,хп),7 = 1,2,...,п, непрерывны по первой переменной и имеют частные производные по остальным переменным,

удовлетворяющим условию Липшица (5); И^) - непрерывные функции,

*

удовлетворяющие следующим условиям: 0 < И* < И(¿) < Н +( при 0 < t <

*

Пусть при — Н <t< 0 заданы начальные условия

Х(0 = П( t), 7 = 1,2,...,п, (12)

где П7(t), 7 = 1,2,...,п, - непрерывные функции.

Будем считать, что задача (11), (12) имеет установившееся решение

* * *

х (0 = (х^),...,хп (0), которое в момент времени ^ =0 получает возмущение Х0 =(х0,..., х^). Тогда

XСо) = + 4,' = 1,2,...,п. (13)

Исследуем устойчивость решения системы уравнений (11) при условиях (12) и (13). Вначале рассмотрим случай, когда на функцию И^) налагается дополнительное условие И^) ^ t при всех значениях ¿,0< t < за

исключением одной точки, которую обозначим как t .

Здесь нужно рассмотреть два случая:

* *

1) И^)> t при t < t и И^)< t при t > t . Как и в работе [1], это

условие будем называть условием А ;

* *

2) И^) < t при t < t и И(0 > t при t > t . Это условие будем называть условием В.

Выше было положено ^=0. Введем функции п, ^), определяемые формулами

Xi (t) = Xj (t) + щ (t), i = 1,2,..., n. (14)

В результате такой замены система уравнений (11) принимает вид

dui (t) * * * */чч

—— = a (t, Xi (t) + Ui (t),..., Xn (t) + Un (t)) - aj (t, Xi (t),..., Xn (t)) +

dt

+bi (t,X*(t - h(t)) + u1(t - h(t)),..., X* (t - h(t)) + un (t - h(t))) -

-bj (t,X*(t -h(t)),...,Xn(t -h(t))), i = 1,2,...,n; (15)

условия (13) принимают вид

Uj (0) = Xo, i = 1,2,..., n, (16)

а условия (12) можно интерпретировать как

щ(t) = 0, — < t <0, i = 1,2,...,n. (17)

Обозначим через 5j достаточно малое положительное число, величина

которого определена ниже для каждого пространства R^ , и предположим, что | | u (0) 11j <5 j, u (0) = (u1(0),..., un (0)), j = 1,2,3. Ниже для простоты обозначений индекс j в 5j будем опускать.

Отметим также, что в случаях, когда проводится рассуждение общее для всех пространств Rjl, j = 1,2,3, индексы j опускаются.

Найдем условия, при которых траектория решения задачи Коши (15)-(17) не покидает шар 5(0,5), где 51 <%/2An, 52 <%/2An3/2,53 <%/4An2. Величина х будет определена ниже.

Вначале рассмотрим первый случай.

*

Пусть t е [0,t ). В этом промежутке времени система уравнений (15) имеет вид

du; (t) * *

—— = ai (t, x (t) + uj(t),..., xn (t) + un (t)) -dt

-a;(t,x*(t),...,x*(t)), i = 1,2,...,n. (18)

Устойчивость решения задачи Коши (18), (16) исследуется методами, изложенными в введении данной статьи.

*

Из теоремы 2 следует, что если при t е [0, t )

Лj(A(t)) < -х, х > 0, 0 < t < t*, j = 1,2,3, (19)

X = const, то траектория решения задачи Коши (15)-(17) не покидает шар

X X X B(0,5 ), 61 <-, 62 <-ттт, 5з <-т-, в течение этого промежутка

J 2 An 2 An3/2 4 An2

времени.

Здесь A(t) = {aj(t)}, i, j = 1,2,...,n, aj (t) = DJa; (t, x1 (t),...,xn (t)),

i, j = 1,2,..., n.

Из непрерывности решения x(t ) = ( x1(t ),..., xn (t)) следует, что

траектория решения задачи Коши (15)—(17) не покидает шар B(0,5)

*

в течение промежутка времени [0, t ].

*

При t > t необходимо учесть влияние слагаемых, содержащих задержки.

*

В промежутке времени t е [t , «) система уравнений (15) имеет вид

dui(t) * * *

—— = at (t, x1 (t) + u1 (t),..., xn (t) + un (t)) - ai (t, x1 (t),..., x* (t)) +

dt

+Ь (t, x* (t - h(t)) + u1(t - h(t)), x2 (t - h(t)) + u2(t - h(t)),..., x* (t - h(t)) +

+u* (t - h(t))) - Ьi (t, x* (t - h(t)), x2 (t),..., x* (t)), i = 1,2,..., n. (20)

*

Покажем, что в случае выполнения при t е [t , «) условий

Лj(A(t))+ ||B(t) |j<-х, tе [t*,«О, j = 1,2,3, (21)

X = const > 0, траектория решения системы уравнений (20) при начальном значении

u (t) = (u1(t),..., un (t)), t е [t* - H*, t*], (22)

*

не покидает шар B(0,5). Напомним, что || u(t)||j<5,0 <t <t , j = 1,2,3. Здесь

B(t) = {bjj (t)}, i, j = 1,2,..., n, bjj = Djb (t,x*(t -h(t)),...,x*(t -h(t))), i, j = 1,2,...,n.

Доказательство проведем от противного. Пусть в момент времени 7j,

*

t < 71 < траектория решения задачи Коши (20), (22) покидает шар B(0,5). Представим при t > 7 систему уравнений (20) в следующем виде:

du (t) n ^ ^

-d^- = YP3^i (71, X1T1),..., x* (T1 ))u j (t) +

j=1

n

+Y,DJ'bi (T1,X*(t-h(t)),x2 (t - h(t)),..., x*(t - h(t))) (t- h(t)) +

j=1

+gi (, u (t), u (t - h(t))), i = 1,2,3,..., n. (23)

Структура функций gj(t,u(t),u(t -h(t))), i = 1,2,3,...,n, очевидна. Систему уравнений (23) можно представить в операторном виде

ddP = A(T1 )u (t) + B(T1 )u (t - h(t)) + G (t, u (t), u (t - h(t))), (24)

где

A(T1) = a (T1)}, B(T1) = {bj (T1)},i, j = 1,2,..., n, G(t,T,u(t),u(t -h(t))) = (g1(t,T,u(t),u(t -h(t)),...,gn(t,T,u(t)),u(t -h(t)))T .

Построение вектор-функции G(t,T,u(t),u(t - h(t))) очевидно. Уравнение (24) при t > T1 имеет решение

t

u(t) = /(7lXt-71 )u(T1) + ^eA{T1){t-s)B(T1)u(s - h(s))ds +

T1

t

+

JeA(T1)(t s)G(s,T,u(s), u(s -h(s)))ds. (25)

T1

Переходя к нормам, имеем

Л . (A(T1))(t-T1) I u(t) ||j< e j | |u(T1)| j +

P Л.(A(T1))(t-s)/ii „ „ „ >

+Je jV (B(T1)u(s - h(s))| |j + | |G(s,T1, u, h(s)) | j )ds <

T1

< Л j (A(T1))(t-T1) (t \л (A(Tx))(t-s) h

< e J 1 1 | |u(T1)| j + Je J 1 | |B(T1)u(s -h(s)) 11j ds +

T1

+ e

Л (A(Tx))(t-s)

G(s,Ti,u,h(s)) |1 j ds =

= (0 + ^ (О + у (0, ] = 1,2,3. (26)

Оценим каждое слагаемое в (26) в отдельности. Оценка У ] ( 0 очевидна:

3]( 0 = Л'(А(71Ш_71) 11„(31) 11], ] = 1,2,3. Для оценки ^2 ] ( 0 необходимо предварительно оценить норму

| | 5(7])и(5 - Н(5)) |] в каждом из пространств , ] = 1,2.

3

Замечание. В пространстве слагаемое 32 (^) не оценивается, так как соответствующая оценка была получена в [1].

Так как рассуждения, проводимые в пространствах Я], ] = 1,2,3, аналогичны, то в тех местах, где это возможно, будем проводить общие рассуждения, отмечая принадлежность оценок различным пространствам нижними индексами. Например: , ^2г, ^31, ' = 1,2,3. Очевидно,

| B(Ti)u(t -h(t)) | 1 < ££ Djbi (Ti,x* (t -h(t)),...,xn (t -h(t)))uj(t -h(t))

i=1 j=1

ri n

< £ |uj (t - h(t))| £ Djbi (t, x* (t - h(t)),., Xn (t - h(t))

j=1

i=1

< max £ Djb (Ti,x*(t -h(t)),.,Xn(t -h(t))

1<j<ni=1 v

=| | B(Ti)| 1 | |u(t - h(t ))| 1;

|u(t - h(t ))| 1 =

I B(Ti )u(t - h(t)) | g = £ | £Djbi (Ti, x* (t - h(t)),., x*(t - h(t)))u j (t - h(t)) | | 2<

i=1 j=1

< ££iDjbi (Ti, x*(t - h(t)),., xn (t - h(t )))|2 £ i uj (t - h(t ))|2= i=1 j=1 j=1

= £ | uj (t - h(t))|2 ££ Djbi (, x*(t - h(t)),..., x*(t - h(t)) j=1 i=1 j=1

=| | B(T1) | |2 11 u(t - h(t))| |2.

Известия высших учебных заведений. Поволжский регион В работе [1] было показано, что

I | В(Т)и^ -И(0) | |з< 11В(?!) | |з11 и^ -И(t)) I 3.

Оценим | |В(Т1)и(t-И(t)) ||у в метрике пространств ЯП, у = 1,2.

* *

Здесь нужно рассмотреть два случая: 1) Т = t , 2) Т > t .

Вначале рассмотрим второй случай.

*

Так как И^) < t при t > t , то t - И(t) > 0 при t > Т1. Как отмечалось выше, | |и (t)| |у < 5, у = 1,2 при t е (0, Т[). Следовательно, существует промежуток времени [7],Т[ + Д2Т1] такой, что | | и^ -И(Г» ||у < 5, у = 1,2 при tе[?1,Т1 +Д2Т1].

Так как по предположению | | и(Т1) | |у = 5, у = 1,2, то существует промежуток времени [Т1, Т +Д3Т1] такой, что

| | иЦ) ||у > (1 - а) 11 и(Т1) | | = (1 -а)5, у = 1,2,

где а >0 - вещественное число.

Следовательно, в промежутке времени [71, Т+Д4Т1], где Д4Т1 =ш1и(Д1Т1,..., Д3Т1) справедливо неравенство

| | и(5 -И(5)) ||у< 11и(Т1) ||у< 11 и(5) ||у /(1 -а), у = 1,2,

и имеет место оценка

1 К Л ■ (А^))^-5)

^2 у (t) <| | В(Т1) | |у — ¡е у | |и(5)| уж, у = 1,2. (27)

Нетрудно видеть, что а можно выбрать таким образом, что из (27) вытекает оценка

К Л . (А(Т1)). (t-5)

^2у(t)<(| |В(Т1) ||у +е)¡е | u(s)|^Ж, у = 1,2. (28)

Т1

Приступим к оценке Jз.

Так как ау (0, Ьу (,u(t),и^ -И(t))), I,у = 1,2,.,п, - непрерывные функции, то найдется такой промежуток времени [ТьТ +Д5Т1], в котором

| |С(5,Т1,и,И(5))|у<е11 и(5) ||у, у = 1,2. (29)

Из (29) следует, что

t

J

3j

,Л . (A(T1 ))(t—s)

(t) <e Je ^ ^ | u (s )| Jds, J = 1,2. (30)

T1

Из неравенств (26)-(30) следует, что

Л ] (А(Шг-71)

j t

+

u (t)| | j < e j 1 1 | |u (T1)| j +

Л . (A(Tj))(t-s)

В(71) |] +2е)|е ] 1 | | и(5)|] Ж, ] = 1,2,3. Т1

Отсюда имеем

| | и(О |]< ехр{(Л](А(Т{))+ 11 В(Т1) |] +2е) -7^} || «(71) |] Ж, ] = 1,2,3. (31)

* *

Рассмотрим теперь случай, когда 71 = ^ . Так как 7] = ^ , то и(Т1 - Н(Т1)) = и (^ - Н(**)) = «(0) и | | и(Т1 - Н(Т1)) 11 ] = 11 и(0) |] = 5 = 11 и(Т1) 11 ], ] = 1,2,3. Следовательно, существует промежуток времени [Т^Т +ДТ1], в течение которого

(1 -е1) 11 и (0) | ] < 11 и (5) | ] < (1 + е1) 11 и(0) | ] . Тогда при tе[71, Т1 +ДТ1] имеем

32] (t) <|| В(Т1) 1] \е ] | | и(5 - Н(5)) |] Л <

Т1

< 1 + £1 | |В(Т )| | Г Л] (^»(^ | ( )| | л ■ 1 2 3

<--11|В(Т1) 11 ] 1е ^ | |и(5) 11 .-Ж, ] = 1,2,3,

1 Т1

где £1 - как угодно малое положительное число. Для простоты обозначений положим

||В(Т1) 1] = | | В(Т1) 1] +£, ] = 1,2,3.

1 -£1 Следовательно,

t

J2

, чГЛ ■ (А(Т))^-5)

2](t)<((В(Т1)||] +е)|е ] 1 | |u(s)|] = 1,2,3. (32) Т1

Аналогичным образом доказывается справедливость неравенства (30)

*

при Т1 = t .

Из неравенств (30), (32) следует оценка

| | и(0 |]< ехр{(Л](А(Т1))+ 11 В(Т1) |] +e)(t -Т1)} 11 и(Т{) |], ] = 1,2,3. (33)

*

Из (31), (33) следует, что эта оценка справедлива при Т е ^ ,

Таким образом показано, что при выполнении условий (21) траектория решения задачи (15)-(17) не покидает шар В(0,5) в течение промежутка времени [0,

*

Рассмотрим теперь случай, когда И(у)< t при tе [0,t ). В этом случае

*

при t е [0, t ) возмущение задачи Коши (11), (13) описывается задачей Коши

(15), (16). Исследование устойчивости тривиального решения проводится по

*

аналогии с исследованием первого случая при t > t . Таким образом, при выполнении неравенства

Л у (А(0)+ 11 В(0 | у < -х < 0,t е [0, М, у = 1,2,3,

траектория решения задачи Коши (15), (16) не покидает шар В(0,5).

*

Рассмотрим теперь промежуток времени ^ , В этом промежутке

И^)> t и, следовательно, устойчивость решения задачи Коши (15), (16)

сводится к исследованию системы (18). Это исследование было проведено выше.

Из проведенного выше анализа следует, что при выполнении неравенства

Л у (А(0)+| |В(0 | у <-х < 0, [t *, ~), у = 1,2,3, (34)

*

траектория решения задачи Коши (11), (12) не покидает шар В(х ,5). Остановимся на вопросе об асимптотической устойчивости. Выше было показано, что вопрос об устойчивости решения задачи Коши (11), (12) сводится к исследованию устойчивости системы уравнений (23).

Объединяя проведенное выше исследование и выкладки, проведенные в работах [1-3], можно показать, что при выполнении условий (34) справедливо неравенство

| | и(01у< ехр{(Лу(А(0)+ 11 В(01у) / 2} 11 и(0) |у, у = 1,2,3.

Таким образом, доказано следующее утверждение. Теорема 3. Пусть выполнены следующие условия:

*

1) задача Коши (11), (12) имеет установившееся решение х ^);

2) функции а^ ^, у,..., гП ), Ь ^, г^,..., гП ) непрерывны в шаре

В(х>), г), г >0;

3) функция И(у) непрерывна при t е [0,

4) выполнено неравенство (34) при у = 1,2,3.

*

Тогда установившееся решение х (^ асимптотически устойчиво

в метрике пространства ЯП, у = 1,2,3.

Из этой теоремы непосредственно следует утверждение об

*

асимптотической устойчивости установившегося решения х ^) задачи Коши (11), (12) в предположении, что уравнение И^) = t имеет конечное число решений при t е [0,

В случае, если уравнение h(t) = t имеет счетное множество различных корней, то для доказательства асимптотической устойчивости решения задачи Коши (11)—(12) приходится ввести дополнительное условие.

Теорема 4. Пусть выполнены условия 1-4 теоремы 3 и пусть, кроме того, выполнены условия:

5) уравнение h(t) = t имеет счетное множество корней Z1,...,Zn,•••;

Z1< Z2< - <Z n < •••;

6) справедливы неравенства Zk+1 - Zk - H >0.

Тогда установившееся решение задачи Коши (11)—(12) асимптотически устойчиво.

Доказательство является объединением доказательств предыдущей теоремы и теоремы 3 работы [3].

Замечание. Аналогичные утверждения имеют место и в случае возмущения вектор-функции n(t) = (П1 (t ),•••, Пп (t)) вектор-функцией П (t ) = (П1(t),..., П п (t)).

Замечание. Случай кратных корней не рассматривается, так как все рассуждения проводятся по аналогии с предыдущими.

Заключение

В работе получены достаточные условия устойчивости решений систем нелинейных обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени. Условия представлены в виде систем неравенств, связывающих коэффициенты исходных уравнений.

Библиографический список

1. Бойков, И. В. Устойчивость установившихся решений систем нелинейных неавтономных дифференциальных уравнений с запаздываниями / И. В. Бойков // Дифференциальные уравнения. — 2018. — Т. 54, № 4. — С. 435—457.

2. Бойков, И. В. Достаточные условия устойчивости систем обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени. Часть I. Линейные уравнения / И. В. Бойков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. — 2018. — № 4 (48). — С. 3—19.

3. Бойков, И. В. Достаточные условия устойчивости решений систем обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени. Часть II. Линейные уравнения / И. В. Бойков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. — 2019. — № 1 (49). — С. 63—77.

4. Далецкий, Ю. Л. Устойчивость решений дифференциальных уравнений в банаховом пространстве / Ю. Л. Далецкий, М. Г. Крейн. — Москва : Наука, 1970. — 536 с.

5. Деккер, К. Устойчивость методов Рунге — Кутты для жестких нелинейных дифференциальных уравнений / К. Деккер, Я. Вервер. — Москва : Мир, 1988. — 335 с.

6. Бойков, И. В. Устойчивость решений дифференциальных уравнений : монография / И. В. Бойков. — Пенза : Изд-во ПГУ, 2008. — 244 с.

References

1. Boykov I. V. Differentsial'nye uravneniya [Differential equations]. 2018, vol. 54, no. 4, pp. 435—457. [In Russian]

2. Boykov I. V. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2018, no. 4 (48), pp. 3-19. [In Russian]

3. Boykov I. V. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2019, no. 1 (49), pp. 63-77. [In Russian]

4. Daletskiy Yu. L., Kreyn M. G. Ustoychivost' resheniy differentsial'nykh uravneniy v banakhovom prostranstve [Stability of solutions of differential equations in the Banach space]. Moscow: Nauka, 1970, 536 p. [In Russian]

5. Dekker K., Verver Ya. Ustoychivost' metodov Runge - Kutty dlya zhestkikh nelineynykh differentsial'nykh uravneniy [Stability of Runge-Kutta methods for stiff nonlinear differential equations]. Moscow: Mir, 1988, 335 p. [In Russian]

6. Boykov I. V. Ustoychivost' resheniy differentsial'nykh uravneniy: monografiya [Stability of solutions of differential equations: monograph]. Penza: Izd-vo PGU, 2008, 244 p. [In Russian]

Бойков Илья Владимирович

доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

Boykov Il'ya Vladimirovich Doctor of physical and mathematical sciences, professor, head of sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)

E-mail: boikov@pnzgu.ru

Образец цитирования:

Бойков, И. В. Достаточные условия устойчивости решений систем обыкновенных дифференциальных уравнений с запаздываниями, зависящими от времени. Часть. III. Нелинейные уравнения / И. В. Бойков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2019. - № 2 (50). - С. 3-20. - БО! 10.21685/2072-3040-2019-2-1.