怎样用FDTD仿真整个超透镜？几何相位、传输相位

这个问题太多人问我了，有好多问题我自己也没解决，我只会最基础的，在这里简单写一下。
透射超透镜基本就两种
1.几何相位
这个比较简单，只需要直接在FDTD里面编程就行，我直接把程序粘这了，拿来改改就能用。（这是官网案例不是我写的，贴上去只是因为太多人找我要了）
#TiO2 metalens
#Paper Metalenses at visible wavelengths:Diffraction-limited focusing andsubwavelength resolution imaging

selectall;
delete;
nm=1e-9;
um=1e-6;
W=95nm; #微柱宽
L=250nm; #微柱长
H=600nm; #微柱高
S=325nm; #微柱周期（间隔）
Lambda=532nm; #工作波长
f=1um; #超透镜焦距（测试用）
#f=90um; # 文献中焦距
R=3um; # 超透镜半径
N=10;M=10; # 横向N、纵向M 微柱个数

#循环生成微柱、根据微柱坐标旋转结构。再用判断语句，将结构限制为圆形
for (i=-N:N){
for(j=-M:M){
x=iS; y=jS;

Circle=(x^2+y2<R^2);
#if(Circle){ #判断语句用于将结构限制为圆形
Theta=-180*(f-sqrt(x