2019 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 15. Вып. 4
ПРИКЛАДНАЯ МАТЕМАТИКА. ИНФОРМАТИКА. ПРОЦЕССЫ УПРАВЛЕНИЯ
ПРИКЛАДНАЯ МАТЕМАТИКА
УДК 517.929 МБС 34К20
Устойчивость однородных нестационарных систем дифференциально-разностных уравнений с линейно возрастающим запаздыванием
А. В. Екимов, О. Н. Чижова, У. П. Зараник
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Екимов А. В., Чижова О. Н., Зараник У. П. Устойчивость однородных нестационарных систем дифференциально-разностных уравнений с линейно возрастающим запаздыванием // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2019. Т. 15. Вып. 4. С. 415-424. https://doi.org/10.21638/11702/spbu10.2019.401
Рассматриваются однородные дифференциально-разностные системы с периодическими коэффициентами и линейно возрастающим временным запаздыванием. Эти системы можно представить как модель распространения эпидемии среди населения. Кроме того, системы с линейно возрастающим запаздыванием описывают динамику работы информационного сервера, смесительного бака, процесс образования пробок на кольцевой дороге и т. д. Вводится понятие усредненной системы. Такой подход позволяет свести анализ задачи устойчивости по Ляпунову нулевого решения исходной системы к исследованию нулевого решения усредненной системы. Сформулированы достаточные условия устойчивости стационарной системы. К изучению устойчивости исходной системы применен подход Разумихина. Построена функция Ляпунова. В результате получены новые достаточные условия асимптотической устойчивости нулевого решения нестационарных однородных систем с линейно возрастающей задержкой времени. Эти условия являются обобщением известных результатов для линейных систем с линейно возрастающей временной задержкой.
Ключевые слова: однородная дифференциально-разностная система, линейное запаздывание, асимптотическая устойчивость.
Введение. В последнее время в центре внимания исследователей все чаще оказываются системы дифференциально-разностных уравнений, содержащие переменные, в том числе и неограниченные, линейно зависящие от времени, запаздывания аргумента. Это объясняется большим количеством ситуаций, когда запаздывание в модели нельзя считать постоянной величиной. Например, системы, описывающие взаимодействие объектов, расходящихся друг от друга, содержат транспортное за© Санкт-Петербургский государственный университет, 2019
паздывание аргумента — время, необходимое для передачи информации от одного объекта к другому, которое линейно возрастает при увеличении расстояния между объектами. Системы и уравнения с линейным запаздыванием аргумента появляются в математических моделях радиоактивного распада [1], работы информационного сервера [2] и смесительного бака [3], движения по кольцевой автомобильной дороге [4] и в ряде других случаев. Следует отметить, что многие из этих моделей являются нелинейными, причем разложение правых частей по степеням искомых функций не имеет линейных членов относительно них или же матрица при линейных слагаемых не позволяет использовать теоремы об устойчивости и неустойчивости по линейному приближению [1, 5, 6]. В этих случаях первым, в широком смысле, приближением оказываются системы уравнений с однородными правыми частями, что объясняет необходимость их всестороннего рассмотрения [7-9].
Заметим, что системы нелинейных уравнений, содержащие линейное запаздывание времени, до последнего времени остаются мало изученными.
Важной проблемой, возникающей при их исследовании, является проблема устойчивости, поскольку наличие в системе неограниченного запаздывания может привести к потере этого свойства. Теоремы об устойчивости однородных систем без запаздывания были доказаны в работах [10-12]. Уточнение известных критериев устойчивости по первому, в широком смысле, приближению систем без запаздывания приводится А. Ю. Александровым [13]. Отдельные классы уравнений и систем с однородной правой частью и линейным запаздыванием аргумента описаны в работах [2, 14, 15].
В данной статье рассматривается нелинейная нестационарная система уравнений с однородными периодическими правыми частями и линейным запаздыванием аргумента. Такая система, в частности, может служить математической моделью распространения эпидемии в некоторой популяции. Действительно, некоторые болезни, например ВИЧ, сифилис, герпес и др., имеют спящий период, длительность которого может быть неограниченной.
Для исследования тривиального решения на устойчивость были использованы метод функционалов Ляпунова—Красовского [16] и подход Разумихина [17], что позволило свести задачу к рассмотрению усредненной стационарной системы и получить некоторые достаточные условия асимптотической устойчивости нулевого решения.
Модель распространения эпидемии. Дифференциальные уравнения с однородными правыми частями могут применяться при моделировании самых различных процессов. Рассмотрим модель распространения эпидемии в некоторой популяции. Через x(t) обозначим число заболевших особей в момент времени t. Динамику процесса естественно описать дифференциальным уравнением
X = f (t,x) + g(t,x),
где f (t, 0) = g(t, 0)=0.
Функция f (t, x) характеризует факторы, способствующие распространению вируса. В простейших моделях полагают, что f (t, x) = ax, где a > 0 — постоянная. Пропорциональность f (t, x) числу заболевших существенно упрощает реальную картину, но далеко не всегда адекватно отражает ситуацию. Выберем в качестве f (t, x) следующую функцию:
f (t,x) = a(t)xM (t),
T
в которой a(t) — периодическая функция с периодом Т такая, что ^ / а(т)д,т = а > 0.
о
Постоянную а естественно считать коэффициентом заболеваемости. Постоянная ¡л > 0 характеризует скорость распространения эпидемии. Функция g(t, x) отражает факторы, препятствующие распространению эпидемии (природные факторы, укрепление иммунитета, изоляция заболевших особей и т. д.). Положим, что g(t,x) =
т
b(t)xa(t), где b(t + Т) = b(t); ^ f Ь(т)с1т = Ъ < 0. Постоянная b — коэффициент вы-
о
здоровления, а — постоянная. При ¡л = а получаем дифференциальное уравнение с однородной правой частью. Поскольку мероприятия по борьбе с эпидемией (например, вакцинация), как правило, не дают мгновенного эффекта, то естественно учесть запаздывание. В простейшем случае запаздывание можно считать постоянным. С учетом этого модель можно представить в виде
X = a(t)x»(t) + b(t)xa(t - т).
В дальнейшем возможен переход к переменному запаздыванию. В частности, если учитывать спящий период, то запаздывание можно считать линейно возрастающим, т. е.
x(t) = a(t)x»(t) + b(t)xa(at), 0 < a < 1.
Перейдем к описанию векторной модели распространения эпидемии. Пусть п популяций сосуществуют на общей территории, подверженной воздействию вируса. Вектор состояния системы обозначим через x(t) = (x\(t),... ,xn(t))T. Здесь xi(t) — число заболевших особей г-й популяции в момент времени t. Динамику процесса представим в виде системы дифференциальных уравнений с запаздыванием:
x(t) = A(t)f (x(t)) + B(t)g(x(at)), (1)
в которой A(t), B(t) — (ихп)-непрерывные матрицы с Т-периодическими элементами, описывающие внутрипопуляционные и межпопуляционные факторы, влияющие на распространение эпидемии. Вектор-функции f, g выбираем следующим образом:
f (x) = (xi,..., xn) , g(x) = (xi,..., xn) .
Примем, что
т т
0
Анализ устойчивости нулевого решения системы (1) в терминах параметров усредненной модели позволит выявить существенные факторы, влияющие на динамику эпидемии:
х(г) = А ■ /(х(г)) + В ■ д(х(аЬ)).
Предварительные результаты. Рассмотрим систему дифференциально-разностных уравнений с однородными правыми частями
х(г) = Г (г,х(г)) + С(х(аЬ)), Ь > г0 > 0, 0 <а< 1, (2)
где векторные функции Г(Ь,х), О(х) непрерывно дифференцируемы по своим аргументам и являются положительно однородными функциями порядка ¡л > 1, при этом ¡л — рациональное число с нечетными числителем и знаменателем.
Кроме того, будем считать функцию Е(Ь, х) периодической по первому аргументу, т.е. Е (Ь + Т,х) = Е (Ь,х) при всех Ь €
Известно [2], что нулевое решение уравнения
Х(Ь) = —рх^(Ь) + чх»(аЬ) (3)
асимптотически устойчиво по Ляпунову, если ¡1 = > 1, ¡3 > 0 и | —^| < /3.
Приведем некоторые достаточные условия асимптотической устойчивости для стационарной системы с однородными порядка ¡л > 1 правыми частями
х(Ь) = Е (х(Ь)) + С(х(аЬ)). (4)
Пусть указанные функции удовлетворяют оценкам
а-у\\хГ < \\Е(х)|| < а2\\хГ, Ъ^хГ < ПОДИ < Ь2\\хГ,
где а1, а2, Ьх, Ь2 — положительные постоянные.
Теорема 1 [17]. Если существует положительно определенная, положительно однородная порядка к > 2 функция V(х) € С2(Яп), удовлетворяющая условиям:
• ( ) ■Р(х) = ]¥(х) — отрицательно определенная функция,
) 'Т
I дх
тах < I
V(у)4У(х) ' ^ дх
то нулевое решение системы (4) асимптотически устойчиво по Ляпунову.
Замечание. Для функций, удовлетворяющих теореме 1, выполнены следующие неравенства:
Ь1\\х(аЬ)Г < тх(аЬ))\\ < Ь2\\х(аЬ)Г, С1\\х(Ь)\\к < V(х(Ь)) < С2\х(Ь)\к, дV
сзмтк-1 <
дх
< С4\х(Ь)\к-1, Р1\х(1)\к+^-1 < Ш(х(Ь)) < Р2\х(1)\к+^-1.
Здесь величины С1, С2, сз, С4, р1, р2 — положительные постоянные.
Следствие. Если функции V(х) и Ш(х) удовлетворяют условиям теоремы 1, то существует число Я > 1 такое, что справедливо неравенство
тах №)У - 1.
V (у)^№ (х) V дх
Пример 1. Уравнение (3) удовлетворяет условиям теоремы 1, если выбрать функцию V(х) = х4. Из [10] следует, что тривиальное решение уравнения (3) будет асимптотически устойчиво, если производная от этой функции, в силу уравнения (3), будет определенно отрицательной на решениях данного уравнения, подчиняющихся условию Разумихина, которое в рассматриваемом случае имеет вид
х4(аЬ) < х4 (Ь). (5)
Продифференцируем функцию V(х) = х4, в силу уравнения (3), и получим
—4/Зх^+3 (Ь) + 4^х3х^ (аЬ). (6)
¿V
(з)
Из условия (5) будет следовать, что |х(аЬ)| ^ |х(Ь)|. Теперь можно оценить правую часть равенства (6):
(Ш
< 4(-в + Ы)х»+3.
Очевидно, правая часть неравенства определенно отрицательна и удовлетворяет условиям теоремы 1 при выполнении неравенств в > 0 и |71 < в.
Пример 2. Рассмотрим систему уравнений с постоянными матрицами
¿(г) = Р
( (г) \ (г)
\ ¿1 (г)
+ Я
( ¿1 (аг) \ ¿1 (аг)
V ¿П(аг) /
(7)
Если существует положительно определенная диагональная матрица Л diag(Лl, .. .,Лп) такая, что
/
шах < (Л1 ¿1
Л ¿1)
Лп^п)
Р
\
( х1(г) \
X1(г)
V ¿1 (г)
+ Я
( у1 \
у1
\ уП /
< -¿Ы21 (8)
при условиях Лгу1+ ^ Л, то нулевое решение системы (7) будет асимпто-
г= 1 г=1
тически устойчивым.
Доказательство заключается в проверке условий теоремы 1 для функции
п
V(¿) = £ Лix1+ 1.
г= 1
Дифференцируя функцию V(¿), в силу системы (7), получаем левую часть неравенства (8). Так как оценка добавочного слагаемого (^ +1)(Л Л^1 ... ЛптП )х хЯХ 1(аг) через неравенство Разумихина приводит к выполнению (8), это означает, что производная от функции V(¿), в силу системы (7), определенно отрицательна. Следовательно, нулевое решение системы (7) будет асимптотически устойчивым. Основные результаты. Скалярным случай. Рассмотрим уравнение
ж = -в(г) • ¿3(г) + ^я3(аг),
(9)
где /3(г) — периодическая функция с периодом Т, т. е. /3(г + Т) = (3(г); 7 — постоянная. Введем среднее значение функции (3 = ^ /0 /3(£)сЙ.
Теорема 2. Если выполняется условие |7| < в, то нулевое решение уравнения (9) асимптотически устойчиво по Ляпунову. Доказательство. Рассмотрим функцию
V(г,т(г)) = ¿4(г) - 4т6(г) (в - ¡з(г)) ¿т.
Найдя производную этой функции вдоль решений уравнения (9), получим, что
ду{1,х{1)) м
= 4т3(г) • (-/З(г):с3(г) + 73?(аь)) - 4т6(г) (в - в(г)) -
г
- 4 (Х6(г)) • / в - Р(т) ¿т = -4вx6(г) + 47x3(г) • ¿3(аг) -0 ^ '
- 24x5(t) —(t) ■ x3(t) + Yx3(at)) ■ j (j3 - 3(т)) dr.
Введем величину [3(Ь) = / I в — [3(т)) ¿т, которая, очевидно, будет удовлетворять
о ^ '
условиям [3(Ь + Т) = [3(Ь) и |/3(Ь)1 ^ В < ж. Тогда на множестве |х(Ь)| < Н < ж будут справедливы следующие оценки:
х4(Ь) (1 — 4Н2В) < V(Ь, х(Ь)) < х4(Ь) (1 + 4Н2В) .
Далее определим из условия Разумихина V(а,х(а)) ^ шУ(Ь,х(Ь)), а ^ Ь, т> 1, что
х (аЬ) < тх (¿) ^ _4Я2д ) ,
Введем теперь величину М = откуда следует, что |ж(ой)| ^ М\х(Ь) \.
Тогда можно оценить производную функции V(Ь,х(Ь)) на решениях уравнения (9):
dV{t,x{t)) dt
< 4(-в + y)x6(t) + 24B(f3(t) + YM)x8(t) <
(9) 6 г - 2 1 (10)
< -4x6(t) (в - y) - (6B3+ yM)x (t)
Так как |x(t)| < H, то постоянную H можно считать настолько малой, что величина в квадратных скобках в неравенстве (10) будет положительной. Применяя теорему 1, завершаем доказательство.
Общий случай. Теперь будем исследовать систему (2)
x(t) = F (t,x(t)) + G (x(at)), t > t0 > 0, 0 <a< 1.
T
Введем величину F(x) = ^ f F(t, x(t))dt.
T0
Теорема 3. Если для системы
x(t) = F(x(t))+ G(x(at))
выполнены условия теоремы 1, то нулевое решение системы (2) асимптотически устойчиво по Ляпунову.
Доказательство. Перепишем систему (2) в виде
x(t) = F(x(t)) + G(x(at)) + F(t, x(t)) - F(x(t)). (11)
Рассмотрим функцию
v(t,x(t)) = V(x(t))- (^))})T/ (F(T,x(t))~Hx(t)j)dT, (12)
0
где V(x) — функция из теоремы 1; верхний индекс T обозначает транспонирование.
Введем векторную функцию
г
¥^,¿(1)) = ! (¥(т^г) - Р^г))^ ¿т. (13)
0
Можно заметить, что функция (13) будет Т-периодической, ограниченной и однородной, при этом д^^г)]]1 < ^(г, ¿(г))]] < д^^г)]]1, где величины д1, д2 — положительные постоянные, а функция (12) удовлетворяет оценкам
с1Ыг)\]к (1 - с^Ыг)!!1-1) < v(г,x(г)) < с2Ыг)\]к (1 +с^Ыг)]]1-1)
и, следовательно, является определенно положительной в достаточно малой окрестности нуля ^¿(г)! < Но. Производная функции (12) на решениях системы (11) будет иметь вид
. т
¿М*))|(п) = (Ё}Шй) [Р(х(г)) + С(х(аг))
т г
_Ш\т?(+ I 1
дх(г)
0
Введем матрицу
г
*(*)) + С(хМ))]) / (^(т, х{Ь)) - ВД*))) ¿т - (14)
/ + ОД«*))]-
д V(^¿(г)) - Р^г)) Р(г,х(г)) = I -—-'-¿т. (15)
0
Матрица (15) будет Т-периодической, ограниченной и однородной, при этом выполнено \\Р(г,¿(г))] ^ Ц^г)]]1-1, где I — положительная постоянная величина.
В условиях теоремы 3 можно оценить сверху равенство (14) на множестве V(y) < ЯУ^):
(11) т (16) " (^Щг) *(*))№, х®) +С(х(ат.
Условие Разумихина [17] для функции (12) запишется в виде v(а, ¿(а)) ^ v(г, ¿(г)) при а ^ г или
т ®
У(х(а)) - / (Р(т, х(а)) - Р(х(а))) ¿г <
0
< - (^|)))))Г / (П^М)-Яхй)) ¿т.
0
Учитывая введенные ранее оценки, получим, что
С1\\х(а)\\к - ^ сЛх{Щь ^ + .
Обозначим А = ^ > 0 и выберем величину Н = 1/ "~\/2А, тогда 1 - АН^1 = 1/2, Вестник СПбГУ. Прикладная математика. Информатика... 2019. Т. 15. Вып. 4 421
Пусть В = тах{ ^/Зсг/сх, ^/йсТ/сг} и = Н/В. Тогда при условии ||ж|| < будут выполнены неравенства ||ж(<т)|| ^ Зс2/сх||)|| и У(ж(<т)) ^ ДУ(ж(£)) для значений <т ^ С учетом того, что ^ а также выполнены неравенства ^
Сб\х(Ь)Ук-2 и Е(Ь,х(Ь)) + С(х(аЬ)) ^ (а2 + ВмЬ2)\\х\\м, проведем дальнейшую оценку правой части неравенства (16):
(11) 2
-KiMm^-1 (i-^^ikwr-1)
Поскольку л > 1, найдется величина Н2 > 0 такая, что при \\х(Ь)\\ < Н2 будет выполнено неравенство 1 — > 0.
Обозначим Н = тт{Я0,Я1,Я2}. Тогда при ||ж(г)|| < Д"^/"2(1+0492^-ПР0ИЗ" водная определенно положительной функции (12) на решениях системы (2) будет определенно отрицательной величиной, что доказывает теорему.
Заключение. Исследована нелинейная нестационарная система уравнений с однородными периодическими правыми частями и линейным запаздыванием аргумента. С помощью функционалов Ляпунова—Красовского и метода Разумихина построена усредненная стационарная система, для которой получены достаточные условия асимптотической устойчивости нулевого решения.
Литература
1. Беллман Р., Кук К. Л. Дифференциально-разностные уравнения / пер. с англ. А. М. Звер-кина, Г. А. Каменского; под ред. Л. Э. Эльцгольца. М.: Мир, 1967. 548 с. (Bellman R., Cooce K. L. Differencial-difference equations.)
2. Жабко А. П., Чижова О. Н. Анализ устойчивости однородного дифференциально-разностного уравнения с линейным запаздыванием // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2015. Вып. 3. С. 105—115.
3. Жабко А. П., Чижова О. Н. Гибридный метод анализа устойчивости линейных дифференциально-разностных систем с линейно возрастающим запаздыванием // Вестн. Тамбов. гос. ун-та. Сер. Естественные и технические науки. 2015. Т. 20. № 4. C. 843—850.
4. Zhabko A., Chizhova O., Zaranik U. Stability analysis of the linear time delay systems with linearly increasing delay // Cybernetics and Physics. 2016. Vol. 5. N 2. P. 67—72.
5. Валеев К. Г. Линейные дифференциальные уравнения с запаздыванием, линейно зависящим от аргумента // Сиб. матем. журн. 1964. Т. 5. № 2. С. 75—83.
6. Laktionov A. A., Zhabko A. P. Method of difference transformations for differential systems with linear time-delay // Proceeding of the First IFAC Conference LTDS-98. Grenoble, France, 1998. P. 201-205.
7. Купцова С. Е. Об асимптотическом поведении решений систем нелинейных нестационарных дифференциальных уравнений // Труды Средневолжск. матем. об-ва. 2006. Т. 8. № 1. С. 235-243.
8. Тихомиров О. Г., Темкина Е. В. Асимптотическое положение покоя для систем однородных нестационарных дифференциальных уравнений // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2014. Вып. 3. С. 58-65.
9. Жабко А. П., Тихомиров О. Г., Чижова О. Н. Устойчивость асимптотического положения покоя возмущенных однородных нестационарных систем // Журн. Средневолжск. матем. об-ва. 2018. Т. 10. № 1. C. 13-22.
10. Зубов В. И. Устойчивость движения. М.: Высшая школа, 1973. 272 с.
11. Красовский Н. Н. Об устойчивости по первому приближению // Прикладная математика и механика. 1955. Т. 19. № 5. С. 516-530.
12. Малкин И. Г. Теория устойчивости движений. М.; Л.: Гостехиздат, 1952. 432 с.
13. Александров А. Ю. Устойчивость движений неавтономных динамических систем. СПб.: Изд-во С.-Петерб. ун-та, 2004. 184 с.
14. Александров А. Ю., Жабко А. П. Об асимптотической устойчивости решений нелинейных систем с запаздыванием // Сиб. матем. журн. 2012. Т. 53. № 3. С. 495—508.
15. Гребенщиков Б. Г., Клечин Ю. И. Об устойчивости одной однородной нестационарной системы с линейным запаздыванием // Дифференциальные уравнения. 2003. Т. 39. № 12. С. 1600— 1604.
16. Красовский Н. Н. Некоторые задачи теории устойчивости движения. М.: Физматгиз, 1959.
212 с.
17. Разумихин Б. С. Устойчивость систем с запаздыванием // Прикладная математика и механика. 1956. Т. 20. № 4. С. 500-512.
Статья поступила в редакцию 6 июля 2019 г. Статья принята к печати 7 ноября 2019 г.
Контактная информация:
Екимов Александр Валерьевич — канд. физ.-мат. наук, доц.; a.ekimov@spbu.ru
Чижова Ольга Николаевна — канд. физ.-мат. наук, доц.; o.chizhova@spbu.ru
Зараник Ульяна Петровна — канд. физ.-мат. наук, ст. преподаватель; u.zaranik@spbu.ru
Stability of homogeneous nonstationary systems of differential-difference equations with linearly time delay.
A. V. Ekimov, O. N. Chizhova, U. P. Zaranik
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Ekimov A. V., Chigova O. N., Zaranik U. P. Stability of homogeneous nonstationary systems of differential-difference equations with linearly time delay. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2019, vol. 15, iss. 4, pp. 415-424. https://doi.org/10.21638/11702/spbu10.2019.401 (In Russian)
These systems can be considered as a model for the spread of the epidemic in the population. In addition, systems with linearly increasing delay describe the dynamics of the information server, mixing tank, the process of formation of traffic jams on the ring road, etc. For the study, the concept of an average system is introduced. This approach allows us to reduce the analysis of the Lyapunov stability problem of the zero solution of the original system to the investigation of the zero solution of the averaged system. Sufficient conditions for stationary system stability are formulated. Then the application of Razumihin's approach to the study of stability original system is used. The Lyapunov function is constructed. As a result, new sufficient conditions for the asymptotic stability of the zero solution of nonstationary homogeneous systems with a linearly increasing time delay are obtained. These conditions are the generalization of well-known results for the linear systems with a linearly increasing time delay.
Keywords: homogeneous differential-difference system, linearly increasing time delay, asymptotic stability.
References
1. Bellman R., Cooke K. L. Differentsialno-raznostnyie uravneniya [Differential-difference equations]. Moscow, Mir Publ., 1967, 548 p. (In Russian)
2. Zhabko A. P., Chizhova O. N. Analiz ustoychivosti odnorodnogo differentsialno-raznostnogo uravneniya s lineynyim zapazdyivaniem [Stability analysis of homogeneous differential-difference equation
with linear delay]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2015, iss. 3, pp. 105—115. (In Russian)
3. Zhabko A. P., Chizhova O. N. Gibridnuj metod analiza ustoichivosti linejnuh differencialno-raznostnuh sistem s linejno vozrastaushchim zapazdyvaniem [Hybrid method of stability analysis of linear difference-differential systems with linear increasing decelerating]. Tambov University Reports. Series Natural and Technical Sciences, 2015, vol. 20, no. 4, pp. 843—850. (In Russian)
4. Zhabko A., Chizhova O., Zaranik U. Stability analysis of the linear time delay systems with linearly increasing delay. Cybernetics and Physics, 2016, vol. 5, no. 2, pp. 67—72.
5. Valeev K. G. Lineynyie differentsialnyie uravneniya s zapazdyvaniem, lineyno zavisyaschim ot argumenta [Linear differential equations with linear time-delay]. Siberian Mathematical Journal, 1964, vol. 5, no. 2, pp. 75—83. (In Russian)
6. Laktionov A. A., Zhabko A. P. Method of difference transformations for differential systems with linear time-delay. Proceeding of the First IFAC Conference LTDS-98. France, Grenoble, 1998, pp. 201— 205.
7. Kuptsova S. E. Ob asimptoticheskom povedenii reshenij sistem nelinejnuh nestacionarnuh dif-ferencialnuh uravnenij [On the asymptotic behavior of nonlinear nonstationary differential equation solutions]. Proceeding of Middle-Volga Mathematical Society, 2006, vol. 8, no. 1, pp. 235—243. (In Russian)
8. Tikhomirov O. G., Temkina E. V. Asimptoticheskoe polojenie pokoia dlia sistem odnorod-nyhnestatsyonarnyh differentialnyh uravnenii [Asymptotic quiescent position for systems of homoge-neousnon-autonomous differential equations]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2014, iss. 3, pp. 58—65. (In Russian)
9. Zhabko A. P., Tikhomirov O. G., Chizhova O. N. Ustoychivost asimptoticheskogo polozheniya pokoya vozmuschennyh odnorodnyh nestatsionarnyh sistem [Stability of the asymptotic quiescent position of perturbed homogeneous nonstationary systems]. Proceedings of Middle-Volga Mathematical Society, 2018, vol. 10, no. 1, pp. 13-22. (In Russian)
10. Zubov V. I. Ustoichivost dvizheniya [Stability of motions]. Moscow, Visshaya shkola Publ., 1973, 272 p. (In Russian)
11. Krasovskij N. N. Ob ustoichivosti po pervomu priblizheniyu [On stability by the first approximation]. Journal of Applied Mathematics and Mechanics, 1955, vol. 19, no. 5, pp. 516-530. (In Russian)
12. Malkin I. G. Teoriya ustojchivosti dvizhenij [The theory of motion stability]. Moscow, Leningrad, Gostekhizdat Publ., 1952, 432 p. (In Russian)
13. Alexandrov A. Yu. Ustoichivost dvizhenij neavtonomnyh dinamicheskih sistem [The stability of motion of nonautonomous dynamical systems]. St. Petersburg, Saint Petersburg State University Publ., 2004, 184 p. (In Russian)
14. Alexandrov A. Yu., Zhabko A. P. Ob asimptoticheskoy ustoychivosti resheniy nelineynyh sistem s zapazdyvaniem [On the asymptotic stability of nonlinear system solutions with delay]. Siberian Mathematical Journal, 2012, vol. 53, no. 3, pp. 495-508. (In Russian)
15. Grebenschikov B. G., Klechin Yu. I. Ob ustoychivosti odnoy odnorodnoy nestatsionarnoy sistemy s lineynym zapazdyvaniem [Stability of a homogeneous nonstationary system with linear delay]. Differential Equations, 2003, vol. 39, no. 12, pp. 1600-1604. (In Russian)
16. Krasovskij N. N. Nekotoryie zadachi teorii ustoychivosti dvizheniya [Some problems of the theory of stability of motion]. Moscow, Fizmatgiz Publ., 1959, 212 p. (In Russian)
17. Razumikhin B. S. Ob ustoichivosti sistem s zapazdyvaniem [Stability of systems with delay]. Journal of Applied Mathematics and Mechanics, 1956, vol. 20, no. 4, pp. 500-512. (In Russian)
Received: July 06, 2019.
Accepted: November 07, 2019.
Author's information:
Alexander V. Ekimov — PhD in Physics and Mathematics, Associate Professor; a.ekimov@spbu.ru
Olga N. Chizhova — PhD in Physics and Mathematics, Associate Professor; o.chizhova@spbu.ru
Uliana P. Zaranik — PhD in Physics and Mathematics, Senior Lecture; u.zaranik@spbu.ru