CC BY
42
УДК 353.42
СТОХАСТИЧЕСКИЙ ПОДХОД К ОПТИМИЗАЦИИ КВАНТОВАННЫХ ДИФРАКЦИОННЫХ ОПТИЧЕСКИХ ЭЛЕМЕНТОВ
© 2002 B.C. Павельев
Институт систем обработки изображений РАН, г. Самара
Работа посвящена применению стохастических методов оптимизации к расчету дифракционных оптических элементов с квантованной фазовой функцией.
Введение
В настоящее время дифракционные оптические элементы (ДОЭ) нашли широкое применение для решения задач управления пучками лазерного излучения в таких областях как оптическое приборостроение, технологические лазерные установки, волоконная оптика и т.д. ДОЭ представляет из себя пластинку из оптического материала с нанесенным микрорельефом, высота которого определяет локальное изменение фазы освещающего пучка [1].
Значительное количество современных технологий формирования микрорельефа использует аппроксимацию в общем случае непрерывной функции микрорельефа ДОЭ кусочно-постоянной (квантованной) функцией [2]. Это определяет актуальность построения процедур расчета квантованных ДОЭ (ДОЭ с квантованной функцией микрорельефа). Расчет квантованных ДОЭ с помощью итерационных алгоритмов типа Фьенапа [3, 4] затруднен в силу необходимости решения обратной задачи на каждой итерации. Использование стохастических процедур [5] не предполагает необходимости решения обратной задачи, т.е. поиск функции микрорельефа ДОЭ может осуществляться непосредственно над множеством технологически реализуемых функций, что дает возможность избежать погрешности квантования.
Однако построение стохастической процедуры оптимизации реального ДОЭ требует решения большого числа прямых задач, что в общем случае может привести к нереализуемым на практике вычислительным затратам.
Представляется целесообразным рассмотреть применение стохастических проце-
дур оптимизации ДОЭ специально для случаев, когда решение прямой задачи не требует больших вычислительных затрат (например, в случае радиально-симметричного ДОЭ или в случае если функция пропускания ДОЭ допускает разделение по координатам). Данная работа посвящена применению известной генетической стохастической процедуры нахождения экстремума функции многих переменных к расчету квантованных ДОЭ, фокусирующих излучение одномодового лазера в двумерные фокальные области в виде круга и прямоугольника.
Оптимизация фазы радиальносимметричного ДОЭ с помощью генетического алгоритма
Генетический алгоритм (ГА) - это итерационный стохастический алгоритм, который работает с набором "индивидов" (популяцией). Каждый индивид представляет собой потенциальное решение задачи. К примеру, в случае расчета радиально-симметричного ДОЭ в качестве отдельного индивида может быть определена фазовая функция
ДОЭ ф(г)е [0,2л] такого вида, что соответствующий микрорельеф
к(г)=ф(гЦ(2л(п -1)), (1)
где п - показатель преломления материала, 1 - длина волны излучения, реализуем с помощью имеющихся в наличии исследователя технологических возможностей. В случае наличия технологии изготовления бинарных (двухуровневых) элементов в качестве индивидов целесообразно рассматривать фазовые функции, принимающие в каждой точке одно
из двух значений - 0 или р (или 2р и р). В общем случае расчета ДОЭ с числом уровней квантования М, фазовая функция принимает значения (р(г)е {рт / М): т е [1, М ] Отметим, что если в случае использования "детерминированных" методов оптимизации [3, 4] расчет квантованных ДОЭ с малым числом уровней квантования М (М=2-8) требует дополнительных усложненных модификаций процедуры расчета, то для стохастической оптимизации уменьшение числа уровней квантования означает уменьшение мощности множества допустимых решений, и, следовательно, снижение вычислительных затрат. Перед началом работы алгоритма популяция генерируется случайным образом. Затем каждому индивиду присваивается при помощи некоторого "оценивающего" (целевого) функционала Ф(ф), мера ценности индивида ф в отношении рассматриваемой задачи. Значение Ф(ф) используется алгоритмом для направления поиска. Конкретный вид функционала в нашем случае определяется приложением ДОЭ. На каждой итерации процедуры численность популяции искусственно увеличивается за счет определенным образом вводимой операции "скрещивания" индивидов исходной популяции. Отобранные лучшие индивиды (в смысле введенного функционала Ф(ф)) используются в качестве исходной популяции для следующей итерации алгоритма. Вид функционала Ф(ф(г)) определяется конкретным приложением ДОЭ и его вычисление вовсе не обязательно требует решения прямой задачи теории дифракции. К примеру, если требуется рассчитать элемент, формирующий пучок с заданным модовым распределением, вычисление функционала Ф(ф(г)) потребует разложения светового поля по модам в плоскости непосредственно за оптическим элементом. Рассмотрим подробнее построение генетической процедуры оптимизации радиально-симметричного квантованного ДОЭ радиуса Я. Начальная популяция представляет собой набор
Nгеп сгенерированных случайным образом
одномерных массивов jj}N= ,0 . Каждый массив соответствует одной реализации фазы элемента j(r), взятой вдоль его радиуса:
jj = j(j'A)e{2pm / M }, (2)
где А = R/(N -1), R - радиус ДОЭ.
После генерации начальной популяции (популяции "родителей") организуется цикл, в теле которого последовательно выполняются следующие операции:
- скрещевание - формирование Nchild дополнительных новых индивидов-массивов, значение каждого отсчета которых выбирается как значение соответствующего отсчета случайно выбранного массива-"родителя";
- мутация - стохастическое изменение значений случайно выбранных отсчетов получившихся новых массивов. Значение вероятности мутации p определяет скорость сходимости процедуры;
- селекция - определение следующей популяции "родителей" путем выбора Ngen лучших (в смысле решаемой задачи) индивидов с помощью вычисления функционала F(j)
для всех Ngen + Nchild индивидов.
Таким образом, на каждой итерации необходимо найти (Ngen + Nchild ) значений функционала Ф(<р). Количество итераций N , требуемых для оптимизации фазы ДОЭ в данной работе выбиралось исходя из количества итераций, требуемых для оптимизации тестового функционала, для которого было известно значение глобального экстремума.
B качестве примера, иллюстрирующего работоспособность предлагаемого подхода, был выбран расчет ДОЭ, фокусирующего гауссов пучок в круг, что объясняется широким спектром потенциальных применений таких элементов. Выбирая функционал для оптимизации фазы ДОЭ, необходимо учесть следующие требования:
- адекватность функционала физической задаче,
- минимизация вычислительных затрат, требуемых для его реализации.
Рассмотрим задачу синтеза ДОЭ с радиусом апертуры я для формирования из плоского гауссова пучка
10 (г)= ехр(- 2г Vо2) (3)
фокального распределения интенсивности (рис.1)
і (р)--
\1,, 0 <р< я, I 0, р > я,
я,
с,
N. -1
фр(г ))=£ I1 (Р1 )- 1,
I=0
Рис.1. Постановка задачи синтеза ДОЭ для фокусировки в круг
(р)=Х
і=0
(4)
Я1 - радиус фокальной области. Значение
где
интенсивности 11 может быть получено из соображения сохранения энергии:
цией плоской волны на кольцеобразной апер туре [4]:
ехР(- 0+1/ а2 )ехр(ф}+1 )х
х(+1 ( 1 (Р(+11/)/р-- -1 (/7/)/р)
(7)
= ’А , J1 (•)- Бессель-функция 1го порядка. Следовательно, значения интенсивности в точках р1 могут быть найдены следующим образом:
(5)
(6)
1 (г )=
или
ехр(- г^1о )(Рі)х
х(+1 ((кРігі(,!/)/Р,--(кРігі !рР /
,м константа, описывающая потери гд:ергии, связанные с квантованием фазы по м уровням. Элемент рассчитывается как добавка к внешней сферической линзе.
В качестве минимизируемого функционала ф(р(г)) логично выбрать невязку:
і (р,)--
N-2 .
X ехР(р (, Рі
і=0
(8)
(9)
где N, = Я, /А, - радиальный размер фокальной области в дифракционных пятнах, А, = 1//2Я - размер дифракционного пятна,
р1 = (, ( 0,5),, 1 (р) - фокальное распределение интенсивности, сформированное ДОЭ с фазовой функцией р(г). Комплексная амплитуда светового поля w(р) в фокальной плоскости линзы может быть представлена как сумма амплитуд, образованных дифрак-
где
С’ = ехр(- (+1/ а" )х
х (г/+1 ( (/р1г’+11/)/р1 -( J1 (/ГГ’//)/р)^ (10) Отметим, что коэффициенты С’ не зависят от ф(г ) и могут быть посчитаны до начала работы оптимизационной процедуры. После завершения работы оптимизационной процедуры подсчитывалось значение энергетической эффективности
Я /я
е =11 (р)р^р \1о(г(11)
0 / 0
и среднеквадратичного отклонения интенсивности от среднего значения I
5= и ^ 11 (р)-112рсір , (12)
1 ЦЯ
і о
1 =
Я
11 (р)р&р .
(13)
Для постановки вычислительного эксперимента выбирались следующие параметры: длина волны освещающего пучка 1=10,6 мкм, радиус освещающего пучка а=1,85 мм, фокусное расстояние внешней линзы /=100 мм, радиус апертуры Я=3,7 мм, радиус фокальной
1(9)
Л
о
10 20 30 40
р/Дг
50
Рис.2. Радиальное распределение интенсивности, сформированное рассчитанным
ДОЭ в фокальной плоскости (результат численного моделирования)
области Я1=0,6 мм, число отсчетов фазовой функции N=100, число уровней квантования М=4, количество итераций стохастической процедуры N=6000. На рис.2 представлена расчетная интенсивность в фокальной плоскости элемента. Значение расчетной энергетической эффективности в ходе вычислительного эксперимента составило е = 80,3%, значение среднеквадратичного отклонения интенсивности составило 8 = 8%.
Таким образом, применение стохастической итерационной процедуры с функционалом (6) позволяет рассчитывать четырехуровневые ДОЭ, фокусирующие гауссов пучок в круг с энергетической эффективностью, близкой к теоретическому пределу (81%), и низким значением среднеквадратичного отклонения. Отметим также небольшие вычислительные затраты, требуемые для реализации расчета элемента, что объясняется простотой вычисления функционала (6) после того, как значения коэффициентов (10) предварительно посчитаны.
Оптимизация фазы ДОЭ, фокусирующего гауссов пучок в прямоугольную область
Отметим, что если фазовая функция элемента ср(и,V) и комплексная амплитуда в сечении освещающего пучка А0 (и, V) допускают разделение по переменным, т.е.
\ф(и,v )= ф ( )+ ф ( )
1 Ао (и v)= А1(и)А2 (v) (14)
вычисление двумерного оператора сводится
к вычислению двух одномерных [4] и, таким образом, для оптимизации фазовой функции в этом случае также целесообразно использовать стохастические методы. Рассмотрим применение такого подхода к расчету квантованного ДОЭ, фокусирующего гауссов пучок в прямоугольную фокальную область. Актуальность такой задачи объясняется широким спектром потенциальных приложений [6].
Рассмотрим задачу фокусировки пучка с распределением интенсивности в сечении
10 (и, V)= ехр(- 2(и2 ( V2 )/а2 ) =
= ехр(- 2и2 / а2 )ехр(- 2v2 / а2) в прямоугольник с фокальными размерами сторон 2а и 2Ь с помощью ДОЭ с фазовой функцией Р(и, V), определенной на апертуре
В размером 2А х 2А (рис.3).
Определим фокальную интенсивность следующим образом:
где
1
1 = 1х • 1у ,
\СХ, -а < х < а [ 0, иначе
Гс„
(15)
(16)
(17)
-Ь < у <Ъ
[ 0, иначе ’
С^ Су-константы, выбираемые из соображений сохранения энергия:
1 а
Сх = — 1 ехр(- 2и2 / а2 )и,
2а - А
(18)
10(и,у)
1(х,у)
Рис.3. Постановка задачи фокусировки в прямоугольную фокальную область
1 A
Cy = — I" exp(- 2v2 / о2 ) y 2b V ^ .
(19)
В случае, если ДОЭ с фазовой функцией р(и, V ) помещен в заднюю плоскость Фурье-каскада с фокусным расстоянием линзы /, сформированная интенсивность I (х, у) в
выходной плоскости может быть найдена как квадрат амплитуды Фурье-образа входной комплексной амплитуды [4]:
ik
JJexp(-(u2 + v2)/о2 )х
I (x,y)= 2pf'
х exp(i j(u, v ))exp(- ik (xu + yv)l f )udv
(20)
B случае расчета ДОЭ, включающего в себя линзу, формула (20) перепишется следующим образом:
~( ) P Hexp(-(“2 + v2)/ о2 )expM«, v))х
I (x,y ) = 2PfJD .
х exp (ik (u2 + v2)/ 2f - ik (xu + yv)/ f )dv
(21)
Учитывая (14), можно переписать формулу (21) в следующем виде:
где
~х (х ) =
~(х,У )= ~х (х R (У ),
A
дjikf2pf Jexp(-u2/a2)x
(22)
- A
x exp(i jj (u))exp(iku2 /2f - г'Ахм / f )| , (23)
~y ) =
A
^ikj2pf Jexp(- v2/a2)x
- A
x exp(i j2 (v))exp(ikv2 /2 f - ikyv/ f )
(24)
B этом случае функционал
1 a i~
F(j(u,v))=— ^7х(х, jj (u)) - сх |^+
-a
+ Tb Jl~y (y,j2(v))-Cy|d-y
2b
(25)
-b
характеризует качество формирования заданного распределения интенсивности с помощью ДОЭ с фазовой функцией р(и, V).
Перейдя к дискретной модели с числом элементов дискретизации в плоскости ДОЭ и шагом дискретизации в фокальной плоскости, равным размеру дифракционного
пятна А / = А//(2А), можем заменить (25) функционалом
F(j)= ТТ X |Гх (х, ,j,(u) - С,
Nх n=0
N-,-1
+
+ ^Г X I7» (.j2(v )) - С
N y m = 0
Ny-1
y
, (26)
где
Nx = 2aIА f х, = -a + nA f ,
ym =-b + m^f , Ny = 2b / A f .
~х (х, )= X Cp exp(i(ji (up ^ , (27)
p=0
N-1
~y (ym )= X Cq exp( j2 (vq И , (28)
q=0
Up = -A + pA, vq = -A + qA .
up+i _________
C, = J J ik/ 2pf exp(- u^ a2 )x
up
x exp(- ilanul f + iku2 /2f )u
(29)
vp+i ________
Cqm = J J ik/ 2pf exp(- v21 a2)
vp
x ex p(- ikymvlf + ikv2 /2f
, (30)
А = 2 А / N. Отметим, что коэффициенты СП, Ст не зависят от выбора фазовых функций р1 ,р2 и могут быть посчитаны до начала работы оптимизационной процедуры, впрочем, также как и коэффициенты Сх, Су . Для минимизации функционала (26) в данной работе также как и в предыдущем пункте данной статьи использовалась модификация известной итерационной генетической процедуры [5]. Перед началом работы итераци-
онной процедуры формировалась начальная популяция, представлявшая собой набор
N„„„ сгенерированных случайным образом
gen
массивов
Р =01 р е {(2ят/М): те [1,М]. (31)
Каждый массив содержал последовательно записанные отсчеты функций р1 (и)
и Р 2 ():
[ Р1( А + 7а) 7 < N
[р2 (-А + ( - N)А) ) > N .
Pi =■
(32)
Для увеличения скорости сходимости в начальную популяцию вводилось в качестве начального приближения (т.е. одного из "индивидов" нулевой популяции) геометрооптическое решение р(и, V) задачи фокусировки гауссова пучка в прямоугольник, квантованное по М уровням. Такое решение представляет из себя суперпозицию двух одномерных фокусаторов в отрезки с длинами 2а и 2Ь:
p(U,v)=p1 (u)+^2 (v).
(33)
Для улучшения сходимости генетическая процедура совмещалась с применением метода вариаций на каждой итерации. После завершения итерационной процедуры, следующие критерии были выбраны для оценки качества фокусировки: энергетическая эффективность
N1 -1^-1~ ~ -lNy-1
е = X (хп )~у (Ут ) XX СхСу , (34)
п=0 т=0 I п=0 т=0
и среднеквадратичное отклонение сформированной интенсивности от среднего значения I
f Nx -lNy -1
ч1/2
s
EEEIx ([ fy (Vm )-i]
n=0 m=0
x
f Nx -1Nv-1r b
E ЕІ (x, Ж. kvm )f
-1/2
n=0 m=0
Nx -1Ny -1
, (35)
I =EE'x (x, % (vm )М^,,. (36)
Для постановки вычислительного эксперимента выбирались следующие параметры: длина волны освещающего пучка 1=10,6 мкм, радиус освещающего пучка 0=1,85 мм, фокусное расстояние элемента f=100 мм, размер апертуры 2*A=6,4 мм, размер фокальной области 2a*2b=1,75*3,5 мм, число отсчетов фазовой функции на апертуре элемента N*N=100*100, число уровней квантования M=8, количество итераций стохастической процедуры N.t=100, число индивидов исходной популяции Ngen = 20 , число индивидов, образованных операцией скрещивания N child = 20 . Уменьшение значения функционала невязки с ходом итераций показано на рис.4.
На рис. 5 представлены итерационно рассчитанная фаза ДОЭ (a) и интенсивность в фокальной плоскости элемента: (b) - результат численного моделирования итерационно рассчитанного ДОЭ, (c) - результат численного моделирования геометрооптического приближения.
Значения расчетной энергетической эффективности (34) и среднеквадратичного отклонения интенсивности (35), полученные в ходе вычислительного эксперимента для квантованного геометрооптического приближения и итерационно найденного решения, приведены в таблице.
Таким образом, применение стохастической итерационной процедуры с функционалом (25) позволяет рассчитывать восьмиуровневые ДОЭ, фокусирующие гауссов пу-
ф 1
1
ч s
V
1
1
нэмер итераи^и
n=0 m=0
Рис.4. Зависимость значения невязки от номера итерации
в а
а) б) в)
Рис.5. Итерационно рассчитанная фаза ДОЭ (черный цвет соответствует значению фазы 2л, белый - 0) (а), (Ь) - распределение интенсивности, сформированное итерационно рассчитанным ДОЭ в фокальной плоскости (результат численного моделирования), (с) - распределение интенсивности, сформированное геометрооптическим приближением
чок в прямоугольник с энергетической эффективностью, близкой к эффективности геометрооптического решения, и более чем в 3 раза меньшим значением среднеквадратичного отклонения для выбранного набора физических параметров. Отметим также небольшие вычислительные затраты, требуемые для реализации расчета элемента, что объясняется простотой вычисления функционала (26) после того, как значения коэффициентов (29) и (30) предварительно посчитаны.
Заключение
В данной работе рассмотрено применение методов стохастической оптимизации к расчету квантованных двумерных ДОЭ, фокусирующих лазерное излучение в двумерные области. Приведенные результаты вычислительных экспериментов подтверждают перспективность применения стохастического подхода к оптимизации квантованных ДОЭ.
Таблица. Результаты вычислительного эксперимента
Метод расчета e, % d,%
Г еометрооптический 87 45
Итерационный 88 13
СПИСОК ЛИТЕРАТУРЫ
1. Сойфер В.А. Глава 1. Введение в компьютерную оптику // Методы компьютерной оптики / Под ред. В.А. Сойфера. М: Физ-матлит, 2000.
2. Головашкин Д.Л., Волков А.В., Казанский Н.Л., Павелъев В.С., Сойфер В.А., Усп-ленъев Г.В. Глава 4. Технология создания ДОЭ // Методы компьютерной оптики / Под ред. В.А. Сойфера. М.: Физ-матлит, 2000.
3. Fienup J.R. Iterative method applied to image reconstruction and to computer-generated holograms // Optical Engineering. 1980. V.19.
4. Soifer V.A., Kotlyar V.V., Doskolovich LL. Iterative Methods for Diffractive Optical Elements Computation. London: Taylor & Francis Ltd. 1997.
5. Lu C. Y., Liao H.Z., Lee C.K., Wang J.S. Energy control by linking individual patterns to self-repeating diffractive optical elements // Applied Optics. 1997. V.36. №20.
6. Duparre’ M., Golub M.A., Luedge B., Pavelyev V.S., Soifer V.A., Uspleniev G.V., Volotovskii S. G. Investigation of computergenerated diffractive beam shapers for flattening of single-modal CO2-laser beams // Applied Optics. 1995. V.34. №14.
STOCHASTIC APPROACH TO QUANTIZED DIFFRACTIVE OPTICAL ELEMENTS OPTIMIZATION
© 2002 V.S. Pavelyev
Image Processing Systems Institute of Russian Academy of Sciences, Samara
This paper is devoted to using the stochastic optimization procedure for designing quantized DOEs that focus light into pregiven focal domains.
