二维时域电场积分方程新平均稳定法.docVIP

二维时域电场积分方程新平均稳定法 　　摘 要:对二维理想金属导体柱的时域电场积分方程(TDEFIE)后期振荡问题提出了一种新的时间平均方法,此方法用前一时刻和后一时刻的值来修正当前时刻的电流值,消除了振荡,而且相对以往的平均方法提高了MOT算法的稳定性,对精度的影响很小。这里采用电场积分方程(EFIE)MOT算法的隐式格式(Implicit),在TM,TE高斯平面波激励情况下分别计算了两个例子:圆柱和方柱。又采用一种新的时间基函数BLIFs来结合新的时间平均法,取得了很好的效果。经过Matlab软件的数值分析,可以看出此时间平均法更加稳定和精确。 　　关键词:时域电场积分方程;时间平均法;隐式算法;时间步进算法 　　中图分类号:TN820.1文献标识码:A 　　文章编号:1004-373X(2010)04-174-03 　　 　　A Novel Average Technique for the Stability of Two-demension TDEFIE 　　ZHANG Cheng,LUO Wei,WU Xianliang 　　(School of Electronic Science and Technology,Anhui University,Hefei,230039,China) 　　Abstract:A new time average method using time domain electric integral equation,which adds the front and back of the current electricity to eliminate oscillation,it is presented to calculate two-demension conducting cylinder for improving the stability more than other average approaches and the precision influenced by the method can be almost ignored.An implicit method using EFIE is adopted to calculate two examples apiece inspired by TM and TE gauss plane wave:circle cylinder and square cylinder.Except that,a new temporal basis funcion BLIFs used in the novel average technique acquire better impression.??The results of Matlab simulation prove the stability and precision of the approach. 　　Keywords:time domain electric integral equation;time average method;implicit method;marching-on-in-time 　　 　　0 引 言 　　 　　TDIE是四种数值算法中发展较为滞后的一种,但由于其在电磁散射和辐射中无可比拟的优势,得到越来越多的关注。影响它发展的原因主要是后期不稳定现象。许多学者经过潜心研究对其做出了各种假设,其中具有代表性的有:A.G.Tihuis认为在建模时对时间和空间进行了一系列近似带来误差,从而在时间步进过程中积累,最终导致发散;B.P.Rynne认为目标本身存在某些内谐振频率。以上只是许多猜想中的两个假设,人们基于以上假设对TDIE稳定性的改进做了很多尝试,在一定程度上改善了MOT算法的时间后期稳定性。这里基于时域电场积分方程的隐式算法[1-3],提出了一种??新的时间平均方法,并与时间基函数(BLIFs)相结合[4],进一步的推迟了振荡产生时刻;其后将本文方法与B.P.Rynne的三步平均法和吴振森的先五步后三步的平均方法[5]做了比较,证明此方法的优越性。 　　 　　1 TDEFIE的隐式MOT迭代过程 　　 　　由理想导体表面切向电场为零的边界条件,得到时域电场积分方程(TDEFIE): 　　[??A??t+??φ]tan??=[Ei]tan??(1) 　　其中矢量位为: 　　A=μ4π∫c∫∞z′=-∞??J(ρ ′,t-R/c)Rdz′dc′(2) 　　标量位为: 　　φ=14πε∫c∫∞z′=-∞??qs(ρ ′,t-R/c)Rdz′dc′(3) 　　由电流连续性定理??s#822