УДК 517.929 Вестник СПбГУ. Прикладная математика. Информатика... 2023. Т. 19. Вып. 2
MSC 93D09
Метод Разумихина в обобщенной задаче Мышкиса для систем с распределенным запаздыванием*
А. В. Егоров
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 74)
Для цитирования: Егоров А. В. Метод Разумихина в обобщенной задаче Мышкиса для систем с распределенным запаздыванием // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2023. Т. 19. Вып. 2. С. 148161. https://doi.org/10.21638/11701/spbul0.2023.202
В статье представлены достаточные условия разрешимости обобщенной задачи Мышкиса для системы уравнений с распределенным переменным запаздыванием и постоянным ядром. Получены условия на ядро, при выполнении которых система будет равномерно устойчива при любом допустимом запаздывании. Под допустимым запаздыванием понимается кусочно-непрерывная функция, ограниченная сверху по величине и скорости роста. На примере двух систем показано, как можно проверить найденные условия на практике.
Ключевые слова: система с запаздыванием, устойчивость, распределенное запаздывание, обобщенная задача Мышкиса, метод Разумихина.
1. Введение. В статье рассматривается система дифференциальных уравнений с распределенным ограниченным запаздыванием, зависящим от времени. Продолжается разработка идеи применения модифицированного метода функций Разумихина для решения обобщенной задачи Мышкиса.
Под обобщенной задачей Мышкиса мы понимаем задачу определения параметров системы, при которых равномерная устойчивость наблюдается для любых запаздываний из заданного ограниченного множества пространства кусочно-непрерывных функций. Можно сказать, что решается задача исследования робастной устойчивости, где параметры варьируются в пределах бесконечномерных множеств. Для простейшего скалярного уравнения
в работе А. Д. Мышкиса [1] был получен весьма любопытный результат: в случае, когда ЪН е [0, 3/2], при любом запаздывании т, для которого
все решения уравнения (1) ограничены, что, учитывая линейность, можно трактовать как устойчивость по Ляпунову. Равномерная устойчивость при указанных условиях доказана в работе [2].
Кроме того, в [1] показано, что при ЪН < 3/2 уравнение асимптотически устойчиво, а при ЪН > 3/2 найдется запаздывание, ограниченное константой Н, для ко-
* Исследование выполнено при финансовой поддержке Российского научного фонда (грант № 21-71-00052).
© Санкт-Петербургский государственный университет, 2023
x(t) = -bx(t - т(t))
(1)
0 ^ т(t) ^ h при всех t,
(2)
торого уравнение неустойчиво. Таким образом, 3/2 является точной верхней границей множества, которое мы называем множеством разрешимости задачи Мышкиса. В иностранной литературе (см., например, [3]) данная теорема называется 3/2 stability theorem.
Теорема Мышкиса была замечена далеко не сразу, но с 1970-х годов стали появляться статьи с обобщениями этого красивого и нетривиального результата. В частности, нелинейное уравнение с переменным запаздыванием было исследовано в [2], уравнение с зависящим от времени параметром b рассмотрено в [3], случай, когда кроме запаздывающего члена в уравнении присутствует слагаемое без запаздывания, представлен в [4], уравнению с несколькими сосредоточенными запаздываниями посвящены работы [5, 6], а случаю распределенного запаздывания — статья [7].
Во всех перечисленных работах, как и в статье Мышкиса [1], исследовались скалярные уравнения с вещественными коэффициентами. Нашей целью с самого начала было изучение систем из нескольких уравнений. При внимательном рассмотрении мы пришли к выводу, что методы, представленные в перечисленных выше статьях, не годятся для такого исследования, поэтому был разработан новый подход на основе метода функций Разумихина.
Сначала метод был успешно применен для системы из двух вещественных уравнений, соответствующей скалярному уравнению (1) с комплексным коэффициентом b
был распространен на системы с произвольным числом уравнений и системы, в которых присутствуют как запаздывающий член, так и слагаемое без запаздывания [9, 10]. Однако до настоящего времени мы рассматривали только случай сосредоточенного запаздывания, в этой статье представленный метод будет впервые применен для системы с распределенным запаздыванием.
Стоит сказать, что предложенный подход имеет ряд ограничений. Во-первых, пока не удалось показать ни в одном из рассмотренных случаев, что полученная оценка области разрешимости задачи Мышкиса является точной, как в случае с теоремой о трех вторых. К сожалению, разработанный нами метод не позволяет доказывать обратные утверждения об устойчивости, и эта проблема пока не решена. Во-вторых, метод требует от запаздывания не только ограниченности по величине, но и ограниченности скорости роста единицей. Стоит отметить, однако, что скорость убывания запаздывания может быть любой, в том числе и бесконечной, т. е. запаздывание может иметь скачки вниз, но не может иметь скачков вверх. В-третьих, прямая проверка условий из представленных нами теорем в общем случае является довольно непростой задачей, так как они выражены в терминах пустоты некоторых множеств, описываемых невыпуклыми системами неравенств. Для того чтобы обойти эту сложность, для систем с сосредоточенными запаздываниями будут предложены некоторые аналитические процедуры, а также метод, позволяющий свести проблему к нескольким задачам выпуклого квадратичного программирования, которые, однако, вносят некоторый консерватизм в результаты.
Отметим, что для решения обобщенной задачи Мышкиса может быть также использован метод функционалов Ляпунова — Красовского (для систем с точечным запаздыванием см. [1ИЗ]), но при его применении от запаздывания требуется непрерывная дифференцируемость, а метод функций Разумихина, как правило, не нуждается в этом ограничении. Стоит сказать, что запаздывание, выбранное Мышки-сом для доказательства неустойчивости уравнения (1) при bh > 3/2, есть разрывная функция со скачками вниз.
В настоящей работе рассмотрена система уравнений с распределенным запаздыванием, ядром которого является постоянная матрица. Получены достаточные условия разрешимости обобщенной задачи Мышкиса в терминах пустоты множества решений некоторой системы неравенств. Отдельно исследован частный случай скалярного уравнения с вещественным коэффициентом, для которого найдены явные условия разрешимости обобщенной задачи Мышкиса, которые позволили сравнить представленный нами результат с известным условием [7]. Кроме того, на примере системы из двух уравнений показано, как установленные достаточные условия могут быть проверены на практике.
В статье используются следующие обозначения. Для векторов применяется евклидова норма || • ||, для матриц — индуцированная ею норма. Симметризованная версия для квадратной вещественной матрицы А обозначается через
А
н ЛТ + А 2
Пространство непрерывных функций со значениями в М", заданных на от резке [а, в], обозначается С(0)([а, в], М").
2. Система с распределенным запаздыванием. Рассмотрим систему
о
х(г) = А ! х(г + в) ¿в, (з)
-г (4)
в которой х(г) е М" для всех г, запаздывание т подчиняется условию (2). Если масштабировать время, то систему (3) можно свести к случаю Н = 1. И действительно, легко показать, что следующая система эквивалентна исходной в смысле устойчивости:
о
у(г) = Н2А ! у(г + в) ¿в,
-т(г)
где для всех г
т (Нг) г п
Поэтому далее, не уменьшая общности, будем считать, что запаздывание в системе (3) ограничено единицей. Добавив к этому условию требование ограниченности скорости роста запаздывания, которое необходимо в силу особенностей применяемого нами метода, сформулируем строгое определение допустимого запаздывания.
Определение 1. Кусочно-непрерывную функцию т(г), г е [—1, то), будем называть допустимым запаздыванием, если она удовлетворяет таким условиям: 1)
0 < т(г) < 1, г > —1;
2)
т (г2) — т (г1)
г2 — г1
< 1, г2 >г1 > —1.
В качестве начальных будем рассматривать непрерывные функции, заданные на отрезке [— 1,0], начальный момент времени для определенности будем считать
неотрицательным. Через х(Ъ, ¿о, ф) обозначим решение системы (3), удовлетворяющее начальному условию
х(^о,ф)= - ¿о), Ъ € [¿о - Мо]-
Определение 2. Систему (3) с допустимым запаздыванием т будем называть равномерно устойчивой, если для любого е > 0 найдется 5 > 0, такое что \\х(Ъ,Ъо,ф)\\ < е для всех ¿о ^ 0 ¿о и всех начальных функций ф, удовлетворяющих условию вир^^ о] \\ф(в)\\ < 5.
Теперь сформулируем основное определение.
Определение 3. Будем говорить, что матрица А решает обобщенную задачу Мышкиеа для системы (3), если система с этой матрицей равномерно устойчива при любом допустимом запаздывании.
3. Метод Разумихина. В классической формулировке метод Разумихина требует для равномерной устойчивости системы существования положительно определенной функции V, производная которой в силу системы неположительна на всех непрерывных функциях из множества
^(ф(в)) <v(ф(0)), в € [-1, 0]}.
Применим модификацию этой теоремы. Идея заключается во введении искусственного запаздывания N > 1. При этом па интервале [-N -Я + 1] задаются произвольные начальные функции, а на интервале (-Я + 1,0] они доопределяются с помощью уравнения (3). Ниже приведена строгая формулировка этой модификации теоремы Разумихина. Результат базируется на теореме со с. 255 книги [14]. Мы адаптировали ее для случая нестационарной системы, взяв в качестве функции Разумихина простейшую положительно определенную функцию \\ • \\2. Воспользовавшись линейностью системы, ограничимся нормированными функциями.
Теорема 1. Система (3) с допустимым запаздыванием равномерно устойчива, если существует N > 1 такое, что для любого Ъ > 0 и любой функции ф €
о
(0)А J ф(в) ¿в < 0. (4)
Здесь
Х,м = | Ф € С(о)([-Я, 0], М")
-г (г)
\ф(0)\\ = 1, \ф(в)\\ < 1, в € [-Я, 0), в
ф(в)= А I ф(в) ¿в, в € (-Я +1,0)|.
в, Ив
в-т (г+в)
Далее будем использовать эту теорему при N = 3.
4. Вспомогательные результаты. Чтобы применить модифицированную теорему Разумихина для построения достаточных условий разрешимости обобщенной задачи Мышкиеа, введем два вспомогательных результата.
Сначала построим оценку сверху для выражения в левой части неравенства (4).
Лемма 1. Пусть т — допустимое запаздывание. Для любого г > 0 на множестве з выполнено неравенство
рт(0)А J ф) ¿в < т(г) рт (0)АНр(0) +
2
(г)1|А2||
-т (4) 1
г3 (г)||А3||
У (1 — в2) шах {0, х(1 + т(г)(в — 1), — т(г)*} ¿в,
0 где
(1) = !И2|| + ут(0)(А2)нф) (1) = \\АЦ+^(0)АМ-т
V ||А3|| ' ||А3||
Доказательство. Зафиксируем г ^ 0. Если т (г) = 0, результат очевиден. т( г) > 0
в
(1 ~сШ
в-т (4+в)
I ¥>(*)<** (5)
из определения множества та интервале от £ до 0, где £ — некоторая точка из [—т(г), 0]
о в
р(0) — ф) = а! У р(в) ¿в¿в. (6)
е в-т(4+в)
Отсюда после домножения слева па рт(0)А, симметризации первого слагаемого и ин-
[—т(г), 0]
0 о о в
У рт(0)Аф) ¿£ = т(г) рт(0)АНр(0) — рт(0)А2 I I I ф) ¿s¿в¿£.
-Т(4) -т(4) е в-т(4+в)
Изменение порядка интегрирования в последнем его слагаемом приводит к формуле
о о в
У рт(0)АФ) ¿£ = т(г) рт(0)АНр(0) — I (в + т(г)) У рт(0)А2р(в) ¿в ¿в. (7)
-т (4) -т (4) в-т(4+в)
Чтобы получить оценку сверху, нужно оценить снизу произведение рт(0)А2р(в) для в е [—2,0]. Сделаем это тремя способами, а после соберем из трех оценок наилучшую.
Оценка 1. Простейшее неравенство можно получить из условия ||р(в)|| < ||р(0)||, в е [—3,0), го определения множества Х4з:
рт(0)А2ф) > —м0)||||А2||м*)|| > —||а2||. Оценка 2. Из равенства (6) получаем
о в
рт(0)А2ф) = рт(0)(А2)Нр(0) — рт(0)А3 J у ф) ^¿в > рт(0)(А2)Нр(0) —
Я в-т(4+в)
о
о в о
-\\А3\\1 I \ф(п)\\ ¿пйв>фТ (0)(А2 )нф(0)-\\А3\\1 т(г + в) ¿в.
я в-т(г+в) я
Так как допустимое запаздывание ограничено единицей, имеем неравенство
ФТ(0)А2ф(в) > фТ(0)(А2)Нф(0) + в\\А3\\ = — \А2\ + X1 + ^ \\А3\\.
Оценка 3. Проинтегрируем равенство (5) на интервале от -т(¿) до в и применим процедуру, аналогичную той, что и при построении оценки 2. Однако стоит учесть, что значение в может быть как больше, так и меньше, чем -т(¿):
я в
ФТ(0)А2ф(в) = ФТ(0)А2ф(-т(¿)) + фТ(0)А3 ! у ф(п) ¿^¿в >
-т(г) в-т(г+в)
>ФТ(0)А2ф(-т(¿)) - \\А3\
в
> \А2\ + (у$-|в + т(г)|) \\А3\\.
J т (г + в) ¿в
-т(г)
Невозможно однозначно сказать, какая из трех оценок лучше, поэтому объединим их в одну:
ФТ(0)А2ф(в) > -\\А2\\ + \\А3\\ шах{0, хV + в, у^ - 1в + т(
В итоге из формулы (7) выводим неравенство
о о
-Тт\ I п /|2|
ФТ(0)А у ф(4) ¿4 < т(г) фТ(0)АНф(0) + \\А2\\ ! (в + т(г))т(г + в) ¿в -
-т(г) -т(г)
о в
-\\А3\\ I (в + т(г)) I ша^0, х<1) + в,у^ - 1в + т(*)|} ¿s¿в.
-т(г) в-т(г+в)
в-
т(г + в) ^ -т(г). Так как выражение в последнем интеграле полученной выше оценки неотрицательно, можно заменить нижний предел интегрирования в - т(г + в) на
-т (г):
т (г)
Т Н 2
ФТ(0)А у ¿4 < т(г) фТ(0)АНф(0) + \\А2\\ у в ¿в -
-т(г) о
о в
-\\А3\\ I (в + т(г)) I ша^0, х<1) + в, у£
- |в + т(г)| [• ¿в¿в.
-т(г) -т(г)
После изменения порядка интегрирования и простейших преобразований получим неравенство
то
ТН
ФТ(0)А у ф(4) ¿4 < т(г) фТ (0)АНф(0) +
г2(г)\\А2\\
-т(г)
т (4) в-т (4)
— ||А3||| в I шах{0, х« + в, у(14 — в + т(г)|} ¿в ¿в < т (г) рт (0)АНр(0) +
о -т(4)
- (4)
2
которое преобразуется к искомому с помощью замены переменной интегрирования. □ Можно заметить, что полученная оценка монотонна по х^ и у^. Так как необходимо ориентироваться на худший возможный случай, можно взять наименьшие (1) (1) (1) (1) значения для х^ и т. е. х^ = =0. Но если хотим получить менее консервативный результат, можно вывести некоторое ограничение на р(0) и р(—т(г)), которые
(1) (1) задают х^ и у^^-
Лемма 2. Пусть т — допустимое запаздывание. Для любого г > 0, такого что т(г) > 0, на множестве Х4,3 выполнено неравенство
рт(0)р(—т(г)) > 1 — т(г)||А|| +
1
+ т2(г)||А2||У(1 — в) шах{0, х^ + т(г)(в — 1), у(24 — т(г)*} ¿в,
о
где
(2) \\А\\-<рт(0)Ан<р(0) (2) \\А\\-<рт(0)А<р(-т№ х<р \\А21| ' ||А2||
Доказательство. Из формулы (6) получаем равенство
то тв
рт(0)р(—т(г)) = 1 — J ! рт(0)Ар(в) ¿в ¿в. (8)
-т(4) в-т(г+в)
Так как требуется оценка снизу, оценим произведение рт(0)Ар(в) сверху при в е [—2, 0]
Оценка 1:
рт(0)Ар(в) < ||А||.
Оценка 2:
то тв
рт(0)А<р(в) = рт(0)АНр(0) — рт(0)А2 J ! р(п) ¿ц ¿в <
я в-т(4+в)
< рт(0)АНр(0) — в||А2|| = ||А| + Рт(0)АНр(0) — ||А| — в||А2||. Оценка 3:
Я в
рт(0)Ар(в) = рт(0)Ар(—т(г))+ рт(0)А2 I ! р(п) ¿^¿в <
-т(4) в-т(4+в)
< рт(0)Ар(—т(г)) + ||А2|| |в + т(г)| = ||А| + рт(0)Ар(—т(г)) — ||А| + |в + т(г)| ||А2||.
Из полученных трех оценок составляем одну:
ФТ(0)Аф(в) < \\А\\ - \\А2\\ шах {0, х™ + в, у$ -|в + т(*)|} . Подставив ее в уравнение (8), получим неравенство
о в
(0)у(-т(t)) = 1 - J J (0)Ay>(s) dsdO >
-T (t) в-т (t+в) 0 0 в
> 1 -\\A\\ j т (t + в) de + \\A2\\ j j maxjo, x(2) + s, y(2) -¡s + т (t)|} dsdû.
)dû
-T (t) -т (t) в-т (t + в)
Так как в + т(t) < т(t + в) < 1,
(0)<р(-т(t)) >
о в
[ i „(2) , „ „(2)
> 1 - т (t)WAW + \\A2W J J maxjo, x™ + s, y£t -¡s + т (t)|} ds de =
-T (t) -T (t) 0
= 1 - т(t)\\A\\-\\A2\\ y s ma^0,x^ + s, y(2,t - s - т(t)} ds.
I2" I в шах|0, х(2) + в, у(2) -т(г)
Замена переменной интегрирования в на т(г)(в - 1) приводит к требуемому неравенству. □ 5. Основной результат.
А
решении (т, 7, /) € (0,1] х М" х М" следующей системы неравенств пусто:
1Т1 = 1, /Т / < 1,
1
Н, , _II л2ц ^2|
2ytAhy + т\\A2W -^WA^W /(1 - s2)max{0, + т(s - 1), y^l - тЛ ds > 0,
о ^ J
YTM - 1+ тWAW-т2\\A2\\ /(1 - s) max{0, x(2) + т(s - 1), y™ - т^ ds > 0,
0
(9)
где
(1)
1И2||+7Т(Л2)Н7 (1) _ ||A2||+7^42M
||A3|| ' У™ ||A3JJ (2) _ ||A||-7rAg7 л (2) _ ||A||-7^4M
||A2|| > Ь'7,M ||A2||
Доказательство. Рассмотрим произвольное допустимое запаздывание
т
нению условий теоремы 1 при N = 3. Предположим от противного, что существуют г ^ 0 и ф € Хг,3 такие, что
0
(0)A У y(s) ds> 0,
-T(t)
т(г) > 0
2рт (0)АН р(0) + т (г)||А21| —
1
— т2(г)|А3\у1 (1 — в2) шах{0, х(1 + т(г)(в — 1), у™ — т(г)*} ¿в> 0,
о
а следовательно, третье неравенство в системе (9) выполнено для т = т(г), 7 = р(0), М = Р(—т (г))
Так как р е Х4,3, по лемме 2 получаем четвертое неравенство в системе (9) при тех же т, 7 и м- Первые два неравенства являются прямым следствием из определения множества
□
Покажем на примерах, как условия полученной теоремы могут быть проверены на практике.
6. Скалярное уравнение. Рассмотрим простейший вариант системы (3) — скалярное уравнение
то
х(г) = а J х(г + в) ¿в, (ю)
-т (4)
аеМ
тогда и только тогда, когда а е [—4, 0]. Применим к этому уравнению для сравнения полученное достаточное условие.
В скалярном случае система из теоремы 2 принимает вид
2а + та2 — т2 |а|3 Г(1 - в2)тах{о, Д- + т(5 - 1), - тЛ > 0,
о I |а| N ]
7/х — 1 + т\а\ — т2а? Г(1 — в) тах <0, —---\- т(а — 1), т—г — — — та > д,а > 0.
о | а| а | а| а
ат точно малым, а ^М близким к единице. При а = 0 система решений не имеет. Остается рассмотреть случай а < 0. Учитывая, что знак 7 не влияет на разрешимость системы, положим 7 =1. Получаем упрощенную систему:
М < 1,
2 + та + т2а2 Г(1 - я2) тах 1 0, Д- +г(й - 1), 1 М - та 1 ¿8 < 0, , .
о I |а| |а| ) (И)
1 — /х + та + т2а? Г(1 — а) тах < 0, — + т(а — 1), ^ — та > ¿а < 0.
о I |а| |а| )
Очевидно, что в рассматриваемом случае выражение в левой части второго неравенства всегда превосходит выражение в левой части третьего неравенства. А заметив, что левая часть второго неравенства не убывает с увеличением м5 придем к выводу, что система (11) эквивалентна в смысле разрешимости неравенству
2 + та + т2а2 J( 1 - в2) тах |о, р- + т(в - 1), -гв| ¿в < 0.
о
Упростим подынтегральное выражение за счет того, что третий член в скобках неположителен, и сделаем замену переменной интегрирования:
1
2 + та + т2а2 J в(2 - в) тах |о, ~ Тв| <
о
Если а ^ —2, то неравенство неразрешимо, так как все слагаемые в левой части неотрицательны. Предположим теперь, что а < —2, и рассмотрим следующие случаи: 1) т € (0, 2/|а|];2) т € (2/|а|, 1].
В первом случае получаем кубическое неравенство относительно т:
2 +та - ^т2а - ^т3а2 <
Простое исследование показывает, что при а < —2 левая часть монотонно убывает по т и достигает наименьшего значения 2/|а| в точке т = 2/|а|. Таким образом, неравенство не имеет решений.
В случае т е (2/|а|, 1] имеем неравенство
8 4 2 +та- — -—— <0, 3а 3та2
левая часть которого монотонно убывает на указанном интервале. Следовательно, в качестве «наихудшего» т нужно взять 1. После домножения на 3а2 приходим к неравенству
За3 + 6а2 — 8а — 4 < 0.
Если найти наименьший корень полинома 3а3 + 6а2 — 8а — 4, то это будет наимень-
а
Таким образом, заключаем, что а € [—2.785,0] решает обобщенную задачу Мышкиса.
Этот результат существенно слабее, чем точная оценка, полученная детальным исследованием поведения решений уравнения в работе [7], но она лучше, чем можно было определить, например, стандартным преобразованием уравнения (10) к виду
о г
х(Ь) = ат(£)ж(£) — а J J ж(£) ^ ¿я =
-Т(г) г+э
о г о
= ат(¿)ж(£) — а2 J J ж(£ + п) ¿п ¿я
-т(г) г+э -т(О
и применением классического метода Разумихина с квадратичной функцией. Легко проверить, что этот подход приведет к оценке а € [—2, 0].
7. Система из двух уравнений. Рассмотрим систему из двух уравнений, соответствующую скалярному уравнению с комплексным коэффициентом:
о
т = — а) / *(*+5)гЬ, (12)
—в а
-Т(г)
в котором а € М, в > 0. В этом случае
\\Ак\\ = \\А\\к = ^^+^)\ к = 0,1,2,..., 1Т Ан 7 = а, 1Т (А2 )н 7 = а2 — в2, 7Т А/л = ар + вя, 1Т АV = (а2 — в2)р + 2авд.
Здесь
т т I 0 1
Р = 1 Я = 7 ( _1 0 )
Заметим, что \\/\\2 = р2 + я2- Таким образом, мы исключили векторные переменные / и 7, заменив их скалярными переменными р и я-
р2 + я2 < 1,
^ + г||А||-г2||А||2/(1-52)тах|о, Ц- ■ ^
р — 1+ т\\А\\ —т2\А\\2 /(1 — в) тах(0,
о
+ т(в — 1), у{р/д — тв \ ¿в > 0
+ т (в — 1), уР2) — тв) ¿в > 0,
А где
(1) = (а2 - 132)р + 2а.13д (2) = _1__ар + /Зд
Ур1И|Г 1И||3 ' Ур||А|| ||А||2 '
Сразу можно заметить, что при а > 0 система разрешима, и далее рассматривать только случай а ^ 0.
Учитывая, что в > 0, при я ^ 0 с ростом |я| множество решений системы (13) монотонно сжимается. Потому достаточно рассмотреть я € [0,1). Можно произвести замену д = у7в? —р2, (I € [0,1). При этом с ростом (I значения у^ -и у^ -
убывают, а следовательно, множество решений системы (13) расширяется. Таким образом, достаточно рассмотреть случай в, = 1, д = у!— р2, отказавшись от первого неравенства в системе. Остается исследовать систему с двумя параметрами р € ( —1,1) и т € (0,1]:
^ + т||А||-т2||А||2/(1-52)тах{о, ^ + у^-г*} ^ > 0,
р-1+т\\А\\-т2\\А\\2}(1-8) шах {о, + ¿8 > 0,
в которой у^к) = у{к) г—-, к = 1,2. р,\/ 1-р2
Далее можно заметить, что при р ^ а/\\А\\ с увеличением значения р множество решений системы расширяется, поэтому есть смысл рассматривать толькор ^ а/\\А\\.
Дальнейшая аналитическая проверка условий затруднена, хотя присутствие в системе всего двух ограниченных параметров р и т позволяет использовать методы трехмерной визуализации при конкретных заданных значениях а и в-
Кроме того, можно несколько ослабить условия, т. е. построить более простую систему, содержащую все решения исходной. Например, можно положить у^2 = 0
т
т
р
Наиболее простые условия будем иметь, если положим у^1 = 0. Тогда система
т
1
^ + г|И||-г2||А||2|(1-52)тах|о, + ^ > 0.
о
Проверка этого условия не представляет сложности, так как, если вычислить интеграл явно, задача сведется к паре полиномиальных неравенств не выше третьей степени (одно для случая т > 2а2/||А||3, а другое для случая т < 2а2/||А||3).
В частности, при в = 0.5 параметр а такой, что матрица
а в —в а
решает обобщенную задачу Мышкиса для системы (12), может быть выбран из отрезка [—2.649, —0.134].
Полученный результат можно сравнить, например, с теоремой о равномерной экспоненциальной устойчивости из работы [15], где изучается скалярное уравнение с комплексным параметром. При в = 0.5 представленное там условие принимает вид а е [—1.866, —0.134]. Отметим, что при постоянном запаздывании, равном 1,
а е [—4.221, —0.131]
8. Заключение. Представлены достаточные условия разрешимости обобщенной задачи Мышкиса для системы с распределенным запаздыванием. Условия выражены в терминах пустоты множества решений некоторой системы неравенств. Для двух частных случаев произведена проверка этих условий.
В дальнейшем хотелось бы улучшить результат, так как, например, для скалярного уравнения полученная оценка области разрешимости обобщенной задачи Мышкиса достаточно далека от точной, тогда как для систем с сосредоточенным запаздыванием, которые были рассмотрены ранее, оценки близки, а в некоторых случаях даже совпадают. Другая проблема, которую предстоит решить, — автоматизация проверки найденных условий. Опять же, для систем с сосредоточенным запаздыванием удалось свести задачу к набору задач выпуклого программирования, а для системы с распределенным запаздыванием — пока нет.
Литература
1. Мышкис А. Д. О решениях линейных однородных дифференциальных уравнений первого порядка устойчивого типа с запаздывающим аргументом // Математический сборник. 1951. Вып. 28. № 3. С. 641-658.
2. Yorke J. A. Asymptotic stability for one dimensional differential-delay equations // Journal of Differential Equations. 1970. Vol. 7. P. 189-202.
3. Yoneyama T. The 3/2 stability theorem for one-dimensional delay-differential equations with unbounded delay // Journal of Mathematical Analysis and Applications. 1992. Vol. 165. P. 133-143.
4. Amemiya Т. On the delay-independent stability of a delayed differential equation of 1st order // Journal of Mathematical Analysis and Applications. 1989. Vol. 142. P. 13-25.
5. Малыгина В. В. Об устойчивости решений некоторых линейных дифференциальных уравнений с последействием // Известия высших учебных заведений. Математика. 1993. Вып. 5. С. 72-85.
6. Krisztin Т. On stability properties for one-dimensional functional differential equations // Funkcialaj Ekvacioj. 1991. Vol. 34. P. 241-256.
7. Малыгина В. В. О точных границах области устойчивости линейных дифференциальных уравнений с распределенным запаздыванием // Известия высших учебных заведений. Математика. 2008. Вып. 52. № 7. С. 15-23.
8. Egorov A. On the stability analysis of equations with bounded time-varying delay // 15th IFAC Workshop on Time Delay Systems. Sinaia, Romania, 2019. P. 85-90.
9. Egorov A. The generalized Myshkis problem for a linear time-delay system with time-varying delay // Stability and Control Processes (SCP 2020), Lecture Notes in Control and Information Sciences — Proceedings / eds.: N. Smirnov, A. Colovkina. Cham, Switzerland: Springer, 2022. P. 223-231.
10. Egorov A. Robust stability analysis of time-varying delay systems // 17th IFAC Workshop on Time Delay Systems. Montreal, Canada, 2022. P. 210-215.
11. Fridman E. An improved stabilization method for linear time-delay systems // IEEE Transactions on Automatic Control. 2002. Vol. 47. N 11. P. 1931-1937.
12. Wu M., He Y., She J.-H., Liu G.-P. Delay-dependent criteria for robust stability of time-varying delay systems // Automatica. 2004. Vol. 40. P. 1435-1439.
13. Zhang C.-K., He Y., Jiang L., Wu M., Zeng H.-B. Stability analysis of systems with time-varying delay via relaxed integral inequalities // Systems & Control Letters. 2016. Vol. 92. P. 52-61.
14. Kolmanovskii V., Myshkis A. Introduction to the theory and applications of functional differential equations. Dordrecht, Netherlands: Springer Science + Business Media, 1999. 648 p.
15. Сабатулина Т. Л. Об экспоненциальной устойчивости дифференциального уравнения с комплексным коэффициентом и переменным распределенным запаздыванием // Теория управления и математическое моделирование. Ижевск, 2022. С. 119-121.
Статья поступила в редакцию 20 марта 2023 г. Статья принята к печати 25 апреля 2023 г.
Контактная информация:
Егоров Алексей Валерьевич — канд. физ.-мат. наук, доц.; alexey.egorov@spbu.ru
Razumikhin approach in the generalized Myshkis problem for systems with distributed delay*
A. V. Egorov
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Egorov A. V. Razumikhin approach in the generalized Myshkis problem for systems with distributed delay. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Seienee. Control Processes, 2023, vol. 19, iss. 2, pp. 148-161. https://doi.org/10.21638/11701/spbul0.2023.202 (In Russian)
The paper gives sufficient conditions for the solvability of the generalized Myshkis problem for a system of equations with a distributed time-varying delay and a constant kernel. Conditions on the kernel which guarantee the uniform stability of the system for any admissible delay are obtained. The admissible delay in this paper is a piecewise continuous function bounded from above in magnitude and growth rate. The applicability of the obtained conditions is illustrated by two examples.
Keywords: time-delay system, stability, distributed delay, generalized Myshkis problem, Razumikhin approach.
* This research was supported by the Russian Science Foundation (project N 21-71-00052).
References
1. Myshkis A. D. O resheniiakh lineinykh odnorodnykh difFerentsial'nykh uravnenii pervogo poriadka ustoichivogo tipa s zapazdyvaiushchim argumentom [On solutions of linear homogeneous differential equations of the first order of stable type with a retarded argument]. Matematicheskii sbornik [.Mathematical Collection], 1951, iss. 28, no. 3, pp. 641-658. (In Russian)
2. Yorke J. A. Asymptotic stability for one dimensional differential-delay equations. Journal of Differential Equations, 1970, vol. 7, pp. 189-202.
3. Yoneyama T. The 3/2 stability theorem for one-dimensional delay-differential equations with unbounded delay. Journal of Mathematical Analysis and Applications, 1992, vol. 165, pp. 133-143.
4. Amemiya T. On the delay-independent stability of a delayed differential equation of 1st order. Journal of Mathematical Analysis and Applications, 1989, vol. 142, pp. 13-25.
5. Malygina V. V. Ob ustoichivosti reshenii nekotorykh lineinykh difFerentsial'nykh uravnenii s posledeistviem [On stability of solutions of some linear differential equations with aftereffect], hvestiia vysshikh uchebnykh zavedenii. Matematika [Russian Mathematics], 1993, iss. 5, pp. 72-85. (In Russian)
6. Krisztin T. On stability properties for one-dimensional functional differential equations. Funkcialaj Ekvacioj, 1991, vol. 34, pp. 241-256.
7. Malygina V. V. O tochnykh granitsakh oblasti ustoichivosti lineinykh difFerentsial'nykh uravnenii s raspredelennym zapazdyvaniem [On the exact boundaries of the stability domain of linear differential equations with distributed delay], hvestiia vysshikh uchebnykh zavedenii. Matematika [Russian Mathematics], 2008, iss. 52, no. 7, pp. 15-23. (In Russian)
8. Egorov A. On the stability analysis of equations with bounded time-varying delay. 15th IFAC Workshop on Time Delay Systems. Sinaia, Romania, 2019, pp. 85-90.
9. Egorov A. The generalized Myshkis problem for a linear time-delay system with time-varying delay. Stability and Control Processes (SCP 2020), Lecture Notes in Control and Information Sciences — Proceedings. Eds N. Smirnov, A. Golovkina. Cham, Switzerland, Springer Publ., 2022, pp. 223-231.
th
Delay System-s. Montreal, Canada, 2022, pp. 210-215.
11. Fridman E. An improved stabilization method for linear time-delay systems. IEEE Transactions on Automatic Control, 2002, vol. 47, no. 11, pp. 1931-1937.
12. Wu M., He Y., She J.-H., Liu G.-P. Delay-dependent criteria for robust stability of time-varying delay systems. Automatica, 2004, vol. 40, pp. 1435-1439.
13. Zhang C.-K., He Y., Jiang L., Wu M., Zeng H.-B. Stability analysis of systems with time-varying delay via relaxed integral inequalities. Systems & Control Letters, 2016, vol. 92, pp. 52-61.
14. Kolmanovskii V., Myshkis A. Introduction to the theory and applications of functional differential equations. Dordrecht, Netherlands, Springer Science + Business Media Publ., 1999, 648 p.
15. Sabatulina T. L. Ob eksponentsial'noi ustoichivosti differentsial'nogo uravneniia s kompleksnym koeffitsientom i peremennym raspredelennym zapazdyvaniem [On the exponential stability of a differential equation with complex coefficient and time-varying distributed delay]. Teoriia upravleniia i matemati-cheskoe modelirovanie [Theory of management and mathematical modeling]. Izhevsk, 2022, pp. 119-121. (In Russian)
Received: March 20, 2023.
Accepted: April 25, 2023.
Author's information:
Alexey V. Egorov — PhD in Physics and Mathematics, Associate Professor; alexey.egorov@spbu.ru