Метод фазового пространства в анализе статистики флуктуаций оптических импульсов в нелинейных диспергирующих средах Текст научной статьи по специальности «Физика»

С.М.Широков

МЕТОД ФАЗОВОГО ПРОСТРАНСТВА В АНАЛИЗЕ СТАТИСТИКИ ФЛУКТУАЦИЙ ОПТИЧЕСКИХ ИМПУЛЬСОВ В НЕЛИНЕЙНЫХ ДИСПЕРГИРУЮЩИХ СРЕДАХ

1. Введение

Процессы распространения сверхкоротких световых импульсов в нелинейных диспергирующих средах привлекают в последние годы возрастающий интерес в связи с возможностями их использования для передачи с особо высокими скоростями, создания элементов оптических компьютеров, новых оптических приборов и т.п. [1-3]. По мере уменьшения длительности импульсов в таких оптических системах все более важное значение приобретает учет флуктуации источника излучения, снижающих степень его когерентности. Знание статистики флукгуаций необходимо для выбора алгоритмов обнаружения, различения и фильтрации сигналов в указанных системах, а также оценки их характеристик [4, 5]. Хотя задачи анализа флуктуирующих сигналов в нелинейных средах исследованы уже довольно подробно, для их функций распределения известны лишь отдельные результаты статистического моделирования на ЭВМ и качественные оценки применительно к некоторым упрощенным моделям среды передачи [3, 6-9]. Ниже предлагается общий подход к решению указанных задач, основанный на представлении модели оптической среды (канала передачи) в гамильтоновой форме относительно некоторых обобщенных фазовых координат.

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

Процессы нелинейного самовоздействия сверхкоротких световых импульсов в нелинейных диспергирующих средах описываются волновыми уравнениями, которые в рамках дисперсионного приближения второго порядка по методу медленно меняющихся амплитуд после соответствующей нормировки обычно допускают представление в виде нелинейного эволюционного уравнения шредингеровского типа (НУШ) [1,3]

ду д1 (р I о /1 ч

/—+а—-+/рф+к0|ф| <р = 0 (1)

Зг[ дх

с начальным условием

Ф(0,г)=ф), (2)

где г[, х - нормированные координаты (пространственная и временная, отсчитываемая от Центра движущегося импульса); ф(л^) - комплексная огибающая импульса излучения, нормированная на начальное пиковое значенеие; j?(t) - входное воздействие; а, (3, к0 -параметры дисперсии, затухания и нелинейности (ко<0 соответствует нормальной, а ко>0

- аномальной дисперсии).

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

С учетом флукгуаций излучения входное воздействие (2) представляется в виде

5(т)=и(т)ц(т)+^) (3)

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

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

3. Метод фазового пространства

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

М/(г1,т)=ф(т1т)^ (4)

относительно которой оно приобретает вид

где

ду У | 12

'тг^тт+^ОМ ч>=о >

£7Г| 01

к(т1)=к0е_2Рт1 .

(5)

(6)

В качестве фазовых координат здесь целесообразно выбрать вещественные функции -квадратурные компоненты х(т1,т), у("п,т) комплексной огибающей

у(л,т)=х(т|,т)+о>(тът) .

(7)

После подставновки (7) в (5) получается система из двух вещественных уравнений, которые можно представить в виде

дх ЪН ду ЪН

дц 5у' а-п Ъх где Н(х,у) - гамильтониан, определяемый выражением

(8)

Н(х

дх

\2 Г ^2 ду

К/ 2 2\2 "Т1* + У )

(к

(9)

Аналогичные комплексные уравнения можно записать, выбрав в качестве фазовых переменных у(г|,т) и у*(п,т) [10].

Полученные уравнения (8), в отличие от аналогичных уравнений, рассматриваемых в статистической механике [11], описывают распределенную динамическую систему с бесконечномерным фазовым пространством и поэтому содержат в правой части не обычные, а вариационные производные (по функциям континуальной переменной т ). При случайном входном воздействии с заданным функционалом плотности вероятности №0{х,у) его эволюция в системе, описываемой гамильтоновыми уравнениями (8), подчиняется уравнению Лиувилля [11], которое в рассматриваемом фазовом пространстве также должно быть записано с использованием вариационных производных:

—+{#, Ф} = 0 , дц

(Ю)

где

ЪН 8(У 8Н 5IV Ъу 8х 5х 5у

(к

(П)

скобки Пуассона.

Уравнение (10) можно получить также из функционального уравнения Фоккера-Планка-Колмогорова [12], положив в нем коэффициенты диффузии равными нулю.

4. Решение с использованием дискретизации по времени

Непосредственное использование уравнений в вариационных производных типа (10) и описываемых ими функционалов плотностей вероятностей для анализа характеристик нелинейных оптических систем затруднительно, так как такие функционалы допускают строгое определение лишь в классе обобщенных функций. При практическом анализе и моделировании на ЭВМ функции х(г|,т) и у{г\,х), как правило, дискретизуются по х с некоторым шагом Дт и представляются векторами отсчетов

где N - число отсчетов,

*к(л)=*(п> (¿-аддт),

Л(Л)=3<1{к-Щ2)Дт); (12)

Полученные в результате такого конечномерного представления случайные векторы X, дописываются обычной 2Лг-мерной плотностью вероятности ЩХ,У,ц), а скобки Пуассона в (10) приобретают вид

дН дГ дН дм>

*=о \дукдхл дхк дук)

(13)

где гамильтониан

После подстановки (14) в (13) и простых преобразований уравнение Лиувилля (10) принимает вид

ттЫ*«-2**! ^(уы-ЪЫ+У*)2

= о (15)

*=о Ц(Дт)

Решая его при начальном условии

»{ХЛЪ)-*^, (16)

ще м0(Х, У) - плотность вероятности квадратурных компонент входного воздействия, методом характеристик, получаем общее выражение 27У-мерной плотности вероятности векторов отсчетов Х(г|), Дт]) в произвольной точке ц :

^ХХЛ) = у>0[Г^Х,Г)] (17)

Здесь /л (Х,У) - оператор нелинейного отображения вход-выход, определяемого уравнением (5) при заданном г\ с учетом дискретизации (12); * - знак

сопряженного оператора. При выводе формулы (17) учтено, что НУШ вида (5) относится к классу уравнений с унитарной нелинейностью [13], поэтому обратный оператор может быть заменен сопряженным.

Как в теоретическом, так и в практическом отношении большой интерес представляет вопрос о существовании стационарных распределений, не изменяющихся при преобразовании (17), а также равновесного распределения, к которому стремится ЦХ, У,л)

при неограниченном увеличении Оц. Согласно выводам статистической механики [11] свойством стационарности обладает любая начальная плотность вида

*0(Х,Г) = ыхт

(18)

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

Ненулевые состояния статистического равновесия в рассматриваемой системе возможны только при отсутствии потерь (/£=0) и характеризуется равномерным распределением по эквиэнергетической поверхности (микроканоническим распределением Гиббса) с заданным значением гамильтониана ЩХ, У)=Н0,

где С - нормирующая константа.

В условиях флуктуаций входного воздействия флуктуирует и Н0 с некоторым распределением ЦЯ0), которое может быть найдено по известному распределению ™о(Х, У) . При этом равновесное распределение следует определять по формуле

В отличие от задач статистической механики, где изменения Н0 обусловлены термодинамическим взаимодействием с внешней средой, в данном случае распределение Ц#0) не обязательно является экспоненциальным и поэтому мр(Х, У) не совпадает с каноническим распределением Гиббса. Для реальных начальных распределений и>о(^,)0 расчет пр(Х,У) описанным методом, как правило, приводит к распределениям, не выражающимся через элементарные функции и допускающим только приближенные или численные представления.

Равновесное распределение, как легко видеть из (19) и (20), определяется только функционалом Н{ХуУ) и не зависит от начальной формы импульса и(х) в (3), т.е. характеризует статистику сигнала на таком большом расстоянии г| , на котором в результате флуктуаций и дисперсионного расширения импульсов содержащаяся в них информация потеряна.

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

Рассмотренный метод обеспечивает определение в явном виде 2Лг-мерной плотности вероятности отсчетов флуктуирующего сигнала в любой точке оптического канала по заданному распределению на входе и поэтому при достаточно большом N применим для анализа как медленных, так и быстрых флуктуаций при любой форме импульса. Это позволяет записать в общем вице отношение правдоподобия и синтезировать оптимальные алгоритмы обнаружения и различения сигналов [4, 5].

Непосредственный численный расчет плотности вероятности по формуле (17) требует решения эволюционного уравнения (5), поэтому его целесообразнно осуществлять на ЭВМ рекуррентно, в соответствии с последовательной процедурой решения указанного уравнения (например, методом расщепления).

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

= СЬ[Н0-Н(Х,У)]

(19)

(20)

амплитуды, скорости и фазы солитонов, как это принято в [6-9]), а их форма остается детерминированной.

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

5. Анализ флуктуации параметров с использованием автомодельных приближений

Уравнения эволюции параметров импульса в замкнутой форме могут быть выведены из НУШ (1) или (5); если принять допущение о сохранении формы импульса в процессе его распространения в среде передачи. Такое допущение справедливо для импульса гауссовской формы в линейном режиме и импульса секанс-гиперболической формы с определенными начальными параметрами в нелинейном односолитонном режиме. В остальных случаях указанные автомодельные решения можно рассматривать как естественные приближения к точному решению НУШ. Автомодельное решение НУШ вида (5) гауссовской формы получено в [14] и определяется выражением

Ч/Длд) = Жл) ехр| - —-—+/¿>(т[)т7 [ (21)

I 2л(т1) J

где Жл)> а(л)» Ь(г|) - параметры импульса (комплексная амплитуда, полуширина и коэффициент линейной частотной модуляции - ЛЧМ), изменяющиеся вдоль волокна.

Аналогичное решение в форме секанс-гиперболических импульсов (солитонов), выведенное вариационным методом в [15] для М мод, при М= 1 можно представить в виде

Ч/Дл,х) = ^(л)5ес/[^]ехр[/г»(т1)т2] . (22)

Как показано в [15], эволюция параметров моделей (21) и (22) описывается одинаковыми по виду уравнениями, отличающимися только числовыми коэффициентами

и Н2 :

(На

— = АсшЬ ;

<Ь\ (23)

4г| \аА ) а3

где О=\А(0)\2а(0) - энергетический параметр.

Для импульсов секанс-гиперболической формы ^4/л2 , \у1=2/ж1 >5 а для импульсов гауссовской формы с тем же значением 0 получается ^=1, ^=1/2^/2; если же

фиксирована энергия импульса, то |Л2=*/8>/2.

Уравнения эволюции параметров (23), как и исходное НУШ (5), могут быть

представлены в гамильтоновой форме:

(24)

дц ' 9л дР

где в качестве фазовых координат выбраны

р = Ь, Я = а2, (25)

а гамильтониан определяется выражением

ца 2ц2к£>

#(/?,*) = 4ар2«?+---—р (26)

9 V?

Соответствующее уравнение Лиувилля для плотности вероятности Мр,я) имеет вид

дм> дм>

-+8арц--

ду дд

ща

4ар-—+—=г Я Я^Я .

^ = 0. (27)

др

Решив его методом характеристик, переходя от я к более удобному для анализа и измерений параметру - интенсивности 1=0/^5 и объединяя оба параметра в вектор 2Н1,Ъ\ можно представить связь распределения вероятностей ЦД^) на произвольном расстоянии т| с начальным распределением м>0(2)=м<7,0) в виде

= ехр| 12а}б (£) | ^ м»0 [ (2)] (28)

где 1\2) - нелинейное отображение в двумерном пространстве параметров, определяемое уравнениями (23) с учетом указанной подстановки, Р{1{2) - ьая ветвь обратного отображения.

В отличие от 2Лг-мерного распределения (17), двумерная плотность, определяемая (28), может быть сравнительно просто рассчитана на ЭВМ. Процедуре расчета целесообразно придать рекуррентный характер, соответствующий выбранному алгоритму решения уравнений эволюции параметров (23). При этом плотность вероятности вектора отсчетов параметров 2^=[/(ААп), ¿К&Дг])] в каждой последующей точке определяется по плотности в предшествующей точке согласно формуле

((А+иДл 1

12а ^¿(ОкИО] <29>

где - отображение

в , определяемое формулами Рунге-Кутга или другими аналогичными разностными соотношениями; £=0,1,2,...

6. Результаты моделирования и выводы

Для рассмотренных методов анализа распределений вероятностей были разработаны соответствующие алгоритмы моделирования на ЭВМ. На рис.1 и 2 в качестве примеров показаны полученные диаграммы эволюции плотностей вероятностей Ц/,п) пиковой интенсивности импульсов секанс-Гийерболической формы для двух значений параметра нелинейности: к=2 и к=4 "соответственно. Начальное распределение интенсивности принято экспоненциальным, что соответствует гауссовской модели флуктуации квадратурных компонент, а для коэффициента ЛЧМ - равномерным, при постоянной энергии всех реализаций импульса. Значения интенсивности / на диаграммах нормированы на математическое ожидание т1.

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

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

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

Непосредственно из полученных диаграмм видно, что в рассматриваемых условиях интенсивность импульса в нелинейной диспергирующей среде флуктуирует по вероятностному закону, явно отличающемуся от экспоненциального, т.е. статистика Флуктуации существенно отличается от гауссовской. Проверка соответствия Распределений квадратурных компонент гауссовскому закону по критерию хи-квадрат Дает следующие результаты: для флукгуаций на входе получено х2=7,74, на выходе, при 11=0,5 и к=4 х2=615. При этом уровню значимости а=0,1 соответствует х2а=24,1. Это

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

Аналогичное моделирование для случая, когда нет ограничения, касающегося постоянства энергий всех реализаций, показывает, что деформация начального экспоненциального распределения при этом происходит значительно медленнее и без образования вторичных пиков. При тех же условиях на дальности 11=0,5 получено Х2=25,1, т.е. распределение флукгуаций, хотя и отличается от гауссовского, но не столь значительно, как в первом случае.

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

Рассмотренные методы расчета распределений вероятностей использованы при анализе помехоустойчивости, дальности, скорости передачи и других характеристик солитонных систем оптической связи и позволяют уточнить приведенные в [4, 5] их оценки.

ЛИТЕРАТУРА

1. Сисакян И.Н., Шварцбург А.Б.//Квантовая электроника, 1984, т.11. N 9, стр.1703-1721.

2. Зуев MA., Колчанов И.Г., Шварцбург А2>.//Компьютерная оптика: Сб.МЦНТИ, М.,1992, вып.12, стр.56-60.

3. Ахманов CA., Выслоух ВА., Чиркин A.C. Оптика фемтосекундных лазерных импульсов. M.: Наука, 1988.

4. Широков С.Л/у/Электросвязь. 1992, N 3, стр.3-8.

5. Широков С.АГ//Компьютерная оптика: Сб.МЦНТИ, М., 1993, вып. 13, стр.59-64.

6. Bava G.P., Chione G., Maio /.//Electr. Letters, 1984, V.20, No.24, pp.1002-1003.

7. Кандидов В.П., Шленов СЛ.//Известия вузов. Радиофизика, 1984, т.27, N 9, стр.1158--1167.

8. Иванов A.B., Матвеев Л.#.//Вестник МГУ, 1987, т.28, N 5, стр.82-85.

9. Выслоух ВА., Иванов A.B., Чередник И.В.ЦИзвестия вузов. Радиофизика, 1987, т.30, N 8, стр.980-990.

Ю.Тахтаджян JIA., Фаддеев ЛД. Гамильтонов подход в теории солигонов. М.: Наука, 1986. U.Боголюбов H.H., Боголюбов H.H. (мл.) Введение в квантовую статистическую механику. М.: Наука, 1984.

И.Гардинер К.В. Стохастические методы в естественных науках. М.: Мир, 1986. ХЪ.Маслов В.П. Комплексные марковские цепи и континуальный интеграл Фейнмана. М.: Наука, 1976.

^Anderson DJ/ Phys.RevA: Gen.Phys., 1983, v.27, No.6, pp.3135-3145. \5.Shirokov £Af.//Proc.ISFOC-93: Boston-St.Petersburg, 1993, pp.206-211.