Файл: Ersoy O.K. Diffraction, Fourier optics, and imaging (Wiley, 2006)(ISBN 0471238163)(427s) PEo .pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 910

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

5

Fresnel and Fraunhofer Approximations

5.1INTRODUCTION

Fresnel and Fraunhofer approximations of the scalar diffraction integral allow simpler Fourier integral computations to be used for wave propagation. They also allow different input and output plane window sizes. However, they are valid only in certain regions, not very close to the input aperture plane. The valid regions for the Rayleigh–Sommerfeld integral, Fresnel, and Fraunhofer approximations are shown in Figure 5.1.

The Rayleigh–Sommerfeld region is observed to be the entire half-space to the right of the input diffraction plane. The Fresnel and Fraunhofer regions are parts of the Rayleigh–Sommerfeld region. Approximate bounds indicating where they start will be derived below as Eqs. (5.25) and (5.27) for the Fresnel region and Eq. (5.41) for the Fraunhofer region. However, these bounds need to be interpreted with care. See Chapter 7 for further explanation.

The term far field usually refers to the Fraunhofer region. The term near field can be considered to be the region between the input diffraction plane and the Fraunhofer region.

This chapter consists of seven sections. Section 5.2 describes wave propagation in the Fresnel region. The FFT implementation of wave propagation in the Fresnel region is covered in Section 5.3. The fact that the Fresnel approximation is actually the solution of the paraxial wave equation is shown in Section 5.4. Wave propagation in the Fraunhofer region is discussed in Section 5.5.

Diffraction gratings are periodic optical devices that have many significant uses in applications. They also provide excellent examples of how to analyze waves emanating from such devices in the Fresnel and Fraunhofer regions. Section 5.6 discusses the fundamentals of diffraction gratings. Sections 5.7, 5.8, and 5.9 highlight Fraunhofer diffraction with a sinusoidal amplitude grating, Fresnel diffraction with a sinusoidal amplitude grating, and Fraunhofer diffraction with a sinusoidal phase grating, respectively.

Diffraction, Fourier Optics and Imaging, by Okan K. Ersoy

Copyright # 2007 John Wiley & Sons, Inc.

63

64

FRESNEL AND FRAUNHOFER APPROXIMATIONS

Rayleigh–Sommerfeld

Integral Region

A

Fresnel Region

Fraunhofer Region

Figure 5.1. The three diffraction regions.

5.2DIFFRACTION IN THE FRESNEL REGION

Let the input wave field be restricted to a radial extent L1:

Uðx; y; 0Þ ¼ 0 if

p

ð5:2-1Þ

x2 þ y2 > L1

Similarly, the radial extent of the observed wave field U(x,y,z) at the output plane is confined to a region L2 so that

q

Uðx0; y0; z0Þ ¼ 0 if x02 þ y02 > L2

ð5:2-2Þ

r01 in Eq. (4.6-3) is given by

r

 

¼

z 1

þ ð

x0

 

x

Þ

2 þ ðy0 yÞ2

#

1=2

5:2-3

Þ

 

 

 

 

 

 

 

z2

 

 

01

"

 

 

 

ð

where (x,y,0) are the coordinates of a point on the input plane and ðx0; y0; zÞ are the coordinates of a point on the observation plane. With the two restrictions discussed above, the following is true:

½ðx0 xÞ2 þ ðy0 yÞ2&max ðL1 þ L2Þ2

ð5:2-4Þ

The upper limit will be used below. If

jzj L1 þ L2;

ð5:2-5Þ

the term ði=jlÞðz=r012 Þ in Eq. (4.7-6) can be approximated by 1=jlz. However, more care is required with the phase. kr10 can be expanded in a binomial series as

kr10 ¼ kz þ 2kz ½ðx0 xÞ2 þ ðy0 yÞ2& 8kz3 ½ðx0 xÞ2 þ ðy0 yÞ2&2 þ . . .

ð5:2-6Þ


DIFFRACTION IN THE FRESNEL REGION

65

The third term has maximum absolute value equal to kðL1 þ L2Þ48jzj3Þ. This will be much less than 1 radian if

z

3

 

kðL1 þ L2Þ2

ð

5:2-7

Þ

 

j j

 

8

 

With this constraint on jzj, the phase can be approximated as

 

 

 

 

 

k

 

 

 

 

kr10 kz þ

 

½ðx0 xÞ2 þ ðy0 yÞ2&

ð5:2-8Þ

2z

The approximations made above are known as Fresnel approximations. The region decided by Eqs. (5.2-5) and (5.2-7) is known as the Fresnel region. In this region, Eq. (4.7-6) becomes

Uðx0

; y0

; zÞ ¼ jlz

ðð

Uðx; y; 0Þejlz½ðx0

xÞ þðy0

yÞ &dxdy

ð5:2-9Þ

 

 

 

ejkz

1

 

p

2

2

 

 

 

 

 

1

 

 

 

 

 

 

this is a 2-D convolution with respect to x and y and can be written as

Uðx0; y0; zÞ ¼ Uðx; y; 0Þ hðx; y; zÞ

ð5:2-10Þ

where the impulse response is given by

 

 

 

 

 

 

 

ejkz

 

p 2

2

 

 

hðx; y; zÞ ¼

 

ej

lz

ðx

þy

Þ

ð5:2-11Þ

jlz

The corresponding transfer function is given by

 

 

 

Hðfx; fy; zÞ ¼ e jkze jplzð fx2þfy2Þ

ð5:2-12Þ

The quadratic terms in Eq. (5.2-9) can be expanded such that Eq. (5.2-9) becomes

Uðx0

; y0

; zÞ ¼ jlz ej2zðx0

þy0Þ

ðð

U0ðx; y; 0Þe jlzðx0xþy0yÞdxdy

ð5:2-13Þ

 

 

 

ejkz

 

k 2

2

1

 

2p

 

 

 

 

 

 

 

 

 

1

 

 

 

 

where

U0ðx; y; 0Þ ¼ Uðx; y; 0Þej

k

ðx2þy2Þ

ð5:2-14Þ

2z


66

FRESNEL AND FRAUNHOFER APPROXIMATIONS

Aside from multiplicative amplitude and phase factors, Uðx0; y0; zÞ is observed to be the 2-D Fourier transform of U0ðx; y; 0Þ at spatial frequencies fx ¼ x0=lz and fy ¼ y0=lz.

EXAMPLE 5.1 Determine the impulse response function for Fresnel diffraction by starting with Eq. (4.4-21) of the ASM method.

Solution: When z r in Eq. (4.4-21), the second and third terms on the right hand side can be approximated by 1. Also expanding ½z2 þ r2&1=2 in a Taylor’s series with two terms kept, Eq. (4.4-21) becomes

ejkz

hðx; y; zÞ ¼ jlz ejlpzðx2þy2Þ

which is the same as Eq. (5.2-11).

EXAMPLE 5.2 Determine the wave field in the Fresnel region due to a point source modeled by Uðx; y; 0Þ ¼ dðx 3Þdðy 4Þ.

Solution: The output wave field is given by

Uðx; y; zÞ ¼ Uðx; y; 0Þ hðx; y; zÞ

dðx

3Þdðy

 

1

ejkz 1þ

x2þy2

 

4Þ

 

2z2

jlz

1

ejkz 1þ

x 3 2þðy 4Þ2

 

 

 

 

 

 

ð Þ 2z2

 

 

 

 

 

jlz

 

 

 

 

 

EXAMPLE 5.3 A circular aperture with a radius of 2 mm is illuminated by a normally incident plane wave of wavelength equal to 0.5 m. If the observation region is limited to a radius of 10 cm, find z in order to be in the Fresnel region.

Solution: With L1 ¼ 0:2 cm and L2 ¼ 10 cm, Eq. (5.2-4) gives

z 10:2 cm

If a factor of 10 is used to satisfy the inequality, we require z 1:02 m

The second condition given by Eq. (5.2-7) is

z3 2pð0:102Þ4m3 4 10 6

z 5:4 m

Using a factor of 10, we require

z 54 m


DIFFRACTION IN THE FRESNEL REGION

67

It is observed that the second condition given by Eq. (5.2-7) yields an excessive value of z. However, this condition does not have to be true for the Fresnel approximation to be valid [Goodman]. In general, kr01 is very large, causing the quadratic phase factor ejkr10 to oscillate rapidly. Then, the major contribution to the diffraction integral comes from points near x x0 and y y0, called points of stationary phase [Stamnes, 1986]. Consequently, the Fresnel approximation can be expected to be valid in a region whose minimum z value is somewhere between the two values given by Eqs. (5.2-4) and (5.2-7). As a matter of fact, the Fresnel approximation has been used in the near field to analyze many optical phenomena. The major reasons for this apparent validity are the subject matter of Chapter 7.

EXAMPLE 5.4 Determine the Fresnel diffraction pattern if a square aperture of width D is illuminated by a plane wave with amplitude 1.

Solution: The input wave can be written as

 

x

y

Uðx; y; 0Þ ¼ rect

 

rect

 

 

D

D

The Fresnel diffraction is written in the convolution form as

 

ejkz

D=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ð ð ej

p

2

 

 

 

 

 

 

2

 

 

 

 

 

 

Uðx0; y0; zÞ ¼

 

 

½ðx x0Þ

þðy y0

Þ

&dxdy

 

 

 

lz

 

 

jlz

 

 

 

 

D=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ejkz

2

1

D=2

ej

p

ðx x0Þ2 dx3 2

 

1

D=2

 

ð

lz

 

ð

¼ j

 

 

 

4

 

 

 

 

 

 

 

 

 

 

5 4

 

 

 

6

 

D=2

 

 

 

 

 

 

 

 

7 6

 

 

D=2

Consider the integral

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ix ¼ p1lz

Dð=2

ej

p

ðx x0Þ2 dx

 

 

 

 

lz

 

 

 

 

 

 

 

 

D=2

 

 

 

 

 

 

Let

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n ¼

2

 

ðx x0Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lz

 

Then, Ix becomes

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

xe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 2

 

 

 

 

 

 

 

 

 

Ix ¼

p2

 

ð ej2v dv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xb

 

 

 

 

 

 

7

ð

 

Þ

ej

p

ðy y0

Þ2 dy3

 

5:2-15

 

lz

 

 

 

 

 

5

 

 

 

ð5:2-16Þ