Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебники / НЕЛИНЕЙНАЯ ДИНАМИКА СЛОЖНЫХ.pdf
Скачиваний:
640
Добавлен:
30.03.2022
Размер:
28.89 Mб
Скачать

http://profbeckman.narod.ru/

Вторая константа Фейгенбаума отношение между шириной ветви и шириной одной из её подветок (кроме тех, которые ближайшие к изгибу). Не зависит от точной формы f.

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

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

между точками бифуркации имеют одно и то же предельное значение.

20.6 Устойчивость циклов одномерных отображений

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

Прежде всего отметим, что нелинейная система в хаотическом режиме является, тем не менее не полностью непредсказуемой. Белые полосы или «окна» на серой вуали хаотического будущего указывают на локальные состояния порядка с отрицательными показателями Ляпунова; в этом состоит смысл концепции движения к порядку через флуктуации. Положительное значение показателя Ляпунова означает, что орбиты разбегаются экспоненциально быстро, а система чувствительна к начальным условиям. В области rс<r<rf=4 (граница фенитности) показатель Ляпунова возрастает при увеличении r. Поскольку при r<rпоказатель λ<0, то в этом интервале значений хаотичного поведения системы не ожидается. Отрицательные пики при r>rсоответствуют периодическим окнам.

Рис. 22. Зависимость показателя Ляпунова от управляющего параметра =r/4.

Зависимость показателя Ляпунова от r для логистического отображения имеет вид

= lim∑ ln| ′( )|

(19)

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

Рассмотрим поведение графика (r) при большом увеличении.

 

http://profbeckman.narod.ru/

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

 

xn+1=4x(1-xn).

(20)

Здесь 4 =r

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

 

=

 

(large)

 

(21)

 

 

 

 

Экспонента определяется как

 

= lim

ln| ′(

)|

(22)

 

 

 

 

 

 

 

При

>0

соседние

траектории расходятся друг от друга при больших n, что

соответствует хаосу. Однако, если траектории сходятся к фиксированной точке или предельному циклу, они будут сближаться, если <0.

На рис.22 можно легко отличить области хаотического >0 от областей, которые стремятся к фиксированной точке или предельному циклу ( <0). В некоторых точках (первая при =0,75), где показатель Ляпунова достигает 0 и затем снова становится отрицательным. Это бифуркации удвоения периода. Именно в точке удвоения периода система находится на пределе хаоса, но затем, когда период удваивается, становится нехаотичной. Однако в конце периода удвоения режима, при 0,822, пересекает ось и система переходит в хаотический режим. Для в диапазоне значений, превышающих точку, в которой сначала положительно, существует много областей, где показатель отрицателен (острова устойчивости), так что поведение равновесно или является предельным циклом. Рис. 22 наглядно демонстрирует, что при изменении масштаба, за

исключением мелких деталей, не изменяется. Хаос ( >0) возникает при значении

,

лежащем между 0,892486 и 0,892487. Скорость разбегания траекторий

| |=

 

|

|

 

 

.

20.7 Топологическая энтропия

Помимо показателя Ляпунова, для анализа устойчивости отображения и для исследования процессов перехода порядок-беспорядок-порядок используется топологическая энтропия. Её, как меру сложности системы, применяют к характеристики логистического и синус-отображений карте, для описания вариабельности сердечного ритма, к интепрерации электроэнцефалограмм пациентов с эпилепсией и в других областях. Эта мера сложности является информационно-теоретической и называется топологической (перенормированной) энтропией. Исторически она происходит из S- теоремы Климонтовича (S обозначает самоорганизацию). Перенормированная энтропия зависит от изменения управляющего параметра и свидетельствует об относительной степени порядка в системе.

Напомним, что энтропия Шеннона вводится как

S f X f X ln f X dX ,

(23)

где =0 или 1.

Не останавливаясь сейчас на свойствах топологической энтропии, отметим, что она задаётся как разность энтропий состояний 0 и 1.

~

S S f1 X S f0 X , (24)

Перенормированная энтропия выступает как мера сложности. На рис. 23а показано поведение энтропии в области периода 2, где параметр управления r лежит между 3 и 4. В этой области относительная степень порядка возрастает до точки накопления периода, т. е. r(2)c~3.569 поскольку система переходит из состояния равновесия в новое стационарное состояние в соответствии с процессом самоорганизации. Следуя самоорганизации системы в этой области значений r, относительная энтропия монотонно уменьшается. Точка накопления 2-периода движется до тех пор, пока не будет достигнута последнее

http://profbeckman.narod.ru/

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

Рис. 23. Бифуркационная диаграмма (верхняя панель) и ренормированная энтропия (нижняя панель) для 2-периода (а), 3-периода (б) и 5-периода логистического отображения.

Открытая неравновесная динамическая система немонотонно приближается к стационарному состоянию через удвоение периода и отступает от стационарного состояния посредством слияния зон. Соответственно, относительная энтропия уменьшается до точки накопления и увеличивается после её прохода. Стоит отметить внезапные изменения в значениях перенормированной энтропии, которые сигнализируют о существовании окон порядка в хаотической области. На рис. 23в показано поведение энтропии в окне периода 5, которое является одним из автомодельных окон бифуркационной диаграмме логистического отображения. Из-за автомодельности поведение, подобное периоду 3, отражается энтропией: она уменьшается при приближении параметра r к r(5)c~3,743, а затем начинает увеличиваться в соответствии с уменьшением относительной степени порядка. Есть повороты в относительной степени порядка для каждой из трёх точек накопления периода, представляющих стационарное состояние неравновесной динамической системы, обладающей фрактальной структурой.

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

20.8 Синус-отображение

Рассмотрим свойства отображения

 

 

 

 

 

 

 

 

xn+1= sin( xn), 0 1,

0 x 1

(25)

Рис. 24. Синус-отображение.

 

Производная Шварца

 

 

 

 

2

 

3

 

2

 

 

(Sf )(x)

 

1

 

tg

 

( x)) 0

 

 

2

 

 

 

 

 

 

 

 

(26)

предсказывает хаотическое поведение. И действительно,

в этом отображении хаос наступает при =0,98.

При малых есть всего одна неподвижная точка.

Она появляется, если x> sin( x) для всех х или эквивалентно, когда

d

sin x 1, т.е.

dx

 

 

 

 

 

 

 

 

<1/ . Зависимость х( ) представлена на рис. 25. Видно, что при х=0, xn+1=f (xn) xn так

 

 

 

 

 

 

 

http://profbeckman.narod.ru/

что при х=0 неподвижная точка стабильна пока <1/ , маргинальна при =1/ и

неустойчива при >1/ .

 

 

 

 

 

 

 

Рис. 25. График зависимости х( ) от .

Разложение sin(πx) в ряд Тейлора в точке х=0

 

 

 

 

 

и

1

1 . Тогда n+1 (1+ )( n-

даёт x

 

6

 

 

 

 

 

 

 

 

n3), т.е. обычную форму субкритической бифуркации вил, которая всё же отличается от транскритической бифуркации при x=0 для логистического отображения, ибо синусотображение имеет симметрию, что приводит к появлению дополнительной устойчивой

неподвижной точки для отрицательного x. Условие устойчивости для второй

неподвижной

точки

 

 

 

предсказывает существование

области

стабильности, которую

можно найти численно. Условие устойчивости

, где c ≈

|

å cos( å[ ])|< 1

 

 

 

этой точке. Это

 

. I( [0.72]= −1)

 

 

удвоения бифуркации в

0.7200. Поскольку

 

 

 

, то можно найти период

< <

 

 

означает,

что можно получить устойчивый период орбиты 2,

который

затем подвергнет себя периоду удвоения. Можно рассчитать точку (µ, x), при которой n-

кратная итерация отображения ( ) есть устойчивая точка и одновременно ( ) = −1. Неподвижная точка синус-отображения, при которой ненулевой х становится нестабильным: х=0,64574, =0,719962. Точки, которые становятся неустойчивыми при последующих итерациях приведены в табл.1.

Табл. 1. Неподвижные точки синус-отображения, теряющие устойчивость при последовательных итерациях.

Табл. 2. Точки бифуркации, соответствующие локализации устойчивых неподвижных точек

n

Период

n

dn

n

1

2

1,80031

0,2777733

 

2

4

4,04565

-0,107204

-2,59069

3

8

4,55592

0,042518

-2,52139

4

16

4,64516

-0,016962

-2,50668

9

512

4,66920

0,000173

-2,50291

Данные табл. 1 и 2 позволяют вычислить первую и вторую постоянные Фейгенбаумаlim n n 1

 

 

n

n 1

n

(27а)

 

 

 

dn 1

и

n

 

 

 

 

dn .

 

(27б)

 

 

 

 

Здесь dn – наименьшее расстояние между точками в 2n орбитой (уменьшается по

геометрической прогрессии).

. Î

. Í

≈ 4.6605

 

 

Например, для m=3 имеем

n

n

 

. Î

.

 

.

 

 

Как следует из таблицы значения an

и dn различны,

однако параметры

и

стремятся при n к постоянным значениям, причём к тем же, к каким стремится и логистическое отображение.

http://profbeckman.narod.ru/

Рис. 26. Бифуркационная диаграмма синус-отображения: а расчёт , исходный масштаб, [0,1]; б - увеличенный маcштаб.

Бифуркационная диаграмма синус-отображения аналогично по структуре диаграмме f: x x2+c. Она включает переход к хаосу путём удвоения периодов, с окнами устойчивости для нечётных периодов. В квадратичном отображении бифуркационная диаграмма конечна и заканчивается значением параметра, выше которого все орбиты выходят на бесконечность. Однако в синус-отображение продолжается с резкими полосами в хаотическом режиме, раньше переходит в хаос. Эта область перемежается окнами с чётными периодами; наиболее известным из которых является 4-цикл, окна устойчивости шире, чем в логистическом отображении. Хаотический режим расширяется, а затем заканчивается на двух циклах, каждая сторона которых раздваивается. Картина является самоподобной.

Период удвоения бифуркации объясняет орбиты любого периода, кроме периода 3-орбиты. Решим уравнение = ( ) и ′( ) = 1 для

нахождения периода 3-орбиты, которая стабильна только при [0.9378, 09425]. Устойчивый период 3-орбиты также имеет соответствующую неустойчивую орбиту, которая видна на рис. 27.

Рис. 27. Точка, при которой генерируется период 3-орбиты (отмечены устойчивые точки).

Из бифуркационной диаграммы (рис.2.1) видно, что the sine отображение становится хаотичным при приближении r к 1. Для характеристики хаоса используется экспонента Ляпунова, которая для синус-отображения рассчитывается по формуле

= lim

∑ ln| Õcos( )|

(28)

 

 

Зависимость ( ) представлена на рис.28. Снова эта картинка очень похожа на логистическую.

Области, в которых граф находится над пунктирной линией ( =0), являются точками, в которых карта становится хаотичной. Топологическую энтропию кусочно-

монотонного отображения интервала f можно выразить так: h

( f ) lim

1

ln e

, где en

 

top

n n

n

 

 

 

число критических точек отображения fn на отрезке. Топологическая энтропия синус-

отображения при =1 задаётся

h

( f ) lim

1

ln 2n 1 ln 2.

Снова тоже значение, что и

 

 

top

1

n n

 

 

 

 

 

для максимума топологической энтропии для логистического отображения. В более общем виде отображение (x)= sin(x) представляется в виде двух функций +sin( ),

http://profbeckman.narod.ru/

которые графически задаются (+ , ). Бифуркационная диаграмма предполагает, что массив пар кривых C(n)={ n(+x,x)|x R} есть непрерывный ряд всей картины.

Рис. 28. Бифуркационные диаграммы: а xn+1= sin( xn); б xn+1=rxn(1-xn)

Рис.29. Показатели Ляпунова для синус-отображения: а – исходный масштаб, б – увеличенный масштаб.

Рис.30. Полная бифуркационная диаграмма синус-отображения для [-2 ,2 ].

xn+1= sin(xn)+c,

http://profbeckman.narod.ru/

Рис.31. Пары кривых C(n) (n=1, 2, 3, 4), наложенные на бифуркационную диаграмму синус-отображения.

Мы рассмотрели стандартное синусотображение. Коротко остановимся на модифицированном синус-отображении f: xsin x + c.

Рис. 32. Бифуркационная диаграмма f: x sin ( x)+c.

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

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

Имеет место отклонение от поведения, характерного для квадратичного отображения, в широкой хаотической области. Здесь, как и в синус-отображении, хаотическая область имеет окна периодичности. (самые большие окна открыты с одной стороны.) Эта область заканчивается другой квадратичной бифуркационной диаграммой, повернутой на 180° к первой. Поведение диаграммы на интервале [2, 4] идентично поведению на [0, 2] только сдвинутым на две единицы выше. Интервал [0, 2] характерен для остальной части значений параметра и может использоваться как элементарная ячейка. Полная диаграмма проходит по диагонали через начало координат через пространство параметров от отрицательной бесконечности до положительной бесконечности. Ниже приведены некоторые более простые схемы бифуркации. Каждая гладкая и немонотонная функция проходит путь удвоения периода до хаоса и в хаотическом режиме всегда есть окна с нечетной периодичностью.

В более общем виде модифицированное синус-отбражение представляют в виде

(29)

где значения xn ограничены интервалом [1/α, 1-1/α] с 2<α<∞. В Ур.29 максимальная точка возникает при xn=0,5, а ее значение +c, а минимум имеет место при xn=1/α и его значение

sin( / )+c. Таким образом 1-1/ =+c

 

 

1

 

 

 

 

 

 

 

 

 

 

sin

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решения этого уравнения

 

 

 

 

 

2

 

 

 

и c

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 sin

 

 

 

 

 

 

1 sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив в Ур.29, получим

(30)

(31)

Кубическое отображение xn+1=a-bxn+xn3
Положительное значение показателя Ляпунова указывает на то, что динамическая система хаотична. Можно показать, что показатели Ляпунова для модифицированного синус-отображения всегда положительны и равны или превышают значения для стандартного синус-отображения. Это показывает, что модифицированное отображение более хаотично.
На рис. 34 нет белых полос и вся площадь почти вся покрыта точками.
20.8 Хаос в кубических одномерных отображениях

= . [sin( ) − 1]+

http://profbeckman.narod.ru/

(32)

Рис. 33. Показатели Ляпунова: а стандартное, б модифицированное синусотображение.

Рис. 34. Бифуркационные диаграммы: а стандартное, б модифицированное синусотображение.

Показатель Ляпунова – критерий чувствительности к начальным условиям нелинейной динамической системы. Для дискретных во времени систем он имеет вид.

= lim∑ ln| ( )|

(33)

(34)

Кубическое отображение - бимодальное одномерное отображение, представляет собой итерированную дискретную динамическую систему, которая демонстрирует хаотичное поведение.

Отображение

f(x)=mx(1-x2) (35)

имеет две критические точки.

Отображение с двумя критическими точками в симметричном случае может описывать динамику численности некоторых генетических групп.

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

http://profbeckman.narod.ru/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Графики

 

 

 

зависимости

 

 

f(x)

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

различных итерациях приведены на рис. 35.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 35. Зависимость f(x) для m=1,0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,5, 2,0, 3,0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим кубическое отображение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x)=mx(1-x2); f(x)=mx(1+x)(1-x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(36)

 

 

 

 

 

 

 

 

 

Если -1х0 или 1х, то f(x)0, если х -1

или 0х1, то f(x)0. Точки максимума и

минимума 1/

 

и 1/

 

 

; Значения максимума и минимума отображения

2m

 

и

2m

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

соответственно. Когда m=3, f

 

 

 

 

 

 

 

 

 

и f

 

 

 

 

 

 

 

 

 

. В интервале -1х1 и при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

3

3

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

m 0, max f (x)

 

,

т.е. m3.

 

Таким

образом,

 

кубическая функция

 

f(x)=mx(1-x2)

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отображает [-1,1] в

2

 

,

2

, где параметр m заключён в [0,3]. Фиксированные точки

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x)

задаются f(x)=x и есть 0,

 

m 1

 

. Теперь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

имеем

df

m 1 x2 ;

 

df

 

 

 

 

 

 

 

 

 

3 2m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

dx

 

x

m 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 35. Бифуркационная диаграмма

 

 

 

 

 

 

 

 

кубического отображения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неподвижные

точки: х=0 – аттрактор

 

 

 

 

 

 

 

 

(стабилен)

для m<1 и репеллер (нестабилен)

 

 

 

 

 

 

 

 

для

m>0.

Производная

df

 

 

 

 

 

 

 

 

 

 

3 2m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

m 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неподвижные точки x m 1 обе устойчивые точки для 1<m<2 и неустойчивые точки m

для 2<m<3.

На рис. 35 показана бифуркационная диаграмма кубического отображения, из которой следует:

1) В наборе функций {Pni}, n=0, 1, 2,..., (i=1 или 2), где P0i(m)=xci и Pn+1i(m)=f(Pni(m)),

критические значения хс1=1/

 

 

есть

P1

(m)

1

 

 

; P1

m

 

2m

;

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

3

1

 

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4m2

 

 

 

 

3

 

 

4m2

4

 

4

 

 

4m2

 

2

 

 

 

 

 

2m2

1

 

 

 

 

 

 

 

2m

1

 

1

 

 

 

 

 

m

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

 

27

 

 

 

 

 

27

 

 

 

 

 

 

 

1

m

 

 

 

 

27

1

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

; P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

3

3

 

3

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

C

 

C

 

 

 

C

 

 

 

C

 

 

 

 

 

 

 

C

C

(37)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и т.д

http://profbeckman.narod.ru/

Аналогично те же полиномы с отрицательным знаком существуют для

критического значения хс2 1/ 3

Рис. 36: Скелет бифуркационной диаграммы для кубической карты.

(ii) В пределах периодической области (т. е. до значения параметра m=2,3) разного порядка изображения первой критической точки кубического отображения генерируют траектории, которые находятся выше

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

Misiurewicz.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 37. Бифуркация вил: а –

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

суперкритическая,

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

субкритическая.

 

 

 

характера функции f(х) (Ур. 35) при изменении m.

 

 

 

 

 

 

Рассмотрим

изменение

 

 

 

 

 

 

 

 

 

 

 

df

 

 

m(1 3x2 ) ;

 

df

x(1 x2 ),

d 2 f

6mx

;

 

df

 

 

 

m .

Для

 

m=1,

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

dm

 

 

dx2

 

 

 

dx

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

df

 

 

 

1,

df

 

 

 

0,

 

d 2 f

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x 0

dm

 

x 0

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Бифуркация типа вил происходит при х=0 для m=1. Для 0<m1 есть только одна ветвь притягивающуяся к неподвижной точке при х=0. После m=1 эта неподвижная точка теряет свою стабильность, распадаясь на две ветви, притягивающие неподвижные точки

 

m 1

. Как только знак кубического члена

в Ур.1 изменяется на положительный,

 

 

m

 

 

 

 

 

возникает субкритическая бифуркация вил (рис. 37).

 

 

 

 

 

 

 

 

В этом случае неподвижные точки 0 и

 

 

1 m

 

. Для 0<m1 одна ветвь аттрактор

 

 

 

 

 

 

 

m

неподвижной точки х=0 и две новых точки ненулевые неустойчивые точки 1 m m

появляются. После m=1 неподвижная точка х=0 теряет стабильность и другие две неподвижные точки исчезают.

Бифуркационная диаграмма относительно двух критических точек симметрична.

Найдя период удвоения бифуркации точки, рассмотрим одну критическую точку 1/ 3 . Некоторые удвоения бифуркации точек происходят для других критических точек также.

 

df

m 1 3x2 ,

df

x(1 x),

d 2 f

6mx;

df

 

 

 

 

 

3 2m . Если m<2, то

 

 

 

dx2

 

 

 

 

 

 

dx

dm

 

 

 

 

dx

 

x

m 1

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

неподвижные точки x

m 1

 

неустойчивы. Область стабильности неподвижной точки

 

m

http://profbeckman.narod.ru/

x

 

m 1

 

df

 

 

 

 

1, также

df

 

 

 

 

0,

d 2 f

 

 

 

 

 

0 .

 

 

есть 1<m<2. При m=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

m

 

m 1

 

 

dm

 

m 1

 

 

dx2

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

x

m 1

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

m

Бифуркация удвоения периода происходит при m=m1=2, который есть значение первого периода удвоения. Чтобы найти второе бифуркационное значение, рассмотрим

вторую итерацию f2(x) отображения. Неподвижные точки f2(x) задаются f2(x)=x m4x9- 3m4x7+3m4x5-(m4+m2)x3+m2x=0. Его решение

0, ±

, ±

,

±

,

±

 

 

 

 

 

 

 

 

 

 

 

(38)

 

 

 

 

 

 

 

 

 

 

Рис.

38.

 

 

Период

удвоения

 

 

 

 

 

 

 

 

бифуркации при 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из приведённых 9 решений,

 

 

 

 

 

 

 

 

первые 3 – неподвижные точки f,

 

 

 

 

 

 

 

 

которые становятся нестабильными ниже

 

 

 

 

 

 

 

 

m=2, а другие 6 решений – неподвижные

 

 

 

 

 

 

 

 

точки f2(x). Но когда мы определяли

 

 

 

 

 

 

 

 

кубическое

 

 

 

 

 

 

 

 

 

 

 

отображение

 

 

 

 

 

 

 

 

 

f : 1,1

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

, то

пренебрегли

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

неподвижными

точками

x

m 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

m2 4

Новые

неподвижные точки

f

(x)

которые

нас

интересуют x

 

 

 

,

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

2

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х12(1)=

(4) = −

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4) = −

,

 

 

 

 

 

 

 

 

 

(39)

 

 

 

 

 

 

 

 

 

 

 

При значениях параметра меньше m=2,

траектории имеют тенденцию изменяться между двумя значениями x11(1) и х12(1) (верх рис.

39) и х11(2) и х12(2) (низ рис. 39).

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

( )( )

=12

 

( )

 

( )

( )( )

( )

и

( )( )

=

 

( )

 

=

( )( )

( )

(40)

( )

(1)

 

12

 

 

 

 

( )

 

( )

 

 

D

 

D

 

D =

(2)

 

D

 

 

(2)

D

 

D

D

 

D

 

(1)

 

 

и f(x

 

(2)

, и f(x11

(2)

 

 

и f(x12

(2)

 

(2)

 

 

 

 

 

где f(x11 )=x

 

 

=x11

 

 

)=x12

=x11

 

 

 

 

 

 

Производная f(2) одинакова для обеих фиксированных точек x11

(1) и x12

(1) так же как

х11(2) и х12(2) – верхняя и нижняя части двух циклов. Обе эти неподвижные точки или устойчивы или обе неустойчивы, и имеют одинаковую степень стабильности или нестабильности. При дальнейшем увеличении m, производная f(2) уменьшается и

http://profbeckman.narod.ru/

неподвижная точка становится стабильной, причём наклон f(2)<1 для этих фиксированных точек. 2-Цикл неподвижных точек остаётся стабильным до m=m2=2,23606797790250... где производная f(2) при двух циклах неподвижны точек равна -1. Следовательно для значения параметра несколько больше m2, 2-цикла неподвижные точки становятся неустойчивыми. При значении параметра больше m2 траектории располагаются в 4-цикле. Процесс удвоения циклов продолжает продуцировать циклы 23, 24..., пока не будет достигнута предельная точка m=m (аккумулирующая точка удвоения периода). В этой точке периодичность 2 и итерации отображения становятся апериодическими. Этот путь в хаос

через удвоение периода называется сценарием Фельгенбаума. Для расчёта периодических точек f4, необходимо решить 9 одностепенных уравнений типа f4(x)=x. Что даст оставшиеся точки бифуркации:

 

 

 

 

 

Рис. 39. Расстояние между

 

 

 

 

 

успешными

бифуркациями

на

 

 

 

 

 

бифуркационной диаграмме.

 

 

 

 

 

 

Основываясь

на

этих

 

 

 

 

 

бифуркационных

можно

мы

 

 

 

 

 

рассчитать константу Фейгенбаума ,

 

 

 

 

 

которая есть отношение длины между

тремя

успешными точками

бифуркации

cходимости в направлении m-оси.

Тогда

lim

k

4,6692016091029...

Природа

универсальна, одна

и та

же для многих

k

 

 

 

 

 

 

 

отображений. Пусть {mn}-последовательность точек бифуркации. Если m1, m2 известны,

 

 

 

 

 

 

 

m2 m1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

то используя

, то

m

 

 

m

2

, и

m

 

m

 

m

 

 

 

 

 

m , повторяя этот

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

4

 

 

2

1

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аргумент

получим

m

 

mn mn 1

m .Рассмотрим

последовательность {m

 

,n},

 

 

 

 

 

 

 

,n

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn mn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

m ,

где mn – экспериментальное значение

точки бифуркации, Здесь

 

,n

 

 

1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim m ,n

m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

=

 

= 4.542933389479799… ,

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− 3 = 4.64120337856012361… ,

=

 

 

= 4.6631839258902891… ,

 

 

=

 

 

= 4.667909342695353… ,

 

 

 

 

 

 

 

=

= 4.66892502419740… ,

 

 

 

=

 

 

 

 

= 4.66914233815390…

 

 

 

 

 

 

была

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последовательность

аккумулирующих

точек

 

 

рассчитана для некоторых

значений n.

= 2.30040567344817444… ,

 

 

,

= 2.3022834623012582211… ,

 

 

 

 

 

,

 

 

 

 

 

 

 

,

= 2.302193902688090932… ,

 

 

,

= 2.3022834626823954807… ,

 

 

 

 

 

,

= 2.3022793351198851045… ,

 

 

,

= 2.3022834626998726653… ,

 

 

 

 

 

,

= 2.322832728438919278… ,

 

 

,

= 2.3022834627006714427… ,

 

 

 

 

 

,

= 2.3022834539969183689… ,

 

 

,

4 = 2.3022834627007119534…

 

 

 

 

Последовательность сходится к значению 2.3022834627007 – точке

аккумулирования.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Оно

определяется как отношение измеренной

ширины

одной

вилки к таковой к

http://profbeckman.narod.ru/

следующей генерации. Кубическое отображение имеет две критические точки, но для

нахождения достаточно рассмотреть итерации от одной критической точки

 

xc

1/

 

3 .

Пусть

xm – максимум

f, а dn

расстояние от

xm

 

до ближайшей точки в 2n-цикле.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оттношение

dn

стремится к при n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.

 

 

 

 

40.

 

 

 

 

 

Часть

 

 

 

 

бифуркационной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

диаграммы кубического отображения (di

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

расстояние от xm

 

 

до ближайшей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точки 2n-цикла.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Суперстабильный цикл периода 2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

содержит

 

 

одну

критическую

 

 

точку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xc

1/

 

3 в

 

 

 

которой

f

'

 

 

 

 

 

 

 

 

0 .

 

dn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

расстояние между хс и его ближайшим элементом. Элемент,

ближайший к xc

1/

 

 

 

есть

3

 

n-1

 

тая

 

 

 

 

 

 

 

 

 

 

2n 1

 

1

 

 

1

 

,n 1,2,3...Для

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2

 

)

 

 

итерация

x 1/ 3

и

dn f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нахождения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dn f

2n 1

1

сначала находят

суперпритягивающую

 

точку

периода

2-цикла.

 

Если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{х1,х2} – 2-цикл для, то x1=mx2(1-x22) и x2=mx1(1-x12);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1x2=m2x1x2(1-x12)(1-x22) m2(1-x12)(1-x22)=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(41)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если х1 - супперпритягательная точка 2-цикла то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(f2)'(x1)=f'(x1)f'(x2)=0; m2(1-3x12)(1-3x22)=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(42)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 41. Локализация первых четырёх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

изображений

критических

точек

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

построенные

на

 

 

 

бифуркационной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

диаграмме.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

 

или

х1

 

1

или х2

 

 

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

 

 

 

 

 

 

 

х1

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 f (x1) f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

4m2

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81 0 Так

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

1

 

 

1

 

 

 

1 8m

 

 

 

54m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

как f – суперпритягивающая неподвижная точка, когда m=3/2, то (m-3/2) есть фактор его 4-

 

 

 

 

 

 

 

3

3

 

2

 

3

3

3

3

 

ой степени

полинома

m

 

8m

 

12m

 

36m 54 0 8 m

 

m

 

m

 

 

m

 

 

 

0

 

 

 

 

 

 

 

 

 

 

3

 

 

 

2

 

 

 

 

2

2

 

2

 

2

 

Это даёт m

 

, что означает что притягивающая точка для периода 2-цикла находится

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при m 3 .

2

Другие суперпритягивающие неподвижные точки для других периодов 4,8,16,...

приведены ниже.

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

http://profbeckman.narod.ru/

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

Табл.3. Суперпритягивающие неподвижные точки для разных периодов.

Тогда lim dn 2,502...Природа универсальна, одинакова для многих отображений.

n dn 1

Рис. 42. Точки Misiurewicz для каждого периода бифуркации удвоения существует.

Первая Misiurewicz точка

есть точка пересечения

P21(m)=f2(xc1) и P22(m)=f2(xc2),

через которую проходят все последующие высокие итерации. Первую точку Misiurewicz находят решая уравнение P21(m)- P22(m)=0. Между m и первой

точкой Misiurewicz находится большое число отдельных полос которые сливаются вместе при изменении параметра и наконец становятся отдельной полосой. Эта полоса слияния процесса имеет место только когда нестабильная неподвижная точка попадает в аттрактор. Полоса снижает порядок 2n до 2n-1 с n=1,2,3... имет место только когда точки ударяют хаотический или аттрактор или полосу порядка 2n.

Табл. 4. Последовательности {A ,n},

A

,n

 

bn bn 1

b

(b

n

– точки Misiurewicz)

 

 

 

 

1

n

 

 

 

 

 

 

 

 

 

На рис. 42 пунктир означает нестабильные линии, показывающие неустойчивые орбиты, созданные при каждой бифуркации. Справа от точки аккумуляции находится зона, сформированная хаотическими полосами), она существенно отличает от зоны слева от неё, но и имеет много общего, что в левой зоне. В правом экстремуме для m=3 есть только одна зона распространяющаяся на весь интервал от 0 до 1. Это – 1-хаотическая полоса. При уменьшении m полоса сужается. При m=b1 (первая точка Misiurewicz) две полосы раскалываются на четыре пары, что составляет 2 (=21) хаотические полосы. При m=b2 (указана) две полосы раскалываются на четыре пары которые составляют 4(=22) хаотические зоны и т.д. Период удвоения каскада хаотических полос завершается при

http://profbeckman.narod.ru/

точке аккумулирования m =b =3,56994567...Этот регион назван хаотическим регионом С,

а bi – точки Misiurewicz.

Рис. 43. Бифуркационная диаграмма кубического отображения для

2,3 m 3

Интересно отметить, что приведенная выше последовательность также сходится к значению

2.302283462700 ..., т.е. стандартной точке накопления.

Бифуркацию седло-узел часто называют касательной бифуркацией одномерного отображения. Эта бифуркация – одновременное появление устойчивой (узел) и неустойчивой (седло) неподвижных точек, когда параметр проходит через критическое значение, например mt. При изменении m происходит слияние устойчивой и неустойчивой точек, после чего они исчезают. Касательные бифуркации (седло-узел) ответственны за создание периодических окон,

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

Рис. 44. Для m=2,46 f имеет 6 цикл. Чёрные точки s аттракторы 6-цикла,а чёрные точки u указывают нестабильный 6-цикл.

При касательной бифуркации

df ( n) x * 1, где х* – неподвижная точка n-ой dx

итерации отображения, т.е функции fmt(n)(x) – касательная к диагональной линии, т.е. у=х.

Бифуркационная диаграмма демонстрирует вхождение периодического окна внутрь хаотического региона. Из бифуркационной диаграммы рис. 43 видно, что первое наиболее заметное периодическое окно в кубическом отображении есть окно периода 6, которое возникло благодаря касательной бифуркаци. Если отбражение бимодально, то периода 6 окно отражается периода 3 окнами как две полосы. Третья иттерация отображения f3(x) кубического отображения есть ключ к нахождению параметра значения где рождается периода 6 цикл. Любая точка p' в периода 6 цикле удовлетворяет p=f3(p). Так как f3(x) есть 27-кратный полином, который нельзя решить точно для неподвижных

точек, то это делается численно. Из 27 решений уравнения 12 реальны: они отмечены буками s и u на рис. 44: 6 – в верху, 6 – внизу).

Рис.45. Касательная бифуркация в логистическом отображении.

При m=2,44 кривая удаляется по диагонали, поэтому некоторое промежуточное значение между m=2,44 и m=2,46 кривой f3(x) должно стать касательной

к диагонали. При критическом значении m=mt стабильный и нестабильный циклы периода 3 сливаются и аннигилируют в касательной бифуркации. Касательная бифуркация происходит при m=mt=2,450440964487875... Период 6 появляется при m=mt, при котором

http://profbeckman.narod.ru/

f3 прикасается к диагонали в 6-цикле {x1',x2',x3',x4',x5'x6'}, где 3-цикл на верху и 3-цикл внизу.

При mt=2,450440964487875...,

df 3

 

1 и

d2 f

3

 

0 При дальнейшем

dx

dx2

 

 

 

x xi '

 

 

x xi ' .

 

 

 

 

 

 

увеличении m имеет место ряд бифуркаций удвоения периода, что приводит систему в хаотический режим, в котором содержится бесконечное число окон, и бифуркации седлоузел происходят бесконечное число раз.

Показатели Ляпунова играют важную роль в описании поведения динамической системы. Это – метод количественного определения хаотического поведения. В нелинейной системе хаотическое поведение порождает некую случайность и потерю информации, которая объясняет её сложное поведение. Показатель измеряет среднюю скорость расхождения или сходимость орбит, начиная с ближайших исходных точек.

Рассмотрим в xn=f(хn-1) n-итераций от двух близко расположенных начальных точек

х0 и х0+ . Тогда en x0

 

f n x

0

f n x

 

 

, т.е.

x

 

1

ln

 

f n x0

f n x0

 

, где

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

x0 1 f ' x f ' x1 ... f '(xn n

1

 

1

ln

 

f ' x

 

 

ln

 

f ' x

 

 

... ln

 

f ' x

 

 

 

 

 

 

 

 

 

 

 

n

 

 

0

 

 

 

 

1

 

 

 

 

n 1

 

 

(43)

 

 

(

) =

 

 

ln| ′( )|

 

 

 

 

 

 

 

(44)

Рис. 46. Соответствие показателя Ляпунова бифуркацтонной диаграмме.

Переходя к пределу n получим выражение для показателя Ляпунова

 

x0 lim

1

log

df n (x0 )

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

n n

 

 

 

 

(45)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если <0, то

lim f

n x

lim f

n x , а

 

 

n

 

 

0

n

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если

>0,

то

 

 

 

f n x

f n x

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

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

Для кубического отображения период удвоения периода орбит заканчивается в точке в точке аккумулирования m =2,3022834627 – точке, соответствующей орбите периода 2 . Показатель Ляпунова при точке аккумулирования не положителен, что говорит о там, что аттрактор периодический. Для значения m>m экспонента положительна, отражая хаотическое поведение соответствующих аттракторов. Снова Для m<m показатель Ляпунова отрицателен за исключением точек бифуркации (в них =0).

Табл. 5 Значения показателя Ляпунова от параметра m (начиная с точки 1/√3).

http://profbeckman.narod.ru/

Кубическое отображение – простейший представителем бимодальных отображений, демонстрирующих структуры типа

«перекрестка» (crossroad area) и

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

их концах, а также бесконечное множество критических точек коразмерности два (точки МакКая и Ван-Зейтца).

Рис. 47. Зависимость параметра b от параметра а в кубическом отображении.