计算电磁学之FDTD算法的MATLAB语言实现.docx

SouthChinaNormalUniversity

课程设计实验报告

课程名称： 计算电磁学
指导老师：
专业班级： 2014级电路与系统
姓名：
学号：

FDTD算法的MATLAB语言实现

摘要：时域有限差分(FDTD)算法是K.S.Yee于1966年提出的直接对麦克斯韦方程作差分处理,用来解决电磁脉冲在电磁介质中传播和反射问题的算法。其基本思想是:FDTD计算域空间节点采用Yee元胞的方法,同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区,这两个区域是以连接边界相连接,最外边是采用特殊的吸收边界,同时在这两个边界之间有个输出边界,用于近、远场转换;在连接边界上采用连接边界条件加入入射波,从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件,尽量消除反射波在吸收边界上的非物理性反射波。
本文主要结合FDTD算法边界条件特点,在特定的参数设置下，用MATLAB语言进行编程，在二维自由空间TEz网格中，实现脉冲平面波。
关键词：FDTD；MATLAB；算法

1绪论
1.1课程设计背景与意义
20世纪60年代以来，随着计算机技术的发展，一些电磁场的数值计算方法逐步发展起来，并得到广泛应用，其中主要有：属于频域技术的有限元法（FEM）、矩量法(MM)和单矩法等；属于时域技术方面的时域有限差分法(FDTD)、传输线矩阵法(TLM)和时域积分方程法等。其中FDTD是一种已经获得广泛应用并且有很大发展前景的时域数值计算方法。时域有限差分（FDTD）方法于1966年由K.S.Yee提出并迅速发展，且获得广泛应用。K.S.Yee用后来被称作Yee氏网格的空间离散方式，把含时间变量的Maxwell旋度方程转化为差分方程，并成功地模拟了电磁脉冲与理想导体作用的时域响应。但是由于当时理论的不成熟和计算机软硬件条件的限制，该方法并未得到相应的发展。20世纪80年代中期以后，随着上述两个条件限制的逐步解除，FDTD便凭借其特有的优势得以迅速发展。它能方便、精确地预测实际工程中的大量复杂电磁问题，应用范围几乎涉及所有电磁领域，成为电磁工程界和理论界研究的一个热点。目前，FDTD日趋成熟，并成为分析大部分实际电磁问题的首选方法。

1.2时域有限差分法的发展与应用

经过四十多年的发展，FDTD已发展成为一种成熟的数值计算方法。在发展过程中，几乎都是围绕几个重要问题展开的，即数值稳定性、计算精度、数值色散、激励源技术以及开域电磁问题的吸收边界条件等。数值稳定和计算精度对任何一种数值计算方法都是至关重要的。其中激励源的设计和引入也是FDTD的一个重要任务。目前，应用最广泛的激励源引入技术是总场/散射场体系。对于散射问题，通常在FDTD计算空间中引入连接边界，它将整个计算空间划分为内部的总场区和外部的散射场区，如图1.1。利用Huygens原理，可以在连接边界处引入入射场，使入射场的加入变得简单易行。

图1.1总场/散射场区

2FDTD法基本原理
2.1Maxwell方程和Yee氏算法
根据电磁场基本方程组的微分形式，若在无源空间，其空间中的煤质是各向同性、线性和均匀的，即煤质的参数不随时间变化且各向同性，则Maxwell旋度方程可写成：

式中，E是电场强度，单位为伏/米（V/m）；H是磁场强度，单位为安/米（A/m）；ε表示介质介电系数，单位为法拉/米（F/m）；μ表示磁导系数，单位为亨利/米（H/m）；σ表示介质电导率，单位为西门子/米（S/m）；σm表示导磁率，单位为欧姆/米（Ω/m）。
K.S.Yee在1966年建立了如图2.1所示的空间网格，这就是著名的Yee氏元胞网格。

图2.1Yee氏网格及其电磁场分量分布

电场和磁场被交叉放置,电场分量位于网格单元每条棱的中心,磁场分量位于网格单元每个面的中心,每个磁场(电场)分量都有4个电场(磁场)分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足,而且还允许旋度方程在空间上进行中心差分运算,同时也满足了法拉第电磁感应定律和安培环路积分定律,也可以很恰当地模拟电磁波的实际传播过程。

2.2时域有限差分法相关技术
2.2.1数值稳定性问题
FDTD方程是一种显式差分方程，在执行时，存在一个重要的问题：即算法的稳定性问题。这种不稳定性表现为在解显式方程时，随着时间步数的继续增加，计算结果也将无限制地增加。Taflove于1975年对Yee氏差分格式的稳定性进行了讨论，并导出了对时间步长的限制条件。数值解是否稳定主要取决于时间步长Δt与空间步长Δx、Δy、Δz的关系。对于非均匀媒质构成的计算空间选用如下的稳定性条件：

上式是空间和时间离散之间应当满足的关系，