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

.页眉. .页脚. S outh China Normal University 课程设计实验报告课程名称: 计算电磁学指导老师: 专业班级: 2014 级电路与系统姓名: 学号: FDTD 算法的 MATLAB 语言实现摘要: 时域有限差分(FDTD) 算法是 于1966 . .页脚. 程作差分处理, 用来解决电磁脉冲在电磁介质中传播和反射问题的算法。其基本思想是:FDTD 计算域空间节点采用 Yee 元胞的方法, 同时电场和磁场节点空间与时间上都采用交错抽样;把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区, 这两个区域是以连接边界相连接, 最外边是采用特殊的吸收边界, 同时在这两个边界之间有个输出边界,用于近、远场转换;在连接边界上采用连接边界条件加入入射波, 从而使得入射波限制在总场区域;在吸收边界上采用吸收边界条件,尽量消除反射波在吸收边界上的非物理性反射波。本文主要结合 FDTD 算法边界条件特点, 在特定的参数设置下,用 MATLAB 语言进行编程,在二维自由空间 TEz 网格中,实现脉冲平面波。关键词: FDTD ;MATLAB ;算法 1 绪论 课程设计背景与意义 20世纪 60年代以来,随着计算机技术的发展,一些电磁场的数值计算方法逐步发展起来,并得到广泛应用,其中主要有:属于频域技术的有限元法( FEM )、矩量法(MM) 和单矩法等; 属于时域技术方面的时域有限差分法(FDTD) 、传输线矩阵法(TLM) 和时域积分方程法等。其中 FDTD 是一种已经获得广泛应用并且有很大发展前景的时域数值计算方法。时域有限差分( FDTD )方法于 1966 年由 提出并迅速发展, 且获得广泛应用。 用后来被称作 Yee 氏网格的空间离散方式,把含时间变量的 Maxwell 旋度方程转化为差分方程,并成功地模拟了电磁脉冲与理想导体作用的时域响应。但是由于当时理论的不成熟和计算机软硬件条件的限制,该方法并未得到相应的发展。 20世纪 80年代中期以后,随着上述两个条件限制的逐步解除, FDTD 便凭借其特有的优势得以迅速发展。它能方便、精确地预测实际工程中的大量复杂电磁问题,应用范围几乎涉及所有电磁领域,成为电磁工程界和理论界研究的一个热点。目前, FDTD 日趋成熟,并成为分析大部分实际电磁问题的首选方法。 时域有限差分法的发展与应用经过四十多年的发展, FDTD 已发展成为一种成熟的数值计算方法。. .页脚. 过程中,几乎都是围绕几个重要问题展开的,即数值稳定性、计算精度、数值色散、激励源技术以及开域电磁问题的吸收边界条件等。数值稳定和计算精度对任何一种数值计算方法都是至关重要的。其中激励源的设计和引入也是 FDTD 的一个重要任务。目前,应用最广泛的激励源引入技术是总场/ 散射场体系。对于散射问题,通常在 FDTD 计算空间中引入连接边界,它将整个计算空间划分为内部的总场区和外部的散射场区,如图 。利用 Huygens 原理,可以在连接边界处引入入射场,使入射场的加入变得简单易行。图 总场/散射场区 2 FDTD 法基本原理 Maxwell 方程和 Yee 氏算法根据电磁场基本方程组的微分形式,若在无源空间,其空间中的煤质是各向同性、线性和均匀的, 即煤质的参数不随时间变化且各向同性,则 Maxwell 旋度方程可写成: 式中,E是电场强度,单位为伏/米(V/m );H是磁场强度,单位为安/米(A/m ); ε表示介质介电系数,单位为法拉/ 米( F/m );μ表示磁导系数,单位为亨利/ 米(H/m );σ表示介质电导率,单位为西门子/米(S/m );σ m表示导磁率,单位为欧姆/米( Ω/m)。.页眉. .页脚. 在1966 年建立了如图 所示的空间网格,这就是著名的 Yee 氏元胞网格。 Yee 氏网格及其电磁场分量分布电场和磁场被交叉放置, 电场分量位于网格单元每条棱的中心, 磁场分量位于网格单元每个面的中心, 每个磁场(电场) 分量都有 4 个电场(磁场) 分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足, 而且还允许旋度方程在空间上进行中心差分运算, 同时也满足了法拉第电磁感应定律和安培环路积分定律,也可以很恰当地模拟电磁波的实际传播过程。 时域有限差分法相关技术 数值稳定性问题 FDTD 方程是一种显式差分方程,在执行时,存在一个重要的问题:即算法的稳定性问题。这种不稳定性表现为在解显式方程时,随着时间步数的继续增加, 计算结果也将无限制地增加。 Taflove 于1975 年对 Yee