Онлайн-доступ к журналу: http: / / mathizv.isu.ru
Серия «Математика»
2018. Т. 25. С. 159-169
УДК 517.929.4 МЭС 34К20
Б01 https://doi.org/10.26516/1997-7670.2018.25.159
Об устойчивости решений дифференциальных уравнений нейтрального типа с распределенным запаздыванием*
Т. К. Ыскак
Институт математики им. С. Л. Соболева СО РАН, Новосибирский государственный университет, Новосибирск, Российская Федерация
Аннотация. Рассматривается один класс систем линейных неавтономных дифференциальных уравнений нейтрального типа с распределенным запаздыванием. Матрица перед производной неизвестной вектор-функции с запаздыванием постоянна, матрица перед неизвестной вектор-функцией имеет непрерывные Т-периодические элементы, ядро интегрального оператора состоит из непрерывных функций, Т-пе-риодических по аргументу Цель работы заключается в исследовании асимптотической устойчивости нулевого решения с использованием метода модифицированных функционалов Ляпунова - Красовского. Метод функционалов Ляпунова -Красовского является развитием второго метода Ляпунова. Достоинством этого метода является простота формулировок и сведение исследования асимптотической устойчивости к решению хорошо обусловленных задач. Кроме того, метод модифицированных функционалов Ляпунова - Красовского позволяет получить оценки на решения линейных систем дифференциальных уравнений с запаздывающим аргументом. Заметим, что использование модифицированных функционалов Ляпунова -Красовского также позволяет получить оценки решений нелинейных дифференциальных уравнений и оценки на множество притяжения. Ранее система периодических дифференциальных уравнений нейтрального типа рассматривалась в работах Г. В. Демиденко и И. И. Матвеевой, в которых были получены достаточные условия асимптотической устойчивости нулевого решения и указаны оценки решений этой системы. Система линейных периодических дифференциальных уравнений с распределенным запаздыванием рассматривалось автором статьи. Для этой системы также были получены достаточные условия асимптотической устойчивости нулевого решения и указаны оценки решений. В настоящей работе получены достаточные условия асимптотической устойчивости нулевого решения системы нейтрального типа с распределенным запаздыванием в терминах матричных неравенств и уста,-
* Работа выполнена при финансовой поддержке РФФИ, грант 18-31-00408.
новлены оценки решений системы, характеризующие экспоненциальное убывание на бесконечности.
Ключевые слова: асимптотическая устойчивость, функционал Ляпунова - Кра-совского, распределенное запаздывание, уравнение нейтрального типа.
1. Введение
Начало теории функционально-дифференциальных уравнений было положено в 50-х годах прошлого века в работах [5; 6]. Вопросам устойчивости функционально-дифференциальных уравнений посвящено множество монографий и статей, отметим лишь несколько из них [1;2; 4; 5; 7; 8; 10].
В данной работе рассматривается следующая система нейтрального типа с распределенным запаздыванием
t
jt(y(t)+Dy(t-t)) =A(t)y(t)+j B(t,t-s)y(s)ds, t> 0, (1.1)
t-r
где r > 0 — запаздывание, D — квадратная матрица порядка n, A(t) — квадратная матрица порядка п с непрерывными Т-периодическими элементами, B(t,£) — квадратная матрица порядка п с непрерывными элементами, Т-периодическими по первой переменной, т. е. B(t,£) = B(t + T,o.
Цель работы заключается в исследовании асимптотической устойчивости нулевого решения и получении оценки решений исследуемой системы, которая характеризует скорость убывания при t —> оо. При исследовании асимптотической устойчивости нулевого решения будем использовать следующую модификацию функционала Ляпунова - Кра-совского, введенную в [3; 9]
V(t,y) = (H(t) (y(t) + Dy(t — г)), (y(t) +Dy(t-r)) )
ri t
+ J J {K(t-s)y(s),y(s))dsdr]+J {M(t-s)y(s),y(s))ds. (1.2)
0 t-7] t-T
В работах [3 ; 9 ; 11] исследован случай дифференциального уравнения нейтрального типа с запаздыванием. В [12] рассмотрена система (1.1) в случае D = 0.
Автор выражает глубокую благодарность д.ф.-м.н. Г. В. Демиденко, к.ф.-м.н. И. И. Матвеевой и к.ф.-м.н. М. А. Скворцовой за внимание и ценные советы.
ОБ УСТОЙЧИВОСТИ РЕШЕНИЙ УРАВНЕНИЙ НЕЙТРАЛЬНОГО ТИПА 161 2. Основная часть
Рассмотрим начальную задачу для системы (1.1) при t > 0:
t
jt (y(t) + Dy(t - т)) = A(t)y(f) + j B(t, t - s)y(s)ds,
t-r (2.1)
y(t)=<p(t), t € [-т,0], I y(+0) = p(0).
Теорема 1. Пусть существуют гладкая Т-периодическая матрица H(t) = H*(t) такая, что H(t) > 0 и матрицы K(s) = K*(s), M(s) = M*(s) € С1 [0, г] такие, что
K(s) > 0, ^j-K(s) < 0, M{s) > 0, ^-M(s) <0, s € [0,т].
(Л/О (Л/О
Обозначим через q(t) минимальное собственное число матрицы Q(t) = H"2(t) (Qn(t) - Qi2(t)Q^(t)Q*12(t)
t-T
QiÁt, s)Qsi(t, s)Q*13(t, s)ds^j H~l2(t),
где
Qn(t) = - (jtH(t) + H(t)A(t) + A*(t)H(t) + M(0) + tK{0) ) ,
Qu(t) = -^=(H(t)A(t) + M(0))D + V^K(0)D, v r
Q13(t,s) = —H(t)B(t, t — s), Q22 = —(M(r) - D*M(0)D) - D*K(0)D > 0,
T
Q33(t,s)=K(t-s). Выберем число k > 0 такое, что
^-K(s) +kK(s) < 0, M(s) + kM(s) <0, s € [0,т].
Uib Cvb
Тогда для решения задачи (2.1) верна следующая оценка
V(t,y) < exp I - J j(s)ds
х (<Я(0)(р(0) + D<p(-t)), МО) + D<p(-t)))
т О О
+ JJ(К(-з)ф),ф)) dsdV + J(М(-з)ф),ф)) , (2.2)
О где
7(t) = min {q(t),k}
Доказательство. Пусть у it) — решение начальной задачи (2.1). Рассмотрим модифицированный функционал Ляпунова - Красовского (1.2) на решении y(t). Его производная будет иметь вид
= l^jH{t)+H{t)A{t) + A*{t)H{t)+rK^)+M{^ y(t),y(t) + (y(t)t ( jH{t)D + A*H(t)Dj y(t -r)^
+ (y(t),H(t) J B(t, t — s)y(s)ds^J
+ ( (jH(t)D + A*H(t)Dj y(t - r),y(t + ( (D*JtH(t)D ~ MW) V(t ~ 'h V(f ~ r)
+ (v(t - t), D*H(t) J B(t, t - s)y(s)ds^J + /tf(i) J B(t,t-s)y(s)ds,y(t)\
t-T t
+ (D*H{t) j B(t,t-s)y(s)ds,y(t-r)\
\ t-T '
t T t
+ J (K(t-s)y(s),y(s))ds + J J (J^K(t-s)y(s),y(s)^ dsdri
t-T 0 t-17
ОБ УСТОЙЧИВОСТИ РЕШЕНИИ УРАВНЕНИИ НЕЙТРАЛЬНОГО ТИПА t
+ J (ftM(t-s)y(s),y(s))ds.
t-T
Перепишем данное тождество в терминах квадратичной формы 4 ' ^(t) \ ( y(t)
^-V(t,y) + í (c(t,s)(Vy(t-T)),(Vy(t-r))\ds dt ¿\ \y(s) ) \y{s) I/
t t
0 t—T¡ t
+
t-T
t
где
Cn(t) C12(t) C13(t,s) C(t,s) = — I C*12(t) C22(t) C23(t,s) C*13(t,s) eut, s) C33(t,s)
Cn(í) = i (^-H(t) + H(t)A(t) + A*(t)H(t) + Ttf(O) + M(0) t \at
Cu{t) = l{jtH{t)D + A*H^D) ' Ci3(M) = H{t)B{t,t-s),
C22{t) = -(D*^-H{t)D-M{r) t \ at
C23(t,s) = D*H(t)B(t,t-s),
C33(t,s) = -K{t-s).
Нетрудно убедиться в справедливости следующего тождества:
/ / 21 \ / \ \ / ( zi+ Dz2\ f zi+Dz2
MSKSJH0(m)U j'U
где
/ On(t) Ol2{t) 0\3(t, s)
0(M)= o*l2(t) o22 0
\0*l3(t,s) 0 033(t,s) Ou(i) = -i ( jH(t) + H{t)A{t) + A*(í)tf (í) + M(0)) - K{0),
On(t) = -(H(t)A(t) + M(0))D + K(0)D, т
013(t,s) = —H(t)B(t, t — s),
022 = —(M(r) - D*M(0)D) - D*K(0)D, т
033(t,s) = K(t-s).
Действительно,
/ / / -D 0 \ / zi + DZ2 \ ( I -D 0 \ / Zl + DZ2
= {C(t,s) 0 / 0 \\ z2 , 0 / 0 \\ z2
\ \0 0 IJ \z3 J \0 0 IJ \z3
Следовательно,
(c\t,s) ^2 j , ^2
/ 0 0 \ / / -D 0 \ / Zl + Dz2 \ ( Z1 + DZ2
-D* I 0 C(t,s) 0 / 0 \\z2 , I z2
0 0 IJ \0 0 IJ \z3 J \z 3
Перемножив матрицы, получим указанное тождество.
В силу того, что О22, Озз(ij s) > 0, имеет место следующее представление
On(t) o12(t) o13(t,s)\ / zA ( zx
out) 022 0 \[z2\,[z2
0*13(t,s) 0 Озз(i, s) J \z3J \z3
= ((On(i) - o12(t)o^o*12(t) - 013^,8)0^1^8)0*13^8))г1,г1)
+(022\Q22Z2 + 0*l2(t)zг), (022z2 + 0*l2(t)zA +(033\t, s)(033(t, s)z3 + OUt, s)z 1), (033(t, s)z3 + 0*13t8(t, s)zA-Следовательно, из (2.3) имеем
d_
dt
I
V(t,y)<- J {[0n(t)-012(t)0^(t)0Ut)
t-T
-013(t, s)0^(t, s)OUt, s)](y(t) + Dy(t - t)), (y(t) + Dy(t - r)))ds
т t t
+
0 t-r] t-т
T t t
/ / (jtK{t - dsdrl+ / (ftM(t-s)y(s),y(s)^ ds,
ИЛИ
fv(t, у) < —(Q(t)H^(t) (y(t) + Dy(t - т)), m{t){y{t) + Dy(t - т))>
т t
+
О t-rj t-т
т t t
J J (jtK(t-s)y(s),y(s)^ dsdr] + J ^M(t-s)y(s),y(s)^ ds.
Из определения <?(£), к и 7^) следует справедливость следующего неравенства
Домножив обе части неравенства на ехр ( J 7(s)ds) и проинтегрировав,
^ о '
получим оценку (2.2).
□
Замечание 1. Нетрудно увидеть, что функция 7(t) является Т-перио-дической.
Отметим, что выполнение условий теоремы 1 (точнее условий на матрицы Q221 К (s), M (s)) влечет принадлежность спектра матрицы D единичному кругу. Согласно критерию Ляпунова спектр матрицы D лежит в единичном круге тогда и только тогда, когда существует решение L в классе положительно определенных эрмитовых матриц дискретного уравнения Ляпунова
L - D*LD = С, С = С* > 0.
Введем обозначения
к = шах \\Н~\0\\1/2
1
т V г \ 2
х | 2||Я(0)||(1 + Ц-СЦ2) + 11 \\K(s)\\dsdr] + I \\Mmd8
0 0 о
/ Т ?
а = ехр I шах ( 4 [ - [
1 ?€[о,т] I Т .} 2 У 2 ^ \ о о
В следующей теореме мы установим оценки решения начальной задачи (2.1), являющиеся аналогами оценок из [3;9;11].
Теорема 2. Пусть
а) выполнены условия теоремы 1,
б) I — такое минимальное число, что ||£>г|| < 1,
т
«) л = I ^>0.
о
1. Если \\01\\ ехр(Д/т) < 1, то для решения задачи (2.1) справедлива оценка
ШН <
V 3=0
+ тах{||_0||еЛг,..., ) шах
J ге[-т,о]
2. Если ||-Сг|| ехр(Д/т) = 1, то для решения задачи (2.1) справедлива оценка
( 1
V з=0
+ тах{||_0||еЛг,..., ||-Ог_1||е(-г_1-)Лг, 1} | шах ^(¿)|е"А4.
) *е[-г,о]
3. Если ||-Сг|| ехр(Д/т) > 1, то для решения задачи (2.1) справедлива оценка
4 з=о
+ Рг||-1 тах{1, ||£>||,... , Ц^"1!!}) тах ехр Л.
) *€[—т,0] V 1т /
Доказательство этой теоремы проводится по аналогии с доказательством теоремы 6 из [11].
Следствие 1. Пусть выполнены условия теоремы 2, тогда нулевое решение системы (1.1) асимптотически устойчиво.
3. Заключение
В работе исследована асимптотическая устойчивость нулевого решения системы (1.1). В теореме 2 и следствии 1 получены достаточные
условия асимптотической устойчивости в терминах матричных неравенств и указаны три оценки решений системы (1.1) в различных случаях. Данные оценки характеризуют экспоненциальное убывание решений на бесконечности.
Список литературы
1. Азбелев Н. В., Симонов П. М. Устойчивость решений уравнений с обыкновенными производными. Пермь : Изд-во Перм. ун-та, 2001. 230 с.
2. Беллман Р., Кук К. Дифференциально-разностные уравнения. М. : Мир, 1967. 548 с.
3. Демиденко Г. В., Матвеева И. И. Об оценках решений систем дифференциальных уравнений нейтрального типа с периодическими коэффициентами // Сиб. мат. журн. 2014. Т. 55, № 5. С. 1059-1077.
4. Колмановский В. В., Носов В. Р. Устойчивость и периодические режимы регулируемых систем с последействием. М. : Наука, 1981. 448 с.
5. Красовский Н. Н. Некоторые задачи теории устойчивости движения. М. : Физматгиз, 1959. 212 с.
6. Мышкис А. Д. Линейные дифференциальные уравнения с запаздывающим аргументом. М.-Л. : Гостехиздат, 1951. 256 с.
7. Хейл Дж. Теория функционально-дифференциальных уравнений. М. : Мир, 1984. 424 с.
8. Эльсгольц Л. Э., Норкин С. Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. М. : Наука, 1971. 296 с.
9. Demidenko G. V. Stability of solutions to linear differential equations of neutral type // J. Anal. Appl. 2009. Vol. 7, N 3. P. 119-130.
10. Hayes N. D. Roots of the Transcendental Equation Associated with a Certain Difference-Differential Equation //J. London Math. Soc. 1950. Vol. 25. P. 226-232. https://doi.Org/10.1112/jlms/sl-25.3.226
11. Yskak Т. K. Sufficient conditions for the asymptotic stability of solutions to one class of linear systems of neutral type with periodic coefficients // Динамические системы. 2015. Т. 5(33), № 3-4. С. 177-191.
12. Yskak Т. Stability of solutions to systems of differential equations with distributed delay // Funct. Differ. Equ. 2018. Vol. 24, N 1-2. P. 97-108.
Тимур Кайратулы Ыскак, аспирант, Институт математики им. C.J1. Соболева СО РАН, Российская Федерация, 630090, Новосибирск, пр-т Академика Коптюга, 4; ассистент, Новосибирский государственный университет, Российская Федерация, 630090, Новосибирск, ул. Пи-рогова, 2; тел.: 8(913)9196662 (e-mail: [email protected])
Поступила в редакцию 10.08.18
On the Stability of Solutions of Neutral Differential Equations with Distributed Delay
T. K. Yskak
Sobolev Institute of Mathematics SB RAS, Novosibirsk State University, Novosibirsk, Russian Federation
Abstract. We consider one class of systems of linear nonautonomous differential equations of neutral type with distributed delay. The matrix in front of the derivative of the unknown vector-function with delay is constant, the matrix in front of the unknown vector-function has continuous T-periodic elements, the kernel of the integral operator consists of continuous functions, which are T-periodic with respect to argument t. The aim of the work is to study the asymptotic stability of the zero solution using the method of modified Lyapunov-Krasovskii functionals. The method of Lyapunov-Krasovskii func-tionals is the development of the Lyapunov second method. The advantage of this method is the simplicity of formulations and the reduction of the study of asymptotic stability to solving well-conditioned problems. In addition, the method of modified Lyapunov-Krasovskii functionals allows to obtain estimates of solutions to linear systems of delay differential equations. Note that the use of modified Lyapunov-Krasovskii functionals also allows to obtain estimates of solutions to nonlinear differential equations and estimates of the attraction set. Previously, a system of periodic differential equations of neutral type was considered in the works by G.V. Demidenko and I.I. Matveeva, where sufficient conditions of the asymptotic stability of the zero solution were obtained and estimates of solutions to this system were established. A system of linear periodic differential equations with distributed delay was considered by the author of this paper. For this system it was also obtained sufficient conditions of the asymptotic stability of the zero solution and established estimates of solutions. In the present paper, we obtain sufficient conditions of the asymptotic stability of the zero solution to the neutral type system with distributed delay in terms of matrix inequalities and establish estimates of solutions to the system characterizing the exponential decay at infinity.
Keywords: asymptotic stability, Lyapunov - Krasovskii functional, distributed delay, neutral type equation.
References
1. Azbelev N.V., Simonov P.M. Ustoychivost' resheniy uravneniy s obyknovennymi proizvodnymi [Stability of solutions to equations with ordinary derivatives]. Perm, Izdatelstvo Permskogo universiteta, 2001, 230 p. (in Russian)
2. Bellman R., Cooke K.L. Differential-difference equations. London, Academic Press, 1963, 482 p.
3. Demidenko G.V., Matveeva I.I. On estimates of solutions to systems of differential equations of neutral type with periodic coefficients. Sib. Math. J., 2014, vol. 55, no. 5, pp. 866-881. https://doi.org/10.1134/S0037446614050061
4. Kolmanovskiy V.V., Nosov V.R. Ustoychivost' i periodicheskie rezhimy reguliruemyh sistem s posledeystviem [Stability and periodic regimes of controlled systems with aftereffect]. Moscow, Nauka Publ., 1981, 448 p. (in Russian)
5. Krasovskii N.N. Stability of Motion. Applications of Lyapunov's Second Method to Differential Systems and Equations with Delay. Stanford, Stanford University Press, 1963, 188 p.
6. Myshkis A.D. Lineynye differentsial'nye uravneniya s zapazdyvayushchim argumentom [Linear differential equations with retarded argument]. Moscow, Leningrad, Gostehizdat Publ., 1951, 256 p. (in Russian)
7. Hale J. Theory of Functional Differential Equations. New York, Springer, 1977, 366 p. https://www.doi.org/10.1007/978-l-4612-9892-2
8. El'sgol'ts L.E., Norkin S.В. Vvedenie v teoriyu differentsial'nyh uravneniy s otklonyayushchimsya argumentom [Introduction to the theory of differential equations with deviating argument]. Moscow, Nauka Publ., 1971, 296 p. (in Russian)
9. Demidenko G.V. Stability of solutions to linear differential equations of neutral type. J. Anal. AppL, 2009, vol. 7, no. 3, pp. 119-130.
10. Hayes N.D. Roots of the Transcendental Equation Associated with a Certain Difference-Differential Equation. J. London Math. Soc., 1950, vol. 25, pp. 226-232. https://doi.Org/10.1112/jlms/sl-25.3.226
11. Yskak Т.К. Sufficient conditions for the asymptotic stability of solutions to one class of linear systems of neutral type with periodic coefficients. Dinamicheskie Systemy, 2015, vol. 5(33), no. 3-4, pp. 177-191.
12. Yskak T. Stability of solutions to systems of differential equations with distributed delay. Fund. Differ. Equ., 2018, vol. 24, no. 1-2, pp. 97-108.
Timur Yskak, Postgraduate, Sobolev Institute of Mathematics SB RAS, 4, Acad. Koptyug avenue, Novosibirsk, 630090, Russian Federation; Assistant, Novosibirsk State University, 2, Pirogov st., Novosibirsk, 630090, Russian Federation; tel.: 8(913)9139196662 (e-mail: [email protected])
Received 10.08.18