matlab光学仿真——楔板

《matlab光学仿真——楔板》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、5、用平行光（点光源+准直镜）照射如图7所示不同形状楔板，观察干涉图样（1）程序代码如下：Clear%清除原有变量Lambda=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）程序代码如下：clear%清除原有变量Lambda=500*(1e-9);%设置波长r=5;%设置介质的尺寸t=3;ni=1000;%分割份数y=linspace(-0.005,0.005,ni);%y方向上进行分割fork=1:nih(k)=t-(r-sqrt(r

3、^2-y(k)^2));%对应点的介质厚度Delta=2*h(k)+Lambda/2;%对应点的光程差Phi=2*pi*Delta/Lambda;%对应点的相位差B(k,:)=4*cos(Phi/2).^2;%对应点的光强endfigure(gcf);%以下为形成图像过程，将光强与图像辉度相对应NCLevels=250;Br=(B/4.0)*NCLevels;image(0,y,Br);colormap(gray(NCLevels));title('二维强度分布');（3）程序代码如下：clear%清除原有变量lambda=500*1

4、e-9;R=5;H=2;%设置波长和介质尺寸ni=2000;%设置分割份数x=linspace(-0.005,0.005,ni);%在x，y方向上进行分割y=linspace(-0.005,0.005,ni);fori=1:niforj=1:nir=x(i).^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

5、(gcf);%以下为形成图像过程，将光强与图像辉度相对应NCLevels=250;Br=(B/4.0)*NCLevels;image(x,y,Br);colormap(gray(NCLevels));title('二维强度分布');（4）程序代码如下：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)

6、;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);NCLevels=250;Br=(B/4.0)*NCLevels;image(x,y,Br);colormap(gray(NCLevels));title('二维强度分布');