MATLAB编程用两种方法模拟光学实验Word文件下载.docx

1、正文： 大学教学课程中引入计算机模拟技术正日益受到重视，与Basic、C和Fortran相比，用MATLAB软件做光学试验的模拟，只需要用数学方式表达和描述，省去了大量繁琐的编程过程。下面来介绍利用MATLAB进行光学模拟的两种方法。（一）衍射积分方法： 该方法首先是由衍射积分算出接收屏上的光强分布，然后根据该分布调制色彩作图，从而得到衍射图案。1单缝衍射。把单缝看作是np个分立的相干光源，屏幕上任意一点复振幅为np个光源照射结果的合成，对每个光源，光程差=ypsin,sin=ys/D,光强I=I0（cos）2+（sin）2，其中=2/=ypys/D编写程序如下，得到图1clear lam=5

2、00e-9;a=1e-3;D=1;ym=3*lam*D/a;ny=51;ys=linspace（-ym,ym,ny）;np=51;yp=linspace（0,a,np）;for i=1:ny sinphi=ys（i）/D; alpha=2*pi*yp*sinphi/lam; sumcos=sum（cos（alpha）; sumsin=sum（sin（alpha）; B（i,:）=（sumcos2+sumsin2）/np2;endN=255;Br=（B/max（B）*N;subplot（1,2,1）image（ym,ys,Br）;colormap（gray（N）;subplot（1,2,2）pl

3、ot（B,ys）;2杨氏双缝干涉（gnsh） 两相干光源到接收（jishu）屏上P点距离r1=（D2+（y-a/2）2）1/2, r2=（D2+（y+a/2）2）1/2,相位差=2（r2-r1）/,光强I=4I0cos2（/2）编写程序如下（rxi），得到图2lam=500e-9a=2e-3;ym=5*lam*D/a;xs=ym;n=101;ys=linspace（-ym,ym,n）;n r1=sqrt（ys（i）-a/2）.2+D2）; r2=sqrt（ys（i）+a/2）.2+D2）; phi=2*pi*（r2-r1）./lam;）=sum（4*cos（phi/2）.2）;图2 杨氏双缝干

4、涉的光强分布Br=（B/4.0）*Nimage（xs,ys,Br）;plot（B,ys）3光栅（gungshn）衍射 公式（gngsh）：I=I0（sin/）2（sin（）/sin）2 =（a/）sin=（d/）sin 编写程序如下：得到图3N=2;a=2e-4;D=5;d=5*a;ym=2*lam*D/a;n=1001;图3 黑白光栅衍射光强分布 alpha=pi*a*sinphi/lam; beta=pi*d*sinphi/lam;）=（sin（alpha）./alpha）.2.*（sin（N*beta）./sin（beta）.2; B1=B/max（B）;NC=255;Br=（B/max

5、（B）*NC;colormap（gray（NC）;plot（B1,ys）;（二）傅立叶变换（binhun）方法（fngf）：在傅立叶变换光学中我们（w men）知道夫琅和费衍射场的强度分布就等于屏函数的功率谱。因此我们可以直接（zhji）将光屏进行傅立叶变换，再处理（chl）得到衍射图样。物体图像的生成可以直接由矩阵运算生成，也可利用Windows下的画图工具，生成一幅黑白图像，并调用命令函数imread（）输入图像，输入的图像是一个巨大的二维矩阵，利用MATLAB函数库中的fft2（）命令对该矩阵进行二位离散傅立叶变换，得到图像的频谱，该频谱是一个复数矩阵，然后用取模函数abs（）对该复数矩

6、阵取模，得到振幅谱矩阵，利用函数fftshift（）对取模后的矩阵进行频谱位移，这是因为变换后的二维矩阵的直流分量位于图像的周边角，该函数交换矩阵的1、3象限和2、4象限，使直流分量移到频谱中心，从而使FFT频谱可视效果与实际图像相吻合。最后利用imshow（）函数将图像显示出来。编写程序如下：cleara=imread（E:1.bmp）;grid onfigure（1）imshow（a,）afft=fft2（a）;aabs=abs（afft）;aabss=fftshift（aabs）;figure（2）imshow（aabss,）colormap（gray）;colorbarfigure（3

7、）plot（aabss）figure（4）meshc（aabss）maxx1=max（max（aabss）;输入黑白或灰度的衍射屏图像，得到输出的衍射图样和光强分布。1简单孔径的衍射（a）圆孔衍射：可以看到光强分布随孔径大小的变化（b）矩孔衍射（ynsh（c）三角（snjio）孔衍射：2单缝衍射：3黑白（hibi）光栅（gungshn）衍射： 缝宽不变，光栅（gungshn）常数变大： 缝宽变小，光栅（gungshn）常数不变：可以明显地观察到干涉条纹受衍射因子的调制，以及（yj）缺级现象4正弦光栅衍射： 我们先设计程序制得一张正弦光栅，用正（余）弦函数作图如下：clear allxm=10*

8、pi;ys=xm;xs=linspace（-xm,xm,500）;B=cos（xs）+1;Br=B/2*N将该图片保存成灰度BMP文件。调用程序得到衍射图样，只有三个衍射斑，与理论一致。（三）利用傅立叶方法还可以得到其他（qt）有趣的衍射花样，并研究其规律性 每个方孔几何（j h）尺寸相同（xin tn）可以（ky）看出，基本衍射图样还是方孔衍射（ynsh）但是在任意两个方孔连线的中垂线方向会出现等间距暗纹，三角形摆放则会出现类似蜂巢的结构。 两个方向上的光栅衍射叠加， 如果线条不平行则会在每个 中垂线方向都向外发展。 可以看到，圆孔衍射的图样受到 椭圆孔衍射因子的调制一些（yxi）字母和符号

9、（fho）也可以（ky）产生（chnshng）很有趣的衍射图样：-参考文献：1：胡守信、李柏年，基于（jy）MATLAB的数学实验，科学出版社2：张铮、杨文平、石博强、李海鹏，MATLAB程序设计与实例应用，中国铁道出版社内容总结（1）MATLAB编程用两种方法模拟光学实验摘要：利用MATLAB软件编程实现了用衍射积分的方法对单缝衍射、杨氏双缝干涉、黑白光栅衍射的计算机模拟（2）以及用傅立叶变换方法对简单孔径衍射、黑白光栅及正弦光栅夫琅和费衍射的模拟（3）本文介绍了通过MATLAB软件编程实现用衍射积分和傅立叶变换实现夫琅和费衍射计算机模拟的方法（4）计算机模拟为衍射实验的验证提供一条简捷、直观的途径（5）B1=B/max（B）（6）编写程序如下：