РОССИЙСКАЯ АКАДЕМИЯ НАУК

СИБИРСКОЕ ОТДЕЛЕНИЕ

ИНСТИТУТ ВЫЧИСЛИТЕЛЬНЫХ ТЕХНОЛОГИЙ

На правах рукописи

Штырина Ольга Владимировна

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ СОЛИТОННЫХ

ОПТОВОЛОКОННЫХ ЛИНИЙ СВЯЗИ

05.13.18 Математическое моделирование,

численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Новосибирск 2006

Работа выполнена в Институте вычислительных технологий Сибирского отделения РАН

Научный руководитель: доктор физико-математических наук, Михаил Петрович Федорук

Официальные оппоненты: доктор физико-математических наук, профессор Александр Михайлович Блохин кандидат физико-математических наук, доцент Елена Геннадьевна Шапиро

Ведущая организация: Новосибирский государственный университет, г. Новосибирск.

Защита состоится 8 сентября 2006 года в 1500 часов на заседании диссертационного совета Д 003.046.01 при Институте вычислительных технологий СО РАН по адресу:

630090, г. Новосибирск, проспект Академика Лаврентьева, 6, конференцзал ИВТ СО РАН.

С диссертацией можно ознакомится в специализированном читальном зале вычислительной математики и информатики научного отделения СО ГПНТБ (проспект академика М.А.Лаврентьева, 6).

Автореферат разослан августа 2006 г.

Ученый секретарь диссертационного совета доктор физико-математических наук, профессор Л. Б. Чубаров

Общая характеристика работы

Актуальность темы. Беспрецендентное увеличение пропускной способности волоконно-оптических линий связи обусловлено увеличением спроса на телекоммуникационные услуги. Действительно за последние 25 лет суммарная скорость передачи информации возрасла на пять порядков и достигла величины 1 Тбит/с. Однако, даже такие скорости передачи информации не смогут удовлетворить все возрастающие потребности пользователей Интернета, поскольку их число непрерывно растет. Так, если в 1998 году в мире было около 25 млн. человек пользователей Интернета, то в конце 2005 г. число пользователей составило уже 970 млн. человек. Поэтому весьма актуальной является задача дальнейшего увеличения пропускной способности оптических линий связи, которое возможно либо путем расширения спектрального диапазона передачи данных, либо путем увеличения спектральной эффективности передачи данных, которая определяется как отношение скорости передачи информации в одном частотном канале к расстоянию между соседними частотными каналами.

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

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

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

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

Разработан новый эффективный численный алгоритм для построения солитонных решений в рамках усредненного обобщенного уравнения Шредингера. Найдены новые солитонные решения для важных в практическом отношении дисперсионных конфигураций волоконных линий связи.

Предложена конфигурация оптического 2R регенератора сигналов на основе полупроводникового насыщающегося поглотителя. Определены параметры конкретных конфигураций многоканальных оптоволоконных линий связи с дисперсионным управлением и периодически встроенными в них оптическими регенераторами, и установлено, что оптические сигналы в подобных линиях связи способны передаваться на расстояния порядка 8-10 тысяч километров с пропускной способностью 40 Гбит/сек в одном частотном канале и спектральной эффективностью 0.2 бит/Гц/с.

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

Основные положения выносимые на защиту:

• новый эффективный численный алгоритм решения обобщенного интегрального уравнения Шредингера с использованием рекуррентного метода вычисления интегрального оператора;

• оригинальный экономичный численный алгоритм для решения нелинейного усредненного уравнения Шредингера в спектральной области, основанный на теореме о свертке и методе стабилизирующей поправки;

• новые солитонные решения нелинейного усредненного уравнения Шредингера в спектральной области для практически важных конфигураций волоконно-оптических линий связи с дисперсионным управлением (распределенной дисперсией);

• методология эффективного уширения непрерывной волны накачки в широкополосных Рамановских усилителях на основе эффекта модуляционной неустойчивости;

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

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

1. CLEO Europe-2005, 2005, Munich, Germany.

2. International Conference on Computational Mathematics, 2004, Novosibirsk, Russia;

3. Nonlinear Guided Waves and Their Applications, 2004, Toronto, Canada;

4. МНСК "Интеллектуальный потенциал Сибири", 2003 г., Новосибирск, Россия;

5. International Conference on Computational Mathematics, 2002, Novosibirsk, Russia;

6. МНСК "Студент и научно-технический прогресс", 2001-2003 гг., Новосибирск, Россия;

Материалы диссертационной работы использовались при выполненении грантов РФФИ (гранты № 99-02-16688, № 03-02-16496 и № 06-01междисциплинарного интеграционного проекта СО РАН № 31, грантов Министерства образования и науки РФ (проекты № 1.80.01 и № 1.13.06) и инновационного гранта INTAS-03-56-203.

Публикации. Основные результаты диссертации опубликованы в работах [1-12]. При выполнении работ [1-8], опубликованных совместно с научным руководителем и другими соавторами, О.В. Штырина принимала участие в постановке задачи, разработке численных алгоритмов, обсуждении полученных результатов, подготовке и представлении статей и докладов на конференциях. Ею выполнена программная реализация разработанных численных алгоритмов, проведены расчеты тестовых задач и значительный цикл вычислительных экспериментов. Кроме того, Штыриной О.В. в [7,8] проведена серия расчетов, а также сравнение и анализ полученных результатов с экспериментальными данными.

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения и списка литературы.

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

Глава 1 состоит из 4 параграфов, в которых формулируется постановка задачи, обсуждаются математические модели и алгоритмы их численной реализации излагаются результаты численного моделирования по исследованию процесса уширения непрерывной волны накачки в волоконных световодах вследствие развития модуляционной неустойчивости CW-волны. На основе этих результатов формулируются основные правила, необходимые для эффективного уширения непрерывной волны накачки в широкополосных Рамановских усилителях на основе эффекта модуляционной неустойчивости.

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

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

Обобщенное интегральное нелинейное уравнение Шредингера для комплексной огибающей A амплитуды электромагнитного поля имеет вид:

параметр дисперсии групповой скорости, 3 дисперсия третьего порядка, коэффициент керровской нелинейности, эффективный коэффициент, учитывающий затухание и усиление сигнала. Величины 2, 3, и представлены как функции от z, чтобы учесть изменения этих параметров при переходе от одного типа световода к другоn му. Коэффициент нелинейности определяется формулой =, где n2 нелинейный показатель преломления, 0 несущая длина волны, 0 = радиальная частота несущего сигнала, cl скорость света, Aef f эффективная площадь собственной моды световода, R (t) = (1 fR ) (t) + fR hR (t) является функцией запаздывающего отклика среды. Здесь (t) дельта-функция, и fR парциальный вклад рамановского отклика. Для функции рамановского отклика hR использовалось приближение демпфированного осциллятора [1]:

где 1 =12.2 фс, 2 =32 фс.

Для сигналов c длительностями TF W HM 10 фс, уравнение (1) может быть существенно упрощено. В этом случае подынтегральную функцию |A(z, t t )|2 можно разложить в ряд Тейлора:

Определим первый момент нелинейной функции отклика как и заметим, что R(t)dt = 1, тогда уравнение (1) сводится к обобщенному нелинейному уравнению Шредингера:

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

При решении уравнения (1) основная трудность состоит в построении эффективного численного алгоритма вычисления интегрального оператора Таким образом, используя аналитическое выражение (2) для функции hR мы приходим к вычислению интеграла Сделаем замену переменных s = t t, тогда интеграл (4) перепишется в следующем виде Далее легко показать, что Вычисление интегралов Ni (A, t) проводится по следующей схеме где tj = jh, j [ N + 1,..., N ], tj [a, a].

Легко видеть, что для вычисления интегрального оператора (4) потребуется O(N ) операций, где N число точек на интервале интегрирования.

В четвертом параграфе на основе вычислительного алгоритма, предложенного в третьем параграфе, выполнено численное моделирование процесса уширения непрерывной волны накачки за счет развития модуляционной неустойчивости. Для существенного расширения спектрального диапазона передачи данных и увеличения числа информационных каналов необходимо создать волоконный рамановский усилитель с равномерной полосой усиления порядка нескольких сотен нанометров [2]. Для эффективного уширения и выравнивания полосы усиления непрерывной волны накачки (CW-волны) предлагается использовать эффект модуляционной неустойчивости. В этом параграфе представлено детальное исследование процесса нелинейного уширения рамановской волны накачки, инициированного развитием модуляционной неустойчивости. Выполнена серия расчетов по оптимизации параметров уширения в волоконном Truewave (TW) световоде. На рисунке 1 показана зависимость периода колебаний от входной мощности накачки. Сплошная кривая результаты численных расчетов, пунктирная соответствует теоретическому закону убывания периода колебаний T Pin. Рисунок 2 иллюстрирует полученную численно зависимость ширины спектра от хроматической дисперсии на выходе из TW-световода длиной 25 км и входной мощностью накачки Pin = 0.9 Вт, а на рисунке 3 показан характерный вид спектра уширенной накачки для величины хроматической дисперсии D = 0.03 пс/нм/км. Изолинии спектральной ширины в плоскости (длина световода, дисперсия) представлены на рисунке для выходной мощности накачки Pout = 0.25 Вт.

Глава 2 также состоит из 4 параграфов.

В первом параграфе дается краткий обзор литературы и описываются основные свойства, так называемых, дисперсионно-управляемых солитонов (ДУ-солитонов) в волоконно-оптических линиях связи с периодически меняющейся дисперсией.

Солитоны с дисперсионным управлением являются новым типом оптических носителей информации. Численное моделирование и данные лабораторных экспериментов позволили выявить ряд особых свойств ДУ-солитонов, существенно отличающих их от фундаментального солитона нелинейного уравнения Шредингера с постоянными коэффициентами (см., например, [3]).

Во втором параграфе формулируется математическая модель для построения ДУ-солитонных решений в рамках усредненного нелинейного уравнения Шредингера в спектральной области. Рассмотрим нелиOscillation period (ps) нейное уравнение Шредингера с периодическими коэффициентами здесь где g(z) эффективный коэффициент, описывающий затухание и усиление оптических сигналов. В предположении, что нелинейная длина LN L много больше дисперсионной длины LD, уравнение (5) можно усреднить и получить следующее уравнение [4]:

средняя дисперсия, G(, ) есть нелинейный интегральный оператор:

Averaged spectral power (dB) с матричным элементом T123, который является функцией величины Функция R(z) определяется из обыкновенного дифференциального уравнения Rz = d(z) d. Уравнение (6) имеет стандартную форму типичную для моделей, описывающих четырехволновое взаимодействие с квадратичным законом дисперсии. Поэтому специфические свойства модели определяются зависимостью матричного элемента T от [5].

Будем искать солитонное решение уравнения (6) в следующей форме:

Тогда уравнение, описывающее форму ДУ-солитона () принимает вид:

В третьем параграфе предложен численный алгоритм реализации уравнения (7), основанный на теореме о свертке и методе стабилизирующей поправки Петвиашвили. Идея метода состоит в аппроксимации матричного элемента T () соответствующим набором функций. Такая аппроксимация позволяет применить алгоритм быстрого вычисления свертки и уменьшить число операций до величины порядка M N log2 (N ), где M зависит от конкретной аппроксимации T ().

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

где Rn некоторые коэффициенты.

Способ получения тригонометрического разложения заключается в нахождении коэффициентов ряда Фурье функции T () на интервале интегрирования. Матричный элемент T () имеет очевидную симметрию, поэтому интегральный оператор представляется в симметричной форме. Интегрируя по 1, и, вводя новые симметричные переменные 2 = 1 + и 3 = 2 +, получим интеграл:

где = 2 +(+ 1 + 2 )2 (+ 1 )2 (+ 2 )2. Факторизуя матричный элемент, получим следующую симметричную форму, которая удобна для быстрых вычислений:

Для вычисления этого интеграла можно использовать быстрый алгоритм для свертки и корреляции.

В четвертом параграфе построены конкретные примеры солитонных решений уравнения (7) для случаев мало-масштабного и крупномасштабного дисперсионных управлений. Здесь мы используем выражения "крупно-масштабные"или "мало-масштабные", если длина дисперсионной компенсации соответственно больше или меньше расстояния между оптическими усилителями.

Мы ограничимся здесь только одним характерным примером построения солитонных решений для случая мало-масштабного дисперсионного управления. Распределение мощности ДУ-солитона (сплошная линия) и фундаментального солитона с равной амплитудой (прерывистая линия) представлены на рисунке 5 при различных J. Здесь J =, где Za расстояние между усилителями, L длина компенсации дисперсии. Заметим, что при больших J, форма ДУ-солитона близка к форме фундаментального солитона. Этот результат полностью согласуется с теорией, развитой в [6].

Power Power Глава 3 состоит из 4 параграфов.

В первом параграфе описаны основные способы оптической регенерации сигналов и сформулирована задача оптимизации волоконнооптических линий связи в терминах коэффициента ошибки (Bit-Error Rate). Известно, что наиболее экономичным способом увеличения спектральной эффективности передачи данных в многоканальных линиях оптической связи является увеличение скорости передачи данных одного частотного каналам. В работе [7] отмечается, что увеличение скорости передачи информации в четыре раза обеспечивает снижение стоимости единицы передаваемой информации примерно в 2.5 раза.

К числу наиболее перспективных способов дальнейшего увеличения пропускной способности одного частотного канала до скоростей передачи 40 Гбит/c и более относится технология оптической регенерации сигналов.

Во втором параграфе приведена принципиальная схема 2R оптического регенератора и сформулирована математическая модель полупроводникового насыщающегося поглотителя, который является основным элементом данного устройства. Принцип работы насыщающегося поглотителя (SA) заключается в поглощении мощности входящего в него оптического сигнала, если она окажется ниже некоторой пороговой мощности Psat. При мощностях, больших Psat, коэффициент пропускания SA быстро возрастает и асимптотически приближается к единице. В таких условиях маломощное излучение усиленного спонтанного шума и фоновое дисперсионное излучение подавляются SA. Использование SA в сочетании с узкополосным оптическим фильтром и сильно нелинейным волоконным световодом (HNF) позволяет подавить шумы в единичных битах. Функция потерь (t) в насыщающемся поглотителе (SA), зависящая от времени и мощности входного сигнала, описывается следующим уравнением:

где P (z, t) = |A (z, t)| распределение мощности входного сигнала, 0 = 3 дБ постоянные потери, z const фиксированное расстояние, Psat является пороговой мощностью насыщения, и соответствует времени спада импульса.

Тогда передаточная функция T (t) = 1 (t, P (z, t)), и действие SA на сигнал описывается следующим образом:

Здесь Pin и Pout соответственно входная и выходная мощности оптического сигнала.

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

Четвертый параграф посвящен автосолитонным режимам распространения одиночных оптических импульсов в симметричной волоконно-оптической линии связи. Здесь же рассмотрены примеры оптимизации двух практических многоканальных оптических линий связи с симметричными конфигурациями со скоростью передачи данных B = 40 Гбит/с в одном частотном канале. Распространение оптических сигналов вдоль волоконных световодов определяется уравнением (3).

Результаты численного моделирования показали, что любой входной импульс (в пределах некоторой области значений параметров) эволюционирует в устойчивое асимптотическое решение (оптический автосолитон). На рисунке 6 показан фазовый портрет сигнала в плоскости (параметр фазовой модуляции, ширина сигнала). Рисунок 7 представляет эволюцию входного гауссовского импульса на выходе из оптического регенератора.

C -0. Для оценки коэффициента ошибки обычно используется непрямой статистический метод оценки концепция Q-фактора оптической линии связи [8]. В случае бинарного ООК-формата "включено-выключено"в предположении, что статистика нулевых и единичных битов подчиняется нормальному закону распределения, величина Q-фактора системы равна где µi и i при (i = 0, 1) математическое ожидание и среднеквадратичное отклонение для нулевых и единичных битов соответственно.

Тогда здесь erf c модифицированная функция ошибок. Отметим, что коэффициенту ошибки BER 109 соответствует значение Q 6. Мы использовали в роли критерия качества передачи сигнала величину Qфактора, и дистанцию передачи данных определяли как расстояние, для которого величина Q 6. Для вычисления дистанции передачи данных в расчетах, представленных ниже, использовались псевдослучайные последовательности с длинами в 210 1 214 1 бита, и дистанция передачи данных вычислялась как расстояние, при котором Q 6. Далее дистанция передачи данных выбиралась как наименьшее расстояние по всем частотным каналам.

На рисунке 8 представлен характерный пример результатов многочисленных численных расчетов. В расчетах рассматривалось 4 частотных канала с расстоянием = 200 ГГц между соседними каналами.

На этом рисунке показан результат оптимизации линии PSCF + RDF+PSCF+EDFA в плоскости параметров (средняя дисперсия линии, входная пиковая мощность). Параметры соответствующих волоконных световодов приведены в тексте диссертации. Видно, что соответствующим выбором параметров оптического регенератора, входной пиковой мощности импульсов и средней дисперсии линии можно достичь дистанции передачи данных свыше 10000 км. Данная система без оптических регенераторов демонстрирует дистанцию распространения около км.

Оптимизационные расчеты были выполнены также для линии TL + RTL + TL + EDFA. Здесь при оптимальных параметрах системы удалось достичь дистанции распространения более 8000 км.

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

Список опубликованных работ по теме диссертации 1. Штырина О.В., Турицын С. К., Федорук М.П. Исследование волоконно-оптических линий связи с оптической регенерацией сигналов // Квантовая электроника.–2005.–Т.35.–№ 2.–С.169-174.

2. Федорук М.П., Штырина О.В. Математическое моделирование и оптимизация волоконно-оптических линий связи с дисперсионным управлением и оптической регенерацией сигналов // Вычислительные технологии.–2004.–Т.9.–С.150-158.

3. Fedoruk M.P., Shtyrina O.V., Turitsyn S.K. Autosolitons in dispersionmanaged systems with in-line saturable absorbers // Optics Letters.– 2004.–V.29.–№21.–P.2464-2466.

4. Shtyrina O.V., Fedoruk M.P. Numerical simulation and optimization of bre links with optical regenerators // Труды Международной конференции по вычислительной математике. – Новосибирск (МКМВЧ.2.–С.932-935.

5. O.V.Shtyrina, S.B.Medvedev and M.P.Fedoruk Dispersion-managed soliton for path-averaged model of optical ber communication line // Proceedings of the International Conference on Computational Mathematics.–2002.–Part 2.–P.697-703.

6. Medvedev S.B., Shtyrina O.V., Musher S.L., Fedoruk M.P. Pathaveraged optical solitons in double-periodic dispersion-managed systems // Physical Review E.–2002.–V.66.–P.066607-1-0.66607-6.

7. Ellingham T.J., Fedoruk M.P., Shtyrina O.V., Ania-Castanon J.D., Turitsyn S.K. Design of a ber scheme for nonlinear Raman pump broadening through modulation instability // Conference Digest for CLEO Europe.–2005, Munich, Germany.–P.CJ-1-WED.

8. Ellingham T.J., Ania-Castanon J.D., Shtyrina O.V., Fedoruk M.P., Turitsyn S.K. CW Raman pump broadening using modulational instability // Nonlinear Guided Waves and Their Applications.–2004.– None.–MC42.

9. Штырина О.В. Солитоны с управляемой дисперсией для усредненной модели в оптоволоконных линиях связи. // Материалы Международной научной студенческой конференции "Интеллектуальный потенциал Сибири". Новосибирск.–2003.–С.13.

10. Штырина О.В. Математическое моделирование генерации суперконтинуума в волоконных световодах // Материалы XLI Международной научной студенческой конференции "Студент и научнотехнический прогресс", секция "Математика". Новосибирск.– 2003.–С.128.

11. Штырина О.В. Солитоны с управляемой дисперсией для усредненной модели в оптоволоконных линиях связи. // Материалы XL Международной научной студенческой конференции "Студент и научно-технический прогресс", секция "Математика".

Новосибирск.–2002.

12. Штырина О.В. Численное моделирование оптоволоконных линий связи в средах с периодически меняющейся дисперсией // Материалы XXXIX Международной научной студенческой конференции "Студент и научно-технический прогресс", секция "Математика".

Новосибирск.–2001.–С.106.

Литература [1] Agrawal G.P. Nonlinear Fiber Optics. New York: Academic Press, 2001.

[2] Дианов Е.М. От тера-эры к пета-эре // Вестник РАН – 2000 – Т. – №11.– С.1010-1015.

[3] Turitsyn S.K., Shapiro E.G., Medvedev S.B., Fedoruk M.P., Mezentsev V.K. Physics and mathematics of dispersion-managed optical solitons // Comptes Rendus. Physique.–2003.–V.4.–P.145-161.

[4] Gabitov I., Turitsyn S.K. Averaged pulse dynamics in a cascaded transmission system with passive dispersion compensation // Opt.

Lett.– 1996.–V. 21.–P.327-329.

[5] Medvedev S.B., Turitsyn S.K. Hamiltonian averaging and integrability in nonlinear systems with periodically varying dispersion // JETP Lett.–1999.–V.69.– P. 499-506.

[6] Turitsyn S.K., Turitsyna E.G., Medvedev S.B., Fedoruk M.P.

Averaged model and integrable limits in nonlinear double-periodic Hamiltonian systems // Phys. Rev. E.–2000.–V. 61.– P.3127-3132.

[7] Величко М.А., Наний О.Е., Сусьян А.А. Новые форматы модуляции в оптических системах связи // Lightwave Russian Edition.–2005.– №4.–С.21-30.

[8] Agrawal G.P. Fiber-Optic Communication Systems. Second edition.

New York: John Wiley & Sons, Inc., 1997.

Штырина Ольга Владимировна

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ СОЛИТОННЫХ

ОПТОВОЛОКОННЫХ ЛИНИЙ СВЯЗИ

диссертации на соискание ученой степени кандидата физико-математических наук 630090, г. Новосибирск, пр. Ак. Лаврентьева, 6. тел. 330-72-02.