基于matlab的衍射系统仿真.doc

工程光学综合性练习二题目： 基于matlab的衍射系统仿真综合练习大作业二一、要求3-4人组成小组，对下面给出的各题目利用Matlab等工具进行仿真。练习结束时每组提交一份报告及仿真程序。在报告中应注明各仿真结果所对应的参数，如屏与衍射屏间距、孔径形状尺寸等。二、仿真题目1. 改变观察屏与衍射屏间距，观察观察屏上发生的衍射逐渐由菲涅耳衍射转为夫琅和费衍射 1) 原理图：S点光源发出的波长lam=500纳米 S点发出光线经过单缝， 缝宽a； 单缝到衍射屏的距离L2) Matlab代码clear;clcl=10;%l=input (单缝到衍射屏的距离L=);a=0.2;%a=input(单缝的宽度(mm)a=);lam=500e-6;%lam=input(波长(nm);x=-1:0.001:1;%接收屏边界y=x./sqrt(x.2+l2);z=a.*y/lam;I=1000*(sinc(z).2; %计算接受屏某点光强subplot(2,1,1)%绘制仿真图样及强度曲线image(2,x,I)colormap(gray(3)title(单缝衍射条纹)subplot(2,1,2)plot(x,I)title(光强分布)3)初始仿真图样（d=10）4)改变d之后的图样(d=1000)5) 变化规律根据衍射屏以及接受屏的相对位置不同，由此产生菲涅尔衍射和夫琅禾费衍射的区别，根据我们模拟的情况得到菲涅尔衍射和夫琅禾费衍射的明显不同是夫琅禾费衍射条件下：中央有一条特别明亮的亮条纹，其宽度是其他亮条纹的两倍；其他亮条纹的宽度相等，亮度逐渐下降。2. 改变孔径形状、尺寸，观察图样变化 1）原理图 矩孔衍射：透镜焦距：1000mm;照射光波长：500nm;孔高：a（mm）；孔宽：b（mm）；圆孔衍射：圆孔直径：r（mm）；照射光波长：500nm;照射光波长：500nm;2）matlab代码矩孔衍射：focallength=1000;lambda=500;a=2.0;b=2.0;resolution=64;center=(resolution)/2;A=zeros(resolution,resolution);for i=1:1:resolutionfor j=1:1:resolutionif abs(i-center)a*10/2 & abs(j-center)b*10/2A(j,i)=255;endendendE=ones(resolution,resolution);k=2*pi*10000/focallength/lambda;imag=sqrt(-1);for m=1:1:resolutionx=m-center;for n=1:1:resolutiony=n-center;C=ones(resolution,resolution);for i=1:1:resolutionp=i-center;for j=1:1:resolutionq=j-center;C(j,i)=A(j,i)*exp(-imag*k*(x*p+y*q);endendE(n,m)=sum(C(:);endendE=abs(E);I=E.2;I=I.(1/3);I=I.*255/max(max(I);L=I;I=I+256;CM=pink(255).(2/3);gray(255);Colormap(CM);edge=(resolution-1)/20;X,Y=meshgrid(-edge:0.1:edge);x=linspace(-edge,edge,resolution);y=linspace(-edge,edge,resolution);subplot(1,2,1);surf(x,y,L);axis(-edge,edge,-edge,edge,0,255);caxis(0,511);subplot(1,2,2);image(x,y,I);axis(-edge,edge,-edge,edge,0,511);view(2);axis square;圆孔衍射：clearlmda=500e-9; %波长r=1.2e-3; %f = 1; %焦距N = 19;K = linspace(-0.1,0.1,N) ;lmda1 = lmda* ( 1 + K) ;xm = 2000* lmda* f;xs = linspace(-xm,xm,2000) ;ys = xs;z0 = zeros( 2000) ;x,y= meshgrid( xs) ;for i = 1: 19s = 2*pi*r*sqrt(x.2 + y.2)./(lmda1( i) ) ;z = 4* ( besselj( 1,s)./( s + eps) ).2; %光强公式z0 = z0 + z;endz1 = z0 /19;subplot( 1,2,1)imshow( z1* 255) ; %平面图xlabel( x)ylabel( y)subplot( 1,2,2)mesh( x,y,z1) %三维图colormap(gray)xlabel( x)ylabel( y)zlabel( 光强)3）仿真图样：矩孔衍射：a=1，b=2a=2，b=2可知：矩孔在一个维度上展宽一定倍数将导致衍射图样在相同维度上缩短相同倍数，同时能量会更向中心亮斑集中。圆孔衍射：r=1.5r=1.2圆孔半径增大，衍射图样中的同心圆环半径相应变小，更集中，而能量则更向中心亮斑集中。进一步分析各种形状的衍射孔可得如下结论：1.衍射孔在哪些方向上有几条对称轴，则衍射图样亦在哪些方向上有几条对称轴；2.衍射孔的面积越小则衍射图样上的能量分布越分散，越大则能量越集中。3. 多缝的夫琅和费衍射，使用平行光照明，观察衍射图样随缝宽、缝距、缝数量的变化1) 原理图S点光源发出波长lam=500纳米经L1准直透镜产生平行光。缝宽 a缝距 d缝数量 N2) matlab代码 clearlam=500e-9; %定义波长lam=500纳米N=2; %单缝个数Na=2e-4; %缝宽af=5; %衍射屏距接受屏距离zd=5*a; %光栅常数dxm=1*lam*f/a;%接受屏边界y0=xm;n=1001;x0=linspace(-xm,xm,n);%构建接受屏矩阵for i= 1: n %计算接受屏上第i个点的光强sinphi=x0(i)/f;alpha=pi*a*sinphi/lam;beta=pi*d*sinphi/lam; %相位差2*bataB(i,:)=(sin(alpha)./alpha).2.*(sin(N*beta)./sin(beta).2;B1=B/max(B);endNC=255; %绘制衍射仿真图样及强度曲线Br=(B/max(B)*NC;subplot(1,2,1)image(y0,x0,Br);colormap( gray(NC) );subplot(1,2,2)plot(B1,x0);3) 初始仿真图样4) 改变参数后的仿真图样改变缝宽a a=0.0001 a=0.0003改变缝距d d=4*a d=6*a d=7*a缝数量N N=3 N=4 N=65) 变化规律根据仿真图样可以得出以下结论：改变缝宽a 增加缝宽a，亮纹间距减小；反之，亮纹间距增加。改变缝距d 增加缝距d，亮纹数量增加，亮纹间距减小；反之，亮纹数量减少，亮 纹间距增加。缝数量N 增加缝数量N，亮纹变的越细亮；反之，亮纹宽度增加。4. 多缝的夫琅和费衍射，使用平行光照明，观察衍射图样随点光源位置（光源上下移动）的变化1) 原理图与3.相同。（点光源S上下移动，即为准直透镜后平行光与衍射屏入射夹角的变化）2) matlab代码clearlam=500e-9; %定义波长lam=500纳米N=2; %单缝个数Na=2e-4; %缝宽af=5; %衍射屏距接受屏距离zd=5*a; %光栅常数dxm=1*lam*f/a;%接受屏边界r=0; %入射角ry0=xm;n=1001;x0=linspace(-xm,xm,n);%构建接受屏矩阵for i= 1: n %计算接受屏上第i个点的光强sinphi=x0(i)/f;alpha=pi*a*(sinphi-sin(r)/lam;beta=pi*d*(sinphi-sin(r)/lam; %相位差2*bataB(i,:)=(sin(alpha)./alpha).2.*(sin(N*beta)./sin(beta).2;B1=B/max(B);endNC=255; %绘制衍射仿真图样及强度曲线Br=(B/max(B)*NC;subplot(1,2,1)image(y0,x0,Br);colormap( gray(NC) );subplot(1,2,2)plot(B1,x0);3) 初始仿真图样4) 改变参数后的仿真图样（即增加入射角r） 5) 变化规律根据仿真图样可以得出以下结论：移动点光源位置，衍射条纹发生平行移动。具体规律为，点光源向上移动，则条纹向下平行移动，及条纹移动方向与点光源移动方向恰好相反。5. 多缝的夫琅和费衍射，使用扩展光源照明，前后、上下移动光源，观察衍射图样的变化 （1）原理图与4类似。（扩展光源即点光源无数个点光源移动后同时作用，产生的衍射图样的合成）（2）matlab 代码clearlam=500e-9; %定义波长lam=500nmN=2; %单缝个数Na=2e-4; %缝宽af=5; %后面的透镜焦距d=5*a; %光栅常d=5；xm=1*lam*f/a; %设定接受屏边界y0=xm;n=1001;x0=linspace(-xm,xm,n); %构建接受屏矩阵for i= 1: n %依次计算接受屏上的第i个点光强 for j=1:n %依次计算第j个入射点 r=pi*(j-501)/5000000; %各个入射点的角度不同 sinphi=x0(i)/f; alpha=pi*a*(sinphi-sin(r)/lam; beta=pi*d*(sinphi-sin(r)/lam; %相位差beta B(i,j)=(sin(alpha)./alpha).2.*(sin(N*beta)./sin(beta).2; end B2(i) = sum(B(i,:); %进行求和 B1=B2./max(B2);endNC=255; %绘制衍射仿真图样及强度曲线Br=B1*NC;subplot(1,2,1)image(y0,x0,Br);colormap( gray(NC) );subplot(1,2,2)plot(B1,x0);（3）初始仿真图样初始的光源对L1透镜的张角为rad。（4）改变参数，即上下移动光源光源向上移动少许距离后的衍射图样。光源向下移动少许距离后的衍射图样。（5）变化规律由图样可知，当光源向上移动时，图样向下移动，光源向下移动时，图样向上移动。同时，可以知道，扩展光源的光强分布与点光源的有所不同，主要体现在清晰度的变化。三、感想这次大作业要求对课本上的衍射部分有更深层次的理解，通过仿真实验，对菲涅尔衍射和夫琅禾费衍射的特点的理解加深了，对不同衍射孔的夫琅禾费衍射图样特点有了更感性的认识，对光源与衍射图样的关系也有了更深刻的掌握，并进一步熟悉了Matlab在工程光学方面的强大的仿真功能，我们相信对我们以后的学习会有很大的助益。四、参考文献1光学衍射实验的MATLAB仿真A.1004-7077(2010)梁兰菊 田贵才 张裕仕 2基于MATLAB的夫琅和费衍射实验的计算机仿真.张志峰 苏玉玲 蒋逢春 3基于Matlab的光学衍射实验仿真.陈芳芳 钱淑珍 吴平辉毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。 作者签名： 日期： 毕业论文（设计）授权使用说明本论文（设计）作者完全了解*学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。 作者签名： 指导教师签名： 日期： 日期： 注 意 事 项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词 5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它目 录第一章 总 论 11.1项目概况 11.2研究依据及范围 21.3结论 31.4建议 4第二章 项目建设的背景和必要性 52.1项目建设的背景 62.2项目建设的必要性 7第三章 项目服务需求分析 9第四章 项目选址与建设条件 114.1选址原则 114.2项目选址 114.3建设条件 12第五章 建设方案与设计 125.1建设规模与内容 125.2总体规划设计 135.3结构方案 185.4主要配套设备 195.5给水工程 205.6排水工程 225.7电气设计 235.8节能设计 26第六章 项目实施进度和招投标管理296.1项目实施进度 296.2招投标管理 31第七章 环境影响分析 317.1项目主要污染源分析327.2 环境保护措施及治理效果 35第八章 消防、安全与卫生防护 378.1 消防 378.2 劳动安全 388.3 卫生防护 39第九章 组织机构、运作方式与项目实施进度399.1 项目建设组织机构399.2项目运营组织机构 419.3劳动定员 42第十章 投资估算和资金筹措4210.1投资估算4310.2项目所需流动资金 4910.3资金筹措 49第十一章 经济和社会效益评价 5011.1 经济效益评价5011.2 社会效益评价 50第十二章 结 论 50 12.1 主要结论 50 12.2 建 议 51附录：1、*发改局关于*迁建工程项目建议书的批复 2、*村镇建设管理所*迁建项目用地红线