基于MATLAB衍射矩形孔实验模拟.docVIP

基于MATLAB衍射矩形孔实验模拟 　　[摘 要] 根据物理光学中光的矩形孔衍射理论，利用MATLAB工具编写程序进行实验模拟，绘出矩形孔衍射图样和光强分布曲线，直观的展现了矩孔衍射的物理现象，并考虑到实际衍射的情况，对非单色光进行了讨论，结果与理论计算符合。 　　[关键词] MATLAB 矩形孔衍射 　　 　　MATLAB是一种主要用于数值计算即可视化图形处理的工程语言和应用软件。它将数值分析、矩阵计算、图形图像处理、信号处理和仿真等诸多强大的功能集成在较易使用的交互式计算机环境中，为科学研究，工程应用提供了一种功能强、效率高的编程工具。 　　光的衍射是光的波动性的主要标志之一。利用MATLAB模拟实验可以形象直观的演示实验现象，探究光的衍射问题，而且不受实验仪器和实验场所的限制，可以通过改变模拟参数获得不同的仿真结果，并在显示屏上直接显示出来，动态直观地展现各种物理量之间的关系，对于理解光学理论具有积极的作用。 　　1、夫琅禾费衍射实验的理论 　　夫琅禾费衍射装置的光路图如下，其中分别在孔径平面和透镜建立坐标系和，两坐标原点和在透镜光轴上。假设光在衍射屏Σ面上某一波前点Q和观察屏P点的坐标分别是（x1,y1）和（x，y）。在远场条件下满足夫琅禾费近似， 即： 　　则： 　　其中：为的距离；为的距离；为波矢量。 　　图1: 夫琅禾费衍射装置的光路图 　　Fig1: Diagram of Fraunhofer diffraction equipment 　　并假定的是衍射孔受平面波垂直照射，即是常数设为，则光强复振幅的分布满足以下夫琅禾费衍射公式： 　　 （1） 　　式中 　　因此，观察屏上的任意一点的光强分为： 　　 （2） 　　对于不同的情况的衍射屏，只要对衍射孔积分，理论上是可以求出观察屏上任意一点的的光强分布I。 　　2、夫琅禾费矩孔的衍射模拟 　　当衍射孔是长宽分别为a，b，的矩孔时，其主要特征是：衍射亮斑集中分布在相互垂直的两个轴（x轴和y轴）上，并且短边的亮斑的宽度比长边的宽度大，这一点与矩形在两个轴上的宽度关系这好相反。为便于计算在此令： 　　 （3） 　　 （4） 　　则根据夫琅禾费的衍射公式即（1）式可得观察屏上的复振幅为: 　　（5） 　　其中： 　　因此可推出：（6） 　　即简写为： 　　（7） 　　其中 ；是时，即点处的光强。 　　下面我们利用用MATLAB7.10.0(R2010a)模拟观察现象，并验证矩孔衍射的主要特征。 　　模拟时我们取:=1m;a=0.0005m;b=0.008m,=550nm, =1, 程序运行完成后依次得到的衍射图样如图2所示： 　　图2：矩形孔单色光的衍射图样 　　 Figure 2: The rectangular hole diffraction pattern of monochromatic light 　　由图j可以直观清晰的看出衍射图样，图k则形象的反映了衍射光强的三维的分布，图l和图m分别画出的是宽度a方向和长度b方向的光强分布，从图中可以容易的看出衍射宽度和长度光强的分布与上述矩孔的衍射特征相一致。 　　进一步考虑，在实际情况下，纯粹的单色光不易获得，通常都有一定的波谱宽度。为了便于模拟，并接近于实际情况，在围绕中心波长取区间长度，并假定光波波长在该区间内均匀分布，可取中心波长的10%，则现在的波长就是一个变量，我们可以把该区间均匀的分成20份，即取其中的19根谱线，则现在波长就是一个一维向量。计算光强时应把这19 根谱线产生的光强叠加并取平均值，即 　　 （8） 　　则通过MATLAB7.10.0(R2010a)很方便的实现这一思想，程序运行完成后依次得到的衍射图样如图3所示： 　　图3：矩形孔非单色光的衍射图样 　　Figure 3: The rectangular hole diffraction pattern of non- monochromatic light 　　从以上衍射图样中可以看出非单色光的衍射现象相对要弱，边缘波长的衍射现象不明显。 　　3、夫琅禾费单缝的衍射模拟 　　进一步我们可虑矩形孔的一个方向的宽度比另一个方向的宽度要大得多，时， 矩形孔就变成了狭缝。入射光在y方向的衍射效应可以忽略，衍射图样只分布在x轴上。下面我们也很方便的用MATLAB7.10.0(R2010a)进行模拟，画出x方向的光强分布。则单色光和非单色光的单缝衍射图样图4所示： 　　图4：单缝衍射图样（单色光/非单色光） 　　Fig4: Single slit diffraction pattern（monochromatic light / non- monochromati