基于Matlab的光学衍射仿真

1、I基于基于 Matlab 的光学衍射实验仿真的光学衍射实验仿真 摘摘 要要 光学试验中衍射实验是非常重要的实验. 光的衍射是指光在传播过程中遇到障碍物时能够绕过障碍物的边缘前进的现象, 光的衍射现象为光的波动说提供了有力的证据. 衍射系统一般有光源、衍射屏和接受屏组成, 按照它们相互距离的大小可将衍射分为两大类, 一类是衍射屏与光源和接受屏的距离都是无穷远时的衍射, 称为夫琅禾费衍射, 一类是衍射屏与光源或接受屏的距离为有限远时的衍射称为菲涅尔衍射。本文用 Matlab 软件对典型的衍射现象建立了数学模型，对衍射光强分布进行了编程运算，对衍射实验进行了仿真。最后创建了交互式 GUI 界面，用户

2、可以通过改变输入参数模拟不同条件下的衍射条纹。本文对于衍射概念、区别、原理及光强分布编程做了详细全面的介绍关键字关键字：Matlab；衍射；仿真; ;GUI 界面; ;光学实验IIMatlab-based Simulation of Optical Diffraction ExperimentAbstractOptical diffraction experiment is a very important experiment. is the diffraction of light propagation of light in the obstacles encountered in t

3、he process to bypass the obstacles when the forward edge of the phenomenon of light diffraction phenomenon of the wave theory of light provides a strong Evidence. diffraction systems generally have light, diffraction screen and accept the screen composition, size according to their distance from eac

4、h other diffraction can be divided into two categories, one is the diffraction screen and the light source and the receiving screen is infinity when the distance between the diffraction Known as Fraunhofer diffraction, one is diffraction screen and the light source or accept a limited away from the

5、screen when the diffraction is called Fresnel diffraction. In this paper, Matlab software on a typical phenomenon of a mathematical model of diffraction, the diffraction intensity distribution of the programming operation, the diffraction experiment is simulated. Finally, create an interactive GUI i

6、nterface, users can change the input parameters to simulate different conditions of the diffraction pattern. This concept of the diffraction, difference, intensity distribution of programming principles and a detailed comprehensive descriptionKey word: matlab；diffraction; simulation; gui interface;

7、optical experimentIII目目 录录1 绪论.11.1 光学仿真的研究意义.11.2 国内外研究现状.21.3 MATLAB仿真的优越性.21.4 仿真的主要内容.22 衍射.32.1 光的衍射现象.32.1.1 衍射定义.32.1.2 光的衍射现象.32.2 惠更斯费涅耳原理.62.2.1 原理表述.62.2.2 原理的定量表达式.62.3 夫琅禾费原理.72.3.1 夫琅禾费衍射的装置.82.3.2 夫琅禾费矩孔衍射.92.3.3 夫琅禾费单缝衍射.102.3.4 夫琅禾费多缝衍射.112.3.5 多缝衍射图样.122.4 菲涅尔衍射原理.132.4.1 菲涅尔半波带法.1

8、32.4.2 菲涅尔单缝衍射.142.4.3 矩孔菲涅尔衍射.153 夫琅禾费衍射仿真.163.1 夫琅禾费单缝衍射仿真.173.2 夫琅禾费多缝衍射仿真.193.3 夫琅禾费矩孔衍射仿真.204 菲涅尔衍射仿真.274.1 菲涅尔方孔衍射仿真.234.2 菲涅耳单缝衍射仿真.265 交互式 GUI 界面.296 总结.30IV参考文献.31致 谢.33毕业设计（论文）知识产权声明.34毕业设计（论文）独创性声明.35附录 1（GUI 编程）.361 绪论11 绪论绪论1.1 光学仿真的研究意义光学仿真的研究意义在工程设计领域中，人们通过对研究对象建立模型，用计算机程序实现系统的运行过程和得到

9、运算结果，寻找出最优方案，然后再予以物理实现，此即为计算机仿真科学。在计算机日益普及的今天，计算机仿真技术作为虚拟实验手段己经成为计算机应用的一个重要分支。它是继理论分析和物理实验之后，认识客观世界规律性的一种新型手段。计算机仿真过程是以仿真程序的运行来实现的。仿真程序运行时，首先要对描述系统特性的模型设置一定的参数值，并让模型中的某些变量在指定的范围内变化，通过计算可以求得这种变量在不断变化的过程中，系统运动的具体情况及结果。仿真程序在运行过程中具有以下多种功能(l)计算机可以显示出系统运动时的整个过程和在这个过程中所产生的各种现象和状态。具有观测方便，过程可控制等优点; (3)借助计算机的

10、高速运算能力，可以反复改变输入的实验条件、系统参数，大大提高实验效率。因此，计算机仿真具有良好的可控制性(参数可根据需要调整)、无破坏性(不会因为设计上的不合理导致器件的损坏或事故的发生)、可复现性(排除多种随机因素的影响，如温度、湿度等)、易观察性(能够观察某些在实际实验当中无法或者难以观察的现大幅度节省实验所耗费的人力物力，特别是在一些重复实验工作强度较大且对实验器材、实验环境等要求较苛刻的情况下，如在大型激光仪器的建造过程中，结合基准实验的仿真计算结果可为大型激光器的设计和优化提供依据。仿真光学实验也可应用于基础光学教学。光学内容比较抽象，如不借助实验，学生很难理解，如光的干涉、菲涅耳衍

11、射、夫琅禾费衍射等。国外著名的光学教材配有大量的图片(包括计算和实验获得的图片)，来形象地说明光学中抽象难懂的理论。光学实验一般需要稳定的环境，高精密的仪器，因此在教室里能做的光学实验极为有限，而且也受到授课时间的限制。为了克服光学实验对实验条件要求比较苛刻的缺点，可采用计算机仿真光学实验，特别是光学演示实验，配合理西安毕业设计（论文）2论课的进行，把光学课程涉及的大多数现象展示在学生面前，以加深对光学内容的理解。如利用计算机仿真联合变换相关实验，可以得到清晰的相关峰，而在实验中液晶光阀的分辨率较低，很难得到清晰的相关峰;又如光学菲涅耳衍射与夫琅禾费衍射它们之间的演化规律，清楚地说明二者之间的

12、联系与区别。学生们可以根据对光学原理和规律的理解，自己设置在仿真光学实验中的可控参数，探索和发现光学世界的奥秘，调动学习的积极性。1.2 国内外研究现状国内外研究现状 在科学计算方面，国外的光学实验仿真是在模拟设计和优化光学系统的过程中发展起来的。在这方面，美国走在最前面，其中最具代表性的是劳伦斯利弗莫尔实验室光传输模拟计算软件 Prop92 及大型总体优化设计软件CHAINOP 和 PROPSUI 法上有独到之处，主要体现在其快速傅里叶变换的计算效率很高:软件采用特殊方法能够处理小于计算分辨率的灰尘点的衍射过程以及截止频率小于计算网格分辨最小频率的滤波过程等。另外，该软件图形显示界面友好，运

13、行稳定。我国用于科学研究的光学实验计算机数值仿真软件虽开发较晚，但也己经取得了显著成绩。特别是年，神光一川原型装置 TIL 分系统集成实验的启动为高功率固体激光驱动器的计算机数值模拟的研究创造了条件。目前己基本完成的 SG99 光传输可靠，模拟计算软件的开发，推出的标准版本基本能稳定运行，对 SG99 主要计算模块的验证结果表明 SG99 对能流放大、线性传输、非线性传输的计算是合理可靠的，其中线性传输的计算模块的计算精度与国外同类软件 Fresnel 相当;目前该软件已经应用于神光一主机可行性论证的工作中。在光学教学方面，国外己有相关的配有光盘演示光学实验的教材，该教材主要针对高年级学生和研

14、究生使用。其中不仅详尽的介绍了几何光学、物理光学、光学成像技术及图像处理技术，而且利用现在普遍使用的软件工具 Matlab对它们进行了系统的仿真。也有针对理科和工科低年级学生使用的光学教材，该教材使用 Matchcad 绘制各种逼真的光学仪器，创造出仿真的光学实验室，学生可利用其进行探索和发现性学习，充分调动学生的积极性。还有网络版光学教材，该西安毕业设计（论文）3教材采用进行光学仿真计算，结合 LiveGraPhic3DJaval.1 的动画制作功能在网络上实时演示各种光学实验的结果图。我国光学教材在利用计算机仿真方面相对落后，至今没有同类教材出现。在 2003 年北京举行的网络教育软件展上

15、，有关光学实验的网络教学软件都偏重于理论分析方面，对计算机应用于光学实验的仿真方面未给与充分重视。结合国家十五教材建设计划，在光学实验仿真方面进行大量的研究，各项研究工作将在后续各章中一一介绍。1.3 Matlab 仿真的优越性仿真的优越性 Matlab 是 Mathworks 公司于 1982 年推出的一套高性能的数值计算和可视化软件。它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便、界面友好的用户环境。它还包括了 ToolBox 江具箱)的各类问题的求解工具，可用来求解特定学科的问题。其特点是:(l)可扩展性:Matlab 最重要的特点是易于扩展，它允许用户自行建立指定功能

16、的 M 文件。对于一个从事特定领域的工程师来说，不仅可利用 Matlab 所提供的函数及基本工具箱函数，还可方便地构造出专用的函数，从而大大扩展了其应用范围。当前支持 Matlab 的商用 Toofbox(工具箱)有数百种之多。而由个人开发的 Toolbox 则不可计数。(2)易学易用性:Matlab 不需要用户有高深的数学知识和程序设计能力，不需要用户深刻了解算法及编程技巧。(3)高效性:Matlab 语句功能十分强大，一条语句可完成十分复杂的任务。如ffi 语句可完成对指定数据的快速傅里叶变换，这相当于上百条 C 语言语句的功能。它大大加快了工程技术人员从事软件开发的效率。据 Mathwo

17、rks 公司声称，Matlab 软件中所包含的 Matlab 源代码相当于 70 万行 C 代码。由于 Matlab 具有如此之多的特点，在欧美高等院校，Matlab 已成为应用于线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具;在研究单位、工业部门，Matlab 也被广泛用于研究和解决各种工程问题。当前在全世界有超过 40 万工程师和科学家使用它来分析和解决问题。1.4 仿真的主要内容仿真的主要内容 本课题主要培养学生进行光学设计以及计算机仿真的综合能力。光的衍射现象是光学重要物理现象之一，在大学物理课程学习中占有重要的地位，用计算机对光衍射

18、现象的模拟是对其物理本质更好的理解和补充。本课题使用Matlab 软件结合所学的物理光学中光的衍射原理，对夫琅禾费衍射实验和菲涅西安毕业设计（论文）4尔衍射的光强分布进行编程运算，包括了单缝，多缝以及矩空缝宽条件下，并输出计算得到的衍射图样分布，对实验现象进行仿真。最后做成了用户可以通过改变不同的输入参数条件下就模拟出不同的衍射实验的 GUI 交互式界面。设置的计算参数观察仿真图样的变化规律，给出物理光学理论解释。本课题涉及到光学知识，计算机仿真等知识内容的综合运用。2 衍射52 衍射衍射2.1 光的衍射现象光的衍射现象2.1.1 衍射定义衍射定义当波传播过程中遇到障碍物时，波就不是沿直线传播

19、，它可以到达沿直线传播所不能达到的区域。这种现象称为波的衍射现象（或绕射现象） （原因是波阵面受到了限制而产生的） 。(1) 理解衍射现象的两个要点:光波的波面可以看作是连续分布的次波源;次波源发射的次波满足相干条件,观察场中衍射光强的重新分布是次波相干叠加的结果(2)衍射现象的特点: 光束在什么方向受限制,衍射图样就沿什么方向扩展.光束被限制得越厉害,衍射图样越扩展,衍射效应越强./a 1/1000 时，衍射现象不明显；1/100 /a a,则沿着 y 轴较沿 x 轴的暗点间距为密，在 x轴和 y 轴各点的光强度，要根据它们的坐标进行计算，从上面的分析我们不难明白，强度为零的地方是一些和矩孔

20、边平行的直线，亦即平行于 x 轴和 y 轴的直线，在两组正交暗线形成的一个个矩形格子内，各有一个亮斑。可以看出，中央亮斑的强度最大，其他亮斑的强度比中央亮斑要小得多，所以绝大部分光能集中在中央亮斑内。中央亮斑可认为是衍射扩展的主要范围，它的边缘在 x和 y 轴上分别由条件 决定。若以坐标表yxbasinsin和示，则有 fbyfax00,可见，衍射扩展与矩孔的宽度成反比，而与光波波长成正比。当孔宽时，衍射扩展趋于零，衍射效应可以忽略，所得结果与几何光学的结果一致。所以，在几何光学可以看成是波长的极限情况。02.3.3 夫琅禾费单缝衍射夫琅禾费单缝衍射如果矩孔一个方向的宽度比另一个方向的宽度大得

21、多，比如 ba，矩孔就变成了狭缝。单缝的夫琅禾费衍射，由于入射光在 y 方向的衍射效应可以忽略，衍射图样只分布在 x 轴上。显然，单缝衍射在 x 轴上的衍射光强分布公式也是 在衍射理论中通常称为单缝衍射因子。矩孔衍射的相对强度是两个单缝衍射因子的乘积。根据前面的讨论，可知在单缝衍射图样0II中，中央亮纹是在下式决定的两个暗点范围内： fax0这一范围集中了单缝衍射的绝大部分能量。在宽度上，它是其他亮纹的两倍。在单缝衍射实验中丝测径仪来精确测定金属丝或者纤维丝的直径。因为直径为 a 的细丝和不透明屏上的距离为 a 的单缝可看成是一对互补屏，所以应用了巴俾涅原理很容易找到细丝衍射图样和单缝衍射图样

22、的关系。在单缝衍射的讨论中，已经知道，衍射条纹的间距（相邻两暗纹之间的距离西安毕业设计（论文）11 faxe因此，直径为 a 的细丝的衍射条纹间距也有上式表示。在实际测量中，只要测量出细丝的衍射间距间距，便可以由上式计算细丝的直径。目前已把细丝测量仪的生产过程做连续的动态监测。 2.3.4 夫琅禾费多缝衍射夫琅禾费多缝衍射多缝夫琅禾费衍射装置如图所示，图中 S 是与图面垂直的线光源，位于透镜的焦面上；G 是开有多个等宽等间距狭缝（缝宽为 d）的衍射屏，多缝的1L方向与线光源平行。多缝的衍射图样在透镜的焦面上观察。假如多缝的取向2L是方向，那么很显然，多缝衍射图样的强度分布只沿着 x 方向变化，

23、衍射条1y纹是一些平行于 y 轴的亮暗条纹。 多缝衍射图样的强度分布同样应该用夫琅禾费衍射公式进行计算，这时积分区域是多个狭缝露出的波面。不过，我们也可以利用上节得到的结果来简化计算，无须逐个缝进行积分运算。在方向上两个相距为 d 的平行等宽狭缝在 P 点产1x生的复振幅有一位相差，而单个 P 点产生的振幅为sin2d 因此，P 点光强为 式中,是单缝在点的光强2202sin2sinsinNaaII200EI 0P西安毕业设计（论文）12度。上式便是 N 缝衍射的强度分布公式。容易看出，当 N=2 时，上式就是双缝衍射的强度公式。式中包括两个因子：单缝衍射因子和多束光干涉因子，2)sin(aa

24、表明多缝衍射也是衍射和干涉两种效应共同作用的结果。单缝衍射因子只与单内引入的振幅和位相的变化）有关而多光束干涉因子来源于狭缝的周期性排列，与单缝本身的性质无关。因此，如果有 N 个性质相同的缝在一个方向上的周期排列起来，或者 N 个性质相同的其他形状的孔径在一个方向上周期性排列起来，它们的夫琅禾费衍射图样的强度分布式中就将出现这个因子。这样，只要把单个衍射孔径的衍射因子求出来，将它乘上多光束干涉因子，便可以得到这种孔径周期排列的衍射图样的强度分布。这个规律对于求多个周期排列的孔径的衍射是很有用的。2.3.5 多缝衍射图样多缝衍射图样多缝衍射图样中的亮纹和暗纹位置可通过分析多光束干涉因子和单缝衍

25、射因子的极大值和极小值条件得到。当时，它有极大mdsin.2, 1, 0m值，其数值为。这些极大值称为主极大。当2N Nmmdsin1,.2 , 1,.;2, 1, 0Nmm时，它有极小值，其数值为零。不难看出，在两个相邻主极大之间有 N-1 个零值。相邻两个零值之间（）的角距离，相邻两个主极大与相邻一个1m零值之间的角距离而是，所以主极大的半角宽度为 表明缝数 N 越大，主极大的缝宽越小。此为，在相邻两个零值之间也有一个极大值。这些极大值叫做次极大，它们的强度比主极大要弱得多。可以证明，次级大的强度与它离开主极大的远近有关，但主极大旁边的最强的次极大，其强度也是只有主极大强度的 4%左右。显

26、然，次极大的宽度也随 N 增大而减小，当 N 是一个很大的数目时，它们将于强度零点混成一片，成为衍射图样的背影。对应 4 个缝的干涉因子的曲线，这时在两相邻主极大之间有 3 个零点，2 个次极大。可以看出，与双缝衍射的情况类似，各级主极大的强度也受到单缝衍射因子的调制。各级主极大的强度为西安毕业设计（论文）13202sinaaINIm它们是单缝衍射在各级主极大位置上产生的强度的倍。其中零级主极大的强2N度最大，等于。如果对应于某一级主极大的位置，那么该02IN0sin2aa级主极大的强度也降为 0，该级主极大就消失了，我们知道这就是缺级。缺级的规律如上述，还可以看出，各级主极大的相对强度与缝数

27、 N 无关，它只依赖于缝距 d 与缝宽 a 之比。2.4 菲涅尔衍射原理菲涅尔衍射原理菲涅尔衍射是在在菲涅尔近似成立的距离上观察到的衍射现象。相对于观察夫琅禾费衍射而言，观察菲涅尔衍射是在离衍射屏比较近的地方。衍射屏上圆孔直经为 2cm，光波波长 600nm，这时为满足菲涅尔近似，要求观察屏到衍射屏的距离大于 25cm，而菲涅尔衍射的一般装置中，S 是点光源，K 是开有某种形状孔径的衍射屏（也可以是一个很小的不透明屏） ，M 是观察屏，在距离衍射屏不太远的地方。通常光源离衍射屏的距离都要比衍射屏上的孔径大得多，为处理简明起见，可以认为光源发出的光波垂直照射在孔径上。在某些特别需要精确的情况下，

28、可以不用这一假设，但处理方法完全相同。2.4.1 菲涅尔半波带法菲涅尔半波带法考察单色平面波垂直照射圆孔衍射屏的情形，我们利用菲涅尔波带法来决定点的光强度，光强度位于通过圆孔中心 C 且垂直于圆孔平面的轴上。假0P0P设单色平面波在圆孔范围内可以按照如下方法：以为中心，以0P.为半径分别做出一系列球面，每个球面都与相交成圆，而11,2zz则被划分一个环带，在这些环带中，两相邻带的相应点到点的光程差为半个0P波长，这些环带因此叫做菲涅尔半波带或菲涅尔波带。显然，点的复振幅就0P西安毕业设计（论文）14是波面上所有波带发出的子波在点产生的复振幅的叠加。由惠更斯-菲涅尔0P原理得知，各个波带在点产生

29、的振幅正比与该带的面积，反比于该带到的0P0P距离，因此，第 j 个波带（圆心 C 所在的为第 1 波带，向外依次为第 2， ，j,波带）在点产生的振幅可以表示为0P 2cos1jjjrACE式中 C 是比例常数，是 j 个波带到点的距离，是第 j 个波带的面积。这jr0PjA样一来，各波带在差为半波长，它们发出的子波到达点产生位相差为 。因0P此，若把奇数波带在点产生的复振幅的位相为零，则偶数波带在点产生的0P0P复振幅的位相就是 ;相邻波带产生的复振幅分别为一正一负。这样，个波带在产生的总振幅总和为0P 321.EnEEEE )(022n22111为偶数为奇数）（nEEEEEEnn 另一方

30、面，对于一定的圆孔大小和光波波长，波带数 n 取决于点的距离0P，即不同的点对应不同的波带数 n。因此，当把观察屏沿光轴 C平移1z1z0P0P时，同样可以看到点忽明忽暗地交替变化。利用菲涅尔衍射的计算公式可以0P证明，点的光强随的变化是大所致） 。因此表明这时点的复振幅等于第0P1z0P1 个波带的复振幅的一半，强度为第 1 波带产生的强度的 1/4。由此可见，当圆孔包含的波带的数目很大时，圆孔的大小不再影响点的光强度。这实际上也0P是从光的直线传播定律出发所得出的结论。所以我们可以说：从波动概念和从光的直线传播概念得出的结论，当圆孔包含的波带的数目很大时开始吻合。2.4.2 菲涅尔单缝衍射

31、菲涅尔单缝衍射衍射装置如图所示。单缝宽度为 a，缝长，缝长方向平行于轴。 当选取坐1y西安毕业设计（论文）15标原点 C 通过单缝中心时，观察屏上的复振幅分布为： 这就是单缝的菲涅尔衍射的公式。它表示，单缝菲涅尔衍射同样可以利用菲涅尔积分和科纽蜷线来计算。在科纽蜷线图上，上式上式大括号里面的两个复数差也有一个矢量表示，矢量起点在 终点在 。由于1122zaxw1222zaxw1122zawww对于一个特定的装置，它是常数，与 x 无关，所以不管考察观察屏上 x 坐标为何值的点，这个矢量两端点之间的曲线长度相等。这样一来，当矢量两端点在科纽蜷线上 w=0 附近（两端点位置取决于 x 值，当 x=

32、0 时，两端点对称位于原点两边） ，一般地矢量长度较短。不过，矢量实际的长短变化变化与缝宽 a 很有关系，不能一概而论。2.4.3 矩孔菲涅尔衍射矩孔菲涅尔衍射设矩孔在方向的宽度为 a，在方向的宽度为 b。选取矩孔中心为坐标1x1y原点，由式，得到矩孔衍射公式11112122211122211112222222212exp2expexp,zbyFzbyFzaxFzaxFiEdyyyzikdxxxzikziikzyxEbbaa该式表明，矩孔衍射图样的振幅（强度）分布是两个互相垂直的单缝衍射图样3 夫琅禾费衍射仿真163 夫琅禾费衍射仿真夫琅禾费衍射仿真光的衍射现象是光的波动性的重要表现之一.波动

33、在传播过程中,只要其波面受到某种限制,如振幅或相位的突变等,就必然伴随着衍射的发生. 然而,只有当这种限制的空间几何线度与波长大小可以比拟时,其衍射现象才能显著地表现出来. 所有光学系统,特别是成像光学系统,一般都将光波限制在一个特定的空间域内,这使得光波的传播过程进光学滤波器设计等具有非常重要的意义. 然而,由于光波波长较短,与此相应的复杂形状衍射屏的制作较困难,并且实验过程中对光学系统及环境条件的要求较高. 因而在实际的实验操作和观察上存在诸多不便. 计算机仿真以其良好的可控性、无破坏、易观察及低成本等优点,为数字化模拟现代光学实验提供了一种极好的手段. 本文探讨利用 MATLAB 软件实

34、现对任意形状衍射屏的夫琅禾费衍射实验的计算机仿真.夫琅禾费衍射实验装置如图 由基础光学可知,任意衍射屏的夫琅禾费衍射可借助两个透镜来实现. 如图1 所示,位于透镜 L1 物方焦平面上的点源 S 所发出的单色球面光波经 L1 变换为一束平面光波,照射在衍射屏 AB 上. 按照平面波理论,衍射屏开口处的波前向各个方向发出次波,方向彼此相同的衍射次波经透镜 L2 会聚到其像方焦平面的同一点上. 满足相长干涉条件时,该点为亮点;满足相消干涉条件时, 该点为P暗点. 所有亮点和暗点的集合构成了该衍射屏的夫琅禾费衍射图样. 其次,从傅里叶光学角度, 任意衍射屏在单位振幅的单色平面波垂直照射下, 其夫琅禾费

35、衍射光场复振幅即衍射屏透射系数的傅里叶变换, 而衍射图样实际上就是衍射屏的空间频谱强度分布.因此,可以用两种方法实现夫琅禾费衍射实验的仿真:1) 直接计算法. 通过推导给定衍射屏的夫琅禾费衍射图样强度分布公式, 得到观察屏上强度分布与位置的关系, 然后利用绘图函数将其光强度分布曲线和衍射图样绘出;2) 傅里叶变换法. 将衍射屏作为输入图像,经过二维傅里叶变换运算,得到衍西安毕业设计（论文）17射屏的频谱分布,即衍射图样. 前者适合于一些形状简单,且可以直接通过数学推导得到衍射光场强度分布公式的孔径(如单缝、双缝、矩形孔、多缝、光栅等) ;后者适合于任何形状的孔径. 如图所示为采用 MATLAB

36、 软件编写夫琅禾费衍射仿真实验程序的计算机流程图仿真方法 直接计算法3.1 夫琅禾费单缝衍射仿真夫琅禾费单缝衍射仿真一束单色平行光通过宽度可调的狭缝，射到其后的光屏上. 当缝宽足够小时，光屏上形成一系列亮暗相间的条纹，这是由于从同一个波前上发出的子波产生干涉的结果. 当光源到衍射屏的距离和光屏到衍射屏的距离都是无穷大时，即满足远场条件时，我们称这种衍射为夫琅禾费衍射. 所以夫琅禾费衍射中入射光和衍射光都是平行光. 为了模拟单缝衍射现象，我们把单缝看成一排等间隔光源，共NP个光源分布在 AB区间内，离 A 点间距为 yp，则屏幕上任一点S处的光强为NP个光源照射结果的合成. 如图 2.27 所示

37、，子波射线与入射方向的夹角称为衍射角，时，子0波射线通过透镜后，必汇聚到O点，这个亮条纹对应的光强称为主极大. NP个光源在其他方向的射线到达S点的光程差，应等于它们到达平面AC的光程差，即，其中sinyp sinysD为 S 点的纵坐标，则与A点光源位相差为ys西安毕业设计（论文）18 Y s A B C O 图 2.27 单缝衍射的模拟实验设单缝上NP个光源的振幅都为 1，在x，y轴上的分量各为，合振幅cossin，的平方为：. 又光强正比于振幅的平方，所以相对于O点22COSaCOSa 主极大光强也为程序模拟了单缝衍射现象，这里取波长 =500nm，缝宽 a=1mm，透镜焦距D=1m，运

38、行结果如图所示. 西安毕业设计（论文）19 单缝衍射仿真图分析图 2.28 中的衍射条纹，我们可以看出所有亮暗条纹都平行于单缝，O点光强为最大，这都和理论推导结果相一致3.2 夫琅禾费多缝衍射仿真夫琅禾费多缝衍射仿真 设每个狭缝的宽度均为 a ,相邻两狭缝间不透明部分的宽度为 b , 则缝间距(光栅常量) 为 d=a+b. 同样取 为衍射角, f 为透镜 L2 的焦距,衍射屏上透光的总狭缝数为 N ,则当平面光波垂直照射衍射屏时,沿 方向的衍射光波在P 处的合振动的相对强度为 式中, 称为单缝衍射因子,称为缝间的1sinaa sindv 22/sinuu干涉因子. 同样,若到的距离为 ,则由此

39、可以得到相对衍射光强P0P x度与之间的关系. x西安毕业设计（论文）20 3.3 夫琅禾费矩孔衍射仿真夫琅禾费矩孔衍射仿真假设矩形孔沿 x,y 方向的边长分别为 a,b,衍射光波的方向用二维衍射角和1西安毕业设计（论文）22来表示,则衍射光波在透镜 L2 的后焦平面上会聚点(x,y) 点的合振动20P的相对强度为式中 上式表明, 矩形孔衍射的相对强度是两/sin,/sindaa 0/PIPI个单缝衍射因子的乘积.从式(1)(2)(3)可以看出,确定了仿真程序中的可调参数为入射光波长 ，衍射屏的缝宽 ab,透镜的焦距 f 后,只要求出接收屏上每一点的相对衍射光强度值,就可以绘出衍射图样. 矩孔

40、衍射仿真图西安毕业设计（论文）224 菲涅尔衍射仿真234 菲涅尔衍射仿真菲涅尔衍射仿真当衍射物的尺寸比光波长大得多时，标量衍射理论是有效的。在光学系统设计时，光信息处理和传输等众多领域，标量衍射理论有着重要的应用。然而，基于惠更斯-菲涅尔原理的衍射积分的计算较于困难，为此需要对衍射积分进行近似处理并采用数值计算方法。当所研究的衍射场局限在旁轴区域时，菲涅尔近似在大多数情况下的衍射场数值计算问题已经进行了大量的研究，但未见有高性能的仿真算法的介绍。以 Matlab 为计算平台，以菲涅尔衍射积分为基础，采用子波叠加概念，针对方形孔径衍射设计了一种高性能的仿真算法，并给出了相应的程序和仿真结果。依

41、据数值模拟结果，考察了菲涅耳衍射场的光强分布对菲涅耳数的依赖关系,以及菲涅耳衍射过渡到夫琅禾费衍射的问题.文中给出的方孔菲涅耳衍射的仿真算法与采用菲涅耳正弦余弦积分的算法比较,在相同的参数和相同的高密度采样(例如:输入面做 3000 x3000 采样) 情况下,计算时间要少几倍到几十倍. 并且稍作修改即可用于具有任何透射率函数的矩形孔衍射计算和分析,只要该透射率函数可表示为分离变量形式. 运用文中所给出的算法,通过对大量的输出结果的分析,得出在菲涅耳近似下衍射光场的光强分布仅由菲涅耳数决定的结论. 并且,当菲涅耳数小于 0.08 时,菲涅耳衍射过渡到夫琅禾费衍射。4.14.1 菲涅尔方孔衍射仿

42、真菲涅尔方孔衍射仿真将方形孔径在平行于长度和高度方向分割成 N N 个微小单元，当 N 足够大时，每个单元可视为一个次级点源。所有点源在观察屏上点合成复oyxP,0振幅，可由式将积分求和得到选用 Matlab 为计算平台，在观察屏上取适当大小的正方形区域，并进行 MM 采样，采样点阵的坐标用二维数组存储；用二维数组存储输00YX 和yxII 和出面上各采样点对应的经由式计算，每次循环均采用数组运算，一次算出一个源点对所有目标点的贡献。即4 菲涅尔衍射仿真25 22/2022NNtzxxxixtteNLIzyYyeNLINNttjx/222202上式形式上完全相同的装置，从而的转置。因此只需计算

43、00XY 是XII 必是y一项，这可以减少一半的计算量。于是TXIIIx22z1本算法的优点是充分利用了 Matlab 高效的数组运算能力，在高密度采样下运行效率可大大提高输出面上的采样密度仅影响结果的平滑度。MATLAB 程序 方孔衍射图取各种不同的值时，输出的衍射仿真图和光强沿 x 轴分布图（横轴以孔径半宽度 L 为单位）输入面（孔径）上采样密度通常应大于输出面（观察屏）上的采样密度，且 z 值越小输入面上的采样密度应越大。作为比较另外采用菲涅尔正弦积分法设计了仿真算法，结果发现两者输出完全相同，但后者的时间、空间代价均大于前者。方孔衍射的菲涅尔数的定义为由数值计算的结果可以看到，当菲涅尔

44、数等于 0.3 时，采用夫琅禾费衍射公式和采用菲涅尔衍射公式计算出的衍射图样已经基本相同，光强分布的差别主要在极大值和极小值的取值，并且当菲涅尔数进一步减少时这种差别单调地也减少，最终菲涅尔衍射过渡到夫琅禾费衍射。图是菲涅尔数为 0.3 时的菲涅尔光强分布和夫琅禾费衍射光强分布，在中心点处最大，即。当 0maxII菲涅尔数等于 0.08 时，中心处最大光强的相偏差略大于 1%，因此可以认为菲涅尔数小于 0.8 时，菲涅尔衍射过渡到夫琅禾费衍射。也就是说，在试验中若采用夫琅禾费近似计算，欲确保精度，试验装置的菲涅尔数应小于 0.08。文中4 菲涅尔衍射仿真25给出的方孔菲涅尔衍射的仿真算法与采用

45、菲涅尔4 菲涅尔衍射仿真26正弦余弦的算法比较，在相同的参数相同的高密度采样（例如:输入面做3000 3000 采样）情况下，计算时间要少到几倍懂啊几十倍，并且稍作修改即可用于具有任何投射率函数的矩形孔衍射计算和分析，只要该投射率函数的可表示为分离变量形式。运用文中所给出的算法，通过大量的输出结果的分析，得出在菲涅尔近似下衍射光强的光强分布仅由菲涅尔数决定的结论。并且，当菲涅尔数小于 0.08 时，菲涅尔衍射过渡到夫琅禾费衍射。4.2 菲涅耳单缝衍射仿真菲涅耳单缝衍射仿真 (l)菲涅耳单缝衍射的光强分布 如图 6-2-6 所示，弧长 As 与单缝缝宽 d的关系为 rRRrds2 其中，R 为光

46、源与单缝的距离,r 为单缝与接收器的距离。在菲涅耳单缝衍射中，特定缝宽在屏幕上不同位置处将得到不同的光强。利用这种方法，借助MATLAB 编程可计算出菲涅耳单缝衍射在接收屏上沿垂直于缝方向的光强分布。用 MATLAB 计算菲涅耳单缝衍射的光强分布并进行彷真。解题分析：一束单色平行光，通过宽度可调狭缝，射到其后的接受屏上。接收屏上任一点尸的光强为 4 菲涅尔衍射仿真28 02 02 2sin,2cosdvvvdvvvvvV 是 P 点对应靠扭曲线上的弧长。对于平面波 （d 为缝宽）022rdv 在计算时，如果我们保持狭缝的位置固定，而计算观察平面的所有各点上的扰动，那么就起来将非常的麻烦，现在用

47、一种近似计算方法，将 S-O-P 直线（O 为狭缝 x 轴方向中心处）固定。这样，当在狭缝上下通过小的位移时，使原来 O 点相对于一个新的数值，这些新数值代入中又对应新的一点 P21,vvpI的光强值，从而可求出整个接收屏面的光强值。这种近似方法中的误差是可以忽略的，只要移动距离比光屏到狭缝的距离小得多的话。改方法更加适合平面波的情况。 程序运行上述程序，可分别取 d 为 1e-3，2e-3，3e-3 和 3e-3。在取不同的 d 值时，可适当调节最大光强和最大灰度的对比度，以得到衍射条纹图的最佳效果 单缝衍射图4 菲涅尔衍射仿真285 交互式 GUI 界面295 交互式交互式 GUI 界面界

48、面GUI 是 Graphical User Interface 的简称，即图形用户界面，通常人机交互图形化用户界面经常读做 “goo-ee” ，准确来说 GUI 就是屏幕产品的视觉体验和互动操作部分。 GUI 的广泛应用是当今计算机发展的重大成就之一，他极大地方便了非常专业用户的使用人们从此不再需记硬背大量的命令，取而代之的是可以通过窗口、菜单、快捷方式来方便地进行操作。而嵌入式GUI 具有几个方面的基本要求：轻型、占用资源少、高性能、高可靠性、便于移植、可配置等特点。本课题的 GUI 设计流程：1.创建窗口控件布局并美化 2.编写回调函数 3.编写算法核心程序 4. 加载算法核心程序到回调函数中 5.调试并完善 衍射仿真的 GUI 界面6 总结306 总结总结随着计算机运算能力的飞速发展计算机仿真技术作为虚拟实验