菲涅尔多缝衍射的数值计算.doc

 第 24 卷 第 1 期 2011 年 2 月 大 学 物 理 实 验 P H YSICAL EXP ER IM EN T O F COLL E GE Vol . 24 No . 1 Feb . 2011 文章编号 :100722934 (2011) 0120055204 菲涅尔多缝衍射的数值计算 郭小花 ,王志坚 ,令维军 ( 天水师范学院 ,甘肃 兰州 741001) 摘 要 : 本文直接由菲涅尔2基尔霍夫衍射积分公式出发 ,通过数值积分的方法 ,计算菲涅尔单缝 、 双缝以及多缝衍射的光强分布 ,绘制出其相应的光强分布曲线 。同时给出了计算菲涅尔多缝衍射的光 振动分布的具体表达式 。研究结果对于理解菲涅尔衍射现象具有重要的理论参考意义 。 关 键 词 : 菲涅尔衍射 ;复振幅 ;光强 中图分类号 : O436 . 1 文献标识码 : A 光的衍射现象是光波动性的主要特性之一 。 以惠更斯22菲涅耳原理作为理论基础可以研究光 波的衍射现象 ,利用菲涅尔2基尔霍夫衍射积分公 式可以得到衍射光波场的光振动特点以及光强分 布 。但是光强分布的解析解在通常情况下需要通 过适当的近似处理 ,进行衍射积分计算 ,其过程繁 琐且复杂 ,难以得到准确的解析解 。对于夫琅和 费衍射 ( Fra unhof e r diff ractio n ,即远场衍射 ,光源 和接收屏幕距离衍射屏均为无穷远) ,通过简单的 计算可得到比较精确的衍射光强分布 ; 而对于菲 涅尔衍射 ( Fre snel diff ractio n , 即近场衍射 , 光源 和接收屏幕距离衍射屏均为有限远或其中之一为 有限远) ,一般是通过半波带法和矢量图解法等近 似方法可以定性或半定量地得到轴上一点的衍射 光强[ 1 ] ,轴外各点的衍射光强分析非常复杂 。目 前常见的是对于菲涅尔圆孔[ 225 ] 和圆屏[ 6 ] 等特殊 形状衍射屏[ 7 ,8 ] 的轴外光强的分析研究 , 而定 量 研究菲涅尔多缝衍射是波动光学中较困难的问题 之一 ,也较少见 。 本文直接从菲涅尔2基尔霍夫衍射积分公式 出发 ,运用 Mat la b 进行数值计算 ,通过数值积分 的方法 ,研究菲涅尔单缝 、双缝以及多缝衍射的光 强分布 ,并绘制其光强分布曲线 。 观察屏上的直角坐标系为 o x y , Z 为观察屏到单 缝的距离 。若单色平面光波垂直地入射到单缝上 , 则观察屏上一点 P 的复振幅 E ( x , y) [ 9 ] 为 : 图 1 单缝衍射坐标系 a/ 2 1 0 ) 2 i ( kz23π/ 4) ik ( x2x / 2 z ∫2 E ( x , y) = e λ e 2 d x0 ( 1) z a/ 式中 k = 2 π/λ, 改变式 1 的 x 值 , 利用 Mat la b 算 ( ) 出观察屏上的复振幅分布后再取模方或其共轭就 可得到观察屏上的相对光强分布 。 令波长λ = 632. 8 nm ,缝宽α = 0. 608 0 mm ,观 察离 衍 射 屏 的 垂 直 距 离 Z 分 别 取 0. 1460 m , 0. 073 0 m ,0. 048 7 m ,0. 036 5 m ,1. 460 4 m ,2. 0 m , 可得到如图 2 (a) ～ (f) 所示的光强分布曲线 。 由图 2 可知 , 在近场衍射中 , 衍射光主要分布 在与缝平行方向上 , 光强分布中出现了多个峰值 , 随距离 Z 的减小 , 峰值数增加 , 其分裂出的峰值个 数由 K = a2 / 4 Zλ[ 10 ] 决定 。随着距离 Z 的增大 , 衍 射光逐渐向缝的正交方向集中 , 光强分布也从“多 峰”变为“单峰”, 图形 ( 如图 2 ( a) ) 已开始向夫琅 禾费衍射靠近 , 如果再增大观察屏与衍射屏的距 1 单缝衍射 如图 1 所示 , 设 ox 0 y0 是位于单缝平面的直角 坐标系 , 单缝缝宽为 a , 且 y0 轴与单缝中线重合 , 收稿日期 : 2010209215 图 2 随距离 Z 变化的菲涅尔单缝衍射的光强分布曲线 随距离 Z 的逐渐增大 , 条纹越见清晰 , 光强极小值 接近于 0 , 当距离 Z 增大到 5 . 0 m 时 , 整个强度分 布曲线的包迹与单 缝衍 射 强度 分布 曲线 形 式一 样 , 即过渡到了标准的双缝夫琅和费衍射图样 。 2 双缝衍射 若用双缝衍射屏替换单缝衍射屏 , 设双缝对 称地分居 y0 轴两侧 , 双缝的缝宽都为 a , 缝距 ( 指 两狭缝中心的距离) 为 d 。并设有一束单色的平面 光波垂直照射双缝 ( 波长λ = 632 . 8 nm