Пространственный детерминированный хаос: модель и демонстрация явления в вычислительном эксперименте

Методические заметки

Изв. вузов «ПНД», т. 13, № 1-2, 2005 УДК 535:530.182 + 519.713

ПРОСТРАНСТВЕННЫЙ ДЕТЕРМИНИРОВАННЫЙ ХАОС: МОДЕЛЬ И ДЕМОНСТРАЦИЯ ЯВЛЕНИЯ В ВЫЧИСЛИТЕЛЬНОМ ЭКСПЕРИМЕНТЕ

И.В. Измайлов, А.В. Лячин, Б.Н. Пойзнер, Д.А. Шергин

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

Введение

В рамках изучения явлений нелинейной динамики можно выделить три ветви исследований: временных процессов, в том числе динамического хаоса (если говорить о кольцевых оптических системах, то это, например, [1-6]), пространственно-временных процессов (например, [1, 3, 7-11]) и чисто пространственных (то есть статичных) распределений динамических переменных системы [10], в том числе «пространственного» хаоса. Как известно, задачами, связанными с явлением динамического (временного) хаоса, традиционно занимается радиоэлектроника. Пространственно-временной хаос (турбулентность) проявляется в системах различной природы и рассматривается в рамках соответствующих специализированных отраслей науки.

Явление «пространственного» хаоса характерно для оптических, радиофизических (и, видимо, геофизических [12]) систем. Именно здесь типичными являются случаи, когда поперечный размер нелинейной среды больше длины волны или соизмерим с ней. Тогда аналогами всевозможных процессов, протекающих во времени

в колебательных системах (например, устройствах радиодиапазона), являются пространственные распределения величин каких-то значимых характеристик, скажем, показателя преломления среды, амплитуды, фазы волны. Предложенная аналогия становится более явной, если ввести понятие наблюдателя, который регистрирует эти величины в точках пространства, последовательно перебирая точки по некоторому алгоритму. Системы, в которых возможен пространственный детерминированный хаос (ПДХ), до сих пор мало изучены.

Цель данной работы: 1) показать, что моделью ПДХ может служить дискретное отображение (ДО), получаемое из обыкновенных дифференциальных уравнений (ОДУ) в статическом режиме; 2) в качестве примера получить описание нелинейного кольцевого интерферометра (НКИ) на языке ДО (для случаев двухчастотного оптического поля в приближении больших его потерь и отсутствия дифракции) и продемонстрировать явление пространственного детерминированного хаоса в нелинейном кольцевом интерферометре.

1. Переход от обыкновенных дифференциальных уравнений к дискретным отображениям

Временные процессы в т-мерной динамической системе описывают с помощью системы т эволюционных уравнений

тп ^ = N

где тп - время релаксации, q(t) = {д1 (Ь), д2(Ь),..., дт(Ь)} - вектор динамических переменных. Будем описывать процессы в динамической системе с запаздыванием в многоконтурной обратной связи системой т эволюционных уравнений

Тп ^ =

где тп - диагональная (т х т)-матрица характерных времен релаксации; т - вектор длины в либо матрица размера тхв времен запаздывания. Тогда Q(t,т) = {д^—т^)} либо Q(t,т) = {^(Ь — т^)} - упорядоченное множество из запаздывающих аргументов системы уравнений (2), которое удобно записывать в виде матрицы размера т х в, г е [1; т], ] е [1; в]. Очевидно, что можно принять Q(t, 0) = q(t).

В статическом режиме и при разрешимости уравнения N(q(t), Q(t, т)) = 0 относительно q(t) модель (2) сводится к рекуррентному соотношению, от которого в приближении т = т, У г е [1; в] можно перейти к модели, описывающей состояние системы в моменты времени Ьп = Ьо + пт:

qп = ¥^п-1).

Выражение типа (3) является т-мерным ДО [5]. Но в силу статичности режима qn = qп-l и, следовательно,

q = гы.

Можно ли трактовать выражение типа (4) - с позиции упоминавшегося выше наблюдателя - как ДО? Можно, но как ДО с пространственной эволюционной переменной, у которого размерность ] < т, причем количество итераций ограничено

и не превышает т — 1. А область допустимых значений начальных условий ДО удовлетворяет требованию совместности системы уравнений относительно проекций вектора q, в качестве которой интерпретируется (4).

Поясним это на примере. Пусть все ¥г являются функциями скалярного аргумента и для Уд и Уг верно Fi(q) = ¥(д), величина д1 явно зависит только от дт, а дг - от дг-1, г = 1. Тогда из (4) получим

дг = ¥

Здесь предполагается, что д1 есть начальное состояние ДО (но не динамической системы (2)). Очевидно, если не вводить ограничения на число итераций, то всегда имеет место предельный цикл (так как д1 = ¥(дт)).

При малых значениях т затруднительно квалифицировать (5) как ДО. Если же т достаточно велико, то из-за требования статичности режима динамической системы (2) на первый план выходит непроизвольность д1. Однако в системах высокой размерности (большие т) в ряде случаев область допустимых значений начальных условий ДО практически квазинепрерывна [13], и эта непроизвольность д1 несущественна. А если предположить незамкнутость системы уравнений (4), то ограничения на д1 нет вовсе. Заметим, что возможно т = о.

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

Попытаемся сформулировать определение последнего, но предварительно напомним, что, согласно [14], динамический (детерминированный) хаос - нерегулярное, апериодическое изменение состояния (движение) динамической системы, обладающее основными свойствами случайного процесса.

Тогда можно предложить - в качестве рабочей - следующую дефиницию: пространственный (статический) детерминированный хаос - статическое состояние динамической системы, такое что последовательность подмножеств Рг множества динамических переменных Р, выделяемая наблюдателем по некоторому алгоритму, является нерегулярной, апериодической, обладающей основными свойствами случайного процесса. Причем подмножества Рг равномощны, не совпадают (Рг = Р^ при г = ]), покрывают множество динамических переменных Р (Р = игРг), и пересечения Рг П Рг+1 равномощны.

Эта последовательность может иметь конечную длину (конечная последовательность). Указанный алгоритм (обусловливающий вид, мощность и конкретное «наполнение» подмножеств, порядок их следования) эвристически обусловлен природой динамической системы и структурой взаимной зависимости значений динамических переменных (например, вида (5)). Если алгоритм выбрать затруднительно,

то в зависимости от мощности Р можно предложить два варианта алгоритма «по умолчанию».

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

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

Заметим, что переход от (2) к (3) возможен не только в случае статического режима динамической системы. Он осуществляется и в известном - благодаря работам К. Икеды - приближении мгновенного отклика нелинейной среды, справедливом при |тп| ^ |т|. Но в отличие от случая статического режима перейти к выражению (4) нельзя.

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

2. Модель динамики нелинейного фазового набега поля в кольцевом интерферометре

Развитые выше положения удобно иллюстрировать на примере НКИ, который исследовался теоретически и экспериментально в работах [1-3, 6, 15].

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

Основным отличием рассматриваемого НКИ от модели К. Икеды [1] является то, что в цепь обратной связи интерферометра введен линейный элемент О, обеспечивающий поточечное преобразование поля (сдвиг, наклон, поворот, растяжение, сжатие). Таким образом, траектория луча, начинающаяся в точке с координатами (х', у') после обхода резонатора попадет в точку (х, у) = О(х', у').

Пусть на вход НКИ поступает сумма двух квазимонохроматических полей с амплитудами а(г^), Ь(г,Ь) и с частотами ю ± круговой поляризации различных

Рис. 1. Схема нелинейного кольцевого интерферометра и ход лучей в НКИ при повороте светового поля на А = 120о в плоскости х0у: а - траектории лучей 1, 2, 3, замыкающиеся после трех обходов; б - проекция траекторий лучей 1, 2, 3 на плоскость х0у

у

Ф

2а

I

У

ф

(при ш > О) либо одинаковых (при ш < О) направлений вращения (рис. 2):

ех (г, г) = а(г, г) соб[(ш + О)г + ф(г, г) + ^(г, г)] + + ь(г, г) сов[(ш - о)г + ф(г, г) - ^(г, г)],

еу(г, г) = а(г, г) ят[(ш + о)г + ф(г, г) + ^(г, г)] -- ь(г,г) в1п[(ю - о)г + ф(г, г) - ^(г,г)].

Здесь ш (либо О при ш < О) имеет смысл средней частоты, а 2О (2ш при ш < О) - частотный интервал между составляющими поля. Чтобы отразить специфику спектра рассматриваемого оптического поля, мы оперируем параметром бихроматичности д = О/ш [16].

Тогда в приближении медленно меняющихся амплитуд, фаз, модуляции положения плоскости поляризации, времени запаздывания ге поля, медленной релаксации нелинейной части показателя преломления, а также в приближении больших потерь (либо в приближении одного прохода) и в пренебрежении дифракцией поля можно получить описание динамики нелинейного фазового набега и в НС НКИ [16]:

Тп (г) = Пе (г)Аи (г, г) - и (г,г) + / (г,г),

/(г, г) = КаЬ(г, г, гп) + рКаЬ(г', г - т, г) + [у(г', г)/а] х х{Ка(г,г, г',г - т) сов[(1 + д)шт + Ф(г,г) - ф(г',г - т) + ^(г,г) - ^(г',г - т)]+ +кь(г, г, г', г - т) соб[(1 - д)шт + ф(г, г) - ф(г', г - т) - ^(г, г) + <ф(г', г - т)]}.

Здесь к = 2пД - волновое число; и (г, г) - нелинейный фазовый набег в НС на частоте ш; т = т(г' ,г) = ге (г' ,г) + и (г' ,г - ге (г' ,г))/ш; ге (г, г) = п0 (г)1/с + г0 (г, г) -

со-О

О)

со+О

# =£2/со

Рис. 2. Спектр бихроматического излучения

эквивалентное время запаздывания (величина времени распространения светового поля, обусловленного наличием линейной части По(г) показателя преломления НС и наличием времени запаздывания Ь0(г, Ь) в НКИ (I - длина НС, с - скорость света)); у(г', Ь) - удвоенный амплитудный коэффициент потерь излучения за один проход через НКИ; о - коэффициент растяжения пучка в элементе О; Бе(г) - коэффициент диффузии молекул НС, нормированный на время релаксации тп(г); р = 0 для приближения больших потерь, но р = [у(г', Ь)/о/2]2 для приближения одного прохода; «смешанный» (КаЬ) и «парциальные» (К а, КЬ) параметры нелинейности

КаЬ(г,Ь, гп) = (1 — Е)п2(гп)1к[а2(г,Ь) + Ь2(г,Ь)], Ка(г, Ь, г', Ь — т) = (1 — К)п2(г)1к а(г, Ь)а(г', Ь — т), КЬ(г, Ь, г', Ь — т) = (1 — К)п2(г)1к Ь(г, Щг', Ь — т);

п2 (г) - параметр нелинейной рефракции НС.

В пренебрежении диффузией молекул НС из этой модели можно получить описание динамики нелинейного фазового набега и в НС НКИ в «точечном приближении»

т (Шт=тт (Ь) + г

^ — иг\Ь) * Лу ¡г = ¡г(Ь) = КаЬф) + рКаЬг-ц(Ь — т) + [уг-1(Ь)/о] х х{КаМ-1(М — т) ео8[(1 + д)ют + фг (Ь) — фг-1(Ь — т) + — — т)] +

+КЬг,г-1(Ь,Ь — т) еов[(1 — д)ют + фг(Ь) — фг-1(Ь — т) — + %-1(Ь — т)]}.

где т = тг-1(Ь) = Ьег-1(Ь) + иг-1(Ь — Ьег-1(Ь)) / ю; р = 0 в приближении больших потерь, но р = [уг-1(Ь)/о/2]2 в приближении одного прохода;

КаЬг,з (Ь) = (1 — Щп2з ¡к [а2(Ь) + Ь2(Ь)],

Каг,г-1(Ь,Ь — т) = (1 — К)п2г ¡к аг(Ь)аг-1(Ь — т), КЬг,г-1(Ь,Ь — т) = (1 — К)п2г ¡к Ьг(Ь)Ьг-1 (Ь — т).

Термин «точечное приближение» означает, что все множество точек поперечного сечения светового пучка в НКИ - в зависимости от вида крупномасштабного преобразования поля элементом О в контуре обратной связи - разбивается на бесконечное число независимых друг от друга подмножеств. Для случая поворота поля на 120о в поперечной плоскости пучка, изображенного на рис. 1, б, подмножества состоят из трех точек - вершин равностороннего треугольника. Независимость подмножеств понимается в том смысле, что процессы (нелинейные и интерференционные) в различных подмножествах протекают независимо. То есть между полем и участками нелинейной среды (а также между самими полями) физическое взаимодействие отсутствует, если они соответствуют различным подмножествам.

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

собой цепочки точек, «связываемые» последовательным смещением светового луча в поперечной плоскости пучка (благодаря действию элемента G).

Иными словами, луч света, проходя через НС и контур обратной связи НКИ в точке i (например, i = 1, 2, 3 на рис. 1, б), приобретает фазовый набег Ui и испытывает временную задержку tei. Из-за наличия элемента G луч попадает в точку i + 1. Здесь, «складываясь» с одним из входных лучей интерферометра, он воздействует на темп изменения величины нелинейного фазового набега Ui+i. Отметим, что в силу замыкания лучей на рис. 1 значение индекса i + 1 = 4 следует положить равным i + 1 = 1. Именно так набег Ui в точке i влияет на набег Ui+i в точке i + 1. Данный тип точек, согласно терминологии [15], называют транспозиционными точками, а m - порядком транспозиции. Соответственно будем говорить о цепочках транспозиционных точек (ЦТТ). Иллюстрация к понятию ЦТТ дана на рис. 1, б. При такой организации обратной связи траектория луча замыкается после m обходов НКИ. Согласно принятому способу нумерации транспозиционных точек, в случае замкнутой ЦТТ под записью i +1 подразумевается операция (i mod m) +1, где символ i mod m означает остаток от деления i на m. Физически это означает, что луч из m-й точки попадает в первую.

Понимаемое указанным образом точечное приближение позволяет из общих соотношений в [16] получить систему ОДУ (6), описывающую динамику нелинейного фазового набега для одной ЦТТ. Нетрудно заметить, что модель (6) является частным случаем (2).

Подобно выполненному переходу от (2) к (4) преобразуем модель (6) - в предположении, что dU/dt = 0 (статический режим), а величины ai, bi, ф^ Yi, tei постоянны во времени, - к виду

Ui = Kabi,j + pKab—ii + [Yi-i/o]x x{Kai,i-i  cos[(1 + я)(Ф-1 + Ui-i) + фi - ф— + % - %-i]+ (7) +Kbi , i-i  cos[(1 - я)(Ф—1 + Ui-i) + фi - ф— - % + %i-i]},

где Фi = rntei, Kabi,j = (1 - R)n2j lk [a2 + b2], Kai, i-i = (1 - R)n2i lkaia—i, Kbi,i-i = (1 - R)n2i lkbibi-i. По своей математической форме выражение (7), как известно [5], является одномерным ДО.

В случае однородности (в пределах одной ЦТТ) оптических свойств НС НКИ (n2 = n2j) и амплитуд входного поля (a = ai, b = bi) верно соотношение: Kab = Ka + Kb. Удобно ввести суммарный параметр нелинейности K и долю Qa интенсивности компоненты с частотой (1 + д)ш по правилу K = Kab = (Ka + Kb), Qa = Ka/K. Тогда Ka = K  Qa, Kb = K(1 - Qa). В случае однородности остальных оптических свойств НКИ (Ф = Ф^ y = Yi) и входного поля (% = 0, фi = 0) нетрудно получить ДО

Ui+i = K {1 + p + Y{Qa cos[(1 + д)(Ф + Ui)] +

+(1 - Qa) cos[(1 - д)(Ф + Ui)]}/o}.

Напомним, что в приближении больших потерь p = 0, в приближении одного прохода p = (y/o/2)2; в случае замкнутых ЦТТ i е [1, m], и запись i + 1 обозначает операцию (i mod m) + 1, то есть Ui = Um+i; в случае незамкнутых ЦТТ Ui = K, а

г € [1, +то), если ЦТТ бесконечна. Очевидно, что ДО (8) является частным случаем ДО (5).

В случае монохроматического излучения на входе НКИ (д = 0) и при р = 0, о =1 из (8) получим следующее ДО [17, 18]:

иг+1 = К [1+ у ес8

3. Демонстрация явления пространственного детерминированного хаоса в нелинейном кольцевом интерферометре

Проиллюстрируем соображения, выдвинутые в разделе 1, на примере модели динамики нелинейного фазового набега в НКИ (8). Поскольку ДО (8) - частный случай ДО (5), то для него справедливы утверждения, сделанные для (5) при малом числе т динамических переменных системы (2); большом т, т = то, (не)замкнутости системы (5), интерпретируемой как система уравнений.

Тип ЦТТ, зависящий от пространственного преобразования светового поля, осуществляемого линейным элементом О (см. рис. 1, а, рис. 3), задает структуру взаимной зависимости значений динамических переменных, упоминаемую в дефиниции ПДХ, и определяет, какой именно из перечисленных случаев реализуется. Для простейших преобразований поля в О подмножества Рг = {иг}.

Пусть ЦТТ незамкнута. Тогда, согласно смыслу модели (6) на языке ОДУ, неизбежно имеет место установление (релаксационно-колебательное либо асимптотическое) статического режима, устойчивого во времени. А моделирование на языке

Простейшие типы преобразования пучка элементом О

Поворот

А = 2пМ/т

Сдвиг 5

Сжатие 1/о

Растяжение о

Тип цепочки транспозиционных точек

Замкнутая конечная

Незамкнутая конечная

Незамкнутая бесконечная

т = —оо

Формирующиеся идеализированные структуры

т

зо

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

ДО (8) того же самого физического процесса в НКИ с ростом дискретной эволюционной переменной позволяет ожидать появление неустойчивого (и потому хаотического) поведения, то есть ПДХ.

Действительно, соответствующая ДО (9) зависимость иг от г при увеличении параметра нелинейности К в некоторых интервалах его значений демонстрирует весьма сложное поведение (рис. 4). При этом в поперечной плоскости лазерного пучка в НКИ формируется сложное распределение нелинейного фазового набега и (г) (см. рис. 3).

Нетрудно убедиться, что ПДХ имеет место, выполнив комплекс расчетов и построений: дискретного пространственного распределения (эквивалент временной реализации системы), фазового

портрета, спектра Фурье, ляпуновских характеристических показателей (ЛХП), фрактальной размерности. Они приведены на рис. 5 для модели (8) в случае незамкнутой ЦТТ. Из рис. 5 очевидно, что аттракторы в ДО могут иметь положительные ЛХП (Л = 0.34457) и дробную фрактальную размерность (Дз = 0.97994), а также характерный для детерминированного хаоса фазовый портрет и спектр Фурье.

Рис. 4. Зависимость от параметра нелинейности К значений нелинейного фазового набега и в пяти первых точках незамкнутой ЦТТ для статического режима НКИ при сдвиге элементом О оптического поля (вдоль оси Ох в плоскости хОу). Цифры обозначают номер точки в ЦТТ

Рис. 5. Характеристики ПДХ, описываемого ДО (8) при незамкнутой ЦТТ: дискретная пространственная реализация (а), фазовый портрет (б), спектр Фурье для амплитуды (в) и фазы (г) для параметров модели: К = 5, у = 0.6, Ф = 0, д = 0, =0

Рис. 6. Характеристики ПДХ, описываемого ДО

В случае же замкнутой ЦТТ оперирование моделью (6) на языке ОДУ позволяет выявить не только устойчивые, но и неустойчивые статические состояния, а также их бифуркации. При этом, одному и тому же набору параметров могут соответствовать несколько (не)устойчивых состояний. Очевидно, что если считать ЦТТ ДО бесконечной, то в математически идеальной ситуации при помещении начального условия ДО (8) на одну из ветвей бифуркационной диаграммы ДО даст периодическую последовательность значений Щ с периодом, не превышающим числа ш точек в ЦТТ. Эта ситуация обсуждалась выше, когда затрагивался вопрос о малости величины ш.

С физической точки зрения ясно, что структуры, соответствующие неустойчивым статическим состояниям, формироваться в НКИ не смогут - в отличие от устойчивых. Но переносится ли характеристика устойчивости статического состояния как целого (состояния системы (6)) на устойчивость соответствующего режима в отображении (8)? И обладает ли режим в отображении сложностью (оцениваемой хотя бы на уровне визуального восприятия)? С целью выяснения был предпринят расчет для замкнутой ЦТТ, состоящей из ш = 126 точек. На рис. 6 - для устойчивого (Л = -0.00807, Б0 = 0.17549) и неустойчивого (Л = 0.89427, Б0 = 0.96147) состояний - представлены: дискретные пространственные реализации, фазовые портреты, спектры Фурье для амплитуды и фазы. Результаты расчета позволяют ответить утвердительно на заданные вопросы.

Заключение

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

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

Получена модель динамики нелинейного фазового набега в нелинейном кольцевом интерферометре на языке дискретных отображений. На примере последнего иллюстрированы особенности ПДХ. Понятие цепочек транспозиционных точек раскрыто в двух аспектах: 1) зависимости структуры цепочки от пространственного преобразования светового поля, производимого элементом О в нелинейном кольце-

вом интерферометре; 2) влияния структуры цепочки на (не)замкнутость системы алгебраических уравнений, задающей ДО, и (бес)конечность их числа итераций m — 1, а также на вид оптической структуры U(r, t), формирующейся в поперечном сечении пучка. Для случаев (не)замкнутых цепочек выполнено моделирование на языке дискретных отображений и построены: дискретные пространственные распределения, фазовые портреты, спектры Фурье, рассчитаны ляпуновские характеристические показатели, значения фрактальной размерности. Показано, что режим в отображении и для замкнутой цепочки транспозиционных точек обладает сложностью.

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

Авторы признательны рецензенту за критику, стимулировавшую переоценку результатов, сокращение текста в раздела 1 и уточнение формулировок.

Работа поддержана грантом Федерального Агентства по образованию Ми-нобрнауки РФ (Программа 2005 г.: «Развитие научного потенциала высшей школы». Подпрограмма 3. Раздел 3.3), регистрационный № 60321.

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

1. Ikeda K. Multiple-valued stationary state and its instability of the transmitted light by ring cavity system // Opt. Comm. 1979. Vol. 30, № 2. P. 257-260.

2. Ахманов С.А., Воронцов М.А. Нестабильности и структуры в когерентных нелинейно-оптических системах, охваченных двумерной обратной связью // Нелинейные волны: динамика и эволюция: сб. ст. М.: Наука, 1989. С. 228-237.

3. Розанов Н.Н.Оптическая бистатбильность и гистерезис в распределенных нелинейных системах. М.: Наука, 1997. 336 с.

4. Измайлов И.В., Калайда В.Т., Магазинников А.Л., Пойзнер Б.Н. Бифуркации в точечной модели кольцевого интерферометра с запаздыванием и поворотом поля // Изв. вузов. Прикладная нелинейная динамика. 1999. Т. 7, № 5. С. 47-59.

5. Кузнецов С.П. Динамический хаос. Курс лекций. Учебное пособие для студентов вузов, обучающихся по физическим специальностям. М.: Физматлит, 2001. 296 с.

6. Балякин А.А. Исследование хаотической динамики кольцевого нелинейного резонатора при двухчастотном внешнем воздействии // Изв. вузов. Прикладная нелинейная динамика. 2003. Т. 11, № 4-5. C. 3-15.

7. Берже П., Помо И., Видаль К. Порядок в хаосе. О детерминистском подходе к турбулентности. М.: Мир, 1991. 368 с.

8. Мельников Л.А., Конюхов А.И., Рябинина М.В. Динамика поперечной поляризационной структуры поля в лазерах // Изв. вузов. Прикладная нелинейная динамика. 1996. Т. 4, № 6. С. 33-53.

9. Weiss C.O. et al. Generation of optical vortices in laser field // Appl. Phys. B 1999. Vol. 68. P. 151-168.

10. Измайлов И.В., Пойзнер Б.Н. Варианты реализации нелинейно-оптического устройства скрытой передачи информации // Оптика атмосферы и океана. 2001. Т. 14, № 11. С. 1074-1086.

11. Melnikov L.A., Konukhov A.I., Veshneva I.V. et al. Nonlinear dynamics of spatial and temporal patterns in lasers and atom optics: Kerr-lens mode-locked laser, Zeeman laser and Bose-Einstein atomic condensate // Izv. VUZ. Applied nonlinear Dynamics. 2002. Vol. 10, № 3. P. 40-62.

12. Рыскин М.И., Иванов А.В. Нелинейная динамика в науках о Земле // Изв. вузов. Прикладная нелинейная динамика. 2003. Т. 11, № 6. С. 138-148.

13. Измайлов И.В., Пойзнер Б.Н., Денисов П.Е. Равносильность: от обоснования понятия до анализа бифуркационного поведения. Томск: ТГУ, 2003. 46 с. Вестник Том. гос. ун-та. Бюллетень оперативной научной информации. 2003. № 15. Октябрь 2003.

14. Заславский Г.М., Кириченко Н.А. Хаос динамический // Физическая энциклопедия / Гл. ред. А.М. Прохоров. М.: Сов. энциклопедия. Т. 5, 1998. С. 397-402.

15. Аршинов А.И., Мударисов Р.Р., Пойзнер Б.Н. Механизмы формирования простейших оптических структур в нелинейном интерферометре Физо // Изв. вузов. Физика. 1995. № 6. С. 77-81.

16. Измайлов И.В., Магазинников А.Л., Пойзнер Б.Н. Моделирование процессов в кольцевом интерферометре с нелинейностью, запаздыванием и диффузией при немонохроматическом излучении // Изв. вузов. Физика. 2000. № 2. С. 29-35.

17. Шергин Д.А., Измайлов И.В. Дискретные отображения как средство описания детерминированного пространственного хаоса // Сб. тез. 9-й Всерос. науч. конф. студентов-физиков и молодых ученых: В 2-х т. 2003. Екатеринбург-Красноярск: АСФ России, 2003. Т. 2. С. 90-93.

18. Шергин Д.А., Измайлов И.В., Пойзнер Б.Н. Дискретные отображения как язык описания пространственного детерминированного хаоса // Современные проблемы физики и высокие технологии: Матер. Междунар. конф. (29 сентября -4 октября 2003 г., г. Томск). Томск: Изд-во НТЛ, 2003. С. 186-189.

Томский государственный Поступила в редакцию 30.08.2004

университет После доработки 29.03.2005

SPATIAL DETERMINISTIC CHAOS: THE MODEL AND DEMONSTRATION OF PHENOMENON IN COMPUTING EXPERIMENT

I.V. Izmailov, A.V. Lyachin, B.N. Poizner, D.A. Shergin

The concept of spatial deterministic chaos is justified. An attempt to give its set-theoretic definition is undertaken. Transition from the ordinary differential equations to discrete maps without use of an approximation of the instantaneous response is realized for mathematical description of spatial deterministic chaos. The developed theoretical theses are applied for deriving a dynamics model in terms of discrete maps of nonlinear phase shift in a ring interferometer. In case of the model discrete realizations, phase portraits Fourier's spectrums illustrating peculiarities of spatial deterministic chaos in the ring interferometer are constructed. A concept of discrete maps undergoing an evolution is introduced.

Измайлов Игорь Валерьевич - родился в 1976 году в Алма-Ате. Поступил на радиофизический факультет Томского государственного университета (1993). Аспирант кафедры квантовой электроники и фотоники (1999). Область научных интересов: моделирование процессов в нелинейных системах различной природы, применение нелинейной динамики для обработки и скрытой передачи информации.

Лячин Александр Владимирович - родился в Томске (1979). Окончил с отличием радиофизический факультет Томского госуниверситета (2002). Аспирант кафедры квантовой электроники и фотоники радиофизического факультета. Область научных интересов - моделирование и изучение процессов в нелинейных динамических системах, динамический хаос и его приложение в системах обработки и передачи информации. Автор более 40 публикаций. E-mail: 765-alex.rff@elefot.tsu.ru

Пойзнер Борис Николаевич - родился в Томске (1941), окончил радиофизический факультет Томского государственного университета. Защитил кандидатскую диссертацию по теории колебаний и волн (1970), доцент кафедры квантовой электроники и фотоники ТГУ. Читает лекции по нелинейной оптике, физике, физике лазеров, принципам управления лазерным излучением, основам синергетики. Область научных интересов: квантовая электроника, применение нелинейной динамики в оптике и материаловедении, прикладная наукометрия, культурологическая теория образования. Имеет много статей по указанной тематике. Инициатор подготовки и редактор семи библиографических указателей (в том числе «Синергетика и сопредельные науки», «Университетское образование и его социальная роль», «Интеллигенция в российском обществе и университете», «Психика и интеллект обучаемого»). Действительный член Всероссийского общества библиофилов. E-maik: Pznr@elefot.tsu.ru

Шергин Денис Александрович - родился в 1982 году, в настоящее время заканчивает обучение в магистратуре радиофизического факультета Томского государственного университета. Область научных интересов - применение методов нелинейной динамики в оптике и криптографии; современные телекоммуникационные системы. Имеет 5 публикаций.