本模拟通过在单色光的双光束干涉的图样上进行更改频率后多次叠加模拟了干涉光的谱宽
lam=520e-9;
theta=pi/2;
aph=pi/2;
A1=1;
a=1;
A2=a*A1;
xm=0.000002;
ym=xm;
n=1001;%采样率
xs=linspace(-xm,xm,n);
ys=linspace(-ym,ym,n);
[xs,ys]=meshgrid(xs,ys);
E=0;
for lam=520e-9:1e-9:820e-9%更改波长范围
E1=A1.*exp(-1i*xs*cos(aph)*2*pi/lam);
E2=A2.*exp(1i*xs*cos(theta-aph)*2*pi/lam);
E=E1+E2+E; %反复叠加
end
I=abs(E)^2;
figure(1)
pcolor(xs,ys,I);
shading flat;
colormap gray
figure(2)
xs1=xs(1,:);
Iplot=I/max(max(I));
Iplot1=Iplot(1,:);
plot(xs1,Iplot1);
xlabel('空间坐标x'), ylabel('合成光强度');
title('合成光强空间分布函数');
两光束干涉,改变干涉光光谱宽度,干涉光最小为520nm,取样率为1nm。
520nm单色光时