XX 大学XXXX 学院 1 工程光学综合练习 -----圆孔、矩孔的菲涅尔衍射模拟 XX 大学XXXX 学院 2 圆 孔 和矩孔 的菲涅尔衍射模拟 一、原理 由惠更斯-菲涅尔原理可知接收屏上的 P 点的复振幅可以表示为 ᵃ̃(ᵄ) = ᵃ ∬ᵃ̃(Q) exp(ᵅᵅᵅ)ᵅᵯK(θ)dσ (1) 其中ᵃ̃(Q)为衍射屏上的复振幅分布，K(θ)为倾斜因子。根据基尔霍夫对此公式的完善，有 C = 1ᵅᵰ K(θ) = 1 + ᵅᵅᵆᵰ2≈ 1 设衍射屏上点的坐标为（x1, y1），接收屏上点的坐标为（x, y），衍射屏与接收屏间距离为 z1，当满足菲涅尔近似条件时，即 ᵅ8ᵆ13 [(ᵆ − ᵆ1)2 + (ᵆ − ᵆ1)2]ᵅᵄᵆ2≪ ᵰ (2) 此时可得到菲涅尔衍射的计算公式 ᵃ̃(x, y) = exp(ᵅᵅᵆ1)iλᵆ1∬ᵃ̃(ᵆ1, ᵆ1)∞−∞exp { ᵅᵅ2ᵆ1[(ᵆ − ᵆ1)2 + (ᵆ − ᵆ1)2]} ᵅᵆ1ᵅᵆ1 (3) 把上式指数项中的二次项展开，并改写成傅里叶变换的形式，可以写成 ᵃ̃(x, y) = exp (ᵅᵅᵆ1)iλᵆ1exp [ ᵅᵅ2ᵆ1(ᵆ2 + ᵆ2)] f {ᵃ̃(ᵆ1, ᵆ1) exp [ ᵅᵅ2ᵆ1(ᵆ12 + ᵆ12)]} (4) 上式为菲涅尔衍射的傅里叶变换表达式，它表明除了积分号前面的一个与x1、y1 无关的振幅和相位因子外，菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。 相对于夫琅和费衍射而言，菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中，输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场以及光强分布，考虑到这三个量都是二维分布,而且 Matlab 主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布,并分别以矩阵的列数和行数来对应平面的直角坐标值(x, y)以及(x1, y1)。 二、圆孔菲涅尔衍射 用 MATLAB 分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同，考虑到运算量的要求，采样点数不能过多，所以每个屏的 x 和 y 方向各取200 到 300 点进行运算。根据式（4），选取合适的衍射屏和接收屏尺寸和相距的XX 大学XXXX 学院 3 距离，模拟结果如下： 取典型的He-Ne 激光器波长λ=632.8nm，固定衍射屏和接收屏尺寸和相距的距离，分别取不同的圆孔半径，得到以下三组衍射图样，其圆孔半径分别为12mm，20mm，50mm 图 1（r=12mm） 图 2（r=20mm） XX 大学XXXX 学院 4 图 3(r=50mm) 三、矩孔的菲涅尔衍射 步骤与上述相同，仅需改变...
