（等倾、等厚干涉）MATLAB在迈克尔逊干涉仪中的应用

理论回顾：

迈克尔逊干涉仪，是1881年美国物理学家迈克尔逊和莫雷合作，为研究“以太”漂移而设计制造出来的精密光学仪器。它是利用分振幅法产生双光束以实现干涉。通过调整该干涉仪，可以产生等厚干涉条纹，也可以产生等倾干涉条纹。主要用于长度和折射率的测量，若观察到干涉条纹移动一条，便是 $M_2$的动臂移动量为 $\frac{\lambda }{2}$，等效于$M_1$与$M_2$ 之间的空气膜厚度改变 $\frac{\lambda }{2}$。
干涉光强合成公式：根据波的叠加理论，两束同频率单色光在空间某一点光矢量的大小为 $E_1=A_1\cos (\omega t+{\phi }_1)$和$E_2=A_2\cos (\omega t+{\phi }_2)$，$A_1$、$A_2$ 分别是两个光矢量的振幅,${\phi }_1$、${\phi }_2$分别是初相。若两个光矢量的方向相同，合成光矢量为 ，振幅和初相分别为$$\begin{array}{l}{A}_0=\sqrt{{A_1}^2+{A_2}^2+2A_1{A}_2\cos ({\phi }_2-{\phi }_1)}\\ {\phi }_0=\arctan \frac{A_1\sin {\phi }_1+A_2\sin {\phi }_2}{A_1\cos {\phi }_1+A_2\cos {\phi }_2}\end{array}$$在一定时间内观察到的平均光强 $I$与光矢量$A_0$ 的平方的平均值成正比,$$I=k{A_0}^2=k\left[{A_1}^2+{A_2}^2+2A_1{A}_2\cos ({\phi }_2-{\phi }_1)\right]$$

$k$是比例系数, $$I=I_1+I_2+2\sqrt{I_1{I}_2}\cos \Delta \phi$$ $$\Delta \phi ={\phi }_2-{\phi }_1$$ $$I_1=I_2=I_0$$ $$I=2I_0(1+\cos \Delta \phi )=4I_0{\cos }^2\frac{\Delta \phi }{2}$$ $\Delta \phi =\frac{2\pi }{\lambda }\delta$ $$I=4I_0{\cos }^2(\frac{\pi }{\lambda }\delta )\delta =r_2-r_1$$

计算模型：

迈克尔逊干涉仪的实验原理图如下：

任一圆环半径$R=x^2+y^2$其中，在屏幕前加上焦距为f的透镜，则
$\tan \theta =\frac{\sqrt{x^2+y^2}}{f}$，$\theta =\arctan \frac{\sqrt{x^2+y^2}}{f}$
$\delta =2dcos\theta =2dcos(\arctan (\frac{\sqrt{x^2+y^2}}{f}))$
当$M_1$与${M_2}^{\prime }$之间存在一定角度时：

根据这些关系式编写出计算屏幕上各点的光强程序：
打开Matlab的m文件编辑器，在里面输入如下程序代码:

Lambda = input('输入光的波长（单位为nm）：');
Lambda = Lambda * 1e-9; 
d = input('输入M1和M2’之间的距离（单位为mm）：');
d = d * 1e-3 ; 
f = 0.1; %透镜的焦距：m
I0 = 1;
temp=input('请输入M1与M2’之间的夹角(单位：pi) :');
if temp==0
    Screen_length = 0.015;
else
    Screen_length = 0.00015;  %定义干涉仿真范围：m
end
[Screen_x,Screen_y]=meshgrid(linspace(-Screen_length,Screen_length,800));
interference_r = abs(Screen_x +1i*Screen_y);
g= abs(Screen_x);
I_delta=(2*pi*(d+g.*tan(temp))/Lambda).*cos(atan(interference_r./f));
I = 4*I0*(cos(I_delta)).^2;
I = I./max(max(I));                      %光强分布归一化
I = I*255;                                %光强归一,扩大显示
image(Screen_x(1,:),Screen_y(:,1),I);    %设置x和y的像素，显示数值
colormap(gray(255));                            %也可以采用gray
colorbar;
xlabel('空间坐标x'),ylabel('空间坐标y');
title('迈克尔逊干涉光强空间分布');

输入完成后，按 “保存”按钮，选择合适的名字，并以“.m”为扩展名保存该文件。存盘完成后，我们就可以在命令窗口输入该文件的名字运行这段代码，选择如下参数：

输入光的波长（单位为nm）：589.3
输入M1和M2’之间的距离（单位为mm）：0.23
请输入M1与M2’之间的夹角(单位：pi) :0

结果显示：

输入：

输入光的波长（单位为nm）：589.3
输入M1和M2’之间的距离（单位为mm）：0.23
请输入M1与M2’之间的夹角(单位：pi) :pi/100

结果显示：