杨氏双缝干涉matlab代码

以下是一个简单的杨氏双缝干涉matlab代码：

% 定义常量
lambda = 632.8e-9; % 波长
D = 1e-3; % 双缝间距
L = 1; % 屏幕到双缝距离
w = 0.2e-3; % 单缝宽度
N = 1000; % 屏幕像素数
dx = 10e-6; % 屏幕像素宽度

% 定义坐标系
x = linspace(-N/2,N/2-1,N)*dx;
[X,Y] = meshgrid(x);

% 计算双缝的干涉模式
k = 2*pi/lambda;
theta = atan(X/L);
sin_theta = sin(theta);
I0 = (w*sin(k*D*sin_theta/2)./(k*D*sin_theta/2)).^2;
I0(theta==0) = 1; % 避免除零错误

% 显示干涉图案
figure;
imagesc(x,x,I0);
xlabel('x (m)');
ylabel('y (m)');
colormap('gray');
colorbar;

代码中首先定义了一些常量，包括波长、双缝间距、屏幕到双缝距离、单缝宽度、屏幕像素数和屏幕像素宽度。然后定义了坐标系，并计算了双缝的干涉模式。最后使用`imagesc`函数显示干涉图案。