CC BY
40
ОПТИМИЗАЦИЯ ДОЭ, ФОКУСИРУЮЩЕГО ГАУССОВ ПУЧОК В ПРЯМОУГОЛЬНУЮ ФОКАЛЬНУЮ ОБЛАСТЬ
В. С. Павельев Институт систем обработки изображений РАН Самарский государственный аэрокосмический университет
Аннотация
Данная работа посвящена применению известной генетической итерационной процедуры нахождения экстремума функции многих переменных к расчету квантованного ДОЭ, фокусирующего Гауссов пучок в прямоугольник.
Введение
В [1] показано, что возможность сведения вычисления двумерного оператора распространения света к вычислению одномерного в радиально-симметричном случае позволяет успешно применять стохастические методы к оптимизации квантованной фазы радиально-симметричных ДОЭ (дифракционных оптических элементов).
Если фазовая функция элемента <р(и, V) и амплитуда в сечении освещающего пучка А0 (и,V) = ^10(и,V) допускают разделение по переменным, т.е.
где
\ф(и,v) = (и) + (v)
[ A0 (и, v) = Ai (и )A2 (v) '
(1)
вычисление двумерного оператора сводится к вычислению двух одномерных и, таким образом, для оптимизации фазовой функции в этом случае также целесообразно использовать стохастические методы. Данная работа посвящена применению такого подхода к расчету квантованного ДОЭ, фокусирующего гауссов пучок в прямоугольную фокальную область. Актуальность такой задачи объясняется широким спектром потенциальных приложений [2], а также тем, что значительное количество технологий изготовления ДОЭ связано с квантованием его фазовой функции [3] и уменьшение числа уровней квантования значительно удешевляет процесс изготовления.
Постановказадачи
Рассмотрим задачу фокусировки пучка с распределением интенсивности в сечении
10 (и, V) = ехр (-2(и2 + V2 )/а2 ) =
= ехр (-2м2 / а2) ехр ((V2 / а2)
в прямоугольник с фокальными размерами сторон 2а и 2Ь с помощью ДОЭ с фазовой функцией <р(и, V), определенной на апертуре Б размером
2А х 2А (рис. 1).
Определим заданную фокальную интенсивность следующим образом:
C , -а < х < a
Ix = \ х •
0, иначе
I =\СУ -b <y<b
y ' 0, иначе
(3)
Cx, Cy -константы, выбираемые из соображений сохранения энергии:
1 A
Cx = — J exp(-2u2/а2)du,
2а - A
Су = 2т J exp(-2v2/а2)) .
(4)
(5)
I0(u,v)
I(x,y)
V
Рис. 1. Постановка задачи фокусировки гауссова пучка в прямоугольную фокальную область
В случае, если ДОЭ с фазовой функцией ((м,V)
помещен в заднюю плоскость Фурье-каскада с фокусным расстоянием линзы f, сформированная
интенсивность /(х, у ) в выходной плоскости может
быть найдена следующим образом [4]:
ik
I (x,y ) =
:JJ exp (-(u2 + v2)/а2 )exp (((u, v ))x
2nfD
x exp (-ik (xu + yv ))f )dudv| . (6)
В случае расчета ДОЭ, включающего в себя линзу: ik
I (x,y ) =
:JJ exp (-(u2 + v2)/а2 )exp (((u, v ))x
I = Ix•I„
(2)
2nfD
x exp (ik (u2 + v2 )/2 f -ik (xu + yv))f )dudv| . (7)
Учитывая (1), можно переписать формулу (7) в следующем виде:
1{х,у) = ~* (х). !у (),
(8)
где
I (х) =
л
фк/2п/ | ехр (-и 7ст2 )>
<ехр((¡(и))ехр(ки2/2/- 1кхи//)с"| , (9)
^ (У) =
фк/2п/ | ехр(-V1/ст2)
х ехр (<р2 (V)) ехр ((/2/ - 1куу / /) с) В этом случае функционал 1 а
ф((р(и,V)) = — { У*(*,((")) -С*|с* +
(10)
+¿1 к; (у-^2 (V)) - Су^у
(11)
характеризует качество формирования заданного распределения интенсивности с помощью ДОЭ с фазовой функцией (р(и,V) = ((и) + (2(у) .
Построение процедуры оптимизации квантованной фазовой функции ДОЭ
Перейдя к дискретной модели с числом элементов дискретизации в плоскости ДОЭ и шагом дискретизации в фокальной плоскости, равным размеру дифракционного пятна Д/ = Л//(2А), можем заменить (11) функционалом
1 нх-1
Ф(() = Т- Е \!х(*„,((и))-Сх\ +
N * п=о
1
Ту -1
+^ Е у> (у» (00) - Су
Т У т=0
(12)
где = 2а/Д/ хп = -а + пДу, ут = -Ь + тД/,
Ту = 2Ь/Д / ■
I (хп ) =
N-1
Е С"п ехр (((ир ))
р=0
N-1
Е С» ехр ((2 (V, ))
1у (у» ) =
ир = - А + рД, Vч = - А + ,Д,
"р+1
Ср =| Щ2к/ехр (-и2!ст2)
"р
х ехр (-/ кхпи!/ +1ки2 /2/) с1и ,
(13)
(14)
(15)
С» = | фк/2п/ ехр (-V2!ст2)>
с ехр
(-1кутV/ + кV /2 / ) СУ.,
Д = 2 А / N ■ Отметим, что коэффициенты С"р, С» не зависят от выбора фазовых функций и могут
быть посчитаны до начала работы оптимизационной процедуры, впрочем, также как и коэффициенты
С* , Су .
Для минимизации функционала (12) в данной работе использовалась модификация известной итерационной генетической процедуры [1, 5].
Перед началом работы итерационной процедуры формировалась начальная популяция, представлявшая собой набор Nёеп сгенерированных случайным образом массивов
>2 N -1 >]=0
{( }2Л'0-1,( £ {(2пт /М) : т е [1,М]} . (17)
Каждый массив содержал последовательно записанные отсчеты функций ( (и) и (2 (у):
= ((-А + уД) ] < N
( ((-А +(-N)Д)> N
(18)
Для увеличения скорости сходимости в начальную популяцию вводилось в качестве начального приближения (т.е. одного из «индивидов» нулевой популяции) геометрооптическое решение ((и, V) задачи фокусировки гауссова пучка в прямоугольник, квантованное по М уровням. Такое решение представляет собой суперпозицию двух одномерных фокусаторов в отрезки с длинами 2а и 2Ь:
((и,у) = ( (и ) + ( (V) . (19)
Далее, в цикле последовательно выполняются следующие операции:
- скрещевание - формирование дополнительных новых индивидов-массивов, значение каждого отсчета которых выбирается как значение соответствующего отсчета случайно выбранного массива-«родителя»;
- мутация - стохастическое изменение значений случайно выбранных отсчетов, получившихся новых массивов. Выбор вероятности мутации определяет скорость сходимости алгоритма;
- селекция - определение следующей популяции «родителей» путем выбора Nёеп лучших (в
смысле решаемой задачи) индивидов с помощью вычисления целевых функций для всех индивидов.
Для улучшения сходимости генетическая процедура совмещалась с применением метода вариаций на каждой итерации.
После завершения итерационной процедуры, следующие критерии были выбраны для оценки качества фокусировки: энергетическая эффективность
р
2
«= ZZ Ix (x (ym VEZ cxC,, (20)
и=0 m=0 / и=0 m=0
и среднеквадратичное отклонение сформированной интенсивности от среднего значения 1
( Nt -1 Ny-1
HZZ[ 1 (Хп Ж (Уm )- 1 \
( Nx-1 Ny-1
Y
zz[ h (x. )i; (y.)]
^ .= 0 m=0
_ Nx-1 Ny -1
1 = ZZ 4 (Х. )y (Ут )/NxNy •
(21) (22)
Результаты вычислительного эксперимента Для постановки вычислительного эксперимента выбирались следующие параметры: длина волны освещающего пучка Х=10,6 мкм, радиус освещающего пучка ст=1,85 мм, фокусное расстояние 7=100 мм, размер апертуры 2*А=6,4 мм, размер фокальной области 2а*2Ь=1,75*3,5 мм, число отсчетов фазовой функции на апертуре элемента Л*Л=100*100, число уровней квантования М=8, количество итераций стохастической процедуры Лй=100, число индивидов исходной популяции N¡,¡¡„=20, число индивидов, образованных операцией скрещивания Л/сЫа=20.
Уменьшение значения функционала невязки с ходом итераций показано на рис. 2.
0,16
0,14
0,12
61 81 101 номер итерации
Рис. 2 Зависимость значения невязки от номера итерации
На рис. 3 представлены итерационно рассчитанная фаза ДОЭ (а) и интенсивность в фокальной плоскости элемента: (б) - результат численного моделирования итерационно рассчитанного ДОЭ, (в) -фаза геометрооптического приближения, (г) - результат численного моделирования геометроопти-ческого приближения.
Значения расчетной энергетической эффективности (20) и среднеквадратичного отклонения интенсивности (21), полученные в ходе вычислительного эксперимента для квантованного геометрооп-тического приближения и итерационно найденного решения, приведены в таб. 1.
Рис. 3. Итерационно рассчитанная фаза ДОЭ (а),
(б) - распределение интенсивности, сформированное итерационно рассчитанным ДОЭ
в фокальной плоскости (результат численного моделирования), (в) - фаза геометрооптического приближения, (г) - распределение интенсивности, сформированное геометрооптическим приближением
a
б
в
г
Таблица 1
Метод расчета e, % S,%
Геометрооптиче-ский 87 45
Итерационный 88 13
Таким образом, применение стохастической итерационной процедуры с функционалом (12) позволяет рассчитывать восьмиуровневые ДОЭ, фокусирующие гауссов пучок в прямоугольник с энергетической эффективностью, близкой к эффективности геометрооп-тического решения, и более чем в 3 раза меньшим значением среднеквадратичного отклонения для выбранного набора физических параметров. Отметим также небольшие вычислительные затраты, требуемые для реализации расчета элемента, что объясняется простотой вычисления функционала (12) после того, как значения коэффициентов (15) предварительно посчитаны.
Заключение
В данной работе рассмотрено применение известной генетической итерационной процедуры нахождения экстремума функции многих переменных к расчету квантованного ДОЭ, предназначенного для фокусировки гауссова пучка в прямоугольник.
Приведенные результаты вычислительных экспериментов подтверждают перспективность рассматриваемого подхода.
Благодарности
Автор выражает благодарность С.Г. Волотов-
скому за помощь в постановке вычислительного
эксперимента.
Литература
1. Павельев В.С. К расчету квантованных радиаль-но-симметричных ДОЭ // Компьютерная оптика, 2001. №22 (принято к публикации).
2. M. Duparre', M.A. Golub, B. Luedge, V.S. Pave-lyev, V.A. Soifer, G.V. Uspleniev, S.G. Volotovskii Investigation of computer-generated diffractive beam shapers for flattening of single-modal CO2-laser beams // Applied Optics, 1995, v. 34, N 14, р. 2489-2497.
3. A.V. Volkov, N.L. Kazansky, V.S. Pavelyev, V.A. Soifer, G.V. Uspleniev Methods For Computer Design Of Diffractive Optical Elements // Chapter 4. Technology of DOE fabrication. In Victor A. Soifer (Ed.): A Wiley Interscience Publication. John Wiley&Sons, Inc., New York, 2002, p. 267-345.
4. V.A. Soifer Methods For Computer Design Of Dif-fractive Optical Elements // Chapter 1. Introduction to Diffractive Optics. In Victor A. Soifer (Ed.): A Wiley Interscience Publication. John Wiley&Sons, Inc., New York, 2002, p. 1-54.
5. C.Y. Lu, H.Z. Liao, C.K. Lee, J.S. Wang Energy control by linking individual patterns to self-repeating diffractive optical elements // Applied Optics, 36, 1997, N.20, p. 4702-4712.
