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

格将在某处被截断。这要求在网格截断处不能引起波的明显反射,因而对向外传播的波而言,就像在无限大的空间传播一样,一种行之有效的方法是在截断处设置一种吸收边界条件。使传播到截断出的波被边界吸收而不产生反射。下面只给出 Engquist-Majda 吸收边界条件,采用 Mur 差分格式,其总体虚假反射在 1%~5% 之间。一维一阶近似情形, x=0 边界: )]0()1([)1()0( 1 1uuxtc xtcuu nnn???????????( 14) 二维二阶近似情形, x=0 边界: )]1,1(),1(2)1,1()1,0(),0(2)1,0([ .)()(2 )( )],1(),0([ 2 )],0(),1([),1(),0( 2 2 1111????????????????????????????????????jWjWjWjWjWjW xtcy xtcjWjWxtc x jWjWxtc xtcjWjW nnnnnn nn nnnn三维二阶近似情形, x=0 边界: ( 15) )]1,,1(),,1(2)1,,1()1,,0(),,0(2)1,,0([)()(2 )( )],1,1(),,1(2),1,1(),1,0(),,0(2),1,0([ .)()(2 )( )],,1(),,0([ 2 )],,0(),,1([),,1(),,0( 2 2 2 2 1111????????????????????????????????????????????????????kjWkjWkjWkjWkjWkjWxtcz xtc kjWkjWkjWkjWkjWkjW xtcy xtckjWkjWxtc x kjWkjWxtc xtckjWkjW nnnnnn nnnnnn nn nnnn 5.仿真 m文件见附件,程序表现的是使用二维 FDTD 算法对 TE 波的仿真。( 16)