MATLAB光学仿真-楔板

1、5、用平行光(点光源+准直镜)照射如图7所示不同形状楔板，观察干涉图样(1)二%清除原有变量%设定波长%设定倾角%微元个数%竖直方向分割%水平方向的对应坐标%对应点的光程差%对应点的相位差%光强与相位差相联系%以下为形成图像过程，将光强与图像辉度相对应程序代码如下：ClearLambda=500*(1e-9);theta=1*(1e-6);ni=500;ds=linspace(0,0.0005,ni);fork=1:niy(k)=ds(k)/sin(theta);Delta=2*ds(k)+Lambda/2;Phi=2*pi*Delta/Lambda;B(k,:)=4*cos(Phi/2).入

2、2;endfigure(gcf);NCLevels=250;Br=(B/4.0)*NCLevels;image(0,y,Br);colormap(gray(NCLevels);title(二维强度分布);(2)(4)程序代码如下:clearLambda=500*(1e-9);r=5;t=3;ni=1000;y=linspace(-0.005,0.005,ni);fork=1:nih(k)=t-(r-sqrt(r入2-y(k)入2);Delta=2*h(k)+Lambda/2;Phi=2*pi*Delta/Lambda;B(k,:)=4*cos(Phi/2).入2;endfigure(gcf);

3、NCLevels=250;Br=(B/4.0)*NCLevels;image(0,y,Br);colormap(gray(NCLevels);title(二维强度分布);%清除原有变量%设置波长%设置介质的尺寸%分割份数%彷向上进行分割%对应点的介质厚度%对应点的光程差%对应点的相位差%对应点的光强%以下为形成图像过程，将光强与图像辉度相对应(3)程序代码如下：clearlambda=500*1e-9;R=5;H=2;ni=2000;x=linspace(-0.005,0.005,ni);y=linspace(-0.005,0.005,ni);fori=1:niforj=1:nir=x(i).

4、入2+y(j).入2;delta=2*H+2*(R-sqrt(R入2-r)+lambda/2;phi=2*pi*delta/lambda;B(i,j)=4*cos(phi/2).入2;endendfigure(gcf);NCLevels=250;Br=(B/4.0)*NCLevels;image(x,y,Br);colormap(gray(NCLevels);title(二维强度分布);%清除原有变量设置波长和介质尺寸%设置分割份数%在乂，y方向上进行分割%点距介质中心的直线距离相应点的光程差相应点的相位差%相应点的光强%以下为形成图像过程，将光强与图像辉度相对应x10程序代码如下：clearlambda=500*1e-9;R=5;H=2;ni=2000;T=0.005;omega=2*pi/T;%a2=axes(Position,0,1,0,1);%图片在figure中的位置前两个代表x方向的坐标，后两个是y方向的坐标x=linspace(-T,T,ni);y=linspace(-T,T,ni);fori=1:niforj=1:nir=sqrt(x(i).入2+y(j).入2);delta=2*H+sin(omega*r)*1e-6;phi=2*pi*delta/lambda;B(i,j)=4*cos(phi/2).入2;endendfigure(gcf);NCLevel