Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

CALCULATION OF THE FRESNEL SYNTHESIZED HOLOGRAM

Zhumaliev K.M. 1 Ismanov Yu.Kh. 1 Alymkulov S.A. 1
1 Institute of Physics named after academician Z.Z. Zheenbaev of the National Academy of Sciences of the Kyrgyz Republic
The task of calculating the light fields in the Fresnel zone is considered. It is shown that in the general case this problem reduces to the calculation of the three-dimensional Kirchhoff integral. However, since this problem is practically not solvable analytically, it is proposed to use an approximation to calculate holograms in this zone, which reduces the Kirchhoff integral to the Fresnel transform. Therefore, to calculate the light fields in the so-called middle zone (Fresnel zone), the three-dimensional Kirchhoff integral is reduced to a two-dimensional Fresnel integral. In this paper, we propose an optimal method for such a transition, which allows minimizing the errors that arise. In addition, a discretization scheme has been proposed for the amplitude and phase components of the light wave, which makes it possible to effectively discretize the Fresnel integral. It is also shown that the discrete Fresnel transform can be reduced to a discrete Fourier transform, which allows the use of fast Fourier transform algorithms in the calculations of synthesized Fresnel holograms. When synthesizing holograms, it was assumed that the object is described by a fairly smooth function | g1(x,y) |, which is a reflection coefficient of intensity. Such a function can be restored at its discretization points by an interpolating operation using a certain function β(x, y).
hologram synthesis
wave amplitude
wave phase
discretization
Fresnel hologram
interpolation

Внешний вид любого объекта задают такие его характеристики, как его отражающие и рассеивающие свойства. Вид этих характеристик можно задать с помощью распределения коэффициентов отражения света, падающего на поверхность тела, по интенсивности G(x, y, z) и амплитуде g(x, y, z). Функциональная зависимость для амплитуды имеет вид

gumal03.wmf (1)

Модуль амплитудной составляющей gumal04.wmf и ее фазовая составляющая α задают фазовое изменение света, который падает после отражения от произвольной точки (x, y, z) рассматриваемого объекта. Функциональная связь между значениями G и g задается соотношением

gumal05.wmf (2)

где «*» представляет собой комплексно сопряженные величины.

Для распределения амплитуд и фаз светового поля вдоль поверхности исследуемого объекта можно записать следующее выражение gumal06.wmf. Распределение светового поля в этом случае для некоторой рассматриваемой поверхности задается в виде интегрального соотношения Кирхгофа

gumal07.wmf

Если выполняется условие

gumal08.wmf (3)

здесь φмакс – наибольшая величина угла, значение его задается в радианах, под которым виден объект с расстояния d, θ – некоторый коэффициент, задающий наибольшую фазовую ошибку, величина которой равна π/θ, то интегральное соотношение Кирхгофа можно преобразовать в интеграл Френеля [1–2]:

gumal09.wmf (4)

Соотношение (4) можно использовать для синтеза голограмм, которые так и называют – синтезированные голограммы Френеля.

Целью данной работы является рассмотрение процесса синтеза голограмм Френеля.

Дискретное представление голограммы Френеля

Рассмотрим дискретное представление голограммы Френеля:

gumal10.wmf (5)

Для ограниченного объекта (–Xмакс, Xмакс; –Yмакс, Yмакс), знание точек отсчета функции gumal11.wmf позволяет полностью ее восстановить методом интерполяции:

gumal12.wmf (6)

где

gumal13.wmf

Δξ = λd/2Xмакс, Δη = λd/2Yмакс. (7)

То есть функция F(ν, μ) полностью определена заданными точками отсчета [3–4]:

gumal14.wmf (8)

Подбор фазы может быть может быть осуществлен, например, путем освещения голограммы

gumal15.wmf (9)

волновым фронтом сферического вида с соответствующим радиусом кривизны

gumal16.wmf (10)

Посмотрим теперь, как найти отсчеты gumal17.wmf. Из (5) получаем

gumal18.wmf (11)

Объект можно описать зависимостью коэффициента отражения по интенсивности от координат точек поверхности gumal19.wmf. Точки дискретизации функции gumal20.wmf позволяют, так как функция достаточно гладкая, полностью ее восстановить методами интерполяции с помощью некоторой функции β(x, y) [5–7]:

gumal21.wmf, (12)

здесь Δx, Δy – шаги дискретизации вдоль координат x и y .

gumal22.wmf

gumal23.wmf (13)

здесь α1(x, y) – функция, пропорциональная фазовой составляющей отражающего коэффициента рассматриваемого объекта.

gumal24.wmf (14)

Добавим взаимно уничтожающие друг друга фазовые множители в дискретизированном виде, как под знак суммы, так и под знак интеграла в (13)

gumal25.wmf

gumal26.wmf

gumal27.wmf

gumal28.wmf (15)

Введем дополнительные условия:

gumal29.wmf gumal30.wmf (16)

Указанные дополнительные условия показывают, что погрешность в описании профиля объекта и отклонение формы волнового фронта от сферической формы были пренебрежительно малы, и потому можно было записать [8–10]

gumal31.wmf

gumal32.wmf

gumal33.wmf (17)

Интегральное соотношение в правой части (17) можно преобразовать следующим образом:

gumal34.wmf (18)

gumal35.wmf (19)

Сделав подстановку (18) и (19) в (15), получаем

gumal36.wmf (20)

Процесс суммирования по k и i в (20) осуществляется в пределах [–Xмакс/Δx, Xмакс/Δx] и [–Yмакс/Δy, Yмакс/Δy] соответственно.

gumal37.wmf можно определить как маскирующую функцию, так как вне пределов некоторого интервала по ν и μ она обращается в нуль. В случае интерполяции объекта в соответствии с теоремой отсчетов

gumal38.wmf (21)

Откуда

gumal39.wmf (22)

Из (22) можно видеть, что наибольшие значения r и s, для которых вычисляется сумма, задаются соотношениями [11–13]

gumal40.wmf (23)

Cогласно (7)

gumal41.wmf (24)

Таким образом, имеем

gumal42.wmf (25)

С целью устранения в (25) размерных величин, введем следующие обозначения:

gumal43.wmf (26)

Окончательно получаем

gumal44.wmf (27)

Соотношение (27) описывает дискретное преобразование Фурье. Вычисление такого преобразования требует помимо наличия матрицы g1(k, i), которая определяет распределение амплитуды света вблизи объекта, еще знание величин εν и εμ, которые задаются условиями (26) и характеризуют относительные размеры объекта, наблюдаемого из произвольной точки плоскости регистрации голограммы [14–15].

Как видно из сказанного, расчет синтезированной голограммы Френеля сводится к определению матрицы {F(r, s)} по матрице отсчетов объекта g1(k, i) и аналогового интерполирования полученных точек дискретизации в соответствии с (8)–(10).

Выводы

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