单缝衍射的matlab分析报告.docx

单缝衍射的MATLAB 分析学院：精密仪器与光电子工程学院专业：生物医学工程班级:1班：单缝衍射的MATLAB 分析摘要:在光的衍射概述和发展历史的基础上，说明了单缝衍射的图样 特点，介绍了夫琅禾费衍射和菲涅耳衍射，几种实现夫琅禾费衍射的 方法和原理及光强、条纹分布特点。并利用衍射公式的近似对基尔霍 夫衍射公式进行了推导，从理论上得出了夫琅禾费单缝衍射的光强公 式，利用Matlab软件进行了光强分布的图样仿真，并用实验采集到的图样对理论和仿真的结论进行了验证，计算结果与实验结果得到了很好的吻合。关键字：单缝衍射 夫琅禾费单缝衍射光强分布 条纹分布一、光的衍射概述1. 光的衍射现象物理光学中，光的衍射现象是指光波在空间传播遇到障碍时，其传播方向会偏离直线传播，弯入到障碍物的几何阴影中，并呈现光强的不均匀分布的现象。通常将观察屏上的不均匀的光强分布称为衍射图样。光的衍射是光的波动性的主要标志之一。光波遇到障碍物以后会或多或少地偏离几何光学传播定律的现象。几何光学表明，光在均匀媒质中按直线定律传播，光在两种媒质的分界面按反射定律和折射定律传播。但是，光是一种电磁波，当一束光通过有孔的屏障以后，其强度可以波及到按直线传播定律所划定的几何阴影区，也使得几何照明区出现某些暗斑或暗纹。1.1衍射现象的基本问题1. 已知照明光场和衍射屏的特征，求屏幕上衍射光场的分布;2. 已知衍射屏及屏幕上衍射光场的发布，去探索照明光场的某些特性；3. 已知照明光场及屏幕上所需的衍射光场发布，设计、计算衍射屏的结构和制造衍射 光学元件。1.2衍射现象的分类根据光源、衍射物（衍射屏）和衍射场（观察屏）三者之间的位置确定1. 夫琅和费衍射：（远场衍射）光源和衍射场都在衍射物无限远处的衍射。2.菲涅耳衍射：（近场衍射）光源和衍射场或二者之一到衍射物的距离比较小时的衍射。1.3衍射现象及单缝衍射图样让一个足够亮的点光源S发出的光透过一个圆孔2,照射到屏幕K上，并且逐渐改 变圆孔的大小，就会发现：当圆孔足够大时，在屏幕上看到一个均匀光斑，光斑的大小就是 圆孔的几何投影，随着圆孔逐渐减小，起初光斑也相应的变小，而后光斑开始模糊，并且在 圆斑外面产生若干围绕圆斑的同心圆环，当使用单色光源时，是一组明暗相见的同心环带， 当使用白色光源时，是一组色彩相间的彩色环带；此后再使圆孔变小，光斑及圆环不跟着变 小，反而会增大起来。光源与照衍射屏衍射空间接收幕屏耳（工5） e 1） e（x.y ）白光衍射图样单色红光衍射样1.4衍射的应用光的衍射决定光学仪器的分辨本领；气体或液体中的大量悬浮粒子对光的散射，衍射也起重要的作用。衍射应用大致可以概括为以下四个方面:1、光谱分析：如衍射光栅光谱仪。2、结构分析:衍射图样对精细结构有一种相当敏感的"放大"作用故而利用图样分析结构。如X射线结构学。3、成像:在相干光成像系统中，引进两次衍射成像概念，由此发展成为空间滤波技术和光学信息处理。如光瞳衍射导出成像仪器的分辨本领。4、波阵面再现：一种全新的两步无透镜成像法，也称为波阵面再现术，这是全息术原理中的重要一步。2.衍射现象的发展过程大约1818年前，一直没有人注意到有可能根据波动理论说明衍射效应。1818年，菲 涅耳的著作问世。他在论文中证明，应用慧更斯作图法，结合干涉原理，能够解释衍射现象。菲涅耳的分析后来由基而霍夫在1882年给出了完善的数学描述。衍射问题是光学中遇到的最为困难的问题之一。在衍射理论中那种在某种意义上可以 认为是严格的解，是很少有的。直至1896年，才由索末菲给出了第一个解。他在一篇重 要论文中讨论了一个完全导电的半无限平面屏对平面波的衍射。此后，对少数其它衍射问题（二维）也求得了严格解。由于在数学上的困难，在大多数有实际意义情况下，必须采用近似方法。这些方法中惠 更斯-菲涅耳理论是最富成效的，它适用于处理光学仪器中所遇到的大多数光学衍射问题。二, 单缝衍射原理1. 惠更斯一菲涅耳原理最早成功地用波动理论解释衍射现象的是菲涅耳，他将惠更斯原理用光的干涉理论加以 补充，并予以发展。图2-2单色点光源S对P点的光作用惠更斯原理是描述波动传播过程的一个重要原理其主要容是女口图2-1所示的波源S, 在某一时刻所产生波的波阵面为2,则2面上的每一点都可以看作是一个次波源，它们发出 球面次波，其后某一时刻的波阵面2即是该时刻这些球面次波的包络面，波阵面的法线面 的法线方向就是该波的传播方向。惠更斯原理能够很好地解释光的直线传播，光的反射和折 射方向，但不能说明衍射过程及其强度分布。菲涅耳在研究了光的干涉现象后，考虑到次波来自于同一光源，应该相干，因而波阵面 2上每一点的光振动应该是在光源和该点之间任一波面上的各点发出的次波场叠加的结 果。这就是惠更斯一菲涅耳原理。利用惠更斯一菲涅耳原理可以解释衍射现象：在任意给定的时刻，任一波面上的点都 起着次波波源的作用，它们各自发出球面次波，障碍物以外任意点上的光强分布，即是没有 被阻挡的各个次波源发出的次波在该点相干叠加的结果。根据惠更斯一菲涅耳原理，图2-2所示的一个单色光源S对于空间任意点P的作用， 可以看作是S和P之间任一波面2上各点发出的次波在P点相干叠加的结果。假设波面£上 任意点Q的光场复振幅为E（Q），在Q点取一个面元d8,则d8面元上的次波源对P点光 场的贡献为：式中，C是比例系数；r = QP，K（0 ）称为倾斜因子，它是与元波面法线和QP的夹 角0 （称为衍射角）有关的量，按照菲涅耳的假设：当0 = 0时，K有最大值；随着0的增 大，K迅速减小；当0 N兀/2时K = 0o因此，途中波面上只有ZZ'围的部分对P点光振 动有贡献。所以P点的光场复振幅为：E（P）= CE（Q ）竺 K G）姑rz这就是惠更斯-菲涅耳原理的数学表达式，称为惠更斯-菲涅耳公式。当S是点光源时，Q点的光场复振幅为：E（Q ）= AeikR式中，R是光源到Q点的距离。在这种情况下，E（Q）可以从积分号中提出来，但是由于 K（0 ）的具体形式未知，不可能由（2-1）式确切地确定E（P）值。因此，从理论上来讲，这 个原理是不够完善的。2. 实现夫琅禾费衍射的几种方法无论是在实验室中或者别的什么地方都不可能将光源和衍射场放在无限远，实际接收 夫琅禾费衍射的装置有以下四种：1.焦面接收装置（以单缝衍射为例，下同）把点光源S放在凸透镜L1的前焦平面上，在凸透镜L2的后焦平面上接收衍射场，见图 2-5。£PeP0图2-5焦面接收装置D 2件是：光源离狭缝很远，即R >>，其中，R是光源到狭缝的距离，P为狭缝宽度的人P 2 一半；接收场距狭缝足够远，即z >>，其中，z为衍射场距狭缝的距离。观察点P人P2在Z >>的条件下，只要求其满足傍轴条件即可，而这一般都是满足的。图2-6为远场人接收光路，假设一束平行光垂直入射到狭缝上。3、象面接收装置（一）衍射屏处于透镜的后方，如图2-7所示。S在光轴上，£代表点光源的象面，S'为S的象点。理论上已经证明了£面上呈现的图样为夫琅禾费衍射图样，即屏上任一点P的复 e振幅与角度e的函数关系符合夫琅禾费衍射的积分形式。图2-7象面接收装置（一）衍射屏处于透镜的前方，如图2-8所示。p点是场点/p的共轭点，S也在光轴上。如果光路逆转自右向左，S'变为点光源，衍射屏便处于透镜的后方了，£'面上的衍射图样就 同象面接收装置（一）£面上的情况，z'相应地取代z，所以实际呈现在图2-8的£面上 的衍射图样可由物面上设想的共轭衍射图样导出，二者为物象关系。3. 夫琅禾费衍射光强强度的计算现在我们用惠更新-菲涅耳原理来解释上述现象。如图2-13所示。为了清楚起见，图 中狭缝的宽度BB'已经放大。平行光束垂直于缝的平面入射时，波面和缝平面重合（垂直于 图面）。将缝的面积分为一组平行于缝长的窄带，从每一条这样的窄带发出次波。其振幅正 比于窄带的宽度dx，设光波的初位相为零，b为缝BB的宽度，AJb，而宽度dx的窄条 上次波的振幅为A0dx/b，则狭缝处各窄带所发次波的振动可用下式表示：dE = A# cos s t这些次波都可认为是球面波，各自向前传播。现在，首先对其中沿图面与原入射方向成 。角（称为衍射角）的方向传播的所有各次波进行研究。在入射光束的平面波面日日上各 次波的位相都相等，通过透镜乙2后在焦平面FF上的同一点P处叠加。要计算P点的合振 幅，必须考虑到各次波的位相关系，这取决于由各窄带到P点的光程如何。现在作平面BD 垂直于衍射方向BD，根据BD面上各点的位相分布情况即可决定在P点相遇的各次波的 位相关系。我们知道，从平面BD上各点沿衍射方向通过透镜而达到P点的光程都相等。这 就只要算出从平面BB到平面BD的各平行直线段之间的光程差就可以了。M 为沿着衍 射角。进行的任一条路程，令丽=x,则MN = x sin 0，这就是从M和从B两点所发次波 沿平行于MN方向到达平面BD时的光程差。得BD面上N点的光振动的表达式为dE = ”0赤 cos(半xsin0 -st) b 人其复振幅为：为简化计算起见，上式中假设各次波到达P点时有相同的振幅（不考虑振幅与光程与 反比的关系以及华侨因数）。根据惠更斯一菲涅耳原理，将上式对整个缝宽（从x=0到x=b） 积分。最后可得沿着衍射角。方向传播的所有次波在观察点P叠加起来的合振幅:sin（誓 sin 6） "oFT -sin 6入令u =（兀bsin6）/人，通常称（sinu）/u为u的sinc函数，并写成sincu，故P点的光强为I = I sin c2u4.夫琅禾费衍射图样的光强分布当光屏放置在透镜L2的焦平面上时，屏上出现衍射花样，光强的分布可由上式决定。不同的衍射角e对应于光屏上不同的观察点。首先来决定衍射花样中光强最大值和最小值的 位置。即求出满足光强的一阶导数为零的那些点：d sin2 u2sin u （u cos u 一 sin u）（）=0du u 2u 3由此得sin u = 0, u = tgu分别解以上两式，可得出所有的极值点。1、单缝衍射中央最大值的位置：由sin u = 0 ,解得满足u° = «b sin6°）人=0的一些衍射方向，即I sin6 o = 0 （中间最 大值的位置），也就是在焦点P0处，I = A2，光强为最大。这里，叠加的各个次波位相差p00为零，所以振幅叠加互相加强。2、单缝衍射最小值的位置:由sin u = 0，解得满足"广2兀焙sin七人=2质的一些衍射方向，即 sin七=kb（k = ±1,±2,±3）（最小值位置）时，AP为零，屏上这些点是暗的。三、光强分布和条纹分布分析1. 基尔霍夫衍射公式的近似对于一些极简单的衍射问题，也因为被积函数形式复杂而得不到解析形式的积分结果。为此，必须根据实际条件进一步做近似处理。1.1傍轴近似在一般的光学系统中，对成像起主要作用的是那些与光学系统光轴夹角极小的傍轴光 线。对于傍轴光线，如图3-3开孔的线度和观察屏上的考察围都远小于开孔到观察屏的距 离。因此，下面的两个近似条件都成立：1、cosG,广）» 1,于是KG）q 1 ；2、r 总 z1。于是，（3-11）可以化简为E（P）=上jj E（Q）e，krd81 X指数中的r未用z1代替，这是因为指数中r所影响的是次波场的相位，r的微小变化都 会引起相位的很大变化。图3-3孑L径的衍射1.2距离近似根据衍射现象，离衍射孔不同距离处，衍射图样是不同的。一种是菲涅耳衍射或近场 衍射，指的是光源和接收屏与衍射屏的距离均为有限远，或其中之一是有限远的情形；另一 种是夫琅禾费衍射或远场衍射，指的是光源和接收屏与衍射屏的距离均为无限远的情形。1、菲涅耳近似 如图3-3所示，设Q = r，则由几何关系有-气)2 + (y-yi )2 = Z1 +i X - XG-X1)2 + (y -y .J2z 21G-X 1)2 + (y -y 1)2z 21因为在指数上的相位因子决定了函数的周期性，每当相位因子改变兀时，指数函数反2+号，这种变化是不可忽略的，相位因子中只有远小于兀的项才可忽略。于是，当乙大到满,max << 兀k 匕-气* +(y-y»81 云1时，上式第三项及以后的各项都可略去，即为X 2 + y 2xx + yy X 2 + y 2"2 ZZ2 Z这一近似称为菲涅耳近似，在这个区域观察到的衍射现象叫菲涅耳衍射。在菲涅耳近 似下，P点的光场复振幅为(x- x h +(y-y1 2 z22、夫琅禾费近似当观察屏离孔的距离很大，满足Q + y 2)k 11max << 兀2 z1时r = z +挡工-欢1 +当12ZZ这一近似称为夫琅禾费近似，在这个区域观察到的衍射现象叫夫琅禾费衍射。在夫琅禾费近似下，P点的光场复振幅为： ()好 1 ik认 .-ikT!,Ex, y =一e2z. JJ E(x , y )ez.dx dy人 z1 11 11s2. 夫琅禾费单缝衍射光强分布如果只考虑单色平行光垂直入射到开孔平面上的的夫琅禾费衍射，则通常采用图3-4 所示的夫琅禾费衍射装置。单色点光源放置在透镜的前焦平面上所产生的平行光垂直入射开孔矿由开孔的衍射，在透镜、的后焦平面上可以观察到开孔酒夫琅禾费衍射图样后焦平面上的各点的光场复 振幅由(3-19 )式给出。图3-4夫琅禾费衍射装置若开孔面上有均匀的光场分布，令E(x, y1 )= A =常数。又因为透镜贴近孔径，一 riA ik f+z1牝f。所以后焦平面上的光场复振幅可为E(x, y)= Cjj e-ik(xx1+yy1 >f dx dy , C =-e s对于夫琅禾费单缝衍射，水平缝宽为a,垂直缝宽为b,则b » a，沿y方向的衍射效应不明显，只在x方向有亮暗变化的衍射图样。由(3-20 )式，衍射屏上P点的光场复振幅E (p )=e V CLe J2 = 2Cf sin1 ikx - a、 ikx2sin=Ca 'kxa ' 一巨 kxa 2f-sin a=E 0 a式中E0 = Ca是观察屏中心点p0处光场复振幅。相应P点的光强度为"以J式中，I02, a=绑=KE ¥ & sin 9 M,9为衍射角。在衍射理论中，通常称 2 J人J(in 2 a. a 2)为单缝衍射因子，10为中央明纹中心处光强度，a为单缝边缘光线与中心光线 的相位差。根据(3-22 )式可的单缝衍射的光强分布特征:1、在中央P点,9 = 0，我们使a为一很小的趋于零的角，对a求极限，则主兰=1 , a，为中央主极大光强。2、3、在a = k兀时(k = 1>2),振幅A = A - sn以=0，光强I = 0为暗纹的光强。0 asinGk'2)当a = 3兀,2时，同理有112=0.04710为一级次极大光强。4、当a= 5k/2时，同理有12-sin 喝 2)2=0.01610为二级次极大光强。5、当a= 7官2时，同理有13 = A0 -2=0.00810为三级次极大光强。3. 夫琅禾费单缝衍射光强分布31光强分布的极值点当光屏放置在透镜L2的焦平面上时，屏上出现衍射花样，光强的分布可由(3-22 )式 决定，也就是说，衍射光场在P点的强度大小主要由因子a或sin2 a.'a2决定。不同的衍射角。对应于光屏上不同的观察点。1、当a =。或9 = 0时，该点光强度取最大值:I(P) = I(P) = Imax，称为主极大值；2、当a = ±馈或a sin 9 = 土瓜,(k = 1,2,3)时，该点强度取极小值:I (P) = 0 = I .;min3、相邻两个极小值之间存在一个极大值，由于因子sin2 a.a2 < 1，该极大值的强度总是小于主极大值，故称之为次极大值。求出满足光强的一阶导数为零的那些点:解得d (sin2a'da2sin a (a cos a - sin a) =0k a ；sin a = 0, a = tan a于是得到夫琅禾费单缝衍射的次极大值位置满足关系为：a = tan a。对于这一超越方程其根为a =±1.43 兀,±2.46 兀,±3.47兀对应的sin 0值为sin 0 = ±1.43一,±2.46 ,±3.47 a a a3.2条纹的角宽度和线宽度1.暗纹的角宽度暗纹出现的位置在a = 土兀，±2兀，土3兀，的地方，对应的sin0值为sin 0=±X ,± 王,± 全，aa aX任何两相邻暗纹间的衍射角的差值为履=±一，即暗纹是以以点为中心等间隔左右对称分 a0布的。2.亮纹的角宽度和线宽度以相邻暗纹的焦距离作为亮纹的角宽度。在傍轴条件下，即0很小时，暗纹的衍射角kX位置简化为0 r sin0 =±一，k =1,2,3。中央主极大明纹宽度由k = ±1的暗纹的衍射 a角所确定，则主极大亮纹角宽度和线宽度分别为A0 =竺，心=f-A0 =如(1)aa在两相邻暗纹间存在的次极明纹，其角宽度和线宽度分别为 A0Xf AxAO = 一 = , Ax = f - AO = = （2）a2a23.3条纹分布的影响因素从以上结论可以看出，夫琅禾费单缝衍射图样的强度随衍射角度按着函数关系sin 2 a. a 2变化；相邻暗条纹中心的角间距相等，因而所有次极大值亮纹的角宽度相等，但 主极大值亮纹的角宽度为次极大值的两倍;相邻次极大值亮纹中心不等间距，随着衍射级次 的增大，相邻次极大值亮纹中心的间距趋于恒定。由（1）和（2）式可以看出，条纹分布的影响因素有：1、狭缝的宽度对条纹分布的影响对于给定的波长，亮纹的宽度是与狭缝的宽度成反比的，在波前上对光束限制越大， 也就是缝宽越小的时候，衍射场越弥散，衍射图样铺开的越宽；反之，当缝宽很大，光束几 乎自由传播时，A。r 0，这表明衍射场基本上集中在沿直线传播的原方向上，在透镜焦面 上衍射图样收缩为几何光学像点。2、入射光的波长对条纹分布的影响在保持缝宽不变的条件下，波长越长，衍射效应越明显；波长越短，衍射效应越可忽 略；但是，实际中，由于波长的数量级比较小，所以对衍射条纹的影响不是非常明显。3、接收屏的距离对条纹分布的影响当其他条件不变的情况下，若接收屏沿轴向向前移动时，接收屏上的条纹间距变小， 衍射图样展开围减小；若沿轴向向后移动，接收屏上的条纹间距变大，衍射图样展开围也增 大。,单缝衍射的实验模拟根据理论的推导，利用衍射积分法对单缝夫琅禾费衍射进行Matlab仿真。衍射积分法 是由衍射积分算出接收屏上的光强分布然后根据该分布调制色彩作图从而得到衍射图案。1.夫琅禾费单缝衍射程序和夫琅禾费单缝衍射的仿真程序：波长500nm,缝宽1mm，屏距1m。clearlam=500e-9;a = 1e-3;D=1;ym=3*lam*D/a;ny=51;ys=linspace(-ym,ym,ny);np=51;yp=linspace(0,a,np);for i = 1:nysinphi=ys(i)/D;alpha=2*pi*yp*sinphi/lam;sumcos二sum(cos(alpha);sumsin二sum(sin(alpha);B(i,:) = (sumcos人2+sumsin人2)/np人2;endN=255;Br=(B/max(B)*N;subplot(1,2,1)image(ym,ys,Br);colormap(gray(N);subplot(1,2,2)plot(B,ys);2. 不同缝宽条纹光强分布的程序和图样当入射光波长为710.nm，狭缝到接收屏的距离为1m，狭缝的缝宽分别为0.5mm,0.1mm, 0.05mm，进行仿真： clear;%清除标量l=710e-9;%红光波长 d=0.5,0.1,0.05*1e-3;% 设置缝宽矩阵n = -1:0.0001:1;%自变量向量fai二n*pi/180;% 角度向量FAI,D = meshgrid(fai,d);%绘制三维网格，建立矩阵u=2*pi*D.*sin(FAI)/l;% 中间变量u(u =0)=eps;%若零改为小量i=(sin(u)./u).人2;% 光强figure%创建图形窗口plot(fai,i,'-','LineWidth',2) % 画曲线title('单缝夫琅禾费衍射不同缝宽的光强曲线,'FontSize',19,'Fontname'，黑体')%标题xlabel('衍射角 fai','FontSize',19,'Fontname',黑体') 标记横坐标ylabel(相对光强 i','FontSize',19,'Fontname'，黑体') 标记纵坐标text(0.012,0.93,'红线 representd=0.05');% 标注text(0.012,0.83,'绿线 representd=0.1');% 标注text(0.012,0.73,'蓝线 representd=0.5');% 标注grid on%绘制网格c=255;%颜色大小figure%创建图形窗口image(i*1000) % 画图colormap(gray(c)%形成线性灰度色图title('单缝夫琅禾费衍射不同缝宽下的条纹分布,'FontSize',19,'Fontname'，黑体')%标题 xlabel('至上而下：d=0.5 d=0.1 d=0.05','FontSize',19,'Fontname',黑体') 标注图样:从图中可以看出，当缝宽改变时，条纹宽度随着缝宽的增大而减小，反之亦然。可以得出结论：当只有缝宽改变，其他条件不改变的情况下，缝宽也与条纹宽度成反比。参考文献：1郁道银，谈恒英.工程光学M机械工业，2011年：物理光学部分2伟民.单缝衍射图样的计算机仿真与分析J.工科物理,1999,(01).3厉江帆，姜宗福，黄春佳,贺慧勇，朱江峰.夫琅禾费衍射公式的一般形式J.大学物理.2003,(11): 9-14凯华,钟锡华.光学M.大学,1982:65-66.5耿顺山.基于MATLAB的单色光单缝夫琅禾费衍射现象的模拟J.电脑学习.2008,(01):20-216春红.夫琅禾费单狭缝衍射条纹分析J.职业技术学院学报.2003,3(03):23-26.