Преобразование Абеля в задачах синтеза градиентных оптических элементов

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

В.В. Котляр, А.С. Мелёхин Институт систем обработки изображений РАН, Самарский государственный аэрокосмический университет

Аннотация

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

Введение

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

Флетчер [2] и Морган [3] получили разные виды интегрального уравнения для обобщенной задачи Лунеберга. Если параллельный пучок лучей падает на обобщенную линзу Лунеберга, то все лучи сходятся в точке на оси внутри или снаружи сферической поверхности линзы. Решение для обобщенной линзы Лунеберга с помощью преобразования Абеля, а также для других типов обобщения задачи Лунеберга найдены в [4-7]. В [6] решена задача расчета сферического градиентного аксикона. В данной работе показана связь между уравнениями Моргана [3] и Флетчера [2], а также, с помощью производной от преобразования Абеля, из решения для обобщенной линзы Лунеберга получен показатель преломления для обычной линзы Лунеберга.

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

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

Для среды с показателем преломления п = п(г), зависящим только от радиальной координаты г в сферической системе координат, имеет место лучевой инвариант, который является обобщением обычного закона преломления Снелиуса [9]:

<• ч 2 60 ,

п(г)г — = п = еоп$1 61

Уравнение (1) выполняется вдоль траектории луча. На рис. 1 показан ход луча (и обозначения) в обобщенной линзе Лунеберга. Все лучи, параллельные оси г, должны сходиться в точке на оси и на расстоянии / от центра линзы.

Радиус линзы г = 1, п(1) = 1, П - расстояние от оси до выбранного луча, г* - минимальное расстояние луча от центра линзы.

-*—7

п г*

еХ а 1 гК г

Г

п

Рис. 1. Ход луча в обобщенной линзе Лунеберга Из рис. 2 видно, что приращение 61 вдоль траектории луча связано с приращением по радиусу 6г и углу 6в следующим образом:

61 =

+ (г60)2 . Из уравнения (2) следует:

- = , г2 +Г -

60 V 160

2

а/

Ссг -- 1 Сг 2 2 2 г* гуп г - И сСп 2 2 2 пуп г - И

гсе

/ 1 пСг - п * гап =

/ п(г) 1 /22 ,2 г* п^п г - И 2 2 2 1 гт! п г - И

се/ 1 пСг + гСп 1 С

/ 2 2 2 ь пгу п г - И ь (пг)^(пг)2 - И2

Рис. 2. Участок траектории луча в градиентной среде Подставим выражение (3) в уравнение (1), получим:

А = Ир

&в И

Из уравнения (4) следует интегральное уравнение для поиска показателя преломления в произвольной сферической градиентной среде:

' аг *

1

=-И! а* .

г1 г^/п2(г)г2 - И2 И *

Интегральное уравнение для обобщенной линзы Лунеберга Интегральные уравнения для расчета показателя преломления обобщенной линзы Лунеберга были получены в [2, 3]. Требуется найти п(г), чтобы параллельные оптической оси лучи сходились на оси в точке на расстоянии / от центра сферы с радиусом г = 1 (/> 1). Из рис. 1 видно, что г1 = г*, п(г*)г* = И .

Т а* . ( й

Так как г— = бш^ (см. рис. 2), то лучевой инвари-

а

ант (1) можно записать в виде п(г)г бш^ = И .

Граничные условия равны: при г = 1, п(1) = 1, тогда = бш*= И; г2 = 1; *1 -*2 = д*=п а *

2

$хпу = —; а = в-у ; д*=|£-*|^ = агссоБ И +

п У.<

/ ' ^2 'У 2 2

Все последние равенства нетрудно получить из геометрии на рис. 1.

Тогда уравнение (5) в частном случае для обобщенной линзы Лунеберга можно записать в виде:

Сг

2агссоБ И + агсБ.

,п( у}

г* г^п 2(г)г 2 - И 2

2И

где И = г п(г ).

Это уравнение было получено Морганом [3]. Аналогично, можно получить уравнение Флетчера [2]:

п

I-

1 п

дп

^п2 (г)г2 - И2

2F

где п = п(г ).

Покажем как перейти от уравнения (7) к уравнению (6):

= 1

¿4

^^И2

агссоБИ

И

Последнее равенство в уравнении (8) следует из табличного интеграла [10].

Прямое и обратное преобразования Абеля Чтобы решить уравнение (7) с помощью преобразования Абеля приведем некоторые предварительные сведения. В [11] дана пара (прямое и обратное) преобразований Абеля с верхним бесконечным пределом:

Р(г) = 2|

Г(х)хс1х

I 2 2

г ух - г

ОТ

г(х) = - -1

ТТ

ар(г) Сг

Сг

Заметим, что повторное применение преобразования (9) к уравнению (10) приводит к выражению для производной от преобразования Абеля:

СГ(х)_

1-р/ ч ОТ - '

СР(г^о „г Сх

-Сх

С"

= 2г I

г-Ц-г 2

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

:(х) = -1Ах |

ТГ Г1V *

Г(г)Сг

Если приравнять правые части выражений (10) и (12), то можно получить полезные выражения для производной от преобразования Абеля:

1 ОТ А1-

Г(г)Сг

Сх'х«/г 2

л/Т2-

=11

V

СГ г с

От| 1---Сг

аг

ОТ

А1-

:(х)хСх

= г 1-

г

аг(х) Сх

Сх

Ь'гг^х2 - г 2 гТх2^

22 х -г

Для прямого и обратного преобразований Абеля в конечных пределах в [13] даны следующие выражения:

Г(х)Сх

Г(г) = 1-

гл/г-х

2

2

од=-

агсг) р(0)+| ^Г

>/х 0 V х - г

6г

Заметим, что в уравнении (15) преобразование Абеля представлено в виде свертки и поэтому для его обращения можно использовать преобразование

Фурье, с учетом того, что функция у г-г является

инвариантной по отношению к преобразованию Фурье [14]:

1 +<» -хй

1 г е *

I

6х =

1

^ЛдД дй'

В интегральных уравнениях (6) и (7) пределы интегрирования отличаются от пределов интегрирования в уравнениях (15) и (16), поэтому приведем пару преобразований Абеля в конечных пределах, удобных для решения уравнения (7):

Р(г) = 2|

1 £(х)х6х

/х

, 2 2 г \Х - г

ад=—

ж

6Р(г) 6

1 -6г

Г 6г_

, 2 2 х-^г - X

Р

л1г 2 - х2 л/1—

Проверим, что уравнение (19) является обращением преобразования (18). Действительно:

Р(г) = |

1 ^ х)6х2

Л.- г 2

1 1 II

— 6г '6х 2 6г

хг^

- Р(!)|

6х2

л/(! - х2)(X2 - г2) ^

I— 6г' I

л Лт' У

6х1

6г

- р(1) I

д/

ах2

>/

= - — -Ы—6г'-яР(1)1 = Р(г) .

п [ х6г' 1

В уравнении (20) предпоследнее равенство по лучено на основе табличного интеграла [10]:

I

6х

-у/

= -агс81п-

г - г

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

I

6х

д/

= п

На основе преобразований Абеля (18), (19) можно получить две разные формы для производной от преобразования Абеля. Первая форма для производной приведена ранее в работах [4-6]. Для ее получения возьмем по частям интеграл в преобразовании Абеля вида:

р(г) =\-г^ = Г =

7 2 2 х - г

1

= Д[х) х]6

Д(1)агссо8 (г) г

- 11 агссов I —

х\]х2 -

^х) + х — 6х

6х

При получении уравнения (23) использовалось соотношение:

_6_

6х

1 I г — агссоБ] —

г I х

1

хд/х2 - г"

Дифференцируя выражения (23) по г, получим:

6^х)

х"х f

_6Г(г) 6г

6х

х6х

V х2 - г 2 л/Т-Т

Если исходить из преобразования Абеля в виде выражения (18), то вместо (25) можно получить другой вид производной от преобразования Абеля. Выполняя интегрирование по частям в уравнении (18), получим:

Б(г) = 2Д1)л/1 - г 2 - 2^

2 2 6^

х - г — 6х . 6х

Дифференцируя уравнение (26) по г, получим:

6Р

'1 V

1

6Г(г ) 6г

2г

6х

-6х

Г2 2 1

г ух - г \1 "

Производная от преобразования Абеля в виде выражения (25) будет использована при получении решения для обычной линзы Лунеберга [ 1].

Решение интегрального уравнения для обобщенной линзы Лунеберга

Будем исходить из уравнения Флетчера [2]. Перепишем уравнение (7) в другом виде:

61п п

1 д/п 2 (г)г 2 - п 2

= ОД

2

2

г

г

2

г

г

1 ( h где f(h) =—arcsin I —

2h I f

Введем обозначения:

lnn(r) = F(r), n(r)r = p, F(r) = F(p). Тогда вместо (28) запишем:

dF

if dp 77 J

dp

1

/—2-= -f(h) .

yjp - h2 T

С учетом того, что Б(1) = 1пп(1) = 1п1 = 0, уравнение (29) с точностью до обозначений совпадает с уравнением (19), и поэтому для его обращения можно использовать уравнение (18). Тогда получим:

И

F(p)=П f

, arcsin

1 1 f

dh

y[h2

-P

Из уравнения (30) для показателя преломления обобщенной линзы Лунеберга, представляющей собой радиально-симметричную градиентную сферу, фокусирующую параллельный пучок лучей в точку (рис. 1), будем иметь [6, 7]:

n(p) = exp

П f

, arcsinI — |dx

1 1 f

np v*2-p2

где р = п(г)г .

Решение для обычной линзы Лунеберга с помощью производной от преобразования Абеля Интеграл в уравнении (31) не берется в аналитических функциях даже при / = 1 (случай обычной линзы Лунеберга [1]). Однако производная от функции показателя преломления (31) может быть записана в элементарных функциях. Действительно, используя выражение для производной для преобразования Абеля (25), получим:

С | . х |, .1 агсБт— /

dn(p) dp "

n

Tip

n

Tp

n(p) np

1 f

d* I f

dx

<Jx

p

p

aicsin— d* f

-p2)(( -x2) Vbp 1 . (p2 + f -2 '

T

---arcsin

4 2 I-2

arcsin— f

p - f2

1

p

Заметим, что уравнение (32) в [4] получено в другом виде. Далее используя уравнение (32) можно получить решение для классической (обычной) линзы Лунеберга. Для этого в уравнении (32) положим /= 1 (фокусировка параллельного пучка лучей на по-

верхности линзы), тогда вместо уравнения (32) получим:

dln n(p) = 1 dp 2p

1 -

1

Интегрируя уравнение (33), получим:

ln n(p) = ln

щ

1 - p , p = n(r)r .

При выводе уравнения (34) использовался табличный интеграл [ 10]:

dx 1 ,

= = — ln

2 a

xyja2 -

x

V2 2 a - x

Из уравнения (34) нетрудно получить функцию показателя преломления для линзы Лунеберга единичного радиуса [ 1]:

n2(r) = 2 - r 2 .

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

В [8] Микаэляном получено и решено интегральное уравнение для бесконечного цилиндрического ГОЭ, собирающего все лучи, исходящие из точечного источника на оптической оси, в точки на оптической оси, расположенные периодически. Для конечного ГОЭ длиной Ь задачу Микаэляна можно сформулировать как задачу Лунеберга для сферического ГОЭ. Требуется найти п = п(г), г - радиальная координата в цилиндрической системе координат, который фокусирует параллельный поток лучей, падающих перпендикулярно на поверхность цилиндрического ГОЭ, в точку на оси на противоположной стороне ГОЭ (см. рис. 3).

i r

k

h г \ z

0 n(r) L

Рис. 3. Ход луча в цилиндрической линзе Микаэляна

Для цилиндрической задачи лучевой инвариант, аналогичный выражению (1), имеет вид закона преломления Снелиуса [9]:

n

где h - расстояние от оси рассматриваемого луча при z = 0.

В уравнении (37), которое выполняется вдоль луча, учтено, что при z = 0 угол в для всех лучей равен нулю. Из рис. 4 видно, что имеет место связь между приращениями:

x

p

dr = I 1 dz V cos2 в

-1 = ¡n^l -1

' n2

Рис. 4. Приращение вдоль траектории луча Из уравнения (38) следует интегральное уравнение относительно п(г) для линзы Микаэляна:

= -|6г = -L .

II

n 2

n 2

-1

В [8] Микаэлян подбором нашел решение уравнения (39) и заметил при этом, что уравнение (39) нельзя отнести к определенному типу интегральных уравнений. Но это не совсем так. Покажем, что уравнение (39) можно решить с помощью преобразования Абеля.

Перепишем уравнение (39) в виде:

i

dr

Jn2(r )-

р

L

р

где р = п(И). И введем обозначения: п(г) = х, г = m( х). Тогда уравнение (40) можно переписать в виде:

6т

, п(0) '

П I

ГГ •>

-dx

i

dx

22 x - р

L

Пр

С учетом того, что m(n(0)) = г = 0, уравнение (41) отличается от уравнения (19) только обозначением. Тогда, применяя обратное преобразование (18) к уравнению (41), получим:

ар

2L n

m(x) =— f -=-

П x ^р2-x2

2Ll = — ln

П

2L l

= ln

П

р + у[р2 - x2

n

n(0) + д/n 2(0) - x 2

Из уравнения (42) следует:

n

= ex

Из уравнения (43) нетрудно получить решение Микаэляна [8] :

n

n(r) = -

chl п

\ 2L

Интегральное уравнение для обобщенного цилиндрического ГОЭ

Как решения Моргана [3], Флетчера [2] и Фло-реса [4-6] обобщают линзу Лунеберга [1], также можно попытаться обобщить решение Микаэляна [8]. В [15], например, рассмотрены составные цилиндрические ГОЭ, являющиеся обобщением линзы Микаэляна и предназначенные для фокусировки параллельного пучка лучей в произвольное радиально-симметричное распределение интенсивности на задней поверхности ГОЭ.

Получим уравнение для обобщенной линзы Микаэляна. Требуется найти функцию n = n(r) для цилиндрического ГОЭ, фокусирующего падающий параллельно оси z пучок лучей в точку на оси за оптическим элементом на расстоянии f (см. рис. 5).

h , r

n

0 L L+f

Рис. 5. Ход луча в обобщенной линзе Микаэляна

Рассуждая аналогично, вместо уравнений (39) и (40) получим уравнение для обобщенной линзы Микаэляна:

dr

Wn2

L n

где И - радиальная координата луча во входной плоскости, И] - в выходной плоскости (см. рис. 5).

Из рис. 5 также можно получить еще одно уравнение, связывающее координаты И и И :

n(h )sine(hj ) = sin у =

h

f+hh2

п(И1 )81п0(И1) = п(И) .

Уравнение (46) - это закон преломления Сне-лиуса, примененный к выходной поверхности, а уравнение (47) - это лучевой инвариант (37). Исключая из уравнений (46) и (47) угол 0, получим:

n2(hx) - n 2(h) =-2^ h2 + f

С помощью замены переменных: п(г) = х, п(И) = р, п(И1) = р1, г = т(х), уравнение (45) можно привести к уравнению Абеля: 6т (х) ,

= L

22 руx -р

р

р

x

x

где:

где:

Р1 =

р2 +

m1

m12 + /2

m1 = m(р1 ) = И1 .

В отличие от стандартного уравнения Абеля (19) в уравнении (49) верхний предел интегрирования не является постоянной величиной, а зависит от нижнего предела интегрирования, поэтому применять для решения уравнения (49) стандартную формулу обращения преобразования Абеля (18) нельзя.

Обращение преобразования Абеля с двумя переменными пределами интегрирования

Общий случай

Попытаемся с помощью преобразования Фурье найти формулу для обращения преобразования Абеля с двумя переменными пределами интегрирования. Будем исходить из преобразования Абеля вида:

Г(р) = 21

Г(х)хСх

22 рл/ х -р

где р = р (р). Введем замену переменных: х = у,

р2 = в, рр =Д, тогда уравнение (51) перепишется в виде интеграла свертки:

ад=Г-Су =ОТ

НУ-в -От

гес/^ ГУСу

где: Г(в) = Г(р), Г(у) = Г(х), В =в-2в1, А = —^

Заменим в уравнении (52) все функции на их Фурье-образы:

ф(ю) = 1 Б(в)еСв,

ОТ ОТ

/(у) = 1 Б(#)е-гу Л Л

гей У-А 1 = 2В1япс(Вп)ег(у-А)п Сп,

В

1ОТ

е

1

г( У-в

д/|У -в ^^-ОТ АМ

-С.

Тогда вместо уравнения (52) получим:

ф|) = 1 1 8(#)[2В 81пс(ВП)]е,у(^+п+!)е-гАп.

• е-гве-1в<я А&п&СвСу = =1118<£^П£1^£

2п 7

= я 8(-п-^ + ю) апС!

д(п,7 +Ю)=1[2В(Р)Бтс(5(Р)>1]^

гА

ср,

б - дельта функция Дирака.

Частный случай

Уравнение (57) в общем виде не позволяет воспользоваться преимуществом интеграла свертки, но это можно сделать для частного случая зависимости функций В(в) и А(в) от в.

Предположим, что верхний и нижний пределы в первом интеграле (52) связаны простой формулой:

в = в + 2М ,

тогда В(в) = М, А(в) = в+ М и интеграл (58) вычисляется в явном виде:

+ ю) =

= 2М вшсМ п е-Мг| 1 е-'в(п+!+ш) Св =

-ОТ

= [2М вшсМ п]е-Мп 8(л +7 + ю). Тогда уравнение (57) можно переписать: „ч, ч пЧ-п-Ф-МбшСМп]

ф(ю)=я—дг^ •

• е-гМп 8(а + п +!)пС =

Л [2М¡апсМп]гМп , = 8(ю) 1 --, =-ап =

= Б(ю)и(ю, М).

Уравнение (61) позволяет легко найти функцию Б(ю), а также, с помощью преобразования (54), и искомую функцию /(у):

г( у) = 1

ф

егую1,

где:

М ОТ [2МБ.псМп]е-гМп с и(ю, М) = 1 --, 1 -Сп =

=е

¡Мю

М 2

1

п + ю

Л

3м

2 _ ей

,-Л 1 ^С!-

VI л + М| М 7171

_Мд/ Л+ М| м

2 2

Интеграл (63) можно свести к табличным интегралам Френеля (заменой переменной 7 = р2):

I» 2

и(ю,М) = 2 1 егюр Ср = 2И1(ю,М) =

е

ОТ

-ОТ

ОТ

1

где табличный интеграл Френеля имеет вид: C(x) = JC0S[§ t2 jjdt,

S (x) = j sinfc t2 jjdt.

f(x) = _fj F(p)e- _ pdpda.

В исходных обозначениях выражение (62) для обращения преобразования Абеля (51) при условии (59) имеет вид:

е-г

е , ,

Интегральное уравнение для цилиндрического градиентного аксикона В отличии от обобщенной цилиндрической линзы, показатель преломления которой ищется из уравнений (47) и (48), для аксикона, который все лучи параллельные оптической оси и отстоящие от неё на расстоянии И, собирает в точке на оптической оси, отстоящей от второй поверхности ГОЭ на расстояние _ДИ) (см. рис. 6).

h _> , r

n

0 L fh

Рис. 6. Ход луча в градиентном аксиконе Так как все лучи, выходящие из ГОЭ аксикона, должны иметь одинаковый угол Т, то вместо уравнения (46) для аксикона получим:

n2(hl) - n2(h) = Sin2T .

Тогда пара уравнений для нахождения показателя преломления цилиндрического аксикона будет иметь вид:

6т(х) , д ——6х Т

I 16с _ =-Т,

22 руx -р

Р

р2 - р2 = sin2T .

Чтобы воспользоваться формулой (66) для решения уравнений (68) и (69), необходимо использовать связь обозначений

dm

= 2f (x) x,

F( p) = --,

р

sin2 T

M =

Тогда формула обращения преобразования Абеля (66) для аксикона будет иметь вид:

6т

2х6х ^ Щю, М) . (

Интегрируя уравнение (71), получим:

» п(0) . п(0) .

т(х) =-Т I I егйу 6у I е~Шр 6р =

' П Сл, ЧЛ ' 7 3

Uj

= -L j da

(eia"2(0)- e'**2) V1(a, n(0)) ia U1(a, T) '

где:

V3 sin T

U1 (a, T) = j e'ap dp

Sin T 2

n

V1 (a,n(0)) = j e~iap dp.

Функции и1 и VI также как и функция и(, М), сводятся к интегралу Френеля. Переходя к исходным обозначениям, окончательно вместо (72) получим:

=iL j

da

a

ian2

Vx

Uj

.

Решим уравнение (75) для «слабого» аксикона. Пусть п 2(г) = п 2(0) -Ап 2(г), где Ап 2(г) << п 2(0). Тогда уравнение (75) можно упростить:

ад

- Т = 116е'"2

где:

R

U1

Из уравнения (76) следует, что: n 2 (r) = n 2(0) + r

LW (n(0), T) n(r) n 2(0) + [[ (n(0), T)]

где:

W

-ад

Литература

1. R.K. Luneberg Mathematical Theory of Optics // Brown U. Press, Providence, R.I., 1944.

2. A. Fletcher, T. Murphy, A. Young Solution of two optical problems // Proc. R. Soc. London Ser. A, v. 223, p. 216-225, 1954.

r

z

x

3. S.P. Morgan General solution of the Luneberg lens problem // J. Appl. Phys., 1958, v. 29, p. 13581368.

4. J.R. Flores, J. Sochacki, M. Sochacki, R. Staronski Quasi-analytical ray tracing through the generalized Luneberg lens // Appl. Opt., 1992, v. 31, no. 25, p. 5167-5170.

5. J. Sochacki, J.R. Flores, C. Gomex-Reino New methods for designing the stigmatically imaging gradient-index lenses of spherical symmetry // Appl. Opt., 1992, v. 31, no. 25, p. 5178-5183.

6. J.R. Flores Gradient-index with spherical symmetry // J. of Modern Opt., 1999, U. 46. №. 11. Р. 1513-1525.

7. J.M. Gordon Spherical gradient-index lenses as perfect imaging and maximum power transfer devices // Appl. Opt., 2000, v. 39, no. 22, p. 3825-3832.

8. Микаэлян А.Л. Применение слоистой среды для фокусирования волн // Доклады академии наук СССР, 1951. Т. LXXXI. С. 569-571.

9. G.I. Greishuk, S.T. Bobrov, S.A. Stepanov Optics of diffractive and gradient-index elements and systems // SPIE Press, Bellingham, 1997.

10. Прудников А.П., Брычков Ю.А., Марычев О.И. Интегралы и ряды // М., Наука, 1981.

11. Гельфанд И.М., Гиндикин С.Г., Граев М.И. Избранные задачи интегральной геометрии // М., Добросвет, 2000.

12. Оптическая голография // Под ред. Г. Колфилда, М., Мир, 1982, т. 1.

13. Диткин В.А., Прудников А.П. Интегральные преобразования и операционное исчисление // М., Физматлит, 1961.

14. Корн Г., Корн Т. Справочник по математике // М., Наука, 1968.

15. Котляр В.В., Мелёхин А.С. Расчёт составного градиентного оптического элемента формирующего заданное распределение интенсивно-сти // Компьютерная оптика, ИСОИ РАН, Самара, 2001, вып. 21, с. 92-95.