任意正多边形小孔夫琅禾费衍射成像探讨.docVIP

任意正多边形小孔夫琅禾费衍射成像探讨 宋易知 北京师范大学物理学系 X 关注成功！ 加关注后您将方便地在 我的关注中得到本文献的被引频次变化的通知！ 新浪微博 腾讯微博 人人网 开心网 豆瓣网 网易微博 摘????要： 对夫琅禾费衍射积分进行了理论推导, 基于夫琅禾费衍射积分公式的一般形式, 推导出了一种对任意正多边形小孔夫琅禾费衍射情况的计算方法, 并利用Matlab软件进行计算机模拟, 做出了只用1个函数文件即可计算出所有可能的多边形小孔夫琅禾费衍射成像的尝试. 关键词： 正多边形小孔; 夫琅禾费衍射; Matlab; 作者简介：宋易知 (1998-) , 女, 四川乐山人, 北京师范大学物理学系2014级本科生. 收稿日期：2017-06-05 Discussion on the Fraunhofer diffraction imaging of arbitrary regular polygon hole SONG Yi-zhi Department of Physics, Beijing Normal University; Abstract： Based on the general form of Fraunhofer diffraction integral, the calculative method to analyze arbitrary regular polygon Fraunhofer diffraction was put forward.Using the Matlab software, computer simulation was carried out, attempting to get all possible Fraunhofer diffraction imaging of the regular polygon hole with just one function file. Keyword： regular polygon hole; Fraunhofer diffraction; Matlab; Received： 2017-06-05 1 理论基础 光的衍射现象能体现光的波动性, 即光在直线传播中遇到圆孔、小球、直边等障碍物时, 不被挡住反而绕到其后面传播, 进入光强本应为零的暗区[1-2].惠更斯-菲涅耳原理以波动理论解释光如何传播, 指出波阵面上每一点均为次波源, 这些次波的包络面即为新波面. 设衍射平面上有1束波长为λ的平面波垂直入射, 衍射平面到与之平行的观察平面之间的距离为z, 衍射平面上的P (x0, y0) 点与观察平面上的Q (x, y) 点间的距离为r.U0 (x0, y0) 和U (x, y) 分别表示P和Q点的光场, , 则有 其中C为比例系数, K (θ) 为倾斜因子[当θ增大时, K (θ) 缓慢减小], Σ为平面波的1个波面[3]. 索末菲使用恰当的格林函数, 从光的波动方程得到瑞利-索末菲衍射公式及, K (θ) =cosθ的结论. 而cosθ精确值为, 则 其中, 距离. 在通常的光学处理系统中, z应远远大于衍射小孔和观察区域的最大线度, 且从P到Q点的距离应远远大于光的波长, 即rλ, 相当于用平面波的振幅替代球面波前在观察平面上产生的复振幅[3].易推知θ应为小量, 则有 虽然 (3) 式分母中的r可以直接用z代替, 但不能忽略的是, 当λ很小时, k值很大, 指数中的r的微小误差可能引入远大于2π的相位变化, 使之不可以直接用z代替.对其做二项式展开[4] 保留展开式的前两项, 有 于是得到了菲涅耳衍射积分 从二项式展开式中弃去高次项所引入的误差越小, 则菲涅耳近似的精度越高.要求该误差所表示的最大相位变化, 即距离z满足, 此区域称为菲涅耳衍射区[5]. 考虑另一近似以便于计算.采用夫琅禾费近似, 即, 那么在整个孔径上近似等于0, 使得x0和y0的平方项的影响可忽略.因此在夫琅禾费衍射区内有 即在夫琅禾费衍射区内观察平面上的夫琅禾费衍射积分公式的一般形式. 2 计算原理 2.1 计算思路 如图1所示, 由于正多边形的中心对称性, 对于正L边形的小孔, 只需要计算出1个小三角形的衍射场分布, 将其进行L-1次坐标旋转变换, 随后将其全部相加, 就能获得任意正多边形小孔的夫琅禾费衍射成像.此前的衍射成像模拟程序都局限在某个特定形状小孔上, 而在这种思路的指导下, 本文做出了只用1个函数文件计算出所有可能的多边形小孔夫琅禾费衍射成像的尝试. 图1 正多边形由多个小三角形组成 ??下载原图 2.2 等腰三角形孔的夫琅禾费衍射 设衍射屏是等腰三角形孔, 如图2所示, 其顶角为θ, 底边长为a.单位