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

О свойствах скейлинга идентичных связанных логистических отображений с двумя типами связи без шума и под воздействием внешнего шума Текст научной статьи по специальности «Физика»

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

Аннотация научной статьи по физике, автор научной работы — Кузнецов А. П., Кузнецов С. П., Седова Ю. В.

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

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

Похожие темы научных работ по физике , автор научной работы — Кузнецов А. П., Кузнецов С. П., Седова Ю. В.

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

About scaling properties of identical coupled logistic maps with two types of coupling without noise and under influence of external noise

In this paper the influence of noise in system of identical coupled logistic maps with two types of coupling dissipative and inertial is discussed. The corresponding renormalization group analysis is presented. Scaling property in the presence of noise is considered, and necessary illustrations in «numerical experiment style» are given.

Текст научной работы на тему «О свойствах скейлинга идентичных связанных логистических отображений с двумя типами связи без шума и под воздействием внешнего шума»

О СВОЙСТВАХ СКЕЙЛИНГА ИДЕНТИЧНЫХ СВЯЗАННЫХ ЛОГИСТИЧЕСКИХ ОТОБРАЖЕНИЙ С ДВУМЯ ТИПАМИ СВЯЗИ БЕЗ ШУМА И ПОД ВОЗДЕЙСТВИЕМ ВНЕШНЕГО ШУМА

А.П. Кузнецов, С.П. Кузнецов, Ю.В. Седова

В работе обсуждается свойство скейлинга в системе идентичных связанных логистических отображений с двумя типами связи - диссипативным и инерционным. Представлен соответствующий ренормгрупповой анализ. Обсуждается свойство скейлинга в присутствии шума и даны необходимые иллюстрации в стиле «численного эксперимента».

Введение

Связанные системы с удвоениями периода являются одним из популярных и глубоко изученных объектов нелинейной динамики (см., например, монографии [1, 2] и цитированную там литературу). Определенным направлением исследований таких систем является изучение критических явлений в идентичных подсистемах у порога хаоса, то есть явлений, допускающих ренормгрупповое (РГ) описание по аналогии с известным анализом Фейгенбаума. Так, в [3] было предложено применить РГ-анализ к связанным идентичным системам и было введено представление о двух типах связи, инвариантных относительно РГ-преобразования. В работе [4] продолжено обсуждение РГ-анализа и дано описание элементарных объектов на плоскости параметров в терминах бифуркаций. Весьма объемный РГ-анализ представлен в [5], обсуждение устройства плоскости параметров синфазных движений можно найти и в [6, 7], квазипериодические режимы и бифуркация Неймарка - Сакера обсуждались в [8, 9], а вопросы, связанные с мультистабильностью - в [10].

Важнейшим выводом теории, основанной на ренормгрупповом анализе, является представление о скейлинге, то есть о самоподобном устройстве пространства параметров. Отметим, однако, что первые работы (например, [3]) в этом направлении были выполнены в середине 1980-х годов, и иллюстрации скейлинга были даны всего для нескольких бифуркационных точек и нескольких «уровней» удвоений. В настоящее время развитие компьютерной техники позволяет дать более «весомые» иллюстрации скейлинга, которые будут представлены в настоящей работе. Кроме

того, в начале 1980-х годов был выполнен ряд исследований с применением метода ренормгруппы по воздействию шума на классическую одиночную систему с удвоениями периода (см., например, [11-12]). Что касается связанных систем с шумом, то подобное исследование в настоящее время известно для своего рода альтернативной ситуации односторонней связи существенно неидентичных подсистем [13, 14]. Поэтому представляется важным обсудить соответствующие вопросы для случая идентичных связанных подсистем со взаимной связью. Нам кажется полезным в этом плане воспроизвести также простейший вариант РГ-анализа в одиночной системе с шумом, поскольку именно его «технология» и сопоставление с РГ-анализом идентичных связанных систем делает понятным свойство скейлинга в связанных системах с шумом. Последнее будет проверено в численных экспериментах в различных сечениях пространства параметров для разных типов связи между подсистемами.

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

Рассмотрим систему двух связанных отображений вида

Xn+1 = g (Xn) + £ф (Xn, Vn) , Vn+1 = g (Vn) + (Xn, Vu) ,

где g(x) = 1 — 1.5276x2+0.1048x4 + ... - функция, представляющая собой неподвижную точку ренормпреобразования Фейгенбаума-Цвитановича g (x) = ag (g (x/a)); a = —2.5029... - константа Фейгенбаума; e - малый параметр; y(x,y) - функция связи, удовлетворяющая условию y(x, x) = 0. Вблизи диагонали на плоскости (x, у) имеем

u = (x + у)/2, v = (x — у)/2, |v| « 1, ф (u) = lim ф (x,y),

y^x y — x

Un+1 = g (Un) , Vn+1 = g' (Un) Vn + еф (un) Vn.

Двукратная итерация приведенных уравнений дает

Un+2 = g (g (un)),

Vn+2 = g' (g (un)) g' (un) + e [g' (g (un)) ф (un) + g' (un) ф' (g (un))] , и после перенормировки масштаба u — u/a, v — v/a получаем un+2 = ag (g (un/a)), Vn+2 = {g' (g (un/a)) g' (un/a) + +e [g' (g (un/a)) ф (un/a) + g' (un/a) ф' (g (un/a))]} v = {g' (un) + еф1 (un)} v.

Это уравнения такого же вида, как исходные, но с новой функцией связи ф1 (u) = g' (g (u/a)) ф (u/a) + g' (u/a) ф' (g (u/a)).

Рассмотрим теперь задачу на собственные значения

vф (u) = g' (g (u/a)) ф (u/a) + g' (u/a) ф' (g (u/a)).

Численное решение выявляет две существенные собственные функции

ф1 = д' (и)/и + (о^ раН), ы = а = -2.502907... Ф2 = д (и) , У2 = 2.

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

Модель с чисто диссипативной связью выглядит следующим образом [3]:

Здесь X, У - динамические переменные; X - параметр, отвечающий за удвоения; е1 - параметр связи.

Специальным выбором коэффициентов можно также реализовать модель с чисто инерционной связью [3]:

Обозначения те же, е2 - параметр связи.

Кроме того, имеется собственное число линеаризованного уравнения ренорм-группы бр = 4.6692 ..., фигурирующее в теории Фейгенбаума и связанное с отклонением от критической точки по управляющему параметру. Таким образом, перенормированный оператор эволюции за большое число шагов вблизи критической точки представляется в виде

Отсюда следует свойство универсальности: для любых двух слабо связанных систем с удвоениями периода оператор эволюции вблизи диагонали выражается универсальным соотношением и определяется тремя параметрами (С, С1, С2). Универсальность означает, что для получения полной картины явлений достаточно изучить модели (1), (2).

Кроме того, справедливо свойство скейлинга: в точке пространства параметров (Сб-1, С1а_1, С22-1) реализуются динамические режимы такого же типа, как и в точке (С, С1, С2), но с удвоенным временным масштабом, причем начальные условия для динамических переменных получаются пересчетом на фактор ар = -2.5029 ... В терминах систем (1), (2) свойство скейлинга формулируется так: на плоскости (£1,2, X) вблизи критической точки (0, 1.4011552...) при изменении параметра связи £1 в адис = 2 раза, параметра связи е2 в аИн = ар = -2.5029079... раз и изменении параметра X, отвечающего за удвоения, в бр = 4.6692016... раз получается режим динамики, подобный исходному, но с удвоенным временным масштабом.

Хга+1 = 1 — "ХХ'П + е{Х (XI — Уга2) , Уп+1 = 1 — ХУ2 + £1Х У2 — х2).

(1)

Хп+1 = 1 — XXI + Х£2 • 0.088 • (X2 — У2) + £2 (Уп — Хп) , Уп+1 = 1 — ХУ2 + Х£2 • 0.088 • (Уп2 — X2) + £2 (Хп — Уп) .

(2)

и ^ д (и) + Сбк,

V ^ д' (и) V + С1акф1 (и) V + С22кф2 (и) V.

-0.5 £! 0.5

Рис. 1. Скейлинг на карте динамических режимов системы двух связанных логистических отображений с диссипативной связью (1). Цифрами обозначены области существования циклов основных периодов

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

Свойство скейлинга на языке карт режимов иллюстрируют рис. 1 и 2. Цифры на картах обозначают периоды основных циклов. Исходные значения параметров связи £1,2 меняются в пределах от —0.5 до 0.5, а управляющий параметр X естественно изменять в диапазоне от 0 до 2.0. Для обоих случаев связи хорошо видна

1При компьютерном построении карт динамических режимов было сделано от 2000 до 8000 «предварительных» итераций для выхода на аттрактор (в зависимости от рассматриваемого масштаба структуры), далее определялся период установившегося режима по элементам массива, число которых соответствовало выбранному числу цветов на карте.

Рис. 2. Скейлинг на карте динамических режимов системы двух связанных логистических отображений с инерционной связью (2)

самоподобная, повторяющаяся структура на плоскости параметров, которая демонстрирует все большее сходство при приближении к критической точке. Заметим, что исследуемая система, вообще говоря, характеризуется мультистабильностью [10], и плоскость параметров можно мыслить как совокупность «листов» различных карт. Мы выбираем начальные условия для синфазного листа. Свойство скейлинга будет выполняться на каждом из листов (это, в определенной мере, тоже проявление универсальности, присущей системам, допускающим РГ-описание).

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

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

плоскости параметров: близкому к нулю значению Л соответствует белый цвет; отрицательным значениям - оттенок тем более темный, чем больше по абсолютному значению Л; положительному ляпуновскому показателю соответствует черный цвет. Также белым цветом обозначены точки, в которых расходится итерационный процесс. Как оказывается, на таких картах на плоскости параметров можно визуализировать и анализировать образования, подобные наблюдаемым на картах режимов в динамических моделях без шума. С другой стороны, этот метод выявляет более тонкие и мелкие детали устройства плоскости параметров. Наконец, в отличие от метода карт динамических режимов, он может быть использован и для систем с шумом. (В системах с шумом нельзя выявить режимы с определенным периодом, и карты типа рис. 1, 2 построить невозможно.)

Скейлинг на картах ляпуновского показателя иллюстрируют рис. 3, 4 - масштаб по осям координат будет изменяться так же, как и в случае карт динамических режимов. Однако при построении ляпуновской карты на каждом новом уровне удво-

-0.5 £! 0.5

Рис. 3. Скейлинг на карте ляпуновского показателя системы двух связанных логистических отображений с диссипативной связью (1)

-0.5 е2 0.5

Рис. 4. Скейлинг на карте ляпуновского показателя системы двух связанных логистических отображений с инерционной связью (2). В силу того, что константа скейлинга аИн отрицательна, для наблюдения самоподобия при переходе к следующему уровню иерархии увеличенный фрагмент необходимо зеркально отражать слева направо

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

2. Ренормгрупповой анализ в системе Фейгенбаума с шумом

Метод РГ для систем без шума стал классическим и вошел во многие учебники [15, 16]. Для систем с шумом, однако, он известен в гораздо меньшей степени. Более того, известны его очень сложные варианты (например, с использованием метода прореживания и интегралов по траекториям [16]). Поэтому мы кратко изложим основные принципы РГ-анализа в одиночной системе с удвоениями периода с шу-

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

Хп+1 = до (Хп) + уЦо (Хп) Ьп,

где д0 (х) = 1 — \сх2, Хс = 1.401155189092 ...; и0 (х) = 1; у - параметр, характеризующий интенсивность шума; Ьп - последовательность статистически независимых случайных величин с нулевым средним и фиксированным среднеквадратичным значением о.

Применим это отображение дважды. Полагая параметр интенсивности шума у малым, в первом порядке малости по этому параметру получаем

Хп+2 = до (до (Хп)) + У [д0 (до (Хп)) ио (Хп) Ьп + ио (до (Хп)) Ьп+1] .

Фигурирующее в квадратных скобках выражение будем интерпретировать как новую случайную величину, обозначив ее и (х) 'Ь,п. При этом множитель и1(х) подберем так, чтобы величина 'Е,'п обладала таким же среднеквадратичным отклонением о, как и

Поскольку Ьп и Ьп+1 статистически независимы, среднеквадратичные значения обоих членов суммы просто складываются. Отсюда следует, что

и22 (х) = а2 [(до (до (х)))2 и2 (х) + и2 (до (х)) Итак, уравнение для двукратной итерации мы свели к такому же виду, как исходное:

Хп+2 = д1 (Хп) + уи (Хп) Ьп,

но с новыми функциями

д1 (х)= адо (до (х/а)), и1(х) = ^ д (до (х)))2 ио2 (х) + ио2 (до (х))

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

дк+1 (х) = адк (дк(х/а)), ик+1(х) = ^ (д'к (дк(х)))2 и1(х) + и1 (дк(х))

Согласно теории Фейгенбаума, последовательность дк сходится к пределу -функции д, представляющей собой неподвижную точку функционального уравнения Фейгенбаума - Цвитановича

д (х) = ад (д (х/а)),

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

и к (х) = у^л/фЩ,

то приходим к задаче на собственные функции и собственные значения

Ф (x) = а2 (.д' (д (ж)))2 Ф (x) + Ф (д (x)) .

Ее можно решить численно, если известна большим собственным числом отвечает ^ показанной на рис. 5.

Рис. 5. График собственной функции Ф(х)

функция д и константа а. Решение с наи-= 6.619036513... и собственной функции,

Таким образом, в одиночной системе с удвоениями периода с шумом скейлинг тоже возможен, но при этом появляется новый универсальный фактор ^ = 6.619036513... («константа шума» или иначе «константа перенормировки интенсивности шума»), на который следует уменьшать амплитуду шума при переходе с одного уровня удвоений на другой.

3. Скейлинг в идентичных связанных системах с удвоениями периода

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

• о выборе моделей, которые будут достаточно универсальными,

• о величине константы, ответственной за перенормировку шума.

Если обратиться к представленному выше РГ-анализу в идентичных связанных системах без шума и в одиночных системах с шумом, то можно обратить внимание, что методология и того и другого имеет общую идею - а именно, рассматривается линейное возмущение к оператору, отвечающему удвоениям в одиночной системе. Это позволяет выдвинуть гипотезу о том, что эти факторы в определенной мере независимы, а значит, универсальная модель, для которой будет наблюдаться скейлинг, может быть выбрана в форме (1), (2) с добавленным аддитивным образом шумом. При иллюстрациях скейлинга следует для перенормировки амплитуды шума использовать константу ^ = 6.619036513... Проверим эти утверждения.

Введем в отображениях (1) и (2) в каждую подсистему случайные, статистически независимые последовательности in и r\n, моделирующие воздействие флук-туаций, с одинаковой интенсивностью шума у. В результате получим такие системы двух связанных квадратичных отображений с диссипативной

Xn+1 = 1 - X + £i\ (ХП - Y2) + Yin,

(3)

Yn+1 = 1 - Y + ei\ [Y,2 - ХП) + Yn

и инерционном связью

Хп+1 = 1 - X + \£2 • 0.088 • X - У*) + £2 (Уп - Хп) + у1и, Уп+1 = 1 - У + ^£2 • 0.088 • У - Х2) + £2 (Хп - Уп) + УЛп,

в присутствии шума.

Если амплитуда шума мала, и исследуется поведение на больших масштабах времени, то конкретная форма распределения вероятности для и цп будет несущественна, и поведение системы с шумом будет иметь универсальный характер (этот факт подтвержден и для других критических ситуаций, допускающих анализ в терминах метода РГ, см. [17, 18]). Что касается численных исследований, проведенных в данной работе, то и цп определялись как случайные последовательности, равномерно распределенные в интервале [-0.5,0.5].

При помощи компьютерного моделирования было показано, что фактор ^ действительно позволяет получить точные иллюстрации скейлинга в различных сечениях пространства параметров «интенсивность шума - параметр удвоений - параметр связи» на достаточно глубоких уровнях иерархической структуры. Например, для случая инерционной связи на плоскости (£2, А,) хорошо видно, что с приближением к критической точке последовательность рисунков все с большей точностью повторяет свою структуру (рис. 6).

Рис. 6. Иллюстрация скейлинга для системы связанных квадратичных отображений с инерционной связью (4) в присутствии шума с амплитудой е = 0.01. При переходе от одного фрагмента к другому интенсивность шума уменьшалась в раз

В случае диссипативной связи карта ляпуновского показателя на плоскости параметров (е1, е) отображения (3) при фиксированном критическом значении параметра X = Хс = 1.4011552 ... показана на рис. 7. На исходном фрагменте значения параметра связи е1 находятся в диапазоне от 0 до 0.1, а значения интенсивности шума е - в диапазоне от 0 до 0.5. Для наблюдения свойства скейлинга необходимо пересчитывать масштаб по оси е1 в адис раз относительно е1 = 0, масштаб по оси е - в ^ раз относительно е = 0 и пересчитывать в 2 раза значение ляпуновского показателя, изменяя при этом цветовую гамму.

Рис. 8 представляет собой иллюстрацию скейлинга на карте ляпуновского показателя в плоскости (е, X) отображения (3) для значения параметра диссипативной связи е1 = 0.1. Начальные значения для амплитуды шума е меняются в пределах

0.0 £! 0.1

Рис. 7. Скейлинг на карте ляпуновского показателя для системы двух связанных квадратичных отображений с диссипативной связью (3) на плоскости «параметр связи - амплитуда шума». Параметр удвоений X фиксирован: X = Хс = 1.4011552 ...

0.0 е 0.5

Рис. 8. Скейлинг на карте ляпуновского показателя для системы двух связанных квадратичных отображений с диссипативной связью (3) на плоскости «амплитуда шума - параметр удвоений». Параметр связи фиксирован: £1 =0.1

от 0 до 0.5, а значения параметра удвоений X - в пределах от 0 до 2.0. Для реализации скейлинга надо изменять отклонение по параметру е в ^ раз относительно е = 0, отклонение по параметру X - в Ьр раз относительно X = Хс. Также при переходе к новому уровню скейлинга надо уменьшать £1 в адис раз и изменять правило градации цветов.

Оба рисунка, 7 и 8, для случая диссипативной связи с высокой точностью иллюстрируют самоподобие фрагментов на все более мелких масштабах, что подтверждает наше утверждение о том, что фактор ^ ив случае связанных квадратичных отображений отвечает за перенормировку амплитуды шума.

В случае инерционной связи скейлинг на карте ляпуновского показателя на плоскости параметров (е2, е) для отображения (4) при фиксированном критическом значении параметра X = Хс показан на рис. 9. На большом исходном фрагменте

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

0.3 е2 '0.3

Рис. 9. Скейлинг на карте ляпуновского показателя для системы двух связанных логистических отображений с инерционной связью (4) на плоскости «параметр связи - амплитуда шума». Параметр удвоений X фиксирован: X = Хс = 1.4011552.. .

значения параметра связи £2 находятся в диапазоне от —0.3 до 0.3, а значения интенсивности шума £ - в диапазоне от 0 до 0.3. Скейлинг реализуется при пересчете масштаба по оси £2 в аИн раз относительно £2 = 0, а по оси £ - в ^ раз относительно £ = 0. Рис. 10 представляет собой иллюстрацию скейлинга на карте ляпуновского показателя в плоскости (£, X) отображения (4) для значения параметра инерционной связи £2 = 0.01. Начальные значения для амплитуды шума £ взяты в области от 0 до 0.5, а значения параметра удвоений X - в области от 0 до 2.0. Для визуального наблюдения скейлинга надо изменять отклонение по параметру £ в ^ раз относительно £ = 0, а отклонение по параметру X - в Ьр раз относительно X = Хс. Также при переходе к новому уровню скейлинга надо уменьшать £2 в аИн раз.

Рис. 9 и 10 для случая инерционной связи свидетельствуют о большой степени соответствия между получающимися в результате изменения масштаба фрагмента-

-0.5 е 0.5

Рис. 10. Скейлинг на карте ляпуновского показателя для системы двух связанных логистических отображений с инерционной связью (4) на плоскости «амплитуда шума - параметр удвоений». Параметр связи фиксирован: е2 = 0.01

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

Заключение

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

сечениях пространства управляющих параметров и, в определенной мере, не зависит от типа шума. Константой перенормировки интенсивности шума оказывается универсальная постоянная ^ = 6.61903..., вычисленная ранее для одиночной системы.

Работа выполнена при поддержке РФФИ-ННИО (грант № 04-02-04011), Фонда содействия отечественной науке, аналитической ведомственной программы Минобрнауки РФ и CRDF «Развитие научного потенциала высшей школы» (CRDF BRHE REC-006 SR-006-X1/BF5M06 Y3-P-06-07) и гранта Президента РФ (МК-4162.2006.2).

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

1. Анищенко В.С., Вадивасова Т.Е., Астахов В.В. Нелинейная динамика хаотических и стохастических систем // Изд-во Сарат. университета. 1999. 367 с.

2. Mosekilde E., Maistrenko Y., Postnov D. Chaotic synchronization. Applications to living systems // World Scientific Series on Nonlinear Science, Series A. 2002. Vol. 42. 430 p.

3. Кузнецов С.П.Универсальность и подобие в поведении связанных систем Фей-генбаума // Известия вузов. Радиофизика. 1985. Т. 28, № 8. C.991.

4. Kook H, Ling F.H., Schmidt G. Universal behavior of coupled nonlinear systems // Phys. Rev. A. 1991. Vol.43, № 6. P. 2700.

5. Kim S.-Y., KookH. Period doubling in coupled maps // Phys. Rev. E. 1993. Vol. 48, № 2. P. 785.

6. Schult R.L., Creamer D.B., Henyey F.S., Wright J.A. Symmetric and nonsymmetric coupled logistic maps // Phys. Rev. A. 1987. Vol. 35, № 7. P. 3115.

7. Kim S.-Y., Kook H. Critical behavior in coupled nonlinear systems // Phys. Rev. A. 1992. Vol. 46, № 8. P. 4467.

8. Reick C., Mosekilde E. Emergence of quasiperiodicity in symmetrically coupled, identical period-doubling systems // Phys. Rev. E. 1995. Vol. 52, № 2. P. 1418.

9. Rech P. C., Beims M.W., Gallas J.A.C. Neimark-Sacker bifurcations in linearly coupled quadratic maps // arXiv:nlin.CD/0408010. Vol. 1, № 5. Aug 2004.

10. Астахов В.В., Безручко Б.П., Гуляев Ю.В, Селезнев Е.П. Мультистабильные состояния диссипативно связанных фейгенбаумовских систем // Письма в ЖТФ. 1988. Т. 15, № 3. С. 60.

11. Crutchfield J.P., Nauenberg M., Rudnik J.Scaling for external noise at the onset of chaos // Phys. Rev. Lett. 1981. Vol. 46, № 14. P. 933.

12. Shraiman B., Wayne C.E., Martin P.C. Scaling theory for noisy period-doubling transitions to chaos // Phys. Rev. Lett. 1981. Vol. 46, № 14. P. 935.

13. Kapustina J.V., Kuznetsov A.P., Kuznetsov S.P., Mosekilde E. Scaling properties of bicritical dynamics in unidirectionally coupled period-doubling systems in presence of noise // Phys. Rev. E. 2001. Vol. 64, № 6. 066207 (12 pages).

14. Гуляев Ю.В., Капустина Ю.В., Кузнецов А.П., Кузнецов С.П. О свойствах скей-линга при воздействии одной системы с удвоениями периода на другую при наличии шума // Письма в ЖТФ. 2001. Т. 27, Вып. 22. С. 58.

15. Кузнецов С.П. Динамический хаос // М.: Физматлит. 2006. 356 с.

16. Шустер Г. Детерминированный хаос // М.: Мир. 1988. 240 с.

17. Fiel D. Scaling for period-doubling sequences with correlated noise // J. Phys. A: Math. Gen. 1987. Vol.20. P.3209-3217.

18. Choi S.-Y., Lee E.K. Scaling behavior at the onset of chaos in the logistic map driven by colored noise // Phys. Lett. A. 1995. Vol. 205. P. 173.

Саратовский филиал ИРЭ РАН Поступила в редакцию 19.04.2006

Саратовский государственный университет После доработки 15.07.2006

ABOUT SCALING PROPERTIES OF IDENTICAL COUPLED LOGISTIC MAPS WITH TWO TYPES OF COUPLING WITHOUT NOISE AND UNDER INFLUENCE OF EXTERNAL NOISE

In this paper the influence of noise in system of identical coupled logistic maps with two types of coupling - dissipative and inertial - is discussed. The corresponding renormalization group analysis is presented. Scaling property in the presence of noise is considered, and necessary illustrations in «numerical experiment style» are given.

A.P. Kuznetsov, S.P. Kuznetsov, J.V. Sedova

Седова Юлия Викторовна - родилась в 1979 году. Окончила Саратовский государственный университет (2001). Кандидат физико-математических наук (СГУ, 2004), научный сотрудник Саратовского филиала Института радиотехники и электроники РАН. Научные интересы - влияние шума на динамические системы, дискретные отображения, связанные системы, динамический хаос. Автор 10 статей в отечественной и зарубежной печати.

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