圆孔矩孔的菲涅尔衍射模拟matlab实现工程光学 (2)

工程光学综合练习——－——圆孔、矩孔的菲涅尔衍射模拟

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

一、原理

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

其中为衍射屏上的复振幅分布,为倾斜因子.根据基尔霍夫对此公式的完善，有

设衍射屏上点的坐标为（ｘ1, y1）,接收屏上点的坐标为（x，y),衍射屏与接收屏间距离为ｚ1，当满足菲涅尔近似条件时，即

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

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

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

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

）。

１

二、圆孔菲涅尔衍射

用MATＬＡB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同,考虑到运算量的要求，采样点数不能过多,所以每个屏的x和ｙ方向各取2０0到３00点进行运算.根据式（4）,选取合适的衍射屏和接收屏尺寸和相距的距离,模拟结果如下:

取典型的Ｈe－Nｅ激光器波长λ=6３2．8nm,固定衍射屏和接收屏尺寸和相距的距离,分别取不同的圆孔半径,得到以下三组衍射图样,其圆孔半径分别为

12mm,20mm，50mｍ

图错误!未定义书签。（r＝１２mｍ)

图错误!未定义书签。（r＝２0mｍ)

图错误!未定义书签。（r=5０mm）三、矩孔的菲涅尔衍射

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

图错误!不能识别的开关参数。（a=b=１5mｍ）

图错误!不能识别的开关参数。（a=b=２0ｍm)

图 错误!不能识别的开关参数。(a=ｂ=30mm ）

四、MAT ＬＡＢ程序

％所有长度单位为毫米

ｌa ｍd ａ=632.8e-6; k=２*pi/la ｍda ；

z ＝10０００00;

%先确定衍射屏

Ｎ=3０0; %圆屏采样点数

ａ=1５；

ｂ=15；

[m ，n ］=mes ｈgr ｉd(lins ｐace （-N ／2，Ｎ／２

-１,N)）;

q=exp(j*k*(m 。^2+ｎ。＾2)/2/z ）；

s ｕｂp ｌo ｔ(２,2,1）； ％圆孔图像画在2行2列的第一个位置

ima ｇesc(I) %画衍射屏的形状

col ｏr ｍa ｐ（[0 0 ０； 1 1 1］） ％颜色以黑白区分 ａxis i ｍag ｅ

ｔi ｔl ｅ('衍射屏形状’)

若为圆孔，方框内替换为以下程序 r=12;a=1;b=1; I=zeros(N,N); [m,n]=meshgrid(linspace(-N/2,N/2-1,N)); D=((m-a).^2+(n-b).^2).^(1/2); i=find(D<=r); I(i)=1; %孔半径范围内透射系数为1

L=300;

M=３00；％取相同点数用于矩阵运算［x,ｙ］=meshgｒｉd(linspace（—L/２，L/2,M）)；

h=exp(j*k*z）*exp(（ｊ*k*(x。^２+y.^2）)／(2*z))/（j＊ｌamda*z);%接收屏

H=fｆtｓhｉft（fft2（ｈ))；

B=fｆｔｓhiｆt(fｆｔ2(I）)；%圆孔频谱

G=H.*B; %公式中为卷积,空间域中相卷相当于频域中相乘

Ｕ= fftshift(iｆft2（Ｇ)）; ％求逆变换,得到复振幅分布矩阵

Br=（U／ｍax(U))；%归一化

suｂｐｌot(2，2,2);

imshｏw（aｂs(U)）;

axｉs imagｅ;

colｏrｍap（ｈot)

％ｆiguｒe,ｉmｓhow(C);

titｌｅ(’衍射后的图样’）；

subploｔ（2,2，3);

ｍｅsh(x，ｙ,abｓ(U)）；％画三维图形

subploｔ（2，2,４)；

plｏt（ａbs（Ｂr)）