单缝衍射的matlab分析.doc

单缝衍射的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点取一个面元dδ，则dδ面元上的次波源对P点光场的贡献为：

式中,C是比例系数；，称为倾斜因子,它是与元波面法线和的夹角(称为衍射角）有关的量，按照菲涅耳的假设：当时,K有最大值；随着的增大，K迅速减小;当时K＝0。因此，途中波面∑上只有范围内的部分对P点光振动有贡献。所以P点的光场复振幅为：

这就是惠更斯－菲涅耳原理的数学表达式，称为惠更斯－菲涅耳公式。
当S是点光源时，Q点的光场复振幅为：

式中，R是光源到Q点的距离。在这种情况下,可以从积分号中提出来，但是由于的具体形式未知，不可能由（2—1）式确切地确定值。因此,从理论上来讲，这个原理是不够完善的.

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

S
D

L1
L2

图 2-5 焦面接收装置

D
Z

图 2-6 远场接收装置
 2、 远场接收装置
当满足远场条件时，狭缝前后也可以不用透镜,而直接获得夫琅禾费衍射图样。远场条件是:① 光源离狭缝很远，即，其中，R是光源到狭缝的距离，为狭缝宽度的一半；② 接收场距狭缝足够远，即，其中，z为衍射场距狭缝的距离。观察点P在的条件下,只要求其满足傍轴条件即可，而这一般都是满足的。图2-6为远场接收光路,假设一束平行光垂直入射到狭缝上。
3、象面接收装置(一）
衍射屏处于透镜的后方,如图2—7所示。S在光轴上，代表点光源的象面,为S的象点。理论上已经证明了面上呈现的图样为夫琅禾费衍射图样,即屏上任一点的复振幅与角度的函数关系符合夫琅禾费衍射的积分形式.

D
Z

S

图 2-7 象面接收装置（一)

S

D
L

图 2—8 象面接收装置（二）
4、象面接收装置（二）
衍射屏处于透镜的前方，如图2-8所示。点是场点的共轭点，S也在光轴上。如果光路逆转自右向左，变为点光源，衍射屏便处于透镜的后方了，面上的衍射图样就同象面接收装置(一)面上的情况，相应地取代z,所以实际呈现在图2—8的面上的衍射图样可由物面上设想的共轭衍射图样导出,二者为物象关系。
3.夫琅禾费衍射光强强度的计算
现在我们用惠更新—菲涅耳原理来解释上述现象。如图2—13所示.为了清楚起见，图中狭缝的宽度已经放大。平行光束垂直于缝的平面入射时，波面和缝平面重合（垂直于图面）。将缝的面积分为一组平行于缝长的窄带，从每一条这样的窄带发出次波.其振幅正比于窄带的宽度，设光波的初位相为零,b为缝的宽度，，而宽度的窄条上次波的振幅为，则狭缝处各窄带所发次波的振动可用下式表示:

这些次波都可认为是球面波，各自向前传播。现在，首先对其中沿图面与原入射方向成θ角（称为衍射角）的方向传播的所有各次波进行研究.在入射光束的平面波面BB’上各次波的位相都相等，通过透镜后在焦平面FF上的同一点P处叠加。要计算P点的合振幅，必须考虑到各次波的位相关系,这取决于由各窄带到P点的光程如何。现在作平面BD垂直于衍射方向，根据BD面上各点的位相分布情况即可决定在P点相遇的各次波的位相关系。我们知道，从平面BD上各点沿衍射方向通过透镜而达到P点的光程都相等。这就只要算出从平面到平面BD的各平行直线段之间的光程差就可以了.为沿着衍射角θ进行的任一条路程，令=x，则，这就是从M和从B两点所发次波沿平行于MN方向到达平面BD时的光程差。得BD面上N点的光振动的表达式为

或
其复振幅为：
为简化计算起见，上式中假设各次波到达P点时有相同的振幅(不考虑振幅与光程与反比的关系以及华侨因数）。根据惠更斯—菲涅耳原理,将上式对整个缝宽（从x=0到x=b）积分。最后可得沿着衍射角θ方向传播的所有次波在观察点P叠加起来的合振幅:

令，通常称为u的函数，并写成,故P点的光强为

4。夫琅禾费衍射图样的光强分布
当光屏放置在透镜L2的焦平面上时,屏上出现衍射花样，光强的分布可由上式决定。不同的衍射角θ对应于光屏上不同的观察点。首先来决定衍射花样中光强最大值和最小值的位置。即求出满足光强的一阶导数为零的那些点：

由此得
分别解以上两式，可得出所有的极值点。
1、单缝衍射中央最大值的位置：
由，解得满足的一些衍射方向，即 (中间最大值的位置)，也就是在焦点P0处，,光强为最大。这里，叠加的各个次波位相差为零，所以振幅叠加互相加强。
2、单缝衍射最小值的位置:
由,解得满足的一些衍射方向，即
（最小值位置）
时，AP为零,屏上这些点是暗的。

三、光强分布和条纹分布分析
1.基尔霍夫衍射公式的近似
对于一些极简单的衍射问题，也因为被积函数形式复杂而得不到解析形式的积分结果。为此，必须根据实际条件进一步做近似处理。
1。1 傍轴近似
在一般的光学系统中，对成像起主要作用的是那些与光学系统光轴夹角极小的傍轴光线.对于傍轴光线，如图3-3开孔的线度和观察屏上的考察范围都远小于开孔到观察屏的距离。因此，下面的两个近似条件都成立:
1、 ,于是；
2、 .
于是，（3—11）可以化简为

指数中的r未用代替，这是因为指数中r所影响的是次波场的相位，r的微小变化都会引起相位的很大变化。

y
x
Σ Q r P
O

图 3-3 孔径的衍射
1。2距离近似
根据衍射现象,离衍射孔不同距离处，衍射图样是不同的。一种是菲涅耳衍射或近场衍射，指的是光源和接收屏与衍射屏的距离均为有限远，或其中之一是有限远的情形；另一种是夫琅禾费衍射或远场衍射，指的是光源和接收屏与衍射屏的距离均为无限远的情形.
1、 菲涅耳近似 如图3—3所示,设，则由几何关系有

因为在指数上的相位因子决定了函数的周期性，每当相位因子改变时,指数函数反号,这种变化是不可忽略的，相位因子中只有远小于的项才可忽略。
于是，当大到满足

时，上式第三项及以后的各项都可略去，即为

这一近似称为菲涅耳近似，在这个区域内观察到的衍射现象叫菲涅耳衍射.在菲涅耳近似下，P点的光场复振幅为

2、夫琅禾费近似 当观察屏离孔的距离很大，满足

时
这一近似称为夫琅禾费近似，在这个区域内观察到的衍射现象叫夫琅禾费衍射。在夫琅禾费近似下，P点的光场复振幅为:

2.夫琅禾费单缝衍射光强分布
如果只考虑单色平行光垂直入射到开孔平面上的的夫琅禾费衍射，则通常采用图3-4所示的夫琅禾费衍射装置。
单色点光源放置在透镜的前焦平面上，所产生的平行光垂直入射开孔Σ，由开孔的衍射,在透镜的后焦平面上可以观察到开孔Σ的夫琅禾费衍射图样，后焦平面上的各点的光场复振幅由（3-19）式给出.
f
P

M

C

图 3-4 夫琅禾费衍射装置
若开孔面上有均匀的光场分布，令。又因为透镜贴近孔径，。所以后焦平面上的光场复振幅可为
,
对于夫琅禾费单缝衍射，水平缝宽为a，垂直缝宽为b,则,沿y方向的衍射效应不明显，只在x方向有亮暗变化的衍射图样。由（3—20）式，衍射屏上P点的光场复振幅为

式中是观察屏中心点处光场复振幅。
相应P点的光强度为
式中，,。在衍射理论中，通常称为单缝衍射因子，为中央明纹中心处光强度，为单缝边缘光线与中心光线的相位差。
根据(3—22）式可的单缝衍射的光强分布特征:
1、 在中央P点, ,我们使为一很小的趋于零的角，对求极限，则,，为中央主极大光强。
2、 在时，振幅，光强为暗纹的光强。
3、 当时，同理有为一级次极大光强。
4、 当时，同理有为二级次极大光强。
5、 当时，同理有为三级次极大光强.
3.夫琅禾费单缝衍射光强分布
3.1光强分布的极值点
当光屏放置在透镜L2的焦平面上时，屏上出现衍射花样,光强的分布可由(3—22）式决定，也就是说，衍射光场在P点的强度大小主要由因子或决定。不同的衍射角θ对应于光屏上不同的观察点.
1、 当时,该点光强度取最大值：
，称为主极大值；
2、 当或时，该点强度取极小值:
；
3、 相邻两个极小值之间存在一个极大值，由于因子，该极大值的强度总是小于主极大值，故称之为次极大值。
求出满足光强的一阶导数为零的那些点：

解得
于是得到夫琅禾费单缝衍射的次极大值位置满足关系为：。对于这一超越方程其根为

对应的值为

3.2 条纹的角宽度和线宽度
1。暗纹的角宽度
暗纹出现的位置在的地方，对应的值为

任何两相邻暗纹间的衍射角的差值为,即暗纹是以点为中心等间隔左右对称分布的。
2。亮纹的角宽度和线宽度
以相邻暗纹的焦距离作为亮纹的角宽度.在傍轴条件下，即很小时，暗纹的衍射角位置简化为，.中央主极大明纹宽度由的暗纹的衍射角所确定,则主极大亮纹角宽度和线宽度分别为
， （1）
在两相邻暗纹间存在的次极明纹，其角宽度和线宽度分别为
, （2）
3.3 条纹分布的影响因素
从以上结论可以看出,夫琅禾费单缝衍射图样的强度随衍射角度按着函数关系变化；相邻暗条纹中心的角间距相等，因而所有次极大值亮纹的角宽度相等,但主极大值亮纹的角宽度为次极大值的两倍;相邻次极大值亮纹中心不等间距，随着衍射级次的增大，相邻次极大值亮纹中心的间距趋于恒定。
由（1)和（2)式可以看出,条纹分布的影响因素有：
1、 狭缝的宽度对条纹分布的影响
对于给定的波长,亮纹的宽度是与狭缝的宽度成反比的,在波前上对光束限制越大，也就是缝宽越小的时候，衍射场越弥散，衍射图样铺开的越宽；反之，当缝宽很大,光束几乎自由传播时，，这表明衍射场基本上集中在沿直线传播的原方向上，在透镜焦面上衍射图样收缩为几何光学像点。
2、 入射光的波长对条纹分布的影响
在保持缝宽不变的条件下，波长越长,衍射效应越明显;波长越短，衍射效应越可忽略；但是，实际中，由于波长的数量级比较小，所以对衍射条纹的影响不是非常明显.
3、 接收屏的距离对条纹分布的影响
当其他条件不变的情况下，若接收屏沿轴向向前移动时，接收屏上的条纹间距变小，衍射图样展开范围减小;若沿轴向向后移动,接收屏上的条纹间距变大，衍射图样展开范围也增大。

四、 单缝衍射的实验模拟
根据理论的推导，利用衍射积分法对单缝夫琅禾费衍射进行Matlab仿真。衍射积分法是由衍射积分算出接收屏上的光强分布，然后根据该分布调制色彩作图,从而得到衍射图案。
1。夫琅禾费单缝衍射程序和图样
夫琅禾费单缝衍射的仿真程序 ：波长 500nm,缝宽1mm，屏距1m。
clear  
lam=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:ny 
  sinphi=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; 
end 
N=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
[4］ 赵凯华,钟锡华.光学［M].北京大学出版社，1982:65-66。
［5］ 耿顺山。 基于MATLAB的单色光单缝夫琅禾费衍射现象的模拟［J］。 电脑学习。 2008,（01）：20-21
［6] 赵春红。夫琅禾费单狭缝衍射条纹分析［J]。泰州职业技术学院学报。2003，3（03):23-26。

- 19 -