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

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