Изв. вузов «ПНД», т. 15, № 6, 2007 УДК 517.9
ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ СИНХРОНИЗАЦИИ ДВУХЧАСТОТНЫХ КВАЗИПЕРИОДИЧЕСКИХ КОЛЕБАНИЙ
В.С. Анищенко, С.М. Николаев
Приводится электронная схема автономного генератора двухчастотных квазипериодических колебаний и результаты экспериментальных исследований эффекта синхронизации квазипериодических колебаний внешним двухчастотным сигналом.
Введение
Квазипериодическими колебаниями называют колебания вида
х(*) = / [Ф1(*), Ф2(*),..., ЦиШ,
где фг(í) = г = 1, 2,..., п. Функция х(Ь) имеет период 2п по каждому из
аргументов фг(£):
х(фг + 2п) = х(фг).
Все частоты ю являются независимыми в том смысле, что между ними не существует рациональных соотношений. В фазовом пространстве п-мерному квазипериодическому колебанию отвечает п-мерный эргодический тор.
Исследуем двухчастотные квазипериодические автоколебания, то есть колебания с двумя независимыми частотами, возникающие в автономной нелинейной диссипативной динамической системе. Математически, модель такой системы была введена нами в [1]. В работе [2] представлены результаты численного исследования динамики генератора двухчастотных колебаний описанного в [1]. Установлено, что введенная в рассмотрение модель генератора [2] реализует режим устойчивых двух-частотных колебаний, последовательность бифуркаций удвоения двумерного тора и переход к хаосу [2]. С использованием модели генератора [1] установлено новое явление - эффект синхронизации квазипериодических колебаний в виде захвата числа вращения на двумерном торе [2].
Целью настоящей работы является экспериментальная проверка всех численных результатов, полученных в работах [1,2]. Введем в рассмотрение электронную схему генератора квазипериодических колебаний, покажем ее соответствие математической модели, предложенной в [1], и исследуем экспериментально явление синхронизации двухчастотных квазипериодических колебаний. В итоге демонстрируется полное соответствие экспериментальных и численных результатов.
1. Автономный генератор квазипериодических колебаний
Рассмотрим в качестве исходной известную модель генератора Анищенко-Астахова [3], схема которого представлена на рис. 1, а. Схема включает классический генератор ван дер Поля, в котором введена дополнительная нелинейная инерционная обратная связь. Уравнения генератора представляют собой трехмерную динамическую систему с тремя параметрами
х = тх + у — хх — йх3,
У = —х, (1)
X = —дх + дФ(х).
Первые два уравнения системы (1) описывают генератор ван дер Поля. В этом легко убедиться, положив х = 0 и используя Ф(х) = х2. Как было показано, система (1) реализует переход к хаосу в соответствии с теоремой Шильникова [4] при условии, что нелинейная функция Ф(х) является асимметричной относительно переменной х и задается в виде [3]
{ 1, х > 0,
Ф(х) = I (х)х2, I (х) = { (2)
[ 0, х < 0,
или, например,
Ф(х) = ехр х — 1.
Рис. 1. Схемы: а - генератора Анищенко-Астахова, б - инерционного каскада дополнительной обратной связи генератора Анищенко-Астахова, в - инерционного каскада дополнительной обратной связи генератора квазипериодических колебаний
Условия асимметрии (2) обеспечивают существование в системе (1) особого решения в виде петли сепаратрисы седло-фокуса и, как следствие, реализацию режима спирального хаоса. С физической точки зрения, это достигается за счет самосогласованного воздействия на основной усилитель со стороны обратной связи, заданной третьим уравнением в системе (1). При малых амплитудах сигнала х(Ь) это воздействие незначительно, и система (1) генерирует предельный цикл. С ростом параметра возбуждения т интенсивность колебаний х(Ь) растет, сигнал обратной связи г(Ь) нарастает тоже и более активно начинает управлять коэффициентом усиления основного усилителя. Система реализует последовательность бифуркаций удвоения периода циклов и переход к хаосу.
С целью обеспечения незатухающих двухчастотных колебаний в систему (1) необходимо ввести элемент, характеризуемый собственной частотой, отличающейся от резонансной частоты контура генератора. Одним из возможных способов является использование колебательного контура в цепи дополнительной обратной связи. Необходимо сделать так, чтобы сигнал обратной связи г(Ь) включал колебания независимой частоты, которые будут модулировать коэффициент усиления и обеспечивать квазипериодические автоколебания. Эта идея получила практическое воплощение [1].
Рассмотрим схемы, представленные на рис. 1, б, в. На рис. 1, б показана схема инерционного каскада дополнительной обратной связи генератора Анищенко-Астахова. Каскад представляет собой нелинейный преобразователь и ЯС-фильтр, описываемые одномерным дифференциальным уравнением (третье уравнение в системе (1)). На рис. 1, в представлена схема нового инерционного каскада, который вместо ЯС-фильтра включает диссипативный колебательный контур некоторой резонансной частоты. Уравнения, описывающие схему рис. 1, в, имеют вид
где у - параметр затухания; д - параметр, представляющий нормированную резонансную частоту нового фильтра. Нетрудно убедиться, что уравнения (3) описывают диссипативный колебательный контур в режиме вынужденных колебаний
Экспериментальная реализация генератора квазипериодических колебаний, соответствующая схеме генератора Анищенко-Астахова (см. рис. 1, а), с учетом предложенных изменений (см. рис. 1, в) представляет собой ЯС-генератор с дополнительной обратной связью, состоящей из квадратичного детектора и ЬСЯ-контура. При этом сигнал с дополнительной обратной связи модулирует коэффициент усиления усилителя ЯС-генератора.
На рис. 2. представлена электронная схема генератора квазипериодических колебаний. ЯС-генератор включает ЯС-цепь, состоящую из элементов Е1, С1, К2, С2,
2. Экспериментальная реализация генератора
а = ф,
Ф = —УФ + уФ(х) — да,
(3)
а + у а + да = уФ(х).
(4)
Рис. 2. Электронная схема генератора квазипериодических колебаний
и операционный усилитель и1А, коэффициент усиления которого определяется под-строечным резистором ЕР1 и сопротивлением Ед полевого транзистора Q1. Дополнительная обратная связь состоит из усилителей И1Б и И2Б, однополупериодного квадратичного детектора, в качестве которого используется полевой транзистор Q2, и ЬСЯ-контура, включающего элементы Ь\, Сз и ЕР2. Уравнение ЯС-генератора имеет вид
+__^(2 — ЕР1
ЕС Ез + Ед
)ио —
1
Е2С2
ио = 0,
(5)
где Е = Е1 = Е2 и С = С1 = С2 - элементы ЯС-цепи, и0 - потенциал на неинвер-тирующем входе операционного усилителя И1А. Уравнение дополнительной обратной связи
и + Е2 и 1 + = Е2 /1(ио),
Ь1Сз Ь1
(6)
где и1 - потенциал на выходе ЬСЯ-контура, оказывающий влияние на сопротивление Ед полевого транзистора Ql (см. рис. 2) (аппроксимируется экспоненциальной зависимостью); /1(и0) - функция, описывающая работу однополупериодного квадратичного детектора.
Переходя от уравнений (5) и (6) к уравнениям в первых производных в безразмерном времени Ь = т/(ЕС), получим
X у
х ф
тз — 2)х + у — йх3,
/2 (х)
ф,
= —УФ + 7/1 (х) — дх.
Были введены следующие переобозначения: т = ЕР1, /2(х) = Е1 + Ед, у = ЕСЕР2/Ь1, д = Е2С2/(Ь1С3), параметр й характеризует диссипативную нелинейность усилителя ЯС-генератора и ограничивает амплитуду колебаний в контуре. Как и в случае генератора Анищенко-Астахова (1), в генераторе квазипериодических колебаний (7) с ростом параметра возбуждения т растет амплитуда колебаний и, соответственно, сигнал /2 (х) на выходе цепи дополнительной обратной связи, который начинает оказывать все большее влияние на параметр возбуждения т. В результате система претерпевает бифуркацию мягкого рождения тора, который характеризуется двумя основными частотами, одна из которых определяется ЯС-цепью генератора, а вторая - ЬСЯ-контуром в дополнительной обратной связи. В ходе экспериментов было установлено, что помимо бифуркаций рождения тора и перехода к хаосу через его разрушение, предложенная установка (см. рис. 2.) демонстрирует также последовательность бифуркаций удвоения квазипериодических колебаний [1].
3. Синхронизация квазипериодических колебаний
Проведенные ранее численные эксперименты по внешней и взаимной синхронизации двух генераторов квазипериодических колебаний [2] показали, что область захвата частот модуляции генераторов лежит внутри клюва синхронизации несущих частот. При сближении значений управляющих параметров генераторов в колебаниях системы сначала возникают четыре независимые частоты, затем после захвата несущих частот остаются три независимые частоты, а затем захватываются частоты модуляции и, соответственно, числа вращения. Расчеты полного спектра характеристических показателей Ляпунова показали, что при отсутствии синхронизации система демонстрирует колебания с четырьмя независимыми частотами, затем при частичной синхронизации наблюдаются три независимые частоты, при полной синхронизации - две независимые частоты [2].
Однако остается неясным, от чего зависит взаимное расположение областей захвата частот генератора. Для ответа на этот вопрос будем рассматривать наиболее общий случай, при котором частоты внешнего воздействия и собственные частоты колебаний генератора не равны, а также не совпадают и их числа вращения. Будем менять несущую частоту воздействия, соответственно, будет меняться число вращения. По спектрам мощности колебаний будем исследовать области захвата частот генератора, при этом оставляя без изменений управляющие параметра системы.
Уравнения генератора квазипериодических колебаний под воздействием внешнего возмущения имеют вид
x = m ) x + y — dx3 + kF(t),
Mz)
У = —X' (8)
z = ф,
ф = —Уф + Yf2(x) — gz,
где F(t) = Ao(1 + Kmod sin Qt) sin wt - квазипериодический сигнал от внешнего источника, k - параметр связи. На рис. 3 представлена схема экспериментальной
Рис. 3. Общая схема экспериментальной установки
установки, на которой исследовалось явление синхронизации квазипериодических колебаний. Сигнал от внешнего источника представляет собой квазипериодические двухчастотные колебания. Несущая и огибающая сигнала являются синусоидальными, что соответствует классическому примеру двухчастотных колебаний.
Подадим на генератор внешний амплитудно-модулированный сигнал с частотой огибающей fie = 3.06 кГц и коэффициентом модуляции Kmod = 0.5. При этом значения параметров генератора соответствуют коэффициенту модуляции примерно 0.3, собственная частота несущей f0g = 10.69 кГц, собственная частота модуляции fig = 2.82 кГц. Будем менять частоту несущей внешнего воздействия foe в диапазоне 10 + 11 кГц.
На рис. 4 показана зависимость отношения частот модуляции, несущей генератора и внешнего воздействия от несущей частоты внешнего сигнала. Как и в случае синхронизации двух генераторов, в данном случае при изменении частоты несущей внешнего сигнала сначала захватывается частота несущей генератора, а затем - его частота модуляции. На рис. 5, а-в показаны соответствующие спектры мощности колебаний генератора для трех режимов. Как говорилось выше: при несинхронном режиме наблюдаются четыре независимые частоты; при частично синхронном режиме наблюдаются три независимые частоты (несущая и две частоты модуляции); в случае полной синхронизации - две независимые частоты (несущая и модуляции).
Ос
1.10 1.05 1.0 ^ -.95 0.90
»Л ▲ ь
А.
10.35 10.50 10.65 10.80 10.95
/ое. кГЦ
Рис. 4. Зависимость отношения частоты модуляции 2 и частоты несущей 1 генератора от частоты несущей сигнала внешнего возмущения
Рис. 5. Спектры мощности колебаний на выходе генератора: а - несинхронный режим, /0е = 10.28 кГц; б - частично синхронный режим, /0е = 10.5 кГц; в - синхронный режим, /0е = 10.69 кГц
На рис. 6, а, б представлены области синхронизации соответствующих частот генератора и внешнего воздействия в зависимости от параметра связи, снятые при различных значениях глубины модуляции внешнего воздействия. Как видно из рисунка, взаимное расположение областей захвата несущей частоты и частоты модуляции зависит не только от параметра связи, но и от глубины модуляции. На рис. 7
Рис. 6. Области синхронизации соответствующих сигналов генератора и воздействия при различных значениях глубины модуляции сигнала воздействия: а - Кто^ ~ 0.7, б - Ктоа ~ 0.9. Линиями 1С ограничена область синхронизации частот несущей, линиями 1т - частот модуляции
Рис. 7. Зависимость расположения областей синхронизации от коэффициента модуляции Кто^ при фиксированном параметре связи к; линиями 1С ограничена область синхронизации частот несущей, линиями 1т - частот модуляции
10.6 10.7 10.8 /0е, кГц
представлены области захвата несущих частот и частот модуляции в зависимости от глубины модуляции. Увеличение глубины модуляции не влияет на ширину области захвата несущей частоты, тогда как ширина области захвата частоты модуляции растет.
Следует также отметить тот факт, что в ходе эксперимента не удалось выявить режимы, при которых ширина области захвата частоты модуляции была бы шире области захвата несущей частоты. С дальнейшим ростом коэффициента модуляции области захвата несущей частоты и частоты модуляции совпадают.
Заключение
В результате проведенных исследований предложена и экспериментально апробирована радиотехническая схема генератора двухчастотных квазипериодических колебаний, соответствующая динамической системе, введенной в [1]. Экспериментально исследовано явление внешней синхронизации квазипериодических колебаний модулированным сигналом с различной глубиной модуляции. Подтверждено явление захвата числа вращения на двумерном торе, впервые установленное численно в [2].
В данной работе изучалось явление синхронизации двухчастотных автоколебаний с иррациональным соотношением частот. Синхронизация резонансных структур, возникающих на торе в случае рационального соотношения его частот, является предметом дальнейших исследований.
Работа финансировалась из средств грантов по программам «Развитие потенциала высшей школы» и BRHE.
Библиографический список
1. Анищенко В.С., Николаев С.МГенератор квазипериодических колебаний. Бифуркация удвоения двумерного тора. // Письма в ЖТФ. 2005. Т. 31, вып. 19.
С. 88.
2. Anishchenko V., Nikolaev S. and Kurths /.Winding number locking on two-dimensional torus. Synchronization of quasiperiodic motions // Phys. Rev. E. 2006. Vol. 73.
P. 056202.
3. Анищенко В.С. Сложные колебания в простых системах. М.: Наука, 1990.
4. Shilnikov L. Mathematical problems of nonlinear dynamics: A tutorial // Int. J. of Bifurcation and Chaos. 1997. Vol. 7(9). P. 1953.
Саратовский государственный Поступила в редакцию 15.10.2007
университет
EXPERIMENTAL RESEARCH OF SYNCHRONIZATION OF TWO-FREQUENCY QUASIPERIODIC MOTIONS
V.S. Anishchenko, S.M. Nikolaev
We present the electronic scheme of autonomous generator of two-frequency quasi-periodic motions and the experimental research results of effect of quasi-periodic motions synchronization under the external two-frequency force.
Анищенко Вадим Семенович - родился в 1943 году. Окончил физический факультет СГУ (1966). Защитил диссертацию на звание кандидата физико-математических наук (1970) и доктора физико-математических наук (1986). С 1988 года - заведующий кафедрой радиофизики и нелинейной динамики СГУ. С 1979 и по настоящее время работает в области исследования нелинейной динамики и стохастических процессов в нелинейных системах. Является автором более 300 научных работ, среди которых 7 монографий на русском и английском языках и 4 учебника. Неоднократно читал лекции в ведущих вузах Германии в качестве приглашенного профессора. Член-корреспондент РАЕН, заслуженный деятель науки РФ (1995), Соросовский профессор, лауреат премии Фонда Александра Гумбольдта (1999). Е-таП^а&т^Ьаов.вви.ги^.ги
Николаев Сергей Михайлович - аспирант кафедры радиофизики и нелинейной динамики СГУ. Область научных интересов - нелинейная динамика, хаотические и квазипериодические колебания. E-mail: [email protected]