研究光源宽度对干涉条纹可见度的影响（matlab仿真）

一、干涉条纹的可见度

对于光波来说，干涉现象往往表现为明暗相间的条纹。为了描述干涉图场中的强弱对比，引入可见度的概念，其定义为：

当Imin=0,K=1.0时,条纹最清晰；当Imax=Imin时，K=0，条纹消失。

影响干涉条纹可见度的因素主要有三个：两相干光的振幅比、光源的线宽、光源的单色性。

二、光源的宽度对干涉条纹的影响

一般实验中对干涉条纹性质的分析都是以点、缝或者细线光源为前提，而实际上的光源总是具有一定宽度的。

以杨氏双缝实验为例（见图1）:

图1

设入射光源的宽度为b，于是宽度为b的整个光源在x平面P点处的光强为：

 其中，K为干涉条纹可见度，为：

条纹可见度K和光源宽度b的关系如图2所示：

图2 

三、matlab仿真结果

运用上述数学关系式，在matlab中进行光源具有一定的宽度的双缝干涉实验仿真。设置基本仿真参数为：I0=550nm，d=0.1mm，D=1m，L=80mm。

（1）当K=0.9，即b=0.11mm时

仿真结果：

图3 

结果分析：K=0.9时，能够较清晰地观察到干涉条纹，光强呈余弦规律变化。

（2）当K=0，即b=0.44mm时

仿真结果：

图4

结果分析：K=0时，条纹可见度为0，光强I=2I0。

通过改变光源宽度b的值，可以得到不同的条纹可见度K。从图5可以分析得到，在临界宽度（b=0.44mm）到允许宽度（b=0.11mm）之间，干涉条纹可见度随着光源宽度的增加而减小。

图5   不同可见度K的干涉图样

四、matlab代码

%***********************杨氏双缝干涉实验仿真*************************
clear,clc;
% 基本参数输入及处理*********************************
Lam= input('请输入光的波长lam（单位为nm）：');
Lam = Lam * 1e-9;       %波长单位转为：m
d = input('请输入两个缝的间距d（单位为mm）：');
d = d * 1e-3;                 %缝间隙单位转为：m
b = input('请输入光源的宽度b（单位为mm）：');
b = b * 1e-3;                 %缝间隙单位转为：m
L = input('请输入光源到缝的距离L（单位为mm）：');
L = L * 1e-3;                 %缝间隙单位转为：m
D = input('请输入缝到屏的距离D（单位为m）：');
I0 = 1;                       %两列波的光强

Xmax = 0.01;          
x = linspace(-Xmax,Xmax,101); %坐标取样101份
I = zeros(101,1);
for i=1:101
    r1 = sqrt((x(i)+d/2)^2+D^2);    %在x坐标轴上的情况
    r2 = sqrt((x(i)-d/2)^2+D^2);
    phi = 2*pi*(r2-r1)/Lam;      %两列波在该点的相位差
    K=(sin(pi*b*d/L/Lam))./(pi*b*d/Lam/L);%可见度k
    I(i,:) = 2*I0*(1+K.*cos(phi));     %合成波光强
end

% 双缝干涉条纹图样*************************
I_change = I./4 .*255;               %光强归一化 并转为255灰度值
subplot(2,1,1);
image(x,Xmax,I_change');             %设置x和y的像素，显示数值
colormap(gray(255));                 %妆点灰度色彩
xlabel('空间坐标x');
title('双缝干涉条纹图样');

% 双缝干涉光强空间分布函数图*************************
subplot(2,1,2);
plot(x,I');
axis([-Xmax,Xmax,0,4]);
xlabel('空间坐标x'), ylabel('合成光强度');
title('双缝干涉光强空间分布函数');