Локализация инвариантных компактов системы Lorenz-84 Текст научной статьи по специальности «Математика»

Математика к Математическое

моделирование

Ссылка на статью:

// Математика и Математическое моделирование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2015. № 4. С. 54-65.

Б01: 10.7463/шаШш.0415.0812317

Представлена в редакцию: 24.09.2015

© МГТУ им. Н.Э. Баумана

УДК 517.938

Локализация инвариантных компактов системы Lorenz-84

Рамазанова X. М.1'* * xadi91@yandex.ru

1МГТУ им. Н.Э. Баумана, Москва, Россия

Один из методов построения локализирующих множеств основан на использовании функций, определенных на фазовом пространстве системы - так называемый функциональный метод локализации. В статье с помощью данного метода получены оценки положения инвариантных компактных множеств автономной системы Lorenz-84, используемой при построении некоторых метеорологических моделей. Рассмотрен простейший вариант системы с отсутствующими термическими перегрузками, и общий вариант системы, в которой при некоторых значениях параметров возникает хаотическая динамика. При всех значениях параметров системы описано компактное локализирующее множество.

Ключевые слова: локализация; локализирующее множество; инвариантный компакт; система Ьоге^-84.

Введение

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

Подмножество М фазового пространства динамической системы называют инвариантным, если для любой точки х0 е М любое решение этой системы, проходящее через точку х0, целиком содержится в М (решение рассматривается на максимальном интервале существования). Примерами инвариантных компактов являются: положения равновесия, периодические траектории, аттракторы и репеллеры, инвариантные торы. Инвариантные компактные множества во многом определяют структуру фазового портрета динамической системы.

Один из методов построения локализирующих множеств основан на использовании функций, определенных на фазовом пространстве системы - так называемый функциональный метод локализации [1]. Основные положения метода были опубликованы в 1995 г. при решении задач локализации предельных циклов [2,3]. В окончательном виде суть функционального метода локализации для автономных систем изложена в [4]. Тема локализации инвариантных компактов в литературе обсуждается достаточно активно. Исследования направлены как на развитие метода и его использование для динамических систем других классов, так и на исследование конкретных динамических систем. В исследовании конкретных систем можно выделить два направления. Первое направление исследований по локализации - автономные системы со сложным поведением: системы Лоренца [5,6], ПРТ-системы [7] , Ланфорда [8]. Второе направление исследований - прикладные системы, используемые в биологии, физике [9,10,11]. Система Ьогеи2-84, предложенная Лоренцем в 1984 году, представляет собой простейшую модель общей циркуляции атмосферы в средних широтах [12]. Модель использовалась в различных климатологических исследованиях, а ее поведение активно изучалось. Численные и аналитические исследования системы можно найти в [13]. Бифуркационный анализ положений равновесия системы [14] показал наличие хаоса при определенных значениях параметров модели. Проблема локализации инвариантных компактных множеств системы рассмотрена в [15].

В данной работе рассматривается задача локализации инвариантных компактов системы Ьогеи2-84. Работа организована следующим образом. В разд.1 описан функциональный метод локализации, перечислены простейшие свойства локализирующих множеств. В разд. 2 анализируется простейший вариант системы Ьогеи2-84, когда отсутствуют термические нагрузки. В разд. 3 проводится исследование для общего варианта автономной системы Ьогеи2-84, в результате которого построено компактное локализирующее множество, описанное неравенством. В заключении подводятся краткие итоги работы.

1. Функциональный метод локализации

Рассмотрим нелинейную систему дифференциальных уравнений

X = /(х) , (1)

где

X = (Х1.....хп)теКп, Г(х) = (Мх).....^(х))те Ст(Ип).

Пусть f = 2 г=1 (х)^ — векторное поле на Мп, соответствующее системе (1). Для

произвольной функции уе С(Мп) через ¿/^(х) обозначим производную Ли этой функции по векторному полю /:

¿ж*) = 2"=1 п(*) ^ = </(*)' <р(х)).

Множество

= {х е : = 0}.

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

Положим для функции <р е С(Еп)

Vsup = sup (р(х), yinf = inf (р(х).

Теорема 1 [1,4]. Для любой функции ф(х) е С(Еп) все компактные инвариантные множества системы (1) содержатся в множестве

Пф = (х е Еп: (pinf < (р(х) < (psuv).

Рассматривая различные локализирующие функции, мы получаем некоторый набор локализирующих множеств. Пересечение всех найденных локализирующих множеств также является локализирующим множеством, причем пересечение, как оценка положения инвариантных компактов, представляет собой более сильный результат, чем каждое локализирующее множество в отдельности. Один из подходов к решению задач локализации - выбрать параметрическое семейство локализирующих функций и затем найти пересечение полученного параметрического семейства локализирующих множеств.

Отметим некоторые свойства локализирующих множеств.

1. Пусть ^(х) е С(Еп) и 'ф(х) = h(^(x)), где h: Е — Е - гладкая строго монотонная функция с ненулевой производной. Тогда множества и П^, соответствующие функциям <р и совпадают. В частности, это верно при h(t) = ct + d, с Ф 0.

2. Если функция ^(х) е С®(ЕП) достигает в точке х* своего наименьшего (наибольшего) значения в Еп, то ^inf = ф(х*) (соответственно ^>inf = ф(х*)), а локализирующее множество Пф можно задать неравенством ф(х) < ^sup (соответственно ^nf < ^(х)).

3. Пусть даны функции <pi е СЮ(ЕП), i = 0,1,2, ... Построим множества П0 = Hi = Oi-1 П Oi-i i, i = 1,2, ... , где

fii-U = (*е Еп: VifHi-i) < yt (х) < ^pCty-i)} ; ^¿,sup(^i-i) = sup щ (x); Viffli-i) = inf м (x);

Sy. (ni-1) 4 (ni-i)

(^i-i) = (x е Lf(p(x) = 0}. Тогда П0 2 П1 2 ...2 Hi 2... и пересечение этих множеств содержит все компактные инвариантные множества системы (1). Таким образом, возникает некоторый итерационный процесс построения вложенных локализирующих множеств.

2. Система Lorenz-84 при отсутствии термических нагрузок

Система Lorenz-84 имеет вид [12]:

х = —ах — у2 — z2 + aF(t), у = —у + ху — fixz + G(t), z = —z + fixy + xz.

Она была предложена как модель общей циркуляции атмосферы в средних широтах.

В этой модели х обозначает интенсивность западной циркуляции, у и z - синусные и косинусные составляющие большой бегущей волны. Функции F(t) и G(t) учитывают термические нагрузки, вызванные температурным контрастом между зимой и летом, а также между океаном и земной поверхностью, соответственно. Коэффициенты а и @ — положительные параметры модели.

В своей работе Лоренц отметил, что F и G являются периодическими функциями, периодом которых является год. Однако, в своих расчетах Лоренц считал значение G фиксированным, а значение F — различным для зимних и летних условий. В данной работе ограничимся автономным случаем системы, т. е. полагаем, что термические нагрузки постоянны, т.е. F(t) = /, G(t) = д. В этом случае получаем систему:

гх = —ах — у2 — z2 + а/,

j у = —у + xy — pxz + д, (2)

^z = —z + flxy + xz.

Для начала рассмотрим простой вариант системы, когда параметры f и д равны нулю. Тогда система (2) примет следующий вид:

гх = —ах — у2 — z2, у = —у + ху — fíxz \z = —z + fíxy + xz.

(3)

Система (3) при всех значениях параметров имеет единственное положение равновесия - точку (0,0,0). Оказывается, что нулевое положение равновесия глобально устойчиво.

Для анализа системы (3) используем цилиндрические координаты.

Сделаем следующую замену переменных в исходной системе:

у = р cos z = р sin х = z.

Тогда система примет вид:

= —az — р2, \р = P(z — 1), [ф = pz.

Видно, что в системе выделяется независимая подсистема из первых двух уравнений:

i = —az — p2, (4)

Р = P(z — 1). ( )

Поэтому анализ системы начнем с этой подсистемы.

Система (4) имеет простое поведение. У нее единственное положение равновесия р = 0, z = 0, являющееся глобально асимптотически устойчивым. Действительно, рассмотрим V(p, z) = р2 + z2. Эта функция положительно определена, непрерывно дифференцируема, а ее производная в силу системы V(p, z) = —2р2 — 2az2 < 0 (при а > 0). Следовательно, функция V(p, z) является функцией Ляпунова, а положение равновесия, согласно теореме Ляпунова об асимптотической устойчивости, глобально асимптотически устойчиво.

Приведенные рассуждения позволяют сделать вывод, что в системе существует единственный инвариантный компакт - нулевое положение равновесия. К аналогичным выводам приводит и применение функционального метода локализации. Действительно, рассмотрим многочлен c неопределенными коэффициентами вида:

ф(х, у, z) = Ах2 + By2 + Cz2. Найдем производную этой функции в силу системы (3):

Lf(p = -2Аах2 - 2Ву2 - 2Cz2 - 2Аху2 - 2Axz2 + 2Вху2 + 2Cxz2 - 2Bfixyz + 2Cfixyz.

Выберем случай, когда степень LfV равняется 2, для чего необходимо и достаточно, чтобы А = В = С. Для определенности примем А = В = С = 1. Тогда <р0(х, у, z) = х2 + +у2 + z2, а производная этой функции в силу системы (3) равна

LfVo(x) = -2ах2 - 2у2 - 2z2. Уравнение Lf^o = 0 принимает вид

-2ах2 - 2у2 - 2z2 = 0. Отсюда вытекает, что универсальное сечение S^o состоит из единственной точки х = у = = z = 0, а локализирующее множество имеет вид х2 + у2 + z2 < 0, т.е. состоит из единственной точки - положения равновесия системы. Таким образом, единственным инвариантным компактом системы является ее положение равновесия (0,0,0).

Отметим, что локализирующая функция ф0(х,у,z) оказалась функцией Ляпунова для системы (3).

3. Система Lorenz-84 при наличии термических нагрузок

Теперь рассмотрим общий случай автономной системы (2) (когда /, д Ф 0). Наличие слагаемых f и д заметно усложняет ее поведение.

Выясним положения равновесия системы (2). Для этого рассмотрим систему уравнений

-ах - у2 - z2 + af = 0,

-у + ху - fixz + д = 0, (5)

-z + fixy + xz = 0.

Второе и третье уравнения системы образуют линейную подсистему относительно у и z. Решив ее и подставив решения в первое уравнение, получим уравнение третьей степени:

(1 + (З2)х3 - (ffi2 + f + 2)х2 + (1 + 2f)x + ^-f = 0. (6)

Таким образом, система может иметь одно, два или три положения равновесия. Проведем анализ количества корней уравнения (6) . Введем обозначения

= f(32 + f + 2 r =1 + 2/ r = g2-af 1+ fi2 ; h2 1 + fi2; h3 a(1 + fi2). Тогда уравнение (6) примет вид

x3 + Вхх2 + B2x + B3 = 0 (7)

Заменим в уравнении (7) переменную х на новую переменную у = х + -1 и, упростив полученное выражение, получим приведенное уравнение

у3 + ру + ч = 0,

Вл . г-» -в1 Вл В? . г-,

где р = + В-, Ч = -ф--1^ + Вз.

Число действительных решений кубического уравнения зависит от знака дискрими-

Р 3 $ —

нанта В = (-)3 + (~)-, а именно, при В >0 уравнение имеет одно решение (система имеет одно положение равновесия), а при В <0 уравнение имеет три решения (система имеет три положения равновесия).

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

Теорема 2. Для того, чтобы многочлен второй степени общего вида

ф(х, у, х) = Ах- + Ву- + Сх- + 2Вху + 2Ехг + + 2Кх + Ну + 2Мг

имел производную в силу системы (2) степени не выше двух, необходимо и достаточно, чтобы

В = Е = ^ = 0, А = В = С.

Доказательство. Непосредственным вычислением получаем

Ьф(х) = 2(-А + В + Ер)ху- + 2(-Л + С - рр)хг- + 2(-Вр + Ср + 2¥)хух -

2Ву3 - 2Вуг- - 2Ех3 - 2Еу-х + 2(В + РЕ)х-у + 2(-рВ + Е)х-г + Г, где Г — сумма слагаемых степени не выше двух. Равенство нулю коэффициентов при слагаемых третьей степени приводит к системе уравнений

г -А + В = 0, -А + С-Е = 0, -рв + РС + 2Е = 0, В = 0, Е = 0, В + РЕ = 0, , -РВ + Е = 0.

Решением этой системы является В = Е = F = 0, А = В = С. Теорема доказана. Итак, производная квадратичной функции в силу системы является квадратичной, если эта функция, с точностью до числового множителя, имеет вид

(р(х, у, г) = х- + у2 + г- + 2Кх + 2Ьу + 2Мг. Эту функцию рассмотрим в качестве локализирующей. Ее производная в силу системы:

ф(х,у, х) = 2(-ах- - (К + 1)у- - (К + 1)г- + (Ь + Мр)ху + (М - Ьр)хг + +(/ - К)ах + (д - Ь)у - Мг + (Ка/ + Ьд)).

В результате приходим к задаче

х2 + у2 + г2 + 2Кх + 2Ьу + 2Мг ^ extг, ах2 + (К + 1)у2 + (К + 1>2 - (Ь + М^)ху - (М- Ьр)хг - (8)

-(/ - К)ах - (д - Ь)у + Мг - (Ка[ + 1д) = 0.

Поскольку целевая функция имеет глобальный минимум, нижнее решение задачи <Pinf совпадает с этим минимумом, так что неравенство ф(х) > в данном случае оказывается тривиальным (выполняется во всем фазовом пространстве). В рассматриваемой задаче ограничение описывает поверхность 2-го порядка. Положительно определенная квадратичная функция имеет на этой поверхности конечную точную верхнюю грань, если эта поверхность - эллипсоид. А это равносильно знакоопределенности матрицы квадратичной формы в левой части ограничения. Запишем матрицу квадратичной формы:

/ I + мр М-

а 2 2

Ш =

а

I + МБ --^ К + 1 0

М-ЬВ

0 К + 1 I

2

Согласно критерию Сильвестра, для знакоопределенности квадратичной формы необходимо и достаточно, чтобы все угловые миноры четного порядка были положительны, а миноры нечетного порядка имели одинаковый знак, т. е.

Аз

Д2>0, -3 >0.

2 А

В нашем случае (считая с правого нижнего угла)

Д1= (К + 1); А2= (К + 1)2;

( (В2 + 1)(Ь2 + М2)\ Аз= (К + 1) 1а(К + 1) - ^--Л

Таким образом, получаем следующее условие

(В2 + 1)(1? + М2) К + 1 > —-^--.

4а

Найти точное аналитическое решение задачи (8) вряд ли возможно. Однако можно найти оценку ее решения сверху. Это позволит записать соответствующее локализирующее множество.

Рассмотрим более простой частный случай Ь = М = 0. В этом случае экстремальная задача выглядит следующим образом:

х2 + у2 + х2 + 2Кх ^ ех1г, ах2 + (К + 1)у2 + (К + 1)г2 - (/ - К)ах - ду - Ка[ = 0.

Исключим из целевой функции у2 + г2 с помощью ограничения:

Ка[

( а \ , /а(/ - К) \ д

-ТТгУ + Н"ТГ- + 2Т + кТ!у + к + 1

ех№.

Тогда переменная 2 будет входить только в ограничение. Поэтому ее можно исключить:

((1 -^у+т1

I ах2 + (К + 1)у2 - a(f - i

д Kaf

х + ^—7 У + ^—7 ^ extr, К + 1* К + 1

К)х - ду- Ка[ < 0.

Так как в последней задаче целевая функция линейна по у, то решение задачи достигается на границе области, т.е. неравенство в ограничении можно заменить равенством:

((1 +ш-+2«;

I ах2 + (К + 1)у2 - a(f - i

д Kaf

х + ^—7 У + ^—7 ^ extr, К + 1* К + 1

К)х - ду- Ка[ = 0.

Полученная задача представляет собой поиск максимума квадратичной функции на эллипсе. Ее можно преобразовать в задачу одномерной оптимизации, если задать эллипс параметрически. Для этого преобразуем ограничение, выделив квадраты:

______ а

где

а(*-^)2 + (* + 1)(у- 2^+1))2 = *2,

2 a(f - К)2 д2 a(f + К)2 д2

R2 = 4 J + J + Kaf = 4 J + ■ w

4 4(K + 1) 4 4(K + 1)

Далее введем замену переменных

R f-K R д

х = —cos6 +---—, у = -sin6 +■

^ 2 ' ' 4кП 2(К + 1)"

Тогда получим задачу одномерной оптимизации:

/ а ^Я- Я а Я

(1 ——-) — со8-е + (/ + К) — соз9 + * л ,_=8т9 + Р ^ extг,

V К + 1/ а К + 14кГГ

где

Р = (л+ а \ V - к)- + К(Г д- + КаГ

р = (1 + + к(г -к) + 2(К +1)- + ТП

Здесь параметры удовлетворяют условиям а >0, @ >0, К + 1 > 0. Последнее условие усилим, потребовав К + 1 > а.

Используя неравенства cos26 < 1, pcosQ + qsin6 < ^р2 + q2, получим

/ a ЛИ2 l(f + К)2 ~д2

Vsup <(1 -*т1)т + T^ + CfW +

Подставляя параметр P в выражение, находим

(f-K)2, д2 , д2

где

+jOwTT)+4(FT1)2+K(2f - v+ад

2 Л а (К + 1)3)\у' ' (К + 1)/

В результате получаем семейство локализирующих множеств, которые описываются неравенством

2 2 2 (/ + К)2 д2(а + К + 1) _ ,

(X + К)2 + у2 + 22 < (У 2 ) + +1)2 ) + Н!(К)' К >(а- 1).

Пересечение этого семейства можно описать неравенством

+ *2 < К > (а - 1) + ® + ^ - С + *>2)- (9)

Графическое изображение локализирующего множества (9) в случае а = 0,5, ¡3 = 3,/ = 1,5, д = 0,5 представлено на рис.1.

Рис.1. Ограниченная траектория системы Ьогет-84 и локализирующее множество

Заключение

В данной работе получены оценки положения инвариантных компактных множеств автономной системы Ьогеи2-84, используемой при построении некоторых метеорологических моделей (эта система рассматривалась и как самостоятельная модель и как составная часть более сложной модели). Рассмотрен простейший вариант системы с отсутствующими термическими перегрузками, и общий вариант системы, в которой при некоторых значениях параметров возникает хаотическая динамика. При всех значениях параметров системы построено компактное локализирующее множество. Результат частично перекрывается с результатом Старкова К. Е. [15], но дает дополнительную информацию.

Работа выполнена при финансовой поддержке РФФИ (грант № 15-07-005489).

Список литературы

1. Канатников А.Н., Крищенко А.П. Инвариантные компакты динамических систем. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 231 с.

2. Крищенко А.П. Локализация предельных циклов // Дифференциальные уравнения. 1995. Т. 31, № 11. С. 1858-1865.

3. Крищенко А.П. Области существования циклов // Докл. РАН. 1997. Т. 353, № 1. С. 17-19.

4. Крищенко А.П. Локализация инвариантных компактов динамических систем // Дифференциальные уравнения. 2005. Т. 41, № 12. С. 1597-1604.

5. Krishchenko A. P., Starkov K.E. Localization of compact invariant sets of the Lorenz system // Phys. Lett. A. 2006. Vol. 353, no. 5. P. 383-388. DOI: 10.1016/j.physleta.2005.12.104

6. Li D., Lu J., Wu X., Chen G. Estimating the bounds for the Lorenz family of chaotic systems // Chaos, Solitons and Fractals. 2005. Vol. 23, no. 2. P. 529-534. DOI: 10.1016/j.chaos.2004.05.021

7. Канатников А.Н. Локализация инвариантных компактов ПРТ-системы // Вестник МГТУ. Сер. Естественные науки. 2007. № 1. С. 3-18.

8. Krishchenko A. P., Starkov K.E. Localization of compact invariant sets of nonlinear systems with application to the Lanford systems // Int. J. of Bifurcation and Chaos. 2006. Vol. 16, iss. 11. P. 3249-3256. DOI: 10.1142/S0218127406016768

9. Канатников А.Н., Федорова Ю.П. Локализация инвариантных компактов двумерных непрерывных динамических систем // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 7. С. 159-174. DOI: 10.7463/0713.0583104

10. Starkov K. E. Bounding a domain which contains all compact invariant sets of the Bloch system // Int. J. of Bifurcation and Chaos. 2009. Vol. 19, no. 3. P. 1037-1042. DOI: 10.1142/S0218127409023457

11. Starkov K. E. Bounds for the domain containing all compact invariant sets of system modeling dynamics of acoustic gravity waves // Int. J. of Bifurcation and Chaos. 2009. Vol. 19, no. 10. P. 3425-3432. DOI: 10.1142/S0218127409024864

12. Lorenz E.N. Irregularity: a fundamental property of the atmosphere // Tellus A. 1984. Vol. 36, no. 2. P. 98-110. DOI: 10.1111/j.1600-0870.1984.tb00230.x

13. Masoller C., Schifino A., Romanelli L. Characterization of strange attractors of Lorenz model of general circulation of the atmosphere // Chaos, Solitions and Fractals. 1995. No. 6. P. 357-366. DOI: 10.1016/0960-0779(95)80041 -E

14. Shilnikov A., Nicolis G., Nicolis C. Bifurcation and predictability analysis of a low-order atmospheric circulation model // Journal of Bifurcation and chaos. 1995. Vol. 5, no. 6. P. 1701-1711. DOI: 10.1142/S0218127495001253

15. Starkov K.E. Localization of compact invariant sets of the Lorenz'1984 model // Springer Proceedings in Physics. 2009. Vol. 132. P. 915. DOI: 10.1007/978-3-642-03085-7 225

Mathematics and Mathematical Madelling of the Bauman MSTU, 2015, no. 4, pp. 54-65.

DOI: 10.7463/mathm.0415.0812317

Received: 24.09.2015

Mathematics & Mathematical Modelling

Electronic journal of the Bauman MSTU

© Bauman Moscow State Technical Unversity

Localization of Compact Invariant Sets of the Lorenz'1984 System

Ramazanova Kh. M.

1,*

xadi91@yandex.ru

1Bauman Moscow State Technical University, Moscow, Russia

Keywords: localization, localizing set, invariant compact set, Lorenz'1984 model

Localization of compact invariant sets of a dynamical system is one way to conduct a qualitative analysis of dynamical system. The localization task is aimed at evaluating the location of invariant compact sets of systems, which are equilibrium, periodic trajectories, attractors and repellers, and invariant tori. Such sets and their properties largely determine the structure of the phase portrait of the system. For this purpose, one can use a localization set, i.e. a set in the phase space of the system that contains all invariant compact sets.

This article considers the problem of localization of invariant compact sets of an Autonomous version of the Lorenz-84 system. The system represents a simple model of the General circulation of the atmosphere in middle latitudes. The model was used in various climatological studies. To build localization set of the system the so-called functional localization method is applied. The article describes the main provisions of this method, lists the main properties of the localization sets. The simplest version of the Lorenz-84 system when there are no thermal loads is analyzed, and a common variant of the Autonomous Lorenz-84 system, in which for some values of system parameters chaotic dynamics occurs is investigated. In the first case it is shown that the only invariant compact set of the system is its equilibrium position, and localization function turned out to be a Lyapunov function of the system. For the General version of the system a family of localization sets is built and the intersection of this family is described. Graphical illustration for the localization set at fixed values of the parameters is shown. The result of the study partially overlaps with the result of K.E. Starkov on the subject, but provides additional information.

The theme of localization of invariant compact sets is discussed quite actively in the literature. Research focuses both on the development of the method and its application to dynamical systems of other classes, and on the investigation of specific dynamical systems.

References

1. Kanatnikov A.N., Krishchenko A.P. Invariantnye kompakty dinamicheskih sistem [Compact invariant sets of dynamical systems]. Moscow, Bauman MSTU publ., 2011, 232 p. (in Russian)

2. Krishchenko A.P. Localization of limit cycles. Differentsial'nye uravneniya, 1995, vol. 31, no. 11, pp. 1858-1865 (English translation: Differential Equations, 1995, vol. 31, no. 11, pp. 1826-1833).

3. Krishchenko A.P. Domains of Existence of Cycles. Doklady Akademii nauk. 1997, vol. 353, no. 1, p. 17-19. (English translation: Doklady Mathematics, 1997, vol. 55, pp. 176-178).

4. Krishchenko A.P. Localization of invariant compact sets of dynamical systems. Differentsial'nye Uravneniya, 2005, vol. 41, no. 12, pp. 1597-1604 (English version of journal: Differential Equations, 2005, vol. 41, no. 12, pp. 1669-1676).

5. Krishchenko A.P., Starkov K.E. Localization of compact invariant sets of the Lorenz system. Phys. Lett. A, 2006, vol. 353, no. 5. pp. 383-388. DOI: 10.1016/j.physleta.2005.12.104

6. Li D., Lu J., Wu X., Chen G. Estimating the bounds for the Lorenz family of chaotic systems, Chaos, Solitons and Fractals, 2005, vol. 23, no. 2. pp. 529-534. DOI: 10.1016/j.chaos.2004.05.021

7. Kanatnikov A.N. Localization of Invariant Compact Sets of PRT-System. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki = Herald of the Bauman MSTU. Ser. Natural Sciences, 2007. no. 1, pp. 3-18. (in Russian)

8. Krishchenko A. P., Starkov K.E. Localization of compact invariant sets of nonlinear systems with application to the Lanford systems. Int. J. of Bifurcation and Chaos, 2006, vol. 16, iss. 11, pp. 3249-3256. DOI: 10.1142/S0218127406016768

9. Kanatnikov A.N., Fedorova Yu.P. Localization of invariant compact sets of two-dimensional continuous dynamical systems. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2013, no. 7, pp. 159-174. DOI: 10.7463/0713.0583104 (in Russian).

10. Starkov K. E. Bounding a domain which contains all compact invariant sets of the Bloch systemro Int. J. of Bifurcation and Chaos, 2009, vol. 19, no. 3, pp. 1037-1042. DOI: 10.1142/S0218127409023457

11. Starkov K. E. Bounds for the domain containing all compact invariant sets of system modeling dynamics of acoustic gravity waves. Int. J. of Bifurcation and Chaos, 2009, vol. 19, no. 10, pp. 3425-3432. DOI: 10.1142/S0218127409024864

12. Lorenz E.N. Irregularity: a fundamental property of the atmosphere. Tellus A, 1984, vol. 36, no. 2, pp. 98-110. DOI: 10.1111/j.1600-0870.1984.tb00230.x

13. Masoller C., Schifino A., Romanelli L. Characterization of strange attractors of Lorenz model of general circulation of the atmosphere. Chaos, Solitions and Fractals, 1995, no. 6, p. 357-366. DOI: 10.1016/0960-0779(95)80041 -E

14. Shilnikov A., Nicolis G., Nicolis C. Bifurcation and predictability analysis of a low-order atmospheric circulation model. Journal of Bifurcation and chaos, 1995, vol. 5, no. 6, pp. 1701-1711. DOI: 10.1142/S0218127495001253

15. Starkov K.E. Localization of compact invariant sets of the Lorenz'1984 model. Springer Proceedings in Physics, 2009, vol. 132, p. 915. DOI: 10.1007/978-3-642-03085-7 225