CC BY
44
ДИФРАКЦИЯ СВЕТА НА ПРЯМОУГОЛЬНОМ ОТВЕРСТИИ СО СТАТИСТИЧЕСКИ НЕРОВНОЙ ГРАНИЦЕЙ
В.А. Куделькин, Ю.Л. Ратис Консорциум «Интегра-С», Самарский государственный аэрокосмический университет Институт систем обработки изображений РАН
Аннотация
В работе исследовано влияние нелинейности функции отклика и дифракционных поправок на интенсивность сигнала, выдаваемого оптикоэлектронным датчиком перемещений. Построены вычислительно эффективные аналитические выражения для расчета интенсивности светового потока. Показано, что учет дифракционных явлений позволяет существенно улучшить точность датчиковой аппаратуры.
Введение
При изготовлении и эксплуатации прецизионных оптикоэлектронных и электронно-оптических датчиков перемещений необходимо корректно оценивать нелинейность функции отклика и систематическую приборную погрешность измерений, возникающую из-за дифракции и несовершенства технологии изготовления датчиков [5, 7, 8, 9, 10].
Целью настоящей работы является оценка совместного влияния этих факторов на интенсивность регистрируемого сигнала в случае, когда диафрагма фотоприемника имеет прямоугольную форму.
Для этого рассмотрим задачу о дифракции плоской волны на отверстии, имеющем форму, близкую к прямоугольной. Это означает, что если прямоугольное отверстие имеет ширину a и высоту Ь , то для квазипрямоугольного отверстия со статистически неровной границей выполняются соотношения:
8a
«1,
«1,
(1)
где 8a = 8a(y) и 8b = 8b (x) - случайные функции.
I ! ! ! i
! 1 1
Рис. 1. Схема регистрации дифрагированного света В случае дифракции плоской волны на отверстии произвольной формы амплитуда дифрагированной волны Up в точке наблюдения Р равна [1]:
Ш
U.
=J
2niR
■ df„
(2)
где U - амплитуда падающей волны в окрестности отверстия, к - модуль волнового вектора падающей волны, df - площадь участка волновой поверхности, dfn - проекция вектора df на плоскость, перпендикулярную направлению падающего луча све-та1, а абсолютная величина радиус-вектора R в формуле (2) выражается через координаты (X, Y, Z) точки наблюдения P и координаты точек волновой поверхности (x, y, z):
R = V(X - x)2 + (Y - y)2 + (Z - z)2 . (3)
Соотношения (2)-(3) являются основой для формальной постановки последующих задач.
Дифракция на прямоугольном отверстии Пусть световая волна падает по нормали на идеальное прямоугольное отверстие в плоском экране, лежащем в плоскости xy . Если отверстие имеет
ширину a и высоту b , то координаты точек волновой поверхности удовлетворяют условиям:
-а/2 < x< a/2, • -b /2 < y < b/2, (4)
z = 0,
В этом случае амплитуда U падающей волны в формуле (2) одинакова для всех точек волновой поверхности, лежащей в плоскости экрана (U = U0 = const), соотношение (2) приобретает вид:
kU 2 2 Up = U •} dx J dy ■
2 2
QikJ ( X - x)2 + (Y - y)2 + Z 2
V( X - x)2 + (Y - y)2 + Z2
(5)
и амплитуда дифрагированной волны вычисляется
2
элементарно :
1 Вектор ¿¡[ направлен вдоль внешней нормали к волновой поверхности.
2 Обозначения и определения см. в приложении.
a
Up (X, Y, Z) * U0'е
\E,
! Ec
2i
f l± f £ - x " nZ V 2
+ Ef
f í - y
nZ l 2
Л
+ Еш
' l±f£+xЛ
. nZ l2 „
v v ' j
^ fí + Y ^
nZ v 2
(6)
При переходе от соотношения (5) к выражению (6) мы использовали стандартное разложение величины Я по малому параметру теории дифракции:
V(X - х)2 + (Y - y)2 + Z2 *
* Z +.!(X - x)2 + (Y - У)2
(7)
2 г
В случае наклонного падения волны амплитуда и = и (х, у) зависит от координат точек волновой поверхности, и ее нельзя выносить из-под знака интеграла. В частности, если на отверстии в экране дифрагирует плоская волна, то
и = и0 • ехр(/кО, (8)
амплитуда световой волны в точке наблюдения Р имеет вид:
Up
kU0' eik' 2niZ
2
: J dx'exp
kx + k'
(X - x)2 2Z
(9)
: J dy' exp
kyy + k ^ y 2Z
и выражается через интегралы Френеля:
(k2 + k?) Z ]
Up (X, Y, Z) • e 2i
J J i\ kxX+kyY+kZ-
U0 I x y 2k
(10)
< [ EF (X+) + EF (X-)]' [ EF (Y+) + EF (Y-)],
где
X± =
Y±=J— nZ
£ ±f X - M. _ 2 V k
í±fY KZJ 2
Y
k
(11)
В сферической системы координат
kx = k sin 6 cos ф, ky = k sin 6 sin ф,
kz = k cos 6,
(12)
где 6 - угол между направлением луча и плоскостью - ху , в которой лежит экран с прямоугольным
отверстием. Поэтому k2 + k2 = k2 sin2 6, и при 6 ^ 1
выполняется соотношение:
, 1 • 2 1 — sin2 ( 2
л—
' cos t
(13)
Следовательно, для случая малых углов падения световой волны на экран (6 ^ 1) с учетом соотношения (13) и того, что kz /к = cos6 , получаем:
(k2 + k2) Z _
kX + kY + kZ---y— * kR .
x y 2k
(14)
Таким образом, в приближении малых углов падения световой волны (6 ^ 1) амплитуда дифрагированной волны имеет вид:
ир (X, 7, г) = ^ • ехр\/(к • Я) ]х
2/ \ ] (15)
х [КЕ (X+) + Ее (XX-)] [е, (7+) + Ее (7 )].
Из (15) следует, что интенсивность света в точке наблюдения Р равна:
1р (х ,г, г) = -р (х ,г, г )|2 =
и2 (16)
ио I - |2 I ~ |2 4 '
=-0-\Ер (Х+) + Ер (XX- )| Е (7+) + Ер (7- )| .
Выражения (15) и (16) описывают дифракцию плоской световой волны на прямоугольном отверстии.
Дифракция на отверстии со статистически неровной границей Рассмотрим случай дифракции на прямоугольном отверстии со статистически неровной границей. Для этого представим вьгражение для амплитуды ир световой волны, дифрагированной на идеальном
прямоугольном отверстии в непрозрачном экране в следующем виде:
-Ро( X ,7, г) =
= ' J dx6f£-x|6
2ni J I 2 I
a
x + — 2
U(x,y, z)' J dy62- y^fy + 2
ei^V (X - x)2 + (Y - y)2 + Z 2
(17)
у1(X - х)2 + (7 - у)2 + г2'
Наблюдаемая амплитуда содержит вклад от дифракции на рельефе границы кодирующего элемента:
ир (х ,7, г) = — ? ах ? ау •и (х, у,о) •
•61 2 + 8ra(y)-xJ^lx + £ + 8,a(y)
•6[-2 + S„¿(x) - y' 6 (Vy + 2 + 5d¿(x)
íkV (X - x)2 +(Y - y)2 + Z2
(18)
X - x)2 + (Y - y)2 + Z2
a
2
Для вычисления амплитуды (18) воспользуемся свойствами обобщенных функций:
а + 8l - xUefa - x| + 8ia - x 18/
(19)
где / = а, Ь , 5/ = 5аг, 5а,, 5Ьи, 5Ьа .
Подставляя (19) в (18) и пренебрегая поправками второго порядка, получаем:
8 Up (X, Y, Z) = £8U,
(20)
Где явное выражение для каждой парциальной амплитуды 5 и,, будет приведено ниже.
Воспользуемся основной идеей теории дифракции (см. формулу (7), и представим первое слагаемое в разложении (20) в следующем виде:
к eikZ ik-A
8Uj =--U0e 2 ■
2ni Z
I (X-a/ 2)2
I 2z
2
J dy8ra (y)
ikyy+ik-
(Y - y )
(21)
Для оценки априорной приборной погрешности, связанной с дифракцией на статистически неровной границе диафрагмы, необходимо провести усреднение по статистическому ансамблю датчиков. В соответствии с этим введем величину (SU.). В этом
случае (8ra (y) = const, и этот множитель можно
вынести из под знака интеграла. В результате получаем:
Ь ТТ ikx—+ik\-
k U 0 JkZ x 2 \ 2Z
a, -J (X-a/2)2
<8U.) = -fe-2fe x[E,(Y+) + E,(Y )].
В рамках используемых приближений
kX ~ kY'
(22)
(23)
поэтому целесообразно ввести величину Акх = kX -кх. Тогда из (22) и (23) следует, что:
a .,14 X 2 +a2
<„,) .(Sa,) U .е- ■ e^'« (24)
<[ E, (Y+) + E, (Y-)],
и, аналогично,
I--,, a ,[ 4X2 +a2
IK т J +iAkX —+ik\-
(8U2) .(8a,) J— ■ ^eikZ ■ e 2 [ 8Z .. _ N 2/ X " VnZ 2i (25)
<[ E, (Y+) + E, (Y-)].
Введем характеристическую величину, описывающую флуктуации эффективной ширины щели для заданной длины падающей волны:
(8a(k, kx, X)) =
(8ar)e
-(8ai)
iAky — , X 2
(26)
Тогда из формул (24)-(26) вытекает, что (8Uj) + (8U2) . (8a(k,kx,X))x
'-Z • eikZU°•[E,(Y+) + Ef(Y_)].
(27)
Для второй пары границ отверстия аналогичным образом получаем:
(8U3)+ (8U^ «(8b(k,ky,Y))>
JkZ U0
где
2i [E, (X+ ) + E, (X- )] ,
(8b(k, ky, Y)) =
(28)
. , - iAky— . . iA
(8b„ ) e 2 + (8brf)e
(29)
Соотношения (26)-(29) дают полное решение задачи о нахождении амплитуды световой волны, дифрагированной на прямоугольном отверстии со статистически неровной границей.
Дифракция на узкой щели со статистически неровной границей
Рассмотрим узкую щель со статистически неровной границей. В этом случае поправкой {5из} + (5 и можно пренебречь, и выражение для
полной амплитуды световой волны в точке наблюдения принимает вид:
и
Up (X, Y, Z) = ■[ E, (Y+) + E, (Y-) ]>
x{e' (kR ■[ E, (X +) + E, (X-) ] + , k
+eikZ \8a(k, kx,X))■ '
(30)
Для оценки относительной априорной погрешности датчика перемещений необходимо рассчитать величину отношения 51 р /I р .
Из соотношений (16) и (30) следует, что:
Ip
h
e'(k ■R) + e'kZ
(8a(k, kx, X))
(31)
Ее (X+) + Ее (X-)
В рамках используемых приближений без ограничения общности можно ограничиться рассмотрением случая:
|(5а(к, кх,X))| « 1. (32)
Тогда, удерживая члены первого порядка малости по (5а), получаем:
i=1
ik
2
2
^ = 2Ке^ *
/1 кг-(к•Я) I
(да (к, кх, X ))|
пг ее (х+)+Ее (Х_ )
. (33)
Соотношение (33) позволяет дать априорную оценку минимальной неустранимой приборной погрешности показаний датчика перемещений.
Функция отклика оптикоэлектронного датчика перемещений
Функция отклика оптического датчика ) по определению имеет вид:
Н^,) = | ф(X, 7) 1р (X, 7, г, (34)
где ф(X, 7) - аппаратная функция, а й$п - проекция элемента светочувствительной поверхности датчика на плоскость, нормальную к направлению падения светового луча.
В рамках настоящей работы предполагается, что ф^, 7) = 1 на фотоэлементе датчика, и ф(X, 7) = 0 за его пределами. При этом а, и Ь, - ширина и высота входной диафрагмы фотоэлемента, соответственно, а Xi - величина смещения щели-осветителя относительно щели - диафрагмы фотоэлемента, т. е. искомое перемещение, измеряемое с помощью оп-тикоэлектронного датчика.
Для важного частного случая щелевидной диафрагмы фотоэлемента, входящего в состав датчика перемещений, заведомо выполняется условие Ь, » а,, и, следовательно, для интегралов Френеля можно использовать приближение |ее (7+) + Ее (7_ )| и 2. В
этом приближении функция отклика выражается через полную энергию светового потока, падающего на входную диафрагму фотоэлемента
Н( Xа) = К, (Xd):
Рис. 2. Типичная функция отклика оптикоэлектронного датчика
и 2 ^ +а/2 2
КX) =-2-ь, | \ЕгX )+Ее (XX-) dх. (35)
Xd-а/2
Введем обозначения:
"Ч^; "о=як\
к ( „ а, пг I " 2
"±=^\Xd |; Д" =
к кг
пг к
(36)
В данных обозначениях функция отклика принимает вид:
К (X,) =
Ко
I Е ("о +" - Д")+ее ("о - " - Д")|2 а", пг
(37)
где К, = и ¡Ь,
Для дальнейшего анализа удобно ввести безразмерную функцию отклика (Xd), связанную с оптической функцией отклика (X,) очевидным соотношением:
К(X,) = Ко •к,(X,).
(38)
Легко видеть, что (X,) выражается через вспомогательную функцию V("о,п):
К(X,) = 4[V("о,"+) - V("о,)] , (39)
где "± = "± -Д" , причем:
п
V ("о, п) =1Е ("о +")+ее ("о - ")|2 а", (4о)
-п
Выражение (39) представим в следующем виде:
V ("о, п) = Vо("о, п) -ДV ("о, п), (41)
где
Vо("о, п) = 2п • IЕе ("о + ") + Ее ("о - ")|2, (42)
Ч
ДV("о,п) = 21 Е("о,")"," .
(43)
Подынтегральная функция Е("о,") в формуле (43) выражается через интегралы Френеля:
Е ("о,") = [С ("о +") + С ("о -")]х
ОС8-("о + ")2 -ОС8-("о-")2
+[5 ("о +") + 5 ("о-")]
(44)
8Ш-2("о +")2 -Sin^("о -")2
Несложно показать, что поправочная функция ДV ("о, п) может быть представлена в виде линейной комбинации интегралов специально вида:
X
и
2(3
АГ (^ п) х (45)
а (45)
х {(а, Р) - В1 (а, -Р) - В2 (а, Р) + В2 (а, -Р)},
где а = — , Р = —, а интегралы В; 2 (а, Р) опреде-
2
лены ниже:
В1(а, Р) = | Ж б1п [а/(/ - 2 + 2Р)],
В2(а, Р) = I—б1п [а/(/ - 2 + 2Р)].
(46)
Интегралы В1(а, Р) элементарно выражаются через интегралы Френеля:
В1(а,Р) = С°5[а(Р-1)2] х
(0
х[ 5 ((1 + Р) ) + 5 ((1 -Р) ))-
-5ш[а([3-1)2] ,[С((0(1 + Р)) + с((0(1 -в))].
(0
(47)
Процедура вычисления интегралов В2(а, в) существенно сложнее. Проблема состоит в том, что прямой расчет интегралов (46) по квадратурным формулам Симпсона, Гаусса и т.п. возможен далеко не при любых значениях параметров а ив. Эта проблема является общей проблемой нахождения численных квадратур от быстро осциллирующих функций. Она связана с быстрым накоплением погрешности при численном сложении большого числа знакопеременных слагаемых, имеющих один порядок величины.
В целях преодоления проблемы накопления погрешности вычислений представим интеграл В2 (а, в) виде суммы:
В2(а, в) = А, (р, д) + Ас (р, д), (48)
где р = 4а, д = 4а (в -1), а вспомогательные функции А, (р, д) и Ас (р, д) определены соотношениями (49) и (50).
А, (р, д) = I—Бш(р/2)соБ(д/),
0 1
1Ж 2 Ас (р, д) = I— соб(р/ ) Б1п(д/). 0 /
(49)
(50)
Таким образом, проблема нахождения функции отклика оптикоэлектронного датчика сводится к задаче вычисления интегралов специального вида А, (р, д) и Ас (р, д). Рассмотрению этой проблемы посвящена отдельная работа [12].
Результаты расчетов На рис. 3, 4 представлены результаты численных расчетов интенсивности 1р и суммарного сигнала
Wd (X,) для случая падения зеленого света (X = 0,55 мкм) на бесконечно длинную щель шириной а = 50 мкм при наличии зазора между кодирующей шкалой и фотоэлементом ё = 50 мкм.
Из рисунков хорошо видно, что дифракционная картина носит ярко выраженный френелевский характер. Оптическая функция отклика (X,) нелинейна при малых Xd , что позволяет сделать вывод
об актуальности рассматриваемой в статье задачи, поскольку наибольшие относительные погрешности датчиковая аппаратура выдает именно для малых значений угловых или линейных перемещений, наиболее интересных с прикладной точки зрения.
0,8
0,6
0,4
0,2
-во -40 -20 0 20 40 60
Рис.3. Амплитуда световой волны на поверхности фотоэлемента
я:
/ X
1
0 а/2 а
Рис. 4. Реальная функция отклика оптикоэлектронного датчика
Заключение
Резюмируем вышесказанное следующим образом:
1. Полученные выражения имеют прозрачный физический смысл. Вся информация о влиянии качества механической обработки границы окна кодирующей шкалы оптикоэлектронного датчика на уровень сигнала содержится в функциях (5а(к, кх,X)) и (5Ь(к, ку, У)) .
2. Показано, что при анализе и синтезе оптикоэлек-тронных датчиков перемещений дифракция светового потока на кодирующей шкале является дифракцией Френеля.
0
3. Построены вычислительно эффективные аналитические выражения для расчета интенсивности светового потока и интегрального сигнала (Xа).
4. Проведен численный анализ совместного влияния дифракционных и инструментальных погрешностей на работу оптикоэлектронных датчиков перемещений.
5. Показано, что учет дифракционных явлений позволяет существенно улучшить точность датчи-ковой аппаратуры.
Представляется весьма перспективным использование разработанного формализма для решения задачи синтеза оптикоэлектронных датчиков перемещений с контролируемой нелинейностью оптической функции отклика.
Приложение
Определения и обозначения, используемые в работе
В работе используются стандартные определения интегралов Френеля:
C (z) = } cos (- t2 I dt,
S (z) = J sin 12 ) dt
(п1)
(п2)
быстрый алгоритм расчета которых построен в работах [3, 4, 6, 11].
Введем комплексный интеграл Френеля:
Ее (I) = С(I) + 5(I) = 2Ы (1,2,/п22j , (п3)
где М(а, Ь, г) - вырожденная гипергеометрическая функция.
Интегралы Френеля связаны со вспомогательными функциями /(г) и g (г) стандартными соотношениями:
C(z) =1 + f (z) sin I П z2 I - g (z) cos
S(z) = 2 - f (z) cos z21 - g(z)sin ^ 2 Из соотношений (п1) -(п5) следует, что
(п4) (п5)
Ef (z) = 2+i + [f (z) - ig (z) ]sin -[g (z) + if (z)] cos z 2 j.
z2 I-
(п6)
3 1
nzf(z) = 1 - 2 +...; nzg(z) = —-(nz ) nz
и, следовательно,
Ef (z) ж - —expI i—z1 | +.... f 2 nz 1 ~ 1
(п7)
(п8)
В асимптотике z » 1
Литература
1. Ландау Л. Д., Лифшиц Е.М., Курс теоретической физики. Теория поля. - М. Наука, 1973. - Т. 2. - 504 с.
2. Абрамовитц М., Стиган И. Справочник по специальным функциям. - М. Наука. 1979. - 832 с.
3. Ratis Yu.L., de Cordoba P.F. A code to calculate (high order) Bessel functions based on the continued fractions method. // Preprint University of Valencia FTUV 92-26 and IFIC 92-27, Burjassot, Spain, 1992.
4. Ratis Yu.L., de Cordoba P.F. A code to calculate (high order) Bessel functions based on the continued fractions method // Computer physics communications. 76 (1993) 381.
5. Ратис Ю.Л., Леонович Г.И. Дифракция светового потока на чувствительных элементах электронно-оптических и оптикоэлектронных датчиков механических перемещений // Компьютерная оптика, 1996. -№16. - С. 74-77.
6. Ratis Yu.L., Segura J., de Cordoba P.F. A code to evaluate Modified Bessel function based on the continued fractions method // Computer Physics Communications, v.105 1997. P. 263-272.
7. Ratis Yu.L., Leonovich G.I., Melnikov A.Yu. Light flux diffraction of fiber - optical and optical electronic transducers of mechanical displacement, Proceedings of SPIE, volume 3348 // Computer and Holographic Optics and Image Processing. 1997. - P. 336-343.
8. Ратис Ю.Л., Леонович Г.И., Мельников А.Ю. Быстрый алгоритм расчета интегралов специального вида при решении задач обработки оптических сигналов //Актуальные проблемы анализа и обеспечения надежности и качества приборов, устройств и систем: Тез. докл. МНТК 26-31 мая 1997 г. - Пенза, 1997. -С. 39-40.
9. Ратис Ю.Л., Леонович Г.И., Конюхов Н.Е., Мельников А.Ю. Фотоэлектрические преобразователи перемещений с линеаризованной функцией отклика оптического канала // Актуальные проблемы анализа и обеспечения надежности и качества приборов, устройств и систем: Тез. докл. МНТК 26-31 мая 1997 г. - Пенза, 1997. - С. 39-40.
10. Ратис Ю.Л., Леонович Г.И., Курушина С.Е., Мельников А. Ю. Нелинейные дифракционные искажения оптической функции отклика в кодирующих сопряжениях оптикоэлектронных датчиков // Компьютерная оптика. - Самара-Москва. 1998. - Т. 18.
11. Ratis Yu.L., Bastardo J.L., Abraham Ibrahim S., de Cordoba P.F., Urchueguia J.F.S., Evaluation of Fresnel integrals based on the continued fractions method // Applied Mathematics Letters, 18 (2005) - P.23-28.
12. Ратис Ю.Л. Куделькин В.А. Быстрый алгоритм расчета интегралов специального вида. Направлена в печать.
