我要一段几何相位超透镜相位分布的MATLAB代码

下面是一段简单的 MATLAB 代码，模拟了几何相位超透镜的相位分布：

clear;
clc;

lambda =532e-9; % 波长(单位：米)
f = 50e-3; % 焦距(单位：米)
d = 5e-6; % 光栅周期(单位：米)
z = 1e-3; % 光栅到超透镜的距离(单位：米)

[X, Y] = meshgrid(-2e-3:5e-6:2e-3, -2e-3:5e-6:2e-3); % 创建网格

k = 2 * pi / lambda; % 波数

phase = k * d * (sin(pi * X / d) .^ 2 + sin(pi * Y / d)