Исследование индуцированной шумом синхронизации в системах с дискретным временем

Федеральное агентство по образованию

ГОУ ВПО «САРАТОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИМЕНИ Н.Г.ЧЕРНЫШЕВСКОГО»

Кафедра нелинейной физики

КУРСОВАЯ РАБОТА

Исследование индуцированной шумом синхронизации в системах с дискретным временем

студента 3 курса факультета нелинейных процессов

Костакова Алексея Александровича

Научный руководитель

ассистент ______________________ О.И. Москаленко

Зав. кафедрой,

профессор, д.ф.-м.н. ______________________ Ю.П. Шараевский

Саратов – 2008

Содержание

Введение

Синхронизация колебаний

Цель работы

Синхронизация, индуцированная шумом

Численное моделирование

Вывод

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

Введение

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

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

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

Синхронизация колебаний

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

С развитием теории динамического хаоса было выявлено достаточно различных типов хаотического синхронного поведения связанных динамических систем:

    фазовая синхронизация

    обобщенная синхронизация

    лаг-синхронизация

    перемежающаяся фазовая синхронизация

    перемежающиеся лаг-синхронизация

    перемежающаяся обобщенная синхронизация

    полная синхронизация

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

Цель работы

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

1.1, где

2. 2, где

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

Синхронизация, индуцированная шумом

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

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

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

Механизмы возникновения

Возможны два похожих механизма, приводящих к возникновению режима индуцированной шумом синхронизации:

    Случайный сигнал имеет ненулевое среднее, что фактически переводит систему в нехаотический режим, при котором состояние системы просто ‘следует’ за внешним случайным возмущением .

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

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

Связь обобщенной синхронизации и синхронизации,

индуцированной шумом

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

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

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

Численное моделирование

Описание рассмотренных систем

    Логистическое отображение под воздействием шума:

3, где (1)

Значение управляющего параметра , - параметр связи.

Случайная величина подчиняется нормальному распределению , где , .

Бифуркационная диаграмма для данного отображения имеет вид:

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

    Одномерное отображение вида:

4, где (2)

Значение управляющего параметра , - параметр связи

Случайная величина подчиняется нормальному распределению , где,.

Бифуркационная диаграмма для данного отображения имеет вид:

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

Результаты, полученные с помощью созданной программы

1. Для отображения , где при

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

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

С увеличением параметра связи : точки соответствующие состояниям систем, лежаться на диагональ y=z, что свидетельствует о наличии синхронного поведения в системе.

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

    Для отображения , где , при получаем аналогичные результаты: при синхронизации не наблюдается:

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

Но с увеличением параметра связи ε=0.2 появляется функциональная зависимость, что свидетельствует об установлении режима индуцированной шумом синхронизации.

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100bf4b0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004e2b00001e16000020454d46000001001803000012000000020000000000000000000000000000007f120000771a0000c80000001f010000000000000000000000000000000f030058600400160000000c000000180000000a00000010000000000000000000000009000000100000003b0a000039050000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001800c09518001000000024991800a4961800e4503565249918001c961800100000008c971800089918008a4f3565249918001c9618002000000080e03c671c9618002499180020000000fffffffffc137b01fbe03c67ffffffffffff0180ffff0180ffff0180ffffffff00000000000000000000000064f46e0501000000000000005802000025000000372e9001cc00020f0502020204030204ff0200e1ffac004009000000000000009f01000000000000430061006c0069006200720000000000a0a4170441a135657dbe4caddc137b01d043830050961800d532366720000000010000008c9618008c9618000076346720000000b4961800fc137b016476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000e7298740a48e87400000000057000000010000004c0000000400000000000000000000003b0a00003a05000050000000200035003600000046000000280000001c0000004744494302000000ffffffffffffffff3c0a00003a050000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c02a1003a01040000002e0118001c000000fb020200010000000000bc02000000cc0102022253797374656d003f00003f3f3f3f3f3f3f3f3f3f0800000001003f3f3f3f3f00040000002d010000040000002d01000004000000020101001c000000fb02f5ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0a00000001000400000000003901a10020000600040000002d010000040000002d010000030000000000

С помощью данной программы было найдено, что порог синхронизации индуцированной шумом:

-для первого отображения

-для второго отображения

Ляпуновские экспоненты

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

Для отображений ляпуновский показатель рассчитывается по формуле:

5,

где F(x) – функция, задающая отображение.

Для рассматриваемых систем зависимость ляпуновской экспоненты от управляющего параметра имеет вид:

1. , где

2. , где

Видно, что для логистического отображения (1) ляпуновская экспонента становится отрицательной при  = 1.165, для отображения (2) – при  = 1.151.Таким образом, результаты, полученные при помощи обоих методов диагностики, оказываются приблизительно одинаковыми.

Выводы

Было изучено явление индуцированной шумом синхронизации в системах с дискретным временем. Для диагностики синхронного режима производилось непосредственное сравнение векторов состояния идентичных систем, на которые воздействовал один и тот же источник шума, а также производился расчет условных ляпуновских экспонент. Рассмотрена взаимосвязь индуцированной шумом синхронизации с обобщенной синхронизацией. Была создана программа, иллюстрирующая явление индуцированной шумом синхронизации. С помощью этой программы рассмотрены два отображения. Также для этих отображений получены зависимости ляпуновской экспоненты от управляющего параметра. Полученные результаты хорошо согласуются с результатами работ [1-3].

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

    А.А. Короновский, О.И. Москаленко, А.Е. Храмов “О механизмах, приводящих к установлению режима обощенной синхронизации”, ЖТФ, 76, 2 (2006) 1-9.

    Raul Toral, Claudio R. Mirasso, E. Hernandez-Garcia and Oreste Piro “Analytical and Numerical Studies of Noise-induced Synchronization of Chaotic Systems”, CHAOS, 11, 3 (2001) 665-673.

    A.E. Hramov, A.A. Koronovskii, O.I. Moskalenko “Are generalized synchronization a noise-induced synchronization identical types of synchronous behavior of chaotic oscillators”, Phys. Lett. A, 354, 5-6 (2006) 423-427.

    С.П. Кузнецов Динамический хаос

    Amos Martian, Jayanth R. Banavar “Chaos, Noise, and Synchronization”, Phys. Rev. letters, volume 72, number 10 (1994) 1451-1454

1 Отображение взято из работы [1]

2 Отображение взято из работы [2]

3 Отображение взято из работы [1]

4 Отображение взято из работы [2]

5 Взято из [4]