时域有限差分法FDTD算法的基本原理及仿真.doc

.FDTD算法)时域有限差分法(后上的一篇论文建立起来的,。这种方法通过将在相互交织的网格空,分式而直接在时域求解,通过建立时间离散的递进序列间中交替计算电场和磁场。旋度方程转化为差分形MaxwellFDTD算法的基本思想是把带时间变量的式,模拟出电子脉冲和理想导体作用的时域响应。需要考虑的三点是差分格式、采用一定的网格划分有限差分通常采用的步骤是:解的稳定性、吸收边界条件。建对场内的偏微分方程及各种边界条件进行差分离散化处理,方式离散化场域;立差分格式,得到差分方程组;结合选定的代数方程组的解法,编制程序,求边值问题的数值解。,FDTD方法由Maxwell这样达到了在一定体积内和一段时间上直接将微分运算转换为差分运算,近似,对连续电磁场数据的抽样压缩。方程的旋度方程组为:MaxwellHE??HE??EH???????????(1)mtt??在直角坐标系中,(1)式可化为如下六个标量方程:H??E?H?E??H?E?yx??zE???y?xz??H????x?tz??y?xmtz??y????E?H?H???H?E?E??y??xzE???y?xz??H????(,2)?yt?x??zymt?x??z??H??H?E?E??E?H?yy??xzE???xz??H??????zzmt??x?yt??x?y????算法的基础。上面的六个偏微分方程是FDTD可以写成F(x,y,z)Yee首先在空间上建立矩形差分网格,在时刻时刻,t?nn),ki,Fy,k?z,n?t)?j(?zF(x,y,,t)?F(ix,j?3)(用中心差分取二阶精度:对空间离散:nn),j,kk)?F(i?12F,?F(xy,z,t),(i?12j,????2x?O??x?x?ix??xnn),k2i,j?Fj(i,?12,k)?1(),(?Fx,yz,tF????2y??O?yy??jyy??nn(i,j,k?2k?1)?F12),,F(),,,(?Fxyztij????2zO???zkz??z?z?;..对时间离散:212n?n?1)k,j,(i(i,j,k)F?F(x,y,z,t)?F????2tO???4)(tn?t?t??t:的六个分量,如下图放置和HYee把空间任一网格上的EEHEEEEyEzEzHyHxEyExyox氏网格及其电磁场分量分布图1Yee中,空间上连续分布的电磁场物理量离散的空间排布如图所示。FDTD在由图可见,电场和磁场分量在空间交叉放置,各分量的空间相对位置也适合于方程的差分计算,能够恰当地描述电磁场的传播特性。同时,电场和磁Maxwell旋度方程离散Maxwell抽样时间间隔相差半个时间步,使场在时间上交替抽样,而不需要进行矩阵求逆运从而可以在时间上迭代求解,以后构成显式差分方程,就可以逐步推进地求得以后FDTD算。因此,由给定相应电磁问题的初始条件,:各个时刻空间电磁场的分布。根据这一原则可以写出六个差分方程?t)?,2j,k(i?1/1??),j,2k(i?1/21nn?)j,ki?1?/2,,(i?1,jk)?.E(E?xxt)?/2,j,k(i?1?1?),k?1/2,