时域有限差分法对平面TE波的MATLAB仿真设计

《时域有限差分法对平面TE波的MATLAB仿真设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、word格式时域有限差分法对平面TE波的MATLAB仿真姓名：王云璐学号：2011302021（西北工业大学电子信息学院08041103，陕西西安，710072）摘要：本文分析FDTD算法的基本原理及两种典型边界条件的算法特点，通过MATLAB程序对TE波进行了仿真，模拟了高斯制下完全匹配层中磁场分量瞬态分布。得到了相应的磁场幅值效果图。最后得出用Matlab语言对FDTD算法编程的几点结论。关键词：FDTD；MATLAB；PML；平面TE波1引言电磁场的有限差分一般是在时域上进行的，随着计算机技术的发展和应用，近年来，时域计算方法越来越受到重

2、视。时域有限差分法具有简单、结果直观、网格剖分简单等优点。近些年FDTD发展的十分迅速，在许多方面都有重要应用，包括天线设计，微波电路设计，电磁兼容分析，电磁散射计算，光子学应用等等。时域有限差分(FDTD)算法是K.S.Yee于1966年提出的直接对麦克斯韦方程作差分处理，用来解决电磁脉冲在电磁介质中传播和反射问题的算法。基本思想是：FDTD计算域空间节点采用Yee元胞的方法，同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区，这两个区域是以连接边界相连接，最外边是采用特殊的吸收边界，同

3、时在这两个边界之间有个输出边界，用于近、远场转换；在连接边界上采用连接边界条件加入入射波，从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件，尽量消除反射波在吸收边界上的非物理性反射波。本文主要简述了FDTD算法的基本原理，解的稳定性以及边界条件特点，并用Matlab语言进行对平面TE波进行了编程计算。2FDTD基本原理时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。FDTD空间网格单元上电场和磁场各分量的分布如图1所示。....word格式图1

4、.Yee氏网格及其电磁场分量分布电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元每个面的中心，每个磁场(电场)分量都有4个电场(磁场)分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足，而且还允许旋度方程在空间上进行中心差分运算，同时也满足了法拉第电磁感应定律和安培环路积分定律，也可以很恰当地模拟电磁波的实际传播过程。3解的稳定性FDTD方程是一种显式差分方程，在执行时，存在一个重要的问题：即算法的稳定性问题。这种不稳定性表现为在解显式方程时，随着时间步数的继续增加，计算结果也将无限制地增加。Ta

5、flove等于1975年对Yee氏差分格式的稳定性进行了讨论，并导出了对时间步长的限制条件。数值解是否稳定主要取决于时间步长与空间步长、、的关系。对于非均匀媒质构成的计算空间选用如下的稳定性条件：(1)(1)式是空间和时间离散之间应当满足的关系，又称为Courant稳定性条件。若采用均匀立方体网格：(2)其中，为计算空间中的电磁波的最大速度。4FDTD算法的几个边界条件....word格式边界条件是对FDTD算法应用的另外一个关键问题。时域有限差分法是在计算机的数据存储空间中对连续的实际电磁波的传播过程在时间按空进程上进行数字模拟。在计算电磁场

6、的辐射、散射等问题中，边界总是开放的，电磁场从实际的角度看总是占据无限大空间，由于计算机内存总是有限的，只能模拟有限的空间。也就是说，时域有限差分网格将在某处被截断。如何处理截断边界，使之需要考虑的无限空间的差异尽量的小，这是FDTD算法要解决的重要问题。在当今FDTD算法中，涉及到的边界条件主要有Mur吸收边界条件、廖氏吸收边界条件、超吸收边界条件、完全匹配层(PML)，在本文中主要介绍Mur吸收边界条件和完全匹配层(PML)。4.1Mur吸收边界条件其具体推导过程可参考文献[4]，递推公式如下一阶条件下的递推公式(3)二阶条件下Mur的FD

7、TD递推公式(4)这就是Mur所建议的具有二阶近似的、适用于二维问题的近似吸收边界条件。它在FDTD中有着广泛的应用。4.2PML吸收边界条件1994年，Berengert提出了用完全匹配(PML)来吸收外向电磁波，它将电磁场分量在吸收边界分裂，并能分别对各个分裂的场分量赋予不同的损耗。完全匹配层是一种非物理性的电磁波吸收层，有关性质和在具体条件下PML格式的FDTD递推公式见文献[5]，本文采用PML吸收边界条件。5MATLAB的仿真程序及结果5.1MATLAB程序及相应说明%二维FDTDTE波仿真Clearall；%定义常数pi=3.141

8、5；c=3.0e10；%高斯制下光速f=1.0e15；%频率lambda=c/f；%波长....word格式nmax=400；%时间步数del_s=l