《光学现象的Matlab演示.doc》由会员分享,可在线阅读,更多相关《光学现象的Matlab演示.doc(30页珍藏版)》请在金锄头文库上搜索。
1、苏州大学本科生毕业设计(论文)目 录摘要 1第一章 引言 2第二章Matlab的使用 32.1 命令运行方式 32.2 M文件运行方式 3第三章 波动光学基本理论 53.1 光的衍射 53.1.1夫朗禾费单缝衍射 53.1.2夫朗禾费多缝衍射 83.1.3夫朗禾费圆孔衍射 93.1.4夫朗禾费矩孔衍射 113.2 光的干涉 123.2.1杨氏双缝干涉 123.2.2薄膜干涉(等厚干涉) 153.2.3牛顿环 173.3 光学现象的Matlab模拟 203.3.1夫朗禾费单缝衍射 203.3.2夫朗禾费多缝衍射 203.3.3夫朗禾费圆孔衍射 213.3.4夫朗禾费矩孔衍射 233.3.5杨氏双
2、缝隙干涉 243.3.6薄膜干涉(等厚干涉)243.3.7牛顿环25第四章 结论27参考文献28 光学现象的Matlab演示 指导老师 刘琳苏州大学应用技术学院学院 机电系06级 0616407035 摘要:Matlab强大的矩阵处理与运算功能,丰富的图形绘制能力深受用户的青睐,其应用范围涵盖了当今几乎所有的工业应用与科学研究领域。它集数值分析、矩阵运算、信号处理和图形显示于一体,其丰富的库函数和各种专用工具箱,将使用者从繁琐的底层编程中解放出来;它对科学计算结果迅捷而准确的可视化能力,有助于使用者化抽象思维为形象思维,更好地理解概念。将Matlab软件和光学有机的结合起来,能够直观的建立物理
3、现象,更形象的说明问题。以下主要分析介绍了夫朗禾费单缝衍射,多缝衍射,矩孔衍射,和圆孔衍射以及杨氏双缝干涉,薄膜干涉和牛顿环的干涉等理论,用Matlab软件编写相应程序然后进行计算机模拟,有助于理解和研究衍射和干涉的理论。关键词:Matlab;波动光学;程序设计;计算机模拟Abstract: Matlab is widely used by the image processing and computer vision community, including industrial application and the scientific research .Formidable mat
4、rix processing , the operation function, the rich graph plan ability depth, advanced algorithms and numerical accuracy and access to implementation details are all its superiorities. It made significant contributions to physics. Fraunhofer signal slit diffranction 、multi slits diffranction, circular
5、 aperture diffranction, rectangular aperture diffranction ,Youngs two slits interference ,equal thickness interference and Newton rings are analyzed in this article .Then Matlab is used to compile the corresponding procedures and to simulate by computer. This is useful to understand and study of the
6、 diffraction and interference theory.Key words: Matlab; wave optics; programming; simulating with computer第一章 引言Matlab是美国Mathworks公司于20世纪80年代推出的一种简便的工程计算语言,其应用范围涵盖了当今几乎所有的工业应用与科学研究领域。鉴于Matlab的诸多优点,使它深受工程技术人员及科技专家的欢迎,并很快成为计算机辅助分析、设计、仿真、教学等领域不可缺少的基础软件之一。 光的波动性通常表现为光的干涉,衍射,偏振等,波动光学课程中包含有丰富的重要的光学图像,许多重要
7、的规律都是建立在这些物理现象之上的,Matlab软件,其强大的科学计算和图形图像功能的完美结合使之成为我们进行科学研究和教学的首选。特别是通过人机交互任意改变各参量值, 使我们对光学图样变化与各参量间的关系可以有一个直观感受,从而可以加深对光学现象的认识,更好地理解物理概念。第二章 Matlab的使用 Matlab提供了两种运行方式,即命令方式和M文件方式,两种方式各有特点。2.1 命令运行方式即可以通过直接在命令窗口输入命令来实行计算机或作图功能。例如,要求矩阵A和B的和,其中 首先打开Matlab界面,直接在命令窗口输入下面的命令行A= 2 3;4 9;B=-3 5;-6 7;C=A+B最
8、终显示的运算结果如下图2-1所示 图2-1 Matlab命令运行方式界面2.2 M文件运行方式 在Matlab窗口中单击File菜单,然后依次选择New-File,打开M文件输入运行界面,如图2-2所示。在该窗口中输入程序文件,可以进行调试或运行,与命令行方式相比,M文件方式的优点是可以调试,可以重复应用。 对于前面的矩阵求和问题,在M文件输入运行界面中输入程序,如图2-2所示。然后在Debug菜单中选择Run选项将在命令窗口中输出矩阵C=A+B的值。 图2-2 M文件运行图示M文件有两种格式,即函数式M文件和脚本式M文件。函数式M文件的第一句式以Function语句作为引导的,脚本式M文件就
9、是命令的简单叠加,它与批处理文件很相似。第三章 波动光学基本理论 第3.1节 光的衍射 光在传播路径中,遇到不透明或透明的障碍物,绕过障碍物,产生偏离直线传播的现象称为光的衍射。按光源、衍射屏和接收屏三者之间的相对位置,可将衍射现象分为两种类型菲涅耳衍射和夫琅禾费衍射。光源和观察点距障碍物为有限远的衍射称为菲涅尔衍射。菲涅耳衍射图样是带有衍射条纹的衍射孔的投影像。光源和观察点距障碍物为无限远,即平行光的衍射为夫朗禾费衍射。夫琅禾费衍射图样是带衍射条纹的光源的投影像。3.1.1夫朗禾费单缝衍射由基础光学可知,任意衍射屏的夫琅禾费衍射可借助两个透镜来实现如图3-1-l所示,位于透镜L。物方焦平面上
10、的点源S所发出的单色球面光波经L。变换为一束平面光波,照射在衍射屏AB上。按照平面波理论,衍射屏开口处的波前向各个方向发出次波,方向彼此相同的衍射次波经透镜 会聚到其像方焦平面的同一点P上。满足相长干涉条件时,该点为亮点;满足相消干涉条件时,该点为暗点。所有亮点和暗点的集合构成了该衍射屏的夫琅禾费衍射图样。其次,从傅里叶光学角度,任意衍射屏在单位振幅的单色平面波垂直照射下,其夫琅禾费衍射光场复振幅即衍射屏透射系数的傅里叶变换,而衍射图样实际上就是衍射屏的空间频谱强度分布。BaAXfz 图3-1-1 夫朗禾费衍射实验装置图求接受屏上的衍射强度分布可以通过求解衍射几分公式,得到屏上的复振幅分布,然
11、后再计算光强分布。由于衍射孔径,即光强分布比较小,因而可以认为衍射光是满足近轴条件的,我们可以采用数值积分法。 P点光来自同一方向,倾斜因子相同,满足近轴条件,倾斜因子为常数1,即所有F=1.狭缝上各点的瞳函数有相等的振幅,记作,如果记狭缝中心处O点的相位为,则 (3.1.1) 记,则(3.1.1)式可化为 (3.1.2)从O点到P点的光程记作,则x到P点的光程为,由图3-1-2可见,于是衍射积分公式为 (3.1.3)fB-a/2rPFAa/2xx图 3-1-2 积分法示意图近轴条件下,各次波中心所发出的球面波的振幅相等,即(2)式所示的积分公式中表示球面波振幅因子为常数,记作,积分公式进一步
12、化为 (3.1.4) 其中,为狭缝上Q点附近单位宽度光源发出的沿光轴方向的次波在光轴上的F点(焦点)所引起的复振幅,,为通过整个狭缝的、沿光轴方向传播时在光轴上的焦点所引起的振动,即复振幅。则 为光轴上焦点处的光强。 强度分布如(3.1.5)式 (3.1.5)上所示即为夫朗禾费衍射光强公式。,对于沿光轴方向入射的光,=0,为为单缝衍射因子。用Matlab编程进行模拟,得到的强度分布曲线和模拟图如图3-1-3所示,即为明暗交错的条纹,且缝宽越大,衍射条纹越细。图3-1-3 单缝衍射模拟及强度分布3.1.2夫朗禾费多缝衍射如果单缝沿着衍射屏平移,而衍射装置的其它部分不变,则衍射的强度分布将不会发生
13、改变。因而对于多缝衍射屏来说,其中每一个单缝,即每一个衍射单元在接收屏上所产生的衍射条纹都是相同的。但来自不同狭缝的光由于相干性,相互间进行相干叠加,对于衍射光栅来说,既有来自每一个衍射单元的波列各自的衍射,也有来自不同单元(狭缝)的波列之间的干涉。如果光栅只有两条狭缝,N=2,则衍射光强为 (3.1.6) 用Matlab编程进行模拟,得到的强度分布曲线和模拟图如图3-1-4所示。图3-1-4 多缝衍射模拟及强度分布图3.1.3夫朗禾费圆孔衍射xyzRABO 图3-1-5 夫朗禾费圆孔衍射如图3-1-5所示,Q点发出任意方向光线,与光轴间的夹角为,过中心O作与同方向光线,取坐标系如下:和轴线所在平面为XOZ平面,Z为光轴,过Q作与、垂直的平面,与和X轴分别交与B、A点。则AB与垂直,与YOZ平面的夹角为,A,Q两点发出的次波是等光程的。则任意点Q点发出的次波与中心点O发出的次波间的光程差为在上的投影,即 (3.1.7) 在近轴条件下,焦平面上P点的复振幅为(取实部) (3.1.8)令 则上式可化为 为一阶贝塞尔函数, 光强