基于matlab的衍射系统仿真_.docx

工程光学综合性练习二题目： 基于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上下移动，即为准直透镜后平行光与衍射屏入射夹角r的变化）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） r=蟺10000 r=蟺8000 r=蟺50005) 变化规律根据仿真图样可以得出以下结论：移动点光源位置，衍射条纹发生平行移动。具体规律为，点光源向上移动，则条纹向下平行移动，及条纹移动方向与点光源移动方向恰好相反。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.