Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

СВОЙСТВА ДИСКРЕТНЫХ СИСТЕМ ТРЕТЬЕГО ПОРЯДКА НА ГРАНИЦЕ ИХ ОБЛАСТЕЙ УСТОЙЧИВОСТИ

Нигматулин Р.М. 1 Кипнис М.М. 1
1 ФГБОУ ВПО «Челябинский государственный педагогический университет»
Рассматриваются дискретные системы, описываемые линейными разностными уравнениями третьего порядка с действительными коэффициентами. Исследована граница области асимптотической устойчивости нулевого решения в пространстве коэффициентов уравнения. Полностью описаны свойства корней их характеристических уравнений на границе области устойчивости в трехмерном пространстве параметров. Для всех точек на границе указаны расположения каждого из трех корней относительно единичного круга в комплексной плоскости. Выписаны общие решения разностных уравнений для случаев, когда параметры находятся на указанной границе. Сделаны выводы об асимптотических свойствах траекторий этих систем посредством разделения точек границы на устойчивые и неустойчивые. Указано также разделение точек границы по признаку наличия – отсутствия колебательных (периодических, псевдошумовых) решений. Ставится задача поиска областей частичной устойчивости в пространстве начальных значений.
разностное уравнение третьего порядка
характеристический полином
область устойчивости
асимптотическое поведение решений
1. Баутин H.H. Поведение динамических систем вблизи границ области устойчивости. М.: Наука, 1984. – 176 с.
2. Васильев М.Д. Исследование одной математической модели трeхвидовой конкуренции // Математические заметки ЯГУ. – 2003. – Т. 10, № 2. – C. 33–39.
3. Джури Э. Импульсные системы автоматического регулирования. – М.: Физматгиз, 1963. – 456 с.
4. Кипнис М.М., Нигматулин Р.М. Устойчивость трехчленных линейных разностных уравнений с двумя запаздываниями // Автоматика и телемеханика. – 2004. – № 11. – С. 25–39.
5. Козак А.Д., Новоселов О.Н. Асимптотическое поведение решений линейного однородного разностного уравнения второго порядка // Математические заметки. – 1999. – Т. 66, Вып. 2. – С. 211–215.
6. Кудинов А.Ф. Общее решение разностного уравнения третьего порядка // Вестник ВГУ. Серия: Физика. Математика. – 2009. – № 2. – С. 69–70
7. Садовский П.А. Критические случаи устойчивости математической модели трехвидовой популяции // Вестник МГТУ им. Н.Э. Баумана. Серия «Естественные науки». – 2006. – № 4. – С. 61–71.
8. Elaydi S. An introduction to difference equations. New York: Springer, 2005. – 546 p.
9. Parhi N., Tripathy A.K. On the behavior of solutions of a class third order difference equations // Journal of Difference Equations and Applications. – 2002. – Vol. 8, № 5. – P. 415–426.

Важнейшим свойством дискретной системы, присущим всей системе, а не только отдельным её траекториям, является устойчивость. Известно [8], что исследование устойчивости решений нелинейного разностного уравнения сводится к выяснению расположения корней характеристического полинома соответствующего линеаризованного уравнения. При этом особую сложность представляет изучение критических случаев (называемых также граничной устойчивостью), когда некоторые корни характеристического полинома на комплексной плоскости попадают на единичную окружность [1, 3, 7].

Изучение критических случаев в теории устойчивости напрямую связано с исследованием границы области асимптотической устойчивости в пространстве параметров. Устойчивость в критических случаях изучена для непрерывных систем [1], но мало исследована для дискретных систем. Нашей целью является восполнение этого пробела для случая дискретной системы третьего порядка, описываемой характеристическим уравнением

nigmatul01.wmf (1)

где a, b, c ∈ R. Уравнение (1) называется устойчивым, если все его решения ограничены, и асимптотически устойчивым, если все его решения стремятся к нулю при n → ∞.

В работе [5] изучено более простое, чем (1), уравнение второго порядка nigmatul02.wmf Его область устойчивости в плоскости параметров a, b ∈ R такова: nigmatul03.wmf. В [5] указаны участки границы этой области, на которых возникают различные типы решений: циклы, предельные циклы, псевдошумовые решения. Уравнения и системы третьего порядка в непрерывном случае изучались в [2, 7], а в дискретном – в [9, 6]. В [9] для уравнения (1) получены достаточные условия колебательности – неколебательности решений в виде ограничений на коэффициенты уравнения. Для близкого к (1) уравнения nigmatul04.wmf с запаздываниями m, k ∈ N в [4] получено полное описание области асимптотической устойчивости в пространстве параметров a, b, k, m.

В настоящей работе мы полностью описываем асимптотическое поведение решений уравнения (1) при n → ∞, когда значения коэффициентов уравнения (1) находятся на границе его области устойчивости.

Граница области асимптотической устойчивости уравнения (1)

Известно, что нулевое решение уравнения (1) асимптотически устойчиво тогда и только тогда, когда все корни его характеристического полинома

nigmatul05.wmf (2)

по модулю меньше единицы. С помощью известного алгебраического критерия устойчивости Шура – Кона в [3] приведены необходимые и достаточные условия расположения всех корней характеристического полинома (2) внутри единичной окружности в виде системы ограничений на коэффициенты:

nigmatul06.wmf

Эти же ограничения приводятся в [8] в преобразованном виде:

nigmatul07.wmf

Анализируя эти системы, область асимптотической устойчивости уравнения (1) можно записать в виде, удобном для графического изображения области и ее границ:

nigmatul08.wmf (3)

Область асимптотической устойчивости уравнения (1) в пространстве коэффициентов a, b, c ∈ R изображена на рисунке. Она представляет собой тело, ограниченное гиперболическим параболоидом nigmatul09.wmf и двумя плоскостями

nigmatul10.wmf nigmatul11.wmf

Область асимптотической устойчивости имеет ось симметрии – ось Ob.

pic_14.tif

Область асимптотической устойчивости уравнения (1)

Границу области асимптотической устойчивости уравнения (1) образуют четыре вершины: A(1, –1, –1), B(3, 3, 1), C(–1, –1, 1), D(–3, 3, –1), пять ребер: AC (α ∩ β), AD и CD (α ∩ γ), AB и CB (β ∩ γ), две грани ACD и ABC и гиперболический параболоид γ.

Свойства характеристического полинома и асимптотическое поведение решений уравнения (1) на границе области асимптотической устойчивости

В этом пункте мы указываем свойства корней характеристического полинома (2) и асимптотическое поведение решений уравнения (1) на каждом участке границы области асимптотической устойчивости.

1. В точке A(1, –1, –1) имеем

nigmatul12.wmf

Все корни характеристического полинома P(λ) по модулю равны 1, причем λ = –1 – корень кратности 2. Тогда общее решение уравнения (1) имеет вид

nigmatul13.wmf

и в общем случае является неограниченным.

2. В точке C(–1, –1, 1) имеем

nigmatul14.wmf

Все корни характеристического полинома P(λ) по модулю равны 1, причем λ = 1 – корень кратности 2. Тогда общее решение уравнения (1) имеет вид

nigmatul15.wmf

и в общем случае является неограниченным.

3. В точке B(3, 3, 1) имеем

nigmatul16.wmf

т.е. λ = –1 – корень кратности 3. Общее решение имеет вид

nigmatul17.wmf

и в общем случае неограниченно.

4. В точке D(–3, 3, –1) имеем

nigmatul18.wmf

т.е. λ = 1 – корень кратности 3. Общее решение имеет вид

nigmatul19.wmf

и в общем случае неограниченно.

5. На ребре

nigmatul20.wmf

имеем

nigmatul21.wmf

Все корни характеристического полинома P(λ) действительные и простые, причем nigmatul22.wmf. Общее решение уравнения (1) имеет вид

nigmatul23.wmf

и является ограниченным, при этом в общем случае nigmatul24.wmf не существует.

6. На ребре

nigmatul25.wmf

имеем

nigmatul26.wmf

При –3 < a < 1 многочлен Q(λ) имеет два комплексно сопряженных корня nigmatul27.wmf, по модулю равных 1. Тогда общее решение уравнения (1) имеет вид

nigmatul29.wmf

и является ограниченным при любых начальных условиях, при этом в общем случае nigmatul30.wmf не существует.

7. На ребре

nigmatul31.wmf

имеем

nigmatul32.wmf

При –1 < a < 3 многочлен Q(λ) имеет два комплексно сопряженных корня nigmatul28.wmf, по модулю равных 1. Тогда общее решение уравнения (1) имеет вид

nigmatul33.wmf

и является ограниченным, при этом в общем случае nigmatul34.wmf не существует.

8. На ребре

nigmatul35.wmf

имеем

nigmatul36.wmf

Все корни характеристического полинома P(λ) действительные, при этом λ = –1 – корень кратности 2, а при 1 < a < 3 модуль третьего корня nigmatul37.wmf. Тогда общее решение уравнения (1) имеет вид

nigmatul38.wmf

и в общем случае не является ограниченным.

9. На ребре

nigmatul39.wmf

имеем

nigmatul40.wmf

Все корни характеристического полинома P(λ) действительные, при этом λ = 1 – корень кратности 2, а при –3 < a < –1 модуль третьего корня nigmatul41.wmf. Тогда общее решение уравнения (1) имеет вид

nigmatul42.wmf

и в общем случае не является ограниченным.

10. Во внутренних точках треугольника

nigmatul43.wmf

имеем

nigmatul44.wmf

При c – 2 < a < –c имеем nigmatul45.wmf. Тогда очевидно, что многочлен Q(λ) при (a + 1)2 + 4c > 0 имеет пару действительных корней

nigmatul46.wmf,

по модулю меньших 1, при (a + 1)2 + 4c < 0 имеет пару комплексно сопряженных корней

nigmatul47.wmf

nigmatul48.wmf

nigmatul49.wmf

а при (a + 1)2 + 4c = 0 – действительный корень nigmatul50.wmf кратности 2, по модулю меньший 1. Получаем, что общее решение уравнения (1) имеет вид

nigmatul51.wmf

и является ограниченным, при этом в общем случае nigmatul52.wmf.

11. Во внутренних точках треугольника

nigmatul53.wmf

имеем

nigmatul54.wmf

При –c < a < c + 2 имеем nigmatul55.wmf. Тогда очевидно, что многочлен Q(λ) при (a – 1)2 – 4c > 0 имеет пару действительных корней nigmatul56.wmf, по модулю меньших 1, при (a – 1)2 – 4c < 0 имеет пару комплексно сопряженных корней

nigmatul57.wmf nigmatul58.wmf nigmatul59.wmf

а при (a – 1)2 – 4c = 0 действительный корень nigmatul60.wmf кратности 2, по модулю меньший 1. Получаем, что общее решение уравнения (1) имеет вид

nigmatul61.wmf

и является ограниченным, при этом в общем случае nigmatul62.wmf.

12. В точках, лежащих на гиперболическом параболоиде b = 1 – c2 + ac, в области nigmatul63.wmf, nigmatul64.wmf имеем

nigmatul65.wmf

Характеристический полином P(λ) имеет один действительный корень λ = –c, nigmatul66.wmf и пару комплексно сопряженных корней

nigmatul67.wmf nigmatul68.wmf nigmatul69.wmf

Тогда общее решение уравнения (1) имеет вид

nigmatul70.wmf

и является ограниченным при любых начальных условиях, при этом в общем случае nigmatul71.wmf не существует.

Выводы

На каждом из выделенных участков границы области устойчивости решения уравнения (1) обладают следующими особенностями асимптотического поведения.

В каждой из четырех вершин A, B, C, D в общем случае решения не ограничены и имеют полиномиальный рост, поэтому система (1) неустойчива. На рёбрах AC, AD, BC все решения ограничены, причем решения могут быть чисто периодическими, поэтому система (1) устойчива (не асимптотически). На ребрах AD и BC возможны так называемые псевдошумовые решения. На рёбрах AB и CD в общем случае решения не ограничены и имеют рост линейный по n, поэтому система (1) неустойчива.

Во внутренних точках граней ACB и ACD все решения ограничены, система устойчива (не асимптотически). В общем случае в указанной области решения имеют вид затухающих колебаний. В точках, лежащих на гиперболическом параболоиде, все решения ограничены, поэтому система (1) устойчива (не асимптотически). В общем случае в указанной области решения являются колебательными. Здесь возникают устойчивые циклы и псевдошумовые решения.

Заключение

Мы сделали полный анализ асимптотических свойств систем третьего порядка, когда их параметры находятся на границе области устойчивости. При специальном выборе начальных условий (посредством обнуления констант при неограниченных слагаемых в формуле общего решения) в исследованных областях, где диагностирована неустойчивость общего решения, можно выделить ограниченные или даже сходящиеся к нулю решения. Например, в п. 9 раздела «Свойства характеристического полинома и асимптотическое поведение решений уравнения (1) на границе области асимптотической устойчивости» настоящей статьи, если (a + 2)x0 – (a + 1)x1 – x2 = 0, то в общем решении xn = C1 + C2n + C3(–2–a)n имеем C2 = 0 и nigmatul72.wmf. Поиск областей в пространстве начальных значений, которые дают ограниченные решения (проблема частичной устойчивости), требует отдельного тщательного исследования.

Работа поддержана грантом № 2807 Министерства образования России.

Рецензенты:

Дильман В.Л., д.ф.-м.н., профессор, заведующий кафедрой прикладной математики, ФГБОУ ВПО «Южно-Уральский государственный университет» (национальный исследовательский университет)», г. Челябинск;

Карачик В.В., д.ф.-м.н., профессор кафедры математического и функционального анализа, ФГБОУ ВПО «Южно-Уральский государственный университет» (национальный исследовательский университет)», г. Челябинск.


Библиографическая ссылка

Нигматулин Р.М., Кипнис М.М. СВОЙСТВА ДИСКРЕТНЫХ СИСТЕМ ТРЕТЬЕГО ПОРЯДКА НА ГРАНИЦЕ ИХ ОБЛАСТЕЙ УСТОЙЧИВОСТИ // Фундаментальные исследования. – 2015. – № 9-1. – С. 39-43;
URL: http://www.fundamental-research.ru/ru/article/view?id=38962 (дата обращения: 29.09.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074