利用MATLAB语言进行光学衍射现象的仿真

《利用MATLAB语言进行光学衍射现象的仿真》由会员分享，可在线阅读，更多相关《利用MATLAB语言进行光学衍射现象的仿真（12页珍藏版）》请在装配图网上搜索。

1、利用MATLAB语言进行光学衍射现象的仿真储林华（安庆师范学院物理与电气工程学院 安徽 安庆 246011）指导教师：张杰摘要： 光的衍射是光的波动性的一种重要表现，因此对光的衍射现象的研究，不仅具有重要的理论意 义，而且在光学仪器研制和成像分析等诸多实际应用方面均有重要价值，但是其衍射光强的计算非常复 杂，对实验条件的要求非常高，通常情况下很难得到满意的效果，严重影响了光学的教学。本文从衍射 的相关理论知识出发，首先介绍了惠更斯-菲涅耳原理及其数学表示形式，然后重点讨论了单色光经各 种对称光学衍射元件（单缝，双缝，光栅，圆孔）的夫琅和费衍射情况，并分别给出了它们在焦平面上 的衍射光强计算公式

2、，最后利用科学计算软件MATLAB对光的衍射现象进行了仿真，所得到的图样细 致逼真，使整个物理过程变得直观形象，且与实验所得到的衍射图样进行了比较，两者吻合得很好，从 而为光学的理论分析和实验教学提供了一种新的途径。关键词：光的衍射，光栅衍射，圆孔衍射，Matlab，计算机仿真0 引言光的衍射现象是光具有波动性的重要特征，因此对衍射现象的研究无论在理论上还是在实践中都 有很重要的意义。对光的衍射现象的研究，始于17世纪，当时著名的荷兰科学家惠更斯提出了光是一 种波的假说，并根据波动理论提出了光的传播理论即惠更斯原理1，根据这一原理，他解释了光的 反射定律和折射定律，给出了折射率的意义，光在两种

3、介质中的速度比。到了19世纪，法国年轻的科 学家菲涅耳，根据叠加原理把惠更斯原理进一步具体化，给出了光在传播过程中光强学计算公式，这就 是著名的惠更斯-菲涅耳原理2。但由于在实际应用过程中，障碍物形状的不规则性，导致光强的计算 公式几乎无解析解，只能进行一些数值计算。针对衍射计算中出现的困难，近代的研究人员想到运用科学的计算软件MATLAB，利用其较强的 绘图和图象功能，编写计算程序，使得多种衍射元件（单缝，双缝，光栅，矩孔，圆孔）下的衍射现象 得以在计算机中形象地被模拟仿真。这种做法，条件限制较少，对于衍射的实验教学是一种的补充，起 到了很不错的效果；但值得指出的是，许多前人撰写的论文，总是

4、在系统化，可视化，条理化方面不够 理想，本文将在他们工作的基础上，将此课题进一步做得完美。1惠更斯菲涅耳原理早在十七世纪后期，荷兰科学家惠更斯就提出了光是一种波动的假说，并阐述了关于波面传播的 一种理论，既惠更斯原理。该原理认为，传播中的波面上任何一点都可以认为是一个新的次波源，由于 这些次波源发出的次波是球面波，这些次波的公共包络面就是下一时刻的波面，根据这一原理，他解释 了光的反射定律和折射定律，并给出了折射率的意义，光在两种介质中的速度比。图1波面上各点复振幅的传播du (p) = cu(Q) f )ejkrdsr菲涅耳根据叠加原理把惠更斯原理进一步具体化了，他 假设各次波是球面波，但这

5、些球面只是等相面而不是等幅面， 球面上各点振幅与传播方向有关,这就避免了次波的后向的传 播；同时，他认为，下一时刻空间任一点的振动由各次波到达 该点的振动叠加所决定，此定理称为惠更斯一菲涅耳定理。如 右图(1)所示。菲涅耳假设:Q点所发的次波对P点贡献dU(p) 正比于Q点附近一个小面元的面积ds，正比于Q点的复振 幅U(Q),正比于ejk/r (假定次波是球面波)，以及正比于一个与传播方向有关的函数f (购是r与小面元 法线即波在Q点的传播方向上的夹角)，即：其中C是一个与r，Q，0无关的比例系数，f(0 )叫倾斜因子，它随e增加而缓慢减少于是，按照叠加 原理，有：(1)J du(p) Jc

6、u(Q)f (e ) ejkrdsrss这就是惠更斯一菲涅耳原理的数学表达式，积分表示整个波面S上各点所发的次波传播到P点的作用 的叠加。从上面的表述，我们可以看到菲涅耳的思想比惠更斯有了更大的进步，他着眼于下一时刻空间各 点的振动情况，而惠更斯只着眼于下一时刻波面的形状与位置，因此惠更斯只能定性地描述光的传播方 向，而菲涅耳却能定量地描述衍射后的光强分布。2. 夫朗禾费衍射现象研究光的衍射现象根据光源到衍射屏以及观察屏的距离远近，可以分为近场衍射和远场衍射。如果光源 到衍射屏以及观察屏的距离为有限值，则称为近场衍射(菲涅耳衍射)；如果光源到衍射屏以及观察屏 的距离为无穷远，则称远场衍射，由于

7、这种衍射最先由夫朗禾费在探索天体成像时作了系统的研究，故 亦称为夫朗禾费衍射。一般情况下，利用(1)式进行光强计算时，菲涅耳衍射情况比较复杂，而夫朗 禾费衍射情况比较简单，本文仅讨论后者。2.1圆孔的夫朗禾费衍射在Oxy平面上，有一以O为圆心、R 为半径的圆孔，如图2所示。现用一束平行 于Z轴的光线照射，经圆孔衍射后到达焦平 面上，于是按照(1)式，P点的光强为：u(p) = J cU(Q)f(0)ejkQdsrQ圆孔Q=J 予 cu(Q)f(0)ejkrQpdP却rp=0 屮=0Q设从圆心O到P点的距离为r0，则从圆孔中任一点Q到P点的距离rQ为:r = r + oo sin0 = r +

8、p cosQ sin0Q 00于是，由于程差A = r -r比r小得多，而与波长九可以比拟，因此，可以把振幅中的r以r代替而只 Q 00Q 0考虑相位因子的变化，即：其中:u(Q) f (0 )u (p) = cejk( rQ-r0)p d p 如(2)ejk0为常数，将Tq代入(2)式中得:其中：J1称为一阶贝塞尔函数，X = 2兀Rsin0.认，由此，可知P点的光强可表达为:I (p) = U (p)|2 = C2冗 2R4(3)2 J (x) 21X根据(3)式可以画出焦平面的光强分布图样，如下图(3)所示:(a)(b)图3圆孔衍射图样其中，图(b)是夫朗禾费衍射图样的照片。显然，这是一

9、个圆对称图形，中心为主极大，在F处,R = 0，0 = 0, x = 0,I = I。，其他各点的光强可通过R来表示，由图(2)可知:R = f sin 0xf入2兀R当:x = 1.22兀,II = 0,为第一暗环;R= 0.818 f R 时,x=1.635兀,I/I = 0.0175,；0R= 1.116 f 时，x 二 2.233兀,1匕 二 0,R = 1.619 f/R 时，x 二 3.238兀,11。= 0 ,为第一亮环；为第二暗环;为第三暗环。这就是夫朗禾费圆孔衍射的光强分布，它的中心永远是亮的，并且在中心亮斑处的光能占总光能 的约84%，中心亮斑的半径也是确定的，为0.61九

10、fJR，或角半径为0.612 R。2.2单缝的夫朗禾费衍射图4单缝夫朗禾费衍射示意图设波长为九的平面波射向缝宽AB=b的狭缝，衍射经透镜L会聚在焦平面F上，取XZ坐标如图 (4)所示，根据惠一菲原理，在焦平面上任一点P的复振幅为：u(p)二 J cU(Q)f (Q )eJkrdsr狭缝把狭缝细分为垂直于X轴的许多小面元，面积为ds=Ldx, L是缝的长度，由于是平面波入射,故U(Q)为常数在角度不大的情况下f (6) - lr - r0 + xsin0，由于x r0，故在振幅中的r可视为1只有相位因子中的“ r- r0 - xsin0不可略去，故:e jkr0 b 2u (p) = cLu (

11、Q)J ejkxsin0dxr .0-b 2把积分求出，得：bb2cbu(p) = c/jksin0 (ejk2sin0 -e-j0)=时sin/sin0)令:kfsin0 =兀饮sin0 = P贝y： u(p)二cbsin片帘I = u(p)u* (p) = c2b2 sin2 p / p 2 = I sin2 P 2其中，10为p= 0时的光强，即：衍射斑中心点的光强，p的物理意义是狭缝的边缘与中心的光线在P点产生的相位差，以卩或sin 9作横坐标，以I/。作纵坐标，对(4)式作图可得夫朗禾费单缝衍 射的光强分布，如下图所示：(b)单缝衍射图样图5单缝衍射光强分布及图样计算表明，在-兀 卩

12、 兀之间的主极大集中了 90%的能量，主极大的半角宽为9 - sin9 b，把狭缝改为矩孔，即在x，y两个方向上考虑衍射效应，则光强表达式为：I二Iosin2 B沖2 sin22 其中，a是在y方向上与P对应的量，矩孔衍射图样如图(6)所示,图6矩孔衍射图样2.3双缝的夫朗禾费衍射以上考虑的是一个缝的衍射，如果有两个相邻的缝，由于衍射，通过两个缝的光在观察屏上会相 遇。试验结果告诉我们，在两条缝的衍射光相互交叠的区域，不是简单的呈现光强的叠加，而是出现了 由光强重新分布而产生的明暗相间的条纹，称为光的干涉。显然，光的干涉与光的衍射一样，都是波的叠加原理所必然导致的结果。现在，我们仍然用惠一菲原

13、理来分析两条缝所产生的夫朗禾费衍射的结果，如图(7)所示：图7双缝衍射示意图如图(7),衍射屏上A，B处各有一条宽为b的缝，缝间距为d，由于透镜的作用，这两条缝的衍 射光在焦平面F上的光强分布是完全一样的，但它们的相位分布不同，把坐标原点分别放在A与B的中心，由式：u(p) = c J ejkxsin9dx可得：-b 2u (p) = c JA-b 2ejk ( xa + ；)sin 9dxAu (p) = c J ejk(xB-2)sin9dxBB-b 2令：r = ksin9为单缝中心与双缝中心的光在P点产生的相位差，则:|b 2u (p) = u (p) + u (P) = C I ej

14、kxsinOdx(e + e-)AB-b2=-c b sin B P (ejr + e-)于是，总光强为：I = u (p)u * (p) = I sin 2 p . p 2 (ejr + e - jr )(e- + e)o= 4I sin20P / P2 cos2 r5)这就是夫朗禾费双缝衍射的光强分布的表达式，设d=5b,即：r=5 P，将式(5)作图，可得下图所示的光强分布图样：43210-10-50510 1520-20 -15图 8 双缝衍射光强分布图样显然，这是由sin2 P.P2与4cos2r相乘而得的图样，在有双缝时，原来的单缝衍射图形不是平均地增加到原来的两倍，而是在sin0

15、 =賦d时，增加到原来的四倍，在sin0 = (2k +1)九/2d处，降为零。2.4 多缝夫朗禾费衍射上面关于双缝衍射的讨论，可以进一步推广到更多的缝，如图(9)，有 N 条等间距的缝，缝宽均为 d，则相邻缝的对应程差为A = d sin0，相位差为：g= 2兀d sin0 /九，由式知：u (p) = u (p)(ejr + e-jr) = u (p)ejr (1+ e-2jr) = u (p)(1+ e j)0 0 0图9多缝夫朗禾费衍射示意图如果把坐标原点放在第一个缝的中心，则u (p)e就是它的单缝衍射振幅，而u (p)ej(r-6则是 0 0另一个缝的衍射振幅，根据这个思想，立即可

16、以把上式推广到多缝：u (P) = u ( P ) + u ( P ) HF u ( P )1 一 ejN a1 ejaN12N=u (p)(1+ e- j + e-2 j ff e-(ni)j ) = u (p)ii于是，光强分布为:IN (P) = uN (P)u N (P)=u1(P)u(P)(1 ejNa )(1 e-jNa)(1 eja )(1 e - ja)(6)f/J. IIS1412耳864 2r-2 n 3x/2 j 兀k/20N =4= 10p/2A ； Tk/2 _ 2 1=I sin2 0 . 0 2sin2(Na /2) sin2 (a /2)0设d=5b，N=4，按

17、(5)式画出光强分布如下图所示：图10多缝衍射光强分布图 进而若n=1oo,贝y可得到的衍射图样为：20100008000600040002000-15-10-505101520图 11 光栅衍射图样分析：多缝衍射与双缝衍射都是单缝衍射的结果。它们的区别在于多缝是干涉的结果，使极大值 变细，峰值变高；而在两个极大值之间出现(N1)个光强为零的B值和(N2)个次峰。当N很大时，各 主极大变得十分尖锐，可以用来计量和分光，这种 N 很大的能产生多缝衍射的光学元件称为光栅， d 叫做光栅常数。3. 利用Matlab软件进行衍射现象的仿真3.1 Matlab 的简介这里使用的软件 MATLAB6.5，

18、 它是一个功能十分强大的应用软件，可以在很多学科中得到应用， 它集数值计算，符号计算，数据可视化，系统动态仿真于一体，与其他计算机语言相比它更加灵活，更 加接近科技人员的思维方式，因而编程效率更高。3.2 仿真方法下面我们应用Matlab强大的图像处理与数值计算功能,仿真模拟光的夫朗禾费衍射(各种衍射屏) 图样和强度分布曲线。仿真模拟首先是根据光的衍射光强分布的理论公式I(x, y)及实验参数建立光强 数据矩阵B(x, y),然后运用Matlab的Image命令绘制衍射图像，运用Plot命令绘制光强分布曲线。3.3 圆孔衍射的程序% 圆孔衍射% function ykys (r,f,lambd

19、a)r= 0.03;f=600;lambda=600;x= -20:0.05:20;y=-20:0.05:20;x1,y1=meshgrid(x,y);theta=atan(sqrt(x1.A2+y1.A2)./f);x=1000000*(2*pi*r/lambda).*sin(theta+(theta=0)*eps);I=(2*bessel(1,x)./x).人 2;subplot(2,1,1);mesh(x1,y1,I); axis(-20 20 -20 20 0 1) subplot (2,1,2);subimage (I*255)axis off圆孔衍射图样照片见图3 (b)所示。3.

20、4 光栅衍射程序% 光栅衍射% function gsys (b， d， f， lambda，N)b=0.05;d=0.1;f=600;lambda=600;N=1/2/100; x=-20:0.1:20;theta =atan (x./f); beta=1000000*(pi*b/lambda).*sin(theta+(theta=0)*eps);delta=beta.*(d/b);I=(sin(beta)./beta).人2).*(sin(N.*delta)./sin(delta).人2); subplot(2,1,1);plot(x,I);subplot(2, 1 ,2);colorma

21、p(gray(10)image(I*255);axis off运行后得到的图样分别见图5(b),图8,图11所示。3.5 矩孔衍射的程序% 矩孔衍射% function jkys (a,b,f,lambda)a=0.05;b=0.05;f=600;lambda=600;x= -20:0.05:20;y= -20:0.05:20;X,Y=meshgrid(x,y);theta1=atan (X./f);theta2=atan(Y./f); beta=1000000*(pi*b/lambda).*sin(theta1+(theta1=0)*eps);alpha=1000000*(pi*a/lamb

22、da).*sin(theta2+(theta2=0)*eps);I=(sin(beta)./beta).A2).*(sin(alpha)./alpha).A2);subplot(2,1,1);mesh(X,Y,I);axis( -20 20 -20 20 0 1) subplot(2,1,2);subimage(I*255); 运行后得到矩孔衍射图样见图6 所示。4结束语：经过Matlab软件的仿真，得到各种形状孔屏的夫朗禾费衍射图样，这与实验室做实验得到的 结果符合得很好；而且通过仿真，也更进一步使夫朗禾费衍射的相关知识形象化，可视化，加深了 对这些内容的理解。参考文献：1 章志鸣，光学，北

23、京:高等教育出版社,2000.2 姚启钧，光学教程，北京:高等教育出版社,2002.3 谢嘉宁，陈伟成，matlab在光学信息处理仿真实验中的应用，物理实验，Vol. 24,No. 6,2004:2326.4 谢嘉宁，赵建林，夫朗禾费衍射的计算机仿真，大学物理， Vol. 23,No. 3,2004:5154 楼顺天，于卫，阎华粱，matlab程序设计语言，西安：西安电子科技大学出版社，2000 孙向阳，圆孔衍射与象分辨本领的matlab仿真分析，湖北工业大学学报，Vol. 19, No. 5, 2004： 4749. 王炳武，MATLAB5.3实际教程，北京：中国水利水电出版社，2000.8

24、 赵凯华，钟锡华，光学，北京大学出版社， 1984.9 Grant R. Fowles. Introduction to Modern Optics. New York: Dover Publications, Inc. ,1989.10 Portmann S, Luthi H P. Molekel: An Interractive Molecular Graphics Tool. CHIMIA, 2000.Simulation of Light Diffraction by MATLABLinhua Chu(School of Physics and Electrical Engineeri

25、ng of Anqing Normal College， Anqing 246011)Abstact: The lights diffraction is a very important behavior of the lights fluctuation character. So the study of the lights diffraction phenomenon has not only the important theory value but also the important value on many practical using，such as optics i

26、nstruments development and imagery analysis and so on. However， the calculation of the diffraction light intensity is very complicated and has a very high standard on the experiment condition. Usually, its not easy to get the pleased results, which affects the optics teaching seriously. This thesis

27、is based on the related reason knowledge of diffraction. Firstly， it introduces Huygens-Fresnel principle and the math expression form. Then it emphatically discusses the single frequencys light of Fraunhofer diffraction cases through many kinds of the symmetric element of the ligh ts diffraction (s

28、ingle crack, double crack, grating, rectangular hole, round hole),furthermore, the thesis gives obtain the diffraction light intensity calculation on the focal plane. Finally，it uses the scientific calculation software MATLAB to simulate the lights diffraction phenomenon. The pattern we got is delic

29、ate and realistic， which makes the whole physical process intuitional and visual. And compared it with the diffraction pattern which is be gotten from experiment, they coincide with each other well， which supply a new way for the optics theory analysis and experiment teaching.Key words: lights diffraction，grating diffraction，round hole diffraction，Matlab ，computer simulation