本文对传统时域有限差分法做出了改进,在空间和时间上采用高阶精度的离散格式,使空间差分和时间差分都达到了4阶,得到了新的广义高阶算法。详细推导了高阶FDTD方法的迭代公式,并推导了高阶算法对应的PML吸收边界条件,以及总场/散射场的连接边界。最后将高阶FDTD方法应用到二维金属方柱的电磁散射特性分析中,并将数值计算结果与传统的FDTD方法作了比较和分析,数值结果表明高阶FDTD方法具有更高的计算精度。 学报图1无线长的线源在xy面内辐射的zE分量的等值线图3高阶算法中总场/散射场连接边界设置在广义高阶算法中，需要根据场点是位于总场或散射场来具体修改公式。以一维情况为例。将一维FDTD（2M,4）区域划分成总场和散射场，并令总场边界00zkz，且位于总场边界上的电场节点0()xEk属于总常总场边界散射场区0kM总场区电场节点磁场节点图2总场边界处场点分布在总场区内，从0()xEk一直到0(1)xEkM的电场节点在更新时，需要用到散射场中的磁场信息，而从0(0.5)yHk到0(1.5)yHkM的磁场节点在更新时需要用到散射场中的电场点信息；在散射场中，从0(1)xEkM一直到0(1)xEk的电场节点在更新时需要用到总场中的磁场点信息，而从0(0.5)yHkM到0(0.5)yHk的磁场节点在更新时需要用到总场中的电场点信息。下面给出具体的修正公式，在总场区中电场点的修7月微波学报图3（b）俯视图图3三维时入射波的引入4散射数值算例为了验证本文方法的有效性，对平面波垂直入射到金属圆柱的雷达散射界面（RCS）进行计算。圆柱半径r，空间入射波的波长0.01m。一般算法中，取网格步长为/40，计算区域大小为160160个网格。而高阶算法中，网格步长为一般算法的两倍本文由公司网站滚圆机网站采集转载中国知网整理!   http://www.gunyuanji158.com/，生的RCS仿真分析-数控滚圆机滚弧机张家港电动液压滚圆机滚弧机折弯机即/20，计算区域为8080个网格。图4中给出了传统FDTD方法和高阶FDTD（2，4）方法的计算结果。从图中可以看出两种算法得到的双站RCS曲线图吻合较好，而高阶FDTD方法的计算精度比传统方法要高一些。图4二维金属圆柱的双站RCS5结束语为了提高时域有限差分法的计算效率和精度，本文将FDTD公式中的空间差分和时间差分都达到了4阶，得到了新的广义高阶算法，并将高阶FDTD方法应用到二维金属方柱的电磁散射特性分析中。文中详细推导了高阶FDTD（2,4）算法的迭代公式，并推导了改进算法的PML吸收边界设置和总场边界条件。二维金属圆柱算例的数值结果表明，高阶FDTD理论在空间步长为FDTD方法的2倍时，所获得的计算精度仍较传统FDTD方法的效果好。而且由于在高阶算法中空间网格节点减少，故该方法相比一般算法节省了计算机内存据雷达散射截面(RCS)的计算原理,利用高阶矩量法(HO-MoM)的后处理数据,实现了整体目标中,局部电磁流产生的RCS计算与分析。本文借助HO-MoM能够精确计算电大尺寸目标的优势,对一整架飞机进行了RCS计算。后处理时提取其中若干区域的电磁流分别求电场,进而算出产生的局部RCS。对于局部电磁流产生RCS的分析,将为后续RCS有针对性的减缩提供参考价值。本文考虑到不同区域之间的耦合,并不将局部从整体中剥离,保证了电流连续性。该研究对于隐身目标的散射特性分析具有一定意义。 生的RCS仿真分析-数控滚圆机滚弧机张家港电动液压滚圆机滚弧机折弯机本文由公司网站滚圆机网站采集转载中国知网整理!   http://www.gunyuanji158.com/