球面波干涉报告书.docVIP

目录 1 概述 1 2 理论分析 1 3 MATLAB程序的编写 2 3.1 点光源为单色光 2 3.2 点光源为复色光 5 4 干涉现象的分析 8 4.1 单色光情况下的变化 8 4.2 复色光情况下的变化 9 5 心得体会 9 6 参考文献 10 7 附录程序 10 球面波干涉的计算机模拟 1 概述 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。 两列球面光波的干涉是光波干涉的最简单，也是最重要的例子。设在均匀媒质中有两个同频率简谐振动的相干点光源S1，S2位于同一平面内，光波叠加后将产生干涉条纹，仿真中，可以在不同的观察平面上观察干涉条纹。如图1所示，a平面与点光源S1，S2的连线垂直，b平面与点光源S1，S2的连线平行。由此可将两列波干涉按照观察平面选取的不同分成两个方向进行讨论，下面我们只对a平面的观察屏进行讨论。 两个点光源与观察平面垂直时，即在平面a处观察干涉条纹。这种情况下，假设两点光源间距为d，观察平面取样点数为ny，把直角坐标系的坐标原点取在平面的中心，点光源连线设为x轴，平面a平行于y-z轴，且平面a与y-z轴距离为D，则点光源S1的位置在（-d/2,0,0），点光源S2的位置在（d/2,0,0），要计算点光源S1和点光源S2到观察平面任意点P处的相位差，然后利用公式 （2-5） 来计算观察平面各个点处的光强分布并显示出来。 3 MATLAB程序的编写 3.1 点光源为单色光 点光源为单色光时，只需要设置光源振幅E=1，然后利用如下程序将干涉强度曲线求出来。 D=0.05 d=1/100； Lambda=500/1000000000； E=1； ymax=0.00055；ny=161；y=linspace(-ymax，ymax，ny)； z=y； i=(ny-1)/2； for j=1：ny l1=sqrt(y(j)^2+(D-d/2)^2+z(j)^2)； l2=sqrt(y(j)^2+(D+d/2)^2+z(j)^2)； phi=2*pi*(l2-l1)/Lambda； I(j)=(E/l1)^2+(E/l2)^2+2*E^2/l1/l2*cos(phi)； end M=max(I)； figure(1)； plot(y,I，b)； axis([-ymax，ymax，0，M])； title(单色光强度曲线) xlabel(x)； ylabel(强度) 运行此程序，结果如图2。 图2 单色光干涉强度曲线图 单色光二维干涉条纹用如下程序编写出来： for i=1ny for j=1：ny l1=sqrt(y(i)^2+(D-d/2)^2+z(j)^2)； l2=sqrt(y(i)^2+(D+d/2)^2+z(j)^2)； phi=2*pi*(l2-l1)/Lambda； I(i,j)=(E/l1)^2+(E/l2)^2+2*E^2/l1/l2*cos(phi)； end end nclevels=255； br=I.*255/max(max(I))； figure(2)； image(y,z,br)； xlabel(y(m))； ylabel(z(m))； title(单色光二维干涉条纹) axis([-ymax,ymax,-ymax,ymax])； colormap(pink(nclevels))； 运行上述程序，结果如图3。 图3 单色光二维干涉条纹图 3.2 点光源为复色光 点光源为复色光时，其与单色光干涉的唯一区别也就是