Связь фазовой проблемы в оптике, фокусировки излучения и задачи Монжа–Канторовича Текст научной статьи по специальности «Физика»

СВЯЗЬ ФАЗОВОЙ ПРОБЛЕМЫ В ОПТИКЕ, ФОКУСИРОВКИ ИЗЛУЧЕНИЯ И ЗАДАЧИ МОНЖА-КАНТОРОВИЧА

Н.Л. Казанский '"2, С.И. Харитонов 2, И.Н. Козлова 2, М.А. Моисеев1 1ИСОИ РАН - филиал ФНИЦ «Кристаллография и фотоника» РАН, 443001" Россия, Самарская область, г. Самара, ул. Молодогвардейская" д. 151;

2 Самарский национальный исследовательский университет имени академика С.П. Королева, 443086" Россия, Самарская область, г. Самара, Московское шоссе, д. 34

Аннотация

Статья посвящена использованию вариационных принципов для решения фазовой проблемы в оптике. В работе рассмотрена связь четырёх фундаментальных задач: фазовой проблемы в оптике, обратной задачи фокусировки когерентного излучения, задачи оптимального перемещения масс Монжа-Канторовича и вариационных методов решения модифицированного уравнения Монжа - Ампера. Показано, что решение фазовой проблемы в оптике в рамках асимптотического подхода тесно связано с решением задачи оптимального перемещения масс с неквадратичной функцией стоимости.

Ключевые слова: фазовая проблема в оптике, задача оптимального перемещения масс, уравнение Монжа-Ампера.

Цитирование: Казанский, Н.Л. Связь фазовой проблемы в оптике, фокусировки излучения и задачи Монжа-Канторовича / Н.Л. Казанский, С.И. Харитонов, И.Н. Козлова, М.А. Моисеев // Компьютерная оптика. - 2018. - Т. 42, № 4. - С. 574-587. - Б01: 10.18287/2412-61792018-42-4-574-587.

Введение

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

Для восстановления фазовой информации Герш-берг и Сэкстон [1 - 4] первыми стали привлекать данные об интенсивности, относящиеся не к одной, а к двум плоскостям. Особенно удачным у них оказался итерационный способ нахождения фазы, который заключается в следующем. Допустим, что исходное монохроматическое поле, несущее информацию об исследуемом объекте, было зарегистрировано в двух плоскостях: во входной плоскости и в Фурье-плоскости.

Обозначим поле во входной плоскости Е0(и, V), а в плоскости изображения - Е(х,у). Пусть при регистрации поля получена информация о функциях Е0(и, V) и Е(х, у), а информация о фазовых множителях отсутствует. Задача состоит в построении комплексной функции по заданному модулю и модулю Фурье -образа.

Фазовая проблема в оптике непосредственно связана с решением обратной задачи фокусировки. Решению этой задачи посвящено множество работ [5 - 10].

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

щем: пусть требуется преобразовать когерентное излучение, имеющее в одной плоскости распределение освещённости Е0 (и, V), в пучок света, который в другой плоскости имеет распределение освещённости Е(х, у). Так как мы знаем распределение модуля комплексной амплитуды света в двух различных плоскостях, то можно найти распределение фазы (или эйконала) Ф0 (и, V) во входной плоскости. Теперь у нас есть полная информация о волновом поле. В результате решения фазовой проблемы нам становится известно: чтобы сформировать распределение поля с освещённостью в выходной плоскости Е(х, у), необходимо сформировать во входной плоскости излучение с комплексной амплитудой -у/Е0(и, V) ехр(/'Ф0(и, V)). На

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

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

В своё время инженер и математик Гаспард Монж сформулировал следующий вопрос: какой самый дешёвый способ транспортировки кучи песка с минимальными затратами, если знать, что перемещение отдельной частицы из положения х в положение у имеет стоимость с(х,у). Оказалось, что решение задачи оптимальной транспортировки песка тесно связано с решением классического уравнения Монжа-Ампера. Эта связь изучается в работах [11 - 12]. В этих статьях установлено, что решение классического уравнения Монжа- Ампера можно свести к задаче оптимального перемещения масс с квадратичной функцией стоимости [13 - 14].

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

1. Фазовая проблема в рамках геометрической оптики

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

В рамках скалярной теории дифракции монохроматический свет описывается комплексной функцией и(х, у, 2). Эту функцию иногда называют комплексной амплитудой света. Комплексная амплитуда удовлетворяет уравнению Гельмгольца

ди + к 2и _ 0.

(1)

где Ди - лапласиан, к = ю/с, ю - циклическая частота света, с - скорость света. Представим комплексную амплитуду света в виде

и (х, у, 2) = А (х, у, 2) ехр (/к Ф (х, у, 2)), (2)

где Ф (х, у, 2) - эйконал, А (х, у, 2) - амплитуда световой волны. В результате уравнение сводится к системе

УА2УФ + А2 дФ = 0,

(УФ)2 = 1 = 1 +

ДА

к 2 А'

(3)

Первое уравнение - это уравнение переноса, которое можно интерпретировать как закон сохранения энергии, а второе уравнение назовём дифракционным уравнением эйконала. Отличие дифракционного уравнения эйконала от обычного состоит в наличии члена ДА/(к'А). В асимптотическом пределе вышеприведённые уравнения имеют вид [15]

(УФ)2 = 1,

УА2УФ + А2 дФ = 0.

(4)

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

Для решения этих уравнений введём лучевые координаты. Лучевые координаты появляются как решение характеристической системы уравнений, которая возникает при решении уравнения эйконала методом характеристик. Связь лучевых координат с декартовыми имеет вид

+ дФ о + дФ о х _ и + —-1, у _ V +--;-1,

ди

дv

2 _ ^, ^ =71 — Ф0и

(5)

где Ф0 (и, V) - распределение эйконала в плоскости 2 = 0, I - расстояние вдоль луча. Решение уравнения эйконала в лучевых координатах имеет вид [15 - 16]

Т_Т(и,V) = Ф0 (и,V) + I . (6)

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

и _ и ( х, у, 2), V _ V ( х, у, 2) , I _ Ь (х, у, 2 ) . (7)

В этом случае решение уравнения эйконала в декартовых координатах имеет вид.

Ф( х, у, 2 ) =

_ Т (и (х, у, 2), V (х, у, 2), Ь (х, у, 2)) = (8)

= Ф0 (и (х, у, 2), V (х, у, 2)) + Ь (х,у, 2).

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

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

дФ = Ап дТ А21 дТ + А31 дТ дх (йе/) ди (йе/) дv (йе/) д1 дФ __ А12 дТ + А 22 дТ А 32 дТ ду (йе/) ди (йе/) дv (йе/) д1 дФ _ А13 дТ А 23 дТ + А 33 дТ д2 (йе/) ди (йе/) дv (йе/) д1

(9) (10) (11)

где

А11 _

А12 _

А 32 _

ду д2 ду д2

дv д1 д1 дv

А ду д2 ду д2

ди дv дv ди

дх д2 дх д2

дv д1 д1 дv'

дх д2 дх д2

ди дv дv ди

дх ду дх ду

А 21 _

ду д2 ду д2 ди д1 д1 ди

А 22 _

дх д2 дх д2 ди д1 д1 ди

а а _дхду

дv д1 д1 дv ди д1 д1 ди

А33 _

дх ду дх ду ди дv дv ди

(12)

(13)

(14)

дх дх дх

ди д ~dl

^ ду ду

ди 1v д1

3z 3z 1z

ди 1v д1

(det) =

= S-1 {~Ф2т12 +Фии Ф vvl2 +

+ (Фии Ф2 +Фуу Ф2 + 2Ф vu Фи Ф v ))}+ (15)

+ {{21ФтФvФu + S} -(ФvvФ2 +Фии Ф2 )l +

+ (Фии )l -Фш12 + ФvvФuul2 } S-1. Градиент функции эйконала, записанный в луче-

вых координатах, имеет вид

р =5Ф = 5Ф = йг p =5Ф

Pz =

дх дФ 1"

ди дФ

~дТ

д1 дZ

ду

дФ

ду '

ду д1'

(16)

Анализ полученных выражений показывает, что декартовые компоненты градиента функции эйконала могут быть выражены через производные функций, входящих в определение лучевых координат. Рассмотрим теперь выражение для оператора Лапласа. Оператор Лапласа от функции эйконала имеет простой вид (см. Приложение А)

5 logD (I)

АФ=др^ +dpy_ + др±

дх ду 3z

где

Px =

Pz =

дФ дх дФ

Hz'

ди

дх 1' 3z_ д'

Py =

д1

дФ

ду

дТ

3v

дУ

д1'

(17)

(18)

Я (l ) =

дх дх дх

ди Hv 1l

ду ду_ ду

ди dv ы

1z 1z 1z

ди 1v 3l

(19)

Используя вышеприведённые выражения для градиента и оператора Лапласа функции эйконала в лучевых координатах, получим, что уравнение переноса имеет вид

if+A2! {(i )}=»■

(20)

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

А2 (и,V,I) = А2 (и,V,0)D-1 (I)D(0) . (21)

Осуществляя несложные преобразования, получим также, что решение уравнения переноса в координатах (и, V, ¿) имеет вид

A2 (и,v, z) =

A2 (и, v, 0)

Xu Yv Xv Yu

где X (и, v, z ) = и +Фои—,

Sz

Y (u,v,z) = v + <$0vf-, =V1 -Ф2и -^5v .

(22)

(23)

Следует отметить, что система координат (u, v, z) не совпадает с лучевой системой координат, однако она более удобна для решения фазовой проблемы, так как уравнения плоскостей, в которых определены распределения освещённостей, имеют простой вид z = const. Используя обозначения

A2 (и,v, z) = E(х,у, z) A2 (и,v,0) = E0 (и,v), (24)

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

E(X, Y)(XuYv -XvYu ) = E0 (и,v) , (25)

х = X (и, v, z), у = Y (и, v, z). (26)

Если считать функции E0 (u, v) и E(X, Y) известными, то получается дифференциальное уравнение в частных производных для определения функции распределения эйконала в плоскости z = 0. Это уравнение будем называть модифицированным уравнением Монжа - Ампера. Если известна функция эйконала в некоторой плоскости, легко можно восстановить комплексную амплитуду во всём пространстве.

Во введении было отмечено, что фазовая проблема в оптике тесно связана с обратной задачей фокусировки когерентного излучения. Пусть требуется преобразовать когерентное излучение, имеющее в одной плоскости распределение освещённости E0 (и, v), в пучок света, который в другой плоскости имеет распределение освещённости E(х,у). Так как мы знаем распределение модуля комплексной амплитуды света в двух различных плоскостях, то можно найти распределение фазы (или эйконала) Ф0 (u, v) во входной плоскости. Для того, чтобы сформировать распределение эйконала, на пути входного пучка помещается рефракционный оптический элемент с криволинейной поверхностью. В случае, если входной пучок имеет плоский волновой фронт и выполняются условия применения параксиального приближения, уравнение этой поверхности принимает простой вид

Фр (и, v)

h (и, v) = -

(n -1)

(27)

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

к (и, V) = (п - 1)шоа2п ^0 (и, V) . (28)

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

2. Задача Монжа - Канторовича

Рассмотрим пару сопряжённых преобразований

Прямое лучевое преобразование Т: (и, V) ^ (х,у) имеет вид

х _ X (и, V), у _ У (и, V) . (29)

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

Обратное лучевое преобразование Я: (§, п) ^ (и, V) (сопряжённое преобразование) имеет вид

и _ и (5,п), V _ V (5,п). (30)

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

Можно показать, что для пары преобразований

и _ и (5,п), V _ V (5,п), (31)

X (и, V), п_ У (и, V) (32)

Я—1 )_(XuУv _XvУu )—1(Я), (33)

выполняется следующее соотношение

( ди _ди_ V

^ дп дп где

и _ и (X (и, V), У (и, V)), V _ V (X (и, V), У (и, V)), (34) ^ _ X (и п), V п)), п _ у (и п), V п)). (35)

Значок Я 1 означает, что равенство (33) выполняется, если после нахождения производных применить преобразование (35).

Рассмотрим минимизацию функционала

£ Е0 (и,V)р(и,V,X,У)дм (IV, (36)

р(и,V,X,У)(и _X)2 + (_У)2 + ^2 (37)

при условии

Е0 (и, V)_ Е(X ( и, V ), У (и, V ))х х ((и (и, V)Уv (и, V) _ XV (и, V)Уи (и, V)).

8(XuУv _XvУu) _ 8XuУv + Xu5Уv _5XvУu _XV5Уи _ ^Ю1и (и, V)У + XuЮ2v (и,V) _ УиЮlv (и, V) _XVЮ2и (и, V)

(38)

2.1. Вариация функционала

Рассмотрим условную минимизацию. Построим функционал

Б _|{Е0 (и,V)р(и,V,X,У) +

+Х(и,V)(Е(X,У){XuУv _XvУu)_Е0)}ди(V,

(39)

X (и, V) - неопределённый множитель Лагранжа.

Получим выражение для вариации этого функционала.

Для этого рассмотрим вариацию функций

Xs (и,V) _ X(и,V) + ею! (и, V), У8 (и,V) _ У (и,V) + SЮ2 (и, V).

(40)

Вариация функционала в этом случае имеет вид

5Б _

_ |Е0(и,V)

(X _ и)ю1 (и, V)

,/(и _ X )2 + ( _ У )2 + ^2 (У _ V)ю2 (и,V)

(и (V +

7(и _X)2 + (v _ У)2 + ^2 +/х(и, V)(5Е(X, У_XvУu ) + +Е (X, У) 5 (Xи^ _ XvУu )) аи (V + +|5Х (и, V)(Е(X, У) (XuУv _XvУu) _ Е0)((и (V.

(41)

Условие

5Б 5Х(и, V)

_ 0

(42)

приводит к уравнению

Е(X, У)(XuУv _ XvУu)_ Е0. (43)

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

Вариация 5E(X,У) имеет вид

5Е ( X, У )_

дЕ (X, У )

дX

Ю1 (и, )-

дЕ (X, У )

дУ

Ю2 (u,V)

(44)

_ s

Х-и XV Уи

Вариация 5 (XuУv - XvУu)

(иу _XvУu ).

Запишем полученное выражение в виде (см. Приложения В и Г)

5(ХиУ _ XvУu ) _

sTr ((т (и, V) ) 1Бю(и, V))

(XuУv _ XvУu ).

(45)

(46)

+

Здесь введены обозначения

^Ю1и (и,V)Уv + XuЮ2v (и,V)_УиЮь, (и,V)_XVЮ2и (и,V)

ч XuУv _XvУu

Подставляем выражение для вариации якобиана, получим выражение для функционала в следующем виде

Л

5Б _sJ Е0 ( u, V )

(X _ и )ю1 (и, V)

7(и _ X )2 + (v _ У )2 + Е2 (У _V)Ю2 (и,V)

(и (V +

(47)

^(и _ X )2 + (V _ У )2 + Е 2

+| X (и, V) (5Е (X, У ) (XuYv _ XvУu) +

+Е (X, У ) (XuУv _ XvУu) sTr ((ДТ)_ Дю)) (и (V.

Подставляя в эту формулу выражение для вариации функции Е(X, У), получаем выражение для вариации функционала в виде

5Б _ sj Е0 (и, V)

(X _ и )ю1 (и, V)

7(и _ X )2 + (V _ У )2 + Е 2 (У _V)Ю2 (и,V)

(и (V +

^(и _ X )2 + (V _ У/ )2 + Е 2 +s^X(u,v)(XuУv _XvУu)((Ех (X,У)Ю1 +

+ЕУ (X,У) ю2) + Е(X,У) Тг ((ДТ)_Д>ю)) (и (V.

(48)

Далее, учитывая, что якобиан лучевого преобразования записывается в виде

Е0(и,V)

(и^ _XvУu ) _

Е (X ,У )

(49)

получаем окончательно выражение для вариации функционала в виде

5Б _sj Е0( и, V)

(X _ и )ю1 (и, V)

у/(и _X)2 + ( _ У)2 + Е2 (У _ V )Ю2 (и, V)

(и (V +

^(и _X)2 + (V _ У)2 + Е2

+síX(u,v)^(ЕК (,У)Ю1 (u,v)-

+ (Еу (X,У)Ю2 (и,V)) + +Е0 (и, V)Тг ((ДТ)_Дю)) (и (V. Рассмотрим интеграл

(50)

(51)

Тг ((ДТ (и, V)) 1Дю(и, V) /1 _ |х(и,V)х

х[ (Е^ (X,У)Ю1 (и,V)-

+ Еу (X,У)Ю2 (и,V)) + +Е0 (и,V)Тг ((ДТ)_1Дю)) (и (V. При замене переменных

и _ и (5,п) V _ V (5,п) (52)

выражение для (ДТ)-1(Я)Дю (Я) преобразуется следующим образом

(ДТ )_1(Я )Дю (Я )_

(|,п) 5^1 (|,п)Л

дп

д^2 (|,п) д^2 (^п)

дп

Тг ((ДТ)_1Дю)

(53)

д^1 (5,п) + ап 2 (^п)

дп

Использованы следующие обозначения

□1 (|,п)_Ю1 (и (|,п), V (|,п)),

о 2 (|,п)_Ю2 (и (|,п), V (|,п)). Сделаем в нем замену переменных /1 _|х(и (I, п), V (I, п)) ((Е§ (I, п) О (I, п) -+ Еп(|,п)02 (|,п)) + Е (I, п) х

(54)

(55)

(I, п) д02 (I, п)

д!

дп

(I (п_

//

^(.п)^ (I,п)Е(п)) +дп(^2 (,п)Е(I,п)))]

_/Л^п)^^ (,п)ЕМ))) ]¿п-

(56)

|Л (I,п)Сд|-(^2 (I, п)Е(I, п))((п) ]

д дп

где введено следующее обозначение Л(|,п)_х(u (|,п), V ^,п)). Далее интегрируем по частям

(57)

+

+

+

11 ==^[—(Ъп)Е (%,ц))%-

-I

5Л(%,Ц)

д%

^2 (%,П)Е(%,ц)%

dn +

(58)

+ [1дП(л(%'П)° 2 (%,П)Е (%,ц))-

-|дЛ(П'П)0 2 (%,П)Е .

Первый интеграл и третий интеграл берутся элементарно и, учитывая факт, что на границе Е (%,ц) ^ 0, окончательно получаем

/1 = -I дА(ПП1 (%,ц)Е(%, ц)d%dn-

2 (%,п)Е(Ъп^а. В результате выражение для вариации имеет вид

(59)

de

= |Е0 (и, V)

(X - и )а>1 (и, V)

у/(и - X)2 + ( - У )2 + Р2

(У - V )сй2 (и, V)

у1(и - X )

+ (V - У)2 + Р2

du dv -

(60)

-I

ГдЛ(%,п) (% ) дЛ(%,п) ^

-^1 (%,п)

2 (%,п)

д% дц

хЕ (%,ц)% dn = 0.

Сделаем замену переменных в первом интеграле ¡\ Е (и, V )кЦШ а, („о*

п-У ((%,ц)) " ь (%,п)

ь (%,п) а 2 (%,ц)\\ <з% dn-

(61)

^М^а. (%,п)+дЛ|ц)а 2 (V

хЕ (%,п)а% dп = 0.

дц

В результате получаем, что условие экстремума функционала приводит к условиям

дл(%%,п) = %-и ((%,п))

5% Ь (% ,ц) ' дл(%,п)_п-У ((%,п))

дц Ь (%, а) '

Дополнительное условие

55"

(62)

5Х(и, V)

_0

(63)

приводит к уравнению

Е(X,У)(XuУv -XvУu) _ Е0 . (64)

Следует отметить, что это уравнение представляет ограничение, которое использовалось при записи исходного функционала.

2.2. Обсуждение полученных формул

Недостаток полученных выражений состоит в том, что в уравнения (62) входят функции и(%, ц) и У(%, ц), которые описывают преобразование точек, расположенных в области фокусировки, в точки, расположенные в области оптического элемента. Уравнение (64) содержит функции X(u, V), У (и, V), описывающие преобразование от точек, расположенных в области оптического элемента, в точки, расположенные в области фокусировки. Используя замену координат

%_ X (и, V), ц_ У (и, V), перепишем уравнение (64) в виде

Е(%,ц)-Е0 (и(%, ц), V(%,ц ))(и%Уц -иц¥%)_ 0 .

Теперь во все полученные уравнения входят функции и (%, ц) и У(%, ц), описывающие связь между точками фокальной плоскости (плоскости г = Р и точками в плоскости, в которой расположен оптический элемент (плоскости г = 0). В этом случае функцию Л (%, ц) можно интерпретировать как функцию эйконала в плоскости фокусировки (плоскости г = Р).

Из решения уравнения эйконала можно показать, что связь между эйконалом в плоскости г = 0 и эйконалом в плоскости г = Р имеет вид

Л(%,ц)_Ф0 (и,V) + Ь(%,ц)

при условии, что точки в плоскости г = 0 и плоскости г = Р связаны соотношениями

и _ и(%,ц), V _ V(%,ц).

Нетрудно показать эквивалентность системы уравнений, записанной в плоскости г = 0,

дЛ{Ъц)_%-и ((%,ц)) д% Ь (%,ц) ,

дЛ(%,п)_п-V ((%,ц))

дц Ь (%, а)

Е(%,ц) - Е0 (и (%,ц),У (%,ц))(и%Уц - ицУ% ) _ 0

и системы уравнений, записанной относительно функций и переменных в плоскости г = Р,

9Ф0 (и,V) X(и,V)-и ди Ь1(и,V)

дФ0 (и, V) У (и, V)- V ду Ь (и, V)

Е(X,У)(XuУv -XvУu)- Е0 (и, V)_ 0,

где Ь1 (и, V) _

+

Минимизируя функционал каким-либо известным методом, мы находим функции X(u, v), Y(u, v), X (и, v). Так как функцию Л п) можно интерпретировать как функцию эйконала в плоскости z = F, то функцию X (и, v) можно интерпретировать как функцию эйконала в плоскости z = F, записанную в координатах (и, v). Связь функции X (и, v) с эйконалом в плоскости z = 0 имеет вид

Х(и,v) = Фо (и,v) + L\ (и,v) .

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

Заключение

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

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

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

Работа выполнена при поддержке Российского научного фонда (проект № 18-19-00326) (вычисление вариации в задаче Монжа-Канторовича о перемещении масс во 2 разделе, Министерства науки и высшего образования РФ в рамках выполнения работ по Государственному заданию ФНИЦ «Кристаллография и фотоника» РАН (расчёт Лапласиана и получение уравнения Монжа- Ампера в Приложениях А, Б - соглашение № 007-ГЗ/Ч3363/26), а также грантов РФФИ (№17-47-630164, 16-29-09528) (получение основных уравнений в разделе 1, описывающих фазовую проблему).

Литература

1. Gerchberg, R. Super-resolution through error energy reduction / R. Gerchberg // Optica Acta. - 1974. - Vol. 21. -P. 709-720. - DOI: 10.1080/713818946.

2. Papoulis, A. A new algorithm in spectral analysis and band-limited extrapolation / A. Papoulis // IEEE Transactions on Circuits and Systems. - 1975. - Vol. 22, Issue 9. - P. 735742. - DOI: 10.1109/TCS.1975.1084118.

3. Sheppard, D.G. Iterative multiframe superresolution algorithms for atmospheric-turbulence-degraded imagery /

D.G. Sheppard, B.R. Hunt, M.W. Marcellin // Journal of the Optical Society of America A. - 1998. - Vol. 15, Issue 4. -P. 978-992. - DOI: 10.1364/JOSAA.15.000978.

4. Soifer, V. Iterative methods for diffractive optical elements computation / V. Soifer, V. Kotlyar, L. Doskolovich. -London: Taylor & Francis Ltd., 1997. - 244 p. - ISBN: 07484-0634-4.

5. Гончарский, А.В. Решение обратной задачи фокусировки лазерного излучения в произвольную кривую / А.В. Гончарский, В.А. Данилов, В.В. Попов, А.М. Прохоров, И.Н. Сисакян, В.А. Сойфер, В.В. Степанов // Доклады АН СССР. - 1983. - Т. 273, № 3. - С. 605-608.

6. Doskolovich, L.L. A method of designing diffractive optical elements focusing into plane areas / L.L. Doskolovich, N.L. Kazanskiy, S.I. Kharitonov, V.A. Soifer // Journal of Modern Optics. - 1996. - Vol. 43, Issue 7. - P. 1423-1433. - DOI: 10.1080/09500349608232815.

7. Soifer, V.A. Methods for computer design of diffractive optical elements / V.A. Soifer, V.V. Kotlyar, N.L. Kazanskiy, L.L. Doskolovich, S.I. Kharitonov, S.N. Khonina, V.S. Pavelyev, R.V. Skidanov, A.V. Volkov, D.L. Golo-vashkin, V.S. Solovyev, G.V. Usplenyev, ed. by V.A. Soifer. - New York: John Wiley & Sons, Inc., 2002. - 784 p. -ISBN: 978-0-471-09533-0.

8. Doskolovich, L.L. Analytic design of optical elements generating a line focus / L.L. Doskolovich, A.Yu. Dmitriev, S.I. Kharitonov // Optical Engineering. - 2013. - Vol. 52, Issue 9. - 091717. - DOI: 10.1117/1.OE.52.9.091707.

9. Doskolovich, L.L. Analytical source-target mapping method for the design of freeform mirrors generating prescribed 2D intensity distributions / L.L. Doskolovich,

E.A. Bezus, M.A. Moiseev, D.A. Bykov, N.L. Kazanskiy // Optics Express. - 2016. - Vol. 24, Issue 10. - P. 1096210971. - DOI: 10.1364/OE.24.010962.

10. Досколович, Л. Л. Расчёт эйконала светового поля для фокусировки в набор точек / Л.Л. Досколович, М.А. Моисеев, Е.В. Бызов, С.В. Кравченко // Компьютерная оптика. - 2014. - Т. 38, № 3. - С. 443-448.

11. Канторович, Л.В. О перемещении масс / Л.В. Канторович // Доклады Академии Наук СССР. - 1942. -Т. XXXVII, № 7-8. - С. 227-229.

12. Богачёв, В.И., Задача Монжа-Канторовича: достижения, связи и перспективы / В.И. Богачёв, А.В. Колесников // Успехи математических наук. - 2012. - Т. 67, Вып. 5(407). - С. 3-110. - DOI: 10.4213/rm9490.

13. Glimm, T. Optical design of single reflector systems and the Monge-Kantorovich mass transfer problem / T. Glimm, V. Oliker // Journal of Mathematical Sciences. - 2003. -Vol. 117, Issue 3. - P. 4096-4108. - DOI: 10.1023/A:1024856201493.

14. Sulman, M.M. An efficient approach for the numerical solution of the Monge-Ampere equation // M.M. Sulman, J.F. Williams, R.D. Russell // Applied Numerical Mathematics. - 2011. - Vol. 61, Issue 3. - P. 298-307. - DOI: 10.1016/j.apnum.2010.10.006.

15. Кравцов, А.Ю. Геометрическая оптика неоднородных сред / А.Ю. Кравцов, Ю.И. Орлов. - М.: Наука, 1980. -304 с.

16. Кошляков, Н.С. Уравнения в частных производных математической физики / Н.С. Кошляков, Э.Б. Глинер, М.М. Смирнов. - М.: Высшая школа, 1970. - 710 с.

Приложения

А. Расчёт лапласиана

Связь лучевых координат с декартовыми имеет вид

(А.1)

+ дФ 0 + дФ 0 X = и +--;—! у = V +--;—! г = .

ди дv

Эти уравнения можно разрешить и = и(х,у,г), V = V(х,у,г), ! = Ь(х,у,г) . (А.2) Решение уравнения эйконала имеет вид ф( X, у, 2 ) =

= Т(и (х, у, 2 )V (х, у, х ), Ь (х, у, х ))= (А.3)

= Фо (и (х,у, г) V (х, у, г)) + Ь (х, у, г).

Выражения для производных в лучевых координатах через производные в декартовых координатах имеют вид

( дх ду дг ^ ( дФ^ ( дТ^

ди ди ди дх ди

дх ду дг дФ дТ

ду ду ду ду ду

дх ду дг дФ дТ

[д! д! д! [дг у [д! )

(А.4)

( дФ^ ( дх ду дг -1 ( дТ^|

дх ди ди ди ди

дФ дх ду дг дТ

ду ду ду ду ду

дФ дх ду дг дТ

[дг [д! д! д! у [д!

(А.5)

Обратная матрица имеет вид

' дх ду дг ^ 1

ди ди ди

дх ду дг

дv дv дv

дх ду дг

~д1 ~д! ~д!

А11 -А21 А ^

- А12 А22 -А32 , (А.6)

А13 -А23 А33 у

¿ег = ¿ег

' дх ду дг ^

ди ди ди

дх ду дг

дv дv дv

дх ду дг

~д! ~дИ ~дИ

(А.7)

Выражения для величин Ау (I = 1,3;] = 1,3) имеют

вид

А =ду дг ду дг ^ = дх дг дх дг ду д1 д1 дv дv д1 д1 дv

А13 =

дх ду дх ду ду д1 д1 ду

А =ду дг ду дг ^ = дх дг дх дг ди д1 д1 ди ди д1 д1 ди

А дх ду дх ду

ду д1 д1 ду

А =ду дг ду дг ^ = дх дг дх дг ди ду ду ди ди ду ду ди

А 33 =

дх ду дх ду ди ду ду ди

(А. 8)

(А. 9)

(А. 10)

Выражения для производных в декартовых координатах через производные в лучевых координатах имеют вид

дФ А11 дФ А 21 дФ А 31 дФ 11Ч -= 7—ч--7—ч-+ 7—ч-, (А. 11)

дх Ыеп ди Ыеп ду Ыеп д1

дФ А12 дФ А22 дФ А32 дФ -= -7—ч-+ 7—7--7—7—, (А. 12)

ду ) ди ) ду ) д1

дФ А13 дФ А23 дФ А33 дФ /4 10Ч

-= 7—ч--7—ч-+ 7—ч-, (А. 13)

дх (е^ ди (Ле() ду (Ле() д1

(det) = ^ {2/ФиуФуФи + -

- (ФууФ2 + Фии Ф2 )1 + (Фии + Фуу )) (А. 14) -Ф2у12 +Фуу Фии12 }.

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

дФ дФ дФ

— = Ф0и, — = Ф0у, — = ^ 1 "Ф0„-Ф2у . (А.15)

дх ду дг

Используя эти выражения, найдём выражения для дивергенции векторного поля. Дивергенция векторного поля в декартовых координатах имеет вид

¿IV (р ) = дЕ± +Ъу.+дЕ±,

дх ду дг

дрх _ А11 дрх А 21 дрх + А 31 др,

дх (det) ди (е{) ду (det) д1

дру _ А12 дру А 22 дру А 32 дру

ду

(е^ ди (det) ду (е^ д1

др2 = А13 др2 А 23 др2 + А 33 др2

дг (det) ди (еду (ед1

(А. 16)

(А. 17)

(А. 18)

- +-- +

дх ду дг (Се/)

дх ду

Рх _¥' ^ _¥'Рг _ Пусть

дРх +дру + ?Рг _ 1

Г дРх дРх дРх дх дх дх дх дх дх Л

ди дv д/ ди "ду ~/ ди д ~д/

1 ду ду ду + дРу дРу дРу + ду ду ду

(Се/) ди дv д/ ди дv д/ ди дv д/

дг дг дг дг дг дг дРг дРг дРг

У ди дv д/ ди дv д/ ди дг д/ )

(А. 19)

дг_ д1

дх ду дг (Се/)

д дх д дх д дх дх дх дх дх дх дх

ди д/ "д^ ~д/ ~д/~дЛ ди ¥ ди ¥ ¥

ду ду ду + д_ду_ д_ду_ д_ду_ + ду ду ду

ди ду д/ ди д/ дv д/ д/ д/ ди дv д/

дг дг дг дг дг дг д дг д дг д дг

ди д ~д/ ди ¥ ¥ ди д/ ¥ ¥ ¥¥

д дх д дх д дх дх дх дх дх дх дх

д/ ди ¥¥ ¥¥ ди ду ¥ ди ¥ ¥

1 ду ду ду + дду д_ду_ д_ду_ + ду ду ду

(Се/) ди дv д/ д/ ди д/ дv д/ д/ ди дv д/

дг дг дг дг дг дг д дг д дг д дг

У ди ¥ ¥ ди ¥ ¥ д/ ди ¥¥ ¥¥

дх дх дх дх дх дх

ди ¥ ¥ ди ¥ ¥

1 д ду ду ду д ду ду ду

(Се/) д/ ди дv д/ шв ди дv д/

дг дг дг дг дг дг

ди ¥ ¥ ди ¥ ¥

(А.20)

ИТОГО

дРх + дРу + дрг дх дУ дг д1

дх дх дх

ди ¥ ¥

ду ду ду

ди дv д/

дг дг дг

ди дv д/

(А.21)

В случае

дФ

Рх ' ру

дФ ду '

дг Рг = ~д/

(А.22)

получаем выражение для лапласиана функции эйконала

дФ _

д_

~д1

дРх + дРу+ дРг

дх ду дг

дх дх дх

ди дv ¥

ду ду ду

шв ди дv 5/

дг дг дг

ди дv ¥

(А.23)

Б. Получение уравнения Монжа-Ампера А2 (и,V,/)_ А2 (и,V,0)х

дх дх дх

ди дv д1

ду ду дУ

ди дv д1

дг дг дг

ди дv д1

*(0) ^х(0) *(0)

ди д^ д/

ду (0) ду (0) ду (0)

ди д^ д/

±(0) ±(0) £(0)

ди д^ дЛ'

(Б.1)

Это уравнение можно переписать в виде

А2 (и,V,/)_ А2 (и,V,0)([ х ^]г/) х х( (0)х г (0)] г/ (0)),

+ дФ 0 + дФ 0 . х _ и +--/, у _ V +--/,

ди

дv

г — , — Л /1

дФ 0

дФ0

ди ) У ди ([ги х ^]г/) _ (Уuгv - гuУv)х/ + + (гuXv -Xuгv)у/ +(xuУv -УuXv)г/.

(Б.2)

(Б.3)

(Б.4)

х

Рассмотрим получение выражения ХиУу

V + Фои X = и + -

л/1 -Ф2-Ф2

ф 0у

л/1 -Ф2и -ф\ В результате

У = V +•

-2.

где

Хи = 1 + фОиил^ФО"-®^ -Ф0и2и1 -

Х(1 -Ф 2и -Ф )-1 ^ ,

Уv = 1 + ф0v^л/1 -Фии -ФUV -Ф0у2у1-1

-ф 2и -ф )-1 2,

XV = Фои^л/1—ФОи^Ф^-Фои2у1-1

-Ф ии -Ф Uv )-1 2 , Уи = Ф Ои^л/^-ФОи-Ф^ + Ф 0v2ul-1

2и -Ф Uv )-1 2 ,

-2и8г1-1 =(Ф0иФ0ии +Ф0vФ0vu ),

-2vSzl-1 =(Ф0иФ0uv +Ф0vФ0уу ). После элементарных преобразований

Хи = Хи — х12и21 ^ , Уv = \Уv - Уl2v2—1 ] , XV = [XV — хч2у2— 1 ^ , Уи = [Уи - У^и2—1 ] , ХиУу — ХуУи =

= {фиУу - ХуУи )2l +фуУи - Уу2и )х1 + + (ф2и — хи2у -УЧ } Sz 1.

ИТОГО

дх дх дх

ди ~ду ~д!

дУ_ дУ_ ду

ди ду д1

д2 д2 д2

ди ~ду ~дЧ

Sz = ХиУУ ХУУи ■.

Sz =

— (0) — (0) — (0)

ди ду дГ '

^У (0) ^ (0) ^У (0)

- (0) - (0)

ди дуУ '

дl д2_ дl

(0 )

(Б.5)

(Б.6) (Б.7) (Б.8) (Б.9)

(Б.10)

(Б .11) (Б.1и) (Б.13) (Б.14)

(Б.15)

(Б.16)

и (и V, 1 ) = Аи (и,

дх дх дх -1

ди ду д1

дУ дУ дУ

ди ду д1

д2 д2 д2

ди ду д1

— (0) — (0) — (0)

ди ду дП '

^У(0) ^У(0) &(0)

ди ду дГ '

ь (0) * (0) *(0)

ди ду дГ '

(Б.17)

Аи (и,V,I) =

Аи (и, V, 0)

ХиУу Х уУи

В. Выражение для производных Рассмотрим систему линейных уравнений для определения производных

д®1 {и ф,ч),г (!,п))ди (!,п)

дП д%

+ ck»^фффn-Vфn)_ дУ (!,п)

д¥

+ % , д®1 (и (!,п)) ^(,4) +

дП дц

+ дУ (!,п)

дУ дц

дц

Определитель системы

(дП дУ^

А =

Д1 =

А и =

Итого

д! д!

ди дУ

дп

(д®1 дУ )

1! д!

д®1 дУ

ч ^ дп У

(ди да 1_ ^

~д!

ди да 1

,дп дп У

дю1дУ дю1 дУ д% дп дп д"% '

д®1 ди д®1 ди

дп д% д% дп

дю1 дУ дю1 дУ да1 (и(^п)У(5>п)) = д! дц дц д!

ди ~ и дУ — дЦ_ дУ '

д! дп дп д!

, . дщди дю1 ди

да1 (и(!>п)>У(!,п))= п д! д! дп

дУ

ди дУ — ди_ дУ д! дп дп д!

(В.1)

(В.и)

(В.3)

(В.4)

ХУи

X

По аналогии

дю2 дУ - дю2 дУ дю2 (и(^п),У(^П)) _~д^~дп

ди ~ Ж дУ-дЦ_ дУ '

дп дп

, . дю2 ди дю2 ди

дю2 (и(|,п),У(^п))_ дп дп

дУ

ди дУ ди дУ дп дп

(В.5)

Г. Вспомогательные преобразования Г.1. Формулы для прямого преобразования

ч Гхп X. бт (и, V)_

У Уи

У,

Бю (и, V) _

дю 1 (и, V) дю 1 (и, V)

ди дv

дю 2 (и, V) дю 2 (и, V)

ди

дv

(Г.1)

(БТ )-1(и, V )_

Гх„ х,, Л-1

У У

У и V )

_ (Хиу - ВД )

Г У, -X,, Л

- У X

У и и )

Тг ((БТ (и, V )) 1Бю(и, V ))

ди ди " д^

Г Ю1и (и, V )У, + Хи Ю2v (и, V )-Уи Юlv ( и, V )-XV Ю2и (и, V )Л ХиУ - XvУu

Г.2. Преобразование некоторых выражений при замене координат Далее рассмотрим обратное преобразование Я : п) ^ (и, V)

и _ и(5,п) V _ У (|,п), (Г.4)

такое, что

X (и (£,п),У (|,п)) _ У (и (£,п), У (1,п)) _ п,

БТ (Я)_

-Хи X„ У У

У и V

((и (^п),У (^п)))

Г дХ дУ-дХ дУ дХ ди-дХ ди Л

дп дп дп дп

дУ дУ дУ дУ дУ ди дУ ди

дп дп дп дп

ди дУ ди дУ

дп дп

ди дУ ди дУ

дп дп

Г дУ ди Л

дп дп

дУ ди "д^ "д^

Г дХ дХ Л

дп

5У 5У дп

Г ди диЛ 1

дп

дУ дУ дп

(Г.2)

(Г.3)

(Г.5)

(Г.6)

(БТ )-1(Я )_

Г ди ди Л

дп

дУ дУ дп

-1Л

Г ди диЛ

дп

дУ дУ дп

_ БЯ.

(Г.7)

Бю( Я )_

Гди ди УУ дП1 (|,п) дП1 (|,п)

дУ дУ до2 (1,п) до2 (1,п) дп,) У дп ,

Й (|,п)_ю1 (и (|,п),У (|,п)). Й 2 (|,п)_Ю2 (и (|,п),У (|,п)) .

(Г.8)

(Г.9) (Г.10)

(DT) 1(R)Dra(R) = DRDro(R) =

( dU dUYöU

dV_

( dU dV

dn

dV dn

dU dn dV dn

Л

dV dU

"зё"

dV

(Г.11)

Выражение для якобиана обратного преобразования

dU dV dU dV dn dn

= (( -XVYU )1(R).

Сведения об авторах Сведения об авторе Казанский Николай Львович см. стр. 571 этого выпуска.

Харитонов Сергей Иванович, 1961 года рождения. Доктор физико-математических наук, доцент кафедры наноинженерии, ведущий научный сотрудник лаборатории дифракционной оптики ИСОИ РАН - филиал ФНИЦ «Кристаллография и фотоника» РАН. 1984 г. - окончил физический факультет Самарского государственного университета. 1993 г. - защитил кандидатскую диссертацию на тему «Асимптотические методы дифракционного расчёта фокусаторов лазерного излучения». 2010 г. - защитил докторскую диссертацию на тему «Асимптотические методы расчёта дифракции когерентного электромагнитного излучения на дифракционных оптических элементах». Область научных интересов: дифракционная, квантовая оптика, физика плазмы. В списке научных работ С.И. Харитонова 87 статей, 5 авторских свидетельств и патентов. E-mail: prognoz2007@smail.com.

Козлова Ирина Николаевна, 1984 года рождения. В 2007 году с отличием окончила Самарский государственный аэрокосмический университет имени академика С.П. Королева (СГАУ, ныне - Самарский университет) по специальности «Конструирование и технология радиоэлектронных средств». Кандидат технических наук (2013 г.), доцент кафедры наноинженерии Самарского университета. Область научных интересов: обеспечение надежности микро- и нанообъектов, проектирование устройств наноэлектроники и нанофотоники. В списке научных работ И.Н. Козловой 15 статей, 2 патента. E-mail: inkozlova@ssau.ru .

Моисеев Михаил Александрович, 1986 года рождения, в 2008 году с отличием окончил Самарский государственный аэрокосмический университет имени академика С.П. Королёва (ныне Самарский университет) по специальности 010600 «Прикладные математика и физика». Кандидат физико-математических наук (2011 г.), работает научным сотрудником в Лаборатории дифракционной оптики ИСОИ РАН - филиале ФНИЦ «Кристаллография и фотоника» РАН. Является членом международного научного общества SPIE. М.А. Моисеев -специалист в области решения обратных задач геометрической оптики. В списке научных работ М.А. Моисеева 35 статей и 7 авторских свидетельств. E-mail: mikhail@smr.ru .

ГРНТИ: 29.31.29.

Поступила в редакцию 27 июня 2018 г. Окончательный вариант - 28 июня 2018 г.

THE CONNECTION BETWEEN THE PHASE PROBLEM IN OPTICS, FOCUSING OF RADIATION, AND

THE MONGE-KANTOROVICH PROBLEM

N.L. Kazanskiyl2, S.I. Kharitonov 2, I.N. Kozlova 2, M.A. Moiseev1

1IPSIRAS - Branch of the FSRC "Crystallography and Photonics" RAS, Molodogvardeyskaya str. 151, 443001, Samara, Russia;

2 Samara National Research University, Moskovskoye shosse 34, Samara, 443086, Samara, Russia

Abstract

We discuss the use of variational principles for solving the phase problem in optics. In this paper, we consider the connection between four fundamental problems: the phase problem in optics, the inverse problem of focusing coherent radiation, the Monge-Kantorovich optimal mass trans-

port problem, and the variational methods for solving the equation of a modified Monge-Ampere equation. It is shown that the solution of the phase problem in optics within the framework of the asymptotic approach is closely related to the solution of the problem of optimal mass transport with a nonquadratic cost function.

Keywords: optimal mass transport, phase problem in optics, Monge-Ampere equation.

Citation: Kazanskiy NL, Kharitonov SI, Kozlova IN, Moiseev MA. The connection between the phase problem in optics, focusing of radiation, and the Monge-Kantorovich problem. Computer Optics 2018; 42(4): 574-587. DOI: 10.18287/2412-6179-2018-42-4-574-587.

Acknowledgements: This work was supported by the Russian Science Foundation (Project No. 18-19-00326) (calculation of the variation in the Monge-Kantorovich mass transportation problem, Section 2), by the RF Ministry of Science and Higher Education within the State assignment to the FSRC «Crystallography and Photonics» RAS (calculation of the Laplacian and derivation of the Monge-Ampere equation, Appendices A and B - Agreement No 007-0/^3363/26), and by Russian Foundation for Basic Research (Projects Nos. 17-47-630164, 16-29-09528) (derivation of the basic equations describing the phase problem, Section 1).

References

[1] Gerchberg R. Super-resolution through error energy reduction. Optica Acta 1974; 21: 709-720. DOI: 10.1080/713818946.

[2] Papoulis A. A new algorithm in spectral analysis and band-limited extrapolation. IEEE Transactions on Circuits and Systems 1975; 22(9): 735-742. DOI: 10.1109/TCS.1975.1084118.

[3] Sheppard David G., Hunt Bobby R., Marcellin Michael W. Iterative multiframe superresolution algorithms for atmospheric-turbulence-degraded imagery. J Opt Soc Am A 1998; 15(4): 978-992. DOI: 10.1364/JOSAA.15.000978.

[4] Soifer V, Kotlyar V, Doskolovich L. Iterative methods for diffractive optical elements computation. London: Taylor & Francis Ltd; 1997. ISBN: 0-7484-0634-4.

[5] Goncharsky AV, Danilov VA, Popov VV, Prokhorov AM, Sisakyan IN, Soifer VA, Stepanov VV. Solution of the inverse problem of laser radiation focusing in an certain curve [In Russian]. Doklady of the USSR Academy of Sciences 1983; 273(3): 605-608.

[6] Doskolovich LL, Kazanskiy NL, Kharitonov SI, Soifer VA. A method of designing diffractive optical elements focusing into plane areas. J Mod Opt 1996; 43(7): 14231433. DOI: 10.1080/09500349608232815.

[7] Soifer VA, ed. Methods for computer design of diffractive optical elements. New York: John Wiley & Sons, Inc; 2002. ISBN: 978-0-471-09533-0.

[8] Doskolovich LL, Dmitriev AYu, Kharitonov SI. Analytic design of optical elements generating a line focus. Opt Eng 2013; 52(9): 091717. DOI: 10.1117/1.OE.52.9.091707.

[9] Doskolovich LL, Bezus EA, Moiseev MA, Bykov DA, Kazanskiy NL. Analytical source-target mapping method for the design of freeform mirrors generating prescribed 2D intensity distributions. Opt Express 2016; 24(10): 10962-10971. DOI: 10.1364/OE.24.010962.

[10] Doskolovich LL, Moiseev MA, Byzov EV, Kravchenko SV. Computation of light field eikonal to focus into a set of points [In Russian]. Computer Optics 2014; 38(3): 443-448.

[11] Kantorovich LV. On the translocation of masses. J Math Sci 2006, 133(4): 1381-1382. DOI: 10.1007/s10958-006-0049-2.

[12] Bogachev VI, Kolesnikov AV. The Monge-Kantorovich problem: achievements, connections, and perspectives. Russian Mathematical Surveys 2012; 67(5): 785-890. DOI: 10.1070/RM2012v067n05ABEH004808.

[13] Glimm T, Oliker V. Optical design of single reflector systems and the Monge-Kantorovich mass transfer problem. J Math Sci 2003; 117(3): 4096-4108. DOI: 10.1023/A:1024856201493.

[14] Sulman MM, Williams JF, Russell RD. An efficient approach for the numerical solution of the Monge-Ampere equation. Applied Numerical Mathematics 2011; 61(3): 298-307. DOI: 10.1016/j.apnum.2010.10.006.

[15] Kravtsov AYu, Orlov YuI. Geometrical optics of inhomo-geneous media. Berlin, Heidelberg: Springer-Verlag; 1990. ISBN: 978-3-642-84033-3.

[16] Koshlyakov NS, Smirnov MM, Gliner EB. Equations of mathematical physics in partial derivatives [In Russian] Moscow: "Vysshaya Shkola" Publisher; 1970.

Author's information

The information about author Nikolay Lvovich Kazanskiy you can find on page 573 of this issue.

Sergey Ivanovich Kharitonov (b. 1961) leading researcher of Diffractive Optics laboratory in the IPSI RAS -Branch of the FSRC "Crystallography and Photonics" RAS. Doctor of Physical and Mathematical Sciences. 1984 -graduated from Physics department of Samara State University (presently, Samara National Research University). 1993 - defended his dissertation "Asymptotic methods of calculation of the diffraction of laser radiation focuser". 2010 - defended his doctoral thesis on "Asymptotic methods for calculating the diffraction of coherent electromagnetic radiation in diffractive optical elements". Research interests: diffraction, quantum optics, plasma physics. The list of S.I. Kharitonov scientific papers includes87 articles, 5 patents. E-mail:prognoz2007@smail.com .

Irina Nikolaevna Kozlova (b. 1984) graduated with honors (2007) from Samara State Aerospace University named after S.P. Korolyov (presently Samara National Research University), majoring in Design and Technology of Radio-electronic Devices. She received her Candidate in Technical Sciences (2013) degrees from Povolzhsky State University of Telecommunications and Informatics. She is holding a position assistant professor at nanoengineering department.

Current research interests include reliability of micro- and nanoobjects, design of nanoelectronics and nanophotonics elements. She is co-author of 15 scientific papers, 2 patents. E-mail: inkozlova@ssau.ru .

Mikhail Alexandrovich Moiseev (b. 1986) graduated with honors (2008) from S.P. Korolyov Samara State Aerospace University (SSAU), majoring in Applied Mathematics and Physics. He received his Candidate's Degree in Physics and Mathematics (2011) from Samara State Aerospace University. He is a researcher at the Diffraction Optics laboratory of the IPSI RAS - Branch of the FSRC "Crystallography and Photonics" RAS. He is a SPIE-member. He is coauthor of 23 scientific papers and seven inventions and patents. His current research interests include nonimaging optics. http://www.ipsi.smr.ru/staff/MoiseevM.htm.

Received June 27, 2018. The final version - June 28, 2018.