Научная статья на тему 'Методика визуализации множества Жюлиа с использованием информационных технологий'

Методика визуализации множества Жюлиа с использованием информационных технологий Текст научной статьи по специальности «Математика»

CC BY
337
19
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕЛИНЕЙНОСТЬ / МНОЖЕСТВО ЖЮЛИА / ВИЗУАЛИЗАЦИЯ МАТЕМАТИЧЕСКИХ ОБЪЕКТОВ / НЕПОДВИЖНЫЕ И ПЕРИОДИЧЕСКИЕ ТОЧКИ / ОРБИТА ТОЧКИ / ХАОТИЧНОЕ ОТОБРАЖЕНИЕ / NONLINEARITY / JULIA SET / VISUALIZATION OF MATHEMATICAL OBJECTS / MOTIONLESS AND PERIODIC POINTS / POINT ORBIT / CHAOTIC DISPLAY

Аннотация научной статьи по математике, автор научной работы — Секованов Валерий Сергеевич, Смирнов Евгений Иванович, Бабенко Алена Сергеевна, Селезнева Елена Михайловна, Смирнова Алена Олеговна

В статье рассматривается множество Жюлиа рациональной функции f (z ) = z-(z2+1/2z), включая доказательство ее хаотичности на своем множестве Жюлиа. При этом отображение f обладает существенной зависимостью от начальных условий, транзитивно и множество периодических точек f всюду плотно в L. Согласно определению хаоса по Девани, отображение f хаотично на L. с помощью математических методов доказывается, что это множество есть вещественная ось на комплексной плоскости. Исследуется характер неподвижных точек данной функции и хаотичность ее на своем множестве Жюлиа. Разработан алгоритм построения множества Жюлиа функции f (z ) = z-(z2+1/2z) с использованием информационных технологий. Указаны методические приемы изучения множества Жюлиа данной функции f (z ) = z-(z2+1/2z), что дает возможность обучаемым освоить интеграцию математических методов с ИКТ. Эта методика может позитивно влиять на развитие их мотивации к изучению математики и информатики, способствует развитию креативности и исследовательских компетенций.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Секованов Валерий Сергеевич, Смирнов Евгений Иванович, Бабенко Алена Сергеевна, Селезнева Елена Михайловна, Смирнова Алена Олеговна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A Method of Visualization of Julia Set with Use of Information Technologies

In the article the Julia set of the rational function are considered f (z ) = z-(z2+1/2z), including the proof of its randomness on the Julia set. At the same time the function possesses essential dependence on entry conditions, transitive and the set of periodic points f is dense everywhere in L. According to the definition of chaos according to Devani function f is chaotic on L. By means of mathematical methods it is proved that this set is a material axis on the complex plane. Character of motionless points of this function and its randomness on the Julia set is investigated. The algorithm of creation of the Julia set of the function f (z ) = z-(z2+1/2z) with use of information technologies is developed. Methodical methods of studying of the Julia set of the given function f (z ) = z-(z2+1/2z) are specified that gives the chance to students to master integration of mathematical methods with ICT. This method can positively influence the development of their motivation to study mathematics and informatics, promotes development of creativity and research competences.

Текст научной работы на тему «Методика визуализации множества Жюлиа с использованием информационных технологий»

ТЕОРИЯ И МЕТОДИКА ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

УДК 37

В. С. Секованов, Е. И. Смирнов, А. С. Бабенко, Е. М. Селезнева, А. О. Смирнова, Д. В. Елкин Методика визуализации множества Жюлиа с использованием информационных технологий

7 2 + 1

В статье рассматривается множество Жюлиа рациональной функции 27 , включая доказательство ее

/

хаотичности на своем множестве Жюлиа. При этом отображение обладает существенной зависимостью от начальных

f (z ) = z -

[

ладает с f всюду плотно R L

условий, транзитивно и множество периодических точек всюду плотно в ^ . Согласно определению хаоса по Девани,

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

72 + 1

множестве Жюлиа. Разработан алгоритм построения множества Жюлиа функции 27 с использованием

информационных технологий. Указаны методические приемы изучения множества Жюлиа данной функции

7 2 + 1

f (z ) = z -

юиа. Разработан алгоритм построения множества Жюлиа функции ных технологий. Указаны методические приемы изучения 2 , 1

f(z) = z -

2z , что дает возможность обучаемым освоить интеграцию математических методов с ИКТ. Эта методика может позитивно влиять на развитие их мотивации к изучению математики и информатики, способствует развитию креативности и исследовательских компетенций.

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

_THEORY AND METHODOLOGY OF PROFESSIONAL EDUCATION_

V. S. Sekovanov, E. I. Smirnov, A. S. Babenko, E. M. Selezniova, A. O. Smirnov, D. V. Elkin A Method of Visualization of Julia Set with Use of Information Technologies

f (z ) = z - ^

In the article the Julia set of the rational function are considered 2z , including the proof of its randomness on

f

the Julia set. At the same time the function J possesses essential dependence on entry conditions, transitive and the set of periodic

points^ is dense everywhere in L . According to the definition of chaos according to Devani function^ is chaotic onL . By means of mathematical methods it is proved that this set is a material axis on the complex plane. Character of motionless points of this function and its randomness on the Julia set is investigated. The algorithm of creation of the Julia set of the

z2 +1

function 2z with use of information technologies is developed. Methodical methods of studying of the Julia set of

f (z ) = z -

z2 + 1 2z

the given function 2z are specified that gives the chance to students to master integration of mathematical

© Секованов В. С., Смирнов Е. И., Бабенко А. С., Селезнева Е. М., Смирнова А. О., Елкин Д. В., 2016

methods with ICT. This method can positively influence the development of their motivation to study mathematics and informatics, promotes development of creativity and research competences.

Keywords: nonlinearity, Julia set, visualization of mathematical objects, motionless and periodic points, point orbit, chaotic display.

Множества Жюлиа появились в начале прошлого века в результате итерирования функций комплексной переменной. Исследования показали, что в большинстве случаев множества Жюлиа являются фракталами. В настоящее время интенсивно исследуются фрактальные множества на комплексной плоскости (см., например, [1-9]). Это связано с их использованием при создании математических моделей в физике, экономике и других науках. Полезны комплексные фракталы и для реализации дидактических целей, поскольку математические исследования органически переплетаются с разработкой алгоритмов, реализуемых с помощью современных информационных и коммуникационных технологий, включая параллельное программирование, что дает прекрасную возможность формирования креативности студентов и развития их исследовательских компетенций.

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

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

z2 +1

f (z) = z--, включая доказательство ее хаотичности на своем множестве Жюлиа. Кроме того, в

2 z

работе поэтапно строится алгоритм построения множества Жюлиа с помощью информационных и коммуникационных технологий (ИКТ). Аналогичные подходы изучения множества Жюлиа для

z 2 +1

функции f(z ) = z--нам неизвестны. Мы покажем, что множество Жюлиа для рациональной

2z

z2 ^ 1 z2 — 1 Z 1 í __\

функции f(z) = z--=-=---, (z e С, в e [ü, 2ж\) является вещественной осью OX.

2z 2z 2 2z

Множество Жюлиа для полинома комплексного переменного f(z), обозначаемое J(f), определяется как J(f) = d{z : f(n) (z) ^ да, n ^ да}, где 8 - граница области притяжения бесконечности, а f(n) (z) = f (f(n—1 (z)), n = 1, 2, ....

Покажем, что множество Жюлиа для функции h(z) = z2 есть окружность |z| = 1. На самом деле

h(1)(z) = z2, h(2)(z) = (z2) = z4 = z2 • z2. Так как |z4| = |z2| • |z2| = 1, то |h(2)(z)| = 1. Следовательно,

h(2) (z) находится на единичной окружности радиуса единичной длины с центром в начале координат. Аналогично можно проверить, что точки h(3)(z), z)... также находятся на единичной

окружности. Нетрудно проверить, что последовательность h(n)(z0) = h(h(n—1)(z0)) = z02 , n = 1, 2, ... стремится к да, если |z0| > 1. При |z0| < 1 данная последовательность будет стремиться к 0 так, что в

рассмотренном случае множеством Жюлиа будет окружность единичного радиуса с центром в начале координат.

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

Теорема. Пусть f - полином n-й степени n > 2. Тогда следующие определения множества Жюлиа эквивалентны:

а) множество Жюлиа есть граница области притяжения всех притягивающих неподвижных точек функции f , включая да ;

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

Найдем для функции h(z) = z2 множество Жюлиа, воспользовавшись пунктом б) вышеприведенной теоремы.

Как уже отмечалось, h(n)(z) = z2 . Пусть периодические точки порядка p = 1, 2, 3, ... удовлетворяют уравнению z2 = z . Если z Ф 0, то, сократив на z , получим, что z2-1 = 1 . Следовательно, име-

{2mq |

e2—1 : 0 < q < 2p — 2 > =

í 2rnq . . 2wq ]

= s cos ——- +1 sin ——0 < q < 2 — 2 >. Все эти точки лежат на единичной окружности с цен-

тром в начале координат и распределены на ней равномерно. Поскольку

(h( z)( p) )'

= 2p > 1, то каждая

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

72 + 1 72 — 1 Ъ 1 Рассмотрим функцию / (г) = 7 — -

27 27 2 27

Заметим, что 7 = 1, 72 = —1 - нули функции ^(7) = 72 +1 являются неподвижными притягивающими точками функции / (7) = 7--, поскольку / (7 ) = 7 , / (72 ) = 72 и

27

ч , 472 — 272 — 2 72 + 1

/(7 )=1--4? = -27+- ■

ГХ 1 т

Обозначим вещественную ось через Ь . Заметим, что если 7 = гх е Ь, то /(7) =---е Ь .

2 2гх

Таким образом, точки 7 и /(7) лежат на прямой Ь . Следовательно, если 7 е Ь, то и /(7) е Ь ■

Рассмотрим теперь два отображения, каждое из которых является гомеоморфизмом комплексных

— — . . 7 +1 — 1 . . т + 1)

плоскостей X уу : w = р(7) =- и 7 = р = --—.

7 — 1 w — 1

Пусть I(м>) = w2 . Покажем тогда, что I = р° / ° р 1 (то есть w2 = рр/(р 1 (w)))). Пусть w = w0 .

—1 / ^о+1

Тогда 70 = р (w0 ) = I-. Далее получим

wo -1

z_ = f ( zo) = f (p^K)) =

zl - 1

i)

Л2 Г W + 1 ^

w -1

-1

_o_1 _ У o y

2 • zo 2 (wo + ')'

wo - 1

2

w0 + 2w0 +1 - W0 - 2w0 -1 - w0 + 2w0 -1

—2- - 1

у W0 - 2w0 + 1J = _(w0 - 1) _ =- 2( w0 + Wo

O(wo + ])/ o(wo + ')' (w0 - 1)2 2(wo + 1)i

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(wo -1)2__- 2(wo2 +1), w -1

wo - 1 wo - 1

-(wo2 +1)_ i(wo2 +1)

(wo2 - Ц wo2 - 1

2

.м>1 +1 . ж02 +1 + w02 -1

г-Л— + г —-0—

Далее имеем: (( / (( 1(ж0))) =

= о

-1 = ж02-1 Ж2 + 1 ! ж,2 + 1 - ж,2 + 1 -1 ж2 - 1

"о * жо

Таким образом, ((/(( 1(ж ))) = ж2 для каждого ж е С . Раньше мы установили, что множество

Жюлиа для функции Н(г) = г2 есть окружность единичного радиуса с центром в начале координат комплексной плоскости С . Таким образом, изучение траекторий точки при отображении / сводится к изучению траектории точки при отображении I(ж) = м>2.

Обозначим на расширенной комплексной плоскости через Б - окружность радиуса единица с центром в начале координат.

,-т -1 / ч (ж +1)

Покажем, что отображение г = ( (ж) = --— переводит данную окружность на вещественную

ж -1

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

отображении ( 1 (ж) = (-— . Действительно, пусть ж = г. Тогда ( 1 (г) = --г = 1. Если же

ж -1 г -1

(ж + 1) л (ж + 1) « гт,

ж = 1, то --— = да . И, наконец, если ж = -1, то --— = 0 . Таким образом, при отображении

ж -1 ж-1

, (ж + 1)' ,,

( (ж) = --— трем точкам, лежащим на окружности ж = 1, соответствуют три точки на веще-

ж -1

ственной оси Ь . Ясно, что дробно-линейное преобразование ( 1 (ж) = --— переводит единичную

ж -1

окружность в прямую линию Ь .

Покажем, что функция ( :(ж) = (-— переводит внешность круга Б1 = (ж е Ж : > 1} в верх-

ж -1

нюю полуплоскость р = ( е С : 1тг > 0}, а внутренность круга Б2 = (ж е Ж : |ж| < 1} - в нижнюю полуплоскость Р2 = ( е С : 1т г < 0}.

Согласно свойству дробно-линейного отображения, достаточно найти такие две точки щ еЖ и

— , . . (ж + 1) „ , . . (ж + 1) „ II,

ж еЖ , что ( (ж) = —-— е р, а ( (ж2) = -— е р . Возьмем ж = 3 . Тогда ж > 1 и

ж -1 ж -1

(з +1)'

(-1 (3) = --— = 2г е р . Если же ж2 = ~ , то < 1 и

3 — 1 2

, Л Г1 +1) г -1 (- Г г ==2-1 = -1+2=4-3 г 6 р.

^2) г г г-2 -2 + г 5 5 2

2 2

Предложение 1. Пусть дана функция / (г) =---. Тогда точки г и / (г) лежат в одной и той же

2 2г

полуплоскости относительно прямой Ь .

Доказательство. Докажем сначала, что функция / (г) =---переводит каждую точку, лежащую

2 2г

в верхней (нижней) полуплоскости относительно вещественной оси, в точку, лежащую в верхней (нижней) полуплоскости. Возьмем г = а + Ьг, где Ь > 0 (Ь < 0) . Тогда

I (г) =

а + ¡Ь 1 1 ( а--Ь ) 1 ( а ) 1 .(, Ь

2 2(а + ¡Ь) 2

а + ¡Ь - ■

у а + Ь ) 2 у а + ь ) 2 у

а—;-т ! + _ -

а2 + Ь2) 2

Ь +—--- I. Так как

а2 + Ь2 ]

(* Ь ) п (1 ^ Ь ) Л

Ь + —,-т 1> 0

а2 + Ь2 )

У 2 У а2 + Ь2) у

■)< 01, то / (г)

г) принадлежит верхней (нижней) полуплоскости

в зависимости принадлежности точки г плоскостям р (соответственно р).

Кроме того, функция I (г) отображает комплексную плоскость Z на себя. Действительно, если

г0 е Z, то решая уравнение---= г0, получим, что г2 - 2г0г -1 = 0 . Откуда находим, что

2 2 г

г12 = г0 ± ^г02 +1 . Взяв, например, ^ = г + V2 +1 , мы получим, что I(^ ) = г0 . Таким образом,

каждая точка г0 е Z имеет, по крайней мере, один прообраз.

1 ( \ 2 / ч г + - ч / ч (ж +1)-' Рассмотрим я(ж) = ж , ж = ((г) =-, г = ( (ж) = --—

г - - ж -1

г2 +1 г2 -1 2 1

I (г) = г--=-=---. Тогда согласно лемме 1 (см. [7]) будем иметь

2г 2г 2 2г

((/(и)(( 1 (ж))) = ж2 , поскольку к{"^ж) = ж2".

Предложение 2. Если точка 2 е р с Z. Тогда I) ^ - при " ^ да .

Доказательство. Существует такая точка ж, что |ж| > 1 и г = ( 1 (ж). Поскольку

( 2" +1)

((/1 (ж))) = ж2 , то I(и)(( 1 (ж))= ( 1 (ж2 )= (——„- . Из последнего равенства вытекает, что

о"

ж2 -1

,1т к"±4 = г.

ж2 -1

Предложение 3. Если точка г е р с Z , тогда I("^(г) ^ -- при " ^ да . Доказательство аналогично доказательству предложения 2.

Теорема 1. Прямая Ь является замыканием периодических отталкивающих точек функции

I (г ) = ---.

2 2 г

, / \ 2 г + - _х (ж + 1)г'

Доказательство. Рассмотрим вновь п(ж) = ж , ж = ((г) =-, г = ( (ж) = --—,

г - г' ж -1

г( \ г2 +1 г2 -1 г 1 г + -

I (г) = г--=-=---. Как уже отмечалось, ж = ((г) =- является гомеоморфиз-

2г 2г 2 2г г- -

— — г2 +1 —

мом комплексных плоскостей Z и Ж , а функция ) = г--отображает Z на себя.

Согласно лемме 1 [6] будем иметь (([(и)(( 1 (ж)))= ж2 . Поскольку ^и)(ж) имеет 22 1 непо-

т т А") о2"-1

движных точки, то согласно лемме 2 [7] I ' также имеет 2 неподвижных точки.

Предположим, что точка г е Ь имеет период " для функции I. Существует такая точка ж е Б, что г = ( *(ж). Покажем, что данная точка - отталкивающая для функции I^и)(г) . В силу вышеизложенного имеем: ((/(и)(( 1 (ж)))= ж2 . Тогда I(и)(("1 (ж))=( 1 (ж2 ). То есть

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ап)1 — 1/ \\ № +1)1 ^ у '(р (w)) = —--. Далее находим:

w — 1

(/П'(7)) = (/(п (,,•)„ =

^ 2п ,Л

w + 1

оИ

V w —1 у

= I-

2 —1 (w2n — 1)— 2пУ"—1 (w2n +1)

= I

2 nw2n—1 (— 2)_ .2 п+1 w

(w2n — 1)2

(^ — 1)2 (^2П — 1)2 После преобразований получим:

: (/ п)(р-'(w

>п+1

,.2п—1

w — 1

поскольку точки w и = —1 принад-

лежат

£. Так

у (7 ))| = 1 уп >И

w

как

2п+1 w2 — 1

w е £.

то

0 <

w — 1

< 4.

Следовательно,

> 1. Таким образом, точка 7 - отталкивающая для функ-

ции / п) (7).

Покажем теперь, что множество периодических точек функции /(7) всюду плотно на Ь. Пусть

г е Ь и и - некоторая окрестность точки г . Тогда существует такая точка а е £, что г = р 1 (а). Поскольку множество периодических отталкивающих точек всюду плотно на £ для функции

к(^) = ^^, то существует последовательность {аР }и=1 с £ (рк- период точки арк е £ , к = 1,2,...), сходящаяся к а е £. Согласно лемме 2 (см. [6]) и только что приведенным рассуждениям каждая точка последовательности {гир }и=1 с Ь будет периодическая и отталкивающая для функции

/(7) = 7--. Поскольку р 1 гомеоморфизм £ , то последовательность периодических отталки-

27

вающих точек }и=1 сходится к точке г е Ь, где Ь, гр = р ^(аР). Поэтому существует такой номер п0, что гр е и. Таким образом, множество периодических точек функции /(7) всюду плотно на Ь ■ Покажем, что функция /(7) хаотична на вещественной оси Ь. Установим сначала хаотичность

2

функции к(7~) = 7 на £ . Заметим, что если х = ёв, то )(х) = в2 'в''. Возьмем 5 > 0, е = ^ . Пусть

х е £. Найдем такие п§ е ^ и уп = в

1 в+-

3-2"5

е £, чтобы выполнялось условие - расстояние

Р(х Уп5)<5.

Тогда

Н{п5\х) = ё2п5-в"

^бО

= в

2п5 I ■ в+-

3. 2п5

= ёч

р(к^п5\х),^^"(у^))= 1 >е = 2. Следовательно, ^(7) обладает существенной зависимостью от

начальных условий на £ . Проверим транзитивность отображения к(7) на £ . Пусть и и V - открытые множества в £ . Нетрудно заметить, что существует такое п , что ^п)(и)з £. Следовательно, ^(п)(и) РV непусто и транзитивность к(7) доказана. Всюду плотность периодических точек к(7) на £ отмечалась в начале статьи.

2

2

п

п

г"5 в+п I ■ 1

3

и

и

Пусть теперь 5> 0, а = —, г е Ь . Тогда (р{гх ) принадлежат Б . Обозначим ((г ) = X = и

6

найдем такие " е ^ и уПв = е У 3'2 8), чтобы выполнялось условия р(х, уп )< 8 и р(кМ(х), кМ(уп^ ))= 1 >£ = 2 . Обозначим к("8\х) = ж, й("8)(у ) = ж . Ясно, что точки

ж, ж2 е Б . Тогда /• ('(ж))= (ж1±1>, а ('(ж))= ^2+1> .

ж -1 ж -1

Далее найдем расстояние между точками 1(ж1)) и I(")(c 1(ж)):

2ж - ж]

Поскольку ж, ж е Б ,

I)((-1(ж2 ))-I (")((-1(ж2 )) =

1 -Л

' ж2 + 1 ж + Р

У

ж2 - 1 ж1 - 1 ^

ж2 - 1 ' ж - 1

ж2 -1 |ж1 -1 - 4 и 2|ж - ж > 2, то \I[")(cl (ж2))-IW(c~1(—2 ))= | 2ж1| | ж> 1 и существенная зависимость от начальных условий функции I(г) на Ь установлена. Покажем, что I (г) транзи-тивна на Ь . Пусть и и V - два открытых множества на прямой Ь . Тогда существуют такие два открытых множества Ж и Ж на Б, что выполняются соотношения: и = ( 1(Щ) и V = ( 1(Ж2). Поскольку к(ж) = ж2, то существует такое натуральное число " , что ^и)(ЖБ . Далее, поскольку /"СЖ))=(-1(к(")(Ж1)), то /"(Ж))зЬ. Следовательно, /"(Ж))|^ непусто и транзитивность I (г) установлена. Выше мы показали, что множество периодических точек функции I всюду плотно в Ь .

Таким образом, отображение I обладает существенной зависимостью от начальных условий, транзитивно и множество периодических точек I всюду плотно в Ь . Согласно определению хаоса по Девани (см. [2]) отображение I хаотично на Ь .

г 2 + 1

2. Следующим шагом в изучении множества Жюлиа от функции I(z ) = г--будет построе-

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

Опишем сначала идею алгоритма построения множества Жюлиа для функции

^2 \ ¿2 1 2 1 I (г) = г--=-=---: вычисляем корни уравнения ) = г ^ г2 +1 = 0.

2г 2г 2 2г

Далее исследуем орбиту каждой точки г , принадлежащей комплексной плоскости, которую имитирует монитор компьютера (то есть находим итерацию = I)(г) точки г , где " - натуральное число, стабилизирующее орбиту исходной точки), и оцениваем расстояние от точки г до корней уравнения I (г )= г (точки = - и точки г2 = - -).

Если расстояние от точки до точки 2 = / меньше, чем ее расстояние до точки 22 = —I, то точку 2 окрашиваем в красный цвет. Если расстояние от точки до точки 2 = / больше, чем ее расстояние до точки 22 = —I, то точку 2 окрашиваем в зеленый цвет. Разделяющая красный и зеленый цвета прямая линия окажется вещественной осью, которая будет множеством Жюлиа функции

/ (* ) = 2 — ^

z2 +1

2z

Рис. 1. Инициализация программы

Рис. 2. Масштабирование

Первый этап указан на рисунке 1. На данном этапе устанавливаются высота и ширина экрана, число итераций (iter = 9), масштаб — 2 < Re z < 2, — 2 < Im z < 2 (scall = 2), выделяется массив для хранения корней уравнения и с помощью функции getRoots вычисляются корни уравнения z2 +1 = 0 (кото-

рые будут неподвижными точками функции f (z) = z —

z2 + 1. 2 z

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

На третьем этапе (рис. 3) осуществляется анализ всех точек комплексной плоскости и устанавливается цвет точки исходя из работы функции getColotForPoint() (рис. 4). Именно: рассматривается каждая точка экрана и устанавливается ее цвет, в соответствии с работой функции getColotForPoint(). Функция getColotForPoint() (рис. 4) действует следующим образом: если

расстояние от точки zn до точки z1 = i меньше, чем ее расстояние до точки z2 = —i, то точку z окрашиваем в красный цвет. Если расстояние от точки zn до точки z1 = i больше, чем ее расстояние до точки z2 = —i, то точку z окрашиваем в зеленый цвет. Вещественная ось будет границей,

z2 +1

разделяющей красный и зеленый цвета, и окажется множеством Жюлиа функции f(z) = z--

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2z

(рис. 4).

Отметим еще, что для полной реализации алгоритма используется для вычисления итераций функция: function fx(z:complex, iter: integer) (рис. 5).

Рис. 3. Анализ орбит точек

Рис. 4. Определение цвета точки

Рис. 5. Вычисление итераций функции

f (z ) = z -

z2 + 1 2z

Рис. 6. Нахождение корней уравнения

В полной реализации алгоритма задействована еще функция - function getRoots(iter) (рис. 6). В данной функции используется функция solve (), которая позволяет вычислять корни выше указанного уравнения; реализация функции (входит в библиотеку SymPy на языке программирования Python, аналог функции solve в Mathcad).

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

В заключение отметим, что разработанная методика изучения множеств Жюлиа от функции

z2 +1

/ (^) = 2--дает возможность обучаемым познать интеграцию математических методов с ИКТ,

22

что позитивно влияет на их мотивацию к изучению математики и информатики, развивает креативность и исследовательские компетенции.

Библиографический список

1. Бабенко, А. С., Секованов, В. С. Развитие креативных качеств студентов при изучении метода итераций [Текст] / А. С. Бабенко, В. С. Секованов // Математика в образовании : сб. статей. - Вып. 4 / под ред. И. С. Емельяновой. - Чебоксары : Изд-во Чуваш. ун-та, 2008. - С. 52-54.

2. Божокин, С. В., Паршин, Д. А. Фракталы и мультифракталы [Текст] / С. В. Божокин, Д. А. Паршин. - М. ; Ижевск : НИЦ «Регулярная и хаотическая динамика», 2001.

3. Гринченко, В. Т., Мацыпура, В. Т., Снарский, А. А. Введение в нелинейную динамику: Хаос и фракталы [Текст] / В. Т. Гринченко, В. Т. Мацыпура,

A. А. Снарский. - М. : Изд-во ЛКИ, 2007.

4. Кроновер Ричард М. Фракталы и хаос в динамических системах [Текст] / пер. с англ. ; под ред. Т. Э. Крэнкеля. - М. : ПостМаркет, 2000.

5. Морозов, А. Д. Введение в теорию фракталов [Текст] / А. Д. Морозов. - М. ; Ижевск : Институт компьютерных исследований, 2002.

6. Секованов, В. С. О множествах Жюлиа некоторых рациональных функций [Текст] /

B. С. Секованов // Вестник Костромского государственного университета им. Н. А. Некрасова. - 2012. -Т. 18. - № 2.

7. Секованов, В. С. Элементы теории фрактальных множеств [Текст] : учебное пособие. - 5-е изд. / В. С. Секованов. - М. : Книжный дом «ЛИБРОКОМ», 2013.

8. Смирнов, Е. И. О непрерывности полуаддитивного функционала [Текст] / Е. И. Смирнов // Математические заметки. - 1976. - Т. 19. - № 4.

9. Шредер, М. Фракталы, хаос, степенные законы. Миниатюры из бесконечного рая [Текст] / М. Шредер. - Ижевск : Научно-издательский центр «Регулярная и хаотичная динамика», 2001.

Bibliograficheskij spisok

1. Babenko, A. S., Sekovanov, V. S. Razvitie kreativnyh kachestv studentov pri izuchenii metoda iteracij [Tekst] / A. S. Babenko, V. S. Sekovanov // Matematika v obrazovanii : sb. statej. - Vyp. 4 / pod red. I. S. Emel'janovoj. - Cheboksary : Izd-vo Chuvash. un-ta, 2008. - S. 52-54.

2. Bozhokin, S. V, Parshin, D. A. Fraktaly i mul'tifraktaly [Tekst] / S. V Bozhokin, D. A. Parshin. -M. ; Izhevsk : NIC «Reguljamaja i haoticheskaja dina-mika», 2001.

3. Grinchenko, V. T., Macypura, V. T., Snarskij, A. A. Vvedenie v nelinejnuju dinamiku: Haos i fraktaly [Tekst] / V. T. Grinchenko, V. T. Macypura, A. A. Snarskij. - M. : Izd-vo LKI, 2007.

4. Kronover Richard M. Fraktaly i haos v dinamich-eskih sistemah [Tekst] / per. s angl. ; pod red. T. Je. Kijenkelja. - M. : PostMarket, 2000.

5. Morozov, A. D. Vvedenie v teoriju fraktalov [Tekst] / A. D. Morozov. - M. ; Izhevsk: Institut komp'juternyh issledovanij, 2002.

6. Sekovanov, V. S. O mnozhestvah Zhjulia nekotoryh racional'nyh funkcij [Tekst] / V S. Sekovanov // Vestnik Kostromskogo gosudarstvennogo universiteta im. N. A. Nekrasova. - 2012. - T. 18. - № 2.

7. Sekovanov, V. S. Jelementy teorii fraktal'nyh mnozhestv [Tekst] : uchebnoe posobie. - 5-e izd. / V. S. Sekovanov. - M. : Knizhnyj dom «LIBROKOM», 2013.

8. Smirnov, E. I. O nepreryvnosti poluadditivnogo funkcionala [Tekst] / E. I. Smirnov // Matematicheskie zametki. - 1976. - T. 19. - № 4.

9. Shreder, M. Fraktaly, haos, stepennye zakony. Min-iatjury iz beskonechnogo raja [Tekst] / M. Shreder. -Izhevsk : Nauchno-izdatel'skij centr «Reguljarnaja i ha-otichnaja dinamika», 2001.

i Надоели баннеры? Вы всегда можете отключить рекламу.