matlab仿真杨氏双缝干涉实验的代码

以下是 MATLAB 实现的杨氏双缝干涉实验的代码：

% 设置参数
lambda = 633e-9; % 波长
d = 0.1e-3; % 双缝间距
D = 1; % 屏幕到双缝距离
a = 0.1e-3; % 缝宽
L = 0.5; % 屏幕大小
N = 500; % 屏幕像素数

% 计算干涉图样
x = linspace(-L/2, L/2, N);
[X, Y] = meshgrid(x, x);
r1 = sqrt((X-d/2).^2 + Y.^2);
r2 = sqrt((X+d/2).^2 + Y.^2);
I = (cos(pi*a*r1/lambda/D).^2 + cos(pi*a*r2/lambda/D).^2 + 2*cos(pi*a*r1/lambda/D).*cos(pi*a*r2/lambda/D)) / 2;
I = I / max(I(:));

% 显示干涉图样
figure;
imagesc(x, x, I);
colormap('gray');
axis('equal', 'tight');
title('杨氏双缝干涉图样');
xlabel('位置 (m)');
ylabel('位置 (m)');

解释：

- `lambda`：光的波长。
- `d`：双缝间距。
- `D`：屏幕到双缝距离。
- `a`：缝宽。
- `L`：屏幕大小。
- `N`：屏幕像素数。
- `x`：屏幕上的位置向量。
- `X` 和 `Y`：屏幕上的位置网格。
- `r1` 和 `r2`：屏幕上每个像素到两个缝的距离。
- `I`：每个像素的干涉强度。
- `imagesc`：用灰度图显示干涉图样。

你可以根据需要修改参数，运行代码即可得到杨氏双缝干涉的图样。