Качественный анализ сингулярно возмущенных моделей одного класса оптико-электронных систем

Изв. вузов «ПНД», т. 16, № 4, 2008 УДК 517.929

КАЧЕСТВЕННЫЙ АНАЛИЗ СИНГУЛЯРНО ВОЗМУЩЕННЫХ МОДЕЛЕЙ ОДНОГО КЛАССА ОПТИКО-ЭЛЕКТРОННЫХ СИСТЕМ

Д.В. Глазков

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

Введение

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

Особый класс моделей лазерных систем представляют собой модели, учитывающие воздействие запаздывающей оптической обратной связи. Их сложность обусловлена тем, что этот класс задач описывается в терминах дифференциальных уравнений с запаздыванием, которые обладают бесконечным числом степеней свободы. Кроме того, динамические системы с запаздыванием демонстрируют различные типы неустойчивостей, обусловленные воздействием задержки. Например, для полупроводниковых лазеров с оптической обратной связью наблюдаются различные пути перехода к хаосу: квазипериодический [3], через каскад бифуркаций удвоения периода [4], сценарий Икеды [5], перемежаемость и кризис аттракторов [6,7].

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

колебания Петермана-Тейгера [8,9], регулярные импульсные пакеты [10,11], низкочастотные флуктуации [2,7,12,13], когерентный коллапс [2-4,7,12]. Глубокое понимание механизмов их возникновения открывает новые способы для управления режимами генерации, а значит, и новые возможности в приложениях. В таких областях, как хранение данных или оптические и оптоволоконные коммуникации отражения и связанные с ними сложные сопутствующие эффекты неизбежны. Например, искажения сигнала при передаче данных нередко обусловлены неминуемыми отражениями от торцов волноводов. Поэтому изучение влияния оптической запаздывающей обратной связи на работу лазеров разных типов представляет собой важную прикладную задачу.

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

1. Постановка задачи

Рассмотрим математическую модель динамики одномодового лазера, которая учитывает воздействие отраженного излучения на резонатор и основана на уравнениях Ланга-Кобаяши [14]

((Е

— = v(1+iа)EZ + ув~гю°кЕ (г-Н),

И

^ = Q - Z - (^) \Е\2.

Здесь Е(г) - комплексная амплитуда электрического поля; Z(г) - инверсия носителей; у>0 и (-ю0Н) - сила и фаза обратной связи; ю0 - оптическая частота генерации в отсутствие обратной связи; Q - превышение током накачки первой пороговой величины; V - отношение времен затухания инверсии носителей и фотонов в резонаторе; а - коэффициент уширения линии, отвечающий за нелинейное взаимодействие между амплитудой и фазой поля; Н - время прохода излучения по внешнему резонатору, нормированное в единицах времени затухания инверсии.

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

(1Е

— = v(1+iа)EZ + у/,

/ = iA/ + Л[Е

^ = Q - Z - (^) \Е\2.

В этой системе / (г) - комплексная амплитуда электрического поля на выходе из фильтра, Л - ширина спектра, А - расстройка между частотой излучения уединенного лазера и несущей частотой фильтра.

<

Система (1) при больших значениях параметров изучалась в работах [13,16]. Другие исследования такого рода автору не известны.

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

Методы построения систем с фазовым пространством бесконечной размерности, играющих роль нормальных форм, были предложены в работах [17-21].

Наличие в системе большого параметра приводит к необходимости рассматривать сингулярно возмущенную задачу. Задачи этого типа [22] не могут быть качественно исследованы регулярными методами. Главное преимущество используемой в работе методики состоит в переходе от сингулярно возмущенной задачи к регулярным уравнениям, которые больших параметров уже не содержат.

2. Уравнения Ланга-Кобаяши при

Рассмотрим уравнение (1) при условии Q=e-1, где е - малый положительный параметр. Нормируя переменные следующим образом:

г = ев, Е = -1  в-™3Ег,

Vе

сведем исходную систему к частному случаю следующего векторного уравнения:

^ = Е

Общая методика исследования его динамики была разработана в [18,19]. В соответствии с алгоритмами, представленными в этих статьях, в малой окрестности некоторого периодического решения системы «нулевого приближения» йх/йв = Е(х) строится нормализующая замена, позволяющая получить асимптотические оценки системы (3). Эта замена имеет следующий вид:

У

йв =1 + еф

Здесь (г, т) являются Т-периодичными по т, причем Уо(в) есть периодическое решение системы ОДУ йх/йв = Е(х), ф(£) - скалярная почти периодическая функция.

Отметим, что из (5) можем получить

о

т(в—Не-1) & т(в) — у, где у = Не-1 + ^ ц>(г+т)йт.

-1 -1 +

-н

Стандартная линеаризация йх/йв = Е(х) на У0(в) приводит к системе

ЛТв = А<»>

где А

Подставляя ряды (4), (5) в (3) и приравнивая коэффициенты при одинаковых степенях е, получим цепочку уравнений вида

dVk йт

= А(х)Уй + {к (г, Т, т—у), кеп,

где через (к(г, т, т—у) обозначена неоднородность дифференциального уравнения. На первом шаге из условий периодичности функции VI по т получим систему урав-

нений

{(г(г, т, т—у), И(т)) = 0,

где угловыми скобками ( ) обозначено скалярное произведение

т

(X, Y) = Ц(Х(т), Y(т))dт.

Полученные уравнения (8) играют роль нормальной формы в задаче о локальной динамике уравнения (3) в некоторой окрестности решения Vо(в) системы «нулевого приближения» йх/йв = Р(х). Отыскав решение VI системы (7) в указанном классе функций, можно последовательно найти любое число элементов рядов (4), (5).

В случае уравнений Ланга-Кобаяши при Q»1 система ОДУ вида йх/йв=Е(х) является консервативной и имеет семейство периодических решений, которое в форме вектора с вещественными компонентами может быть записано следующим образом:

V0(в) = (ссо8(юв), с8т(шв), с-2—1)т.

Матрица системы (6) имеет вид

А(в) = ^ (Vо(в)) =

( 0 — ю 0 \

ю 0 0

\ —2 с-1 со8(юв) —2 с-1 8ш(юв) —с2 /

Для системы (6) легко находятся все линейно независимые решения, два из которых можно выбрать периодическими

Ко

( со8(юв) \

8т

( — 8ш(шв) \

К1

С08

К2

V

/

0

\ е-с 3 /

2п с2 _

юс <1.

Очевидно, мультипликаторы системы (6) имеют вид: Ц1=Ц2=1, Цэ=е~

Линейно независимые периодические решения сопряженной к (6) системы выберем таким образом, чтобы (К^И^) = ,г] = 0,1, где - символ Кронекера. Таким условиям удовлетворяют вектор-функции вида

0

0

/ сс^тз) \

Но

8т(мв) 0

( — 8т(мв) ^

Н1

V

СС8(Ю8) 0

/

Рассмотрим свободную постоянную с как функцию «медленного» времени: Ь=е,в. В соответствии с изложенной методикой, собирая коэффициенты при первой степени е, придем к системе

дУо (с(г) дУо дVI ,ОЛ

0 + ^^Ф

дс М

дт

д т

Условие существования 2л/ш-периодического (по т) решения уравнения (9) состоит в ортогональности его неоднородной части функциям Но(т), Н1(т). Выполнив ряд преобразований и переобозначив шф^ф, получим следующую систему интегрально-дифференциальных уравнений, которую мы будем называть квазинормальной формой модели Ланга-Кобаяши в случае

' dc(í) / 1 \ ( о \

- = V I -^тг — с

(И

ф(£) = va

\с(1) 1

#ц) у у8КШоН+_{ф('+т)(Ос(с-г•

1

Заметим, что система (10) заменами

(0

ф(<) = с(£)е*0(4) = Е (I) может быть преобразована к одному комплексному уравнению следующего вида:

(Е

— = v(1+ш) (|Е|-2—1) Е + уе-™0кЕ(г—Н),

которое, очевидно, связано с (1).

Таким образом, при достаточно больших значениях накачки Q оказывается возможно исключить из системы (1) уравнение для инверсии Z, упростив первоначальную задачу.

Стационарные решения системы (10) определяются из трансцендентного уравнения относительно ф

Ф = —у [а сс^юоН+фН) + 8ш(шоЛ,+ф^)] .

Обозначая п=юоЛ,+фЛ,, придем к хорошо известному соотношению для определения числа мод внешнего резонатора [7,23]

П — юоЬ = —уЬ\/ 1+а2 8ш(п + аг^(а)),

которое, имеет, по крайней мере, одно решение. Легко видеть, что, чем больше значение yhV 1+а2, тем больше решений имеет уравнение (12).

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

X2 + 2Ху cos nk (1—e-Xh) + у2 (1—e-Xh)2+

+2(v—уcosnk)[X + y(cosnk—asinnk)(1-e ^h)j =0.

В случае cos nk=1 для моды с максимальным усилением [23] уравнение (13) сводится к следующей простой совокупности:

" X + y(1-e-Ah) = 0 ,

X + y(1-e-Ah) +2(v-y) = 0.

При 0<y<v у каждого из уравнений (14) отсутствуют корни с положительной вещественной частью. Следовательно, мода с максимальным усилением в случае асимптотически большой накачки является устойчивой. Этот результат хорошо согласуется с выводами, сделанными в работе [23].

Мода с минимальной шириной линии достигается при nk=— arctg(a). Достаточным условием ее устойчивости является следующая теорема:

Теорема 1. Пусть nk=— arctg(a) и у < yi = v/\/1 + a2. Тогда все корни уравнения (13) имеют неположительные действительные части.

Критерием устойчивости моды с минимальной шириной линии в случае большого запаздывания h=p,-1, где ц - малый параметр, является следующая теорема:

Теорема 2. Пусть nk=— arctg(a) и у < у2 = (2vV 1 + a2)/(2 + a2). Тогда все корни уравнения (13) имеют неположительные действительные части. В случае Y>Y2 существует корень с положительной действительной частью.

Непосредственной подстановкой несложно проверить, что корней с асимптотикой шц-р, где p>0, у уравнения (13) в случае nk=— arctg(a) нет. Раскладывая X в ряд по степеням малого параметра ц, так что X=Xo+^Xi+p,2X2+..., получаем Re X0=Re X1 =0, а знак Re X2 определяет знак выражения y(2+a2)— 2v\/ 1+a2.

Основные результаты относительно динамических свойств системы (10) получены численным моделированием. При этом исследовалась не сама система (10), а система дифференциальных уравнений с запаздыванием, получаемая заменой ф(t)=d0(t)/dt. Тем не менее, в качестве основных характеристик модифицированной таким образом квазинормальной формы выступали функции c(t), d0(t)/dt.

Численное интегрирование выполнялось с помощью схемы Рунге-Кутты 4-го порядка, адаптированной для систем с запаздыванием. В качестве начальных условий выбирались функции из класса тригонометрических полиномов, заданных на отрезке [—h, 0]. Были обнаружены явления мультистабильности (рис. 1), кризиса аттракторов, последовательного удвоения периода циклов, определены области с хаотической и сложной нехаотической динамикой.

Рис. 1. Иллюстрация явления мультистабильности. При значениях параметров г=100, а=5, rn0h=1, у=97.9, h=0.04 найдено 3 различных устойчивых режима: мода, близкая к состоянию с максимальным усилением, цикл удвоенного периода

Важные особенности динамики квазинормальной формы (10) отражает бифуркационная диаграмма, представленная на рис. 2. Выбраны следующие значения параметров v=100, а=5, rn0h=1. Цифрами обозначены области: 1 - состояние равновесия; 2 - цикл, рождающийся в результате бифуркации Андронова-Хопфа; 3 - сценарий Фейгенбаума - каскад бифуркаций удвоения периода (штриховой линией отмечены первые две бифуркации удвоения); 4 - нефейгенбаумовские переходные режимы (в некоторых случаях интерпретировались как квазипериодические или двухчастотные на основе анализа временного ряда стандартной функцией FFT пакета Maple); 5 - хаос. Штрихпунк-тиром обозначены первые две бифуркации седло-узла, в которых происходит рождение пары мода-антимода. Правее первой из них существуют области мультиста-бильности. Изменение бассейнов притяжения лежит в основе перехода (при возрастании параметров h, у) от все усложняющихся режимов из области 3 к области 2. Область, отмеченная цифрой 4 на верхней части диаграммы имеет сложную струк-

Рис. 2. Бифуркационная диаграмма системы (10) на плоскости (у, К), разделенная на две части с разным масштабом по параметру К. Отметим существование областей мультистабильности

туру и содержит многочисленные участки типа 2 и 5 малых размеров. Отметим, что существование устойчивого состояния равновесия, близкого к моде с максимальным усилением, подтверждается численно. При у>^у1+а2 система (1), вообще говоря, теряет свойство диссипативности. Это объясняется тем, что столь большие значения у физически недостижимы, поскольку величину v-1y следует понимать как ту долю отраженного излучения, которая возвращается в резонатор. Граница области динамического хаоса определялась на основе вычисления корреляционной размерности [24,25]. Расчет старшего показателя Ляпунова [26], в целом, подтверждает полученные оценки.

Подчеркнем связь между установившимися режимами в квазинормальной форме и в исходной модели: устойчивым состояниям равновесия и циклам системы (10) соответствуют устойчивые предельные циклы и торы в (1).

Следует также обратить внимание на то, что известные о когерентном коллапсе факты [2-4,7,12] хорошо согласуются с наблюдаемой динамикой системы (10).

3. Модель с фильтром при V,

Введем в рассмотрение малый положительный параметр е так, что

V = е-1, Q = де-1, и выполним замену переменных

г = ев, Е = —  е-™3Е1, / = —  е-™3¡1. ее

Это приводит к системе, которая является частным случаем (3).

Решение системы «нулевого приближения» представимо в форме

У0(в) = уд (ехр(гюв), сехр(гюв), 0)т.

Матрица линеаризованной на Уо системы имеет вид

/ 0 —ю 0 0 ^д  [еов(юв) — а 8т(юв)] \

ю 0 0 0 уд  ^т(юв) + а еов(юв)]

0 0 0 -ю 0

А

ю 0

\  еоз(юв) —2уд  8т(юв) 0 0 —д /

Система (6) и сопряженная к ней имеют по три (с точностью до константы) периодических решения и им соответствующих единичных мультипликатора. Выберем эти решения следующим образом:

К0(в) = уд  (— 8ш(юв), еоз(юв), 0, 0, 0)т ,

1 т Но(в) =-- (8ш(юв)+аеоз(юв), — еов(юв)+а8ш(юв), 0, 0, 0) ,

К (в) = Н1(в) = (0, 0, еов(юв), вт(юв), 0)т ,

К2(в) = Н2(в) = (0, 0, — вт(юв), еов(юв), 0)т ,

0

0

0

чтобы выполнялось условие (К^, И) = ,1,] =0...2, где б^ - символ Кронекера.

В данном случае замена (4), (5) несколько модифицируется, а именно: V(s,е) = Уо(х1, Т2) + eVl(í, ть Т2) + ть Т2) + ...,

= 1 + еф1(£) + е2% (*) +... , йт2 п

— =1+ £ф2(£)+ е2^2 (*) + ... .

Здесь V0(т1,т2) =  (ехр^шт^, сехр(гшт2), 0)т. Кроме того,

я г

Т2(«)-Т1(«) = еJ [ф2(ет)-ф1(ет)+ .. ]йт = J [ф2(т)-ф1(т)]^т + О(е).

«о го

По аналогии с предыдущим пунктом введем в рассмотрение у1 и у2. Тогда, собирая коэффициенты при первой степени е, получим систему уравнений

д Vо йс(г) дУо д Vо дVl дVl

+ ф1

= AVl + Ф(Уо (т1, т^ , ^(т1-уь Т2-У2)).

Из условий периодичности функции Vl по Т1, Т2 получим следующие уравнения, которые будем называть квазинормальной формой системы (2) в выбранной области пространства параметров,

_ г

ф^) = -с(£)у\/ 1+а2 8т(/(ф1(т)-ф2(т))^т + arctg(а)),

го

¿с({) о г

—— = Л[со8

^ -н го

Л о г

ф2^) = А--— 8ш(шоЛ + / ф^+т-йт +/(ф2(т)-ф1(т))&) .

с(Ч -н го

(Ш  (^

Здесь и далее без ограничения общности считаем, что б1(^о)=б2(^о), где —^—=ф^(£).

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

Подчеркнем связь между установившимися режимами в квазинормальной форме и в исходной модели: устойчивому состоянию равновесия системы (16) соответствует устойчивый предельный цикл в (2).

4. Большая накачка Q^1 в системе с фильтром

Пусть параметр Q достаточно велик. В данном случае считаем, что v=0(1). Положим Q=е-1, где е^ 1. Выполним следующее преобразование переменных:

г = ев, Е = —  е~гюв

Vе

1

—е

Е1, / = —  е-гюз/1.

Полученная в результате система также сводится к уравнению общего вида (3).

Решение системы «нулевого приближения» можно записать следующим образом:

У0(в) = (с1 ехр(гюв), с2ехр(гюв), с-2—1) .

Стандартная линеаризация dx/dг=F(x) на Уо(в) в соответствии с изложенной методикой приводит к системе (6), 2п/ю-периодическая матрица которой имеет вид

/

(Ао(в)) = ^(Уо (в)) =

0

ю

ю

00

00

0

0 -ю 0

ю0

\ —2 с-1 еов(юв) —2 с-1 8т(юв) 0 0 —с2 /

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

В соответствии с изложенной методикой, собирая коэффициенты при первой степени е, придем к следующей системе уравнений:

дУо (1сл(г) дУо ¿¿2(1) , (+) дУо , (+) дУо , дУ1 дУ1

+ —--+ ф1(г^—+ ф2(г)——+ ——+

дс1 ¿г

дс2 ¿г

дт1

дт2 дт1 дт2

АУ1 + Ф(Уо(Т1, Т2), Уо(т1—У1, Т2— У2)).

Условие существования периодического (по т) решения уравнения (17) состоит в ортогональности его неоднородной части функциям Но(т1), Н1(т1), Н2(т2), Н3(т2). Отметим, что, выполнив замену юф ^ ф, можно избавиться от ю.

Следующая система интегрально-дифференциальных уравнений выступает в роли нормальной формы в задаче о локальной динамике уравнений (2) в окрестности решения Уо невозмущенной предельной системы при

= V

— с1(г)) + У (ф2(т) —ф1 (г))йг )с2(г) ,

¿с1

¿г . _ , _

Ф1(г) = сЩ — ^ + 7 в1П(/ (ф2(г) —Ф1 М)^) ,

Г / о 4 \

Л с1(г—Ь) еоЦ юоН + / ф1 (г+т^т + /(ф2(г)—ф^т))¿т) — с2(г) V -н г0 /

( о 1 \ сл(г—ю

ф2(г) = А — Л 8Ш юок + / ф1(г+т)(1т + / (ф2(т)—ф1(т)) ^ —-—.

V -н г0 / с2

^2 (г) ¿г

0

0

0

0

0

0

0

Заметим, что квазинормальная форма (18) при выполнении замен

вд = с1

допускает более компактную запись в терминах исходной системы

- = v

= гА/ + Л[-(¿-Л)е-*юоН - /].

Численный анализ показывает, что динамика (18) оказывается существенно проще, чем (10). В пределах точности вычислений при у^ и 0<Л<10 вне зависимости от значений других параметров удавалось обнаружить лишь один аттрактор - устойчивое состояние равновесия, которому соответствует простой периодический режим в системе (2).

5. Нормализация в неавтономном случае

Полученные результаты распространяются на случай, когда параметры уравнений (1), (2) являются гладкими ограниченными периодическими функциями времени Ь, амплитуда и период колебаний которых не зависят от малого параметра. Сформулируем основной результат для системы (1). Пусть параметры уравнений Ланга-Кобаяши имеют вид [28,29]

а = а(Ь), у = у(Ь), Л = Л(Ь), Q = е-1д(Ь).

Если v=е-1v* (Ь) есть большой параметр, то соответствующая квазинормальная форма уравнений (1) записывается следующим образом:

о

ф

-н

В случае, когда v=v(t) есть величина порядка единицы, нормализованная система имеет вид

( ¿с(ь) а(Ь) о

лГ = v(t)^M^- - с(Ь)) + У(Ь) сс8(шоЛ(Ь) + / ф(Ь+т)Нт)с(Ь-Л(Ь)),

М с

а

ф(Ь) = v(t)a(t)( С2Ут - 1) - у(Ь) 81п(юоЛ(Ь) + / ф(Ь+т)Нт) У .

с

Детальное изучение этой модели представляет собой отдельную задачу.

Исследуем уравнение (20) в случае, когда его параметры представимы в виде

а=ао+а1 81п(ю„ь), у=уо +у1 8ш(шуь), Л=Ло+Л1 8ш(шнь),

где ао>|а1|, уо>|у11, Ло > | Л11. Если частоты юа, , юн соизмеримы, то в системе (20) устанавливается периодический режим. Такой эффект связан с пульсацией эллипса [2, 7, 9], который соотносят с модами внешнего резонатора. При этом возможны

переходы с одной моды на другую. Соответственно такие «медленные» изменения параметров могут приводить к усложнению динамики исходной модели. Однако, как отмечено, например, в [30], физически в модели Ланга-Кобаяши проще всего управлять параметрами Н и Q.

Особый практический интерес представляет изучение двух случаев. Это учет влияния параметрического резонанса и быстро осциллирующих слагаемых на квазинормальную форму. Поскольку у системы «нулевого приближения» есть пара единичных мультипликаторов, то параметры с резонансными слагаемыми имеют вид Q=е-1+q1 8ш(еОдг), у=уо+еу1 8ш(еОуг) и так далее. Из построения нормализованной системы можно установить, что резонансные слагаемые не влияют на вид уравнений (10). То есть реализовать параметрический резонанс в данном случае невозможно.

Наличие быстро осциллирующих составляющих приводит к несколько более содержательным результатам. Пусть параметры системы (1) имеют вид

У = уо+у1 вш(ге-2), Н = Но+Н1 8т(ге-2), а = ао+а1 8ш(ге-2),

V = v0+v1 8ш(ге-2), Q = е-1 [до+д1 эт(ге-2)].

Выполняя стандартные замены в уравнениях (1), усредняя их по «быстрому» време-

12

ни ве 1=ге 2, и действуя далее в соответствии с изложенным алгоритмом нормализации, получим следующую систему с распределенным запаздыванием:

¿с(г) = Vо(^0)—с(г)) +1 } ео8(юоН(т) + / ф(г+ь)йь)с(г—Н(т)) ¿т ,

-Н

dг с(г) п -1 у 1—т

ф(г) = ро(сд^—1) — П i -= «1п(юоН(т) + / ф^+п^п)сс(—НТ1 ¿т.

с2

Здесь

Р0 = vоао+-2-1, Н(т) = Н0+Н1Т.

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

Рассмотрим систему (21) в случае Vl=yl =а1=0. В терминах исходной задачи (1) это ситуация быстро осциллирующего запаздывания. Численный анализ системы (21) показывает, что уже при соотношении Н1/Но=0.01 происходит ощутимый сдвиг вправо по параметру у всех границ между областями, изображенными на рис. 2. При Н1/Н0=0.2 в системе (21) практически всюду в рассматриваемой на рис. 2 области параметров устанавливается стационарный режим, которому соответствует простой периодический режим в усредненной системе (1).

Аналогичные результаты справедливы и для модели с фильтром.

Заключение

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

накачки. Заметим, что аттракторы модельных задач (1), (2) и нормализованных уравнений связаны между собой стандартным образом: устойчивому состоянию равновесия или циклу, наблюдаемому в квазинормальной форме, соответствует устойчивый цикл или тор в исходной системе.

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

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

Автор выражает благодарность С.А. Кащенко и С.Д. Глызину за постановку задачи, помощь и внимание к работе, а также Е.В. Григорьевой за полезные консультации.

Работа выполнена при финансовой поддержке целевой программы «Развитие научного потенциала высшей школы» (проект РНП.2.1.1.630).

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

1. Ханин Я.И. Основы динамики лазеров // М.: Наука, Физматлит, 1999. 368 с.

2. Van Tartwijk G.H.M. and Agrawal G.P. Laser instabilities: A modern perspective // Progress in Quantum Electronics. 1998. Vol. 22. P. 43.

3. Grigorieva E.V. Quasiperiodicity in Lang-Kobayashi model of lasers with delayed optical feedback // Nonlinear Phenomena in Complex Systems. 2001. Vol. 4, № 4. P. 333.

4. Ye J., Li H., McInerny J.M.Period-doubling route to chaos in a semiconductor laser with weak optical feedback // Phys. Rev. A. 1993. Vol. 47, № 3. P. 2249.

5. Fischer I., Hess O., Elsasser W. and Gobel E. High-dimensional chaotic dynamics of an external cavity semiconductor laser // Phys. Rev. Lett. 1994. Vol. 73, № 16. P. 2188.

6. Sacher J., Elsasser W., Gobel E.O. Intermittency in the coherence collapse of a semiconductor laser with external feedback // Phys. Rev. Lett. 1989. Vol. 63, № 20. P. 2224.

7. Sano T. Antimode dynamics and chaotic itinerancy in the coherence collapse of semiconductor lasers with optical feedback // Phys. Rev. A. 1994. Vol. 50, № 3. P. 2719.

8. Tager A.A., Petermann K. High-frequency oscillations and self-mode locking in short external-cavity laser diodes // IEEE J. Quantum Electron. 1994. Vol. 30, № 7. P. 1553.

9. Wolfrum M., TuraevD. Instabilities of lasers with moderately delayed optical feedback // Opt. Commun. 2002. Vol. 212. P. 127.

10. Heil T., Fischer I., Elsasser W., Krauskopf B., Green K., Gavrielides A. Delay dynamics of semiconductor lasers with short external cavities: Bifurcation scenarios and mechanisms // Phys. Rev. E. 2003. Vol. 67. 066214.

11. TabakaA., Panajotov K., Veretennicoff I., Sciamanna M. Bifurcation study of regular pulse packages in laser diodes subject to optical feedback // Phys. Rev. E. 2004. Vol. 70. 036211.

12. Heil T., Fischer I., and Elsasser W. Influence of amplitude-phase coupling on the dynamics of semiconductor lasers subject to optical feedback // Phys. Rev. A. 1999. Vol. 60, № 1. P. 634.

13. Pieroux D., Mandel P. Bifurcation diagram of a complex delay-differential equation with cubic nonlinearity // Phys. Rev. E. 2003. Vol. 67. 056213.

14. Lang R., Kobayashi K. External optical feedback effects on semiconductor injection laser properties // IEEE J. Quantum Electron. 1980. Vol. 16, №1. P. 347.

15. Green K., Krauskopf B. Mode structure of semiconductor laser subject to filtered optical feedback // Opt. Commun. 2006. Vol. 258. P. 243.

16. Alsing P., Kovanis V., Gavrielides A. and Erneux T. Lang and Kobayashi phase equation // Phys. Rev. A. 1996. Vol. 53, № 6. P. 4429.

17. Kaschenko S.A. Normalization in the systems with small diffusion // International Journal of Bifurcations and chaos. 1996. Vol. 6, № 7. P. 1093.

18. Кащенко С.А. Бифуркации цикла в сингулярно возмущенных нелинейных автономных системах // Изв. РАЕН, серия МММИУ. 1998. Т.2, № 4. С. 5.

19. Кащенко С.А. Локальная динамика нелинейных сингулярно возмущенных систем с запаздыванием // Дифференциальные уравнения. 1999. Т.35, № 10. С. 1343.

20. Кащенко С.А. Бифуркации в окрестности цикла при малых возмущениях с большим запаздыванием // Журнал выч. матем. и матем. физ. 2000. № 4.

21. Kaschenko S.A. Bifurcational features in systems of nonlinear parabolic equations with weak diffusion// International Journal of Bifurcation and Chaos. 2005. Vol. 15, № 11.

22. Васильева А.Б., Бутузов В.Ф. Асимптотические разложения решений сингулярно возмущенных уравнений. М.: Наука, 1973. 272 c.

23. Levine A.M., Tartwijk G.H.M., Lenstra D. and Erneux T. Diode lasers with optical feedback: Stability of the maximum gain mode // Phys. Rev. A. 1995. Vol. 52, № 5. P. 3436.

24. Grassberger P., Procaccia I. Estimation of the Kolmogorov entropy from a chaotic signal // Phys. Rev. A. 1983. Vol. 28, № 4. P. 2591.

25. Grassberger P., Procaccia I. Measuring the strangeness of strange attractors // Physica D. 1983. Vol. 9, №1,2. P. 189.

26. Wolf A., Swift J.B., Swinney H.L., Vastano J.A. Determining Lyapunov exponents from a time series // Physica D. 1985. Vol. 16, № 3. P. 285.

27. Глазков Д.В. Простейшие устойчивые режимы в модели Ланга-Кобаяши с большим запаздыванием // Труды XXVII Конференции молодых ученых ММФ МГУ им. М.В. Ломоносова. 2005. С. 27.

28. Колесов Ю.С., Майоров В.В. Новый метод исследования устойчивости решений линейных дифференциальных уравнений с близкими к постоянным почти-периодическими коэффициентами// Дифференциальные уравнения. 1974. Т. 10, № 10. С. 1778.

29. Кащенко С.А., Майоров В.В. Алгоритм исследования устойчивости решений линейных дифференциальных уравнений с последействием и быстро осциллирующими почти периодическими коэффициентами // Сборник статей. Исследования по устойчивости и теории колебаний. 1977. С. 70.

30. Mendez J.M., LajeR., GiudiciM., AliagaJ. andMindlin G.B. Dynamics of periodically forced semiconductor laser with optical feedback // Phys. Rev. E. 2001. Vol. 63. 066218.

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

университет им. П.Г. Демидова После доработки 02.07.2008

QUALITATIVE ANALYSIS OF ONE CLASS OF OPTOELECTRONIC SYSTEMS SINGULARLY PERTURBED MODELS

D.V. Glazkov

Two models of semiconductor laser with delayed optical feedback are studied. We consider singularly perturbed problem because of the large parameter presence. We construct and discuss quasinormal forms of models in trancritical cases.

Глазков Дмитрий Владимирович - родился в 1982 году в Ачинске Красноярского края. Окончил с отличием математический факультет Ярославского государственного университета им. П.Г. Демидова (2004). В настоящее время является ассистентом кафедры математического моделирования ЯрГУ. Область научных интересов - дифференциальные уравнения с запаздыванием, метод нормальных форм, модели динамики лазеров. Автор 5 статей по данной тематике.

E-mail: glazkov_d@mail.ru