УДК 372.851
В. С. Секованов, Е. М. Селезнева, В. А. Ивков
Фрактальные методы в физике и экономике на уроках математики
В данной статье рассмотрены приложения фрактальной геометрии в физике и экономике. Разработана методика исследования множества Жюлиа преобразования перенормировки к (х) _ ( х 2 + ч - 11 для ч = 0 которое совпадает с множест-
ч \ 2 х + ч - 2) '
вом нулей Янга-Ли в термодинамическом пределе. Получено представление функции К0 (х) в виде композиции двух линейных взаимно-обратных функций и одной нелинейной функции ( \_ 2 1 Установлена связь между периодическими
р\х I — X + .
точками функции R0(x) и периодическими точками функции /) = x2 +1 Доказано равенство периодов точек х и <р(х) для 0 P\x) x 4 ■
отображений R0(x) и p(x) соответственно. Установлено взаимно-однозначное соответствие между замыканиями множеств неподвижных отталкивающих точек для функций R0(x) ир(х). Приведены алгоритмы построения множества Жюлиа с помощью системы программирования и в среде MathCad. Указаны применения идей фрактальной геометрии в физике (в частности, гидроаэродинамика, гидрология, метеорология, а также теория турбулентности, связанная с именами А. Н Колмогорова и У. Фриша) и экономике (моделирование динамики финансовых рынков, построение фрактальных структур времени изменения котировок акций на бирже). Отмечен вклад А. Н Колмогорова в развитие идей самоподобия и каскадности энергетических структур в контексте теории турбулентности.
Ключевые слова: множество Жюлиа, перенормировка, неподвижные и критические точки, орбита точки, фазовый переход.
V. S. Sekovanov, E. M. Selezniova, V. A. Ivkov
Fractal Methods in Physics and Economics at Mathematics Lessons
This article describes the application of fractal geometry in Physics and Economics. It is developed a technique to study Julia set
Rt (x) = i i!±jl± 12 = 0
of the renormalization transformation ^2x + q 2 > for which q coincides with the set of Yang-Lee zeros in the thermodynamic limit. Representation is obtained as a function R°(x as the compositions of two linear mutually inverse functions and a nonlin-
p(x)= x2 +--.
ear function 4 The connection between the periodic points of the function R0 (x) and periodic points of the function
p(x) = x 2 + -1.
4 It is proved equal periods of points х and ф(х) for reflections R0(x) and p(x) respectively. It is established a one - to-
one correspondence between the closure of the set of fixed points repelling functions R°(x and p^). The algorithms to construct Julia set by the system of software and environment MathCad. We apply the ideas of fractal geometry in Physics (in particular, Nuclear Power, hydrology, meteorology, and the theory of turbulence associated with the names of A.N.Kolmogorov and U.Frisch) and Economics (modeling the dynamics of financial markets, the construction of fractal structures of time change of stock quotes on the stock exchange. The contribution of A.N.Kolmogorov into the development of ideas of self-similarity and energy cascade structures in the context of the theory of turbulence.
Keywords: Julia set, renormalization, a fixed and critical point, the orbit of the phase transition.
До недавнего времени считалось, что фракталы - это красивые картинки и ничего более. Однако время и интенсивные исследования в различных областях науки развеяли данный взгляд на фракталы. Идеи фрактальной геометрии находят широкое применение в различных областях человеческих знаний. Данной науке посвящено
много сайтов, написаны монографии, проводятся представительные научные конференции. Об интенсивности исследований в области фрактальной геометрии свидетельствуют работы [2 - 4], [6 - 15].
В данной работе мы остановимся на приложении фракталов в физике и экономике.
© Секованов В. С., Селезнева Е. М., Ивков В. А., 2014
Деррида, Де Сезе и Ициксон впервые обнаружили тождественность нулей Янга-Ли в термодинамическом пределе с множеством Жюлиа преобразования перенормировки
( х2 + д - 1 |2 Яд (х) = I --- I . То есть фазовая граница
I 2 х + д — 2 у
Янга-Ли совпадает с фрактальным множеством Жюлиа преобразования [6]. В монографии
[6] рассмотрены в смысле перенормировок фазовые границы Янга-Ли, продолженные в комплексную плоскость. Однако подробных пояснений для построения множеств Жюлиа рацио-
22 х + д — 1
' нет, что
Жюлиа для функции Яд (х) =
х2 + д — 1 2 х + д — 2
д = 0).
Мы имеем: Я0 (х)
( х2 — 1 |2 (х +1)2
V 2х — 2 у
4
Используя замену координат р(х) = + -4
р
1 (РШ)) = 4
V х 5} 16 + 8 + 16
— 1 = £_ + х +1 = Я0 (х)
4 2 4 0Ч/
нальной функции Я„(х) = , „ „
V 2 х + д — 2
вызывает многочисленные вопросы у бакалавров, студентов, магистров и аспирантов, изучающих теорию фазовых переходов и алгоритмы построения множеств Жюлиа.
Как известно [2], [12], множество Жюлиа для рациональной функции комплексного переменного / (х), обозначаемое J (/), определяется как замыкание ее отталкивающих периодических точек. Причем, если / - полином, то
J(/) = д{х е С :11ш /(п)(х) = а>[, (то есть J(/)
является границей множества, орбиты каждой точки которого стремятся к бесконечности [2].
Исследуем подробно построение нулей Янга-Ли в термодинамическом пределе (множество
Поясним подробно вывод формулы (*). Найдем такую линейную замену р(х) = а • х + Ь (а Ф 0),
что р— о р о р = Я0,
где р(х) = х2 + с, Я0(х) = — + х + -1.
4 2 4
_1 х Ь
Р (х) =---•
Заметим, что а а Далее име-
ем:
~К (л \\\ (ах+Ь)2 + с Ь ^ Ь2 с Ь
Р (рЧррх)))^-----= ах + 2Ьх+—+---
а а а а а
. Таким образом,
2 Ь2 с Ь х2 х 1
ах + 2Ьх + — +---= — + — + —.
а а а 4 2 4
Сравнивая коэффициенты многочленов
2 Ь2 с Ь х2 х 1
ах + 2Ьх +---1----и--1---+ —, получаем
а а а 4 2 4
простую систему уравнений: 1
а = — 4
Ь = 1
4
Ы+с—Ь=1
а а а 4
при
Из системы уравнений (**) находим:
1 Ь 1 1 х ( л 1 + 1
а = —, Ь = —, с = — . То есть р(х) = — • х + —, а 4 4 4 4 4
р(х) = х2 + 4 и равенство Я0 = р— о р о р (*)
установлено.
Применив метод математической индукции,
покажем, что Я0(п'=р—1 ор(п'ор (***). Для
(1)=Р-1 о р(1),
получим равенствоЯ0 = р 1 о р ор (*), п = 1 имеем: Я0( = р 1 орор - равенство
где р(х) = х2 +1, р~1 (х) = 4х — 1 то есть 4
Я0(х) = р~1 (р(р(х))). Действительно:
1) р(х)=4+4;
2) р(р(х)) = | х +1| +1 = — + х + —;
уV 4 4 у 4 16 8 16
3)
справедливо.
Предположим, что равенство справедливо при
п > 1. То есть Я0(п'= р— о р(п' о р . Далее имеем:
ро Я0(п '= р(п) ор и
р (п+:) о р = р о р"> о р = р о ро Я0п. Из последнего равенства получаем, что
р-1 о р(п+1) о р = р— о р о р о Я0(п' = Я0 о Я0(п' = Я0(п+1}
Лп )
2
2
, что доказывает наше утверждение, то есть Я0(по р(п) ор (***).
Установим связь между периодическими точ-
п ✓ » (х + 1)2
ками функции К0 (х) —-
и периодически-
Найдем точки периода 2 для функции р(х). Для точек периода 2 функции р(х) имеет место
Г 2 112 1
уравнение: I х +— I +—_ х. Проведя соответ-I 4) 4
ми
точками функции р(х)_ х2 +—. Неподвиж-
ной точкой отображения К0 (х) _
_(х +1)2
ствующие 1 ]2
4 ) 4
вычисления
получим
2 112 1 Г 112 Г 2 51 Л
х +—I + — х _1 х —I I х + х + — 1_ 0.
4
явля-
/^1,1 1
ется точка х _ 1. Найдем р(х) _ — 1 +—_ —.
4 4 2
Нетрудно проверить, что точка р(х) _ — является неподвижной точкой функции р(х)_ х2 +—.
2 4у
Таким образом, точками периода 2 для функции р(х) будут три точки
11
2
2
Найдем р :(х2)_ 4 •
1
х3 _ — /. 2
Г 1 1-1 _
— + /
1 2
Далее имеем:
Г
(4/ - 3 + 1)
,2 У
2
+1
(2/ -1)2 + 1
2
Г-4 - 4/ + 1 + 11
_ (1 + 2/)2 _ 4/ - 3.
Доказательство. Пусть точка х имеет период Таким образом, мы получили, что периоды п для функции К0 _ р- о р о р. То-
точек х2 _-+ / и р 1 (х2)_4/-3 равны 2
гдаК0 (п2 (х)_ х.
для функций р(х)_ х2 + 4 и К0 (х)
соответственно. Аналогично можно проверить (студентам полезно провести необходимые выкладки), что периоды точек х2 _- - / и
р 1 (х2 ) _ -4/ - 3 также равны 2 для функций
( \ 2 1 п / 1 (х +1)2
р[х) _ х + 4 и К0 (х) _-—4 соответственно.
Докажем общее утверждение: точка х отображения К0 _ р- о р о р будет иметь период п тогда и только тогда, когда точка р(х) будет иметь период п для отображения р. Далее замечаем, что
Согласно формуле (***) имеем р-1 (р(п) (р(х))) _ х. Из последнего равенства получаем, что р(п) (р(х)) _ р(х). Таким образом, точка р(х). имеет период п для функции р (х). Докажем обратное утверждение. Пусть точка у имеет период п для функции р. То есть
р(п)(У) _ у. Существует такая точка г е X , что
р(г)_ у. Тогда имеем р ^(р^г))_ у. Из последнего равенства находим, что К {п\г)_р1 (р [п)Ш))_р1 (У)_ г. То есть точка р 1 (у) _ г будет иметь период п для функции К0.
2
2
2
2
2
(Ww)' =
■ 4 • р(
( x — +
14
( (
4 р (
V V
1 > 1
4 у 4
x 1
—+44
1
)
=4
p
(x
v 4
Л Л
+
4
—
= 4 p(
))
x+1 Ii x+1
4 4 4 4 у
-р^Ш
Из последнего равенства заключаем, что неподвижная точка x е X функции R 0 (n) будет отталкивающей тогда и только тогда, когда неподвижная точка p(x) е X будет отталкивающей точкой для функции р(n).
Пусть теперь множество A есть множество всех неподвижных отталкивающих точек функции R0(n) n = 1, 2,..., а множество В - множество всех неподвижных отталкивающих точек р(n) n = 1, 2,... . Тогда между множествами A и В функция р устанавливает взаимно однозначное соответствие, то есть данные множества эквивалентны и мы имеем р 1 (в) = A.
Нетрудно проверить, что функция р будет устанавливать взаимно однозначное соответствие и между замыканиями A и В соответственно множеств A и В , и будет иметь место равенство р— (в)= A. Действительно, пусть y е A. Тогда существует такая последовательность {Уп }Г=1 с A = Р- (в)> что lim Уп = У. То-
гдаУп е Р-(в)
Далее мы имеем lim р(Уп )=р(у), где
{р(Уп )}!=1
с B. Следовательно, Р(У) е в
и
y ер (в) Из последнего включения находим, что и A с р-(в). Пусть теперь x е р- (ß) Тогда p(x) е B. Существует такая последовательность ^п }=j с B = р(Л), что lim xn = р{%).
п ^вд
Следовательно, lim р- (xn) = x, где
р1 ({xn }п=i)с A . Таким образом, x е A , и мы имеем р 1 (B )с A. Поскольку выше было установлено включение A ср"1 (b^ то справедливо равенство р— (в)= A. В силу теоремы (см. [2, с. 227]) J(R0 ) = A и J(р) = B где J(R0 ) - мно-
жество Жюлиа для функции Я0, а J(р) - множество Жюлиа для функции р. Поскольку
(р(х) = — + - линейное отображение, то множество J (р ) = В в нашем случае получается из множества J(Я0) = А с помощью гомотетии (см. [5, с. 46]).
Под заполняющим множеством Жюлиа мы будем понимать множество, орбиты каждой точки которого пойманы (см. [2, с. 218]). Множеством Жюлиа будет граница заполняющего множества Жюлиа (Рис. 3). Алгоритм построения заполняющего множества Жюлиа для функции
/(х) = X2 + с хорошо изучен (см. рис.1 и [12, с. 79]). Построив заполняющее множество J (р), мы без труда построим и заполняющее множество Жюлиа J (Я ) (см. рис.1 и рис. 2).
1
Рис.2. Заполняющие множества Жюлиа для функции р(х)_ х2 +1 (белый цвет)
и функции К0 (х) _
(х +1)2
4
(черный цвет в объедине-
Рис. 1. Алгоритм построения заполняющего множества Жюлиа для функции р(х) _ х 2 + С Нетрудно проверить, что множеством Жюлиа для функции /(х)_ х2 будет окружность радиуса единица с центром в начале координат. Следует отметить, что множеством Жюлиа для
функции /(х)_ х2 - 2 является отрезок [- 2;2] (см. [2, с. 246]). При бесконечно многих значениях параметра С множество Жюлиа функции
/(х) _ х2 + С имеет фрактальную структуру.
нии с заполняющим множеством Жюлиа для функции р(х) _ х +—.
Оставив границу заполняющего множества Жюлиа для функции К0 (см. рис.3), мы получим
множество Жюлиа этой функции, что согласовывается с [6, с. 116(с)].
Ч
Рис. 3. Множества Жюлиа для функций
К0 (х) _
(х +1)2 1
и р(х)_ х2 +1
Замечание. Множество Жюлиа для функции р(х)_ х2 + -4 можно построить и с помощью пакета МаШСАБ (рис. 4, см. [12, с. 86]).
4
МО
Чч)
Рис. 4. Множество Жюлиа для функции р(х) = X2 + , построенное в среде МаШСАБ
Укажем применение идей фрактальной геометрии в других разделах физики и экономики.
Мандельброт отмечает, что метод фрактальной геометрии стал частью математического инструментария гидроаэродинамики, гидрологии и метеорологии. По словам Мандельброта, эффективность метода объясняется уникальной способностью выражать большое количество запутанных, неупорядоченных данных несколькими простыми формулами. Эта способность особенно ярко проявляется в случае мультифрактальности - фундаментального понятия при изучении турбулентности и полезного инструмента на финансовых рынках. Он и другие ученые на протяжении последних нескольких десятилетий использовали понятия фрактальной геометрии для изучения и создания моделей рынков.
Броуновское движение - случайное и хаотическое движение частичек пыли, взвешенных в воде - также является одним из аспектов фрактальной геометрии, имеющим практическое использование. Случайное броуновское движение имеет частотную характеристику, которая может быть использована для предсказания явлений, включающих большие количества данных и статистики. Идеи фрактальной геометрии применяются в теории турбулентности. В классической книге «Фрактальная геометрия природы» [3] основателя теории фракталов Бенуа Мандельброта, переведенной на многие языки и выдержавшей
несколько изданий, говорится: «Кроме того, именно в контексте турбулентности теория каскадов и самоподобия достигла своих прогнозист-ских триумфов между 1941 и 1948 гг. Главными действующими лицами здесь были Колмогоров, Обухов, Онсагер и Фон Вайцзекер, однако традиция связывает достижения этого периода только с именем Колмогорова». Ученик А. Н. Колмогорова А. М. Яглом дает краткое описание турбулентности: «Турбулентность -явление, наблюдаемое в громадном большинстве течений жидкости и газов, встречающихся как в природе, так и в технических устройствах или лабораторных установках. Оно заключается в наличии беспорядочных пульсаций (т.е. хаотических изменений в пространстве и во времени) скорости и , давления Р и других гидродинамических характеристик рассматриваемых течений, делающих соответствующие гидродинамические поля и (х, * ) Р (х, *) и др. резко изменчивыми и крайне нерегулярными» [5].
Уриэл Фриш, известный своими работами по фрактальным моделям однородной турбулентности, в книге «Турбулентность. Наследие А. Н. Колмогорова» отмечает: «Теория динамических систем не только оказалась полезной в ряде случаев при исследовании турбулентности, но и сама развивалась под влиянием этих исследований» [14].
В статье «Локальная структура турбулентности в несжимаемой вязкой жидкости при очень больших числах Рейнольдса», опубликованной в 1941 году, Колмогоров пишет: «С энергетической точки зрения процесс турбулентного перемешивания естественно представлять себе так: пульсации первого порядка поглощают энергию осредненного движения и передают ее последовательно пульсации более высоких порядков; энергия же самых мелких пульсаций рассеивается в тепловую благодаря вязкости. В силу хаотического механизма передачи движения от пульсаций низших порядков к пульсациям более высоких порядков естественно допустить, что в
пределах малых по сравнению с I(1) областей пространства мелкие пульсации высших порядков подчинены приближенно пространственно изотропному статистическому режиму» [1].
Примеры использования фракталов в экономических исследованиях наиболее ярко отражены в работах [4], [7]. Большая часть приложений методов фрактальной геометрии в экономической теории относится к финансовым рынкам. Оказывается, что поведение цен и курсов валют на товарных и фондовых рынках имеет самоподобную структуру, т.е. одна и та же фрактальная модель поведения цены на рынке будет работать одинаково вне зависимости от временного масштаба (ежедневный или ежемесячный график изменения цены будет представлять одинаковые структуры). Исследователю рынка (трейдеру) необходимо только знать основные модели движения цены и понимать характер этих моделей. Одной из популярных идей в анализе финансовых рынков является идея о фрактальной структуре времени изменения котировки акций на бирже [4]. Но эта идея требует дальнейшего развития специалистами в области прогнозирования.
Библиографический список
1. Колмогоров А. Н. Избранные труды. Математика и механика [Текст] / А. Н. Колмогоров. -М. : Наука, 1985. - 476 с.
2. Кроновер Р. М. Фракталы и хаос в динамических системах [Текст] / Р. М. Кроновер; пер. с англ. под ред. Т. Э. Крэнкеля. - М. : Постмаркет, 2000. - 352 с.
3. Мандельброт Б. Фрактальная геометрия природы [Текст] / Б. Мандельброт. - М. : Ин-т компьютерных исследований, 2002. - 656 с.
4. Мандельброт, Б. (Не)послушные рынки: фрактальная революция в финансах [Текст] /
Б .Мандельброт, Р. Л. Хадсон ; пер. с англ. - М. : Издательский дом «Вильямс», 2006. - 400 с.
5. Маркушевич А. И. Введение в теорию аналитических функций [Текст] / А. И. Маркушевич, Л. А. Маркушевич. - М. : «Просвещение», 1977. - 320 с.
6. Пайтген Х-О, Рихтер П. Х. Красота фракталов. Образы комплексных динамических систем [Текст]. - М. : Мир, 1993. - 176 с.
7. Петерс Э. Фрактальный анализ финансовых рынков. Применение теории хаоса в инвестициях и экономике. Хаос и порядок на рынках капитала [Текст] / Э. Петерс. - М. : Интернет-трейдинг, 2004. -304 с.
8. Секованов В. С. Методическая система формирования креативности студента университета в процессе обучения фрактальной геометрии [Текст] / В. С. Секованов. - Кострома: КГУ им. Н. А. Некрасова, 2006. - 279 с.
9. Секованов В. С. О множествах Жюлиа некоторых рациональных функций [Текст] / В. С. Секованов. - Вестник КГУ им. Н. А. Некрасова. - Т. 18. - №2. - 2013.
10. Секованов В. С. Формирование креативной личности студента вуза при обучении математике на основе новых информационных технологий [Текст]. -Кострома: КГУ им. Н. А. Некрасова, 2004. - 231 с.
11. Секованов В. С. Элементы теории фрактальных множеств: учебное пособие с грифом УМО для студентов классических университетов специальности «Прикладная математика и информатика» [Текст] / В. С. Секованов. - Кострома: КГУ им. Н. А. Некрасова, 2005. - 135 с.
12. Секованов В. С. Элементы теории фрактальных множеств [Текст]: учебное пособие / В. С. Секованов. - 5-е издание, перераб. и доп. - М. : Книжный дом «ЛИБЕРКОМ», 2013. - 248 с.
13. Falconer K. Fractal Geometry: Mathematical Foundations and Applications [Текст] / K. Falconer. -New York: John Wiley, 1990. - 367 p.
14. Фриш У. Турбулентность. Наследие Колмогорова [Текст] / У. Фриш. - М. : ФАЗИС, 1998. - 360 с.
15. Шредер М. Фракталы, хаос, степенные законы (миниатюры из бесконечного рая) [Текст] / М. Шредер. - Москва; Ижевск: НИЦ «Регулярная и хаотичная динамика», 2001. - 528 с.
Bibliograficheskij spisok
1. Kolmogorov А. N. Izbrannye trudy. Matematika i mekhanika [Tekst] / А. N. Kolmogorov. - M. : Nauka, 1985. - 476 s.
2. Kronover R. M. Fraktaly i khaos v dinamicheskikh sistemakh [Tekst] / R. M. Kronover; per. s angl. pod red. T. EH. Krehnkelya. - M. : Postmarket, 2000. - 352 s.
3. Mandel'brot B. Fraktal'naya geometriya prirody [Tekst] / B. Mandel'brot. - M. : In-t komp'yuternykh issledovanij, 2002. - 656 s.
4. Mandel'brot, B. (Ne)poslushnye rynki: fraktal'naya revolyutsiya v finansakh [Tekst] / B.Mandel'brot, R. L. KHadson ; per. s angl. - M. : Izdatel'skij dom «Vil'yams», 2006. - 400 s.
5. Markushevich А. I. Vvedenie v teoriyu analiticheskikh funktsij [Tekst] / А. I. Markushevich, L. А. Markushevich. - M. : «Prosveshhenie», 1977. - 320 s.
6. Pajtgen KH.-O., Rikhter P. KH. Krasota fraktalov. Obrazy kompleksnykh dinamicheskikh sistem [Tekst]. -M. : Mir, 1993. - 176 s.
7. Peters EH. Fraktal'nyj analiz finansovykh rynkov. Primenenie teorii khaosa v investitsiyakh i ehkonomike. KHaos i poryadok na rynkakh kapitala [Tekst] / EH. Peters. - M. : Internet-trejding, 2004. - 304 s.
8. Sekovanov V. S. Metodicheskaya sistema formiro-vaniya kreativnosti studenta universiteta v protsesse obucheniya fraktal'noj geometrii [Tekst] / V. S. Sekovanov. - Kostroma: KGU im. N. А. Nekrasova, 2006. - 279 s.
9. Sekovanov V. S. O mnozhestvakh ZHyulia neko-torykh ratsional'nykh funktsij [Tekst] / V. S. Sekovanov. -Vestnik KGU im. N. А. Nekrasova. - T. 18. - №2. - 2013.
10. Sekovanov V. S. Formirovanie kreativnoj lichnosti studenta vuza pri obuchenii matematike na osnove novykh informatsionnykh tekhnologij [Tekst]. - Kostroma: KGU im. N. А. Nekrasova, 2004. - 231 s.
11. Sekovanov V. S. EHlementy teorii fraktal'nykh mnozhestv: uchebnoe posobie s grifom UMO dlya stu-dentov klassicheskikh universitetov spetsial'nosti «Prik-ladnaya matematika i informatika» [Tekst] / V. S. Sekovanov. - Kostroma: KGU im. N. А. Nekrasova, 2005. - 135 s.
12. Sekovanov V. S. EHlementy teorii fraktal'nykh mnozhestv [Tekst]: uchebnoe posobie / V. S. Sekovanov.
- 5-e izdanie, pererab. i dop. - M. : Knizhnyj dom «LIBERKOM», 2013. - 248 s.
13. Falconer K. Fractal Geometry: Mathematical Foundations and Applications [Tekst] / K. Falconer. -New York: John Wiley, 1990. - 367 p.
14. Frish U. Turbulentnost'. Nasledie Kolmogorova [Tekst] / U. Frish. - M. : FАZIS, 1998. - 360 s.
15. SHreder M. Fraktaly, khaos, stepennye zakony (miniatyury iz beskonechnogo raya) [Tekst] / M. SHreder.
- Moskva; Izhevsk: NITS «Regulyarnaya i khaotichnaya dinamika», 2001. - 528 s