РАСЧЕТ ХОДА ПСЕВДОЛУЧЕЙ ЧЕРЕЗ ОПТИЧЕСКИЕ СИСТЕМЫ, ВКЛЮЧАЮЩИЕ ДИФРАКЦИОННЫЕ ЛИНЗЫ, СТРУКТУРА КОТОРЫХ ВЫПОЛНЕНА НА АСФЕРИЧЕСКОЙ ПОВЕРХНОСТИ
Г. И. Грейсух1, Е. Г. Ежов2; С. А. Степанов1 'Государственная архитектурно-строительная академия, г. Пенза 2Институт систем обработки изображений РАН, г. Самара
Методика расчета хода псевдолучей через оптические системы, включающие радиально-градиентные и дифракционные линзы, распространена на системы с дифракционными линзами, структура которых выполнена на асферической поверхности.
Введение
В работах [1-3] описана методика и приведены формулы расчета хода псевдолучей лучей (то есть лучей, траектории которых рассчитываются в приближении заданного порядка малости) через враща-тельно-симметричные оптические системы, содержащие радиально-градиентные элементы и дифракционные линзы (ДЛ), структура которых выполнена на плоской поверхности. В работе [4] указанная методика распространена на случай включения в систему ДЛ, структура которых размещена на сферической поверхности. Данная работа направлена на дальнейшее развитие методики и обобщение ее на оптические системы, содержащие. ДЛ, структура которых выполнена на асферической поверхности.
Расчет хода псевдолуча через такие ДЛ, включает две задачи: прослеживание хода псевдолуча через среду, которая в общем случае может быть ограничена двумя асферическими поверхностями, и определение параметров псевдолуча, дифрагировавшего на структуре ДЛ. Ниже получены формулы расчета, позволяющие решать обе эти задачи.
1. Расчет хода псевдолуча через однородную среду, ограниченную асферическими поверхностями
Ход луча будем описывать в декартовой системе координат, ось 2 которой совпадает с оптической осью. Высоту и наклон луча определим с помощью векторов р и е; при этом вектор р имеет
составляющие [х(г), у(г), 0], а вектор е - составляющие (е х, е у ,0), где е х и е у - направляющие
тангенсы луча, связанные с его направляющими косинусами соотношениями е х = а х/ а 2 и
еу = ау/аг.
При распространении луча между к и (к +1) асферическими поверхностями оптической системы луч на входе в среду зададим векторами р к и е к, а на выходе из среды - векторами рк+1 и е к+1. Если 2к, к+1 есть расстояние вдоль оси 2 от точки входа луча в однородную среду до точки его выхода, то
очевидно, что векторы рк между собой уравнениями
рк+1 =рк + 2к,к+1е к е к+1 =е к
е к и р к+1 , е к+1 связаны
(1)
Пусть ёк - расстояние между вершинами к и (к +1) поверхностей. Тогда
2
к ,к+1= ^к + 2к+1 2к
(2)
где 2к - координата точки пересечения луча с к поверхностью в системе координат, связанной с вершиной этой поверхности и, аналогично, 2к+2 -координата точки пересечения луча с (к +1) поверхностью в системе координат, связанной с вершиной (к +1) поверхности.
Уравнение асферической поверхности в системе координат с началом в вершине этой поверхности запишем в виде
^ ( 2 )= С2 -1 + д/1 -(ср)2 -1 Стз (ср)4 -
о
- — ст5 (ср)6 -... = 0 16 у '
(3)
где с - кривизна поверхности в ее вершине, Стз , СТ5, ... - коэффициенты асферической деформации поверхности.
Тогда координаты 2к лить из выражений
1
и 2к+1 можно опреде-
2к = — ск
+
1 -д/1 -(ск рк)2
1 Стз,к (скрк)4 +
8
+ 77Ст5,к (скрк)6 +...
1 -"\/1 -(ск+1р к+1)2
16
2к+1 =
ск+1
+1 Стз, к+1 (ск+1р к+1)4 +
8
(4)
+"¡т СТ5, к+1 (ск+1р к+1)6 +... 16
Уравнения (1), (2) и (4) являются основой для получения формул расчета хода псевдолуча через однородную среду, ограниченную асферическими поверхностями. При их выводе в качестве первого шага представим векторы р к и е к в виде сумм слагаемых первого, третьего и пятого порядков малости относительно модулей векторов, определяющих высоту и наклон луча во входном зрачке оптической системы:
р _п(1) ,_(3) + р(5) + ,
рк ~Рк + рк +рк + -I 8 к =в к1) +о(3) +„(5)
+ 8), у + „>, у +...
где р®, рк3 и рк5 - составляющие первого, третьего и пятого порядков малости вектора, определяющего положение луча. Обозначения „к1"1 - „к5"1 имеют
аналогичный смысл для вектора, определяющего наклон луча.
Из (2), (4) и (5) нетрудно видеть, что расстояние 7к к+1 можно представить в виде суммы членов
нулевого и четных порядков малости:
_ 7 (0)
2
к ,к+1
_ 2
к ,к+1
г (2) к ,к+1
7 (4)
к ,к+1
(6)
Подставляя (5) и (6) в первое из уравнений (1), получим, что на выходе из среды
р _р(1) +р (3) р к+1 _рк+1 +р к+1
где
р(3) _р(3 рк+1 _рк
7(°)
+ 7
(0)
к ,к+1° к
+ р
.(1) •к
(3)
+...
(7)
+ 7;
(2) „(1)
к,к+1Ь к
-»(5)
(5) + 7(0) „(5)+ 7 (2) „(3)
рк+1 _рк + 7к,к+1„к + 7к,к+1„ к +
(8)
+ 7(4) ет
+ 7к ,к+1Ь к
Для того чтобы найти слагаемые различных порядков малости расстояния гк,к+1, разложим
уравнения (4) в степенные ряды, а результат разложения подставим в (2). При этом воспользуемся уравнениями (7) и (8), введем три инварианта вращения
е1 к, е2 к рк, ез _р к
(9)
и, используя соотношения (5), представим эти инварианты в виде
е1 _ е(2) + е<4) +..;
е2 _ е™ + е24) +... •, (10)
е3 _ е32) + е34) +...
где
,(2)
,(2)
,(2)
„(1) лк
.(!) „(!)
е(4) _ 28а> -„(3) е1 _ 2„ к „ к
--ру „к е(4) _ р(3)
е2 _рк 2
(1) +р(1) „(3)
„/ +
рк „к'
ркХ)
,(4)
^-рк3)
(11)
В результате получим
7(0) _ а 7к,к+1 _ ик
7 (2) _ к ,к+1
_ 1 (ск+1 - ск)
е32) +
+ с
к+1" к 2
ЛЛ е22) + - а^Р
(12)
(13)
7 (4)
к ,к+1
2 (ск+1 - ск )е: +8 [(+а3,к+1 )с!+1 -
Ос 3 ][
8
+ °3, к)
+ ск+1 7 (2)
а к |е24) + 2 | +
+ 7к'к+1( + аке1
(2)
(2)
-1 (1-о \
))с1+11 ¿к [2)
++ст3, к+1
+ 4[ ] + 2е22)е32)-+ 2 + еР)ак2е22)
Уравнения (8), (11) и (12-14) позволяют по известным составляющим первого, третьего и пятого порядков малости параметров луча на входе в однородную среду вычислить соответствующие составляющие параметров луча на выходе из нее, то есть эти уравнения позволяют рассчитать ход псевдолуча пятого порядка через однородную среду, ограниченную двумя асферическими поверхностями оптической системы. При этом каждая из этих поверхностей может представлять собой как поверхность оптического элемента, так и любую другую поверхность, например поверхность предмета или изображения.
2. Расчет хода псевдолуча через дифракционные структуры, выполненные на асферической поверхности Как уже отмечалось во введении, в работе [4] получены формулы расчета хода псевдолуча через ДЛ, структура которой выполнена на сферической поверхности. Распространение этих формул на асферическую поверхность требует лишь получения выражений для величин различных порядков малости единичного вектора нормали к к этой поверхности в точке падения луча, определяемой вектором р.
Единичный вектор нормали к асферической поверхности можно найти, воспользовавшись выражением [5]:
к )2 . (15)
Подставив (3) в (15), получим, что составляющие вектора к описываются соотношениями
' к ^
V ~ I
кг _ 11 +
_ -Нск,
(-р)2 ]
12
(16)
X
где
1 -
ы2 ]2 + -2 стэ ы2
+ Ö СТ5 (сР)4 + -
(17)
Для того чтобы представить составляющие кх,
ку и кг в виде сумм величин различных порядков
малости, вновь используем один из трех инвариантов вращения
e3 =р2 = + е(4) +...
где
,32) = [р(1)
е34) = 2р(1) -р(3).
(18)
(19)
Подставляя (18) в (16) и (17), а затем, раскладывая радикалы в ряд, получим
V ~ J
~.(3) ~.(5)
+ +
У У У
+...
kz = 1 + k
(2)
где
= —c
= —c<
= —c<
+k
x(1)
(4)
+...
(20)
_ У
x(3) ,(3)
x(5)
_ У
,(5)
(2) + b((2) + k}2)
x(1) У(1)
x(3)
У
(3)
+((4) + »(2) +k}4))
x(1) У(1)
, (21)
kp =—_ c~e-z 2
~(4) =— i c2 [4) + 2»(2)e32)
z 2 L3 3
2e(2) 3
3 4 + — c
(2)
(22)
I(2) = 2 (1+ C3 ) 2e32) I(4) = 1 (1+ C3 ) 2ef +
+ 3 (1 + > Z
,(2)
(23)
Соотношения (21-23) совместно с соответствующими формулами работы [4] позволяют произвести расчет хода псевдолуча через ДЛ, структура которой размещена на асферической подложке.
Литература
1. Степанов С.А., Грейсух Г.И. Расчет хода псевдолучей через оптические системы, включающие градиентные и дифракционные линзы // Опт. и спектр. 1996. Т.81, № 4. C. 698-701.
2. Степанов С.А., Грейсух Г.И. Компьютерный расчет оптических систем, в области аберраций высших порядков // Компьютерная оптика. М., МЦНТИ. 1996. В.16. C. 9-12.
3. G.I. Greisukh, S.T. Bobrov, S.A. Stepanov Optics of Diffractive and Gradient-Index Elements and Systems // Bellingham, WA: SPIE Press, 1997. 414 p.
4. Ежов Е.Г., Степанов С.А. Расчет хода псевдолучей через дифракционные структуры, выполненные на сферической поверхности // Компьютерная оптика. М., МЦНТИ. 2000. В. 20. C. 25-28.
5. Корн Г., Корн Т. Справочник по математике // М., Наука. 1977. 832 с.