ON STABILITY WITH DERIVATIVE OF A CLASS OF DIFFERENTIAL EQUATIONS OF NEUTRAL TYPE
- Authors: Balandin A.S1, Malygina V.V1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 1 (2019)
- Pages: 22-50
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/amcs/article/view/2161
- DOI: https://doi.org/10.15593/2499-9873/2019.1.02
- Cite item
Abstract
We consider an important class of differential-difference equations of neutral type, for which we study the asymptotic properties of a solution. We present necessary and sufficient conditions of exponential stability. The conditions have a geometric form of a region in the parameter space. The behavior of a solution at the boundaries of the region is analyzed separately. The loss of stability on the boundary can occur in various ways, namely through the appearance of stationary solutions, and through the appearance of periodic modes. Along with the asymptotic properties of a solution, we study analogous properties of its derivative.
Full Text
Введение. Настоящая работа посвящена исследованию асимптотических свойств одного класса функционально-дифференциальных уравнений нейтрального типа: , (*) где - произвольные вещественные числа, , - заданная локально-суммируемая функция. Это уравнение возникает в прикладных задачах: в работе [1] с его помощью описывается динамика популяции клеток, в работе [2] - движение плоских упругих плит с учетом трения. Но не менее интересно оно с теоретической точки зрения. Уравнение (*) являет собой удачный пример уравнения нейтрального типа: с одной стороны - достаточно простой, чтобы получить эффективные признаки устойчивости тривиального решения, с другой - достаточно сложный, чтобы на нем проявилось все разнообразие асимптотических свойств решений. Это подтверждается значительным количеством чисто теоретических исследований (см., например, [3-10]), в которых изучается уравнение (*). Особое внимание в статье уделяется асимптотическим свойствам производной решения (*). Как во всяком дифференциальном уравнении, информацию о производной следует учитывать и правильно ее использовать. Для уравнений запаздывающего типа (таким будет уравнение (*) при ) асимптотические свойства производной легко выводятся из свойств самого решения. Но при уравнение (*) уже не является разрешенным относительно производной и для него вопрос об асимптотических свойствах производной, оставаясь содержательным, перестает быть тривиальным. 1. Формулы представления решения и его производной. Пусть - множество натуральных чисел, , - множество комплексных чисел, - пространство суммируемых на множестве функций с нормой ( ), - пространство измеримых и ограниченных в существенном на множестве функций с нормой . Символом будем обозначать единичный (тождественный) оператор. Введем в рассмотрение оператор сдвига, действующий в пространстве измеримых (кусочно-непрерывных, суммируемых) на каждом конечном отрезке функций: Легко видеть, что в этих пространствах при любом оператор ограниченно обратим. Изменяя масштаб времени и коэффициенты , , перепишем уравнение (*) в эквивалентной форме, более удобной для дальнейшего изучения: . (1) Как известно [12, с. 35], уравнение (1) с заданным начальным условием однозначно разрешимо в классе локально абсолютно непрерывных функций и его решение имеет вид ([12, с. 84], [6]): , . (2) Функция называется фундаментальным решением, функция - функцией Коши. Представление (2) принято называть формулой Коши. Фундаментальное решение является локально абсолютно непрерывной функцией и однозначно определяется как решение уравнения (1) при с начальным условием . На каждом интервале , , функция абсолютно непрерывна, но имеет разрывы первого рода во всех точках , причем [6]. Как показано в [6], [9], однозначно определяется через равенством , . (3) Получим еще одну полезную формулу, которая связывает функции и . Лемма 1. Для фундаментального решения и функции Коши уравнения (1) справедливо соотношение , . (4) Доказательство. В силу определения фундаментального решения при любом справедливо равенство . Применяем формулу (4) и учитываем перестановочность операторов и : . Следовательно, , что, в силу отмеченной выше обратимости , эквивалентно (4). ▲ Наряду с и введем функцию , которая (однозначно) определяется на по следующему правилу: , . (5) Имея в виду изучение асимптотического поведения решения уравнения (1) вместе с производной, было бы полезно располагать для формулой, аналогичной формуле (2). Заметим, что прямое дифференцирование равенства (2) невозможно, так как функция разрывна и, следовательно, не имеет производной в классическом смысле. Тем не менее удается доказать следующее утверждение. Теорема 1. Пусть - решение уравнения (1). Тогда при любом справедливо представление . (6) Доказательство. Пусть . Используя определение оператора , равенства (5), (4), (2) и свойства фундаментального решения, имеем: . При заданных начальных условиях и определяются уравнением (1) однозначно, следовательно, справедливость представления (6) доказана. ▲ Замечание 1. Если допустить дифференцирование разрывных функций за счет выхода в класс обобщенных функций [13, гл. IV, § 4], то для производной функции Коши получаем . Теперь, если продифференцировать равенство (2) в обобщенном смысле и подставить вместо найденное представление, мы получим другое доказательство формулы (6). Обозначим характеристическую функцию уравнения (1). Из определения функций , , , следует, что они непрерывны или кусочно-непрерывны на . При отрицательных значениях аргументов не нарушая общности их можно доопределить нулем. В работе [9] доказано, что функции и при растут не быстрее экспоненты. Из равенств (4) и (5) следует, что тем же свойством обладают функции и . Следовательно, ко всем четырем функциям применимо преобразование Лапласа: , , (7) , . Приведем в удобной для нас формулировке одно хорошо известное утверждение. Лемма 2 [14, с. 497-498]: Следствие 1. При любом Лемма 3. Пусть . Тогда при любом , где , имеет место равенство Доказательство. В условиях леммы , значит . При в силу следствия 1 при получаем: Аналогично, при имеем , что и требовалось. ▲ 2. Экспоненциально устойчивые решения. Начнем исследование асимптотических свойств решений уравнения (1) с самого сильного типа устойчивости - экспоненциального. Найдем условия на параметры уравнения, при которых введенные выше функции - Y, Z - имеют экспоненциальные оценки. Теорема 2. Пусть и все корни характеристической функции лежат слева от мнимой оси. Тогда существуют и такие, что при всех справедливы оценки: а) , b) , c) , d) . Доказательство. Поскольку , при функция не обращается в нуль. Следовательно, характеристическая функция не имеет нулей в области при достаточно большом . Поскольку - аналитическая функция, в прямоугольнике она может иметь лишь конечное число нулей. Значит, среди них есть корень с наибольшей (отрицательной, в силу условий теоремы) вещественной частью. Следовательно, существует такое число : , что в полуплоскости функция не обращается в нуль. Из (7) заключаем, что при функция является аналитической. По формуле обратного преобразования Лапласа и лемме 2 имеем: (8) откуда следует утверждение а. Для оценки применим формулу (3). Так как , оператор обратим в пространстве функций, ограниченных на полуоси [13, с. 230], причем , следовательно, , и, с учетом установленной выше оценки на Х, а также следующего из определения неравенства , получаем утверждение b: Оценка на следует из формулы (4): . Наконец, оценка на следует из формулы (5): ▲ Следствие 2. Пусть . Тогда из оценки а следуют оценки b, c и d. Следствие 3. Пусть . Тогда оценка а эквивалентна оценке b, а оценка c - оценке d. Остановимся подробнее на одной важной постановке начальной задачи для уравнения (*), отличной от принятой в разделе 1. Будем предполагать, что при отрицательных значениях аргумента решение уравнения доопределено начальными функциями, заданными на отрезке : (9) Если понимать под решением уравнения (9) непрерывное продолжение начальной функции [15-17], то функцию естественно считать частью решения. Тогда логично требовать, чтобы была абсолютно непрерывна, а суммируема на и выполнялись равенства , . С другой стороны, можно не рассматривать начальные функции как часть решения, а считать, что мы доопределяем функции и при отрицательных значениях аргумента некоторыми суммируемыми на функциями соответственно. Здесь равенство уже не обязательно. Отметим, что начальная задача для уравнения (1) во многих работах ставится в виде (9), а в определения устойчивости решения включаются начальные функции [15-18]. Поэтому имеет смысл получить утверждения об асимптотических свойствах решения и его производной с явными оценками начальных функций. Независимо от понимания решения, задачу (9) можно переписать в виде (1) , (10) где роль внешнего возмущения играет функция , определяемая по правилу В такой постановке (9) оказывается частным случаем (1) со специальной правой частью . В любом случае , но если окажется, что при , то, учитывая это, можно получать более тонкие оценки. Теорема 3. Пусть выполнены условия теоремы 2, а . Тогда найдутся такие , что для решения уравнения (10) справедливы оценки: 1. Если , то , . 2. Если , то , . 3. Если , то . Доказательство. Для первое утверждение теоремы вытекает из формулы (2) и пп. a, b теоремы 2: . Для применяем неравенство Гельдера, где : . Докажем второе утверждение. Пусть , . Тогда по формуле (6) имеем . Далее интегрируем это равенство на отрезке и применяем пп. c, d теоремы 2. В итоге при получаем . При применяем неравенство Гельдера, где : Наконец, при имеем ▲ Следствие 4. Пусть выполнены условия теоремы 2, а . Тогда тривиальное решение уравнения (10) экспоненциально устойчиво вместе с производной. Для уравнения (1) можно указать эффективно проверяемые необходимые и достаточные условия экспоненциальной устойчивости решения. Эти условия - в аналитической и геометрической форме - были найдены в работах [4] и [8]. Зададим в пространстве криволинейную поверхность и три плоскости: , , . Пересечение Г и Р обозначим . Поверхности ограничивают в пространстве область D. Из результатов работы [8] следует, что условия теоремы 2 выполнены тогда и только тогда, когда , т.е. область D является областью экспоненциальной устойчивости уравнения (2). Теоремы 2 и 3 усиливают этот результат: в области D имеет экспоненциальную оценку не только решение уравнения (1), но и его производная. Интересен вопрос о том, как меняется асимптотика решения рассматриваемого уравнения при переходе через границу области D. В следующих двух разделах мы изучим поведение решения на плоскости Р, поверхности Г и «ребре» . Поведение решения на плоскостях , наиболее сложно и требует отдельного исследования. 3. Стабилизирующиеся решения. Изучим асимптотическое поведение решения уравнения (1) и его производной в случае, когда неравенство остается справедливым, а точка лежит на плоскости . Теорема 4. Пусть и . Тогда существуют такие , что при всех справедливы оценки: а) , b) , c) , d) . Доказательство. Характеристическая функция уравнения (1) в этом случае имеет вид . Разделяя в уравнении вещественную и мнимую часть, получаем , т.е. на мнимой оси функция имеет единственный корень . Так как , то корень - простой. Применяя теорему Руше [14, с. 454], убеждаемся, что справа от мнимой оси у функции нулей нет. Аналогично тому, как это было доказано в теореме 2, устанавливаем, что при некотором в области функция не обращается в нуль. Следовательно, существует такое число : , что в полуплоскости функция не имеет иных нулей, кроме единственного простого корня . Зафиксируем , и рассмотрим в плоскости прямоугольник , вершины которого находятся в точках , , , , с ориентацией границы против часовой стрелки. Разобьем интеграл от функции по границе этого прямоугольника на четыре интеграла: , (11) и оценим каждый из них. По формуле обратного преобразования Лапласа Применяя установленную при доказательстве теоремы 2 оценку (8), получаем: (12) . Далее, (13) . Аналогично, (14) . С другой стороны, по теореме Коши о вычетах [14, c. 85], . Переходя в равенстве (11) к пределу при , получаем утверждение а. Для оценки применим формулу (4). Поскольку , то оператор обратим, , следовательно, , С учетом установленной выше оценки на и следующего из определения неравенства , получаем утверждение b: Оценка на следует из формулы (4) и утверждения b: И наконец, оценка на вытекает из следствия 2 и утверждения с. ▲ Соответствующий результат получается и для уравнения (10). Теорема 5. Пусть выполнены условия теоремы 4. Тогда при любом найдутся такие , что при любом для решения уравнения (10) справедливы оценки: 1. Если , то . 2. Если , то . 3. Если , то . Доказательство. Первое утверждение следует из формулы (2) и пп. a, b теоремы 4. Доказательство второго и третьего утверждений дословно повторяет рассуждения теоремы 3. ▲ Следствие 5. В условиях теоремы 4 тривиальное решение уравнения (10) является равномерно устойчивым, но не является асимптотически устойчивым. Следствие 6. Пусть выполнены условия теоремы 4, а . Тогда производная тривиального решения уравнения (10) экспоненциально устойчива. Теорема 5 описывает класс уравнений вида (10), для которых решение не является асимптотически (и тем более экспоненциально) устойчивым. При переходе через плоскость уравнение теряет устойчивость за счет появления решений, имеющих на бесконечности ненулевой предел. Однако производная любого решения имеет экспоненциальную оценку: интегральную при и равномерную при . 4. Асимптотически периодические решения. В этом разделе мы рассмотрим асимптотическое поведение решения уравнения (1) и его производной в случае, когда неравенство выполнено, а точка принадлежит поверхности Г. Теорема 6. Пусть и . Тогда существуют и такие постоянные что при всех справедливы оценки: а) , b) , c) , d) , причем во всех оценках а-d . Доказательство. В условиях теоремы 6 характеристическая функция уравнения (1) имеет на мнимой оси два простых комплексно-сопряженных корня , где . Справа от мнимой оси функция не имеет нулей. Аналогично тому, как это было доказано в теореме 2, устанавливаем, что при некотором в области функция не обращается в нуль. Следовательно, существует такое число : , что в полуплоскости функция не имеет иных корней, кроме пары чисто мнимых . Зафиксируем , и рассмотрим в плоскости прямоугольник , вершины которого находятся в точках , , , , с ориентацией границы против часовой стрелки. Рассмотрим интеграл от функции по контуру и представим его в виде (11). Заметим, что полученные при доказательстве теоремы 4 оценки (12)-(14) на интегралы остаются в силе. По теореме Коши о вычетах [14, c. 85] имеем где . Переходя в равенстве (11) к пределу при и используя формулу обратного преобразования Лапласа, получаем . (15) Из формулы (3) имеем Заметим, что где . Следовательно, . (16) Далее, легко видеть, что где . С учетом формулы (4) имеем: . (17) Наконец, полагая , на основе формулы (5) устанавливаем оценку . (18) Теперь, чтобы получить все четыре утверждения теоремы, достаточно обозначить , , и воспользоваться оценками (15)-(18). ▲ Соответствующий результат получается также для уравнения (10). Теорема 7. Пусть выполнены условия теоремы 6. Тогда при любом найдутся и такие , что при всех для решения уравнения (10) справедливы оценки: 1. Если , то . 2. Если , то . 3. Если , то . Доказательство проводится по той же схеме, что и теорема 3: для доказательства первого утверждения используется формула (2) и пп. a, b теоремы 6, для доказательства второго и третьего утверждений - формула (6) и пп. c, d теоремы 6. ▲ Следствие 7. В условиях теоремы 6 тривиальное решение уравнения (10) является равномерно устойчивым, но не является асимптотически устойчивым. Следствие 8. Пусть выполнены условия теоремы 6, а . Тогда тривиальное решение уравнения (10) равномерно устойчиво вместе с производной. Теорема 7 описывает класс уравнений вида (1), решения которых асимптотически сближаются с периодическими функциями. Аналогичным свойством обладают и производные, если оценивать расстояние до периодической функции по норме пространства . Следовательно, при переходе через поверхность Г уравнение (1) теряет устойчивость за счет появления периодических решений. 5. Неустойчивые решения с устойчивой производной. Осталось рассмотреть поведение решения уравнения (2) в случае, когда точка принадлежит . Это означает, что коэффициенты уравнения (1) связаны соотношениями: , причем , а задача (9) принимает вид (19) Покажем, что уравнение (20) не является равномерно (и тем более асимптотически) устойчивым. Выберем в качестве начальных функций , . Подставляя в (19) , убеждаемся, что эта функция является решением. Следовательно, уравнение (19) имеет неограниченные решения. Исследуем для уравнения (19) поведение функций и . Теорема 8. Пусть и . Тогда существуют такие , что при всех справедливы оценки: а) , b) . Доказательство. Характеристическая функция уравнения (1) в этом случае имеет вид . Разделяя в уравнении вещественную и мнимую часть, получаем , т.е. на мнимой оси функция имеет единственный корень . Так как , а , то корень имеет кратность 2. Справа от мнимой оси у функции нулей нет. Аналогично тому, как это было доказано в теореме 2, устанавливаем, что при некотором в области функция не обращается в нуль. Следовательно, существует такое число : , что в полуплоскости функция не имеет иных нулей, кроме двукратного корня . Зафиксируем , и рассмотрим в плоскости прямоугольник , вершины которого находятся в точках , , , , с ориентацией границы против часовой стрелки. Разобьем интеграл от функции по границе этого прямоугольника на четыре интеграла: , (20) и оценим каждый из них. По формуле обратного преобразования Лапласа Пусть . По лемме 3 имеем , откуда следует Далее, . Аналогично, . С другой стороны, по теореме Коши о вычетах, Переходя в равенстве (20) к пределу при , с учетом выбора получаем оценку: . (21) Далее, применяя формулу (5), имеем: откуда следует утверждение а. Утверждение b получается из формулы (5) и оценки а. ▲ Для производной любого решения уравнения (19) имеет место следующее утверждение. Теорема 9. Пусть выполнены условия теоремы 8, а . Тогда найдутся такие , что при всех для производной решения уравнения (19) справедливы оценки: 1. Если , то . 2. Если , то . Доказательство проводится по той же схеме, что и теорема 3: применяем формулу (6) и пп. а, b теоремы 8. ▲ Следствие 9. Пусть выполнены условия теоремы 8, а . Тогда производная тривиального решения уравнения (19) равномерно устойчива. Замечание 2. В процессе доказательства теоремы 8 была установлена формула (21), которая весьма точно характеризует асимптотическое поведение функции Коши. Используя равенство (3), из оценки (21) легко получить формулу асимптотического поведения фундаментального решения. , следовательно, . (2) Таким образом, несмотря на то что уравнение (19) не является устойчивым, с помощью оценок (21) и (22) можно дать хорошее описание асимптотического поведения любого его решения. Замечания и комментарии. Идея изучения асимптотического поведения решений уравнений нейтрального типа вместе с производной выдвигалась и развивалась многими авторами [17, 19-21], но конкретное воплощение этой идеи осуществлялась ими по-разному. Чаще всего реализовался подход, при котором в пространстве решений априорно выбиралась норма, учитывающая и решение, и ее производную. Например, в монографии [17, c. 161] для уравнений нейтрального типа вводится понятие -устойчивости, суть которого в том, что классические определения устойчивости (по Ляпунову, асимптотической и т.д.) применяются к сумме . В задачах вида (9) начальная функция также предполагается непрерывно дифференцируемой (причем , а ). Такой подход позволяет доказать утверждения типа следствий 4 и 8 (ситуация одинакового асимптотического поведения и ), но непригоден для случаев, когда их асимптотика различна. Кроме того, из формулы (6) следует, что в рамках -устойчивости невозможно расширить класс начальных функций до при . Другой пример - работа [20], где норма решения задается равенством результаты формулируются в терминах оценок вида , а и определяются корнями характеристической функции. Сопоставим этот подход с теоремами 3, 5, 7 и 9. Как отмечалось выше, выбор интегральной нормы для функции оправдан и задается пространством, которому принадлежат начальные функции. В [20] предполагается, что , поэтому квадратичная норма для производной выбрана правильно. Но для самого решения такая норма невыгодна, так как для функции можно получать более тонкие поточечные оценки. Другой недостаток выбора общей нормы для и заключен в том, что полученная оценка неизбежно будет диктоваться той функцией, которая растет на бесконечности быстрее, и мы утратим информацию об асимптотическом поведении второй функции. Между тем результаты разделов 4 и 6 этой статьи показывают, что и могут вести себя существенно по-разному: производная стремится к нулю, а решение - нет; производная ограничена, а решение неограниченно возрастает.About the authors
A. S Balandin
Perm National Research Polytechnic University
V. V Malygina
Perm National Research Polytechnic University
References
- Diekmann O., Getto P., Nakata Y. On the characteristic equation and its use in the context of a cell population model // Journal of Mathematical Biology. - 2016. - Vol. 72. - P. 877-908.
- Putelat T., Willis J.R., Dawes J.H.P. Wave-modulated orbits in rate-and-state friction // International Journal of Non-Linear Mechanics. - 2012. - Vol. 47. - P. 258-267.
- Hahn W. Zur Stabilität der Lösungen von linearen Differetial-Differetzen-gleichungen mit konstanten Koeffitienten // Annals of Mathematics. - 1956. - Vol. 131. - S. 151-166.
- Ожиганова И.А. Определение области асимптотической устойчивости для дифференциальных уравнений первого порядка с отклоняющимся аргументом // Труды УДН. - 1962. - Т. 1. - С. 52-62.
- Громова П.С., Зверкин А.М. О тригонометрических рядах, суммы которых являются непрерывными неограниченными функциями на вещественной оси - решениях уравнений с запаздывающим аргументом // Дифференциальные уравнения. - 1968. - № 4. - C. 1774-1784.
- Баландин А.С., Малыгина В.В. Об экспоненциальной устойчивости линейных дифференциально-разностных уравнений нейтрального типа // Известия вузов. Математика. - 2007. - № 7. - С.17-27.
- Junca S., Lombard B. Stability of a critical nonlinear neutral delay differential equation // Journal of Differential Equations. - 2014. - Vol. 256, iss. 7. - Р. 2368-2391.
- Баландин А.С. Об устойчивости одного дифференциального уравнения, не разрешенного относительно производной // Современные методы прикладной математики, теории управления и компьютерных технологий (ПМТУКТ-2017): сб. тр. X Междунар. конф. - Воронеж, 2017. - С. 68-71.
- Баландин А.С. О связи между фундаментальным решением и функцией Коши для функционально-дифференциальных уравнений нейтрального типа // Прикладная математика и вопросы управления. - 2018. - № 1. - С. 13-25.
- Баландин А.С. О разрешимости некоторых функциональных уравнений // Современные методы прикладной математики, теории управления и компьютерных технологий (ПМТУКТ-2018): сб. тр. XI Междунар. конф. - Воронеж, 2018. - С. 60-63.
- Баландин А.С. О разрешимости одного функционального уравнения // Динамические системы в науке и технологиях» (DSST-2018): тез. докладов междунар. конф. - Алушта, 2018. - С. 13-15.
- Азбелев Н.В., Максимов В.П., Рахматуллина Л.Ф. Введение в теорию функционально-дифференциальных уравнений. - М.: Наука, 1991. - 180 с.
- Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. - М.: Наука, 1981. - 544 с.
- Лаврентьев М.А., Шабат Б.В. Методы теории функции комплексного переменного. - М.: Наука, 1987. - 688 с.
- Беллман Р., Кук К. Дифференциально-разностные уравнения. - М.: Мир, 1963. - 548 с.
- Эльсгольц Л.Э., Норкин С.Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. - М.: Наука, 1971. - 296 с.
- Колмановский В.Б., Носов В.Р. Устойчивость и периодические режимы регулируемых систем с последействием. - М.: Наука, 1981. - 484 с.
- Хейл Дж. Теория функционально-дифференциальных уравнений. - М.: Мир, 1984. - 421 с.
- Азбелев Н.В., Симонов П.М. Устойчивость решений уравнений с обыкновенными производными. - Пермь: Изд-во Перм. ун-та, 2001. - 230 с.
- Власов В.В. Спектральные задачи, возникающие в теории дифференциальных уравнений с запаздыванием // Современная математика. Фундаментальные направления. - 2003. - Т. 1. - C. 69-83.
- Курбатов В.Г. Линейные дифференциально-разностные уравнения. - Воронеж: Изд-во Воронеж. гос. ун-та, 1990. - 168 с.
Statistics
Views
Abstract - 66
PDF (Russian) - 46
Refbacks
- There are currently no refbacks.