Исследование дифракционных характеристик фокусатора в кольцо методом вычислительного эксперимента

Положенные в основу расчета фокусаторов в кольцо [1,3] геометрооптические представления не позволяют оценить перечисленные выше дифракционные параметры. При этом для фокусатора в кольцо важно выяснить, какие дифракционные эффекты обусловлены несовершенством геометрооптической фазовой функции, конечностью размера фокусатора, а какие - технологией синтеза элементов компьютерной оптики (ЭКО)[4], т. е. дискретизацией и квантованием фазовой функции. В работе [5] исследована объемная структура светового поля вблизи фокальной плоскости, получены аналитические оценки ширины фокального кольца и энергетической эффективности фокусировки без учета дискретизации ’ и квантования фазовой функции.   В работе   [6]

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

ФОКУСАТОР В КОЛЬЦО

На основе решения обратной задачи фокусировки лазерного излучения в работах [1,3] предложено несколько «видов» фазовой функции фокусатора в кольцо. Используется для фокусировки плоского пучка света в кольцо и пара «аксикон-линза» [2,8,9]. Фазовая функция плоского оптического элемента, действующего аналогично паре «аксикон-линза», имеет вид:

*’(.) = ^о " k^2 + |г-г/ , О s г S а;                (11

здесь г - полярный радиус; к = 2тг/А, А - длина волны света; f - фокусное расстояние; а - радиус фокусатора; ро = const; rQ = const.

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

.2 krr —-^

^г) - % - 4 - -1 * — > Ro = /f+ro (2)

видно, что рассматриваемый элемент эквивалентен комбинации собирающей линзы с фокусом R_q (квадратичная по г фаза) и аксикона [7] (линейная по г фаза). Как следует из работы [1], оптический элемент (1) обладает оптимальной фазовой функцией с точки зрения энергетической эффективности фокусировки в кольцо, поэтому остановимся на исследовании именно этого элемента.

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

— ( 2Rt

.2 2П # о> = М“п ’

откуда

2XR п о

“М

Здесь для каждого значения п - 1, 2, ..., N

возможны одно или два

значения радиуса кольца:

, О s г s а; п n = 1, 2, . . . , N,

а общее количество кольцевых зон

N =

г г2М

(а - r0>2 м-

2AR о

+ 1,

где [^] означает целую часть числа £.

ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ

При исследовании ограничений геометрооптического расчета фокусаторов и влияния эффектов дискретизации на функционирование ЭКО необходимо использование как аналитических методов расчета [5], так и средств вычислительного эксперимента [6,10]. Для решения указанных задач необходимо рассчитать и исследовать создаваемое фокусатором световое поле, определить значения важнейших характеристик фокусаторов - ширину фокальной линии и энергетическую эффективность фокусировки. Такие факторы, как большой объем расчетных данных, многовариантность значений параметров исследуемого светового поля, наличие набора критериев качества, требования естественности восприятия результатов расчета и возможности диалога исследователя с ЭВМ, выявление неизвестных свойств исследуемого объекта, выводят решаемую задачу из рамок расчетных, позволяя охарактеризовать ее как вычислительный эксперимент [6, 10].

Для решения сформулированных задач предлагается использовать агрегированную математическую модель фокусировки лазерного излучения элементами компьютерной оптики, основанную на представлении ЭКО в виде совокупности элементарных модулей кольцевой, линейной или прямоугольной формы и представлении поля от фокусатора в виде суммы базовых решений -аналитически рассчитываемых полей, создаваемых отдельными модулями.

В рамках этой модели считается, что на поверхности подложки ЭКО имеются координаты it = (u, V). Область G, занимаемая ЭКО, разбивается на N неперосекающихся ячеек (модулей) Gn, нумеруемых индексами nel :

LT G = G, G n G ,= 0 (п * п').                 (7)

п                     п п п€1 N

Каждая ячейка соответствует одному положению сканирующего устройства. Значение t^ функции амплитудно-фазового пропускания синтезируемого ЭКО в пределах ячейки полагается постоянным и называется отсчетом этой функции.

Полученная функция комплексного пропускания синтезированного ЭКО является кусочно-постоянной и описывается выражением

Т (3) = Е t • X (5), 3eG, n n

где X(u) - -I1, U®G„

" к 55g .                                            (9)

n

При переходе от ячейки к ячейке значение отсчета изменяется в соответствии с сигналом, поступающим от ЭВМ. В памяти ЭВМ n-му отсчету соответствует последовательность из ш двоичных разрядов.

Введение интерполирующей функции (9) соответствует «равномерной засветке» в пределах одного элемента разрешения амплитудной маски, что имеет место в известных устройствах [11 - 13] . Наибольшее распространение при синтезе фокусаторов получили два типа устройств. При изготовлении ЭКО на устройстве с построчной разверткой [12, 14] фокусатор представляется в виде системы прямоугольников, а на устройстве с круговой разверткой [11, 13] - в виде системы концентрических колец с различными значениями коэффициента амплитудно-фазового пропускания.

При этом поле в фокальной области фокусатора определяется как сумма полей, созданных его модулями:

w( _х. _у, _z) _ 11^ Е t„ J E(U.V)- ехр( Iks)_dudv.         _(i0)

* "е,н с п где E(u, v) - падающий на фокусатор пучок света длины волны X;

f      - фокусное расстояние;

(х, у) - декартовы координаты в фокальной плоскости; z      - расстояние до фокальной плоскости;

к = 2п/Х;

s = /Tf+zT^TTu^xT5^^ .                (и)

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

Термин «агрегировать», согласно словарю иностранных слов (И. Русский язык, 1983, с. 15), означает «объединять, суммировать какие-либо однородные показатели (величины) с целью получения более общих, обобщенных показателей (величин)». Поскольку в рамках разрабатываемой модели фоку- сатор представляется совокупностью (объединением) модулей, а поле от фоку-сатора описывается суммой базовых решений задачи дифракции света на отдельных модулях, и это описание используется для получения дифракционных характеристик всего фокусатора, употребление термина «агрегированная» представляется обоснованным.

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

На начальном этапе исследований важно провести аналитический дифракционный расчет структуры сфокусированного излучения с учетом конечных размеров фокусатора [1,5,15]. Анализ полученных дифракционных соотношений позволяет исследовать ограничения положенного в основу расчета фазовой функции фокусатора геометрооптического подхода, выявив начальные значения физических параметров, при которых происходит разрушение требуемой формы области фокусировки.

На первом этапе вычислительного эксперимента необходимо показать принципиальную работоспособность исследуемого геометрооптического фокусатора с учетом дифракционных эффектов на элементах разрешения ЭКО. Такое исследование целесообразно проводить для реальных физических параметров и при максимально возможном числе М уровней квантования фазовой функции фокусатора (МЫ6).

Для понимания направленности дальнейших исследований следует подчеркнуть, что наиболее качественный фазовый рельеф получается с помощью фотолитографических методов травления [11,16]. Однако из-за высочайших требований к точности совмещения бинарных фотошаблонов, используемых для формирования фазового рельефа по данной технологии, трудно получить большое число градаций фазовой функции. Поэтому следующий важный этап вычислительного эксперимента - это исследование структуры сфокусированного излучения для малого числа уровней квантования фазы фокусатора (М = 2; 4).

Результаты проведенных исследований важно представить в виде, привычном для оптика-экспериментатора: полутоновые распределения яркости на экране телевизионной трубки, графики кривых, изофоты (линии равной интенсивности [17]), изометрии пространственных распределений интенсивности. В результате автоматизированного анализа структуры сфокусированного излучения исследователь должен получить оценки дифракционных параметров фокусатора.

МЕТОДЫ ВЫЧИСЛЕНИЯ ДИФРАКЦИОННОГО ИНТЕГРАЛА

Возможности вычисления на современных ЭВМ интеграла (10) с быстроосциллирующим (ввиду малости А) подынтегральным выражением представляются весьма ограниченными. Имеется возможность существенно уменьшить вычислительные трудности, используя следующее разложение в ряд:

S s s + n

2 C

n

2s n

где s определяется формулой (11);

s = / ( f+Z) 2+(X-U )2+(y-V )2; n                         n           n

( 13)

c = ( u -u) ( 2X-U-U ) + n        n                n

(v-v_r) ( гу-v-vj,

( 14)

где (u,v ) - «центральная* точка модуля G. n n n

Оставляя только первые два члена разложения (12) при аппроксимации фазовой составляющей подынтегрального выражения (10) и только первый член разложения ( 12) при аппроксимации амплитудной составляющей, приходим к приближению, которое называется приближением Френеля [17,18]. Отличие предлагаемого подхода от известных [17,18] заключается в применении приближения Френеля не ко всему оптическому элементу, а к каждому отдельному (относительно малому) модулю ЭКО, что позволяет существенно повысить точность расчета.

Фазовая составляющая подынтегрального выражения (10) более чувствительна к аппроксимации, чем амплитудная. Так, использование двух первых членов при аппроксимации фазовой составляющей ( 10) является более строгим ограничением, чем использование одного члена при аппроксимации амплитудной составляющей [19].

w( х, у, z) = -

• Je( U, V)

G п

• -    „ t ехр(iks )

Г + Z v n n . .2

iA n€is

Nn

• / ike x

expj—2s~~ [ du dv-

( 15)

Исследованию точности приближения Френеля посвящены работы [19, 20]. В них сравнением результатов прямого численного интегрирования по формулам (линейный случай)

w (х, z)

А

E(u) ехр(iks)

UK du s

s

и

w2 (X,z)

xT^z) exP|ik(f+z)

Е(и)ехр

ik( x-u) ² 2( f+z)

( 17)

где

s

(f+z)2+ (x-u)2

( 18)

следующим образом: N N 1     2

т (u.v) =2 X Р=1  1=1

где |u| s U, |v| s V, (Up, vj- центр (р, 1)-го кания t ;

pi rect(х)

u = - и р

V = - V

При использовании для изготовления ЭКО с вращательной симметрией фотопостроителя с круговым сканированием и при аксиальной симметрии освещающего пучка E(u, v)=E(r) для вычисления интеграла (15) применяются методы, описанные в работах [6, 21].

При изготовлении ЭКО на фотопостроителе с построчной разверткой функция комплексного пропускания ЭКО в соответствии с (8) будет выглядеть

, U~U х          , v-vх pi reCE (  SU J reCt ( Sv  J’(19)

модуля ЭКО с функцией  комплексного  пропус-

= $1‘~^sxs^(

1 О, |х| > |;

+ (р - j-) • 6U, р = ТИТ ;

+ (1 - Ь • Sv, 1 =Т1Г ; 22

” ; 4у - -SY. .

N=

Поле w( х, у, z)b фокальной области фокусатора с прямоугольным растром вычисляется следующим образом:

(«У Х)2+ (V,- У)2-

(х-u) 2+ (у-у)2]

2s р 1

du dv,

где s^ = (f+Z)2+ (X-Up)2+ (y-vj2, u j = ~ U + (p-1) Su, p-F

U ! = - U + p . Su,             p

p*:t

= - V + ( 1-1) • 5v,

v = - V + 1- 6-v,

1 1+-

1 = i,n₂.

Для плоского или сферического

освещающих пучков в приближении Френеля

Е( u, v) а А( u, v) • ехр •{ i ^f (u) + ^ (v) I 1                 2

и интеграл в (21) факторизуется, т. е.

w( х, у, z) = -

„   *1   *2 А( u , V ) -t , r+ z V V _____Е_____!____р 1

iX p=i 1=1        s2

exp-l ik s

I I р 1

х)2+

2S

(v^y)2

Следует отметить, что предлагаемый подход (25) точнее, чем методы дифракционного расчета, использующие быстрое преобразование Фурье (БПФ) [22 - 25] за счет более точного вычисления знаменателей фазовой (s вместо f + z) и амплитудной - вместо -^^ составляющих. Однако ос-v s                      '

р 1

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

При выполнении условий приближения Френеля для всего фокусатора (ks»i и f/2a £12) можно существенно ускорить вычисления:

exp ik(f+z)

^w<^x-y.z) = - _ix^L (f^)"¹ •

^H1 ^upx —                                                       (26)

Zr 2 |'Г> X M X-U) 21\ ■

J expji ^ ₁ ( u) + _{2( f+ir}J J du p = 1 u 1

Согласно (26) для вычисления поля в точке (х,у, z)

В1>1Числить h^ интегралов

Vх* - J вкр^^/ч) + гПТ^Г"]/ du u 1

и N2 интегралов v 1 1 ♦ -----

^В1^(у> * Г exp|ib (V) + 2(f+z) I ^dv"                              ⁽²⁸⁾

V 1 1-- 2

Таким образом, для вычисления поля в одной точке требуется найти Nt+ N2 интегралов типа интеграла Френеля, а не N^^, хотя полной факторизации сумм добиться невозможно, так как коэффициент амплитудно-фазового пропускания прямоугольных модулей t свойствами факторизации обладать не обязан. В случае, когда необходимо найти матрицу отсчетов поля ■lw(x , у , z ), к = 1, К ; m = 1, М I в сечении фокальной области z=z , доста-I к ш О                                     J                                                          и точно вычислить не K-M-N^N интегралов, а только K-Ni интегралов А?( xj согласно (27) и M'N2 интегралов В^у^) по формуле (28).

Существует целый ряд методов [18,26-29] вычисления интегралов ( 27), (28) для плоского и сходящегося сферического (при факторизации - цилиндрического) освещающих пучков.

Наиболее простой способ заключается в сведении интеграла в (27) к разности двух значений интегралов Френеля [30], программа вычисления которых имеется в Библиотеке стандартных программ [31]. Данный способ приемлем в случае малого числа элементов дискретизации N *N • Для больших N_t при вычислении интеграла (27) целесообразно применение обычных квадратурных методов [18], но наиболее эффективно использование метода локальной линейной аппроксимации фазовой функции подынтегрального выражения [28, 29]. Последний метод используется и при вычислении двумерного интеграла в (21), когда невозможна факторизация освещающего пучка Е(и, у) или отсутствуют условия для применения приближения Френеля.

ПРОГРАММНЫЙ КОМПЛЕКС ДЛЯ ВЫЧИСЛИТЕЛЬНОГО ЭКСПЕРИМЕНТА

Программный комплекс для вычислительного эксперимента создан на базе пакета прикладных программ обработки изображений и цифровой голографии (ППП) [32], комплекса графических программ ГРАФОР [33] и подробно описан б статье [34]. Используемые в качестве входных параметров массивы отсчетов фазовых функций фокусаторов рассчитываются и кодируются средствами комплекса программ синтеза элементов компьютерной оптики [35].

Вычислительный эксперимент производится с терминала ЭВМ ЕС-1061. оснащенной графопостроителем ЕС-7054, полутоновым дисплеем, графическим дисплеем, и основан на совместном применении ППП, средств визуализации световых полей, средств машинной графики и организации диалога с исследователем. Результаты в виде распределения интенсивности отображаются на экране полутонового дисплея. На графическом дисплее отображаются семейства графиков, изофоты или условные трехмерные изображения световых полей. Для получения твердой копии результатов эксперимента используются графопостроитель и многоградационный фотопостроитель.

РЕЗУЛЬТАТЫ ВЫЧИСЛИТЕЛЬНОГО ЭКСПЕРИМЕНТА

При исследовании квантованного фокусатора (1) с дискретизацией (5) (6) были выбраны следующие параметры: фокусное расстояние f = 500 мм; длина волны А = 10,6 мкм; диаметр 2а = 25,6 мм. В табл. 1, на рис. 2, 3 представлены результаты вычислительного эксперимента с фокусаторами в кольцо радиусом ro = 1 мм (рис. 2) и rQ = 2 мм (рис.З) для различного числа уровней квантования фазы (М = 2; 4; 16).

В табл. 1 используются следующие значения:

• 9    - уровень интенсивности (относительно максимальной), по которому определяются границы кольца фокусировки ( р-, р+);

• с    - энергетическая эффективность фокусатора (доля энергии, попадающая в кольцо фокусировки);

• ( Р_» Р ₊ ) - границы кольца фокусировки, определяемые по уровню 9 максимальной интенсивности;

Др = р*- р_ - ширина кольца фокусировки по уровню 9;

&Рд~ теоретическая оценка ширины кольца фокусировки для неквантованного фокусатора (М = ®, N = «), согласно [5];

• е _@ - теоретическая оценка энергетической эффективности неквантованного фокусатора (табл. 2 статьи [5]).

Анализ приведенных результатов вычислительного эксперимента позволяет сделать следующие выводы:

• -    с уменьшением числа уровней квантования происходит расплывание кольца фокусировки (увеличение ширины кольца и уменьшение доли энергии, попадающей в него);

• ~    при этом ширина кольца и энергетическая эффективность практически ^Не зависят от радиуса кольца г ;

• -    относительная яркость центрального всплеска уменьшается с ростом Радиуса кольца фокусировки г ;

• -    экспериментальный радиус центра кольца фокусировки увеличивается с Уменьшением числа уровней квантования;

• -    теоретическая оценка дифракционной ширины кольца фокусировки Др^ Наквантованного фокусатора неплохо согласуется с данными вычислительного эксперимента;

Зависимость энергетической эффек фокусировки Др по уровню е и от числа уровней М квантования фазы в кольцо радиусом r_Q ( f = 500 мм;

Го (мм)	М	е - 0,5    (Др@ - 0,408 мм)					е - о, 2 ^Ре ■
		е	Р_ (мм)	(мм)	Др (мм)	Су 4е	с	р_ ( мм)	(мм)
1, о	16 4 2	0,682 0,798 1,179 0,381 0,958 0,845 0,652 1,259 0,570 0,811 1,193 0,382 0,801 0,699 0,650 1,266 0,289 0,829 1,218 0,389 0,406 0,353 0,632 1,289
2, 0	16 4 2	0,686 1,797 2,189 0,392 0,963 0,844 1,684 2,282 0,578 1,803 2,220 0,417 0,812 0,695 1,687 2,305 0,282 1,863 2,242 0,379 0,396 0,345 1,729 2,330

Рис. 2. Графики распределения нормированной интенсивности i = IM(^ ,0)/I16(rQ, 0) в фокальной плоскости фюкусатора плоского пучка в кольцо с параметрами:

f - 500 мм; Л = 10,6 мкм; 2а - 25,6 мм; г0 - 1 мм для различного числа М уровней квантования фазы: а - М =16;

б - М = 4; в - М = 2

• -    относительная эффективность квантованного фокусатора £/£_е соответствует относительной эффективности 1_м( О, О)/К о, О) плоской линзы с тем же числом градаций фазы [11] (здесь 1_и(0,0) - интенсивность в фокусе плоской линзы с уровнями квантования фазы, а 1(0,0)  - интенсивность в

фокусе безаберрационной линзы с теми же физическими параметрами).

Таблица 1

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

А = 10, 6 мкм; а = 12, 8 мм)

о, 6i:	) мм)	е ■	0, 1	(Лре -	0, 701	мм)
Ар	с,		Р_	р+	Ар	Су
(мм)			( мм)	(мм)	(мм)	'се
0, 608	0, 992	0, 873	0, 617	1, 306	0, 690	0, 989
0, 617	0, 820	0, 721	0, 611	1, 310	0, 700	0, 817
0, 657	0, 414	0, 365	0, 574	1, 338	0, 764	0, 413
0, 588	0, 991	0, 872	0, 643	2, 328	0, 685	0, 988
0, 618	0, 816	0, 717	1, 644	2, 348	0, 704	0, 812
0, 601	0, 405	0, 364	1, 615	2, 382	0, 767	0, 412

Рис. 3. Г рафики распределения нормированной интенсивности i = 1М( $ ,О)/116(го, 0) в Фокальной плоскости фокусатора плоского пучка в кольцо с параметрами:

мм, X ■ 10,6 мкм; 2а - 25,6 мм; г₀ ■ 2 мм для различного числя М уровней квантования фазы: а - М - 16; о-М=4;в-М«2

Таким образом, эффективность квантованного фокусатора с неравномерной дискретизацией можно определить, зная долю энергии, идущей в рабочий дифракционный порядок при данном числе уровней квантования фазы [11], и эффективность с@ неквантованного фокусатора, описываемого своей геометрооптической формулой (табл. 2 статьи [5]):

е = с0 • 1н (0,0)/1 (0,0).                         (29)

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

При фокусировке гауссова пучка с амплитудным распределением ( г2 1

A(r) = Ао ехр^- — j                       (30)

можно воспользоваться тем же самым фокусатором (1),  (5)-(6). Результаты вычислительного эксперимента с фокусатором в кольцо для гауссова освещающего пучка представлены на рис.4. Параметры фокусатора и освещающего пучка: f = 500 ММ; а = 12, 8 ММ; Л = 10, 6 МКМ; (Г = 0,4 а; Г = 1 ММ;

о

М = 2, 4, 16. Сравнение рисунков 2 и 4 показывает, что при фокусировке гауссова пучка увеличивается ширина кольца фокусировки, что обусловлено фактическим уменьшением активной области фокусатора при выбранном значении а. При большом числе уровней квантования происходит подавление центрального всплеска, который сливается с окружающим фоном. Этот результат совпадает с данными теоретического исследования неквантованного фокусатора [9].

При отсутствии кольцевого фотопостроителя можно изготовить фокусатор в кольцо на фотопостроителе с прямоугольным растром типа PHOTOMATION Р-1700. Дифракционный расчет светового поля от оптического элемента такого типа основан на вычислении интеграла (21).

Результаты исследования элемента компьютерной оптики (плоского аксикона) с прямоугольным растром, фокусирующего в кольцо сходящийся

Рис. 4. Графики распределения нормированной интенсивности i = ^(^ , 0)/I₁₆(fq, 0) в фокальной плоскости фюкусатора в кольцо для различного числа уровней квантования при фокусировке гауссова пучка. Параметры:

Г = 500 мм; X ■ 10,6 мкм; 2а - 25,6 мм; б - 0,4 а; г₀ - 1 мм; а - М - 16; б - М - 4; в - М -2

сферический пучок света, представлены на рис. 5 и в табл. 2. Параметры фокусатора: f = 750 мм; X = 0,6328 мкм; 2а = 25,6 мм; r = 1 мм; N-N - число о квадратных элементов дискретизации, но используются из них только те, центры которых попадают в круг радиусом а; М - число уровней градации фазы. Физические параметры фокусатора выбраны равными физическим парамет--рам плоского аксикона с равномерной дискретизацией по радиусу, который рассматривался в статье [6]. Ширина кольца фокусировки Ар=б4 мкм выбрана в соответствии с оценками работ [5, 6]. Сравнение рис. 2 и табл. 2

Таблица 2

Зависимость энергетической эффективности с фокусатора сходящегося сферического пучка в кольцо с прямоугольным растром дискретизации от числа уровней М квантования фазы

N 00 16 4 2 128 0, 472 0, 470 0, 378 0, 327 256 0, 758 0, 741 0, 602 0, 584 статьи [6] с рис.5 и табл.2 настоящей статьи свидетельствует о значительной схожести результатов. Но для энергетической эквивалентности необходимо, чтобы число элементов дискретизации N вдоль каждой оси

Рис. 5. Графики распределения нормированной интенсивности i = 1М(У, 0)/!^ (гд, 0) в фокальной плоскости фокусатора в кольцо с прямоугольной дискретизацией растра. Параметры:

f * 500 мм; Л - 0,6328 мкм; 2а - 25,6 мм; r0 = 1 мм; N, = N2 » 128; а - М - 16; б - М - 4; в - М = 2

прямоугольного растра было в 2 раза большим, чем количество колец при равномерной дискретизации вдоль радиуса. Об этом же свидетельствует большее количество и яркость паразитных колец на рис.5 по сравнению с рис. 2 статьи [6]. В то же время проведенный анализ показывает, что для исследования матричных фокусаторов с радиально-симметричной фазовой функцией может быть использован математический аппарат, предложенный в работах [6, 21], вместо методов (21) -(28), применяемых для анализа ЭКО с прямоугольным растром, с соответствующей заменой N2 квадратных модулей на N/2 кольцевых, что позволит существенно сократить время вычислений.

ЗАКЛЮЧЕНИЕ

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