[发明专利]分析导体瞬态电磁散射特性的时域高阶Nystrom方法

说明书

技术领域

本发明属于分析导体目标瞬态电磁散射特性的时域积分方程方法，具体是一种分析导体瞬态电磁散射特性的时域高阶Nystrom方法。

背景技术

雷达目标电磁散射特性的获取与分析是电磁问题中的一个非常重要研究领域，目标的电磁散射波是雷达探测、遥感观测以及地质勘测邓众多应用的信息来源，散射特性的定量分析是这些应用系统在设计和工作时的主要依据。雷达目标的形状和体积等物理量都是通过对雷达散射截面等参数进行计算得出的。因此，对于各种目标散射特性的研究在这些应用领域具有特别重要的意义。

近年，瞬态电磁散射特性的分析越来越引起科研学者和工程人员的关注。相比于其它方法，时域积分方程方法非常适合于理想电导体的瞬态电磁散射特性的分析。出现最早、研究最多、且最为成熟的就是基于时间步进的时域积分方程方法(S.M.RaoandD.R.Wilton,“Transientscatteringbyconductingsurfacesofarbitraryshape,”IEEETrans.AntennasPropag.,vol.39,no.1,pp.56–61,1991.)。随着宽频带和非线性电磁散射和辐射系统的仿真与设计需求的增加，一种对于网格鲁棒性的时域求解技术显得非常重要。

但是，对于导体瞬态电磁散射特性的分析，因为传统的基于RWG基函数的时域积分方程，要求离散的三角形单元共内边，这极大限制了该时域方法在某些实际问题中的应用。而时域高阶Nystrom方法所用的矢量插值基函数定义在离散曲面三角形单元内的插值点处，没有共内边的要求，对离散网格有鲁棒性的优点。

发明内容

本发明的目的在于提供一种分析导体瞬态电磁散射特性的时域高阶Nystrom方法，步骤如下：

第一步，建立导体表面时域积分方程，即根据理想导体表面切向连续的边界条件，在金属表面可以建立时域电场积分方程和时域磁场积分方程，入射电场和磁场分别为已知激励，通常使用调制高斯平面波作为入射场，散射场可以用待求的表面瞬态未知电流来表示；

第二步，对导体表面时域积分方程采用三角基函数进行时间上的离散，并采用二阶曲面三角形单元进行空间上的离散；

第三步，在时间上采用加辽金测试，空间上点匹配，形成待求解的矩阵方程，未知电流为导体瞬态面电流；

第四步，矩阵方程的求解以及瞬态电磁散射参数的计算。

本发明与传统的基于RWG的时域积分方程方法相比，其显著优点是：对离散网格单元具有鲁棒性，即不需要离散三角形单元共内边。

附图说明

图1是曲三角形单元映射到局部空间(u,v)示意图。

图2是曲面三角形网格不共内边的导体球示意图。

图3是导体球在不同频率点处的双站雷达散射截面(RCS)，(a)：频率为90MHz(b)：频率为120MHz(c)：频率为150MHz(d)：频率为180MHz(e)：频率为210MHz。

具体实施方式

针对导体结构，本发明采用时域积分方程方法可以分析其瞬态电磁散射特性。时域高阶Nystrom方法与传统的基于RWG基函数的时域积分方程方法相比，具有对离散网格鲁棒性的优点。因为传统的基于RWG基函数的时域积分方程，要求离散的三角形单元共内边，这极大限制了该方法在某些实际问题中的应用。而时域高阶Nystrom方法所用的矢量插值基函数定义在离散曲面三角形单元内的插值点处，没有共内边的要求。

下面结合附图对本发明作进一步详细描述。

结合图1，本发明基于分析导体瞬态电磁散射特性的时域高阶Nystrom方法，步骤如下：

第一步，令电磁波照射到导体结构上，在导体表面上产生表面感应面电流J，根据理想导体的电场边界条件，即金属表面的总场切向分量为0，得到导体目标的时域积分方程TDIE，如下