毕业设计（论文）-基于FDTD的E和H型源激励的矩形波导的截止频率和动态特性研究.doc

摘 要 时域有限差分法是由有限差分法发展出来的数值计算方法。自1966年Yee在其论文中首次提出时域有限差分以来，时域有限差分法在电磁研究领域得到了广泛的应用。主要有分析辐射条线、微波器件和导行波结构的研究、散射和雷达截面计算、分析周期结构、电子封装和电磁兼容的分析、核电磁脉冲的传播和散射以及在地面的反射及对电缆传输线的干扰、微光学元器件中光的传播和衍射特性等等。 由于电磁场是以场的形态存在的物质，具有独特的研究方法，采取重叠的研究方法是其重要的特点，即只有理论分析、测量、计算机模拟的结果相互佐证，才可以认为是获得了正确可信的结论。时域有限差分法就是实现直接对电磁工程问题进行计算机模拟的基本方法。在近年的研究电磁问题中，许多学者对时域脉冲源的传播和响应进行了大量的研究，主要是描述物体在瞬态电磁源作用下的理论。另外，对于物体的电特性，理论上具有几乎所有的频率成分，但实际上，只有有限的频带内的频率成分在区主要作用。 文中主要分析了源处不同性质的激励源对波导截止频率特性的影响，得到了在二维情况下对于矩形波导的场分量瞬态图以及截止频率特性图。详细模拟了TE和TM模式场结构分布以及脊波导中的电磁波传播的场分量瞬态分布。 关键词：时域有限差分 波导 截止频率 场结构特性 目录 摘 要 4 目录 5 第一章 绪论 6 1.1 课题背景与意义 6 1.2 时域有限差分法的发展与应用 7 1.3 时域有限差分法的特点 9 1.4 本课题的主要研究内容 12 第二章 时域有限差分法的基本原理 13 2.1 Maxwell方程和Yee氏算法 13 2.2 FDTD的基本差分方程 15 2.3时域有限差分法相关技术 17 2.3.1数值稳定性问题 17 2.3.2数值色散 18 2.3.3 FDTD离散网格的确定 19 2.4吸收边界条件 19 2.4.1一阶和二阶近似吸收边界条件 20 2.4.2 二维棱边及角顶点的处理 23 2.4.3三维棱边及角顶点的处理 26 2.4.3 Berenger完全匹配层 26 2.5激励源设计 30 2.5.1激励源类型 30 2.5.2激励源的设置 31 2.6 FDTD计算所需时间步的估计 32 第三章 FDTD法分析矩形波导的截止频率 34 3.1 波导理论 34 3.2 傅立叶变换算法的特点 34 3.3 矩形波导的截止频率 35 第四章 波导的场结构特性 43 4.1波导的场结构特性 43 4.2脊波导的场结构特性 48 第五章 结论 50 参考文献 51 攻读学位期间的研究成果 55 第一章 绪论 课题背景与意义 1873年Maxwell建立电磁场基本方程以来，电磁场理论的应用和发展已经有一百多年的历史了。现代技术的许多方面都与电磁场，尤其是与高频电磁场有关，复杂的高频电磁场系统的分析和综合，以及高频电磁场与复杂目标相互作用的分析和计算，都成为现代技术发展的重要课题。只有一些典型几何形状和结构相对简单的问题才有可能求得严格的解析解。当代电磁场工程中高频电磁场问题的主要特点是电磁系统的高度复杂性。虽然对很多典型问题的解析分析仍然能帮助人们加深对电磁规律的认识，但作为工程问题的解决，解析方法往往无能为力。能够较广泛发挥作用的，唯有各种数值方法。数值方法数值方法其中主要有：属于频域技术的有限元法（FEM）、矩量法(MM)和单矩法等；属于时域技术方面的时域有限差分法(FDTD)、传输线矩阵法(TLM)和时域积分方程法等。此外，还有属于高频技术的几何衍射理论(GTD)和衍射物理理论(PLD)等。各种方法都具有自己的特点和局限性，在实际中经常把它们相互配合而形成各种混合方法。时域有限差分（FDTD）方法于1966年由K.S.Yee[] 提出并迅速发展，且获得广泛应用。K.S.Yee的运算简化为乘以的代数运算即可，由此把一个四维电磁场问题简化成三维空间的场问题。因此，这种频域分析法对分析波导结构的稳态特性是非常有效的。然而，随着微波与毫米波技术的发展，仅仅研究波导结构的稳态特性是远远不够的，还必须研究其瞬态特性。若在频域中研究波导结构的瞬态特性，由于要对每一个频率点进行单独计算，如果入射波具有较宽的频谱，那么计算工作量会很大，将四维电磁场问题化为三维场问题的频域分析法的优越性就不明显。而且，频域分析一般假设在单模传输的条件下，如果遇到波导不均匀问题，在不均匀区附近往往会激励起高次模，单模传输不能反映瞬态电磁波传输的真实过程。此外，频域分析以线性叠加原理为基础，原则上只对线性波导系统有效，对于非线性结构，频域分析法就显得十分困难。而在时域中分析波导结构的瞬态特性，尽管遇到的是四维电磁场问题，用计算机进行数值模拟需要较大的内存与较多的计算时间，但这对于计算机技术高速发展的今天