СОДЕРЖАНИЕ
Введение........................................................ 5
Глава 1. Критические явления в автономных системах с удвоениями периода в присутствии шума........................................ 16
1.1 Логистическое отображение как универсальная модель однопараметрической динамики (обзор основных свойств)...... 17
1.1.1 Бифуркационное дерево.............................. 17
1.1.2 График ляпуновского показателя..................... 19
1.1.3 Уравнение РГ, константы Фейгенбаума................ 20
1.1.4 Скейлинг на бифуркационном дереве и на графике ляпуновского показателя.................................. 23
1.2 Шум в системах с удвоениями периода...................... 26
1.2.1 Влияние шума на динамику отображений с квадратичным экстремумом. Уравнение РГ с учетом случайного внешнего воздействия (краткий обзор)..................... 26
1.2.2 Скейлинг в присутствии шума. Иллюстрация универсальности свойства самоподобия по отношению к виду шумового воздействия..................................... 29
1.2.3 Методы оценки универсальных констант шума для унимодальных отображений (краткий обзор)................. 30
1.2.4 Оригинальный численный метод, его содержание и результаты............................................... 36
1.3 Случай двухпараметрических систем - трикритическая динамика................................................... 38
1.3.1 Кубическое отображение как типичное двухпараметрическое отображение (обзор основных свойств)............ 38
1.3.2 Свойства скейлинга кубического отображения в отсутствие и в присутствии шума............................ 44
1.3.3 Трикритическая динамика с шумом на примере нели-
2
нейного осциллятора под действием случайных импульсов. 50
1.4 Случай трехпараметрических систем - динамика типа Т\3),
5и Е........................................................... 58
1.4.1 Отображение четвертой степени как типичное трехпараметрическое отображение. Классификация критических ситуаций коразмерности три (обзор основных свойств)................................................... 58
1.4.2 Трехпараметрический анализ: карты динамических режимов и карты ляпуновского показателя. Скейлинг в отсутствие и в присутствии шума........................... 61
Выводы.................................................. 78
Глава 2. Критические явления в однонаправлено связанных системах с шумом........................................................ 81
2.1 Бикритическая точка. Скейлинг в отсутствие внешнего воздействия.................................................. 82
2.2 Влияние внешнего шума на бикритическую динамику 88
2.3 Численный метод поиска констант шума для двумерных отображений............................................... 92
2.4 Ренормгрупповой анализ бикритичности в присутствии шума........................................................ 97
2.5 Иллюстрации свойства скейлинга для бикритического случая в присутствии шума................................... 104
Выводы........................................................ 108
Глава 3. Сложная динамика, критические явления и влияние шума на связанные системы с удвоениями периода......................... 110
3.1 Неидентичные связанные логистические отображения 112
3.2 Неидентичные связанные осцилляторы Дуффинга под внешним гармоническим воздействием....................... 123
3.3 Идентичные связанные логистические отображения с двумя типами связи под воздействием внешнего шума.............. 127
3
%
Выводы......................................................... 142
Заключение......................................................... 144
Литература......................................................... 147
Список публикаций по теме диссертации............................ 154
*
*
Л
4
ВВЕДЕНИЕ
Актуальность работы В численных и натурных экспериментах для огромного количества самых разнообразных динамических систем, начиная с простейших моделей с дискретным временем (отображений) и заканчивая распределенными системами, наблюдается переход к хаосу через каскад бифуркаций удвоения периода. С появлением работ Фейгенбаума [1-3] стало ясно, что проблема описания перехода к хаосу имеет не только качественный, но и количественный аспект. Системы, описываемые различными математическими моделями (дискретные отображения, дифференциальные уравнения, уравнения с запаздыванием и т.д.), различной физической природы (радиофизические, оптические, биологические и др.), если они относятся к одному определенному классу универсальности, вблизи порога хаоса демонстрируют свойства самоподобия, характеризующиеся одинаковыми универсальными масштабными константами. Теоретическим инструментом для анализа таких свойств является метод ренормализационной группы (РГ), впервые привнесенный в нелинейную динамику М. Фейгенбаумом в контексте перехода к хаосу через удвоения периода и затем развитый многими другими авторами в приложении к разным классам универсальности [4-9].
Чтобы распространить концепцию универсальности и масштабного самоподобия (скейлинга) на реальные системы или наблюдать скейлинговые закономерности в физическом эксперименте, необходимо учитывать влияние неизбежно присутствующего шума. Проблема изучения воздействия шумов на различные системы традиционно является одной из центральных в радиофизике, не только не утратившей своей актуальности, но привлекающей все большее внимание в связи с обнаружением новых явлений в системах с детерминированным хаосом и шумом [10, 11]. Таким образом, логичной и своевременной является постановка вопроса, касающегося влияния шума на критическую динамику, связанную с различными классами универсальности.
5
В этом контексте решение проблемы требует соответствующей разработки РГ анализа. Для типа критичности, изученного Фейгенбаумом, такой подход был развит в работах Кратчфилда с соавторами [12] и Шраймана с соавторами [13]. Как показано этими авторами, чтобы наблюдать каждый следующий уровень удвоения периода, необходимо уменьшать интенсивность шума на фактор /1/г = 6.61903... - универсальную константу, ответственную за скейлинговые закономерности, связанные с влиянием шума в фей-генбаумовском классе универсальности. Были предприняты аналогичные исследования и получены константы скейлинга для перемежаемости [14-17] и квазипериодичности [18] в диссипативных системах, для удвоений периода [19] и разрушения КАМ-тора [20] в гамильтоновых системах. Стоит отметить, что константа Кратчфилда и Шраймана Д/т, по-видимому, имеет более фундаментальное значение, поскольку появляется в задачах о квантовом хаосе как константа перенормировки постоянной Планка при изучении свойств квантовой системы, демонстрирующей переход в классическом пределе к системе с удвоениями периода [21-23].
В контексте многопараметрического анализа перехода к хаосу, критичность определенного типа может иметь место на некоторых поверхностях, кривых или точках в пространстве параметров. РГ анализ - неоценимый инструмент для поиска, изучения и классификации критических ситуаций. Их коразмерность* определяется числом соответствующих неустойчивых собственных векторов линеаризованного уравнения РГ. В последнее время было обнаружено множество типов критического поведения, связанных с удвоениями периода, каждый из которых допускает описание в духе теории Фей-генбаума [9]. Если сценарий Фейгенбаума является типичным однопараметрическим феноменом, то эти новые типы поведения наблюдаются при вариации двух и более параметров. Часть из них встречается в одномерных отображениях, а некоторые характерны лишь для многомерных систем. Таким
* Коразмерность - минимальное число параметров, которое должна иметь динамическая система, чтобы данная критическая ситуация могла реализоваться в типичном случае.
6
образом, возникает проблема исследования воздействия шума на системы, демонстрирующие эти многопараметрические варианты перехода к хаосу.
Объектами исследования будут служить одномерные отображения и системы связанных отображений. К настоящему времени отображения стали популярными и эффектными моделями для исследования различных радиофизических систем, таких как, например, нелинейный колебательный контур [24]; модель Улама, предложенная для объяснения механизма стохастического ускорения космических частиц Ферми [25] и т.д. Благодаря методу сечений Пуанкаре отображения применяются и для анализа радиофизических систем, описываемых дифференциальными уравнениями (генератор с инерционной нелинейностью [26] и др.). Как и традиционно для радиофизических методов, отображения стали универсальным инструментом исследования систем различной физической природы*. Так, например, хорошо изучено отображение Икеды, описывающее кольцевой оптический резонатор со средой с фазовой нелинейностью, возбуждаемый лазером [27]. (К отображению Икеды в рамках метода медленно меняющихся амплитуд приводит и задача о колебаниях осциллятора Дуффинга, возбуждаемого короткими импульсами [28, 29]. Б свою очередь, при определенных приближениях отображение Икеды сводится к одномерному “отображению косинуса”, которое является простым, физически мотивированным отображением со многими экстремумами [28,29].) Одномерные отображения привлекательны и в биологии в задачах о динамике популяций. Более того, во многом из биологии пришли вопросы о сложной динамике простейшего отображения с квадратичным экстремумом (логистическое отображение), послужившие стимулом для создания теории Фейгенбаума. Модели в виде отображений, и, в частности, связанных отображений, привлекают внимание и при анализе химических систем (см., например, работы [30, 31] и цитированную в них литературу.)
* Следует особо отметить значительный вклад нижегородской научной школы как в теорию дискретных отображений, так и в конкретные приложения таких систем.
7
В работе будут использоваться одномерные отображения достаточно общего вида, допускающие разложение в ряд Тейлора. В зависимости от числа членов, оставленных в разложении, возникают одно-, двух- и трехпараметрические модели - логистическое отображение, кубическое отображение, “квартичное” отображение (отображение четвертой степени). Эти модели оказываются достаточно универсальными, что находит свое выражение в картине колебаний, бифуркаций и перехода к хаосу, описываемых ими. “Объекты” на плоскости параметров таких отображений (точки сборки, линии складок, линии удвоений, каскады бифуркаций, известные как “crossroad area” композиции бифуркаций [32] и т.д.) легко узнаваемы теми, кто исследовал сложную динамику нелинейных колебательных контуров, различных осцилляторов и других подобных систем. Сценарии перехода к хаосу для такой многопараметрической системы отображений, как и в теории бифуркаций, могут быть классифицированы по числу существенных параметров (коразмерности), причем, первую “строку” в такой классификации занимает однопараметрический сценарий Фейгенбаума [1-3]. При этом увеличение числа учтенных членов ряда Тейлора приводит к увеличению числа экстремумов отображений, возникновению точек перегиба, что является весьма существенным в классификации сценариев (типов критического поведения) на пороге хаоса. Во второй и третьей главах мы обратимся к моделям в виде связанных логистических отображений, а в третьей главе будет рассмотрена также и дифференциальная система в виде связанных неавтономных нелинейных осцилляторов. Заметим, наконец, что шум в дискретные отображения вводится в виде аддитивной добавки. Если говорить о физической модели, где это можно сделать наиболее явным и наглядным образом, то можно представлять себе нелинейный осциллятор, возбуждаемый импульсами, амплитуда которых является слегка модулированной некоторым случайным образом.
Цель работы состоит в исследовании влияния шума и выявлении связанных с ним свойств универсальности и скейлинга для динамических сис-
8
тем, которые характеризуются наличием перехода к хаосу через последовательность бифуркаций удвоения периода.
Научная новизна работы
1. Предложен численный метод, позволяющий определять скейлинго-вые константы шума для отображений с удвоениями периода, принадлежащих различным классам универсальности, и допускающий обобщение на многомерный случай. Метод имеет то преимущество, что собственно уравнение РГ решать не требуется; при этом получаемая точность приемлема для компьютерных (а, тем более, натурных) иллюстраций свойств скейлинга на большом количестве уровней иерархии.
2. Даны иллюстрации масштабного самоподобия (скейлинга) в присутствии шума для одномерных отображений с одним, двумя и тремя управляющими параметрами. Продемонстрирован скейлинг на бифуркационных деревьях, графиках ляпуновского показателя, картах ляпуновского показателя.
3. Развит ренормгрупповой анализ влияния шума на бикритическое поведение системы двух связанных логистических отображений с однонаправленной связью. Найдена новая универсальная постоянная цв = 2.713695..., определяющая правило пересчета интенсивности шума во второй подсистеме.
4. Показано, что пространство управляющих параметров связанных неидентичных логистических отображений характеризуется наличием критических точек коразмерности два (так называемый С - тип критичности), которые представляют собой концевые точки фейгенбаумовских линий. Также отмечено определенное сходство между сложным устройством дискретной (два связанных неидентичных квадратичных отображения) и потоковой (два связанных неидентичных осциллятора Дуффинга под внешним гармоническим воздействием) систем.
5. Представлены иллюстрации скейлинга на картах ляпуновского показателя для двух связанных идентичных квадратичных отображений под воз-
9
действием шума с двумя разными типами связи: диссипативным и инерционным. Показано, что в этом случае связанных подсистем константой перенормировки интенсивности шума является универсальная постоянная Кратч-филда и Шраймана Ц/г, которая была предложена для одного автономного отображения [12,13].
Достоверность научных выводов работы подтверждается согласованностью и воспроизводимостью всех данных, а также хорошим совпадением результатов, полученных методом ренормализационной группы и компьютерными иллюстрациями масштабного подобия (скейлинга).
Основные положения, выносимые на защиту
1. Результаты работы выявили количественные закономерности, связанные с воздействием шума на системы с хаотической динамикой вблизи различных типов точек перехода к хаосу через удвоения периода. Системы, описываемые двух- и трехпараметрическими отображениями с шумом, вблизи соответствующих критических точек перехода порядок-хаос демонстрируют свойство самоподобия (скейлинга) на бифуркационных деревьях, графиках и сечениях пространства параметров в виде карт ляпуновского показателя с константами \ьт = 8.2439, р$ = 10.0378, |1£= 11.5937, ответственными за пересчет амплитуды шума.
2. Результаты полученного в работе ренормгруппового анализа для системы двух связанных отображений с однонаправленной связью в присутствии шума показывают, что характерному для такой системы случаю бикритической динамики отвечает новая универсальная постоянная рд = 2.713695..., определяющая правило пересчета интенсивности шума.
3. Устройство пространства управляющих параметров двух неидентичных связанных логистических отображений характеризуется наличием критических точек коразмерности два, известных как С — тип критичности. Для идентичных связанных логистических отображений с шу-
10
мом пространство управляющих параметров (параметр нелинейности, константа связи, амплитуда шума) демонстрирует свойство скейлинга. Научно-практическая значимость работы. Рекомендации по использованию научных выводов
Результаты работы, связанные с возможностью вычисления констант шума с помощью предложенного численного метода, открывают перспективы по исследованию влияния флуктуаций и демонстрации скейлинговых свойств для типов критической динамики, как связанных с переходом к хаосу через последовательность удвоений периода, так и для других известных сценариев.
Результаты ренормгруппового анализа бикритичности дают возможность количественного описания влияния шума на связанные системы различной природы с однонаправленной связью. Можно ожидать, что для оставшихся за рамками настоящей работы типов критичности на пороге перехода порядок-хаос также возможно ренормгрупповое описание, учитывающее влияние шума.
Установленное свойство скейлинга в системах с шумом и найденные константы шума позволяют оценивать число удвоений, для которых традиционная картина перехода к хаосу окажется разрушенной. В то же время можно рекомендовать проведение специальных экспериментов, нацеленных на проверку обнаруженных закономерностей, которые требуют, однако, использование специальных генераторов шума с регулируемым его уровнем.
Результаты исследования системы неидентичных связанных отображений и связанных нелинейных осцилляторов выявили тонкую картину бифуркаций в области сосуществования двух разных сценариев перехода к хаосу -через каскад бифуркаций удвоения периода и через разрушение квазиперио-дического режима. Обнаружено, что такой переход в связанных системах может ассоциироваться с типом критичности, известным, как С - тип. В силу универсальности, присущей критической динамике, следует ожидать, что в системах различной природы с подобным типом связи может наблюдаться
11
аналогичная картина бифуркаций, включая соответствующую тонкую картину на глубоких уровнях иерархии в непосредственной окрестности критической точки.
Результаты работы использованы в учебном процессе на факультете нелинейных процессов Саратовского госуниверситета.
Апробация работы и публикации Результаты работы были представлены на школе-конференции “Нелинейные дни в Саратове для молодых-98” (Саратов, 1998), VI международной конференции студентов и аспирантов по фундаментальным наукам “Ломоносов-99” (Россия, Москва, 1999), Международной школе для молодых ученых и студентов по оптике, лазерной физике и биофизике (Россия, Саратов, 1999), школе - конференции “Нелинейные дни в Саратове для молодых-99” (Саратов, 1999), VII Всероссийской школе -семинаре “Волновые явления в неоднородных средах” (Красновидово Московской области, 2000), всероссийской научной школе - конференции “Нелинейные дни в Саратове для молодых-2000” (Саратов, 2000), международной межвузовской конференции “Современные проблемы электроники и радиофизики СВЧ” (Россия, Саратов, 2001), четвертом международном симпозиуме по классической и небесной механике (Россия, Великие Луки, 2001), IX Международном симпозиуме “Nonlinear dynamics and complex systems” (Республика Беларусь, Минск, 2001), шестой международной школе “Хаотические автоколебания и образование структур” (Россия, Саратов, 2001), всероссийской научной школе - конференции “Нелинейные дни в Саратове для молодых-2001” (Саратов, 2001), конференциях молодых ученых и аспирантов научно-образовательного центра “Нелинейная динамика и биофизика” СГУ (Саратов, 2002,2003,2004), Международной конференции по синхронизации хаотических и стохастических колебаний (Приложения в физике, химии, биологии и медицине) “Synchro-2002” (Россия, Саратов, 2002), всероссийской научной школе - конференции “Нелинейные дни в Саратове для молодых-2003” (Саратов, 2003).
12
По теме диссертации имеется 17 публикаций (5 статей в российских и зарубежных реферируемых журналах, 7 статей в сборниках трудов научных конференций, 5 тезисов докладов).
Результаты работы получены в рамках грантов Российского фонда фундаментальных исследований (№ 96 - 15 - 9692, 97 - 02 - 16414, 00 - 02 - 17509, 03 - 02 - 16074, в том числе персонального гранта для молодых исследователей № 01 - 02 - 06388), Американского фонда гражданских исследований и развития CRDF (№ REC - 006), Эйлеровой программы (совместно с группой профессора А. Пиковского, университет Потсдама, Германия), Министерства общего и профессионального образования (№ 97 - 0 - 8.3 - 8.8), Федеральной целевой программы “Интеграция” (№ 696.3, проект А0057) и фонда некоммерческих программ “Династия” при содействии Международного центра фундаментальной физики в Москве.
Личный вклад автора В совместных работах автором выполнено программирование всех задач и проведены численные эксперименты. Формулировка поставленных задач, выбор моделей, а также объяснение и интерпретация полученных результатов проведено совместно с научным руководителем и соавторами. В работе [97] автору принадлежит п.5.
Стру»стура и объем работы Диссертация состоит из введения, 3-х глав и заключения. Работа содержит 158 страниц текста, включая 70 рисунков, 3 таблицы, список литературы из 109 наименований на 11 страницах.
Краткое содержание работы
Во введении обсуждается актуальность работы, формулируется цель исследования, научная новизна и основные положения, выносимые на защиту.
В первой главе в п. 1.1 дан обзор основных свойств логистического отображения - базовой модели для изучения перехода к хаосу по сценарию Фейгенбаума. Обсуждается устройство бифуркационного дерева и графика ляпуновского показателя для дискретных отображений, дано определение скейлинга, представлено краткое изложение метода ренормализационной
13
- Київ+380960830922