请你给我一段菲涅尔圆孔衍射的MATLAB代码

菲涅尔圆孔衍射是一种光学现象，通常用于展示波前通过小孔时的分布变化。在MATLAB中，我们可以使用复数数组和傅里叶变换来模拟这个过程。以下是一个简单的示例代码，它创建了一个圆形光源并通过一个圆孔进行衍射：

% 定义参数
lambda = 500; % 波长，这里以纳米为单位
diameter = 1e-4; % 圆孔直径，同样以米为单位
radius = diameter / 2;
grid_size = 200; % 显示区域的像素大小
dx = 1 / grid_size; % 每个像素的宽度

% 创建空间坐标网格
[x, y] = meshgrid(-grid_size*dx:dx:grid_size*dx);

% 创建圆形光源函数
source_func = @(x, y) exp(1i * 2*pi/lambda * x); % 简化模型，假设为平面波

% 创建圆孔函数
aperture_func = ones(grid_size, grid_size);
aperture_func(sqrt(x.^2 + y.^2) < radius) = 0;

% 对光源进行傅里叶变换，然后与圆孔函数相乘
transformed_func = ifft2(source_func .* fftshift(aperture_func));

% 计算并显示衍射图案
disp_pattern = abs(transformed_func).^2;
imagesc(x, y, disp_pattern);
colormap(gray);
colorbar;
title('Fresnel圆孔衍射图案');

%