CC BY
9УДК 517.925/926
Сингулярно-гиперболический аттрактор отображения многомерного цилиндра
В. Н. Белых, Д. А. Гречко
Волжский государственный университет водного транспорта,
Национальный исследовательский Нижегородский государственный университет им. Н. И. Лобачевского
Нижний Новгород, 603950 E-mail: [email protected], [email protected]
Аннотация. В работе рассматривается многомерное отображение с одной кусочно-гладкой периодической нелинейностью. Получены условия, при которых отображение имеет аттракторы, лежащие в полнотории. Даны критерии, определяющие колебательный и вращательный тип аттракторов. Доказана теорема, указывающая область параметров, при которых аттракторы как колебательные, так и вращательные, являются сингулярно-гиперболическими. При этих условиях динамическое поведение траекторий отображения представляет собой пример динамического хаоса.
Ключевые слова: динамика систем, нелинейное отображение, аттракторы, гиперболичность, бифуркации
Singular hyperbolic attractor of a multidimensional cylinder map
V. N. Belykh, D. A. Grechko
Volga state university of water transport, Lobachevsky university, Nizhny Novgorod.
Abstract. We consider a multidimensional map with a piecewise-smooth periodic nonlinear function of form F: (x,yi) ^ (x + S — ag(x) + ^"=1 Vi, Xi(—big(x) + yi)), i = 1,. ..,n, where x £ S1,y £ Rn, g(x) = g(x + 2n), a, S, Xi, bi are parameters, \g'(x)\ > h. We divide nonwandering trajectories into oscillating and rotating type via the rotation number r = Xnk'. We prove the existence of
absorbing domain D containing the attractors of F. Namely, the next statement is true.
Theorem 1. Let 0 < Xi < X* < 1, i = 1~n. Then the solid torus D = {||y|| < A, x £ S1} is the absorbing domain D, FD с D.
In order to specify the attractors of the map F we introduce two auxiliary 1D maps of the circle F± : x ^ x + S — ag(x) ± y, and prove the next
Theorem 2. 1. If both maps F + and F- have only oscillating attractors then the map F has only oscillating attractor as well. 2. If both maps F + and F- have only rotating attractor then the map F also has only rotating attractor.
We studied the hyperbolic properties of the map F attractors which imply that in each point of an attractor there exist unstable (stable) cone invariant under F (F-1 respectively) with local extension (contraction) property in it. For this purpose we consider the variation equation along attractor trajectories (£, щ) ^ ((1 — ag')£ + Y1 n=1 ni, Xi(—bg'£ + ni)),i = 1,n. Introducing the unstable cone Ku = {£, ni \ ni = ai£, \аА < X,i = 1, n} we obtain a condition when the map £ ^ £ inside the cone is extension. Then finding the condition for value x = X* defining the cone Ku we prove the main theorem.
© В. Н. БЕЛЫХ, Д. А. ГРЕЧКО
Theorem 3. Each attractor (oscillating or rotating) of the map F is singularly hyperbolic in a parameter region ah > 2 + x* + £ > 0.
Finely we note that for this explicitly given region of parameters indipending on whether the attractor is oscillatory or rotatory the attractor is hyperbolic. The dynamical behavior of the attractor trajectories exhibits the example of rear dynamical chaos.
Keywords: dynamical systems, nonlinear maps, attractors, hyperbolicity, bifurcations. MSC 2010: 34D45
Введение
Задача об условиях существования аттракторов и бифуркаций, формулируемых в терминах правых частей динамических систем, возникла при анализе состояний равновесия предельных циклов Пуанкаре и гомоклинических петель сепаратрис в конкретных системах. В период последних почти 50 лет эта задача приобрела широкий интерес в связи с возникновением концепции странного аттрактора как образа динамического хаоса.
Странные аттракторы как притягивающие инвариантные множества целых неустойчивых траекторий можно разделить на 3 типа: гиперболические, структура которых не меняется во всех точках интервала параметра, характеризующего деформации динамической системы; сингулярно-гиперболические, структура которых меняется только в точках бифуркаций; квазистранные — странные не на интервале, а на точечном, как правило, канторовом множестве параметров.
В классе динамических систем с дискретным временем (отображений) примерами гиперболических аттракторов служат диффеоморфизмы Аносова [1], примеры сингулярно-гиперболических — аттрактор Лоренца для модельных отображений [2, 3], аттракторы Лози [4] и Белых [5], примерами квазиаттракторов — аттрактор Эно [6] и другие. Условия существования приведенных аттракторов выписываются аналитически в силу конкретно заданного вида отображений.
В настоящей работе рассматривается многомерное отображение с одной непрерывной, кусочно-гладкой, периодической нелинейностью F : S1 х Rn ^ S1 х Rn, записанное в нормальной форме вида:
{
х = х + 5 - ag{x) + J2 Vi; . — ,n ^
Vi = Xi{-big{x)+ Vi),
где a,ó,bi,\i,i = 1,n — положительные параметры, x Е S 1,y = colomn(y1,...,yn) Е Rn переменные, g(x) — 2^-периодическая, непрерывная, кусочно-гладкая функция с разрывами производной в точках экстремума, удовлетворяющая условию
max \g(x)\ = 1;
п 2п
/ g(x)dx = 0; (0.2)
J о
\g'(x)\ > h, x Е S1.
Используются общепринятые обозначения: точка (x,y) = (x(k + 1),y(k + 1)) есть образ точки (x,y) = (x(k),y(k)). Случай отрицательных значений параметра 6 приводится к исследуемому заменой (x, y) ^ (—x, -y). Отображение (0.1) задано в цилиндрическом фазовом пространстве S1 х Rn. Следовательно, неблуждающие траектории F делятся на колебательные и вращательные. А именно, пусть O* = {x*(k),y*(k),k £ Z} неблуждающая траектория отображения F в накрывающем фазовом пространстве Rn+1. Введем аналог числа вращения Пуанкаре вдоль этой орбиты
г = lim ^. (0.3)
Траекторию O* будем называть колебательной, если г = 0 и вращательной, если г = 0.
И вообще, аттрактор отображения F будем называть колебательным (о-аттрактором), если он не содержит вращательных траекторий, и вращательным (^-аттрактором), если не содержит колебательных. Смешанным аттрактор (m-аттрактор) будет при наличии в нем и вращательных, и колебательных траекторий.
Такие отображения возникают в задачах хаотической динамики [7] и служат математической моделью дискретной системы фазовой автоподстройки частоты с цифровым фильтром высокого порядка [8].
Заметим, что отображение (0.1), записанное в общем виде
{
x = f(x) + LT y; yi = -bg(x) + By),
где Ь = ео1отп (Ьх,..., Ьп), Ь = ео1отп (Ьх,... ,Ьп), В — п х п — матрица с действительными собственными значениями Ах,... , Ап, f (х) = х + 8 — ад(х), линейным невырожденным преобразованием у ^ Ну приводится к нормальной форме вида (0.1).
При использовании свойств ограниченности и отделённости от нуля разрывной производной периодической нелинейности находятся достаточные условия существования сингулярно-гиперболического странного аттрактора.
1. Динамика отображения в вырожденном случае
Рассмотрим простейший случай отображения (0.1) при Аг = 0,г = г, п. При этом по переменным уг происходит абсолютное сжатие, так что за одну итерацию уг = 0,г = г,п и отображение сводится к отображению окружности на себя
X = х + 8 — ад(х) = f (х). (1.1)
Мы опускаем "простой" случай, когда это отображение взаимнооднозначно (1 — ад' > 0) и число вращения является инвариантом.
Очевидно, что все траектории этого отображения неустойчивы при условии
1 - ад' < -1 или 1 - д' > 1 (1.2)
для разных участков монотонности функции д(х).
Эти условия, в этом случае принимаемые за "гиперболичность", в силу (0.2) имеют вид
ак> 2 (1.3)
В результате, мы получаем простое утверждение
Лемма 1. При Хг = 0, г = 1,п при условии (1.3) отображение (0.1) имеет сингулярно-гиперболический аттрактор.
В качестве примера рассмотрим тестовую функцию
{
x, N< V,x + V(mod2*);
, V < x < 2n - v.
n—v '
В этом случае условие (1.2), очевидно, принимает вид
а> 2v, (1.5)
и следовательно отображение (1.1), (1.4) при условии (1.5) имеет сингулярно-гиперболический аттрактор.
Замечание 1. Сингулярность аттрактора, порождаемая разрывами производной в критических точках, приводит к тому, что при увеличении параметра 8 от нуля происходят бифуркации исчезновения колебательных и рождение вращательных траекторий, однако аттрактор при этом остается сингулярно-гиперболическим.
В случае тестовой функции (1.4), при условии
8 = 0; а + V <п; а> 2v (1.6)
траектории сингулярно-гиперболического аттрактора имеют только колебательный тип, а при
8 > а; а > 2v (1.7)
— только вращательный.
2. Диссипативность отображения Г
Для простоты будем рассматривать случай, когда Ьг = Ь,г = 1,п. При 0 < Хг < 1 многомерное отображение Е имеет поглощающую область. Введем норму
1 п
" " £ \Уг\- (2.1)
n
i=1
Справедлива следующая
Теорема 1. Пусть
0 <Xi <A* < 1, i = 1~n. (2.2)
Тогда отображение F диссипативно по переменной y = colomn(y1,...,yn), так что все его траектории из внешности полнотория
( A*b }
D = {||y|| < ,x е S1 }, (2.3)
попадают внутрь него и остаются в нём навсегда.
Доказательство. Найдем область фазового пространства, для которой при любых x е S1 норма ||y|| убывает при каждой итерации F, т.е. выполняется условие
ЦуЦ < М, x е S1. (2.4)
В силу (0.1) (при bi = b) получаем
1 n n ||y|| = -Е Ail - bg(x) + yi\ <-J2(b + Ы). (2.5)
i=1 i=1
Неравенства (2.4), (2.5) выполняются при
A*b
M > , (2.6)
откуда следует, что во внешности полнотория D траектории действительно приближаются к D, попадая в D за конечное число итераций. Неравенство (2.6), кроме того, означает FD С D. □
Следствие 1. Аттракторы отображения F расположены в полноторим D, и их исследование достаточно проводить в области (2.3).
3. Существование аттракторов отображения F
В отображении (0.1), рассматриваемом в области D, оценим сумму, входящую в первое уравнение
l£yil <nM < An = г (3.1)
Введем одномерные отображения сравнения
F+(-) вида
X = x + 8 — ag(x) ± y = f (x) ± y, (3.2)
мажорирующие координаты x траекторий отображения O = {x(k),y(k),k е Z +}. При этом любой образ x, определенный отображением (0.1), ограничен в силу неравенств
f (x(k)) — Y < x(k + 1) < f (x(k)) + y. (3.3)
Неравенства (3.3) позволяют эффективно использовать отображения сравнения (3.2) при оценке качественного поведения траекторий отображения F.
Теорема 2. 1. Если оба одномерных отображения (3.2) имеют только о-аттрмкторгы, то отображение F имеет только о-аттр>актор.
2. Если оба отображения (3.2) имеют 'только ф-аттракторы, т.е. при 6 > a + y, то отображение F также имеет только ф-аттрактор.
Доказательство. 1. Очевидно, что траектории о-аттракторов отображений окружности FO± = {x±(k), k £ Z} ограничены.
x1 < x± < x2, k £ Z; ,
1 < 2 (3.4)
x2 — x1 < 2п.
Тогда в силу неравенств (3.3) x — координаты траекторий аттрактора отображения O* = {x*(k),y*(k),k £ Z} удовлетворяют неравенству (3.4), а y*(k) £ D, т.е. это о-аттрактор отображения F.
2. Оба отображения сравнения (3.2) имеют только вращательный аттракторы при условии
x = f (x) ± y > x, Vx £ S1, (3.5)
которое справедливо при 6 > a + y.
□
Следствие 2. Пусть 6 = 0 и a = a1 соответствует взаимной однозначности отображений (3.2), 1 — a1g'(x) > 0, Vx £ S1. Тогда при a > a1} с увеличением параметра a о-аттрмкторы отобрмжений Fа следовательно, и отобрмже-ния F, становятся т-аттракторами при бифуркации рождения вращательных траекторий.
Такие изменения легко проследить на примере F с тестовой функцией (1.4).
Пример 1. Пусть g(x) есть тестовая функция (1.4). Анализ отображений F+(_) в этом случае приводит к следующему утверждению.
Теорема 3. При условии
(6 — Y )(п — v) v — a + 6 — y — —п — --—--
a (3 6)
(6 + Y )(п — v) [ )
— v + a + 6 + y < п — --—--
a
Отображение F имеет о-аттрактор в области D\I, где
I = |x \ f (v) — y < x < f (—v)+ y },
Доказательство. Второе условие есть нервенство f (—v) < f (x2), где x 2 — координата неподвижной точки отображения F x2 > v.
Первое условие — неравенство f (v) > f (x1), где x1 - координата неподвижной точки отображения Fx1 < —v. Эти условия есть условия того, что интервал I лежит внутри интервала (x1,x2), обеспечивая отсутствие ф-траекторий.
Рис. 1. Аттрактор колебательного типа.
Рис. 1 иллюстрирует интервал, на котором расположены х-координаты о-аттрактора отображения Е.
Рис. 2. Аттрактор смешанного типа.
При 8 > 0 с увеличением параметра а первым нарушается второе неравенство в (3.6). При этом у отображения могут появиться ^-траектории.
При условии f (—V) — 7 > х2 (рис. 2) отображение Е имеет т-аттрактор, содержащий одновременно и о-траектории, и ф-траектории. Рис. 3 иллюстрирует случай 8 > а + 7, когда аттрактор Е не имеет о-траеторий. □
4. Гиперболичность аттракторов
Пусть О = {х(к),у(к),к £ Z} есть неблуждающая траектория отображения Е. Свойства устойчивости (гиперболичности) траектории О определяются уравнением в вариациях вида
п
£ (к + 1) = [1 — аБ (к)]£(к) + ^ Пг(к)
г=1
п(к + 1) = Аг[—ЬБ (к)£(к) + Пг],
=1 i = 1,n, (4.1)
где 5 = д'(х(к)).
Для удобства обозначим р = 1 — аБ (к), д = —ЬБ (к), £ = £ (к + 1), П = п(к + 1), и перепишем (4.1) в виде отображения Т
п
£ = Р£ + ^ Пг
г=х г = 1,п, (4.2)
П= АМ(к) + Пг}, Введем в каждой точке (х(к),у(к)) £ О одинаковые конусы
Ku = {£, Vi \ Vi = N <X,i =1,n}, (4.3)
Найдем образ образующей линии
I = [щ = С е = Т,П}, (4.4)
Подставляя (4.4) в (4.3) получаем Е1 в параметрическом виде
п
i=1 г = 1,п. (4.5) fji = Xi(q +
Из (4.5) получаем образ Fl
Vi = (4.6)
где
k(q + ai)
а = —-, i =1,n. (4.7)
P + 2^i=i ai
Первое уравнение в (4.5) есть одномерное отображение переменной £.
Лемма 2. Отображение £ ^ £ в точках конуса Ku является растягивающим, т.е.
n
\P + ail > 1 + £, £> 0, (4.8)
i=1
при условии
ah > xn + £ + 2 (4.9)
Доказательство. Неравенство (4.8), очевидно, переписывается в виде двух неравенств
1 — ag'(x(k)) + ^^ ai < — 1 — е;
i=1 n
1 — ag'(x(k)) + ^^ ai > 1 + е
i=1 (4.10)
(x(k)) + а
i=1
Используя условие \g'\ > h в (0.2) и границы конуса Ku \ai\ < х, взамен (4.10) получаем неравенство (4.9) и неравенство ah > £ + хп, справедливое при условии (4.9). □
Теорема 4. При условии
1 п А*bn . .. , ah> 2 + £ + ---, \g'\ >h, £> 0,
1 - А* (4.11)
А* = min Ai
i£[1,n]
аттракторы отображения F сингулярно-гиперболические.
Доказательство. Формула (4.7) задает отображение аг ^ аг, г = 1,п координат направляющих векторов прямых (4.4) и (4.6). Перепишем (4.6) в виде
cii
kg
+ (-
k
-)®i
Р + ЕГ=1 аг р + ЕГ=1 аг
В силу неравенства (4.10) и условия \д(х)\ < 1 в (0.2) имеем
kiQ
Р + Е
n
i=1 ai
< A*b,
k
n
p + E i=1 ai
< k
1 , n.
Тогда область, где \аг\ < \аг\, г = 1,п определяется неравенствами
k*b + k*\ai\ < \а\, i =1,n.
Следовательно, область
ы <
X*b
(4.12)
(4.13)
(4.14)
(4.15)
есть притягивающая область отображения (4.12). В результате, полагая % = в неравенствах (4.15), определяющих конус Ки, получаем его инвариантность, ТКи £ Ки. Тогда по лемме 2, где условие (4.9) заменено на условие (4.11), отображение Т в инвариантном конусе Ки является растягивающим. В силу свойств отображения (4.12) в качестве устойчивого конуса Кя может быть выбрана внешность конуса Ки.
□
Заключение
В работе проведено исследование многомерного отображения Е вида
п
(х, уг) ^ (х + 8 — ад(х) + ^2 уг, Хг(—Ьгд(х)+ уг)), г = 1,п
г=1
с одной кусочно-гладкой периодической нелинейностью д(х) .
Получены условия, при которых отображение Е имеет аттракторы, лежащие в полнотории
г \*Ь }
О = {||у|| < , х £ Б1 }.
Даны критерии, определяющие колебательный и вращательный тип аттракторов. Доказана теорема, указывающая область параметров, при которых аттракторы как колебательные, так и вращательные, являются сингулярно-гиперболическими. При этих условиях динамическое поведение траекторий отображения представляет собой пример динамического хаоса.
Список цитируемых источников
1. Аносов, Д. В. Геодезические потоки на замкнутых римановых многообразиях отрицательной кривизны. Тр. МИАН СССР 90, 3-210 (1967).
2. Afraimovich, V. S., Bykov, V. V. & Shilnikov, L. P. On attracting structurally unstable limit sets of Lorenz attractor type. Trans. Mosc. Math. Soc. 44, 153-216 (1982).
3. Robinson, C. Homoclinic bifurcation to a transitive attractor of Lorenz type. Nonlinearity Anal. 23, 1255-1268 (1989).
4. Lozi, R. Un attracteur de Henon. J. Physique 39, 9-10 (1978).
5. Belykh, V. Chaotic and strange attractors of a two-dimensional map. Matematicheski Sbornik 186:3, 3-18 (1995).
6. Henon, M. A two-dimensional mapping with a strange attractor. Communications in Mathematical Physics50:1, 69-77 (1976).
7. Заславский, Г. М., Чириков, Б. В. Стохастическая неустойчивость нелинейных колебаний. Успехи физических наук №105, 3-39 (1971).
8. Белых, В. Н. Модели дискретных СФС и их исследование. В кн. Системы фазовой синхронизации. Под ред. В.В.Шахгильдяна, Л.Н.Белюстиной. М.: Радио и связь, 1982, 161-162.
9. Anosov, В., Sinai, J. Some smooth dynamical systems. Russian Math. Surveys 22:5, 107172 (1967).
10. Arnold, V., Afraimovich, V., Il'yashenko, Yu.& Shilnikov, L. Theory of bifurcations. VINITI Modern Problems of Mathematics 5, 5-218 (1986).
11. Белых, В. Н. Странный аттрактор. В кн. Большая российская энциклопедия, Т. 31, 285-286 (2016).
12. Sataev, E. Invariant measures for hyperbolic maps with singularities. Russian Math. Surveys 47:1, 147-204 (1992).
13. Belykh, V., Ukrainsky, B. Hyperbolic attractor of a continuous piecewise smooth 2-dimensional map. Sbornik Nauchnyh Trudov "Modelirovanie i optimizacia slozhnyh sistem" N. Novgorod: VGAVT, 1997.
14. Belykh, V., Mordvinkina, I. and Ukrainsky B. Multidimensional Lurie systems and Henon maps: Smale's horseshoes and bifurcations. Тезисы докладов Международной конференции «Динамика, бифуркации и странные аттракторы» посвященной памяти Л. П. Шильникова (1934-2011) N. Novgorod, 2013.
15. Belykh, V., Komrakov, N. and Ukrainsky B. Hyperbolic attractors in a family of multidimensional maps with cusp-points. Proc. of int. conf. "Progress in nonlinear science" dedicated to the 100-th anniversary of A. Andronov. N. Novgorod, 2013.
16. Белых, В. Н. Элементарное введение в качественную теорию и теорию бифуркаций динамических систем. Соросовский образовательный журнал №1, 115-121 (1997).
Получена 05.11.2018
