圆孔矩孔的菲涅尔衍射模拟工程光学.docx

XX大学

XX大学XXXX学院

PAGE

PAGE1

工程光学综合练习

圆孔、矩孔的菲涅尔衍射模拟

圆孔和矩孔的菲涅尔衍射模拟

一、原理

由惠更斯-菲涅尔原理可知接收屏上的P点的复振幅可以表示为

???(??)=??? ???(Q)exp(??????)K(θ)dσ(1)

??

??

其中???(Q)为衍射屏上的复振幅分布，K(θ)为倾斜因子。根据基尔霍夫对此公式的完善，有

C=1

????

K(θ)=1+????????≈1

2

1设衍射屏上点的坐标为（x1,y），接收屏上点的坐标为（x,y），衍射屏与接收屏

1

间距离为z1，当满足菲涅尔近似条件时，即

??

??????8??3[(?????1)2+(?????1)2]2

??????

???(2)

1

此时可得到菲涅尔衍射的计算公式

???(x,y)

=exp(??????1) ∞

?( )

????

[(?????

)2+(?????

)2]}????

????

(3)

iλ??

? ????1,??1 exp{2?? 1

1 1 1

1 ?∞ 1

把上式指数项中的二次项展开，并改写成傅里叶变换的形式，可以写成

???(x,y)=exp(??????1)exp[

????

(??2

+??2)]f{???(??

,??

)exp[

????

(??2

+??2)]}(4)

iλ??1 2??1

1 1 2??1 1 1

上式为菲涅尔衍射的傅里叶变换表达式，它表明除了积分号前面的一个与x1、y1无关的振幅和相位因子外，菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。

相对于夫琅和费衍射而言，菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中，输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场以及光强分布，考虑到这三个量都是二维分布,而且Matlab主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布,并分别以矩阵的列数和行数来对应平面的直角坐标值(x,y)以及(x1,y1)。

二、圆孔菲涅尔衍射

用MATLAB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同，考虑到运算量的要求，采样点数不能过多，所以每个屏的x和y方向各取200到300点进行运算。根据式（4），选取合适的衍射屏和接收屏尺寸和相距的

XX大学XXXX学院

距离，模拟结果如下：

取典型的He-Ne激光器波长λ=632.8nm，固定衍射屏和接收屏尺寸和相距的距离，分别取不同的圆孔半径，得到以下三组衍射图样，其圆孔半径分别为12mm，20mm，50mm

图1（r=12mm）

图2（r=20mm）

3

PAGE

PAGE4

XX大学XXXX学院

固孔形状 衍射回样

300

300

25D

.200

150

100

50

100

200

300

.-·.·. 0.25

1.5

1

o.s

忒- -

-2OJ -瓦|0

2[IJ

02

0.15

01

0.05

0。L· ·

10J 2[0 扣0

三、矩孔的菲涅尔衍射

图3(r=50mm)

步骤与上述相同，仅需改变与衍射屏形状对应的矩阵。这里选择矩孔的长宽相等，分别为15mm，20mm，30mm，其衍射图样及强度分布如图4、5、6

300

250

20口

150

100

5口

衍射焊形状

衍射后的图祥

10口 2[D

30□

...?.··

...?.··=.·.········

:i...·

.._

;

0.2

0.15

0.1

0.05

100

2□□

3口口

1.5

1

0.5

盐

勾 0

-200 -2DJ

图4（a=b=15mm）

XX大学

XX大学XXXX学院

PAGE

PAGE5

图5(a=b=20mm)

图6(a=b=30mm)

四、MATLAB程序

%所有长度单位为毫米lamda=632.8e-6;k=2*pi/lamda;z=1000000;

%先确定衍射屏

N=300; %圆屏采样点数a=15;

b=15;

[m,n]=meshgrid(linspace(-N/2,N/2-1,N));

I=rect(m/(2*a)).*rect(n/(2*b));q=exp(j*k*(m.^2+n.^2)/2/z);

若为圆孔，方框内替换为以下程序

r=12;a=1;b=1;

I=zeros(N,N);