CC BY
17Нелинейная динамика. 2014. Т. 10. № 3. С. 265-277. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru
УДК: 517.9+534
М8С 2010: 37D20, 37D45, 37Ь30, 37Ь45
Гиперболический хаос в системах с параметрическим возбуждением паттернов
стоячих волн
В. П. Круглов, А.С.Кузнецов, С.П.Кузнецов
Показана возможность реализации аттракторов типа Смейла-Вильямса с разной кратностью растяжения угловой координаты п = 3, 5, 7, 9,11 у отображений, описывающих эволюцию параметрически возбуждаемых паттернов стоячих волн на нелинейной струне за период модуляции накачки при попеременном возбуждении мод с отношением длин волн 1 : п.
Ключевые слова: параметрические колебания, струна, аттрактор, хаос, показатель Ляпунова
1. Введение
Недавно предложен подход к построению систем со структурно устойчивым гиперболическим хаосом, основанный на попеременном возбуждении в нелинейной среде длинноволновых и коротковолновых паттернов, пространственная фаза которых за характерный
Получено 05 сентября 2014 года После доработки 18 сентября 2014 года
Работа выполнена при поддержке гранта РФФИ №14-02-31162 (В. П. К. и А. С. К.).
Круглов Вячеслав Павлович kruglovyacheslav@gmail.com Кузнецов Сергей Петрович spkuz@rambler.ru
Саратовский филиал Института радиотехники и электроники им. В. А. Котельникова РАН 410019, Россия, г. Саратов, ул. Зеленая, д. 38
Саратовский государственный университет им. Н.Г.Чернышевского 410012, Россия, г. Саратов, ул. Астраханская, д. 83
Кузнецов Алексей Сергеевич alek.art@mail.ru
Саратовский филиал Института радиотехники и электроники им. В. А. Котельникова РАН 410019, Россия, г. Саратов, ул. Зеленая, д. 38
временной период изменяется в соответствии с растягивающим отображением окружности. Этот подход был реализован для систем, сконструированных на основе модификации уравнения типа Свифта-Хоенберга [1, 2], а также для параметрически возбуждаемых стоячих волн на нелинейной струне [3, 4]. Было показано, что в этих системах реализуются аттракторы типа соленоида Смейла-Вильямса.
|.Л = 3 1к = 5 \к = 7
Рис. 1. Преобразование области в форме тора за один шаг эволюции в дискретном времени (вверху) и предельный объект-аттрактор в виде соленоида (внизу) при кратности растяжения по угловой координате к = 3, 5 и 7.
Наиболее распространенный в литературе и учебных курсах вариант построения соленоида Смейла-Вильямса заключается в том, что в трехмерном пространстве рассматривается область в виде тора, которая за один шаг эволюции претерпевает продольное растяжение, поперечное сжатие и складывание двойной петлей, оказываясь расположенной внутри исходного тора. С увеличением числа шагов количество витков трансформируемого объекта стремится к бесконечности, и возникает аттрактор-соленоид, имеющий в поперечном направлении структуру типа канторова множества. Такие аттракторы могут иметь место также и в пространстве размерности выше трех.
Очевидное обобщение конструкции состоит в том, что растяжение вдоль витков можно производить не в два, а в большее число раз (рис. 1). В частности, в работах [1, 3, 4] речь идет о трехкратном растяжении, причем роль угловой координаты, отсчитываемой вдоль витков соленоида, вложенного в бесконечномерное фазовое пространство, играет пространственная фаза паттернов, возбуждающихся на последовательных этапах эволюции нелинейной среды. В настоящей работе будут рассмотрены варианты систем с параметрическим возбуждением стоячих волн, отличающиеся разной кратностью растяжения вдоль витков аттрактора Смейла-Вильямса на одном шаге его построения.
2. Модельное уравнение и механизм функционирования системы
Уравнение в частных производных, описывающее распространение волн в одномерной среде в линейном приближении, имеет вид
рЦ = сЦ (2.1)
сМ2 дх2
и может использоваться во многих физических ситуациях, одна из которых относится к механическим колебаниям струны. В этом случае переменная y(x,t) отвечает поперечному смещению струны в точке с координатой x в момент t, р — линейная плотность струны (масса на единицу длины), G — сила продольного натяжения струны. Для однородной струны при постоянном натяжении величина с = y/G/p определяет скорость распространения волн.
В классическом опыте Мельде [5, 6] при периодическом изменении во времени силы натяжения струны с частотой 2ио, где ио отвечает одной из мод, возбуждаются параметрические колебания на этой моде (рис. 2). Начальная стадия параметрического возбуждения (пока колебания можно полагать линейными) описывается уравнением вида (2.1), где коэффициент G является функцией времени (скажем, G = Т0(1 + asin2uot), 0 < a = const < 1). Колебания силы натяжения выступают в качестве накачки — источника энергии, поддерживающего параметрические колебания.
2ш0
Рис. 2. Опыт Мельде (1859 год): колебания камертона на частоте 2шо вызывают периодическое изменение силы натяжения струны, благодаря чему происходит возбуждение параметрических колебаний на частоте шо.
Пусть в уравнении (2.1) сила натяжения колеблется по закону
G(t) = То • [1 + A(t) sin 2uot + B(t) sin 2nu0t], (2.2)
где коэффициенты A, B изменяются во времени с некоторым периодом Т, поочередно оказываясь по величине большими или малыми, близкими к нулю. Конкретно, мы задаем
A(t) = a cos2 (nt/T), B (t) = b sin2(nt/T), (2.3)
где неотрицательные постоянные величины a, b подчинены условию a + b < 1.
Распределение массы по струне будем считать слабонеоднородным, зависящим от пространственной координаты по закону
p(x) = р0(1 + е sin mk0x), (2.4)
где ко = Wo/co, Со = у/То/ро и m — некоторое целое число.
Кроме того, добавим в правую часть уравнения член вида —(а + ¡3u2)dy/dt, где параметр а отвечает за линейную, а в — за нелинейную диссипацию при колебаниях струны. Присутствие нелинейной диссипации необходимо для стабилизации параметрической неустойчивости на некотором уровне амплитуды колебаний, однако для рассматриваемой здесь системы весьма существенным является еще и то обстоятельство, что кубическая нелинейность обеспечивает возникновение третьей и других нечетных гармоник при колебательно-волновых движениях.
Используя нормировку переменных и параметров, такую, что Со = 1, ко = ио, в = 1, приходим к системе, описываемой уравнением в частных производных следующего вида:
/л \ д2у , 2. ду
(1 + ¿sinткох)—т = -{а + у - 7у+ dt2 dt
+ (l + a cos2 ^ sin2w0í + 6 sin2 ^ sin2??,w0í]> (2.5)
[ T T ) дх2
Будем считать наложенными граничные условия периодичности:
у(х + L, t) = у(х, t); (2.6)
длину L примем равной целому числу длин волн с волновым числом ко: L = 2nN/ko. Добавление в уравнение линейного члена вида —jу в системе с периодическими граничными условиями обеспечивает затухание колебательно-волновых возмущений с нулевым волновым числом.
Механизм функционирования системы состоит в следующем.
Предположим, что при накачке на частоте 2uo в системе возбуждена стоячая волна на частоте Uo с волновым числом ко, у которой расположение узлов и пучностей характеризуется фазовой постоянной 9, так что в грубом приближении
у ~ cos u0t sin(k0x + 9). (2.7)
Амплитуда волны стабилизируется на конечном уровне благодаря нелинейной диссипации, и из-за ее присутствия колебательно-волновое движение будет иметь также составляющие на третьей, пятой, седьмой и других нечетных гармониках. Пространственная фаза гармоники на частоте пио будет, очевидно, даваться величиной n9, где п — нечетное число.
На следующем этапе временной эволюции амплитуда накачки на частоте 2ио становится малой, и колебания на частоте ио затухают. Амплитуда накачки на частоте 2пио, с другой стороны, с какого-то момента оказывается достаточной для развития параметрической неустойчивости стоячей волны с частотой пио и волновым числом пко. Эта волна формируется из начального возмущения, отвечающего n-й гармонике волны, существовавшей на предыдущей стадии, так что пространственная фаза будет определяться величиной п9.
Затем вновь возникает ситуация, когда интенсивная накачка имеет место на частоте 2ио. Затравочное колебательно-волновое возбуждение с волновым числом ко обеспечивается теперь благодаря комбинации возмущения у ~ sin(n^x + п9), оставшегося от предыдущей стадии процесса, и пространственно-периодического возмущения, обусловленного неоднородным распределением массы р(х) и пропорционального е sin ткох.
Чтобы присутствовала составляющая, находящаяся в пространственном резонансе с параметрически возбуждаемой модой с волновым числом ко, величину т нужно выбрать на единицу больше или меньше, чем п. В самом деле, из тригонометрического тождества
sin ткох s'm(nkox + пв) = ^ [cos((m — п)кох — пв) — cos((m + п)кох + пв)] (2.8)
видно, что при т = п + 1 резонансным является первый член в правой части (2.8), которому отвечает волновое число ко. Поэтому новая величина пространственного фазового сдвига 9 получается как результат применения п-кратно растягивающего отображения окружности
вида 9' = —ив + const. Если же m = n — 1, то резонансным является второй член, и отображение имеет вид в' = ив + const.1
Знак минус в отображении соответствует случаю, когда направление обхода вокруг витков соленоида инвертируется на каждом шаге построения, а знак плюс — случаю, когда направления совпадают. И в том, и в другом случае для превышающих единицу чисел и (рассматриваться будут только нечетные n) это отображение с хаотической динамикой, характеризуемое положительным показателем Ляпунова Л = ln n.
При наличии сжатия в фазовом пространстве по остальным направлениям, для отображения, описывающего изменение состояния системы за период модуляции накачки, на этом основании можно предположить присутствие аттрактора типа соленоида Смейла-Вильямса.
3. Численные результаты
Для численного решения уравнения (2.5) используем схему «крест», имеющую второй порядок аппроксимации по пространственному и временному шагу и являющуюся явной [7]. Решение находится в узлах сетки с шагом по координате Н и по времени т. Аппроксимируя производные разностными соотношениями, на основании (2.5) запишем
[1 + е sin(mk0ih)]
y
j+1 — 2yj + yj-1
a + (yj )2
y
j+1
y
j-1
+ [1 + A(jr) sin(2uojr) + B(jr) sin(2nuojr)]
2t
yj+i— 2yj + yj-i
h
2
+
(3.1)
j
Ж.
Нижний индекс г нумерует узлы сетки по пространственной координате, а верхний индекс
] — по времени. Решение можно искать слой за слоем по времени, сохраняя в памяти ком-
2
пьютера на каждом этапе значения переменной, относящиеся к двум предыдущим слоям.
В работах [3, 4] было принято п = 3, т = 4 и показано, что пространственная фаза попеременно возбуждающихся длинноволновых и коротковолновых паттернов стоячих волн за период модуляции накачки подчиняется троекратно растягивающему отображению окружности с обратным направлением обхода для фазы-образа. При другом выборе п и т оказывается возможным обеспечить иную кратность растяжения, как с прямым, так и с обратным направлением обхода (см. таблицу 1).
2
Таблица 1
п (кратность растяжения) 3 5 5 7 7 9 9 И и
m 4 4 6 6 8 8 10 10 12
прямой обход (П)/ обратный обход (О) О П О П О П О п О
1 Присутствующую в этих соотношениях константу можно устранить выбором начала отсчета фазы.
2Известно, что для однородной струны при постоянном натяжении условием вычислительной устойчивости схемы служит соотношение с0(т/Н)2 < 1; можно полагать, что для нашей задачи схема будет устойчива, когда т2Н-2 [1 + Лат)&т(2шо]т) + Ват)ъ\\\(2пшо]т)] < 1 для всех ]. Это предположение вполне согласуется с результатами расчетов.
Рис. 3. Пространственно-временные диаграммы, иллюстрирующие параметрические колебания системы (2.5) с периодическими граничными условиями при различных т и п при задании параметров согласно (3.2) на четырех полных периодах модуляции накачки.
При численном решении уравнений были заданы следующие параметры:
ш0 = 2^, к0 = 2^, Ь = 1,
а = 0.4, Ь = 0.2, е = 0.2, а = 0.4, 7 = 0.03,
Т = 60 (п = 3), Т = 80 (п = 5), Т = 100 (п = 7), Т = 120 (п = 9), Т = 140 (п = 11).
(3.2)
При численном моделировании в системе наблюдается поочередное возбуждение длинноволновых и коротковолновых структур, пространственная фаза которых от одного периода модуляции к другому претерпевает растягивающее отображение кратности п.
На рисунке 3 представлены пространственно-временные диаграммы, иллюстрирующие процесс функционирования системы для разных величин п и т. Форма струны показана в моменты времени, кратные периоду быстрых колебаний, что делает наглядным визуальное представление эволюции огибающей паттернов стоячих волн. Как можно видеть, в системе поочередно возбуждаются длинноволновые и коротковолновые структуры, притом пространственная фаза их меняется от одного периода модуляции накачки к другому. Фаза колебаний стоячих волн во времени, с другой стороны, жестко привязана к фазе накачки.
В численных расчетах можно проверить, что пространственная фаза стоячих волн за период модуляции накачки действительно претерпевает п-кратно растягивающее отображение окружности. Для построения итерационных диаграмм выполняем решение уравнения (2.5) на большом числе периодов модуляции и в моменты времени t = tn = nT вычисляем пространственные фазы
On = arg[y(0, t) + iyx(0, t)/ko] (3.3)
с представлением результатов в координатах (On,On+i) (рис. 4).
Вид графиков соответствует растягивающим отображениям окружности с числом ветвей п, что согласуется с описанным выше механизмом. В самом деле, расположение ветвей таково, что однократный обход окружности для прообраза, то есть изменение аргумента на 2п, соответствует n-кратному обходу для образа (в прямом или обратном направлении в соответствии с таблицей 1), с изменением фазы на величину ±2пп.
На рисунке 5 показаны портреты аттракторов в стробоскопическом сечении в проекции на плоскость двух переменных — смещения струны в некоторой точке и пространственной производной в этой точке. Поперечная фрактальная структура не различима из-за сильного поперечного сжатия, и аттрактор, хотя и представляет собой соленоид, визуально выглядит просто как замкнутая кривая. Расположение изображающих точек на последовательных шагах итераций отображения Пуанкаре описывается отображением окружности с троекратным и семикратным растяжением (соответственно, для левой и правой картинки: см. диаграммы (a) и (d) на рис. 4).
Для вычисления показателей Ляпунова используем метод Бенеттина [8-10], адаптированный к рассматриваемой распределенной системе. Его содержание состоит в следующем. Наряду с уравнением (2.5) одновременно решаем численно комплект уравнений в вариациях
/л ■ i \д2у , 2. dy ~ду _
(1 +tsm тк0х)—- = -{а + у )-щ - 2УУ-^ ~1У+
_ (3.4)
+ (1 + a cos2 sin 2Lü0t + b sin2 ^ sin 'IniüotX —| [ T T ) dx2
в количестве, равном числу показателей S, которое требуется найти. При этом используются такие же периодические граничные условия, как для уравнения (2.5), и тот же конечно-разностный метод. По истечении каждого очередного периода модуляции T производится ортогонализация векторов возмущения ys, s = 1,...,S, по схеме Грама-Шмидта, что дает набор векторов yS. После этого для дальнейшего интегрирования уравнений за исходные берутся переопределенные нормированные на единицу векторы возмущения y° = yS/||yS||.
Рис. 4. Диаграммы для пространственных фаз при задании параметров в системе (2.5) согласно (3.2): (а) п = 3, т = 4, (Ь) п = 5, т = 4, (о) п = 5, т = 6, (^ п = 7, т = 6, (в) п = 7, т = 8, (£) п = 9, т = 8, п = 9, т = 10, (Ь) п =11, т = 10.
Рис. 5. Портреты аттракторов в стробоскопическом сечении в проекции на плоскость, где по осям отложены поперечное смещение и пространственная производная в фиксированной точке струны для т = 3, п = 4 (а), т = 7, п = 6 (Ь). Остальные параметры заданы согласно (3.2).
Тогда показатели Ляпунова, определенные для отображения за период, получаются из соотношений
NN *—'
где N — полное число периодов, на протяжении которых проводится вычисление.
В таблице 2 приводятся первые три показателя Ляпунова, рассчитанные для систем с разными значениями т и п при задании параметров согласно (3.2). Наибольший показатель, как и ожидалось, во всех рассмотренных случаях близок к величине 1п п. Остальные показатели отрицательные. Они отвечают за приближение траекторий к аттрактору, представляющему собой разновидность соленоида Смейла-Вильямса, который получается при трансформации исходной тороидальной области во вложенную п-кратную петлю, претерпевшую продольное растяжение и достаточно сильное сжатие в направлении поперек витков.
Таблица 2
1п п А1 л2 Л3
п = 3, т = 4, Т = 60 1.0986 1.085 -2.028 -18.9
п = 5, т = 4, Т = 80 1.6094 1.515 -2.150 -31.9
п = 5, л? = 6, Т = 80 1.603 -2.242 -32.1
п = 7, л? = 6, Т = 100 1.9459 1.927 -3.273 -30.6
п = 7, т = 8, Т = 100 1.942 -3.223 -32.1
п = 9, т = 8, Т = 120 2.1972 2.194 -3.529 -31.0
п = 9, т = 10, Т = 120 2.185 -3.518 -30.9
п = 11, т = 10, Т = 140 2.3979 2.338 -5.669 -28.5
п = И, т = 12, Т = 140 2.394 -5.707 -27.9
4. Конечномерная модель и проверка гиперболичности
Для уравнений (2.5) в случае п = 3, т = 4, когда реализуется трехкратное растяжение для угловой переменной, в работе [11] показано, как можно с использованием разложения по модам свести динамику к конечномерному описанию в рамках системы обыкновенных дифференциальных уравнений. (Аналогичным образом можно получить конечномерное описание и для вариантов с большими коэффициентами растяжения угловой переменной, то есть для п = 5, 7, 9,11, но в этом случае надо учитывать гармоники разложения в пространственный ряд Фурье с номерами до п включительно, и система уравнений получается очень громоздкой.)
Поскольку для случая п = 3, т = 4 параметрическое возбуждение имеет место на волновых числах ko и 3ko, естественно искать решение уравнения (2.5) в виде
y(x, t) = u1(t) cos k0x + v1(t) sink0x + u3(t) cos 3k0x + v3(t) sin3k0x. (4.1)
Подстановка этого выражения в (2.5) и разложение обеих частей уравнения в ряд Фурье с учетом гармоник того же порядка приводят к системе уравнений для коэффициентов
üi + \dh = /ъ й3 + У vi = /з, vi + = gi, v3 + \eüi = д3, (4.2)
где
(nt nt \
1 + a.cos2 — sin2woí + b sin2 — sin6woí j ui — 7«i — айi —
- ^[(|г/2 + uius + u2 + ^vf + v\v3 + vl)Ú! + + 4«i«3 - vf)u3+ + (uivi - U3V1 + uiv3)Vi + (uivi + 2uiV3)-V3],
2 1 2 nt ,2 nt \
g 1 = —w0 I 1 + a eos — sin 2woí + b sin — sin owoí I vi — 7^1 —
- m')i - - v\v3 + v2 + |г/2 - И1И3 + u3)i'i + |(г/2 + 4i'iv3 - vf)v3+
+ (uivi - u3vi + uiv3)ui + (-uivi + 2u3vi)uÍ3],
2 1 2 nt ,2 nt \
/3 = —9w0 I 1 + a eos — sin 2woí + b sin — sin owoí I «з — 7И3 — au3
~ 3 [(«i + 4«i«3 - vf )üi + 2(-ttiVi + 2tt3vi)vi + (2u\ + 2v2 + Зг/| + г>| )ü3 + 2«зг>зг>з] ,
(nt nt \
1 + a eos2 — sin 2woí + b sin2 — sin 6woí j V3 — 7V3 — сшз —
- j [2(wivi + 2ttiv3)üi + (г/2 + 4viv3 - vf )i'i + 2u3v3ú3 + (2«f + 2vf + г/| + Зг|)г>3].
Уравнения можно переписать, выразив вторые производные явно:
/1 - \e9¿ .. /з -
Ul = 1 _ 1,2 > = ! _ 1,2 .
4 4 ,
4 4 (4.3)
д 1 - 2е/з .. 93 - 2еh
Vi = -V-3 = -——.
1 - Г'2 1 - \£2
Согласно [11], результаты расчетов с использованием конечномерной модели находятся в хорошем соответствии с теми, что получаются при численном решении уравнения в частных производных (2.5) в случае п = 3, т = 4 [3, 4]. Этим можно воспользоваться для проверки гиперболичности аттрактора с помощью известного из литературы метода.
Идея тестирования гиперболичности на основе «критерия углов» была выдвинута в работах [12] и [13], а в наиболее простой и удобной для проведения вычислений форме методика представлена в работе [14]. При наличии одного неустойчивого направления (один положительный показатель Ляпунова) она состоит в следующем.
1500
М
0-1--
0 в тг/2
Рис. 6. Гистограмма распределения углов пересечения локальных устойчивых и неустойчивых многообразий на аттракторе, полученная обработкой данных численного решения применительно к отображению Пуанкаре системы уравнений (4.3) по методике, описанной в тексте. Параметры шо = 2п, а = 0.4, Ь = 0.2, е = 0.2, а = 0.4, 7 = 0.03, Т = 60. Отсутствие углов, близких к нулю, подтверждает гиперболичность аттрактора.
Сначала на достаточно большом временном интервале проводится решение системы уравнений (4.3), которую формально запишем в виде x = F(x,t), где x — восьмимерный вектор, а функция F определяется видом правых частей уравнений. Затем линеаризованное уравнение для вектора возмущения X = F'(x(t),t)X интегрируется вдоль найденной траектории x(t) с нормировкой на единицу вектора X на каждом шаге n отображения Пуанкаре. В результате получаем набор единичных векторов {хп}. Далее, в соответствии с идеей работы [14], вдоль той же самой опорной траектории проводится интегрирование в обратном
гп
времени линейного уравнения U = — [F'(x(t),t)] T и, где верхний индекс T означает матричное сопряжение. Это дает набор векторов {un}, определяющих ортогональное дополнение к устойчивому подпространству векторов возмущения, которые тоже нормируем на единицу. Теперь для оценки угла а между одномерным неустойчивым подпространством и устойчивым подпространством (семимерным в нашем случае) можно найти угол f3n Е [0,п/2] между векторами хп и un при n-ом прохождении сечения Пуанкаре cos f3n = |un(t) • xn(t)| и положить окончательно 9n = п/2 — f3n.
На основе обработки полученных данных анализируется статистика распределение углов 9. Если обнаруживается отличная от нуля вероятность нулевых углов, то это говорит о нарушении гиперболичности из-за присутствия касаний устойчивых и неустойчивых многообразий. Если же распределение отделено от нуля, то тест подтверждает гиперболичность.
На рисунке 6 показана гистограмма распределения углов между устойчивыми и неустойчивыми подпространствами, полученная численно для указанного в пояснении набора значений параметров системы (4.3). Распределение с очевидностью располагается вдали от нулевых значений углов 9, то есть гиперболичность аттрактора подтверждается.
5. Заключение
Представленное в этой статье исследование расширяет круг параметрически возбуждаемых распределенных систем, демонстрирующих гиперболический хаос, в сравнении с предыдущими работами [3, 4], где были представлены первые модели такого рода. Теперь мы вводим в рассмотрение целое семейство систем с кольцевой геометрией (периодические граничные условия), в которых благодаря периодической модуляции накачки реализуется хаотическая динамика, отвечающая равномерно гиперболическим аттракторам типа соленоидов Смейла-Вильямса с разной кратностью растяжения угловой координаты, определяемых нечетными числами n = 3, 5, 7, 9,11. Динамика состоит в попеременном рождении и затухании пространственных паттернов с отношением волновых чисел 1 : n. Аттракторы типа Смейла-Вильямса возникают в связи с тем, что пространственные фазы этих паттернов на каждом следующем периоде модуляции накачки трансформируются в соответствии с n-кратно растягивающим отображением окружности. Для ситуации трехкратного растяжения гиперболичность аттрактора подтверждена результатами компьютерной проверки, показавшей отсутствие касаний для устойчивых и неустойчивых многообразий орбит на аттракторе для редуцированной конечномерной модели.
Как постулируется в теории колебаний, грубость (структурная устойчивость) служит принципиальным основанием для выделения значимых для практики и подлежащих первоочередному теоретическому анализу колебательных систем [15]. Для любых возможных практических приложений хаоса (коммуникация [16], шумовая локация [17], генерация случайных чисел [18]) кажется естественным отдать предпочтение именно генераторам грубого хаоса, нечувствительным к параметрам, деталям изготовления устройств, техническим флуктуациям и шумам. Для построения генераторов грубого хаоса принцип параметрического возбуждения представляется одним из самых перспективных при практической реализации, поскольку подразумевает использование простых нелинейных элементов или сред и открывает широкие возможности для манипуляции частотами и фазами колебательно-волновых составляющих. Поскольку параметрическое возбуждение и усиление колебаний и волн широко известно и применяется в механике [19, 20], акустике [21, 22], радиофизике [23, 24], нелинейной оптике [25], естественно ожидать, что возможность генерации грубого хаоса в параметрических системах с модуляцией накачки будет представлять интерес для этих областей.
Список литературы
[1] Kuptsov P. V., Kuznetsov S. P., Pikovsky A. Hyperbolic chaos of Turing patterns // Phys. Rev. Lett., 2012, vol. 108, no. 19, 194101, 4 pp.
[2] Kruglov V. P., Kuznetsov S.P., Pikovsky A. Attractor of Smale-Williams type in an autonomous distributed system // Regul. Chaotic Dyn., 2014, vol. 19, no. 4, pp. 483-494.
[3] Исаева О. Б., Кузнецов А. С., Кузнецов С. П. Гиперболический хаос при параметрических колебаниях струны // Нелинейная динамика, 2013, т. 9, №1, с. 3-10.
[4] Isaeva O. B., Kuznetsov A. S., Kuznetsov S. P. Hyperbolic chaos of standing wave patterns generated parametrically by a modulated pump source // Phys. Rev. E, 2013, vol. 87, no. 4, 040901(R), 4 pp.
[5] Мандельштам Л. И. Лекции по колебаниям. Москва: АН СССР, 1955. 511с.
[6] Rowland D. R. Parametric resonance and nonlinear string vibrations // Amer. J. Phys., 2004, vol. 72, pp. 758-766.
[7] Калиткин Н. Н. Численные методы. Москва: Наука, 1978. 512 с.
[8] Benettin G., Galgani L., Giorgilli A., Strelcyn J.-M. Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems: A method for computing all of them // Meccanica, 1980, vol. 15, pp. 9-30.
[9] Шустер Г. Детерминированный хаос. Москва: Мир, 1988. 250 c.
[10] Кузнецов С. П. Динамический хаос. 2-е изд. Москва: Физматлит, 2006. 356 с.
[11] Kuznetsov S. P. Some mechanical systems manifesting robust chaos // Nonlinear Dynamics & Mobile Robotics, 2013, vol. 1, no. 1, pp. 3-22.
[12] Lai Y.-Ch., Grebogi C., Yorke J. A., Kan I. How often are chaotic saddles nonhyperbolic? // Nonlinearity, 1993, vol. 6, no. 5, pp. 779-797.
[13] Anishchenko V. S., Kopeikin A. S., Kurths J., Vadivasova T. E., Strelkova G. I. Studying hyperbolicity in chaotic systems // Phys. Lett. A, 2000, vol.270, no. 6, pp. 301-307.
[14] Kuptsov P. V. Fast numerical test of hyperbolic chaos // Phys. Rev. E, 2012, vol. 85, no. 1, 015203(R), 4 pp.
[15] Андронов А. А., Витт А. А., Хайкин С.Э. Теория колебаний. Москва: Наука, 1981. 916 с.
[16] Дмитриев А. С., Панас А. И. Динамический хаос: Новые носители информации для систем связи. Москва: Физматлит, 2002. 252 с.
[17] Lukin K. A. Noise radar technology // Telecomm. Radio Eng., 2001, vol.55, no. 12, pp. 8-16.
[18] Drutarovsky M., Galajda P. A robust chaos-based true random number generator embedded in reconfigurable switched-capacitor hardware // Radioengineering, 2007, vol. 16, no. 3, pp. 120-127.
[19] Пановко Я. Г. Введение в теорию механических колебаний. Москва: Наука, 1971. 256 с.
[20] Борисов А. В., Мамаев И. С. Динамика твердого тела. Москва-Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001. 384 с.
[21] Руденко О. В., Солуян С. И. Теоретические основы нелинейной акустики. Москва: Наука, 1975. 287 с.
[22] Островский Л. А., Папилова И. А., Сутин А. М. Параметрический генератор ультразвука // Письма в ЖЭТФ, 1972, т. 15, №8, с. 456-458.
[23] Люиселл У. Связанные и параметрические колебания в электронике. Москва: ИИЛ, 1963. 352 с.
[24] Каплан А.Е., Кравцов Ю. А., Рылов В. А. Параметрические генераторы и делители частоты. Москва: Советское радио, 1966. 335 с.
[25] Ахманов С. А., Хохлов Р. В. Параметрические усилители и генераторы света // УФН, 1966, т. 88, №. 3, с. 439-460.
Hyperbolic chaos in systems with parametrically excited patterns of standing waves
Vyacheslav P. Kruglov1, Alexey S. Kuznetsov2, Sergey P. Kuznetsov3
i,2,3Kotel'nikov's Institute of Radio-Engineering and Electronics of RAS, Saratov Branch Zelenaya 38, Saratov, 410019 Russia 1,3Saratov State University Astrahanskaya 83, Saratov, 410012, Russia
1kruglovyacheslav@gmail.com, 2alek.art@mail.ru, 3spkuz@rambler.ru
We outline a possibility of implementation of Smale-Williams type attractors with different stretching factors for the angular coordinate, namely, n = 3, 5, 7, 9,11, for the maps describing the evolution of parametrically excited standing wave patterns on a nonlinear string over a period of modulation of pump accompanying by alternate excitation of modes with the wavelength ratios of 1 : n.
MSC 2010: 37D20, 37D45, 37L30, 37L45
Keywords: parametric oscillations, string, attractor, chaos, Lyapunov exponent
Received September 05, 2014, accepted September 18, 2014
Citation: Rus. J. Nonlin. Dyn., 2014, vol. 10, no. 3, pp. 265-277 (Russian)
