基于FDTD方法研究一维孔径菲涅耳衍射.pdfVIP

第 26卷第 6期 泉州师范学院学报 (自然科学) Vo1．26No．6 2008年 11月 JournalofQuanzhouNormalUniversity(NaturalScience) Nov．2008 基于 FDTD方法研究一维孔径菲涅耳衍射 余燕忠 ，林 梅。 (1I东南大学 毫米波国家重点实验室，江苏 南京 210096I2．泉州师范学院 理工学院，福建 泉州 362000；3．江西工业工程职业技术学院，江西 萍乡 337055) 摘 要：应用时域有限差分FDTD方法研究一维孔径菲涅耳衍射花样的基本特征．首先简介了FDTD的 基本原理和实现方法，接着从理论上分析了一维孔径菲涅耳衍射花样的规律，最后给出应用 FDTD方法的数 值仿真结果，该结果与理论分析相吻合． 关键词：一维孔径 ；菲涅耳衍射 ；衍射花样 IFDTD 中图分类号 ：O414．4 文献标识码 ：A 文章编号 ：1009—8224(20o8)06—0037—04 时域有限差分FDTD(FiniteDifferenceTimeDomain)方法是K．S．Yee于 1966年首次提出来的[1]，该方 法 自提出来就得到迅速发展及广泛应用口]．其主要思想是以Yee元胞为空间电磁场离散单元，对电磁场分量 在空间和时间上采取交替抽样的离散方式，将含时间变量的Maxwell方程转化为差分方程，并在时间轴上逐 步推进地求解空间电磁场．因此，FDTD方法非常适合于动态模拟电磁波的实际传播过程，是一种有效的电磁 数值计算与分析的工具． 光的衍射是光波的基本特征，所以对衍射现象和规律的研究是很有必要的．在光学中，菲涅耳衍射 占有重 要的地位，本文只研究一维孔径的菲涅耳衍射．当然 ，该方法也可以用于研究其他的衍射现象和规律．由于菲 涅耳衍射积分的被积函数中包含一个与衍射孔径坐标有关的二次相位因子，因此，即使对于简单的衍射孔径 (如单缝、直边、圆孔等)，也难以得出解析结果0]．基于这一点考虑，本文采用数值分析的方法来研究其衍射现 象和规律，并应用 FDTD方法研究一维孔径菲涅耳衍射的基本特征． 1 FDTD基本方法简介[2,4] 本文是在二维TM 波情形下仿真菲涅耳衍射过程，故以下只介绍二维TM 波的FDTD方法． 1．1 直角坐标中的FDTD公式 在亘角坐标中 ，二fitTM 泼方程为 fI一一等百--a．,H， {等一警+ (1) a【【x一8y—e8t+…帆… 其中，E为电场强度 ，H为磁场强度 ，￡表示介电系数，表示磁导系数，口表示 电导率， 表示导磁率． 利用Yee元胞对式 (1)进行 FDTD差分离散 (如图 1所示)，可得到 FDTD公式为 图 l 二维TM波的Yee元胞 收稿 日期 ：2008一O6—16 作者简介：余燕忠 (1972一 )，男，福建莆田人，讲师，博士生。主要从事 电磁场数值计算及仿真研究． 基金项 目：福建省tl然科学基金计划资助项 目(A0610027) 38 泉州师范学院学报 (自然科学) 20O8年 11月 H (，+丢)一CP(m)H7(f，+1)一OQ(m)堡 ± ， (2) ，2(汁 1 一CP(m)H~- (件 1 +OQ()墅 盟 ， (3) 舛 (H一÷，)一 (i-÷，J) 门(，J+1)一 (f，一÷) E~．-1(i，)= ( E：( )+(】3( [————_厶__Z —— 一———— —— ]·(4)