1、陕西理工学院毕业论文题 目 矩形环夫琅禾费衍射规律的研究 学生姓名 任强 学号 1210014024 所在学院 物 理 与 电 信 工 程 学 院 专业班级 物理学（物理 1201） 指导教师 潘峰 完成地点 陕西理工学院 2016 年 6 月 6 日陕西理工学院毕业论文第 1 页 共 11 页矩形环夫琅禾费衍射规律的研究任强（陕西理工学院 物理与电信工程学院 物理学专业 1201 班，陕西 汉中 723000）指导老师：潘峰摘要夫琅禾费衍射是指光源和观察幕离障碍物均为无穷远的衍射现象，又称为远场衍射。文章讨论了不同边长及波长的矩形孔夫琅禾费衍射的光强分布，并用 Mathematica 软件对
2、矩形环的夫琅禾费衍射进行了模拟。在软件中建立衍射的数学模型，对光强分布进行计算并仿真，最后以图形的方式直观的输出计算结果。即可得出矩形环夫琅禾费衍射的规律。关键词 光的衍射；矩形环夫琅禾费衍射；Mathematica；计算机仿真引言光在传播的工程中，遇到障碍物时，会发生绕过障碍物偏离直线传播，部分光线会进入障碍物的阴影区域，并在该区域出现明暗相间（光强分布不均匀）的现象，把这种现象称之为光的衍射 1。光的衍射现象经过惠更斯-菲涅尔原理后，通过基尔霍夫衍射公式得到衍射光强的分布规律。由于光的夫琅禾费衍射便于计算，同时在现代光学应用和光学成像系统方面应用较为重要。随着计算机技术的迅速发展，利用计算
3、机对光学现象进行仿真也成为一种可能。Mathematica是一款集数值分析、符号运算、图形处理、系统仿真等功能于一体的科学与工程计算软件，它具有编程效率高、简单易学、人机交互好、可视化功能、拓展性强等优点，利用 Mathematica 编程仿真光学现象只需改变程序中的参数，就可以生成不同实验条件下的光学图像，使实验效果更为形象逼真。在理论研究中，能快速的验证实验理论。本文将利用计算机软件 Mathematica 对光的夫琅禾费衍射规律和现象进行模拟，进一步分析矩形孔下光的夫琅禾费衍射现象 1。1 惠更斯-菲涅尔原理在波前上任何一个未受阻挡的点都可以看作是一个频率（或波长）与入射波相同的子波源；
4、在其后任何点的光振动，就是这些子波叠加的结果。如图 1.1 所示，S 为点波源，M 为从 S 发出的球面波在某时刻到达的波面 M，P 为波场中的某个点。把 M 面分割成无穷多的面元 dM ，把每个面元 dM 看成发射次波的波源，从所有面元发射的次波将在 P 点相遇 2。一般情况下，由各面元 dM 到 P 点的光程是不同的，从而在 P 点引起的振动位相不同，P 点的总振动就是这些次波在这里相干叠加的结果。图 1.1 面积元 dM 发出次波惠更斯菲涅耳原理可以表述如下：陕西理工学院毕业论文第 2 页 共 11 页波前上每一个面元都可看成是新的振动中心，它们发出次波（频率与入射波相同） ；在空间某一
5、点 P 的振动是所有这些次波在该点的相干叠加。如图 1.1 所示。点光源 M 在波面 M上任一点 Q产生的复振幅为，(1.1)()ikRQAEe式中，A 是离点光源单位距离处的振幅，R 是波面 M的半径。在 Q 点处取面元 ，面元发出的子波在 P 点产生的复振幅与在面元上的复振幅、面元大小和倾斜因子 K()成正比。面元 在 P 点产生的复振幅可以表示为，(1.2)()()ikRredcKdK()表示子波的振幅随面元法线与 QP 的夹角 的变化。 （ 称为衍射角）c 为一常数,r=QP。菲涅耳假设：当时 = 0，倾斜因子 K 有最大值，随着增加 ，K 减小，当 /2 时，K =0。对 P 点产生
6、作用的将是波面 M中界于 ZZ范围内的波面 M 上的面元发出的子波。则，(1.3)此即为惠更斯菲涅耳原理的菲涅耳表达式，结合(1.1)式此关系式还可推广为下式，即，(1.4)2 菲涅尔-基尔霍夫衍射积分公式对于单色点光源 S 发出的球面波照明无限大不透明屏上孔径 M 的情况，计算 P 点的场值。若孔径线度比波长大，但比孔径到 S 和 P 的距离小得多。则由亥姆霍兹一基尔霍夫积分定理选取包围 P 点的闭合曲面，它由三部分组成孔径；不透明屏右侧 M1 ；以 P 为中心，R 为半径的部分球面 M2 。则 P 点的场强值， (2.1) 对于 M1 和 M2 面，基尔霍夫假定， 在孔径 M 上， 和 的
7、值由入射波决定，与不存在不透明屏时完全相同。即，E(2.2)(2.3)表示外向法线与从 S 到上某点 Q 的矢量之间夹角的余弦。 在不透明屏右侧 M1 上，假定，(2.4)假定称为基尔霍夫边界条件：对于 M2，在 M2 上，R = r,并 则对 M2 上的积分关系，(2.5)n1cos,iklEeAnliklcos,nl 0cos,1nl222241ikRikRiiMikREedneeid()()ikRikrcAedE()ikrQK124()ikrikrMEedn陕西理工学院毕业论文第 3 页 共 11 页为 M2 对 P 点所张的立体角。在辐射场中，(2.6)而 是有界的则 时，可不考虑 M
8、2 的贡献，R(2.7)将(2.3)和 带入上式，并考虑到 1/r、1/l 比 k 值小的多。则，(2.8)此即为菲涅尔-基尔霍夫衍射公式。图 2.1 光经过孔径 M 的情况3 夫琅禾费衍射公式菲涅尔-基尔霍夫衍射公式在一定条件下及时对一些极简单的衍射问题，也会因为被积函数形式复杂而得不到解析形式的积分结果。为此，必须根据实际情况做进一步近似处理-旁轴近似 3。在一般的光学系统中，对成像起主要作用的是那些鱼光学系统光轴夹角极小的旁轴光线，而通常情况下，衍射屏开孔的线度和观察屏上的考檫范围都远小于开孔到观察屏的距离，因此，可以作以下两个近似 4：(1) 取式(2.8)中 ,因为倾斜因子 k(Q)
9、表示次波的振幅在各个方向上是不同的，其值介于 0 和 1。如果一平行光垂直入射到开孔上，则 ， ，因此 。(2) 式(2.8)中的 r 表示的是 QP 的距离，由于 很大，古可得 r=Z1 这样，式(2.8)就可以简化为，(3.1)图 3.1 矩形孔的衍射当单色平面光波垂直照射矩形孔时，会因为平面光波距离圆孔距离远近的不同产生不同的衍射图样。因此，可将衍射区划分为衍射效应可以忽略的集合投影区域、衍射效应不能忽略的近场衍射区，将其中的远场衍射又称为夫琅禾费衍射 5。 lim0REnike 214()ikriedPcos,ikrenikr2cos,iklrMnAecos()1n()cos/21kc
10、os(,)n11()ikrMEedZP,陕西理工学院毕业论文第 4 页 共 11 页如图 3.1 所示设 ，则有坐标关系可得， (3.2)对该式作二项式展开，有，(3.3)当 z1 大到使得上式及第三项引起的相位变化远远小于 时，即，(3.4)上面的第三项以及以后的各项都可以略去，简化为，(3.5)当观察屏离孔的距离很远时，使得（3.5）中的第四项满足，(3.6)进一步简化为，(3.7)这一近似称为夫琅禾费近似，在这个区域内观察到的衍射现象称为夫琅禾费衍射。在夫琅禾费近似下，P 点的光场复振幅为，(3.8)(3.8)即为夫琅禾费衍射公式 13。4 矩形环夫琅禾费衍射如图 3.1 所示的矩孔：(
11、4.1)取矩孔中心作为坐标原点,则观察屏上的 P 点的复振幅为，(4.2)平面波入射，(4.3)(4.4)对于轴上的点 P0,x=y=0,则其复振幅：(4.5)故，P 点（x,y）的复振幅为，(4.6)P 点的强度，2211()zxyr318()zk2211112xyzzrxy212211xyzzr211()(),ikikzxyxyzMeEed11:,:2baxy21(),abxyikfiklxyceEed2112sin(,)xyikfiklxiklyxyikfdabelcikfcA0cabE20sinxyikfleQPr陕西理工学院毕业论文第 5 页 共 11 页(4.7)(4.8)此即为夫
12、琅禾费矩形环衍射的强度分布公式 14。夫琅禾费矩形环衍射的强度分布公式讨论，x 轴上的点强度分布：此时 ，对应强度分布公式：(4.9)时， （P 0 点）有主极大 I/I0=1， 处，有极小值 I=0其零强度点（暗点）满足条件： (4.10)显然，相邻两个零强度点间的距离与宽度 a 成反比。此外，在两个零强度间有一强度次极大，其位置由，(4.11)决定。y 轴上点的强度分布，同 x 轴情况。中央亮斑可以认为是衍射扩展的主要范围，它的边缘在 x和 y 轴上分别由，(4.12)决定。用坐标表示，则有，(4.13)时，有几何光学结果。5 用 Mathematica 模拟矩形孔衍射图样夫琅禾费衍射的装
13、置入图 5.1 所示，平行光正入射到衍射屏上，衍射后面放一个凸透镜，在凸透镜的后焦平面上放一个接收屏，所接收的图像就是衍射光在无限远处形成的干涉图样（的微缩版）。图 5.2 是衍射屏坐标与衍射角坐标的示意图。根据夫琅禾费衍射理论，衍射场的振幅分布为 15，(5.1)其中 u(x,y)是衍射屏的透射函数，而 和 是方向坐标， 是光的波长。积分遍及整个衍射屏 M。图 5.1 夫琅禾费衍射装置 图 5.2 衍射的角坐标示意图*22020siniiklabIE,sin,isxylffklab020sinI,si1,x2in0dtgsin,sinxyab00ff02(sini)(,)(,)ixyMUuy
14、ed陕西理工学院毕业论文第 6 页 共 11 页他0.2 他0.1 0.1 0.25001001500200衍射强度的分布为，(5.2)把 U(,)分成实部和虚部来计算，(5.3)两个部分分别为，(5.4)(5.5)从理论上讲，按以上两个积分公式计算 U1、U 2 是没有问题的，因为积分函数 NIntergrate可以进行二重积分。但是，实际计算却发现，直线积分速度太慢。如果把积分改为求和，速度比积分要快许多，于是 12，(5.6)(5.7)矩形孔下的透射函数形式，(5.8)从 U1、U 2 的表达式中可以看出，它们都是 、 的函数，而在实际接收的衍射图上，用直角坐标系表示衍射图样强度比较方便
15、 6，将衍射强度改为接收屏上的直角坐标 x、y来表示。设在图 5.2中衍射屏到接收屏的距离为 1，则不难证明，(5.9)(5.10)通过 Mathematica 计算出的衍射图样如下 89， 波长为 0.6m，矩形环为边长为 10m 的正方形。程序没有定义 u 函数，因为取样点全部局限在矩形环内，幅度为 1。(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.3 波长为 0.6m，边长为 10m10m 的矩形环夫琅禾费衍射图样2(,)(,)IU12,i1 (sni)(,)(,)coMxyuxydx2 i,si1 (sini)(,)(,)cojkjkjk xyuxy2 2ii,
16、sijkjkjk1,(,)0jkjkxyuxy他22sin1xisy陕西理工学院毕业论文第 7 页 共 11 页他0.2 他0.1 0.1 0.25001001500200他0.2 他0.1 0.1 0.2100200300400500600700 波长为 0.6m，矩形环为边长为 20m 的正方形。(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.4 波长为 0.6m，边长为 20m20m 的矩形环夫琅禾费衍射图样 波长为 0.6m，矩形环为边长为 1020m 的长方形。(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.5 波长为 0.6m，边长为 1
17、0m20m 的矩形环夫琅禾费衍射图样 波长为 0.8m，矩形环为边长为 1010m 的正方形。(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.6 波长为 0.8m，边长为 10m10m 的矩形环夫琅禾费衍射图他0.2 他0.1 0.1 0.25001001500陕西理工学院毕业论文第 8 页 共 11 页他0.2 他0.1 0.1 0.21他1062他1063他1064他1065他1066他106 波长为 0.6m，矩形环为边长为 1060m 的长方形。(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.7 波长为 0.6m，边长为 10m60m 的矩形
18、环夫琅禾费衍射图样以上是对于矩形孔的衍射图样，再次基础上讨论矩形环的衍射图样如下：图 5.8 夫琅禾费衍射矩形环如图 5.8，设有边长为 2a 的正方形小孔，在孔内有边长为 a 的正方形不透光阴影，用单色平行光垂直照射该孔径，使其光斑投影到距离该孔径 z 的光屏上。孔径的振幅透射函数可以表示为两个矩形函数之差，即，(5.11)000(,)()()()2xyxytxyrecttrectrtaa孔径平面的透射光场为，(5.12)其傅里叶变换为，(5.13)2()2 20(,)4sinc()si()sinc()i()xyjfxyxyUxyafafaffeF此即为孔径的频谱，其 中， ，当 其沿 x
19、方向的频谱分布为，(5.14)22(,)sic()sic()xTfzz把 5.13 式代入夫琅禾费衍射方程，则得孔径的夫琅禾费衍射分布为，zy 00000(,)(,)()()()xyxtyrettretrcta陕西理工学院毕业论文第 9 页 共 11 页(5.15)222()0() 2()2 2 ()(,)(,)4sincsi()sinc()i()i()icisincxykjzjxykjzjxy jfxyxykjz jjxyeUxyUyafafaffee azzzz F ()z其对应光强分布为(5.15)2 2 2(,)16sinc()sic()sic()si()8i2naxyxyIxyazz
20、zzzxya用 Mathematica 画出式(5.15)的图形入下图所示，(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.9 波长为 0.6m，外边长为 2020m，内边长为 1010m 的矩形环衍射图样(a)三维衍射图 (b) 二维衍射图 (c) y 轴光强分布图图 5.10 波长为 0.6m，外边长为 2020m，内边长为 55m 的矩形环衍射图样6 结论(1) 明暗相间的衍射条纹与对应矩形的边平行。(2) 中央亮条纹光强度最大，光强度随着级数的增加，迅速减弱。(3) 由光强分布图可知，中央亮条纹宽度是各次极大亮条纹宽度的 2 倍。(4) 由图 5.3 和图 5.6 比较可知，波长越长衍射效应越显著；波长越短，衍射效应越不明显。(5) 由图 5.3 和图 5.4 比较可知，随着矩形环边长的增加，条纹变窄、变密。从衍射强度曲线也可以看出，随着矩形环边长的增加，其中央明纹的宽度在减小，衍射效果也逐渐的不明显。(6) 由图 5.3、图 5.5 和图 5.7 比较可知，随着矩形环宽度的减小，横轴方向能显示的条纹级数和位他20 他10 10 20510152025他20 他10 10 200.51.01.5
