fdtd算法的matlab程序,基于MATLAB的FDTD算法编程

文章编号　100426410 (2006) 040043204 基于MATLAB 的 FDTD 算法编程 赵　嘉 (广西工学院 计算机工程系, 广西 柳州　545006) 摘　要: 介绍了时域有限差分(FDTD)法的基本原理, 推导了二维 TM 模 Yee 算法的 FDFD 表达式, 并结合算例阐述了基于MA TLAB 编程的基本方法。 关　键　词:MA TLAB; FDTD; Yee 算法; 编程中图分类号: TP321　　　　文献标识码:A 收稿日期: 20060820 作者简介: 赵　嘉(19652) , 女, 广西柳州市人, 广西工学院计算机工程系讲师。 图 1　二维 TM 波YEE cell 0　引言时域有限差分(FDTD) 法是六十年代由 KS1Yee 提出并首先用于求解电磁散射问题, 其主要思路是在 三维空间和时间轴上对场量离散, 并且用中心差分代替偏微分, 将麦克斯韦方程组转化为差分方程, 通过在时间轴和空间轴上采取蛙跳法(leapfrog)逐步推进地求解, 最终求出一定边值与初值条件下的空间场解。随着计算机技术的发展, 近年来FDTD 计算技术也得到了越来越多的应用。对于FDTD 算法的编程求解, 最常用的有VC 和 FORTRUN , 而MA TLAB 作为一种可视化效果很好的科学计算软件, 在FDTD 计算中能充分发挥编程简单、可视化程度高、能显示动态场效果的特点。本文首先推导二维 TM 模 Yee 算法的 FDFD 表达式, 并讨论如何结合FDTD 算法边界条件的特点, 用MA TLAB 语言进行编程的方法和应注意的问题。 1　二维 TM 模Yee 算法 111　算式推导 在自由空间中, 对于二维问题,z= 0, 对于 TM 波, H z= Ex = Ey = 0,MAXW ELL 的两个旋度方程可分解为[1]: Ez y = -Λ0H x t , 　　Ez x =Λ0H y 5t , 　　Ε0Ez t =H y x -H x 5y　　构造二维 TM 波 YEE cell 如图 1 所示: 按 YEE 元胞对上式偏导用中心差商代替, 可得: E z (i, j+ 1 2 , n) - E z (i, j- 1 2 , n) ∃y = -Λ0 H x (i, j+ 1 2 , n+ 1 2 ) - H x (i, j+ 1 2 , n- 1 2 ) ∃t (1) 第 17 卷　第 4 期 　　　　　　　　　　广 西 工 学 院 学 报 　　　　　　　　　　　Vol17　No4 2006 年 12 月 　JOURNAL OF GUAN GX IUN IV ERSITY OF TECHNOLOGY　　　 Dec12006 E z (i+ 1 2 , j, n) - E z (i- 1 2 , j, n) x =Λ0 H y (i+ 1 2 , j, n+ 1 2 ) - H y (i+ 1 2 , j, n- 1 2 ) ∃t (2) 图 2　算例计算区域 　　　Ε0 E z (i, j, n+ 1 2 ) - E z (i, j, n- 1 2 ) t　　　　 = H y (i+ 1, j, n+ 1 2 ) - H y (i, j, n+ 1 2 ) ∃x - H x (i, j+ 1, n+ 1 2 - H y (i, j, n+ 1 2 ) ∃y (3)　　方程中出现了半个网格和半个时间步 , 为便于编程, 可将上面差分式改写为如下FDTD 算式: 　　　　　　　　　H x (i, j, k+ 1) = H x (i,