Введение 3
Постановка задачи 5
Глава 1. Некоторые элементы теории хаоса 6
1.1. Хаотические системы 6
1.2. Аттракторы 7
1.3. Показатели Ляпунова 8
Глава 2. Вычисление старшего показателя Ляпунова по временному ряду методом Вольфа 12
2.1. Восстановление аттрактора по временному ряду методом Такенса 12
2.2. Описание алгоритма 13
2.3. Пример. Аттрактор Лоренца 15
Глава 3. Вычисление старшего показателя Ляпунова для одной модели звездной динамики 18
3.1. Описание модели движения звезды 18
3.2. Реконструкция акттракторов по временным рядам псевдоэкспериментальных данных 20
3.3. Расчет старшего показателя Ляпунова для различных классов движений звезды 25
Выводы 29
Список литературы 31
При исследовании свойств различных динамических систем в настоящее время получил актуальность подход, основанный на анализе сигналов, произведенных системой. Такой метод применим в тех случаях, когда сам процесс математически описать практически невозможно, однако есть возможность наблюдать некоторую величину. Например, электрокардиограмма, запись колебаний земной коры, данные метеонаблюдений, биржевой курс, солнечная активность и т.п. Поэтому такой анализ часто реализуется путем обработки фиксируемых сигналов. Изменение со временем некоторой переменной называется (скалярным) временным рядом, а если измеряется несколько переменных — векторным, метод называется реконструкцией динамических систем, а раздел теории динамических систем — анализом временных рядов [1, 2, 3]. Временные ряды естественным образом возникают в эксперимента как натурных (тогда получается непрерывный ряд y(t)), так и численных (тогда получается дискретный ряд). Конечно, наличие только лишь наблюдаемых данных усложняет процесс изучения системы по сравнению с тем случаем, когда имеется полное ее решение. Поэтому необходимо вводить дополнительные ограничения при реализации метода реконструкции [4].
Скалярным временным рядом xi, ...,xn называется массив из N чисел, которые представляют собой значения некоторой наблюдаемой динамической переменной x(t) с постоянным по времени шагом т , где ti = to + (i — 1)т, xi = x(ti),i = I,..., N. При анализе временных рядов выделяются две основные задачи, первая из которых — задача идентификации, а вторая — задача прогноза. Задача идентификации заключается в нахождении таких параметров системы, породившей данный временной ряд, как размерность вложения, корреляционная размерность, энтропия и др. Размерностью вложения называется наименьшее число независимых переменных, которые однозначно определяют установившееся движение исходной диссипативной распределенной системы. Обычно размерность вложения обозначают как d. Подмножество, к которому с течением времени стягиваются все близлежащие траектории системы — это аттрактор [1, 5, 6]. Размерность вложения аттрактора — это минимальная размерность фазового пространства, в которое без самопересечений может быть помещено гладкое многообразие, целиком содержащее этот аттрактор. Например, двумерный тор вложим только в трехмерное пространство и пространство большей размерности. Корреляционная размерность является оценкой фрактальной размерности аттрактора системы и частным случаем обобщенной вероятностной размерности. Понятие энтропии связано с предсказуемостью значений ряда и всей системы.
Исследование временных рядов основано на идее, что для получения удовлетворительной геометрической картины странного аттрактора, вместо входящих в исходную систему переменных можно использовать векторы задержек наблюдаемой Zi = (xj,x+i,...,x+m_i). Впервые данный подход к анализу временных рядов был математически обоснован в работе Ф.Такенса [7]. То есть наиболее интересным приложением теории динамических систем может являться прогнозирование динамики временных рядов, порождаемых ими. При этом предполагается, что характеристики систем могут быть неизвестны. Таким образом, теоретические исследования, основанные на анализе временных рядов, могут дать мощный инструмент для понимания многих явлений, особенно когда имеющихся данных для построения модели может быть недостаточно.
Исследование динамики реальных процессов и явлений зачастую начинается с накопления экспериментальных данных в виде временных рядов. В последующем значения этих временных рядов должы быть обработаны с целью выявления в них каких-либо важных динамических характеристик, присущих наблюдаемому процессу. Одной из таких характеристик является мера экспоненциальной скорости разбегания траекторий движений, начавшихся в достаточно малой окрестности — старший показатель Ляпунова. Положительные значения старшего показателя Ляпунова отвечают хаотическим типам движений динамических систем, анализ которых особенно затруднителен. Тем самым, определение возможной меры хаоса в экспериментальных данных является весьма важной и интересной задачей.
Основные результаты представленной дипломной работы:
• Изучен и отработан на примере модели Лоренца алгоритм Такенса реконструкции аттрактора динамической системы по одному временному ряду экспериментальных данных.
• Методом Вольфа, по временным рядам псевдоэкспериментальнам данных для системы Лоренца (с определённым набором параметров и начальной точкой), был продублирован расчёт старшего показателя Ляпунова.
• На примере одной модели динамики движения звезды при стационарном ротационно-симметричном потенциале ' была произведена генерация псевдоэкспериментальных данных для ряда ключевых наборов параметров рассмотренной системы.
• По полученным временным рядам была произведена соответствующая реконструкция аттракторов, отвечающих основной динамике движения звезды в этой модели. Показана эквивалентность реконструированных аттракторов с реальными.
• Методом Вольфа по временным рядам псевдоэкспериментальных данных рассчитаны старшие показатели Ляпунова, правильнол отражающие наличие или отсутствие хаотических режимов поведения для предложенной модели движения звезды.
Таким образом, показана возможность идентификации меры хаотического поведения в траекториях ротационно-симметричной модели движения звезды по наблюдаемым данным временных рядов.
[1] Кузнецов С. П. Динамический хаос. М.: Физматлит, 2001. 295 с.
[2] Головко В. А. Нейросетевые методы обработки хаотических процессов // Научная сессия МИФИ-2005. VII Всероссийская научно - техническая конференция «Нейроинформатика -2005»: Лекции по нейроинформатике. М.: МИФИ, 2005. С. 43-91.
[3] Лоскутов А. Ю. Анализ временных рядов. Курс лекций. http : // chaos.phys.msu.ru/loskutov / PDF / Lecturestimeseriesanaly sis.pdf
[4] Анищенко В. С. Динамические системы // Соросовский образовательный журнал, 1997. № 11. С. 77-84.
[5] Кратчфилд Д. П. Хаос //В мире науки, 1987. № 2. С. 18-28.
[6] Зеленый А. С. Об определении порога хаоса одной динамической системы // http : //www.twirpx.com/file/1678632/
[7] Takens F. Detecting strange attractors in turbulence // Lecture Notes in Mathematics, Vol. 898, Springer-Verlag, Berlin, 1980, pp. 366-381; and in Dynamical System in Turbulence, Warlock, 1980, Eds. D. Rand and L. S. Young.
[8] Берже П., Помо И., Видаль К. Порядок в хаосе. О детерминистском подходе к турбулентности / Перевод с франц. Ю. А. Данилова - М.: Мир, 1991. 368 с.
[9] Реферат: Детерминированный хаос. http://www.bestreferat.ru/referat — 95144.html
[10] Ивановский А. Хаос наступает //http : // w — o — s.ru/article/4160
[11] Аттрактор и бабочка //http : / /lenta.ru/articles/2008/04/ 18/lorenz/
[12] Ласкар Ж., Фрошле К., Челлетти А. Измерение хаоса с помощью численного анализа фундаментальных частот. Приложение к стандартному отображению/ Перевод с англ. - А. Г. Арзамасцева - 1992
[13] Wolf A., Swift J. B., Swinney H. L., Vasano J. A. Determining Lyapunov exponents from a time series // Physica 16D. 1985. P. 285-317.
[14] Оселедец В. И., Мультипликативная эргодическая теорема. Характеристические показатели Ляпунова динамических систем // Труды Московского Математического Общества, т.19, 1968. С. 179-210
[15] Мун Ф. Хаотические колебания. / Перевод с англ. Ю. А. Данилова и А. М. Шукурова - М.: Мир, 1990. 312 с
...