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

工程光学综合练习一 干涉仿真 一、要求 3-4 人组成小组，对下面给出的各题目利用 Matlab 等工具进行仿真。练习结 束时每组提交一份报告及仿真程序。在报告中应注明各仿真结果所对应的参数，如相干光源间距、光入射倾角等。 二、仿真题目 一）、对于杨氏双缝干涉，改变双缝的缝宽和缝间距，观察干涉图样变化 1） Matlab 程序及注释 建立如下图所示坐标系 两束平面波干涉采用的接收屏是 x-y平面 在接收屏上坐标 (xs,ys)点处，易求得 ： 光程 r( xs d )2 ys2 D 2 2 2 相位差 phi 2 光强 B 4cos( ) 程序： clear lam=500e-9; d=2e-3;  %设定波长为 500nm %设定双缝宽为 2mm； D=1;  %设定光源中心到接收屏的 距离 xm=5*lam*D/d;ym=xm; n=101;xs=linspace(-xm,xm,n); ys=linspace(-ym,ym,n);  %设定光屛的范围 %把光屛的 x方向分成 101点 %把光屛的 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); phi=2*pi*(r2-r1)/lam ;  %光程 r2 %屏上各点的相位差 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; ylabel;  %创建图形窗口 %画干涉条纹 %y轴方向 %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; % 屏幕上 y的最大范围 y=linspace(-ymax,ymax,ny)% 设定光屛的范围 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; % 光屏上各点相位差 I=I1+I2+2*sqrt(I1*I2)*cos(phi); % 光屏上各点光强 %end %end Br=(I/4.0)*N; % 使最大光强对应于最大灰度级 （白色） figure(1); image(Br); colormap(gray(N)); % 根据光强生成图像 生成图像如图 2-1 所示 将参数 d=2e-6 改为 d=2e-5, 即增大两光源的距离，所得图像如图 2-2 所示 图 2-1 图 2-2 三）、改变下列光波场分布，观察干涉图样变化 S1 S2 S1 S2 d z  d 图1 图2 图3 1、如图 1 所示，两平面光波叠加，改变光波振幅比、两光波夹角，观察在接收屏上的干涉图样变化； 1） Matlab 程序及注释 A1=input( 第一束光的振幅 (m)=); A2=input( 第二束光的振幅 (m)=); jiajiao=input( 夹角 =); lanbda=input( 光波长（ nm） =)/1000000000; ymax=0 ny=161; y=linspace(-ymax,ymax,ny); z=y; i=(ny-1)/2; for i=1:ny j=1:ny I1=A1^2; %计算第一束光的光强 I2=A2^2; %计算第二束光的光强 phi=pi*jiajiao*y(i)/lanbd