Формирование и передача на расстояние изображений с помощью мод Гаусса-Лагерра

ФОРМИРОВАНИЕ И ПЕРЕДАЧА НА РАССТОЯНИЕ ИЗОБРАЖЕНИЙ С ПОМОЩЬЮ МОД ГАУССА-ЛАГЕРРА

С.Н. Хонина

Институт систем обработки изображений Российской Академии Наук

Аннотация

Чтобы передать на расстояние изображение через свободное пространство или по волокну нужно сформировать инвариантное изображение. Под инвариантным изображением мы понимаем поперечное распределение интенсивности светового поля, которое сохраняется с точностью до масштаба и поворота при распространении в пространстве. Это можно сделать с помощью аппроксимации заданного изображения набором мод Гаусса-Лагерра (ГЛ) со специально выбранными номерами. Анализируются условия, которым должны удовлетворять номера мод ГЛ, чтобы формировать стабильные или вращающиеся изображения, а также изображения, сохраняющие вид в Фурье-плоскости,

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

Введение

В последнее время возрос интерес к световым пучкам, сохраняющим поперечное сечение при распространении в некоторой среде. Сформировать такие пучки можно с помощью дифракционных оптических элементов (ДОЭ), согласованных с модами Бесселя [1-3], Гаусса-Лагерра [4,5], Гаусса-Эрмита [4,6]. Определенные композиции этих мод позволяют получать стабильные, периодические или вращающиеся световые пучки [7-10]. Стабильные и вращающиеся (можно считать их стабильными с точностью до поворота) композиции значительно расширяют множество пучков с инвариантными свойствами, позволяя формировать более сложные конфигурации поперечного сечения.

В данной работе предложен алгоритм формирования заданных изображений с помощью аппроксимации его модами Гаусса-Лагерра. Показано, что накладывая определенные ограничения на участвующие в аппроксимации номера мод, можно получить изображения, сохраняющиеся при распространении с точностью до масштаба и поворота. Подробно рассмотрены три группы мод, удовлетворяющие условию вращения светового пучка как единого целого [9,11], условию равенства скоростей [12] и условию равенства собственных чисел Фурье-преобразования [5].

В [14] было предложено формировать вращающиеся изображения с помощью "элементарных" спиральных гауссовых пучков Такие изображения представляют собой контурный рисунок, выложенный с помощью таких сдвинутых "элементарных" пучков, и могут вращаться только с одной скоростью, поворачиваясь за все расстояние от нуля до бесконечности на 90 градусов. Предлагаемый в данной работе подход позволяет формировать

изображения с принципиально любой задаваемой скоростью вращения

1. Формирование изображений с помощью мод Гаусса-Лагерра

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

Я= 0 |я £it

(О

где Lllm

г2 = X +у , (р = 0ГС/£ — X

где Г„(х) = (-1)п—-[£„,„(*)] - обобщенный мно-

с1х

гочлен Лагерра, а - параметр, характеризующий эффективную ширину модовой функции (2). Функции Гаусса-Лагерра ортогональны

2,то

о о

\\/1 = т

где опт = < - символ Кронекера, * - знак

комплексного сопряжения.

Из условия (3) следует выражение для вычисления коэффициентов С^ суммы (1):

Сгт = Ц^^г,<р)1*пт(г,<р)Шгй<р

Аппроксимировать объектную функцию Дг,ср) с помощью суперпозиции мод (г. <р), определенной уравнением (1), можно с помощью минимизации следующего критерия:

М = \&г&<р

о о

Минимизировать критерий можно методом последовательных приближений аналогично итеративному алгоритму сокращения ошибки [14],

В данном случае итеративный алгоритм принимает следующий вид:

N

1

п-йтйп

т\ ехр('а»т (г> <Р) г>

= агя| ¡]Г {г, (р) ехр[/Г(р) (г, р)^ (г, <р)гйгй(р

где р - номер итерации, Т{р\г,(р) - фаза комплексной аппроксимирующей модовой функции (1), Вш > О - положительные произвольно заданные числа, характеризующие вес отдельной моды в смеси мод, агб(-} - аргумент комплексного числа.

2. Условие равенства скоростей

Сформированное таким образом изображение (1) будет распространяться в свободном пространстве, подчиняясь уравнению Гельмгольца [15]:

(г, <р, г) = ехр

N

*ЕЕ

На +

¡кг

2Я

Н лН ' 2гг \

к^У)

где (г) = (2и+ю+1)-агй^/г0),

ехр [-фш

Я(г ) = г

( гЛ

1+Ц

V г

- радиус кривизны параболичес-

кого фронта светового поля, сг(г) = с

7,

эффективный радиус пучка, г0 = — - конфокаль-

ный параметр или Френелевское расстояние, к -волновое число света с длиной волны X.

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

(

преломления пг{г)-п\

2 Л

Г

1-2Д-Г-го /

, где п0 - макси-

мальное значение показателя преломления на оси волокна, г0 • радиус волокна, А - параметр дисперсии показателя преломления, изображение (1) будет распространяется следующим образом [16]:

л-0

"•/Г

а

\ /

„2

V

хехр

'У*,*--

и

где а -

ч * "о

Упт ~ 2/1о ~ 4ог_2 (2и+/л +1)| 2 - коэффициент, пропорциональный фазовой скорости соответствующей моды.

В обоих случаях при соблюдении условия равенства скоростей:

П, :2л+\т\ = р,/> = 0,1,2,-

зависимость от г выносится за знак суммы и мы получим стабильный пучок с сохраняющимся поперечным сечением, представляющим собой несколько другую, чем (1) аппроксимацию заданного изображения F(r,(p):

^

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

На рис.1 показаны некоторые группы индексов (п,т) е£1, , удовлетворяющие условию равенства скоростей (11). Подходящие пары индексов легко находятся по пересечению линий с узлами целочисленной сетки. Так, скорость р= 10 имеют 3 моды: (4,-2), (5,0), (4,2); скорость р=13 имеют 4 моды: (5,-3), (6,-1), (6,1), (5,3); скорость р=24 имеют 9 мод: (8,-8), (9,-6), (10,-4), (11,-2), (12,0), (11,2), (10,4), (9,6), (8,8). Ясно, что чем больше значение скорости, тем большее количество мод может участвовать в аппроксимации (12).

2п+|т=24

т: -13 -12 -11-10 -9 -8 -7 -6 -5 -4 -3 -2

2 3 4

6 7

9 10 11 12 13

3. Условие вращения

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

Условие вращения многомодовых пучков Гаусса-Лагерра независимым способом и в различном виде были получены в [9-11]. В [10] рассматривалось световое волокно, а в [11] - свободное пространство. При определенных допущениях условие вращения в обоих случаях выглядит одинаково:

2(л-л') + И~1

ш\

= В.

т-т

где В - постоянная величина.

Это условие легко сводится к виду, полученному в [9]:

П2:2я+|/я|+7>я = д,

где Т^-В - величина, определяющая угловую скорость и направление вращения; при 7>0 вращение по часовой стрелке, при Т<0 вращение против часовой стрелки, при Т=0 вращение не происходит и условие (14) сводится к условию (11). При этом нужно заметить, что в отличие от целого положительного р в (11), постоянные величины Вид являются в общем случае действительными.

Аппроксимация заданного изображения /-"(^ф) в этом случае выглядит следу ющим образом:

I С1Кт{г,<Р).

В [9] дана достаточно сложная формула подбора индексов пары мод исходя из равенства

2 п + \гп\ + Тт = 2 п + \т'\ + Тт

и конкретного значения Т=Р/{2, представляемого в виде несократимой дроби. Однако, нужно заметить,

что композиция из двух любых мод всегда формирует либо вращающийся, либо стабильный пучок. Данное утверждение вытекает из представления условия вращения в виде (13). Это условие выполнено для любых двух пар индексов (л,ш), вращение не будет происходить при ш=т' (В=оо), или при 2п +|/к| = 2л' + |т'| (6=0), что соответствует

условию равенства скоростей (11).

Чтобы определить значения индексов для более, чем двух мод, предлагается решать уравнение (14) графически. Перепишем (14) в виде:

(2п + (Т + 1)т = д,т>0,

[2п+(Т-\)т = д,т< 0.

На рис.2 показано графическое представление (17) для Т>0, ¡7=16. При Т<0 картина будет симметричной относительно оси п. Аппроксимация (15) для одинакового абсолютного значения |7] будет выглядеть одинаково, но для разных знаков будет различное направление вращения. Значение д/2 соответствует положению узла линий с различным 7 на оси п. При д<0 будет существовать только нижний по отношению к узлу пу чок линий.

Как видно из рис.2, при прочих равных-условиях, начиная с />1 подмножество 02 шире подмножества <Г2,. Действительно, для д=\6 П, включает 5 мод (7—0), а может включать 12 мод (7М) или 14 мод (7-3). Это должно дать лучшее качество аппроксимации (15) по сравнению с аппроксимацией (12). Можно сделать предположение, что подмножества с 7-1 и 7=3 наиболее перспективные, так как включают в себя диагонали и линии, параллельные осям п и т.

Условие вращения (17) можно свести к более простому виду

п = р\т\ + 1,

где все т одного знака, 0 и / - целые. Параметр Т в этом случае принимает только целые нечетные значения и определяется по формуле:

Т = -{2р+\)*&

Рис. 2

п:

Понятно, что условие (18) охватывает меньшее подмножество, чем условие (17). В частности, если сравнить рис. 2 и 3 для Т>\, то видно, что условие (18) включает только левые от оси п ветви.

Заметим, что при выборе всех т=0 параметр угловой скорости Т=ос, что легко видеть из формулы (13). Пучки, состоящие из композиции более двух таких мод не обладают свойством инвариантности поперечного сечения при распространении. Эти моды не имеют углового момента (от=0) и скорости распространения у них различны - в (11) 2л^сомл*.

, 1,/п > 0

где 5§п(ш) = < - таковая функция.

[- \,т < 0

На рис.3 показано графическое представление (18) для некоторых Т при 1=5.

Ш

Л

4. Условие равенства собственных чисел

Кроме того, что функции Гаусса-Лагерра удовлетворяют уравнениям распространения (8) и (10), они являются собственными функциями преобразования Фурье:

— I ]ьпя, (г,97)ехр[- \гр СО8(0 - V)] г с1 г&0 =

0 0

что позволяет использовать их для эффективного ввода излучения в световые волокна [16].

В общем случае, если мы сформировали световой пучок (г, (р), который является суперпозицией мод Гаусса-Лагерра и одновременно аппроксимацией заданного изображения Яг,<р), то в зоне дифракции Фраунгофера или в фокальной плоскости линзы мы получим следующее изображение:

= £ ХЛшС-ЛЛа*).

П- 0

где 3 - обозначение Фурье-преобразования, Л«, = (_"1)и(—0 - собственные числа.

Из уравнения (21) видно, что функция Р£(Г><Р) в общем случае не равна своему Фурье-

образу . Однако, при участии в

суперпозиции (1) только мод с одинаковыми значениями собственных чисел:

_ ,2я-Э|п

= const,

Q3:2и + 3|m| = /+4i,/,s = 0,1,2,...,

для интснсивностеи световых пучков такое равенство может бьпъ достигнуто:

2

|з{р* (r,p)j2 =

У^яд-^щя

С Х^

п.теП)

п.етеО з 2

= F,

"з

Более наглядно распределение одинаковых значений собственных чисел в зависимости от номера функции Гаусса-Лагерра приведено в Таблице 1.

Как видно из таблицы, моды Гаусса-Лагерра с одинаковыми собственными значениями дают достаточно степеней свободы, чтобы сформировать изображение, хорошо аппроксимирующее заданное и не меняющее своей структуры при прохождении Фурье-каскада. Заметим, что функции Гаусса-Лагер-ра с взаимно противоположной "закруткой" ехр[±/'|ю|ф] имеют одинаковые собственные значения.

Совокупность мод Гаусса-Лагерра с одинаковыми собственными значениями можно считать собственной группой, так как она также представляет собой собственную функцию преобразования Фурье. Интересно отметить, что условие равенства скоростей удовлетворяет и условию равенства собственных чисел:

а,по3

что легко видеть из соотношения: Я | =(-1)"Н)'_1" = (-0'(-1Г(-0"5" =СОП51.

ят1п^П€В[ 4 у „ д_у_,

Это понятно и просто из логических соображений, считая, что Фурье-образ изображения соответствует

изображению, прошедшему бесконечное расстояние.

В случае вращения, П2 п £1, = П2 только при 7Ь4/, /-целое. Такой вывод можно сделать, учитывая выражение для количества оборотов, которое успевает совершить световой пучок (15) на расстоянии от 2=0 Д0 7=0С [11]:

Т

К = —. 4

Таблица 1. Собственные значения для собственных функций Гаусса-Лагерра

о Iii!

1 •1 /

2 -? -1 Ш

3 -1 / -i

4 № -i -1 i Щ

5 -1 i я -i -1 i

6 -i -1 i ж -i -1

7 -1 i m -i -1 i И -/ 5И553 iSSiSf

8 № -i -i i -i -1

9 -1 i ш -i -1 i Iii -/ -1 i

5. Численные примеры

При численном моделировании использовались следующие параметры: 128 отсчетов но радиусу г и 128 отсчетов по угловой составляющей сре[0,2я]. Диапазон изменения аргументов r/сте [ОД'Ь], где соотношение R/'гт=15 в большинстве случаев, иногда использовалось Л/сг=10. Параметр ст=0.1 мм, длина волны А.=0.63 мкм, фокусное расстояние /=100 мм. Во всех вычислениях рассматриваюсь функции ГЛ с номерами п,т< Л-15.

На рис. 4 представлены результаты аппроксимации (1) некоторых фигур (рис. 4а) с помощью мод ГЛ.

Для получения аппроксимирующих функций, представленных на рис. 4в, было достаточно двух итераций алгоритма (6)-(7), далее наступай стагнация. Распределение коэффициентов, участвующих в (1) представлено на рис. 46.

Понятно, что в общем случае эти изображения не сохраняются ни при распространении, ни при фокусировке линзой. Для формирования стабильного изображения необходимо соблюдать условие равенства скоростей (11). На рис. 5 представлены результаты аппроксимации (12) фигуры "крест" (рис. 5а Л/а-10) при значении параметра р= 12 (рис. 56,в) и р= 24 (рис. 5г,д). Видно, что чем больше значение р из условия (11), тем больше степеней свободы и лучше аппроксимация. Адаптированный для (12) алгоритм (6)-(7) в данном случае не приводил к стагнации и демонстрироват постоянную сходимость (приведенные результаты соответствуют 11-ой итерации).

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

Из рис. 6 видно, световой пучок (12) сохраняет свою структуру с точностью до масштаба на всем

пути распространения. На рис. 7 представлен аналогичный результат для другого изображения.

Набор изображений, хорошо аппроксимирующиеся при условии равенства скоростей (11), достаточно ограничен. Возможность использования аппроксимации (15), учитывающей условие вращения (13) или (14), позволяет расширить этот набор.

и о в

Рис. 4. (а) заданное изображение, (б) распределение коэффициентов ГЛ, (в) аппроксимация изображения

т

Г

-15 -12 -} о -3 3 с 9 12 15

121

I

I

!5'

I-

т

I I I

-15 -12 -} -О -3

3 г •) 12 15

12

15'

Рис. 5. Исходное изображение

«Л

(О) 2 /

\

*Г

(в) 2= 100 мм (г) 7"-150 мм (л) ^-200 мм (с) г 250 мм Ш >00 мм

Рис. 6. 1'аспреде.гение интенсивности и фазы во входной лшспости

На рис. 8 показаны результаты аппроксимации (15) для некоторых заданных изображений (рис. 8а) после 11 итераций алгоритма (6)-(7) при различных параметрах условия вращения (14): "полукольцо" -Г=-1, (^=18. "наклонная линия" - 7=2. <7=10. "треугольник" -Т=\. <7=22.

Результаты численного моделирования повеления таких световых п\чков при распространении в свободном пространстве и в <|юкусс линзы приведены на рисунках У-II В соответствии с формулой (26). за весь путь от г=0 до г-г "полукольцо" должно совершить четверть оборота против часовой стрелки, "треугольник" - столько же оборотов по ча-

совой стрелке, а "наклонная линия" - половину обо- рота по часовой стрелке.

m

-IS -12 -•> -( -3 > > ' I? 15

>Л>

€

(в) 2=0

jgs1

ЛчЛ

Л

•v

&У'

(г) r-ioa vfnî (,т) 2=150 M.vf (p) z=200 мм (ж) 2=250 mi (j) 2=300 mi

Pue 7. Исходное изображение (a), распределение коэффициентов ГЛ (б), распределение интенсивности и фазы аппроксимации (12) во входной плоскости (в), распределение интенсивности на различных расстояниях от входной

плоскости г

Г"

-15 -12 - -• -:< Ï I- •} 12 15

с

I \ I

! i

{

...__I

' п Ш-15 -12 -t -Î • 3 I• ') 12 ti

/Л

.1 I-

!

"-15 -12 - i -3 0 ^ •) 12 15

*

? ». i

а о в

Рис. Л Исходное изображение (а), распределение интенсивности аппроксимации (15) (б), распределение

коэффициентов ГЛ

(В) 7= 100 ММ (I ) 7~ ] 50 ММ (л) 7-200 мм (С.) 7-250 ММ (Ж) 7 - ПИ' мм

Рис. 9. Распредс, /ение интенсивности и фты во входной пшскости (и). распределение иитепс.инностч и фалы в фокшъпой плоскости мины (о), распределение интенсивности ни ракитных расстояниях от иходной ¡носкости /. (в-

мV

V.

(И) 7-0

&

<ж>

Г\

1

\ '

. Л

V" *

(В)7~50мм (г; 7-1 Си'» ММ и)Г- 150 ММ (г) 7- 20«' ММ (Ж/7 2>0ММ

Рис. ]0. Распределение интенсивности и фат со сходной п. юскости (а), ршпреде. итие иппишшпости и фон.! >. фока льной шиккоапи литы (0), распределение интенсивности пи ршличпы.\ расстояниях от входной а юскости / (в-

V)

Из рис. 9-11 видно, что световые пучки (15) сохраняют свою структуру с точностью до масштаба и поворота на протяжении всего пути. При этом они могут вращаться с различной скоростью и в разных направлениях в зависимости от параметра угловой скорости Г из (14). Скорость вращения пучков в свободном пространстве не постоянная, она замедляется с ростом пройденного расстояния [11]:

Т

V =

1 + (г/г0)2

Интересно заметить, что благодаря симметрии, изображение может повторяться не только при пол-

ных оборотах, но и чаще (см., например, рис. 10а,6 -отрезок повторился на полуобороте).

Ранее было показано, что световые пучки, удовлетворяющие условию равенства скоростей (11) автоматически удовлетворяют и условию равенства собственных чисел (22) или (23), то есть являются собственной групповой функцией преобразования Фурье. Вращающийся пу чок (15) будет собственной функцией преобразования Фурье если угловая скорость Т кратна 4. Однако, из примера, приведенного на рис.10, видно, что вращающийся световой пучок может быть Фурье-инвариантен и с другими значениями Г благодаря своей симметрии.

<а) 2-0

(в) 2=100 мм (г) 2=150 мм (д) 2=200 мм (е) 2=250 мм (ж) 2=300 мм

Рис. 11. Распределение интенсивности и фазы во входной плоскости (а), распределение интенсивности и фазы в фокальной тю скости линзы (б), распределение интенсивности на различных расстояниях от входной плоскости г (в-

ж)

На рис. 12 представлены результаты аппроксимации (24) изображения "кольцо" при соблюдении условия равенства собственных чисел (22) после 20 итераций. Были опробованы две группы индексов {(я, т) (рис. 126-г) и {(пут): Лпт = (рис.

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

Такое изображение должно быть Фурье-инвариантно, то есть повторяться в фокусе линзы, что и подтверждает рис. 13.

Однако, видно, что при распространении в свободном пространстве такое изображение не сохраняется. Интенсивность по мере распространения как бы растекается от центра к периферии. Еще один пример Фурье-инвариантного светового пучка приведен на рис. 14. Для аппроксимации (24) изображения "крест" была выбрана группа индексов {(«, т) : Лпт = 1 .

15 -'2 -I

! 1 3

О

¡4

-15 и

о

Рис. 11 Исходное изображение

(а; -о

! i

^ 'Л

I

(в) z~~-100 АШ

(ж) .Г--300 мм

(г) J50 мм (л) Г---200 мм (с) Г--250 мм

Рис. J3. Распределение интенсиепоспш и фазы но входной /носкости

•U -В • ) J

1НЁ№

f#i

ïktMa-î;

(В) 2-0

!_______

(Д) 2=100 мм

(е) 2~ 150 мм

(ж) z-2u<) мм

{ >.) 2=250 мм

H H

-и! I

"ЦК*!

(и) 2" 300 M\1

Рис. ¡4. Исходное изображение (а), распределение коэффициентов i Л (б), распределение ишпенсинности и фат во входной плоскости (в), распределение интенсивности и фит и фока, шиш /носкости лип tt>i (г), рос преде tenue интенсивности пирих'ш иных расстояниях от uxoûiwû ямскоини i (д-и/

Интересно сравнить различные аппроксимации для одного и того же изображения, например, рис. 4, 5 и 14.

Заключение

В данной работе получены следующие результаты:

- адаптирован и исследован итеративный алгоритм для аппроксимации изображений модами ГЛ;

- показана возможность формирования инвариантных с точностью до масштаба и поворота световых пучков, поперечное сечение которых представляет собой заданное изображение;

- введено понятие групповой собственной функции, как суперпозиции мод ГЛ с одинаковыми

собственными числами, не изменяющей своей структуры после прохождения Фурье-каскада;

Благодарность

Работа выполнена при поддержке Российского фонда фундаментальных исследований

Литература

1. Vasara A.,Turunen J.,Friberg А.Т. Realization of general nondiffracting beams with computer-generated holograms. J.Opt.Soc.Am. A 6

2. Khonina S.N., Kotlyar V.V, Soifer V.A., Shinkaiyev M.V., Uspleniev G.V. Trochoson. Optics Communications, 91

3. Paterson С.,Smith R. Higher-order Bessel waves pro-

duced by axicon-type computer-generated holograms. Optics Communications, 124, 123-130

4. Soifer VA, Golub M. A., Laser beam mode selection by computer-generated holograms, CRC Press, Boce Raton, 1994

5. Хонина C.H., Котляр B.B., Сойфер B.A., Дифракционные оптические элементы, согласованные с модами Гаусса-Лагерра, Компьютерная оптика, Самара-Москва, No. 17, с.25-31

6. Kotlyar V.V., Khonina S.N., Soifer V.A., Generalized Hermite beams in free space, Optik, v. 108, No.l, pp.20-26

7. Chavez-Cerda S., McDonald G.S., New G.H.C. Non-

diffracting beams: travelling, standing, rotating and spiral waves. Optics Commun. 123, 225-233

8. Kotlyar VV, Soifer V.A., and Khonina S.N., An algorithm for the generation of laser beams with longitudinal periodicity: rotating images, Journal of Modern Optics, v.44, No.7, pp. 1409-1416

9. E. Abramochkin, V. Volostnikov, Spiral-type beams. Opt. Comm., v. 102, No.3-4, pp.336-350

10. Котляр B.B., Сойфер В.А., Хонина C.H., Алгоритм расчета ДОЗ для генерации вращающихся модальных изображений, Автометрия, СО РАН, Новосибирск, No.5, с.46-54

11. Котляр В В., Сойфер В.А., Хонина С.Н., Вращение световых многомодовых пучков Гаусса-Лагерра в свободном пространстве, Письма в ЖТФ, т.23, вып. 17, с. 1-6

12. Бахарев М.А., Котляр В.В., Павельев B.C., Сойфер В.А., Хонина С.Н. Эффективное возбуждение пакетов мод идеального градиентного волновода с заданными фазовыми скоростями, Компьютерная оптика, Самара-Москва. No. 17, с. 21-25

13. Е. Abramochkin, V. Volostnikov, Spiral-type beams: optical and quantum aspects. Opt. Comm., v. 125, pp.302-323

14. J.R. Fienup. Phase retrieval algorithms. A comparison. Appl. Opt, 1982, v.21, No. 15, pp.2758-2769.

15. Yariv A. Optical electronics, Kolt, Rinehart, and Winston, New York, 1985.

16. Adams M.J. An introduction to optical waveguides, J. Wiley & Sons, Chichster, 1981.