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

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

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

一、原理

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

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

设衍射屏上点的坐标为(x1, y1),接收屏上点的坐标为(x, y),衍射屏与接收屏间距离为z1,当满足菲涅尔近似条件时,即

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

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

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

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

二、圆孔菲涅尔衍射

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

取典型的He-Ne激光器波长λ=632、8nm,固定衍射屏与接收屏尺寸与相距的距离,分别取不同的圆孔半径,得到以下三组衍射图样,其圆孔半径分别为

12mm,20mm,50mm

图1(r=12mm)