clear,clc;
%波长（典型值：600）
Lambda = input('输入光的波长（单位为nm）：');
Lambda = Lambda * 1e-9;
%缝间隙（典型值：1）
d = input('输入两个缝的间距（单位为mm）：');
d = d * 1e-3;
%缝到屏的距离（典型值：1）
D = input('输入缝到屏的距离（单位为m）：');
%设置空间坐标的边界
Xmax = 3 * Lambda * D/d; 
Ymax = Xmax;
Tmax = Xmax;
%坐标取样1001份
x = linspace(-Xmax, Xmax, 1001);      
y = linspace(-Ymax, Ymax, 1001);
t = linspace(-Tmax, Tmax, 1001);
%生成二维网格
[x, y] = meshgrid(x, y);
%计算复振幅
r1 = sqrt((x+d/2).^2+D^2);   
r2 = sqrt((x-d/2).^2+D^2);
phi = 2*pi*(r2-r1)/Lambda;  
E1 = 1;
E2 = 1 * exp(-1i * phi);
E = E1 + E2;
%合成波光强
I = abs(E) .^ 2; 
I = I ./ 4;
%取y轴的中线处的光强
I0 = I(501, :);     
%合成光强灰度分布图
subplot(2,1,1);     
imshow(I);            
xlabel('空间坐标x');
ylabel('空间坐标y');
title('杨氏双缝干涉合成光强空间分布');
% 合成光强函数分布图
subplot(2,1,2);
%设置坐标轴
axis([-Tmax,Tmax,0,1]);
plot(t, I0);
xlabel('空间坐标x'), ylabel('合成光强度');
title('杨氏双缝干涉合成光强空间分布函数');