使用matlab模拟仿真迈克尔逊干涉仪等厚干涉

clear; d=2;%透镜厚度 lamda=632.8e-006;%氦氖激光器 theta=0.25;%初始入射角 rmax=d*tan(theta/2); N=800; x=linspace(-rmax,rmax,N); y=linspace(-rmax,rmax,N);%设置坐标轴范围 for k=0:15 theta=0.25-0.01*k; %模拟倾角的增加 rmax=d*tan(theta/2); for i=1:N for j=1:N x(i)=(i-2)*2*rmax/(N-1)-rmax; y(j)=(j-2)*2*rmax/(N-1)-rmax; r(i,j)=sqrt(x(i)^2+y(j)^2); delta(i,j)=2*d/sqrt(1+r(i,j)^2/d^2); Phi(i,j)=2*pi*delta(i,j)/lamda; B(i,j)=4*cos(Phi(i,j)/2)^2;%条纹位置 end end figure(gcf); N1=255; Br=2.5*B*N1; image(x,y,Br); colormap(gray(N1)); set(gca,'XTick',[]); set(gca,'YTick',[]); axis equal drawnow pause end