工程光学综合练一 干涉仿真.doc

工程光学综合练习一 干涉仿真 一、要求 3-4人组成小组，对下面给出的各题目利用Matlab等工具进行仿真。练习结束时每组提交一份报告及仿真程序。在报告中应注明各仿真结果所对应的参数，如相干光源间距、光入射倾角等。 二、仿真题目 一）、对于杨氏双缝干涉，改变双缝的缝宽和缝间距，观察干涉图样变化 （1）Matlab程序及注释 建立如下图所示坐标系 两束平面波干涉采用的接收屏是x-y平面 在接收屏上坐标(xs,ys)点处，易求得 ： 光程 相位差 光强 程序： clear lam=500e-9; %设定波长为500nm d=2e-3; %设定双缝宽为2mm； D=1; %设定光源中心到接收屏的距离 xm=5*lam*D/d;ym=xm; %设定光屛的范围 n=101;xs=linspace(-xm,xm,n); %把光屛的x方向分成101点 ys=linspace(-ym,ym,n); %把光屛的y方向分成101点 for i=1:n for j=1:n r1=sqrt((xs(i)-d/2)^2+ys(j)^2+D^2); %光程r1 r2=sqrt((xs(i)+d/2)^2+ys(j)^2+D^2); %光程r2 phi=2*pi*(r2-r1)/lam ; %屏上各点的相位差 B(i,j)=4*cos(phi/2)^2; %屏上各点光强 end end N=225; %确定用的灰度等级为225级 Br=(B/4.0)*N; %使最大光强对应于最大灰度级（白色） subplot(1,2,1) %创建图形窗口 image(xs,ys,Br); %画干涉条纹 xlabel; %y轴方向 ylabel; %x轴方向 title; %屏幕上的干涉图样 运行结果如图1—1 图1—1 将B(i,j)=4*cos(phi/2)^2;改为B(i,j)=8*cos(phi/2)^2;改变强度即改变了缝宽，运行结果如图1—2 将d=2e-3改为d=3e-3即改变缝间距，运行结果如图1—3 图1—2 图1—3 二）、对于杨氏双孔干涉，改变双孔的直径和孔间距，观察干涉图样变化 D=1; %设定光源中心到接收屏的距离 d=2e-6; %设定两光源间距为0.002mm R1=1; R2=1; I0=1; ny=101; lam=5e-7; %设定波长为500nm N=255; I1=I0*R1*R1*R1*R1; I2=I0*R2*R2*R2*R2; ymax=5*lam*D/d; %屏幕上的最大范围 x=y; [x,y]=meshgrid(x,y); %屏幕上的X，Y网格 for i=1:ny for j=1:ny end; end; r1=sqrt((x-d/2).^2+y.^2+(D-d/2).^2); %光源1到接收屏的光程r1 r2=sqrt((x+d/2).^2+y.^2+(D+d/2).^2); %光源2到接收屏的光程r2 phi=(r1-r2)*2*pi/lam;