探地雷达FDTD非均匀网格二维仿真

我有一个三段线性变化的介电常数模型，为了平衡其效率和精度的影响，拟采用变化空间步长的的方式来实现二维仿真，需求是得到网格的剖分方式，精度与效率的结果，有基本均匀网格仿真的代码！