Схемы электронных устройств с гиперболическим хаосом на основе связанных осцилляторов ван дер Поля

УДК 53.043

Д.С. Аржанухина

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

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

Аттрактор, гиперболический хаос, отображение «кот Арнольда», отображение Фибоначчи, осциллятор ван дер Поля

D.S. Arzhanukhina DIAGRAM OF THE ELECTRONIC HYPERBOLIC CHAOS ON THE BASIS OF RELATED VAN DER POL OSCILLATORS

Schemes of an electronic device based on coupled van der Pol oscillators. Oscillator phase behavior during the modulation period for the first system corresponds to the mapping Anosov exhibit chaotic dynamics, and for the second - display Fibonacci modification leading to the emergence of DA-hyperbolic attractor. Due to the property of structural stability, which allows to vary the system, keeping the type of dynamics, such schemes are of interest from a practical point of view. They can be used for transmission and coding, generating random numbers, cryptographic schemes.

Attractor hyperbolic chaos mapping «Arnold cat» display Fibonacci van der Pol oscillator

Введение

Одним из наиболее изученных классов динамических систем с хаотическим поведением являются системы с гиперболическими хаотическими аттракторами [1-7]. Такие системы представляют практический интерес. Благодаря свойству структурной устойчивости (нечувствительности динамики к малому изменению параметров уравнений) и возможности строгого математического описания физическая реализация систем с гиперболическим хаосом является перспективной с точки зрения применение в радиотехнических приложениях. В частности, генераторы хаоса на этой основе могут использоваться в информационно-коммуникационных целях.

В качестве особого класса систем со структурно устойчивым хаосом выступают системы Аносова [8,9]. В этих системах фазовое пространство целиком представляет собой гиперболическое инвариантное множество, и все траектории являются седловыми. Простейшим примером системы Аносова может служить отображение Фибоначчи:

jn+1 +jn _1

Двукратную итерацию (1) называют отображением «кот Арнольда» [10-12]:

поскольку В.И. Арнольд использовал графическую иллюстрацию действия этого отображения на область в виде головы кота.

Данное отображение является консервативной системой и, следовательно, аттракторов для него не существует. Для того чтобы система стала диссипативной и возник аттрактор, отображение модифицируют согласно [8]. Для этого к правой части уравнения добавляется член такого вида, что неподвижная седловая точка исходной системы становится отталкивающей, вдоль линии, отвечавшей ранее устойчивому многообразию, и в ее окрестности возникают две седловые точки, тогда как вдоль неустойчивого направления движение не меняется (рис.1) [5]. В результате получается гиперболический аттрактор, называемый DA аттрактором («Derived from Anosov»).

Рис. 1. Иллюстрация модификации отображения в окрестности неподвижной точки

Примерами систем с гиперболическими аттракторами долгое время были лишь искусственные математические конструкции, в то время как в реальных физических системах присутствуют аттракторы не являющиеся гиперболическими. Однако, в последнее время, появились работы, где приведены примеры физически реализуемых систем с гиперболическими аттракторами [5-7, 13-17]. В

21

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

В первом разделе статьи рассматривается система на основе трех осцилляторов ван дер Поля, которые благодаря модуляции бифуркационного параметра возбуждаются попеременно, а также рассматриваются два случая возможной динамики фаз для этой системы. Во втором разделе приводится схема соответствующего электронного устройства, выполненная с использованием программного пакета МиШ8т. Как будет показано, динамика фаз такой системы соответствует отображению Фибоначчи. Третий раздел посвящен исследованию видоизмененной схемы электронного устройства, динамика фаз в которой отвечает отображению Фибоначчи с БЛ модификацией.

1. Система трех осцилляторов с гиперболическим хаосом

Рассмотрим систему трех связанных осцилляторов ван дер Поля следующего вида [18,19].

dx

+ W x = [B + Acos(2pt / T) - x2]---+є

dt

d (yz) dt

d x

~dF

^ + w2y = [B + Acos2p(t /T -1/3) - y2]dy + є

dz

cos w0t + a— 0 dt

dt d2 z dt2

dt dz

+ Wz = [B + A cos2p(t / T - 2/3) - z2]---------------+є

0 dt

Г d ( xz ) dt d (xy) dt

dx

cos wt + a—

dt dy

cos w0t + a— 0 dt

Здесь x, y и z - динамические переменные осцилляторов ван дер Поля, собственная частота которых равна частоте вспомогательного сигнала ш0. Величина A определяет глубину модуляции параметра, ответственного за возбуждение осцилляторов, а B - то, какую часть периода осцилляторы

проводят выше и ниже порога возбуждения, e - параметр связи, а - дополнительный параметр связи.

Как будет ясно из дальнейшего, случай а=0 отвечает системе с динамикой типа Аносова на аттракторе, а если а достаточно велико, то в системе реализуется гиперболический DA аттрактор.

Поясним механизм функционирования системы, обратившись сначала к случаю а = 0. Предположим, что в некоторый момент осцилляторы x и y имеют фазы фх фу т.е. x ~ sin(ffi01 + px), y ~sin(«01 + py), тогда при возбуждении осциллятор z получит фазу колебаний, определяемую с точностью до константы соотношением

j'z = Фx + Фy

где штрихом обозначена величина фазы, относящейся к новой стадии активности осциллятора z. Аналогичным образом фазы изменяются и на следующих этапах передачи возбуждения, так что за полный период модуляции выполняется соотношение

p'z = Px +py ,

j'x =Фy +p'z =Фx + 2Фy ,

p'y =p'z +p'x = 2Фx + 3py .

Отсюда видно, что изменение фаз за один шаг передачи возбуждения соответствует отображению Фибоначчи, а за два - соответственно, отображению «кот Арнольда». За полный период модуляции преобразование вектора, составленного из величин ф^; и фу, определяется отображением

Г Ф/ Л

v фУ у

Г1

2

3

Л

v ф y у

Если теперь обратиться к случаю аф0, то можно показать, что осциллятор х получит фазу колебаний, определяемую соотношением [19]

Фx = Фz + Ф y + f (Ф y ) (mod 2p) =

где г(ф) = -агсід 2а51Пф . Согласно результатам [19], для аттрактора этого отображения наблюда-1 + 2а соф

ется поперечная канторова структура, характерная для БЛ аттрактора.

2. Система с динамикой, соответствующей отображению «кот Арнольда»

На рис. 2 показана схема электронного устройства, сконструированная с помощью программного пакета «Микшт 10.0», которая соответствует уравнениям (3) для случая а = 0. Схема составле-

на из трех осцилляторов ван дер Поля на основе колебательных контуров, каждый из которых образован катушкой индуктивности Ь=1,7 мГН и емкостью С=4нФ. Таким образом, собственная частота каждого осциллятора составляет примерно 60 кГц. В каждый осциллятор введено отрицательное сопротивление (при помощи операционных усилителей и1, и2, и3) и нелинейная проводимость (при помощи трех пар полупроводниковых диодов с противоположными направлениями пропускания тока Б1 и Б2, Б3 и Б4, Б5 и Б6). Для модуляции параметра, ответственного за возникновение автоколебаний, сигнал на частоте 4 кГц от источника трехфазного напряжения V1 подается на затворы полевых транзисторов Р1, Р2 и Р3, в результате чего их проводимость в цепи сток-исток меняется с соответствующим периодом. Передача возбуждения осциллятору г происходит от осцилляторов х и у, сигналы которых перемножаются умножителем Л7, затем дифференцируются схемой на операционном усилителе и7 и умножаются на вспомогательный сигнал постоянной амплитуды и частоты 60 кГц от источника переменного напряжения V3. Аналогичным образом возбуждение передается осциллятору у от осцилляторов х и г, и осциллятору х от осцилляторов у и г.

Рис. 2. Принципиальная схема устройства с динамикой фаз, описываемой отображением «кот Арнольда»

На рис. 3 показаны графики временной зависимости напряжения на конденсаторах и резисторах С1 и R1 (а), С2 и R2 (б), С3 и R3 (в), полученные при моделировании в среде Multisim с использованием двухлучевого осциллографа. Как можно видеть, осцилляторы возбуждаются поочередно, причем каждый из них генерирует последовательность цугов колебаний, следующих друг за другом с периодом модуляции Т = 0,25 мс.

Чтобы продемонстрировать соответствие динамики фаз системы отображению (4), данные моделирования в среде Multisim, записываются в файл с помощью приложения Grapher. Для каждого осциллятора создается свой файл с данными. Затем файлы обрабатываются внешней программой. Сначала выбираются значения напряжений на резисторе и конденсаторе UR и UC , для первого осциллятора соответствующие моменту времени t = (n — 3)T , для второго осциллятора -моменту времени t = (n — 3)T и для третьего осциллятора - моменту времени t = nT . Затем полученные временные ряды для напряжений нормируются так, чтобы сумма квадратов элементов была для них одинаковой, после чего для каждой пары величин Unc, UR, относящихся к одному и тому же моменту времени,

фазы определяются соотношением pn = arg(UC + iUR), и приводятся к интервалу от 0 до 2p.

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

Рис. 3. Реализации напряжения на конденсаторах

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

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

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

24

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

Другой вариант демонстрации соответствия динамики фаз в схеме отображению «кот Арнольда» состоит в том, чтобы представить полученные данные на плоскости с изображением головы кота. Как было показано выше, в системе (3) отображение для фаз за треть периода модуляции соответствует отображению Фибоначчи, за две трети - отображению «кот Арнольда», а за полный период модуляции - отображению (6). Графическое построение производится следующим образом. Для пары

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

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

Рис. 6. Иллюстрация и действия отображения, описывающего динамику фаз, полученного на основе результатов моделирования в программной среде Multisim

На рис. 7 показаны спектры колебаний трех осцилляторов, полученные при моделировании в среде МиШзт с использованием анализатора спектра. Для каждого осциллятора спектр сплошной, что подтверждает хаотическую природу системы, и сосредоточен вблизи собственной частоты 60 кГц.

л; дБ о

-70

О /5кГц 100

Рис. 7. Спектры колебаний напряжения на конденсаторах С1

3. Система с динамикой, соответствующей модифицированному отображению Фибоначчи

На рис. 8 показана схема электронного устройства с динамикой, соответствующей уравнениям (3) для случая а Ф 0 (рис. 8). Основное отличие этой схемы от схемы в разделе 2 заключается в том, что в схему вводятся дополнительные элементы, обеспечивающие поведение фаз, соответствующее динамике на гиперболическом аттракторе.

Передача возбуждения осциллятору х происходит от осцилляторов у и г, сигналы которых перемножаются умножителем А1, затем дифференцируются схемой на операционном усилителе и4, конденсаторе С4 и резисторе Я13 и умножаются на вспомогательный сигнал постоянной амплитуды и частоты 60 кГц от источника переменного напряжения У3. Затем полученный сигнал суммируется с сигналом от осциллятора г, подвергнутого преобразованию посредством дифференциатора на конденсаторе С8, резисторе Я17 и операционном усилителе и8. Аналогичным образом возбуждение передается осциллятору у от осцилляторов х и 2, и осциллятору 2 от осцилляторов х и у.

Заметим, что сигнал от источника трехфазного напряжения, подаваемый на затворы полевых транзисторов, в этой схеме имеет частоту 2 кГц.

На рис. 9 показаны графики реализаций напряжений на конденсаторах и резисторах С1 и Я1 (а), С2 и Я2 (б), С3 и Я3 (в), полученные при моделировании в среде Миііівіш. Как можно видеть, осцилляторы возбуждаются поочередно, и фаза высокочастотного заполнения относительно огибающей хаотически меняется от одного периода активности к другому.

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

На рис. 10а показан график зависимости значений фаз (рп от фп-\, где на качественном уровне можно

наблюдать структуру, соответствующую гиперболическому ЭА аттрактору. Для сравнения на рис. 10б приведен портрет аттрактора отображения (7).

Рис. 8. Принципиальная схема устройства с динамикой фаз, описываемой модифицированным отображением Фибоначчи

Рис. 9. Реализации напряжения на конденсаторах

Рис. 10. График для отображения фазы, полученный на основе результатов моделирования в программной среде МШ^пл

11с,в

-0.15

II і*,в

0.15

-0.15

II і*,в

0.15

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

Рис. 12. Спектры колебаний напряжения на конденсаторах С1

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

На рис. 12 приводятся спектры сигналов, генерируемых тремя осцилляторами, построенные с помощью анализатора спектра в логарифмическом масштабе. Спектры сосредоточены в окрестности собственной частоты каждого осциллятора 60 кГц. В низкочастотной области спектра присутствуют дискретные компоненты на частоте 2 кГ ц, что объясняется периодической модуляцией параметров с этой частотой.

Заключение

В работе представлены схемы электронных устройств, сконструированных с помощью программного пакета «Multisim 10.0». В первом разделе статьи приводятся модельные уравнения, описывающие систему связанных осцилляторов с попеременной передачей возбуждения, а также рассматриваются два случая возможной динамики фаз для этой системы. Во втором разделе рассматривается схема электронного устройства, для которой поведение фаз описывается отображением «кот Арнольда». Для схемы, представленной в третьем разделе, динамика фаз соответствует поведению на гиперболическом DA аттракторе. Благодаря свойству структурной устойчивости, вытекающему из математической теории гиперболических аттракторов, можно варьировать систему, сохраняя тип динамики, в силу чего генераторы хаоса данного типа представляют особый интерес с практической точки зрения. Они могут использоваться, в частности, для передачи и кодирования информации, генерации случайных чисел, криптографических схем.

Работа выполнена при поддержке гранта РФФИ№12-02-00342.

ЛИТЕРАТУРА

1. Afraimovich V. Lectures on chaotic dynamical systems / V. Afraimovich and S.-B. Hsu // AMS/IP Studies in Advanced Mathematics, Vol. 28. American Mathematical Society, Providence RI, International Press, Somerville, MA, 2003.

2. Гукенхеймер Дж. Нелинейные колебания, динамические системы и бифуркации векторных полей / Дж. Гукенхеймер, П. Холмс. М.; Ижевск: Ин-т компьютерных исследований, 2002. 559 с.

3. Devaney R.L. An Introduction to Chaotic Dynamical Systems / R.L. Devaney. N.Y: Addison -Wesley, 1989.

4. Shilnikov L. Mathematical problems of nonlinear dynamics: a tutorial / L. Shilnikov // Int. J. of Bif. & Chaos. 1997. Vol. 7, № 9. P. 1353.

5. Кузнецов С.П. Гиперболические странные аттракторы систем, допускающих физическую реализацию / С.П. Кузнецов // Изв. вузов. ПНД. 2009. Т. 17. № 4.

6. Кузнецов С.П. Пример неавтономной системы с непрерывным временем, имеющей аттрактор типа Плыкина в отображении Пуанкаре / С.П. Кузнецов // Нелинейная динамика. 2009. Т. 5. № 3. С. 403-424.

7. Кузнецов С.П. Динамический хаос и однородно гиперболические аттракторы: от математики к физике / С.П. Кузнецов // Успехи физических наук. 2011. Т. 181. № 2. С. 121-149.

8. Каток А.Б. Введение в современную теорию динамических систем: пер. с англ. / А.Б. Каток, Б. Хасселблат М.: Факториал, 1999. 768 с.

9. Coudene Y. Pictures of Hyperbolic Dynamical Systems / Y. Coudene // Notices of the American Mathematical Society. 2006. 53 (1). P. 8-13.

10. Кузнецов С.П. Динамический хаос: курс лекций / С.П. Кузнецов. М.: Физматлит, 2001. 296 с.

11.Берже П. Порядок в хаосе. О детерминистском подходе к турбулентности / П. Берже,

И. Помо, К. Видаль. М.: Мир, 1991, 368 с.

12.Шустер Г. Детерминированный хаос / Г. Шустер. М.: Мир, 1988. 240 с.

13. Кузнецов С.П. Хаотическая динамика в физической системе со странным аттрактором ти-

па Смейла - Вильямса / С.П. Кузнецов, Е. П. Селезнев // ЖЭТФ. 2006. Т. 129. Вып. 2. C. 400-412.

14.Belykh V. Hyperbolic Plykin attractor can exist in neuron models / V. Belykh, I. Belykh, E. Mo-sekilde // International Journal of Bifurcation and Chaos. 15. 2005. № 11. 3567-3578.

15.Kuznetsov S.P. Plykin type attractor in electronic device simulated in Multisim / S.P. Kuznetsov // CHAOS. 21. 2011. 043105.

16.Isaeva O.B. Arnold's cat map dynamics in a system of coupled nonautonomous van der Pol oscillators / O.B. Isaeva, A.Yu. Jalnine, S.P. Kuznetsov // Phys. Rev. E 74. 2006. 046207.

17.Kuznetsov S.P. Autonomous coupled oscillators with hyperbolic strange attractors /

S.P. Kuznetsov, A. Pikovsky // Physica. D232. 2007. 87-102.

18.Аржанухина Д.С. Система трех неавтономных осцилляторов с гиперболическим хаосом.

Ч. I. Модель с динамикой на аттракторе, описываемой отображением на торе «кот Арнольда» / Д.С. Аржанухина, С.П. Кузнецов // Изв. вузов. ПНД. 2012. Т. 20. № 6. С. 56-66.

19.Аржанухина Д.С. Система трех неавтономных осцилляторов с гиперболическим хаосом.

Ч. II. Модель с DA-аттрактором / Д.С. Аржанухина, С.П. Кузнецов // Изв. вузов. ПНД. 2013. Т. 21.

№ 2.С. 163-172.

АржанухинаДарья Сергеевна - Darya S. Arzhanukhina -

аспирант кафедры «Динамические системы» student

Саратовского государственного университета Department of Dynamic systems

имени Н.Г. Чернышевского N.G. Chernyshevsky Saratov State University

Статья поступила в редакцию 17.08.13, принята к опубликованию 15.09.13