文章编号:1007-2934(2014)04-基于Matlab的牛顿环白光干涉实验仿真王蕴杰(青海师范大学，青海西宁810008)摘要:根据色度学原理，利用Matlab仿真出了牛顿环白光干涉场的分布图，并给出了程序代码，仿真结果形象直观，效果良好。关键词:干涉;牛顿环;Matlab;计算机仿真中图分类号:O411-3文献标志码:A光的干涉现象对于由光的微粒说到光的波动说的推进起到了重大作用。干涉现象是波的相干迭加的必然结果，其现象证明了光的波动性。利用计算机模拟光的干涉，可以不受仪器、场地的限制，实验效果形象、直观［1］。同时计算机仿真实验可以方便的改变实验参数，从而反映出不同条件下的实验现象和结果，方便了实验教学的开展。由于白光源难于控制和白光的干涉场光强较弱，白光的干涉现象观察有一定的难度。文章中根据色度学原理，应用Matlab仿真了实验室中牛顿环的白光等厚干涉现象［2-3］。板之间形成很薄的空气层。如果有光从上向下照射，那么从空气层的上下表面反射出的两束光线1、2将产生干涉。由于空气层厚度相同的各点位于以透镜顶点O为圆心半径r的圆周上，因此干涉条纹是以O为圆心的若干同心圆。当光源发出的许多波长的光发生干涉时，就会形成彩色的干涉条纹。这种干涉现象的发现者是牛顿，故该装置被称为牛顿环［4］。1．2色度学原理色光的三基色是指红、绿、蓝三色，任何颜色的光可以通过红、绿、蓝三色按照不同的比例合成产生，这是色度学的最基本原理，即三基色原理。因此，可将白光分解为红、绿、蓝三基色，白光产生的干涉条纹可认为是由红、绿、蓝三基色的相干单色光波所产生干涉条纹的非相干叠加。由“CIE(国际照明委员会)1931-ＲGB系统”的规定，白光可被分解为波长分别为700nm、546．1nm、435．8nm的红、绿、蓝三基色光。因此按照色度学原理，可以利用Matlab仿真出牛顿1基本原理1．1牛顿环牛顿环属于等厚干涉现象。它是由一块曲率半径较大的平凸透镜(曲率半径值为Ｒ)和一块平玻璃板放在一起组成，原理见图1。透镜和玻2牛顿环干涉实验仿真设计2．1干涉图像的光强分布对于频率相同的相干光1、2在视场某点p处叠加，其合成光强分布为:(1)I=I1+I2+2槡I1I2cosδ式(1)中，I1和I2分别为两列相干光波到达图1牛顿环收稿日期:2014-04-05基金项目:国家自然科学基金(51167016);青海师范大学2012年教学改革研究项目基于Matlab的牛顿环白光干涉实验仿真98p点时的光强，δ为两列光波的相位差。如果两列光波的波长为λ，考虑到半波损失，由图1可知两列光波的相位差可以用式(2)来表示。2．3牛顿环干涉条纹图像如图2所示为利用Matlab模拟的在透镜曲率半径为5m条件下的干涉条纹，在实验室观测到的干涉条纹如图3所示，比较两图2、图3可以发现:图像的中心接触点位置均为一暗点，其余的条纹分布由内至外均是由紫到红变化，同时经观察可以看到，中心几级与实际情况吻合度是很高，而越到边缘越不一致。造成这个现象的主要原因是:实际图像是波长分布于380～770nm的无数单色光干涉条纹的非相干叠加，而仿真图像为波长分别为700nm、546．1nm、435．8nm的红、绿、蓝三种单色光干涉条纹的非相干叠加，实际应用中必存在着一定差别。δ=2π(2d+λ)(2)λ由几何关系可知:d≈r22(3)2Ｒ式(3)中，d为透镜下空气层的厚度，r为干涉条纹半径，Ｒ为凸透镜的曲率半径。由于两列相干光的光强很接近，因此可以近似认为I1=I2=IO，由式(1)、式(2)可知，两列光波在干涉场p点相遇时的光强为:2π．r2(4)I=4IOsinＲλ牛顿环干涉实验的Matlab仿真将光屏沿水平X方向和竖直方向等分为n×n个点，可分别计算每个点上的光强，各个点距离光屏中心的数值即干涉条纹半径r，由几何关系可知r2=X2+Y2。基于上述原理，在Matlab中设计仿真程序，假设凸透镜的曲率半径为5m，借助Matlab所提供的绘图功能，得到牛顿环等厚干涉条纹的图像，源程序中主要程序段如下:clear;clc;Ｒ=5;wlr=700e-9;wlg=546．1e-9;wlb=435．8e-9;［x，Y］=meshgrid(linspace(-0．005，0．005，400));r2=(x．^2+Y．^2);2．2图2牛顿环干涉图像仿真图deta1=pi*r2/(Ｒ*deta2=pi*r2/(Ｒ*wlr);wlg);wlb);In1=abs(sin(deta1)．^2);In2=abs(sin(deta2)．^2);In3=abs(sin(deta3)．^2)...