Прикладные задачи
^^^^^^^^^^»нелинейной теории колебаний и вслн
Изв. вузов «ПНД», т. 21, № 5, 2013 УДК 517.9:537.86:621.373
АВТОНОМНАЯ СИСТЕМА - ГЕНЕРАТОР ГИПЕРБОЛИЧЕСКОГО ХАОСА: СХЕМОТЕХНИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ЭКСПЕРИМЕНТ
С. П. Кузнецов, В. И. Пономаренко, Е. П. Селезнев
Рассмотрена схема электронного устройства, представляющего собой автономную динамическую систему с гиперболическим аттрактором типа Смейла-Вильямса в отображении Пуанкаре, и выполнено моделирование хаотической динамики в программной среде Ми^т. Генератор гиперболического хаоса реализован в виде лабораторного макета, проведено его экспериментальное исследование и продемонстрировано соответствие наблюдаемой динамики результатам численного и схемотехнического моделирования. Привлекательной для возможных приложений является грубость, или структурная устойчивость, системы обусловленная гиперболической природой аттрактора, то есть нечувствительность характеристик генерируемого хаоса к помехам, шумам, погрешностям изготовления.
Ключевые слова: Хаос, аттрактор, соленоид Смейла-Вильямса, автономная система, эксперимент, операционный усилитель, МиШ8т.
Введение
Равномерно гиперболические хаотические аттракторы, обладающие свойством грубости, или структурной устойчивости, такие как соленоид Смейла-Вильямса или аттрактор Плыкина, были введены математиками нескольких десятилетий назад [1-4]. В то время ожидалось, что они могут служить для описания хаоса и турбулентности, однако впоследствии выяснилось, что встречающиеся в приложениях хаотические аттракторы обычно не попадают в этот класс. Физически реализуемые системы с гиперболическим хаосом были обнаружены (или, скорее, сконструированы) только совсем недавно [5-8]. В связи с этим возникли предпосылки для создания на этой основе реально функционирующих устройств, демонстрирующих структурно устойчивый хаос, например, в электронике с возможным их применением в схемах скрытой коммуникации [9, 10], шумовой локации [11], в криптографии [12, 13], для генерации случайных чисел [14]. Привлекательность структурно устойчивых систем [15-17], с практической точки зрения, обусловлена нечувствительностью характеристик генерируемого хаоса к погрешностям изготовления, шумам, помехам и т.п.
При разработке нетривиальных вопросов, связанных с приложениями, важным моментом будет существование надежной теоретической базы, обеспеченной наличием хорошо развитой математической теории. Например, возможность полного символического описания траекторий на однородно гиперболическом аттракторе позволяет трактовать динамический процесс в контексте идей и представлений теории информации и кодирования, что, по-видимому, важно с точки зрения информационных и коммуникационных приложений.
Первое экспериментальное исследование системы с аттрактором Смейла-Виль-ямса было представлено в работе [18]; также реализованы и изучены экспериментально варианты систем, использующих запаздывающую обратную связь [19, 20]. Несколько схем, отвечающих устройствам с аттракторами типа Смейла-Вильямса и Плыкина, рассмотрено в работах [21, 22] с привлечением схемотехнического моделирования в программной среде МиШзт.
Приходится отметить, однако, что все реальные схемы электронных устройств, где пока удалось осуществить гиперболический хаос, отвечают неавтономным системам, требующим для своего функционирования заданного определенным образом периодического внешнего воздействия. До последнего времени не был решен вопрос о реализации равномерно гиперболического аттрактора в автономной системе, то есть о создании полноценного генератора гиперболического хаоса с демонстрацией его функционирования в лабораторном эксперименте. В этом состоит цель исследования, результаты которого представлены в настоящей статье.
1. Автономная система дифференциальных уравнений с аттрактором Смейла-Вильямса
В работе [6] рассмотрено несколько искусственных примеров автономных систем с гиперболическими аттракторами, одна из которых имеет минимальную размерность, допускающую существование аттрактора Смейла-Вильямса. Здесь она будет принята за основу с некоторой модификацией, имеющей целью упрощение схемотехнической реализации.
Следуя [6], обратимся сначала к двумерной системе типа «хищник-жертва», мгновенное состояние которой задается парой неотрицательных переменных Т\, Г2,
Г1 = 2[1 - Г2 + 1 Г1 - 50(Г1 - 1)2]Г1, Г2 = 2(Г1 - 1)Г2. (1)
В установившемся режиме она демонстрирует периодические автоколебания, причем на плоскости переменных п, г2 изображающая точка посещает окрестность начала координат раз за разом. После каждого такого прохождения происходит сначала возбуждение первой подсистемы (переменная п, «жертва»), потом возбуждение второй (г2, «хищник»), далее затухание первой подсистемы, более медленное затухание второй, и затем цикл повторяется.
Пусть теперь г1 и г2 представляют собой квадраты модуля двух комплексных величин й1 и а2. Запишем уравнения для этих переменных и дополним их членами, описывающими связь по приводит к уравнениям
2
описывающими связь подсистем, пропорциональными, соответственно, а2 и аь что
а1 = [1 - |а2|2 + 1М2 - 50(М2 - 1)2]а1 - 1 £1а^, а2 = (|а112 - 1)а2 - ^аь
Здесь е\ и £2 - параметры связи, которые в дальнейшем изложении фиксированы и приняты равными е\ = 0.04, £2 = 0.4.
В системе с комплексными переменными (2) активизация второй подсистемы происходит в присутствии воздействия со стороны первой подсистемы, благодаря чему аргумент а2 наследует величину аргумента, отвечавшую переменной а\. Затем на стадии затухания вторая подсистема, в свою очередь, обеспечивает затравочное воздействие для первой при очередном прохождении орбиты вблизи начала координат. Поскольку соответствующий член содержит квадрат комплексной переменной а2, передача возбуждения сопровождается удвоением аргумента комплексного числа. Далее процесс повторяется, причем на каждом новом цикле угловая переменная, отвечающая аргументу комплексных величин, умножается на фактор 2. Это соответствует растягивающему отображению окружности - отображению Бернулли Фга+1 = 2Фп, которое характеризуется хаотической динамикой и имеет положительный показатель Ляпунова Л = 1п2 ~ 0.693.
Полагая а1 = х + гп и а2 = у + гь, можно отделить в (2) действительные и мнимые части и перейти к уравнениям в действительных переменных
X = [1 - (У2 + ь2) + 1(х2 + п2) - 50(х2 + п2 - 1)2]х - 2£1(у2 - ь2),
п = [1 - (у2 + V2) + 2(х2 + п2) - 50(х2 + п2 - 1)2]п - £1уь, (3)
У = (х2 + п2 - 1)у - £2х, V = (х2 + п2 - 1)ь - £2п.
Для описания динамики в терминах отображения Пуанкаре введем сечение в четырехмерном фазовом пространстве подходящей гиперповерхностью и определим соответствующее трехмерное отображение на этой гиперповерхности. Его аттрактором будет служить объект в виде соленоида Смейла-Вильямса, поскольку имеется угловая переменная, претерпевающая удвоение на каждой очередной итерации, и, как следует из анализа показателей Ляпунова (см. ниже), это сопровождается сжатием фазового объема по остальным направлениям в пространстве состояний.
На рис. 1 приведены графики зависимости от времени динамических переменных, полученные при численном решении уравнений (3). Можно видеть, что стадии активности одной и другой подсистемы перемежаются со стадиями подавления. Согласно вычислениям, средний временной период повторения стадий для данного режима {Т) = 5.68.
На рис. 2 показаны портреты аттрактора - один на плоскости переменных, отвечающих первой подсистеме (х,п), а другой на плоскости амплитудных переменных двух подсистем (г1, г2) = (х2 + п2, у2 + V2).
На рис. 2, б возбуждение первой подсистемы соответствует уходу траектории от начала координат вправо, возбуждение второй - уходу вверх с последующим подавлением сначала первой, а затем второй подсистемы, и новым возвращением в окрестность начала координат.
Для задания секущей Пуанкаре используем соотношение г2 = у2 + V2 = с0, со = 3, причем учитываем проходы траекторий только в направлении увеличения амплитуды г2.
На рис. 3 показан портрет аттрактора в сечении Пуанкаре в проекции на плоскость переменных второй подсистемы и приведена диаграмма, иллюстрирующая преобразование циклической переменной - аргумента комплексной амплитуды
Рис. 1. Зависимости от времени динамических переменных, полученные при численном решении уравнений (3) для £1 = 0.04, £2 = 0.4
Рис. 2. Портреты аттрактора системы (3) при £1 = 0.04, £2 = 0.4: на плоскости переменных, отвечающих первой подсистеме (а) и на плоскости амплитудных переменных первой и второй подсистемы (б)
Рис. 3. Портрет аттрактора в сечении Пуанкаре в проекции на плоскость переменных первой подсистемы (а) и диаграмма преобразования циклической переменной (аргумента комплексной амплитуды) второй подсистемы при последовательных проходах секущей Пуанкаре (б) для системы (3) при £1 = 0.04,
£2 = 0.4
второй подсистемы при последовательных проходах секущей Пуанкаре. Как можно видеть, оно с замечательной точностью соответствует отображению Бернулли.
Показатели Ляпунова для данного аттрактора, рассчитанные по методу Бенет-тина [23],
А.1 = 0.1219, А, 2 = 0.0000, Аз = -0.6867, \А = -1.0761. (4)
Учитывая, что средний период прохода секущей Пуанкаре {Т) = 5.68, для старшего показателя отображения Пуанкаре получаем Л1 = А,1{Т) = 0.692, что прекрасно согласуется с ожидаемым значением 1п2 = 0.6931.... Второй показатель нулевой, с точностью до погрешности вычислений, и должен интерпретироваться, как связанный с вектором возмущения типа сдвига вдоль опорной траектории. Остальные показатели отрицательные. Размерность по Каплану-Йорке [23] для аттрактора отображения Пуанкаре Оку = 1 + А1/|А3| « 1.18, а для аттрактора потоковой системы на единицу больше.
2. Аналоговое устройство с аттрактором Смейла-Вильямса
При конструировании описанной ниже схемы принята идеология по возможности точного воспроизведения уравнений с использованием элементов, применяемых в технике аналогового моделирования, таких как интеграторы, умножители, сумматоры [24,25], что представляется разумным в качестве первого шага на пути построения реальных генераторов гиперболического хаоса. В дальнейшем, по-видимому, можно будет реализовать более простые схемы на традиционной для радиотехники элементной базе (транзисторы, диоды, конденсаторы, резисторы).
Схема, представленная на рис. 4, содержит операционные усилители Ш-Ш и умножители А1-А10. Все умножители характеризуются коэффициентом передачи К, равным по абсолютной величине 1/10 (то есть при входных напряжениях Ц и и2 на выходе получается напряжение КИ1И2), причем для умножителей А6 и А7, обозначенных на схеме серым цветом, этот коэффициент отрицательный.
С каждой из четырех динамических переменных х, и, у, V ассоциируется интегратор на базе операционного усилителя (соответственно, И1, И2, И3, И4), емкости (С1, С2, С3, С4) и резистора (Я1, Я2, Я3, Я4). Собственно величины х, и, у, V отвечают напряжениям на конденсаторах С1, С2, С3 и С4, соответственно.
Динамика во времени соответствует системе уравнений (3) при измерении времени в единицах т. Это константа с размерностью времени т = КО, определенная через емкость О = О1 = О2 = О3 = О4 и характеристическое сопротивление К = 100^,2,3,4 = 10 кОм. При указанных на схеме номиналах т = 0.22 мс.
Умножители А1, А2, операционный усилитель И5 и резисторы Я13, Я14, Я7, Я8 образуют блок, формирующий коэффициент = х2 + и2 - 1. Соответствующий сигнал, будучи умножен на у и V посредством элементов А8 и А9, подается на интеграторы, отвечающие за третье и четвертое уравнения системы (3). Также сигнал ^2 и его квадрат, получаемый посредством умножителя А5, используются, вместе с выходным сигналом блока на элементах А3, А4, И6, для получения коэффициента = 1 - (у2 + V2) + 1/2 (х2 + и2) - 1/50 (х2 + и2 - 1)2 на выходе операционного усилителя и6. Коэффициенты, присутствующие в выражении для Ц1, определяются соотношением сопротивлений Я5, Я9, Я10, Я15 и Я17. Будучи умножен на х и у посредством элементов А6, А7, сигнал поступает на интеграторы, соответствующие первому и второму уравнениям (3). Члены, описывающие воздействие первой
Рис. 4. Схема устройства, динамика которого описывается уравнениями (3), где за единицу времени принята величина т = 0.22 мс, а коэффициенты связи составляют e1 = 0.04, e2 = 0.4. Динамические переменные x, y, z, v отвечают напряжениям на конденсаторах C1, C2, C3, C4, измеренным в вольтах
подсистемы на вторую, формируются благодаря подаче сигналов x и u через одинаковые резисторы R11 и R12. (Их сопротивление отвечает за величину параметра связи £2, будучи ему обратно пропорциональным.) Член в первом уравнении (3), определяющий воздействие со стороны второй подсистемы, формируется благодаря подаче сигнала, пропорционального разности квадратов переменных y и v с выхода операционного усилителя U7 на соответствующий интегратор через резистор R20. Аналогично член, отвечающий за связь во втором уравнении, обеспечивается выходным сигналом умножителя A10, подаваемым на интегратор через резистор R18. Выбор R20=2 R18 обеспечивает нужное соотношение членов в первом и втором уравнениях. Коэффициент связи £i находится в обратной пропорциональности с величинами R18 и R20.
3. Моделирование динамики устройства в среде Multisim
Удобным современным средством схемотехнического моделирования является программный продукт Multisim компании National Instruments [26,27]. Работая с программой NI Multisim, пользователь составляет и изображает графически на экране компьютера электронную схему в обычном виде, следуя стандартным обозначениям, принятом в радиотехнике и электронике. Радиодетали берутся из имеющейся в распоряжении библиотеки электронных компонентов. Предусмотрено использование обширного набора виртуальных приборов, позволяющих выполнять измерения или отслеживать изменение различных электрических величин, в том числе генераторы сигналов, осциллографы, анализатор спектра. Обращение с ними приближено к работе с реальными приборами в лаборатории: подключив виртуальный прибор к соответствующей точке схемы, можно получить информацию о характеристиках процессов, протекающих при ее функционировании.
На рис. 5 приведены осциллограммы напряжений х, и, у, V, полученные при моделировании динамики в среде МиШБт с использованием виртуального многоканального осциллографа, четыре входа которого подключены к соответствующим узлам схемы рис. 4. Как показывает сравнение с рис. 1, вид зависимостей качественно такой же, как для уравнений (3). Средний период повторения стадий возбуждения и подавления подсистем, как можно оценить из рис. 5, составляет примерно 1.25 мс, то есть 5.68 единиц, что согласуется с приведенным в разделе 1 численным результатом.
На рис. 6 приводятся портреты аттрактора, скопированные с экрана виртуального осциллографа. Для получения портрета на плоскости переменных, отвечающих первой подсистеме, показанного на диаграмме (а), вход осциллографа, управляющий горизонтальной разверткой луча, подключался к выходу операционного усилителя и1, а вход, отвечающий за вертикальное отклонение, к выходу операционного усилителя и2. Портрет на диаграмме (б) представлен на плоскости амплитудных переменных подсистем (х2 + и2, у2 + V2). Для его получения один вход осциллографа подключался к выходу операционного усилителя И5, где напряжение пропорционально 1 — (х2 + и2), а второй - к выходу специально добавленного сумматора (на схеме не показан), на два входа которого подавались напряжения от умножителей А3 и А4. Портреты аттрактора находятся в прекрасном соответствии с рис. 2.
Чтобы построить портрет аттрактора в сечении Пуанкаре и диаграмму, которая иллюстрировала бы характерное для аттрактора Смейла-Вильямса преобразование
О 10 20 30 40 г, мс
Рис. 5. Зависимости от времени динамических переменных, полученные с помощью четырехканаль-ного цифрового осциллографа при моделировании схемы рис. 4 в среде Ми^т
у2 + 1)2
10.0
Рис. 6. Портреты аттрактора, полученные с помощью виртуального осциллографа при моделировании схемы рис. 4 в среде МиШ8т, на плоскости переменных первой подсистемы (х, и) (а) и на плоскости квадратичных амплитуд двух подсистем (б)
циклической переменной, приходится обратиться к записи данных моделирования в среде Multisim в файл с последующей их обработкой специально составленной программой. Для этого используем такое же подключение многоканального осциллографа, как при анализе временных зависимостей, и инструмент Grapher, обеспечивающий возможность записи в файл. Выборка делается с шагом по времени 1 мкс, на три порядка меньше характерного периода {T). Затем файл обрабатывается программой, которая, считывая шаг за шагом записанные данные, вычисляет величину y2 + v2. В момент tn, когда эта величина, изменяясь в сторону увеличения, проходит уровень с0 = 3, определяются значения xn, un, yn и vn. Представляя их в координатах (x, u), строим портрет аттрактора в сечении Пуанкаре, который показан на рис. 7, а, а с помощью набора значений угловой переменной Фп = arg(y(tn) + iv(tn)) получаем график отображения на диаграмме рис. 7, б в координатах (Фп, Фп+1). Сравнение рис. 7 и 3 показывает, что визуально они практически не отличаются.
На рис. 8 показаны спектры генерируемых системой сигналов, полученные с помощью виртуального анализатора спектра в Multisim, при соответствующих установках рабочего диапазона частот и разрешения анализа. Спектры приведены в ло-
Рис. 7. Построенные обработкой данных моделирования в среде МиШ8т для схемы рис. 3: а - портрет аттрактора в сечении Пуанкаре в проекции на плоскость переменных первой подсистемы; б - диаграмма, иллюстрирующая преобразование циклической переменной Ф = + гг) при последовательных проходах секущей Пуанкаре
Рис. 8. Спектры колебаний напряжения на конденсаторах С1 (а) и С3 (б), полученные с помощью виртуального анализатора спектра при моделировании схемы рис. 4 в среде МиШвт
гарифмическом масштабе. График на рис. 8, а соответствует сигналу х (для переменной и спектр выглядит точно так же), а график на рис. 8, б - сигналу у (для переменной V спектр выглядит так же). Непрерывный характер спектра свидетельствует о хаотической природе динамики на аттракторе.
4. Лабораторный макет и его экспериментальное исследование
Экспериментальная установка представляла собой макет, изготовленный по схеме, показанной на рис. 4. В качестве операционных усилителей использовались микросхемы типа АЭ822АК, а в качестве аналоговых умножителей - микросхемы АЭбЗЗАК. Сбор и анализ данных осуществлялся с помощью четырехканального цифрового осциллографа типа М808104А. Осциллограф позволяет записывать в цифровом виде данные с точностью 11 двоичных разрядов при максимальной скорости записи 109 выборок в секунду и осуществлять целый ряд вычислений, в том числе и преобразование Фурье.
На рис. 9 показаны временные зависимости напряжений х, и, у, V, полученные для лабораторного макета. Как видно из сравнения с рис. 1 и 5, они качественно хорошо согласуются с зависимостями, построенными по результатам численного решения уравнений (3) и для модели в среде МиШБт. Средний период повторения стадий возбуждения и подавления подсистем, как можно оценить из рисунка, составляет примерно 1.25 мс в согласии с результатами предыдущих двух разделов.
На рис. 10 представлен спектр переменной х, который можно сравнить с рис. 8, а.
На рис. 11 показаны портреты аттрактора на плоскости переменных, отвечающих первой подсистеме, и на плоскости амплитудных переменных подсистем (х2 + и2, у2 + V2). Первый из них (а) получен в эксперименте непосредственно на экране осциллографа и демонстрирует несомненное сходство с рис. 2, а и 6, а.
Портрет аттрактора на диаграмме (б) построен путем обработки полученного в эксперименте четырехкомпонентного временного ряда с шагом выборки 1 мкс при записи в файл длинной реализации (порядка 106 отсчетов) для напряжений х, и, у, V. Вид аттрактора демонстрирует несомненное сходство с рис. 2, б и 6, б, хотя надо
X
5В
и
У
V
О 10 20 30 40 мс
Рис. 9. Зависимости переменных от времени, полученные при помощи цифрового осциллографа MS08104A
Рис. 10. Спектр переменной х, полученный при помощи цифрового осциллографа MSO8104A. По вертикали - 10 dB/дел, по горизонтали - 10 кГц/дел
а
-6.0
-6.0 -3.0
12.0 х2+и2
Рис. 11. Фазовый портрет аттрактора на плоскости (х,и), полученный при помощи цифрового осциллографа MSO8104A (а) и портрет аттрактора в проекции на плоскость амплитудных переменных первой и второй подсистемы (б), построенный обработкой записанных в файл данных для экспериментального макета
отметить довольно существенное зашумление и поперечное уширение аттрактора, обусловленное, очевидно, техническими флуктуациями, возникающими при функционировании реального электронного устройства.
Для построения портрета аттрактора в сечении Пуанкаре и диаграммы, иллюстрирующей преобразование типа отображения Бернулли для циклической переменной, записанный в ходе эксперимента четырехкомпонентный временной ряд подвергается обработке специально составленной программой аналогично тому, как это делалось с данными моделирования в среде МиШзт. Считывая шаг за шагом записанные данные, вычисляем величину г2 = у2 + V2. В момент ¿п, когда эта величина, изменяясь в сторону увеличения, проходит уровень со = 3, определяем значения хп, ип, уп и vn. Учитывая технические флуктуации, которые могут вести к ложному определению момента прохождения секущей Пуанкаре, при обработке исключаются события, происшедшие в пределах 10 шагов по выборке от предыдущего прохода. Также исключаются события, которые отвечают малым амплитудам первой подсистемы п = х2 + п2 в момент прохода. (Ложные срабатывания этого последнего типа отвечают локальному нарушению монотонной зависимости Г2 от времени на
2.0
-2.0
-4.0
•г ,.
■Л
:>' л*
ч
^ / 1-
1
1
Л,- +
/
■ 1 - я
-4.0 -2.0
0
2.0
а
Рис. 12. Построенный путем обработки экспериментальной реализации портрет аттрактора в сечении Пуанкаре в проекции на плоскость переменных первой подсистемы (а) и диаграмма, иллюстрирующая преобразование циклической переменной Ф = + т) при последовательных проходах секущей
Пуанкаре, (б)
нисходящем участке траектории из-за технических флуктуаций.) Представляя в координатах (х, и) данные, отвечающие зарегистрированным моментам пересечения заданного уровня г2, с исключением ложных событий, строим портрет аттрактора в сечении Пуанкаре (рис. 12, а), а с помощью набора значений угловой переменной Фп = arg(y(tn) + гу^п)) в координатах (Фп, Фп+^ получаем график отображения на диаграмме (б). Сравнение полученных диаграмм с рис. 7 и 3 показывает хорошее качественное соответствие, несмотря на довольно сильное поперечное ушире-ние показанного на рис. 11, б аттрактора из-за технических флуктуаций в реальном эксперименте.
Приведенные результаты показывают, что механизм функционирования, ассоциирующийся с присутствием гиперболического хаоса, в реальном устройстве сохраняется, несмотря на довольно большой уровень флуктуаций. Это обстоятельство можно рассматривать как проявление структурной устойчивости наблюдаемого хаотического аттрактора.
Заключение
В статье предложена схема электронного устройства, представляющего собой первый реальный пример автономной электронной системы, демонстрирующей гиперболический хаос. В отображении Пуанкаре, описывающем изменение состояния системы на последовательных стадиях эволюции во времени, хаотическая динамика ассоциируется с аттрактором типа Смейла-Вильямса. Основанием трактовать аттрактор как соленоид Смейла-Вильямса служит то обстоятельство, что удается указать угловую циклическую переменную, которая за время возврата на секущую Пуанкаре претерпевает двукратно растягивающее отображение окружности (отображение Бернулли), тогда как по остальным направлениям в пространстве состояний имеет место сжатие фазового объема. Следует отметить, что рассмотренная схема отвечает системе минимальной размерности, делающей возможным существование аттрактора Смейла-Вильямса.
Одно из возможных применений такого рода устройств состоит в их использовании для генерации случайных чисел, что считается важной прикладной задачей [14,17]. Для систем с аттрактором типа Смейла-Вильямса бинарная случайная последовательность получается естественным образом как символическая последо-
вательность, кодирующая траекторию на аттракторе. (Ноль или единица отвечает попаданию угловой переменной на каждом периоде в одну или другую половину полного интервала.) Согласно принятой в литературе классификации, математическую модель (3) в этом контексте следовало бы отнести к «генераторам псевдослучайных чисел» (в англоязычной литературе употребляется термин Pseudo-Random Number Generator). Однако, коль скоро реализовано соответствующее физическое устройство, мы должны отнести его к классу «истинных генераторов случайных чисел» (True Random Number Generator), так как в этом случае функционирование протекает в присутствии неизбежных шумов и флуктуаций. Их усиление до макроскопического уровня в силу присущей динамике на гиперболическом хаотическом аттракторе чувствительности к возмущению фазовых траекторий приводит к тому, что в процессе временной эволюции система выбирает траекторию на аттракторе реально случайным образом.
Работа подготовлена при финансовой поддержке РФФИ (грант 12-02-00342).
Библиографический список
1. Смейл С. Дифференцируемые динамические системы // УМН. 1970. Т. 25, № 1. С. 113.
2. Аносов Д.В. и др. Динамические системы с гиперболическим поведением. Динамические системы - 9 // Итоги науки и техн. Сер. Соврем. пробл. мат. Фун-дам. направления. Т. 66. М.: ВИНИТИ, 1991. 248 с.
3. Синай Я.Г. Стохастичность динамических систем // Нелинейные волны / ред. А.В. Гапонов-Грехов. М.: Наука, 1979. С. 192.
4. Shilnikov L. Mathematical problems of nonlinear dynamics: A tutorial // Int. J. of Bifurcation and Chaos. 1997. Vol. 7, № 9. P. 1353.
5. Kuznetsov S.P. Example of a physical system with a hyperbolic attractor of the Smale-Williams type // Phys. Rev. Lett. 2005. Vol. 95. 144101.
6. Kuznetsov S.P., PikovskyA. Autonomous coupled oscillators with hyperbolic strange attractors // Physica D. 2007. Vol. 232. P. 87.
7. Кузнецов С.П.Динамический хаос и однородно гиперболические аттракторы: От математики к физике // УФН. 2011. Т. 181, № 2. С. 121.
8. Kuznetsov S.P. Hyperbolic Chaos: A Physicist's View. Higher Education Press: Beijing and Berlin, Heidelberg: Springer-Verlag, 2012. 336 p.
9. Дмитриев А.С., Панас А.И. Динамический хаос: Новые носители информации для систем связи. М.: Физматлит, 2002. 252 c.
10. Короновский А.А., Москаленко О.И., Храмов А.Е. О применении хаотической синхронизации для скрытой передачи информации // УФН. 2009. Т. 179, № 12. Р. 1281.
11. Lukin K.A. Noise radar technology // Telecommunications and Radio-Engineering. 2001. Vol. 16, № 12. P. 8.
12. Baptista M.S. Cryptography with chaos // Physics Letters A. 1998. Vol. 240. P. 50.
13. Птицын Н.В. Приложение теории детерминированного хаоса в криптографии. М.: Изд-во МГТУ им. Н. Э Баумана, 2002. 80 с.
14. Stojanovski T., Kocarev L. Chaos-based random number generators: Part I and II // IEEE Trans. Circuits and Systems. 2001. Vol. 48, № 3. P. 281; P. 382.
15. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М.: Наука, 1959. 915 с.
16. Рабинович М.И., Трубецков Д.И. Введение в теорию колебаний и волн. М.: Наука, 1984. 432 с.
17. Elhadj Z., SprottJ.C. Robust Chaos and Its Applications. Singapore: World Scientific, 2011. 472 p.
18. Кузнецов C.П., Селезнев Е.П. Хаотическая динамика в физической системе со странным аттрактором типа Смейла-Вильямса // ЖЭТФ. 2006. Т. 129, № 2. С. 400.
19. Кузнецов С.П., Пономаренко В.И. О возможности реализации странного аттрактора типа Смейла-Вильямса в радиотехническом генераторе с запаздыванием // Письма в ЖТФ. 2008. Т. 34, вып. 18. С. 1.
20. Баранов С.В., Кузнецов С.П., Пономаренко В.И. Хаос в фазовой динамике осциллятора ван дер Поля с модулированной добротностью и дополнительной запаздывающей обратной связью // Известия вузов. Прикладная нелинейная динамика. 2010. Т. 18, № 1. С. 11.
21. Kuznetsov S.P. Plykin type attractor in electronic device simulated in MULTISIM // CHAOS. 2011. Vol. 21. 043105.
22. Кузнецов C.П. Схемы электронных устройств с гиперболическим хаосом и моделирование их динамики в программной среде Multisim // Известия вузов. Прикладная нелинейная динамика. 2011. Т. 19, № 5. С. 98.
23. Кузнецов С.П. Динамический хаос. М.: Физматлит, 2001. 296 с.
24. Хоровиц П., Хилл У. Искусство схемотехники. М.: Мир, 1986. Т. 1: 510 с. Т. 2: 592 с.
25. Ремпен И.С., Егоров Е.Н., Савин А.Н., Пономаренко В.И.Операционные усилители. Учебно-методическое пособие. Саратов: Изд-во ГосУНЦ «Колледж», 2004. Ч. 1: 19 с. Ч. II: 16 c.
26. Макаренко В.В. Моделирование радиоэлектронных устройств с помощью программы NI MULTISIM. Электронные компоненты и системы (Киев) VD MAIS, 2008, № 1, 50-56; № 2, 51-57; № 3, 44-51; № 4, 44-51, № 6, 46-53; № 7, 54-59; № 8, 46-56; № 9, 65-69; № 12, 47-52.
27. Варзарев Ю.Н., Иванцов В.В., Спиридонов Б.Г. Моделирование электронных схем в системе Multisim. Таганрог: Изд-во ТТИ ЮФУ, 2008. 81с.
СФ ИРЭ им. В.А. Котельникова РАН Поступила в редакцию 27.02.2013
AUTONOMOUS SYSTEM GENERATING HYPERBOLIC CHAOS: CIRCUIT SIMULATION AND EXPERIMENT
S.P. Kuznetsov, V.I. Ponomarenko, E.P. Seleznev
We consider an electronic device, which represents an autonomous dynamical system with hyperbolic attractor of the Smale-Williams type in the Poincare map. Simulation of chaotic dynamics in the software environment Multisim has been undertaken. The generator of hyperbolic chaos is implemented as a laboratory model; its experimental
investigation is carried out, and good compliance with the observed dynamics in the numerical and circuit simulation has been demonstrated. A distinctive feature of the system is attractiveness for potential applications due to its robustness, or structural stability, i.e. insensitivity of the generated chaos characteristics in respect to noises, technical fluctuations, fabrication imperfections etc.
Keywords: Chaos, attractor, Smale-Williams solenoid, autonomous system, experiment, operational amplifier, Multisim.
Кузнецов Сергей Петрович - родился в Москве (1951). Окончил Саратовский государственный университет (1973). Сотрудник Саратовского филиала Института радиотехники и электроники им. В.А. Котельникова РАН (с 1988). Руководит лабораторией теоретической нелинейной динамики. Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в СГУ (1977) и доктора наук (1987) по специальности радиофизика. Профессор базовой кафедры динамических систем СГУ. Автор монографии «Динамический хаос и гиперболические аттракторы: От математики к физике», монографии «Hyperbolic Chaos. A Physicist's View», монографии «Strange Nonchaotic Attractors» (совместно с А. Пиковским и У. Фойдель), а также учебно-научных монографий «Динамический хаос» и «Нелинейные колебания» (в соавторстве с А.П. Кузнецовым и Н.М. Рыскиным). Опубликовал свыше 200 научных статей по нелинейной динамике, радиофизике и электронике. Под его руководством защищены десять кандидатских диссертаций. Лауреат государственной научной стипендии для ученых России (1994-1996), Соросовский доцент (1998), Соросовский профессор (2000, 2001). Член редакционной коллегии журнала «Известия вузов. Прикладная нелинейная динамика».
410019 Саратов, ул. Зеленая, 38
Саратовский филиал ИРЭ имени В.А. Котельникова РАН E-mail: spkuz@yandex.ru
Пономаренко Владимир Иванович - родился в 1960 году в Саратове. Окончил Саратовский государственный университет (1982). Работает в Саратовском филиале Института радиотехники и электроники им. В.А. Котельникова РАН (с 1984). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1992) и доктора физико-математических наук (2008). По совместительству профессор кафедры динамического моделирования и биомедицинской инженерии СГУ. Имеет более 100 научных публикаций.
410019 Саратов, ул. Зеленая, 38
Саратовский филиал ИРЭ имени В.А. Котельникова РАН E-mail: ponomarenkovi@rambler.ru
Селезнев Евгений Петрович - родился в 1960 году в Саратове, окончил Саратовский государственный университет (1982). После окончания СГУ работает в Саратовском филиале Института радиотехники и электроники имени В.А. Котельникова РАН, в настоящее время зам. директора по научной работе. Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1990) и доктора физико-математических наук (2006). Область научных интересов - радиофизика, экспериментальное исследование нелинейных явлений. Имеет более 150 научных публикаций.
410019 Саратов, ул. Зеленая, 38
Саратовский филиал ИРЭ имени В.А. Котельникова РАН E-mail: evgenii_seleznev@mail.ru