实现一个MATLAB函数来绘制Airy disk： 该函数应以数值孔径和光波长作为两个输入变量。为了简化，将其峰值大小归一化为1。那么如何编写代码

以下是一个MATLAB函数来绘制Airy disk的示例代码：

function AiryDisk(numAper, waveLen)
% numAper: 数值孔径
% waveLen: 光波长

% 定义常数
k = 2 * pi / waveLen; % 波矢量
NA = numAper; % 数值孔径
f = NA / waveLen; % 焦距
r = 0 : 0.001 : 1; % 半径范围

% 计算Airy disk强度分布
airyDisk = ((2 * besselj(1, k * NA * r) ./ (k * NA * r)) .^ 2);
airyDisk(r == 0) = 1; % 避免除以0错误

% 绘制Airy disk图像
figure
plot(r, airyDisk, 'LineWidth', 2)
xlabel('Radius (\lambda / NA)')
ylabel('Normalized Intensity')
title('Airy Disk')
end

该函数首先计算一些常量，例如波矢量，数值孔径和焦距。然后，使用`besselj`函数计算Airy disk的强度分布，并将其绘制为图像。在绘制过程中，我们将半径范围定义为0到1，并在避免除以0错误的情况下对Airy disk的中心进行了归一化。

可以使用以下命令调用函数：

AiryDisk(0.5, 633e-9)

其中，输入变量为数值孔径为0.5，光波长为633纳米的情况。