基于MATLAB的光孔衍射数值计算-湖北汽车工业学院.docVIP

基于MATLAB的光孔衍射数值计算 杨永富 张慧 但汉作 （理学部 光信息科学与技术专业 ） 指导老师 罗时军 刘国营 曾维友 摘要：编写了夫朗禾费矩形孔、圆孔、椭圆孔、单缝、多缝衍射的MATLAB仿真程序，对这些衍射现象进行了模拟。总结了利用MATLAB对物理光学现象进行模拟的一般方法。 关键词：MATLAB 光孔衍射 数值计算 Numerical Calculation?of Optics-hole Diffraction Based on MATLAB Tutor Luo Shijun,Zeng Weiyou,Liu Guoying Yang Yongfu,Zhang Hui,Dan Hanzuo (Department of Science,Major of Optic-Information Science and Technology) Abstract: It prepared simulation program for the Fraunhofer rectangular,cycle and elliptic hole, single and multi-slit diffraction using the language of MATLAB, these diffraction phenomenon was simulated. Summed up the general approach of simulation of physical optics phenomenon using MATLAB. Keyword: MATLAB; Optics-hole Diffraction; Numerical Calculation 在通常的光学操作实验中，实验参数的改变引起的接收屏上的光强变化不明显，难以充分展示实验的全部特征，加上实验需要特定的仪器和场所，给研究带了许多不便。MATLAB具有数值计算、符号计算、数据可视化、系统动态仿真等功能，这使得它在理论研究和工程实践中得到了广泛的应用。利用MATLAB的数值计算和仿真功能对光学现象进行模拟，这样就可以方便快捷地改变实验参数，观察不同实验结果。本文主要以夫朗禾费衍射为例，利用MATLAB对各种光孔的衍射进行了数值计算。 _______________________________ 基金项目：湖北汽车工业学院大学生创新性实验项目（17） 作者简介：杨永富（1985-）,男，湖北丹江口人，理学部T503-1班学生 1夫朗禾费矩形孔和圆孔衍射的模拟 在夫朗禾费衍射近似下，接收屏上点的光场振幅为(1) 为开孔面上处的光场振幅。 假设开孔面上有均匀的光场分布，即，相应的后焦平面上的光场振幅为(2) 其中 1.1夫朗禾费矩形孔衍射 若衍射孔是矩形孔，根据(2)式，透镜焦平面上点光场振幅为(3) 其中，是观察屏中心点处的光场振幅； 分别是矩形孔沿轴方向的宽度；分别为 则在点的光强度为(4) 式中是点的光强度，且有。 根据以上计算分析编写夫朗禾费矩形孔衍射仿真程序： %This is the program to simulate Fraunhofer Rectangular Diffraction clear all %=================================================== %===模块1：参数设置================================= Lambda=500e-9; b=2e-3; a=2*b; f=500; ymax=1000*b; xs=ymax; ny=500; ys=linspace(-ymax,ymax,ny); xs=linspace(-ymax,ymax,ny); %=================================================== %===模块2：根据光强分布的理论公式建立光强数据矩阵=== for i=1:nyfor j=1:nysin_theta1=xs(i)/f;sin_theta2=ys(j)/f; alpha=pi*a*sin_theta1/Lambda; belta=pi*b*sin_theta2/Lambda; A=sin(alpha)/(alpha+eps); B=sin(belta)/(belta+eps); L(i,j)=4*(A*B).^2;end end %==================================================== %===模块3：根据光强数据矩阵绘制图像================== clf; figure(gcf); ncle