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

Изв. вузов «ПНД», т. 17, № 6, 2009 УДК 531.36

СТОХАСТИЧЕСКАЯ ЧУВСТВИТЕЛЬНОСТЬ РАВНОВЕСИЙ И ЦИКЛОВ ОДНОМЕРНЫХ ДИСКРЕТНЫХ ОТОБРАЖЕНИЙ

И.А. Башкирцева, Л.Б. Ряшко, И.Н. Цветков

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

Ключевые слова: Модель Ферхюльста, циклы, стохастическая чувствительность.

Введение

В литературе часто встречается утверждение об аналогии между переходами динамических систем от движений одного типа к движениям другого типа (например, от состояния равновесия к периодическому движению) и известными в статистической физике фазовыми переходами второго рода [1,2]. Идея такой аналогии ведет свое начало от работы Г. Хакена [3] и в дальнейшем развита Ю.Л. Климонто-вичем [4]. Эта аналогия оказалась весьма полезной, так как позволяет использовать методы теории фазовых переходов, например методы скэйлинга и ренормализацион-ной группы [5-7].

Регулярные аттракторы (равновесия и циклы) динамических систем, задаваемых детерминированными дискретными отображениями, являются классическими и наиболее изученными объектами современной теории устойчивости [8]. Прежде всего следует отметить исследования переходов от порядка к хаосу через последовательность бифуркаций удвоения периода в соответствии с законом Фейгенба-ума [9-14]. Анализ устойчивости периодических орбит данных систем и оценки показателей Ляпунова приведены в работах [15-18].

В настоящее время внимание исследователей привлекает изучение реакций таких систем на внешние аддитивные и внутренние параметрические случайные возмущения. Большое число работ посвящено исследованию как переходных процессов, так и установившихся режимов - стохастических аттракторов [19-21]. В работах [22-26] показано влияние внешних возмущений на показатели Ляпунова и на интенсивность сплошного спектра цикла. Анализу скейлинга в цепи удвоения периода циклов при переходе к хаосу для стохастически возмущенных дискретных динамических систем посвящены работы [27-29].

74

Полное вероятностное описание случайных состояний системы требует построения функции плотности распределения. Изменение вероятностного распределения во времени задается соответствующим дискретным функциональным уравнением, решение которого даже в одномерном случае связано с большими сложностями. Здесь, наряду с известными проблемами численного решения, возникает необходимость знания закона распределения действующих случайных возмущений. Такая подробная вероятностная информация, как правило, недоступна. Типичной является ситуация, когда известными являются лишь оценки первых двух моментов. В этих обстоятельствах представляет интерес задача получения асимптотических характеристик первых двух моментов случайных состояний динамической системы для возмущений малой интенсивности. Действительно, для случая малых шумов, когда случайные состояния стохастических аттракторов локализуются вблизи соответствующих аттракторов детерминированных систем, отклонения случайно возмущенных траекторий от невозмущенных детерминированных решений можно описать линейными приближениями. Для систем с непрерывным временем, задаваемых стохастическими дифференциальными уравнениями, такой подход реализован на основе теории квазипотенциала с последующей аппроксимацией с помощью функции стохастической чувствительности [30-33].

Метод функции стохастической чувствительности (ФСЧ) для систем из двух стохастических дифференциальных уравнений впервые был описан в [30], где с помощью этого метода была найдена зона сверхчувствительности циклов брюсселято-ра к случайным возмущениям. Возможности метода ФСЧ для систем Лоренца и Рес-слера в зоне удвоения периода исследованы в [31,32]. В работе [33] представлен общий метод построения ФСЧ для циклов трехмерных систем, использующий технику спектрального разложения. Исследованию возможностей ФСЧ в задачах конструирования систем с циклами, имеющими заданные вероятностные характеристики, и подавления хаоса, посвящена работа [34].

Данная работа посвящена распространению аналогичного подхода на дискретные динамические системы.

В разделе 1 вводится конструкция функции стохастической чувствительности для равновесий и циклов общей одномерной дискретной динамической системы. Здесь для значений ФСЧ получены явные параметрические формулы. В разделе 2 на основе метода ФСЧ проводится анализ стохастически возмущенного квадратичного отображения - модели Ферхюльста. Детально исследована стохастическая чувствительность аттракторов (равновесий и циклов) в зонах структурной устойчивости. Показано, что при переходе от порядка к хаосу в цепи бифуркаций удвоения периода стохастическая чувствительность аттракторов растет как геометрическая прогрессия.

1. Функция стохастической чувствительности

Рассмотрим стохастическую систему

xt+1 = f

где f (x) - достаточно гладкая функция, а ^t - некоррелированный случайный процесс с параметрами

E^t = 0, Е^2 = 1, Щг%к = 0 (t = к),

75

E - знак математического ожидания. Функция a(x) позволяет учесть зависимость интенсивности помех от состояния системы. Пусть xt - решение детерминированной системы (1) при е = 0. Рассмотрим xf - решение (1) с начальным условием x0 = Х0 + е^. При малых е чувствительность детерминированного решения xt к случайным возмущениям определяется величиной

xf - xt vt = lim-------.

s—>o е

Для vt справедлива система

vt+i = atvt + atlt,

где at = df (xt)/dx, at = a(xt), v0 = l. Стохастическая динамика последовательности vt определяется изменением первых моментов mt = Evt, Vt = Ev2. Для mt и Vt справедливы следующие уравнения:

mt+i = atmt,

Vt+i = a2tVt + a2t

с начальным условием m0 = El, V0 = El2.

В данной работе исследуется стохастическая чувствительность стационарного и периодических решений. Пусть в системе (1) при е = 0 имеется устойчивое равновесие xt = x. В этом случае системы (2) - (4) имеют вид

vt+i = avt + alt, mt+i = amt, Vt+i = a2Vt + a2,

где a = df (x)/dx, a = a(x). Неравенство |a| < 1 является необходимым и достаточным условием экспоненциальной устойчивости равновесия x. При этом последовательности mt и Vt стабилизируются следующим образом:

lim mt = 0, lim Vt = w

t—x t—x

a

2

1- a2'

При малом е в системе (1) вокруг устойчивого равновесия x формируется стационарное распределение состояний xf с дисперсией D(xf). Рассмотрим отношение

w = 72 D(xs). е2

Величина еtw является первым приближением D(xf). Связывая в системе (1) значения входа е2 и выхода D(xf) ~ еtw, w играет роль коэффициента стохастической чувствительности равновесия x.

Рассмотрим теперь случай, когда система (1) при е = 0 имеет k-цикл Г - множество точек Г = {xi,..., xk}, связанных соотношениями

f(xi) = xi+i (i = 1,...,k - 1), f(xk) = xi.

76

Считаем, что последовательность Xt определена при всех t с условием периодичности xt+k — Xt. Предполагается, что цикл Г является экспоненциально устойчивым. В этом случае системы (2)-(4) имеют периодические коэффициенты

at+k — at> at+k — at.

Пусть a — a\  a2 ...  ak. Неравенство |a| < 1 является необходимым и достаточным условием экспоненциальной устойчивости цикла. При этом последовательность mt стремится к нулю, а последовательность Vt сходится к периодической последовательности wt

lim (Vt — wt) — 0.

t^-Ж

Элементы последовательности wt могут быть найдены следующим образом: величина wi находится по явной формуле wi — q/(1 — a2), q — ak + a|a‘k_i + ... +

+

wi — a-w-i + a2-i (i — 2,...,k).

Условие |a| < 1 гарантирует у системы (1) при малых шумах существование стационарно распределенного решения. Величины e2wt являются первым приближением дисперсии этого стационарно распределенного решения в окрестности точки xt цикла Г.

Значения w1,... ,wk введенной здесь k-периодической функции wt характеризуют реакцию точек Xi,..., Xk цикла Г на малые случайные воздействия системы (1). Функцию wt будем называть функцией стохастической чувствительности цикла.

Явные формулы для значений функции стохастической чувствительности.

Случай 2-цикла. Для цикла, состоящего из двух точек Xi,X2, имеем

wi

a2 + a2a2 w — a2 + aia2

1 — aja2 1 — a2a2

Случай k-цикла. Пусть цикл Г состоит из точек х-\_,... ,Xk. Тогда для каждой точки Xi функция стохастической чувствительности находится по явной формуле

wi

Jk+i-i + ak+i-2afc+i-i + ... + afrii = i;i=i a

1 — П

2

3

Величины m — mini=i k wi, M — maxi=i k wi характеризуют диапазон изменения стохастической чувствительности состояний k-цикла. Их отношение M/m может служить показателем неравномерности стохастической чувствительности цикла. Величину M удобно использовать в качестве характеристики стохастической чувствительности аттрактора в целом; будем называть ее показателем стохастической чувствительности.

77

2. Система Ферхюльста

Рассмотрим стохастическую систему Ферхюльста

xt+1 = ц—^1 - xt) + £%t,

где §t - независимые случайные величины с параметрами E^t = 0, E^2 = 1, а £ - интенсивность шумов.

При отсутствии случайных возмущений (£ = 0) детерминированная модель (7) имеет на интервале 0 ^ ц ^ 4 различные типы динамики. Зоны порядка чередуются с зонами хаоса (рис. 1). На интервале I = (1, цте), где наблюдаются бифуркации удвоения периода, можно выделить зоны структурной устойчивости с постоянной кратностью циклов, разделенные точками бифуркаций - интервалы Io,h,... На интервале In наблюдаются 2п-циклы с состояниями xn,i(ц) (i = 1,...,2п). Точки аттрактора удобно упорядочить так, чтобы первой была точка цикла —пд, являющаяся ближайшей к 0.5: 1 = argmini=1 _ 2п |—пДц) — 0.5|. Для системы (7) имеем Io = (1,3), I1 « (3, 3.45), I2 « (3.45, 3.54) и т.д. При этом на I0 имеем одну устойчивую точку покоя —o , i (ц) = 1 — 1/ц, на Ii имеем 2-цикл с состояниями

ц + 1 — /

2ц ’ —1,2(Ц) = 2ц .

С ростом параметра ц точки бифуркации накапливаются к значению цте = 3.5699...

Рассмотрим бифуркационную диаграмму системы (7) при £ = 0.001. Под воздействием случайных возмущений траектории покидают аттрактор детерминированной системы (7), формируя вокруг него стохастический аттрактор с соответствующим стационарным вероятностным распределением. В присутствии шума тонкая структура детерминированного аттрактора размывается (рис. 2). Как видно из рисунка, разброс случайных состояний неоднороден. При изменении параметра ц на интервале структурной устойчивости дисперсия заметно меняется, резко возрастая вблизи точек бифуркации. На интервале I1 наблюдается различие в разбросе случайных состояний вокруг точек — 1д(ц),—1,2(ц), соответствующих детерминированному 2-циклу: нижняя ветка размывается сильнее верхней.

— 1,1

Рис. 1. Бифуркационная диаграмма детермини- Рис. 2. Бифуркационная диаграмма в присут-

рованной модели Ферхюльста при отсутствии ствии шума интенсивностью £ = 0.001

возмущений

78

Детальную картину указанных особенностей зависимости разброса траекторий для различных точек аттрактора от параметра ц дает функция стохастической чувствительности. На каждом интервале структурной устойчивости In определен набор функций wn,i(^) (i = 1,..., 2n), где wn,i(^) задает стохастическую чувствительность xn,i i-го состояния 2”'-цикла при ц е In. Значения функции стохастической чувствительности на интервале In находятся по рекуррентным формулам

Wn,i (ц) =

Qn и 2 , an — an, 1 . 1 - an . . an,k, k — 2n, an,i — ц

,2 + a 2 a 2 + 'n,k ' ^n,k^n, k— 1 ' • •  + an,k' ... " an,2,

wn,i+1(ц) — an,iwn,i(M') + 1-

На интервалах Io и Ii для функции стохастической чувствительности соответствующих аттракторов - равновесия жод(ц) и 2-цикла с состояниями ж1;1 (ц), ж1;2(ц) -можно получить явное представление

wo,i

1

1 - (2 - ц)2 ’

. . 2 + 2л/~0 + D ,

W1,1

2 - 2y/D + D

D(2 - D) ’

D — (ц - 3)(ц + 1).

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

Величины wn,i(ц) — (1/e2)Dn,i(ц) характеризуют стохастическую чувствительность аттракторов системы (7) при фиксированной интенсивности е действующих шумов. На рис. 3 сплошными линиями изображены ветви функции стохастической чувствительности wn>i(ц) на интервалах I0,I1,I2,I3 и звездочками отмечены значения wn>i(^, полученные прямым численным моделированием решений системы (7) при шумах различной интенсивности. Как видим, при малых возмущениях теоретические кривые хорошо совпадают с эмпирическими данными. Однако, в силу высокой чувствительности системы вблизи точек бифуркаций, стохастическая чувствительность, полученная с помощью линейного приближения, не совпадает с эмпирическими данными для конечного шума. Эта разница становится более заметной с ростом кратности цикла и увеличением интенсивности возмущений.

В поведении ветвей функции стохастической чувствительности наблюдаются следующие общие закономерности. На интервале структурной устойчивости стохастическая чувствительность ведет себя достаточно регулярно: при приближении к концам интервала (точкам бифуркации) - неограниченно возрастает, выпукла вниз, имеет минимум в центральной части интервала. Стохастическая чувствительность элементов цикла существенно различается. При этом, чем выше кратность цикла, тем больше перепад значений функции стохастической чувствительности (см. рис. 3).

79

Рис. 3. Функции стохастической чувствительности

Рассмотрим функции т и M, задаваемые на каждом отдельном интервале In соотношениями

т

min

i=1,...,2n

wn,i(^') >

M (Ц

max

i=1,...,2n

Wn,i(^)-

Рассмотрим зависимость функций т(ц), M(ц) и показателя неравномерности стохастической чувствительности К = M/m от параметра ц. На каждом интервале In существует значение цп, при котором показатель M(ц) принимает наименьшее значение

Цп = argmin^ M (ц), Mn = min M (ц).

ц€.1и

Рассмотрим величины

тп = т(цп), Кп = —т—г-.

т

Среди всех стохастических 2п-циклов, наблюдаемых на интервале In, цикл, отвечающий значению ц = цп, является наименее чувствительным. Будем называть его стохастическим суперциклом. Следует отметить, что уже в случае 2-цикла стохастический суперцикл отличается от детерминированного суперцикла: ^et = 1 + л/5 ~ ~ 3.236; ^toch ~ 3.185. Наиболее чувствительным в 2п-цикле является состояние

жпд, а наименее чувствительным - следующее за ним xn,2.

80

Из формулы (8) следует, что при всех M Е I выполняется неравенство т ^ 1. На каждом интервале In имеем minMe/n т(ц) =

= 1, причем минимум достигается в точках детерминированного суперцикла. Значения т(ц) в центральной части каждого интервала In мало отличаются от 1. Следовательно, тп ~ 1, Кп ~ ^Д-п . Значения Мп, ^^п для интервалов In структурной устойчивости системы (7) представлены в таблице.

Как видим, в цепи бифуркаций удвоения периода стохастическая чувствительность циклов неограниченно возрастает. Для больших n последовательность Мп ведет себя как геометрическая прогрессия со знаменателем п

П = lim пп = 6.99 ...

При увеличении n растет неравномерность стохастической чувствительности циклов. Если у 2-цикла чувствительность состояний различается в 5 раз, то уже для 32-цикла это различие достигает 104.

В работах [23,27-29] рассматривалось влияние внешнего шума на свойство самоподобия участков бифуркационной диаграммы. Показано, что, если уменьшать интенсивность воздействия в m-f = 6.61903... раз при переходе через точки бифуркации, то для некоторых участков бифуркационной диаграммы наблюдается самоподобие. Используя метод функций стохастической чувствительности, значение m-f можно связать с показателем роста чувствительности цп следующей формулой:

2 2 Mf = па ,

где а = —2.5029... - вторая константа Фейгенбаума.

Введенная конструкция функции стохастической чувствительности позволяет охватить и случай мультипликативных шумов. Для системы Ферхюльста с мультипликативными шумами наблюдается подобный геометрический рост стохастической чувствительности при удвоении периода, при этом предел соответствующей последовательности пп не зависит от вида шума. Похожие свойства были обнаружены для некоторых потоковых систем. Для систем Лоренца и Ресслера авторами была получена [32,33] константа роста чувствительности, близкая к 6.99.

Таблица

n Mn Mn Mn+1 Пп = Mn

0 2.0 1.0 4.7775

1 3.1855 4.7775 6.7042

2 3.4876 32.0296 6.9088

3 3.5522 221.2874 6.9777

4 3.56615 1544.0766 6.9899

5 3.56911 10793.038 6.9928

6 3.56977 75474.078 6.9934

7 3.56990 527822.22 6.9935

Заключение

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

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

81

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

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

Работа выполнена при частичной поддержке грантов РФФИ № 09-01-00026, № 09-08-00048, Федерального агентства по образованию № 2.1.1/2571, ФЦП 02.740.11.02.02.

Приложение

Утверждение. Пусть у системы (1) при е = 0 имеется устойчивый k-цикл Г. Тогда последовательность mt стремится к нулю, а последовательность V сходится к периодической последовательности wt (wt+k = wt), элементы которой могут быть найдены следующим образом: величина wi находится по явной формуле wi = q/(1 — a2), q = ok + a2ko\_ 1 + ... + (ak  ...  a2)2of, а остальные Wi находятся

рекуррентно

Wi = a2-iWi-i + o2-i (i = 2,... ,k).

Доказательство. Последовательные k итераций систем (3)-(4) с учетом периодичности at и ot можно записать в виде

m

V

Неравенство |a| < 1 является необходимым условием экспоненциальной устойчивости цикла Г. С учетом этого из уравнения (11) следует, что для любого k последовательность mik+i стремится к нулю при l ^ ж. Следовательно, в силу системы (3), и вся последовательность mt стремится к нулю при t ^ ж.

Рассмотрим произвольную последовательность Vt, задаваемую системой (4), и разность At = Vt — wt. Из (12) и определения wt следует, что

A

В условиях |a| < 1 из (13) следует, что limt^^ At = 0. Таким образом, утверждение доказано.

82

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

1. Полак Л.С., Михайлов А.С. Самоорганизация в неравновесных физико-химических системах. М.: Наука, 1983.

2. Стенли Г Фазовые переходы и критические явления. М.: Мир, 1973.

3. Haken H. Synergetics - a field beyond irreversible thermodynamics // Lect. Notes in Phys. Berlin: Springer. 1978. Vol. 84. P. 140.

4. Климонтович Ю. Л. Статистическая физика. М.: Наука, 1983.

5. Синай Я. Г. Теория фазовых переходов. М.: Наука, 1980.

6. Вильсон К.Дж. Ренормализационная группа и критические явления // УФН. 1983. Т. 141. Вып. 2. С. 109.

7. Hu B. Intoduction to real-space renormalizatin-group methods in critical and chaotic phenomen // Phys. Rep. 1982. Vol. 91, № 5. P. 233.

8. Elaydi S.N. An introduction to difference equations. Springer, 1999.

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

10. Feigenbaum M.J. Quantitative universality for a class of nonlinear transformations // J. Stat. Phys. 1978. Vol. 19, № 1. P. 25.

11. Feigenbaum M.J. The universal metric properties of nonlinear transformations // J. Stat. Phys. 1979. Vol. 21, № 6. P. 669.

12. Feigenbaum M.J. The transition to aperiodic behavior in turbulent systems // Comm. Math. Phys. 1980. Vol. 77, № 1. P. 65.

13. Hauser P. R., Tsallis C., Curado M. F. Criticallity of the routes to chaos of the 1 — a\x\z map // Phys. Rev. A. 1984. Vol. 30, № 4. P. 2074.

14. Derida B., Gervois A., Pomeau Y. Universal metric properties of bifurcations of endomorphisms // J. Phys. A. 1979. Vol. 12, № 3. P. 269.

15. Пиковский А.С. О стохастических свойствах простейшей модели стохастических автоколебаний // Изв. вузов. Радиофизика. 1980. Т. 23, № 7. С. 883.

16. Huberman B.A., Rudnick J.Scaling behavior of chaotic flows // Phys. Rev. Lett.

1980. Vol. 45, № 3. P. 54.

17. Huberman B.A., Zisook A. B. Power spectra of strange attractors // Phys. Rev. Lett.

1981. Vol. 46, № 10. P. 624.

18. Huberman B.A., Hirisch J. E., Scalapino D. J.Theory of intermittency // Phys. Rev.

A. 1982. Vol.25, № 1. P 519.

19. Анищенко В. С. Стохастические колебания в радиофизических системах. Ч. 1,2. Саратов: Изд-во СГУ, 1986.

20. Неймарк Ю.И. О возникновении стохастичности в динамических системах // Изв. вузов. Радиофизика. 1974. Т. 17, № 4. С. 602.

21. Неймарк Ю.И., Ланда П.С. Стохастические и хаотические колебания. М.: Наука, 1987.

22. Crutchfield J.P., Farmer J., Huberman B.A. Fluctation and simple chaotic dynamics // Phys. Rep. 1982. Vol. 92, № 2. P. 45.

23. Crutchfield J.P., Packard N.H. Symbolic dynamics of noisy chaos // Physica D. 1983. Vol. 7D, № 1-3. P. 201.

24. Gutierrez J., Iglesias A., Rodiguez M.A. Logistic map driven by dichotomous noise // Phys. Rev. E. 1993. Vol. 48. P. 2507.

83

25. Hall P. , Wolf R. C. L. Properties of invariant distributions and Lyapunov exponents for chaotic logistic maps // Journal of the Royal Statistical Society. 1995. Vol. 57. P. 439.

26. Linz S.J., Lucke M. Parametric modulation of instabilities of a nonlinear discrete system // Phys. Rev. A. 1986. Vol. 33. P. 2694.

27. Кузнецов А.П., Капустина Ю.В. Свойство скейлинга при переходе к хаосу в модельных отображениях с шумом // Изв. вузов Прикладная нелинейная динамика. 2000. Т. 8, № 6. С. 78.

28. Кузнецов А.П., Кузнецов С.П, Седова Ю.В. О свойствах скейлинга при воздействии шума в отображении окружности с числом вращения, заданным золотым средним // Изв. вузов. Прикладная нелинейная динамика. 2005. Т. 13, № 5, 6.

С. 56.

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

30. Bashkirtseva I.A., Ryashko L.B. Sensitivity analysis of the stochastically and periodically forced brusselator // Phys. A. 2000. Vol. 278. P. 126.

31. Bashkirtseva I.A., Ryashko L.B. Sensitivity analysis of stochastically forced Lorenz model cycles under period-doubling bifurcations // Dynamic systems and applications. 2002. Vol. 11, № 2. P. 293.

32. Башкирцева И.А., Ряшко Л.Б., Стихин П.В. Стохастическая чувствительность циклов системы Ресслера при переходе к хаосу // Изв. вузов. Прикладная нелинейная динамика. 2003. Т. 11, № 6. С. 32.

33. Bashkirtseva I., Ryashko L. Stochastic sensitivity of 3D-cycles // Mathematics and Computers in Simulation. 2004. Vol. 66. P.55.

34. Bashkirtseva I., Ryashko L. Sensitivity and chaos control for the forced nonlinear oscillations // Chaos, Solitons & Fractals. 2005. Vol. 26. P. 1437.

Поступила в редакцию 4.02.2008 После доработки 6.04.2009

STOCHASTIC SENSITIVITY OF EQUILIBRIUM AND CYCLES FOR 1D DISCRETE MAPS

I.A. Bashkirtseva, L.B. Ryashko, I.N. Tsvetkov

The response problem of equilibrium and cycles for stochastically forced Verhulst population model is considered. Theoretical and empirical approaches are used for stochastically sensitivity analysis. The theoretical approach is based on the firth approximation method and the empirical approach is based on direct numerical simulation. The correspondence between the two approaches for Verhulst population model is demonstrated. The increase of discrete system sensitivity to external noise in the period-doubling bifurcation zone under transition to chaos is shown.

Keywords: Verhulst model, cycles, stochastic sensitivity.

84

Башкирцева Ирина Адольфовна - родилась в 1965 году, окончила Уральский госуниверситет (1987), кандидат физико-математических наук (1997), доцент кафедры математической физики Уральского госуниверситета. Область научных интересов - асимптотический анализ и численное моделирование нелинейных стохастических систем.

620083 Екатеринбург, Ленина, 51 Уральский государственный университет E-mail: irina.bashkirtseva@usu.ru

Ряшко Лев Борисович - родился в 1953 году, окончил Уральский госуниверситет (1975), доктор физико-математических наук, профессор кафедры математической физики Уральского госуниверситета. Область научных интересов -устойчивость, стабилизация и управление стохастическими системами. Автор ряда статей по данной тематике.

620083 Екатеринбург, Ленина, 51 Уральский государственный университет E-mail: lev.ryashko@usu.ru

Цветков Иван Николаевич - родился в 1983 году, окончил Уральский госуниверситет (2006), аспирант кафедры математической физики. Область научных интересов - анализ чувствительности стохастических дискретных систем.

620083 Екатеринбург, Ленина, 51 Уральский государственный университет E-mail: itsvet@e1.ru

85