多缝夫琅禾费衍射matlab代码

夫琅禾费衍射是一种光学现象，它是指平行光线经过一个宽度很宽的孔或者一组孔后，产生的衍射图样。下面是一个简单的MATLAB代码实现多缝夫琅禾费衍射：

%定义常量
wavelength = 0.5e-6; %光的波长
slit_width = 10e-6; %缝宽
slit_distance = 100e-6; %缝距
screen_distance = 1; %屏幕距离
screen_size = 0.1; %屏幕大小

%定义一个屏幕
screen_points = 1000;
x = linspace(-screen_size/2, screen_size/2, screen_points);
y = linspace(-screen_size/2, screen_size/2, screen_points);
[X, Y] = meshgrid(x, y);

%计算每个点的衍射强度
intensity = zeros(screen_points);
for i = 1:screen_points
    for j = 1:screen_points
        %计算光程差
        r1 = sqrt((X(i,j)-slit_distance/2)^2 + Y(i,j)^2);
        r2 = sqrt((X(i,j)+slit_distance/2)^2 + Y(i,j)^2);
        path_diff = r2 - r1;
        %计算衍射强度
        intensity(i,j) = (sinc(slit_width*wavelength*path_diff/(2*pi)))^2;
    end
end

%绘制衍射图像
figure;
imshow(intensity, []);

这段代码会生成一幅多缝夫琅禾费衍射图像。如果需要绘制单缝夫琅禾费衍射，只需要将上述代码中的`slit_distance`设为无穷大即可。

需要注意的是，上述代码中使用了`sinc`函数。该函数是MATLAB中的一个内置函数，它计算的是`sinc(x) = sin(pi*x)/(pi*x)`的值。如果你使用的是旧版本的MATLAB，可能会没有这个函数，需要自己编写。