matlab 编写夫琅禾费双缝衍射

夫琅禾费双缝衍射是光学领域中的一个经典问题，可以用 MATLAB 进行模拟。下面是一份简单的代码：

% 定义常量
lambda = 632.8e-9;  % 波长
d = 0.25e-3;  % 双缝间距
D = 1;  % 双缝到屏幕距离
L = 10e-2;  % 屏幕宽度
N = 2048;  % 屏幕像素数
dx = L / N;  % 像素间距

% 定义双缝函数
u = zeros(1, N);
u(N/2-dx/2:N/2+dx/2) = 1;
u(N/2+d+dx/2:N/2+d-dx/2+N) = 1;

% 计算衍射图样
x = linspace(-L/2, L/2, N);
[X, Y] = meshgrid(x);
r1 = sqrt((X-d/2).^2 + Y.^2 + D^2);
r2 = sqrt((X+d/2).^2 + Y.^2 + D^2);
I = (u).^2 .* (sin(pi*d/lambda*r1)./(pi*d/lambda*r1)).^2 ...
    .* (sin(pi*d/lambda*r2)./(pi*d/lambda*r2)).^2;

% 显示结果
figure;
imagesc(x, x, I);
axis equal tight;
colormap(gray);
xlabel('x (m)');
ylabel('y (m)');

代码中，我们首先定义了夫琅禾费双缝衍射的常量，包括波长、双缝间距、双缝到屏幕距离、屏幕宽度、像素数和像素间距等。接着定义了一个双缝函数 `u`，其中双缝的位置通过数组下标来表示。然后通过计算双缝函数的傅里叶变换，得到衍射图样 `I`，并用 `imagesc` 函数将其显示出来。

运行代码后，可以得到夫琅禾费双缝衍射的衍射图样，如下所示：