衍射微透镜 设计 matlab,亚波长衍射微透镜色散的数值分析

亚波长衍射微透镜色散的数值分析3刘玉玲　卢振武 (中国科学院长春光学精密机械与物理研究所 ,应用光学国家重点实验室 , 长春 　130022) (2003 年 8 月 29 日收到;2003 年 10 月 31 日收到修改稿) 　　采用旋转体时域有限差分法对亚波长衍射微透镜进行数值分析 ,给出了在微透镜焦平面上的电场强度分布和从近场到远场的电场强度分布. 研究了亚波长衍射微透镜的色散 ,给出了入射波长不同于设计波长的色散曲线 , 数值分析结果表明了亚波长衍射微透镜的负向色散特性 ,而且随着波长的减小焦距增加得更快. 关键词 : 亚波长衍射微透镜 , 色散 , 时域有限差分法 PACC: 4225F 3中国科学院长春光学精密机械与物理研究所创新项目(批准号:2002LQ. 04)资助的课题. 11引 言 亚波长衍射微透镜具有体积小、重量轻、集成度高、易于复制等特点被广泛应用于图像识别和处理、光学传感、光通讯、激光医学、空间光学等许多领域. 对于亚波长尺寸的元件 ,标量衍射理论已经不能对其进行精确分析[1 ,2] ,必须采用严格的矢量分析方法. 然而目前的文献大多只是对周期结构的光栅进行矢量衍射分析[3] ,而关于分析非周期结构的亚波长衍射微透镜的文献却非常少. 本文把旋转体时域有限差分法(BOR FDTD) [4] 作为分析方法 ,说明了分析过程、理论公式和 Matlab 程序流程 ,以此为工具对亚波长衍射微透镜进行数值分析 ,给出了在微透镜焦平面上的电场强度分布和从近场到远场的电场强度分布. 研究了亚波长衍射微透镜的色散 ,给出了入射波长不同于设计波长的色散曲线 ,数值分析结果表明了亚波长衍射微透镜的负向色散特性 ,即随着入射波长的减小焦距增加 ,而且随着波长的减小焦距增加得更快. 21分析方法 本文采用 BOR FDTD 来分析旋转对称的亚波长衍射微透镜. 使用该方法的优点是不仅能对旋转对称的微透镜进行精确的模拟 ,而且能将三维问题转 化为二维问题并且只需分析半个微透镜结构 ,因此大大简化了计算 ,节约了时间. 图 1 　BOR FDTD 分析微透镜的几何框图 图 1 为采用 BOR FDTD 算法分析亚波长衍射微透镜的几何框图 ,图中 Z 轴为旋转对称轴. BOR FDTD 算法利用 Yee 氏网格对计算空间进行离散 ,被分析的元件位于网格空间中 ,通过赋予网格点电磁参数的方法来模拟元件的电磁结构. 在 BOR FDTD 计算空间的外边界设置了具有一定厚度的完全匹配层(PML) [5] 用来截断计算空间 ,在 PML 内采用完全匹配吸收边界条件(PML ABC)算法来消除非物理因素引起的电磁波的反射 ,从而模拟电磁波在开放空间的传播 ,这样既可以大大减轻计算负担又可以节约计算时间. 为了便于设置连续平面波源 ,把 BOR FDTD 计算空间划分为总场和散射场 ,被分析的元件位于总场内 ,如图 1 所示. 总场内既存在入射波又存在散射波 ,在散射场内只考虑散射波 ,不计入射波 ,总场 第 53 卷 第 6 期 2004 年 6 月 100023290200453(06)Π1782206 物 　理 　学 　报 ACTA PHYSICA SINICA Vol. 53 ,No. 6 ,June ,2004 ν2004 Chin. Phys. Soc. 与散射场的交界处为连接边界 ,垂直入射的平面波便由连接边界引入总场. 在初始时刻整个计算网格空间的电磁场均设为零 ,而后平面波源被接通 ,随着时间步的增加 ,入射平面波由连接边界向总场区内传播 ,当碰到