Синхронизация хаотически модулированных колебаний Текст научной статьи по специальности «Физика»

Теория колебаний

Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 5 (3), с. 292-299

УДК 621.391.1

СИНХРОНИЗАЦИЯ ХАОТИЧЕСКИ МОДУЛИРОВАННЫХ КОЛЕБАНИЙ © 2011 г. В.В. Матросов, В.Д. Шалфеев

Нижегородский госуниверситет им. Н.И. Лобачевского

matrosov@rf.unn.ru

Поступила в редакцию 30.05.2011

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

Ключевые слова: динамический хаос, синхронизация, фазовые системы, автоподстройка.

Введение

Явление синхронизации хаотических колебаний было открыто около трех десятилетий назад [1, 2]. В последующие годы появился огромный поток публикаций по этой проблеме, тем не менее сегодня можно констатировать, что достаточно удовлетворительной теории синхронизации хаотических колебаний пока не существует. Трудности построения стройной теории синхронизации хаотических колебаний в значительной мере обусловлены размытостью понятия фазы и частоты хаотических сигналов и следующей отсюда неопределенностью того, что целесообразно считать явлением синхронизации хаотических колебаний, что является критерием синхронизации и т.д. Важно также отметить, что развитие теории синхронизации хаотических колебаний сейчас происходит на фоне крайне слабой заинтересованности в построении такой теории со стороны физических и технических приложений и, следовательно, при отсутствии четких целевых установок для развития теории; поэтому, к сожалению, многие задачи по синхронизации хаотических колебаний формулируются исходя из модельных представлений, а не из требований практических приложений.

В начале 90-х годов появилась яркая идея практического использования динамического хаоса для скрытой передачи информации [3]. Идея оказалась весьма привлекательной. Теоретические исследования и макетные разработки, выполненные в 90-х годах, подтвердили принципиальную возможность реализации систем передачи информации с использованием хаотических колебаний в качестве несущих. Однако в

различных макетах для создания в приемнике копии хаотических колебаний, использованных при передаче в качестве несущих (что принципиально необходимо при когерентном приеме), потребовалось выполнить условие высокой идентичности параметров генераторов хаоса в передатчике и приемнике, что привело к низкой помехоустойчивости и слабой криптостойкости разработанных систем. Дальнейшее развитие этого направления привело к заключению о бесперспективности создания конкурентоспособных приемопередающих устройств на основе динамического хаоса, реализующих когерентный способ приема [4]. Одним из вариантов альтернативного подхода к решению проблемы является использование хаотически модулированных колебаний (ХМК) при передаче информации [5]. Переход к ХМК, вообще говоря, может привести к некоторому ухудшению свойств широкополосности хаотического сигнала, однако несомненным преимуществом станет возможность решения задачи о надежной синхронизации ХМК. В [6, 7] показано, что на основе свойств коллективной динамики ансамблей связанных систем фазовой автоподстройки (ФАП) можно построить высокоэффективные генераторы хаотически модулированных колебаний со стабилизированной по опорному колебанию центральной частотой. Рассмотрим далее задачу о синхронизации таких колебаний.

Схема и динамическая модель синхронизации

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

Рис. 1. Схема синхронизации хаотических колебаний генераторов Г! и Г2

двух систем ФАЩ и ФАП2 (рис. 1). Для достижения синхронизации ФАП и ФАП2 связаны через дополнительный частотный дискриминатор ЧД, на котором формируется сигнал Ф(у2-уі) для управления частотой генератора Г2. Математической моделью таким образом связанных ФАП1 и ФАП2 является система уравнений:

йф

йх

= Уи

аУі

йх

Ці — = Уі - в1пФ! - У!-Є1г1,

ах

йФ2

йх

й2-,

= У2,

йУ2

йх

(1)

= 2 г,

Ц 2—^ = у 2 - в1п ф2 йх

- У2 -Є 2 22 -^ Ф(У2 - У1 )

где первые три уравнения описывают динамику ФАПЬ последние три уравнения - ФАП2 с дополнительным управлением, нелинейность частотного дискриминатора аппроксимируется функцией Ф(у)=2ау/(1+а2у2). Здесь ф! и ф2 (у и у2) - текущие разности фаз (частот) колебаний генераторов Г1 и Г2 относительно опорного сигнала 80(т); безразмерные параметры у1 и у2 характеризуют начальные частотные расстройки колебаний генераторов относительно S0(т), еь ц1 и е2, ц2 - параметры фильтров Ф1 и Ф2, 8 -силу связи (управления).

При моделировании динамики (1) параметры Уь е1, Ц1 выбирались такими, что сигнал S1(т) являлся хаотически модулированным со стабилизированной по опорному сигналу S0(т) средней частотой, параметры е1=е2, у^у2, Ц1^Ц2, таким образом, сигнал S2(т) в общем случае отличается от S1(x). При отсутствии связи (8 =0) да-

же при близких параметрах ФАП1 и ФАП2 синхронизации колебаний ФАП1 (рис. 2а) и ФАП2 (рис. 2б) нет (рис. 2в). Введением связи (8>0) можно добиться синхронизации ХМК (рис. 2г), при этом силой связи можно обеспечить уменьшение ошибки синхронизации е до малых величин (рис. 2д), т.е. осуществить с заданной точностью синхронизацию ХМК1 с выхода ФАП1 и ХМК2 с выхода ФАП2.

Динамические режимы и характеристики синхронизации хаотически модулированных колебаний

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

Режим синхронизации хаотически модулированных колебаний. В фазовом пространстве модели (1) режиму синхронизации отвечает некоторое устойчивое многообразие ^ - многообразие синхронизации. В случае синхронизации хаотических колебаний идентичных систем

(в рассматриваемом случае это означает равенство в модели (1) параметров ФАП1 и ФАП2: У1=У2, е1=е2, ц1=ц2 ) проекции ^ на плоскости соответствующих фазовых переменных (ф1,ф2), (у1,у2), (г1гг2) представляют собой отрезки,

размещенные на диагоналях. Проекции ^ на плоскости разностей соответствующих фазовых переменных [(ф 2-ф 1>,С>’2 ->’1)], [(Ф2-Ф1),

(22-21)], [С>2->’1),(^2-^1>] есть точки.

В случае неидентичных систем проекции многообразия ^ на плоскости соответствующих фазовых переменных всегда обладают некоторой «толщиной» (рис. За, б, в), а на плоскостях разностей соответствующих фазовых переменных вместо точки появляется некоторое множество, по виду напоминающее проекции колебательного хаотического аттрактора, у которого фазовые переменные изменяются в небольших пределах (рис. 3г, д, е). Следует обратить внимание на то, что в случае неидентичных систем некоторые проекции на плоскости соответствующих фазовых переменных могут размещаться не в окрестности самих диагоналей, а в окрестности прямых, параллельных этим диагоналям, как на рис. 3 а. В этом случае синхронизация наступает с точностью до констант, определяемых величинами сдвига прямых относительно диагоналей. Значения констант определяют остаточные ошибки синхронизации по соответствующим переменным.

Точность хаотической синхронизации. Толщина многообразия ^ или е - радиус окружности, которой можно окружить проекцию на

плоскостях разностей фазовых переменных, характеризует точность синхронизации. Эта величина имеет особое значение, поскольку, в случае практического использования явления хаотической синхронизации, именно она характеризует качество синхронизации. В задачах передачи информации с когерентным приемом от точности синхронизации зависит точность воспроизведения принимаемой информации.

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

Область захвата в режим хаотической синхронизации. Это область пространства параметров, где режим синхронизации наступает гарантированно при любых или при определенных, характерных для каждой конкретной модели, начальных условиях. Отсюда можно предложить несколько подходов к определению границ области захвата: граница определяется либо условиями глобальной устойчивости многообразия ^ , либо условием принадлежности набора характерных начальных условий бассейну притяжения ^ . Задача об определении области захвата в режим хаотической синхронизации, так же как и в случае определения области удержания режима

2.8

Фг

-0.3

/ 1.0 / 0.6 /

/ / /

-0.7 / -0.5 /

-0.3 ф1 (а)

0.05

§

-0.05

01 Ф-Ф1

(г)

-О?

(б)

1.0

-0.5

0.6

(в)

0.05 0.05

(§) € N -0.05 (Ч») £ -0.05 Ф

0.3

0.1 ф-ф! 0.3

(д)

-0.05 уг-ух 0.05

(е)

Рис. 3. Проекции многообразия синхронизации в случае неидентичных систем. Пунктирная линия воспроизводит 8-окружность, характеризующую точность синхронизации

хаотической синхронизации, должна решаться с учетом точности синхронизации.

Время установления режима синхронизации. Это время переходного процесса от начального состояния до установления в системе режима хаотической синхронизации, т.е. время достижения фазовой траекторией некоторой 8-окрестности многообразия ^.

Важное место в исследованиях хаотической синхронизации занимает изучение сценариев установления и срыва режима синхронизации, т.к. понимание этих механизмов позволяет не только разрабатывать оптимальные алгоритмы расчета динамических характеристик систем синхронизации (областей захвата и удержания, точности синхронизации, времени установления режима синхронизации и т.д.), но и минимизировать ошибки этих расчетов.

Моделирование процессов синхронизации ХМК

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

рассчитана по следующему алгоритму. Изначально в модели (1) значения параметров и начальные состояния ФАП и ФАП2 выбираются одинаковыми. В модели (1) устанавливается режим хаотической синхронизации1. Далее значения параметров системы ФАП1 остаются неизменными, а значения параметров ФАП2 меняются, удаляясь от значений параметров ФАП1. При этом начальные состояния модели (1) для новых значений параметров наследуются с предыдущего шага. Рассогласования параметров ФАП2 и ФАП1 увеличиваются до тех пор, пока не произойдет срыв синхронизации.

Расчет границы области захвата осуществляется в обратном порядке. Теперь изначально значения параметров и начальное состояние модели (1) выбираются таким образом, чтобы ФАП2 и ФАП1 находились в асинхронном режиме. Далее параметры системы ФАП2 меняются в сторону уменьшения рассогласований между значениями параметров ФАП1 и ФАП2, начальные состояния модели (1) при этом наследуются. Границе области захвата соответствуют такие значения параметров, когда в модели (1) вновь установится режим хаотической синхронизации.

Если нанести границы областей удержания режима хаотической синхронизации на карту динамических режимов автономной модели ФАП2 [1, 3], то получим рис. 5, из которого видно, что области синхронизации модели (1) могут

Рис. 4. Области захвата и удержания режима хаотической синхронизации в модели (1) при 8! = е2 = =1, щ = 2.38, у! = 0.625, 5 = 0.5, а = 30: точка Ь иллюстрирует значения параметров ФАП1, линия 1 есть граница области удержания режима хаотической синхронизации (точность синхронизации 8=0.2), линия 2 - граница области захвата, линия 3 -граница области синхронизации с повышенной точностью (точность синхронизации 8=0.01)

Рис. 5. Области существования режима хаотической синхронизации модели (1) в случае є2=1, 5=0.1; 0.5; 0.9, а=30 для точек Ь (у1 =0.625, є1=1, ц1=2.38) и а (у1=0.95, є1=1, ц1=6.5). Здесь границы областей синхронизации для точки Ь проведены сплошными линиями, для точки а - пунктирными линиями

захватывать области существования различных динамических режимов автономной модели ФАП2. Отсюда следует, что при разомкнутом кольце автоподстройки (5=0) ФАП2 может функционировать практически в любом из возможных динамических режимов: от режима генерации регулярных колебаний, синхронизированных с опорным сигналом, до режима генерации хаотически модулированных колебаний. Цепи автоподстройки здесь не только позволяют скомпенсировать существенную разницу между двумя колебательными процессами, но и обеспечивают поддержание режима хаотической синхронизации с высокой точностью на протяжении длительного времени. Заметим, что на установление режима хаотической синхронизации существенное влияние оказывают как характеристики самих колебательных движений, так и параметры цепей автоподстройки. В примере, представленном на рис. 5, хаотические колебания в цепи управления ФАП1 для значений параметров, определенных в точках а и Ь, различаются амплитудами, в частности проекции аттрактора модели (1) на плоскость (ф^) в точке а меньше, чем в точке Ь. Из вышесказанного следует, что области синхронизации ХМК двух однонаправленно связанных систем ФАП1 и ФАП2 увеличиваются при увеличении силы связи 5 и при уменьшении размеров хаотического аттрактора ведущей подсистемы.

Теперь рассмотрим процессы установления режима хаотической синхронизации, когда ФАП1 функционирует в режиме генерации ХМК, а в ФАП2 могут реализовываться различные режимы. Пусть параметры ФАП1 определе-

ны в точке Ь (уі = 0.625, 8! = 1, ц1=2.38), а параметры ФАП2 имеют следующие значения: у2 = =0.68, 82 = 1, ц2 = 2. Тогда в отсутствие связи (5=0) в парциальном фазовом пространстве и1 = ={ф1(mod 2л), у1, х{} модели (1) имеет место колебательный хаотический аттрактор, а в и2= ={ф2(mod 2л), у2, г2} существуют два регулярных аттрактора с различной топологией - колебательный Ь1 (рис. 6а) и вращательный L2 (рис. 7 а) предельные циклы. Определим состояние ФАП2 таким образом, что в фазовом пространстве и2 реализуется цикл Ь1, соответствующий регулярному квазисинхронному режиму. Далее будем увеличивать значения параметра 5, добиваясь синхронизации колебаний ФАП1 и ФАП2. Процесс установления режима синхронизации иллюстрирует рис. 6, на котором представлены: однопараметрическая бифуркационная диаграмма модифицированного отображения Пуанкаре (рис. 6б)2, проекции фазовых портретов и осциллограммы (рис. 6в-ж). Здесь при слабой связи 5=0.05 режим синхронизации отсутствует (рис. 6в). Увеличение силы связи до 5=0.15 приводит к появлению на осциллограмме интервалов синхронных режимов, которые достаточно часто перемежаются асинхронными движениями (рис. 6г). Появление коротких интервалов синхронизации не приводит к заметным изменениям в проекциях фазового портрета аттрактора. Дальнейшее увеличение 5 приводит к увеличению интервалов синхронных движений, эти изменения находят свое отражение и на проекциях фазового портрета (рис. 6д) -изображающая фазовые траектории точка боль-

ше времени проводит в окрестности диагонали. Наконец, при 5=0.232 устанавливается режим синхронизации с точностью 8=0.04 по частоте (УгОО - у1(^)|<8) (рис. 6е). На рис. 6ж представлены проекции фазового портрета на плоскости разности фазовых переменных и соответствующая осциллограмма в увеличенном масштабе. Дальнейшее увеличение 5 приводит к повыше-

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

Процесс установления режима синхронизации с ростом 5 в случае, когда при 5=0 в фазовом пространстве и2 реализуется вращательный цикл

Ь2 (рис. 7а), иллюстрируют диаграмма отображения Пуанкаре, представленная на рис. 7б, проекции фазовых портретов и осциллограммы на рис. 7в, г. Здесь введение слабой связи 0<5<<1 приводит к тому, что колебания на выходе ФАП2 становятся хаотическими. Этому режиму в фазовом пространстве модели (1) отвечают движения по аттрактору, у которых разность между координатами ф2 и ф1 постоянно нарастает, а разность между у2 и у1 хаотически изменяется в широком диапазоне (рис. 7в). Таким образом, в случае слабой связи режим хаотической синхронизации колебаний ФАП1 и ФАП2 отсутствует. При увеличении 5 наблюдаемая выше картина не претерпевает существенных изменений вплоть до значений 5=0.234, когда аттрактор асинхронных режимов разрушается, и фазовые траектории из его окрестности переходят на многообразие синхронизации (рис. 7г). Если теперь начать уменьшать параметр 5, то выход из режима синхронизации происходит через перемежаемость синхронных и асинхронных движений при 5=0.232 .

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

новления режима синхронизации: в случае регулярного квазисинхронного режима переход на режим генерации ХМК происходит мягко, через перемежаемость, в случае регулярных биений - жестко, через разрушение асинхронного аттрактора.

Выводы

Результаты моделирования системы (1) показали, что существует принципиальная возможность подстройки (синхронизации) колебаний ФАП2 под хаотически модулированные колебания ФАП1 с помощью схемы, представленной на рис. 1. На основании проведенного исследования заключаем, что на базе кольца ФАП возможно осуществить не только генерацию хаотически модулированных колебаний в достаточно широкой области изменения параметров, но и их синхронизацию. Процессы синхронизации ХМК характеризуются наличием гистерезисности, т.е. наличием в пространстве параметров областей удержания и захвата с несовпадающими границами.

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

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

Работа выполнена при поддержке гранта РФФИ №10-02-00865, ФЦП «Научные и научнопедагогические кадры инновационной России» на 2009—2013 гг. (контракты № П2308,

№ 02.740.11.0565, № 02.740.11.0075).

Примечания

1. Вообще говоря, начальные состояния ФАП1 и ФАП2 могут быть выбраны неидентичными, но в этом случае необходимо, чтобы начальные условия принадлежали бассейну притяжения аттракторов, определяющих один и тот же динамический режим. Тогда при S>0 режим хаотической синхронизации наступит через некоторый интервал времени Ttrans, длительность которого зависит от соотношения начальных состояний ФАП1 и ФАП2. Время Ttrans всегда конечно, более того, оно сравнительно мало.

2. Отличие модифицированного отображения Пуанкаре от обычного заключается в том, что на диаграмме на рис. бб по оси ординат отложены не координаты точек отображения Пуанкаре, а значения некоторой функции R от этих координат, в частности здесь R(y2,y1)=y2-y1.

Список литературы

1. Fujisaka H., Yamada Т. II Prog. Theor. Phys. 1983. V.69. №1. P. 32-46.

2. Афраймович В.С., Веричев Н.Н., Рабинович М.И. // Изв. вузов. Радиофизика. 1986. Т.29. №9. С. 795-8Q3.

3. Pecora L.M., Carroll T.L. II Phys. Rev. Lett. 199Q. VoI.64. №.8. P. 821-824.

4. Дмитриев А.С., Клецов А.В., Лактюшкин А.М. и др. // Радиотехника и электроника. 2QQ6. Т. 51. №10. С. 1193-12Q9.

5. Шалфеев В.Д., Матросов В.В., Корзинова М.В. // Зарубежная радиоэлектроника. Успехи современной радиоэлектроники. 1998. №11. С. 44-56.

6. Shalfeev V.D., Matrosov V.V., Korzinova M.V. II Int. J. Bifurcation and Chaos. 1999. V.9. №.5. P. 963-973.

7. Матросов В.В., Шалфеев В.Д. Динамический хаос в фазовых системах. Нижний Новгород: ННГУ, 2QQ7. 258 c.

SYNCHRONIZATION OF CHAOTICALLY MODULATED OSCILLATIONS

V. V. Matrosov, V.D. Shalfeev

Synchronization of chaotically modulated oscillations is studied using the phase-locked loop principle. The application of this principle makes it possible to improve synchronization accuracy of chaotic oscillations and to enlarge the existence domains of chaotic synchronization regime in the parameter space of phase-locked oscillators.

Keywords: dynamical chaos, synchronization, phase-locked loops.