LINEAR REPRESENTATION OF THE STABILITY AREAS OF SOME CLASSES OF DIFFERENCE EQUATIONS
- Authors: Kandakov A.A1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 3 (2019)
- Pages: 40-53
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/amcs/article/view/2145
- DOI: https://doi.org/10.15593/2499-9873/2019.3.02
- Cite item
Abstract
The work is devoted to obtaining new stability domains for linear autonomous difference equations, as well as a demonstration of the use of the lined-structure property of D -decompositions in solving this problem. It is proposed to study D -decompositions in the form of a set of plane cross-sections consisting of lines. Using the proposed method, one of the known results is obtained and a new area is found. The prospects of transferring the approach to the study of wider classes of equations are indicated.
Full Text
Введение Рассмотрим линейное автономное разностное уравнение произвольного порядка на дискретной полуоси (1) где Определение 1. Решением уравнения (1) назовем функцию удовлетворяющую равенству (1) для всех . Одним из фундаментальных свойств разностных уравнений является устойчивость решений. Оно отражает непрерывность зависимости решения от начальной функции. В данной работе изучается асимптотическая устойчивость уравнения (1), которая для автономных уравнений совпадает с экспоненциальной исходя из формулы представления решения [1]. Для уравнения (1), в силу его линейности и однородности, устойчивость определяется оценкой решения, а наличие одного такого решения дает гарантию устойчивости остальных, поэтому говорят об устойчивости всех решений как об устойчивости уравнения. Определение 2. Уравнение (1) будем называть экспоненциально устойчивым, если существует такая константа , что для каждого решения при некотором для всех имеем . Теорема 1 [1, с. 77]. Уравнение (1) экспоненциально устойчиво, если и только если все корни его характеристического многочлена принадлежат единичному кругу . При изучении устойчивости автономных уравнений целесообразно стремиться получить признаки устойчивости, выраженные в терминах параметров уравнения и соответствующие области устойчивости в пространстве параметров. Определение 3. Областью устойчивости называется множество точек пространства параметров, соответствующих устойчивым уравнениям. Для решения такого рода задач Ю.И. Неймарк [2] на основе идей И.А. Вышнеградского разработал метод D-разбиений. Суть метода заключается в построении границ в пространстве параметров, при переходе через которые изменяется количество корней характеристического уравнения, находящихся внутри заданной области комплексной плоскости. Областью устойчивости является объединение всех областей разбиения, которым соответствует нулевое число корней характеристического уравнения в единичном круге. Области устойчивости некоторых классов уравнений могут иметь свойства, которые облегчают исследование. Таким свойством является линейчатость. Определение 4. Назовем область линейчатой или обладающей линейчатой структурой, если ее граница является объединением отрезков, полупрямых и прямых. Замечание 1. Если D-разбиение обладает линейчатой структурой, то любая часть пространства, ограниченная поверхностями D-разбиения, тоже линейчата. Поэтому и область устойчивости оказывается линейчатой. Очевидно, что исследовать область устойчивости, обладающую линейчатой структурой, намного удобнее в силу ее простоты. Конечно, этим свойством обладают области устойчивости далеко не всяких уравнений. Однако это встречается чаще, чем кажется на первый взгляд. Подтверждение можно найти, в частности, в работах М.В. Мулюкова [3, 4]. В своей диссертации он эффектно и эффективно использует линейчатую структуру D-разбиений для получения признаков устойчивости квазиполиномов и приводит несколько классов дифференциальных уравнений с последействием, обладающих этим свойством. 1. Об одной известной области Рассмотрим уравнение третьего порядка (2) Область устойчивости уравнения (2) получена, например, в [5, 6]. Теорема 2 [5]. Уравнение (2) экспоненциально устойчиво, если и только если выполнены неравенства Представленная на рис. 1 область является линейчатой. Покажем это, получив ту же область другим способом. Характеристический полином уравнения (2) имеет вид (3) Воспользуемся методом D-разбиения. Положим в (3) : Рис. 1. Область устойчивости уравнения (2) домножим обе части полученного равенства на : затем сделаем замену , : В последнем уравнении разделим действительную и мнимую части: В случае получаем прямые для остальных ситуаций можно делить второе уравнение системы на Обозначим . В итоге поверхности D-разбиения задаются уравнениями (4) Обратим внимание, что все сечения состоят из набора прямых. Таким образом, доказана следующая Теорема 3. D-разбиение и область устойчивости уравнения (2) имеют линейчатую структуру. Теперь легко строить сечения, фиксируя Для этого строим прямые подставляем значение в уравнение и строим третью прямую (рис. 2). Рис. 2. Сечения D-разбиения при разных Замечание 2. Осуществленная ранее замена соответствует повороту системы относительно оси ξ. Поэтому можно с уверенностью сказать, где находится область устойчивости. На рис. 2 она попадает в треугольник, образованный в результате пересечения прямых. Соберем сечения при всевозможных «в стопку» и сформируем область устойчивости уравнения (2) (рис. 3). Повернув изображение (рис. 4), получим область, описанную в теореме 2. Теорема 4. Уравнение (2) экспоненциально устойчиво, если и только если точка попадает в область пространства содержащую точку и ограниченную поверхностями (4). Рис. 3. Набор сечений области устойчивости уравнения (2) Рис. 4. Область устойчивости уравнения (2) в линейчатом виде 2. Область устойчивости уравнения четвертого порядка Разовьем предложенный подход: будем рассуждать с точки зрения линейчатости D-разбиения, построенного для получения новой области устойчивости. Увеличив все запаздывания аргумента в уравнении (2) на единицу, получим уравнение четвертого порядка с тремя запаздываниями: (5) Будем действовать аналогично изложенному выше. Возьмем характеристическое уравнение и положим (6) Домножаем на и делаем замену , : Разделяем мнимую и действительную части: В случае получаем прямые для остальных ситуаций разделим обе части второго уравнения на Имеем: Преобразуем и обозначим . В итоге можно представить D-разбиение в следующем виде: (7) Так же, как для уравнения (2), есть возможность построить D-разбиение пространства параметров уравнения (5) по сечениям при фиксированных . Формулы (7) показывают, что любое такое сечение будет состоять из набора прямых, т.е. справедливо следующее утверждение. Теорема 5. D-разбиение и область устойчивости уравнения (5) имеют линейчатую структуру. Сечения формируются из прямых и прямых вида где в роли выступают корни уравнения . При этом количество корней меняется при различных , а следовательно, меняется количество прямых (рис. 5). г в б а Рис. 5. График уравнения и вид сечений D-разбиения: a - при б - ; в - ; г - Собираем сечения при различных «в стопку» и получаем объемный вид D-разбиения (рис. 6). Рис. 6. Объемное представление D-разбиения Теорема 6. Уравнение (5) экспоненциально устойчиво, если и только если точка попадает в область пространства содержащую точку и ограниченную поверхностями (7). Доказательство. Достаточно выяснить, какие из областей D-разбиения составляют область устойчивости. В этом поможет линейчатость D-разбиения. Корни характеристического уравнения (6) непрерывно зависят от параметров a, b, c. Пусть - корень характеристического уравнения. Если принадлежит поверхностям (7), то . Исследуем, как изменяется модуль корня при пересечении поверхностей. Полагая в (6) , получим множества точек пространства параметров, которым соответствуют характеристические уравнения, имеющие корень с модулем . Производя замены и получим D-разбиение в виде (8) Поверхности (7) являются частным случаем поверхностей (8) при . Зададимся целью определить знак производных от по направлениям пересечения поверхностей (8). Для границ, не зависящих от или s, удобно рассматривать изменение модуля при движении вдоль оси Ob. Для этого достаточно знать знак производных для связанных неявно выражениями (8). Он виден из производных, найденных явно: Для прямых (8), заданных параметрически, поступим иным образом. Найдем значение производной не для всех точек поверхностей, а только для пересечения с плоскостью Ниже покажем, что этого достаточно. Для случая можно явно вычислить корни и в дальнейшем найти производную для соответствующих прямых: Оценим динамику корней характеристического многочлена, попадающих в единичный круг. Если производная уравнения границы положительна в точке , то количество корней, попавших в единичный круг при пересечении границы в направлении Ob, увеличивается (так как рассмотрена обратная зависимость ), соответственно, если производная отрицательна, количество корней уменьшается. Благодаря полученной информации мы можем определить область, в которой все корни характеристического многочлена попадают в единичный круг, она и будет областью устойчивости (согласно теореме 1). Все области D-разбиения рассекаются плоскостью На рис. 7 сечение области устойчивости - окрашенный четырехугольник, а неокрашенные области содержат как минимум один корень характеристического уравнения, не попадающий в единичный круг. Рис. 7. Динамика корней характеристического полинома (6) при Исследуем область пространства, сечением которого является окрашенный четырехугольник. Очевидно, что уравнение, соответствующее точке (0, 0, 0), устойчиво. Следовательно, эта область пространства и является областью экспоненциальной устойчивости уравнения (5). Теорема доказана. Рис. 8. Область устойчивости уравнения (5) по сечениям Полученная область представляет собой своеобразный квазиконус, срезанный углом, однако при этом она обладает линейчатой структурой. Вид области с различных ракурсов представлен на рис. 8 (по сечениям) и рис. 9. Рис. 9. Область устойчивости уравнения (4) Замечание 3. Следует указать, что , - многочлены Чебышёва первого и второго рода. Есть предположение, что многочлены этих классов дают описание областей устойчивости и для уравнений более высоких порядков. Это может существенно упростить решение задач устойчивости. Выводы В работе рассмотрены некоторые разностные уравнения вида (1), при этом использовано свойство линейчатости D-разбиений пространства параметров. Представлен способ исследования областей устойчивости в виде набора сечений. Итогом работы стало получение области устойчивости уравнения с тремя запаздываниями (4). Основываясь на полученном опыте, можно выдвинуть гипотезу о возможности получения областей устойчивости уравнений более высоких или даже произвольных порядков аналогичными способами. Автор выражает свою признательность участникам пермского семинара по функционально-дифференциальным и разностным уравнениям за активный интерес к работе, советы и обсуждения. Отдельная благодарность научному руководителю, старшему научному сотруднику К.М. Чудинову, а также ведущему научному сотруднику В.В. Малыгиной за постановку задачи и всестороннюю поддержку в процессе выполнения работы.About the authors
A. A Kandakov
Perm National Research Polytechnic University
References
- Elaydi S. An Introduction to Difference Equations. - N.Y.: Springer, 2005. - 539 p.
- Неймарк Ю.И. Устойчивость линеаризованных систем. - Л.: Изд-во ЛКВВИА, 1949. - 140 с.
- Мулюков М.В. Устойчивость двупараметрических систем линейных автономных дифференциальных уравнений с ограниченным запаздыванием // Изв. Ин-та матем. и информ. - 2018. - Т. 51. - С. 79-122.
- Мулюков М.В. Устойчивость систем линейных автономных дифференциальных уравнений с ограниченным запаздыванием: дис.. канд. физ.-мат. наук. - Пермь, 2018. - 150 с.
- Баландин А.С., Малыгина В.В. О разрешимости одного класса разностных уравнений // Вычислительная механика: сб. науч. тр. - Пермь, 2006. - № 4. - С. 67-72.
- Нигматулин Р.М., Кипнис М.М. Свойства дискретных систем третьего порядка на границе их области устойчивости // Фундаментальные исследования. - 2015. - № 9. - С. 39-43.
Statistics
Views
Abstract - 64
PDF (Russian) - 39
Refbacks
- There are currently no refbacks.