单缝衍射的matlab分析教程

单缝衍射的MATLAB分析学院：精密仪器与光电子工程学院专业：生物医学工程班级：1班姓名：单缝衍射的MATLAB分析摘要:在光的衍射概述和发展历史的基础上，说明了单缝衍射的图样特点，介绍了夫琅禾费衍射和菲涅耳衍射，几种实现夫琅禾费衍射的方法和原理及光强、条纹分布特点。

并利用衍射公式的近似对基尔霍夫衍射公式进行了推导，从理论上得出了夫琅禾费单缝衍射的光强公式，利用Matlab软件进行了光强分布的图样仿真，并用实验采集到的图样对理论和仿真的结论进行了验证，计算结果与实验结果得到了很好的吻合。

关键字：单缝衍射夫琅禾费单缝衍射光强分布条纹分布一、光的衍射概述1.光的衍射现象物理光学中，光的衍射现象是指光波在空间传播遇到障碍时，其传播方向会偏离直线传播，弯入到障碍物的几何阴影中，并呈现光强的不均匀分布的现象。

通常将观察屏上的不均匀的光强分布称为衍射图样。

光的衍射是光的波动性的主要标志之一。

光波遇到障碍物以后会或多或少地偏离几何光学传播定律的现象。

几何光学表明，光在均匀媒质中按直线定律传播，光在两种媒质的分界面按反射定律和折射定律传播。

但是，光是一种电磁波，当一束光通过有孔的屏障以后，其强度可以波及到按直线传播定律所划定的几何阴影区内，也使得几何照明区内出现某些暗斑或暗纹。

1.1衍射现象的基本问题1.已知照明光场和衍射屏的特征，求屏幕上衍射光场的分布；2.已知衍射屏及屏幕上衍射光场的发布，去探索照明光场的某些特性；3.已知照明光场及屏幕上所需的衍射光场发布，设计、计算衍射屏的结构和制造衍射光学元件。

1.2衍射现象的分类根据光源、衍射物（衍射屏）和衍射场（观察屏）三者之间的位置确定1．夫琅和费衍射：(远场衍射)光源和衍射场都在衍射物无限远处的衍射。

2．菲涅耳衍射：（近场衍射）光源和衍射场或二者之一到衍射物的距离比较小时的衍射。

1.3衍射现象及单缝衍射图样让一个足够亮的点光源S发出的光透过一个圆孔∑，照射到屏幕K上，并且逐渐改变圆孔的大小，就会发现：当圆孔足够大时，在屏幕上看到一个均匀光斑，光斑的大小就是圆孔的几何投影，随着圆孔逐渐减小，起初光斑也相应的变小，而后光斑开始模糊，并且在圆斑外面产生若干围绕圆斑的同心圆环，当使用单色光源时，是一组明暗相见的同心环带，当使用白色光源时，是一组色彩相间的彩色环带；此后再使圆孔变小，光斑及圆环不跟着变小，反而会增大起来。

单色红光衍射图样白光衍射图样1.4衍射的应用光的衍射决定光学仪器的分辨本领；气体或液体中的大量悬浮粒子对光的散射，衍射也起重要的作用。

衍射应用大致可以概括为以下四个方面：1、光谱分析：如衍射光栅光谱仪。

2、结构分析:衍射图样对精细结构有一种相当敏感的“放大”作用,故而利用图样分析结构。

如X射线结构学。

3、成像:在相干光成像系统中,引进两次衍射成像概念，由此发展成为空间滤波技术和光学信息处理。

如光瞳衍射导出成像仪器的分辨本领。

4、波阵面再现: 一种全新的两步无透镜成像法，也称为波阵面再现术，这是全息术原理中的重要一步。

2. 衍射现象的发展过程大约1818年前，一直没有人注意到有可能根据波动理论说明衍射效应。

1818年，菲涅耳的著作问世。

他在论文中证明，应用慧更斯作图法，结合干涉原理，能够解释衍射现象。

菲涅耳的分析后来由基而霍夫在1882年给出了完善的数学描述。

衍射问题是光学中遇到的最为困难的问题之一。

在衍射理论中,那种在某种意义上可以认为是严格的解，是很少有的。

直至1896年，才由索末菲给出了第一个解。

他在一篇重要论文中讨论了一个完全导电的半无限平面屏对平面波的衍射。

此后，对少数其它衍射问题（二维）也求得了严格解。

由于在数学上的困难，在大多数有实际意义情况下，必须采用近似方法。

这些方法中惠更斯－菲涅耳理论是最富成效的，它适用于处理光学仪器中所遇到的大多数光学衍射问题。

二、单缝衍射原理1.惠更斯—菲涅耳原理最早成功地用波动理论解释衍射现象的是菲涅耳，他将惠更斯原理用光的干涉理论加以补充，并予以发展。

惠更斯原理是描述波动传播过程的一个重要原理，其主要内容是：如图2-1所示的波源S，在某一时刻所产生波的波阵面为∑，则∑面上的每一点都可以看作是一个次波源，它们发出球面次波，其后某一时刻的波阵面∑＇即是该时刻这些球面次波的包络面，波阵面的法线面的法线方向就是该波的传播方向。

惠更斯原理能够很好地解释光的直线传播，光的反射和折射方向，但不能说明衍射过程及其强度分布。

菲涅耳在研究了光的干涉现象后，考虑到次波来自于同一光源，应该相干，因而波阵面∑＇上每一点的光振动应该是在光源和该点之间任一波面上的各点发出的次波场叠加的结果。

这就是惠更斯—菲涅耳原理。

利用惠更斯—菲涅耳原理可以解释衍射现象：在任意给定的时刻，任一波面上的点都起着次波波源的作用，它们各自发出球面次波，障碍物以外任意点上的光强分布，即是没有被阻挡的各个次波源发出的次波在该点相干叠加的结果。

根据惠更斯—菲涅耳原理，图2-2所示的一个单色光源S 对于空间任意点P 的作用，可以看作是S 和P 之间任一波面∑上各点发出的次波在P 点相干叠加的结果。

假设波面Σ上任意点Q 的光场复振幅为)(~Q E ，在Q 点取一个面元d δ，则d δ面元上的次波源对P 点光场的贡献为：式中，C 是比例系数；rQP =，()K θ称为倾斜因子，它是与元波面法线和P Q 的夹角θ（称为衍射角）有关的量，按照菲涅耳的假设：当0θ=时，K 有最大值；随着θ的增大，K 迅速减小；当/2θπ≥时K ＝0。

因此，途中波面∑上只有ZZ '范围内的部分对P 点光振动有贡献。

所以P 点的光场复振幅为：()()()ikre E P C E Q K d r θσ∑=⎰⎰这就是惠更斯－菲涅耳原理的数学表达式，称为惠更斯－菲涅耳公式。

当S 是点光源时，Q 点的光场复振幅为：()ikRA E Q e R=式中，R 是光源到Q 点的距离。

在这种情况下，)(~Q E 可以从积分号中提出来，但是由于()K θ的具体形式未知，不可能由(2-1)式确切地确定()E P 值。

因此，从理论上来讲，这个原理是不够完善的。

2.实现夫琅禾费衍射的几种方法无论是在实验室中或者别的什么地方,都不可能将光源和衍射场放在无限远，实际接收夫琅禾费衍射的装置有以下四种：1.焦面接收装置（以单缝衍射为例，下同）把点光源S 放在凸透镜1L 的前焦平面上，在凸透镜2L 的后焦平面上接收衍射场，见图2-5。

图 2-5 焦面接收装置图 2-6 远场接收装置2、 远场接收装置θDZP θ0P∑ θSD∑P θP fL 1L 2当满足远场条件时，狭缝前后也可以不用透镜，而直接获得夫琅禾费衍射图样。

远场条件是：① 光源离狭缝很远，即λρ2>>R ，其中，R 是光源到狭缝的距离，ρ为狭缝宽度的一半；② 接收场距狭缝足够远，即λρ2>>z ，其中，z 为衍射场距狭缝的距离。

观察点P在λρ2>>z 的条件下，只要求其满足傍轴条件即可，而这一般都是满足的。

图2-6为远场接收光路，假设一束平行光垂直入射到狭缝上。

3、象面接收装置（一）衍射屏处于透镜的后方，如图2-7所示。

S 在光轴上，∑代表点光源的象面，S '为S 的象点。

理论上已经证明了∑面上呈现的图样为夫琅禾费衍射图样，即屏上任一点P θ的复振幅与角度θ的函数关系符合夫琅禾费衍射的积分形式。

图 2-7 象面接收装置（一）图 2-8 象面接收装置（二）4、象面接收装置（二）衍射屏处于透镜的前方，如图2-8所示。

P θ'点是场点P θ的共轭点，S 也在光轴上。

如果光路逆转自右向左，S '变为点光源，衍射屏便处于透镜的后方了，'∑面上的衍射图样就θDZP θS '∑ SθP θS '∑ θ'S'∑ P θ'Z 'u vxx 'D L同象面接收装置（一）∑面上的情况，z '相应地取代z ，所以实际呈现在图2-8的∑面上的衍射图样可由物面上设想的共轭衍射图样导出，二者为物象关系。

3.夫琅禾费衍射光强强度的计算现在我们用惠更新-菲涅耳原理来解释上述现象。

如图2-13所示。

为了清楚起见，图中狭缝的宽度'BB 已经放大。

平行光束垂直于缝的平面入射时，波面和缝平面重合（垂直于图面）。

将缝的面积分为一组平行于缝长的窄带，从每一条这样的窄带发出次波。

其振幅正比于窄带的宽度dx ，设光波的初位相为零，b 为缝'BB 的宽度，b A 0，而宽度dx 的窄条上次波的振幅为b dx A 0，则狭缝处各窄带所发次波的振动可用下式表示：00cos A dxdE t b ω=这些次波都可认为是球面波，各自向前传播。

现在，首先对其中沿图面与原入射方向成θ角（称为衍射角）的方向传播的所有各次波进行研究。

在入射光束的平面波面BB’上各次波的位相都相等，通过透镜2L 后在焦平面FF 上的同一点P 处叠加。

要计算P 点的合振幅，必须考虑到各次波的位相关系，这取决于由各窄带到P 点的光程如何。

现在作平面BD 垂直于衍射方向'B D ，根据BD 面上各点的位相分布情况即可决定在P 点相遇的各次波的位相关系。

我们知道，从平面BD 上各点沿衍射方向通过透镜而达到P 点的光程都相等。

这就只要算出从平面'BB 到平面BD 的各平行直线段之间的光程差就可以了。

MN 为沿着衍射角θ进行的任一条路程，令BM =x ，则sin MN x θ=，这就是从M 和从B 两点所发次波沿平行于MN 方向到达平面BD 时的光程差。

得BD 面上N 点的光振动的表达式为02cos(sin )A dx dE x t b πθωλ=-或 2(sin )0i x t A dx dE eb πθωλ-=其复振幅为： 2sin 0i x A dx dE eb πθλ=为简化计算起见，上式中假设各次波到达P 点时有相同的振幅（不考虑振幅与光程与反比的关系以及华侨因数）。

根据惠更斯—菲涅耳原理，将上式对整个缝宽（从x=0到x=b ）积分。

最后可得沿着衍射角θ方向传播的所有次波在观察点P 叠加起来的合振幅：sin(sin )sin P bA A bπθλπθλ=令(sin )/u b πθλ=，通常称(sin )/u u 为u 的sin c 函数，并写成sin cu ，故P 点的光强为20sin P I I c u =4.夫琅禾费衍射图样的光强分布当光屏放置在透镜L 2的焦平面上时，屏上出现衍射花样，光强的分布可由上式决定。

不同的衍射角θ对应于光屏上不同的观察点。

首先来决定衍射花样中光强最大值和最小值的位置。

即求出满足光强的一阶导数为零的那些点：223sin 2sin (cos sin )()0d u u u u u du u u -==由此得 sin 0,u u tgu == 分别解以上两式，可得出所有的极值点。