close

Вход

Забыли?

вход по аккаунту

?

Расчет хода псевдолучей через дифракционные структуры выполненные на сферической поверхности.

код для вставкиСкачать
РАСЧЕТ ХОДА ПСЕВДОЛУЧЕЙ ЧЕРЕЗ ДИФРАКЦИОННЫЕ СТРУКТУРЫ,
ВЫПОЛНЕННЫЕ НА СФЕРИЧЕСКОЙ ПОВЕРХНОСТИ
Е. Г. Ежов1,С. А. Степанов2
Институт систем обработки изображений РАН
2
Пензенская государственная архитектурно-строительная академия
1
Представлен математический аппарат расчета хода лучей в приближении заданного порядка малости через
вращательно-симметричные оптические системы, включающие дифракционные линзы, структура которых размещена на сферической поверхности.
1. Введение
В работах [1-3] описана методика и приведены
формулы расчета хода псевдолучей лучей (т.е. лучей, траектории которых рассчитываются в приближении заданного порядка малости) через вращательно-симметричные системы, включающие радиально-градиентные и дифракционные линзы. В этих
работах предполагалось, что структура дифракционной линзы (ДЛ) выполнена на плоскости. Однако
известно, что размещение структуры ДЛ на сферической поверхности открывает возможность управление комой линзы (вплоть до ее устранения) независимо от величины сферической аберрации, расположения предмета и зрачка [3]. Очевидно, что использование ДЛ на сферической поверхности в составе гибридной оптической системы может существенно расширить ее коррекционные возможности,
поэтому целесообразно распространить методику
расчета хода псевдолучей на такие элементы.
Для получения формул расчета хода псевдолуча необходимо иметь представленные в аналитическом виде точные соотношения, связывающие параметры дифрагировавшего и падающего на ДЛ лучей. Такие соотношения для элементов, структура
которых размещена на поверхностях произвольной
формы приведены в работах [4, 5]. Точные лучевые
соотношения следует разложить в ряды по степеням
входных параметров луча. В результате после разложения и соответствующей перегруппировки параметры луча на выходе из ДЛ будут описываться
степенными рядами, состоящими из слагаемых различных порядков малости относительно параметров
луча на входе в оптическую систему. В результате
ограничения рядов конечным числом слагаемых будут получены выходные параметры псевдолуча, вычисленные в требуемом приближении. С целью упрощения вышеописанной процедуры точные лучевые соотношения целесообразно представить в виде
наиболее удобном для последующего разложения.
2. Точные лучевые соотношения
Требуемые лучевые соотношения получим для
представляющего наибольший практический интерес случая, когда размещенная на сферической поверхности кольцевая структура ДЛ симметрична от-
носительно оси, проходящей через центр сферы и
являющейся оптической осью линзы. Введем две
прямоугольные системы координат. Одна из них −
K −система с началом координат, совмещенным с
вершиной сферической поверхности, и осью Z ,
совпадающей с оптической осью. Вторая − локаль~
ная K −система, ее началом координат является
точка падения луча на сферическую поверхность, а
~
ось Z совпадает с нормалью к поверхности в этой
точке. В K −системе уравнение сферической поверхности, на которой размещена структура ДЛ,
имеет вид [6]
z = (1 c )⎡1 − 1 − (cρ )2 ⎤ ,
⎢⎣
⎥⎦
где
c
(1)
− кривизна сферической поверхности,
ρ = x2 + y2 .
В этой же системе координат направляющие
~
~
косинусы оси Z (составляющие орта k ) описываются выражениями
~
~
~
k x = −cx, k y = −cy, k z = 1 − cz .
(2)
~
~
Что же касается двух других ортов i и j , то
их составляющие легко найти, воспользовавшись
очевидными соотношениями
~ ~
~ ~
~ ~
i ⋅ j = 0, i ⋅ k = 0, j ⋅ k = 0⎫⎪
⎬.
~2
~
⎪⎭
i = 1, j 2 = 1
(3)
Поскольку (3) представляет собой систему пяти
скалярных уравнений с шестью неизвестными, то
~
~
одна из проекций орта i или j может быть выбрана произвольно. Положим, например,
~
jx = 0 ,
(4)
~
~
тогда другие составляющие ортов j и i будут определяться соотношениями
~
~
ky
kz
~
~
(5)
jy = ~
~ , jz = − ~2 ~ 2 ;
k y + kz
k y2 + k z2
25
~
ix =
~
iy =
~
iz = −
(k~ ~j
~~ ~ ~
k z j y − k y jz
)
~ ~ 2 ~2
+ kx
z y − k y jz
~~
k x jz
~~ ~~ 2 ~
k z j y − k y j z + k x2
~~
kx jy
(
)
(k~ ~j
z y
~~
− k y jz
)
2
(~j
2
y
(~j
2
y
(
~
+ j z2
~
+ j z2
)
)
~ ~
~
+ k x2 j y2 + j z2
)
⎫
⎪
⎪
⎪
⎪
⎪⎪
⎬.
⎪
⎪
⎪
⎪
⎪
⎪⎭
α1x = −
α1 y = −
(6)
y
s (ρ s ) + (1 − z s )2
2
α 2x = −
⎫
⎪
⎪
⎪
⎪
⎬,
⎪
⎪
⎪
⎪⎭
x
s ′ (ρ s ′) + (1 − z s ′)
2
(7)
(8)
~′ = α
~ −α
~ +α
~ , α
~′ = α
~ −α
~ +α ⎫
α
1x
2x
1y
2y
x
x
y
y
⎪
⎬ , (9)
2
2
~
~
~
α′z = 1 − α′x − α′y
⎪⎭
Направляющие косинусы, входящие в правую
часть выражений (9), изначально известны только в
K −системе. В частности, направляющие косинусы
нормали к первому (сферическому) фронту записи
определяются выражениями
2
1
3
b3 ρ 2 + b5 ρ 4 − ...
2
8
y
α 2y = −
2
s ′ (ρ s ′) + (1 − z s ′)2
+
где m − порядок дифракции, µ = λ λ 0 , λ - длина
волны света, падающего на ДЛ, λ0 – длина волны
записи (в дальнейшем будем полагать, что m = −1 и
µ = 1 ).
Производные же эйконала записи можно выразить через направляющие косинусы нормалей к
двум фронтам записи, по крайней мере, один из которых является асферическим. В результате соотношения (8) приобретают вид
26
2
(10)
а направляющие косинусы ко второму (асферическому) фронту записи – соотношениями
с помощью которой можно перейти от составляющих любого вектора, определенных в K −системе к
составляющим этого же вектора, определенным в
~
K −системе. Обратный переход производится с помощью транспонированной матрицы Α′ .
Целесообразность использования локальной
~
K −системы обусловлена тем, что именно в этой
системе справедливы простые соотношения, связывающие направляющие косинусы дифрагировавшего и падающего лучей с производными эйконала записи структуры ДЛ:
~′ = α
~ + mµ ∂G0 , α
~′ = α
~ + mµ ∂G0 ⎫
α
x
x
y
y
~
∂x
y ⎪⎪ ,
∂~
⎬
⎪
2
2
~
~
~
α′z = 1 − α′x − α′y
⎪⎭
s (ρ s ) + (1 − z s )
2
α1z = 1 − α12x − α12y
~ ~
~
Составляющие ортов i , j и k образуют
матрицу
~ ~ ⎞
⎛~
⎜ ix iy iz ⎟
~ ~ ~
Α = ⎜ jx j y jz ⎟ ,
⎜
⎟
⎜ k~ k~ k~ ⎟
⎝ x y z⎠
x
+
1
3
b3 ρ 2 + b5 ρ 4 − ...
2
8
α 2 z = 1 − α 22 x − α 22 y
⎫
+⎪
⎪
⎪
⎪
⎪
⎪⎪
+⎬ ,
⎪
⎪
⎪
⎪
⎪
⎪
⎪⎭
(11)
где − s и s ′ – отрезки записи, b3 , b5 , …− коэффициенты асферический деформации третьего, пятого
и т.д. порядков [5].
Таким образом, чтобы найти направляющие
косинусы дифрагировавшего луча в K −системе необходимо вначале направляющие косинусы падающего луча и нормалей к фронтам записи, известные
в этой системе, пересчитать с помощью матрицы Α
~
[см. выражение (7)] в K −систему. Затем, используя
соотношения (9), получить направляющие косинусы
~
дифрагировавшего луча в K −системе и совершить
с помощью транспонированной матрицы Α′ обратный переход в K −систему.
3. Схема расчета хода псевдолуча
Рассмотрим теперь, как можно использовать
вышеприведенные соотношения для расчета хода
псевдолуча через структуру ДЛ, размещенную на
сферической поверхности. Для определенности ограничимся приближением пятого порядка малости.
Высоту и наклон псевдолуча в точке его падения на поверхность ДЛ опишем векторами S и T .
Вектор S имеет составляющие (x, y, 0 ) , вектор T −
(
)
составляющие Tx , T y , 0 , где Tx и T y − направляющие тангенсы луча связанные с направляющими косинусами луча соотношениями
Tx = α x α z ,
Ty = α y α z .
(12)
Векторы S и T являются функциями двух
зрачковых векторов, определяющих высоту и наклон луча во входном зрачке оптической системы, и
могут быть представлены в виде суммы слагаемых
первого, третьего, пятого и т. д. порядков малости
относительно модулей этих зрачковых векторов:
S = S (1) + S (3) + S (5) + ... ⎫⎪
⎬.
T = T(1) + T(3) + T(5) + ...⎪⎭
Введем три инварианта вращения
(14)
( ) + (T )
= 2(T T + T
2
e1( 4)
(1) 2
y
(1) (3)
x
x
(1) (3)
y Ty
( ) + (y )
= 2(x x + y
e3( 4)
2
(1)
(1) 2
(3)
(1)
y (3)
(16)
)
e 2( 4) = x (1) T x( 3) + x (3) T x(1) + y (1) T y(3)
e3( 2) = x (1)
( )
⎫
⎪
⎬ , (17)
(3) (1)
+ y T y ⎪⎭
⎫
⎪
⎬,
⎪⎭
(18)
)
1)
3)
5)
α px = α (px
+ α (px
+ α (px
+ ...⎫
⎪
⎪
1)
3)
5)
α py = α (py
+ α (py
+ α (py
+ ...⎬
⎪
α pz = 1 + α (pz2) + α (pz4) + ... ⎪
⎭
z = z ( 2) + z ( 4) + ... ,
(19)
где
1 ( 2)
ce1 ,
2
z ( 4) =
~ ~
~
~
k x = k x(1) + k x(3) + k x(5) + ...⎫
⎪⎪
~ ~
~
~
k y = k y(1) + k y(3) + k y(5) + ...⎬ ,
⎪
~
~
~
k z = 1 + k z( 2) + k z( 4) + ... ⎭⎪
~
~
~
k y(1) = −cy (1) , k y(3) = −cy (3) , k y(5) = −cy (5) ; (28)
Далее, поскольку направляющие косинусы падающего на ДЛ луча связаны с его направляющими
тангенсами соотношениями
(21)
то после разложения в степенные ряды получим
α x = α (x1) + α (x3) + α (x5) + ...⎫
⎪⎪
α y = α (y1) + α (y3) + α (y5) + ...⎬ ,
⎪
α z = 1 + α (z2) + α (z4) + ... ⎪⎭
где
(26)
~
~
~
k x(1) = −cx (1) , k x(3) = −cx (3) , k x(5) = −cx (5) ; (27)
( )
1 + T x2 + T y2 ⎫
⎪
2
2⎪
α y = T y 1 + Tx + T y ⎬ ,
⎪
α z = 1 1 + T x2 + T y2 ⎪
⎭
(25)
где
1 ( 4) 1 3 ( 2) 2
ce + c e1
. (20)
2 1
8
α x = Tx
( p = 1 − 2) .
Далее в виде сумм величин различных поряд~ ~
ков малости представим составляющие ортов k , j
~
и i . В частности, из выражений (2), (13) и (19) по~
лучим, что составляющие орта k можно записать
как
Разлагая уравнение (1) в степенной ряд и используя соотношения (14), (15), координату z точки
пересечения луча с поверхностью ДЛ можно представить как
z ( 2) =
(24)
Аналогично, разлагая в степенные ряды уравнения (10) и (11), можно представить в виде сумм
величин различных порядков малости направляющие косинусы нормалей фронтов записи:
⎫
⎪
⎬,
⎪⎭
e 2( 2) = x (1) T x(1) + y (1) T y(1)
⎫
⎪⎪
⎬ ; (23)
2
1 ( 4) ⎤ ⎪
− e1 ⎥
2
⎦ ⎭⎪
2
1
1
3
α (z2) = − e1( 2) , α (z4) = − e1( 4) + e1( 2) .
2
2
8
(15)
где
e1( 2) = T x(1)
1
α (y1) = Ty(1) , α (y3) = T y(3) − Ty(1) e1( 2)
2
1
⎡3
α (y5) = T y(5) − Ty(3) e1( 2) + Ty(1) ⎢ e1( 2)
2
⎣8
( )
и, используя соотношения (13), запишем их в виде
e1 = e1( 2) + e1( 4) + ... ⎫
⎪⎪
e2 = e2( 2) + e2( 4) + ...⎬ ,
⎪
e3 = e3( 2) + e3( 4) + ...⎪
⎭
⎫
⎪⎪
⎬ ; (22)
2
1
⎤
− e1( 4) ⎥ ⎪
2
⎦ ⎪⎭
( )
(13)
e1 = T 2 , e2 = T ⋅ S , e3 = S 2
1
α (x1) = Tx(1) , α (x3) = Tx(3) − Tx(1) e1( 2)
2
1
⎡3
α (x5) = Tx(5) − Tx(3) e1( 2) + Tx(1) ⎢ e1( 2)
2
⎣8
~
~
k z( 2) = −cz ( 2) , k z( 4) = −cz ( 4) .
(29)
Аналогично, раскладывая в степенные ряды
уравнения (4)-(6) и используя соотношения (26),
~ ~
представим составляющие ортов j и i в виде
~
⎫
jx = 0
⎪⎪
~
~ ( 2) ~ ( 4)
j y = j y + j y + ...
(30)
⎬
⎪
~ ~ (1) ~ (3) ~ (5)
j z = j z + j z + j z + ...⎪⎭
и
(21)
~
~
~
ix = 1 + ix( 2) + ix( 4) + ... ⎫
⎪⎪
~ ~ ( 2) ~ ( 4)
iy = iy + iy + ...
⎬.
⎪
~ ~ (1) ~ ( 3) ~ (5)
iz = iz + iz + iz + ...⎪⎭
(31)
27
Затем с помощью матрицы (7) находим составляющие различных порядков малости направляющих косинусов падающего луча и нормалей к фронтам записи. В частности, составляющие направляющих косинусов падающего на поверхность ДЛ
~
луча в K −системе можно записать как
~ =α
~ (1) + α
~ (3) + α
~ (5) + ...⎫
α
x
x
x
x
⎪
~ =α
~ (1) + α
~ (3) + α
~ (5) + ...⎪
α
⎬
y
y
y
y
⎪
( 2)
( 4)
~
~
~
α z = 1 + α z + α z + ... ⎪⎭
(32)
и аналогично для направляющих косинусов нормалей фронтов записи
~ =α
~ (1) + α
~ (3) + α
~ (5) + ...⎫
α
px
px
px
px
⎪
⎪
(1)
(3)
( 5)
~
~
~
~
α py = α py + α py + α py + ...⎬
~ = 1+ α
~ ( 2) + α
~ ( 4) + ... ⎪
α
⎪⎭
pz
pz
pz
( p = 1 − 2) .
(33)
Подставляя соотношения (32) и (33) в (9) получаем в виде сумм величин различных порядков малости составляющие дифрагировавшего луча в
~
K −системе
(34)
1.
2.
4.
5.
где
~ ′( q ) = α
~ (q) − α
~ (q) + α
~ ( q ) + ...⎫
α
x
x
⎪
1x
2x
~ ′( q ) = α
~ (q) − α
~ (q) + α
~ ( q ) + ...⎬
α
⎪⎭
y
y
1y
2y
(q = 1, 3, 5) ;
(35)
2
2
(1) 2
y
(1)
x
(3)
x
⎫
⎪⎪
⎬.
(1) ~ (3) ⎪
y αy ⎪
⎭
(36)
Далее, используя транспонированную матрицу
Α′ , находим в виде сумм величин различных порядков малости направляющие косинусы дифрагировавшего луча в K −системе. Наконец на последнем этапе с помощью соотношений, аналогичных
уравнениям (12), переходим от направляющих косинусов к направляющим тангенсам дифрагировавшего луча и тем самым находим в виде сумм составляющих различных порядков малости искомый
вектор T′ , определяющий наклон луча, дифрагировавшего на структуре ДЛ, размещенной на сферической поверхности.
3.
~′ = α
~ ′(1) + α
~ ′ ( 3) + α
~ ′(5) + ...⎫
α
x
x
x
x
⎪
~′ = α
~ ′(1) + α
~ ′(3) + α
~ ′(5) + ...⎪ ,
α
⎬
y
y
y
y
⎪
( 2)
( 4)
~
~
~
α′z = 1 + α′z + α′z + ... ⎪⎭
28
( ) + (α~ ) ⎤⎥⎦
( ) − α~ α~ − α~
~ ′( 2) = − 1 ⎡ α
~ (1)
α
z
x
2 ⎢⎣
~ ′( 4) = − 1 α
~ ′ ( 2)
α
z
z
2
6.
Литература
Степанов С.А., Грейсух Г.И. Расчет хода псевдолучей через оптические системы, включающие
градиентные и дифракционные линзы // Опт. и
спектр. - 1996. - Т.81, № 4. - C.698-701.
Степанов С.А., Грейсух Г.И. Компьютерный
расчет оптических систем, в области аберраций
высших порядков // Компьютерная оптика. – М.:
МЦНТИ, 1996. – Вып.16. - C.9-12.
Greisukh G.I., Bobrov S.T., Stepanov S.A. Optics
of Diffractive and Gradient-Index Elements and Systems. - Bellingham, WA: SPIE Press, 1997. - 414 p.
Welford W.T. A Vector raytracing equation for
hologram lenses of arbitrary shape // Optic Communication. – 1975. Vol.14, № 3. - P.322-323.
Бобров С. Т., Грейсух Г. И., Туркевич Ю. Г. Оптика дифракционных элементов и систем. - Л.:
Машиностроение, 1986. - 223 с.
Корн Г., Корн Т. Справочник по математике.
М.: Наука, 1977. – 832 с.
Документ
Категория
Без категории
Просмотров
4
Размер файла
210 Кб
Теги
структура, сферическая, дифракционной, выполненной, псевдолучей, расчет, поверхности, через, хода
1/--страниц
Пожаловаться на содержимое документа