MATLAB仿真光的干涉（四）

本文以牛顿环为代表介绍光的等厚干涉。

和等倾干涉不同，等厚干涉中平行平面板的两个表面不严格平行。

一个很著名的例子是牛顿环，考虑垂直入射的情况：

光在透镜的上表面（平面）和玻璃平面分别反射，形成干涉。

考虑到从光疏介质到光密介质的反射涉及到半波损失（详见本专栏的光的反射相关文章），两光束的相位差为：半波损失导致的+往返透镜和空气的光程差*

详细的项目代码如下：

R = 1e1;%透镜半径
Lambda = 0.55*1e-6;
x = linspace(-0.01,0.01,1001);
y = x;
[X Y] = meshgrid(x,y);
r = sqrt(X.^2+Y.^2);
I = cos(pi*(r.^2/R+Lambda/2)/Lambda).^2%I = I0^2*(cos(δ/2))^2
NCLevels = 255;%指定调色板
colormap(gray(NCLevels));
Ir = I*NCLevels;
image(x,y,Ir);
title('牛顿环（单位：）m');
axis square;
figure
r(501,:)
I(501,:)
plot(r(501,:),I(501,:),'k')%相当于两条线重合
title('干涉光强分布图');
xlabel('r/（环半径）m');
ylabel('相对光强');

结果：

可以看出：

中心为暗斑，越向外环的宽度越细、越密集。 

从干涉光强分布图中可以更清晰的看出这一点。

参考《高等光学仿真——光波导、激光》