时域有限差分法(FDTD基本原理及仿真
算法)的
时域有限差分法（FDTD算法）
时域有限差分法是1966年K.S.Yee发表在AP上的一篇论文建立起来的，后被称为Yee网格空间离散方式。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解, 通过建立时间离散的递进序列, 在相互交织的网格空间中交替计算电场和磁场。
FDTD算法的基本思想是把带时间变量的Maxwell旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应。需要考虑的三点是差分格式、解的稳定性、吸收边界条件。有限差分通常采用的步骤是：采用一定的网格划分方式离散化场域；对场内的偏微分方程及各种边界条件进行差分离散化处理，建立差分格式，得到差分方程组；结合选定的代数方程组的解法，编制程序，求边值问题的数值解。 1.FDTD的基本原理
FDTD方法由Maxwell旋度方程的微分形式出发，利用二阶精度的中心差分近似，直接将微分运算转换为差分运算，这样达到了在一定体积内和一段时间上对连续电磁场数据的抽样压缩。
Maxwell方程的旋度方程组为：
??H???E?H??E ??E?????mH （1） ?t?t在直角坐标系中，（1）式可化为如下六个标量方程：
??Ex?Hz?Hy?????Ex??Ez??Ey????Hx??H?mx??y?z?t?y?z?t????Ey?Hx?Hz???Hy??????Ey?，?Ex??Ez?????mHy? （2）
?z?x?t?z?x?t???Hy?Hx??Ey?Ex?Ez??Hz?????Ez???????mHz??x?y?t?x?y?t????上面的六个偏微分方程是FDTD算法的基础。
Yee首先在空间上建立矩形差分网格，在时刻n?t时刻，F(x,y,z)可以写成
F(x,y,z,t)?F(i?x,j?y,k?z,n?t)?Fn(i,j,k) （3）
用中心差分取二阶精度： 对空间离散：
?F(x,y,z,t)?xFn(i?12,j,k)?Fn(i?12,j,k)2 ????O?xx?i?x?x???F(x,y,z,t)?yFn(i,j?12,k)?Fn(i,j?12,k)2 ????O?yy?j?y?y???F(x,y,z,t)?zFn(i,j,k?12)?Fn(i,j,k?12)2????O?z z?k?z?z??
对时间离散：
?F(x,y,z,t)
?tFn?12(i,j,k)?Fn?12(i,j,k)2????O?t （4） t?n?t?t??Yee把空间任一网格上的E和H的六个分量，如下图放置:
zExHzExEyEzHxEyoExyEzHyEzEyx
图1 Yee氏网格及其电磁场分量分布
在FDTD中，空间上连续分布的电磁场物理量离散的空间排布如图所示。由图可见，电场和磁场分量在空间交叉放置，各分量的空间相对位置也适合于Maxwell方程的差分计算，能够恰当地描述电磁场的传播特性。同时，电场和磁场在时间上交替抽样，抽样时间间隔相差半个时间步，使Maxwell旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求逆运算。因此，由给定相应电磁问题的初始条件，FDTD就可以逐步推进地求得以后各个时刻空间电磁场的分布。根据这一原则可以写出六个差分方程:
?(i?1/2,j,k)?t2?(i?1/2,j,k)nExn?1(i?1,j,k)?.E(i?1?/2,j,k)?(i?1/2,j,k)?tx1?2?(i?1/2,j,k)?t1?..?(i?1/2,j,k)1??(i?1/2,j,k)?t2?(i?1/2,j,k)1?Hzn?1/2(i?1/2,j,k)?Hzn?1/2(i?1/2,j?1/2,k)Hy[??y
n?1/2（5）
(i?1/2,j,k?1/2)?Hyn?1/2(i?1/2,j,k?1/2)?z]其余的也如法可以写出，每个网格点上的个场分两的新值依赖于该点在前一时间步长时刻的值机该点周围的临近点上另一场量在早半个时间步长时的值。因此任一时刻可一次算出一个点，并行算法可计算出多个点。通过这些运算可以交替算出电场磁场在各个时间步的值。
时域有限差分法(FDTD算法)的基本原理及仿真
