1.一种曲面共形阵列目标电磁特性分析方法,其特征在于,包括:步骤1,根据阵列目标的结构特性将目标整体分解为若干分解部件,对每个分解部件进行表面三角剖分并定义RWG函数,在每个分解部件外部创建外部等效源;步骤2,基于计算平台的硬件配置,创建若干个并行线程,指定其中一个为主线程,剩余的为子线程,根据分解部件的数量和子线程的数量,分配每个子线程负责的分解部件数量和编号;步骤3,对于每个子线程,利用分解部件自身的RWG函数和对应的外部等效源,独立构建该子线程负责的所有分解部件的综合函数,将结果存储至本地内存并上报主线程;步骤4,重复步骤2,根据分解部件自身的RWG函数和照射电磁波的特性参数计算每个子线程负责的所有分解部件的自阻抗矩阵、激励矩阵以及不同分解部件之间的互阻抗矩阵,对计算得到的自/互阻抗矩阵、激励矩阵分别进行压缩获得压缩阻抗矩阵和压缩激励矩阵,将结果存储至本地内存并上报主线程;步骤5,根据压缩阻抗矩阵和压缩激励矩阵构建针对目标整体结构的压缩矩阵方程,求解该方程获得目标上定义的综合函数的电流系数,根据综合函数的定义获得目标感应电流源,进一步根据自由空间理想电流源的辐射特性计算得到目标整体的电磁特性;所述步骤4包括:步骤41,对于分解部件b i ,按照传统MoM的阻抗矩阵计算方法,计算分解部件b i 基于自身RWG函数的自阻抗矩阵,记为Z ii ;步骤42,利用分解部件b i 上定义的综合函数展开系数矩阵P i 对自阻抗矩阵进行压缩,压缩方式为z ii =[P i ] T Z ii [P i ],其中,上标T为矩阵的转置运算符,z ii 为部件b i 压缩后的自阻抗矩阵;步骤43,对于分解部件b i 与部件b j 之间的互阻抗,按照传统MoM的互阻抗矩阵计算方法,计算两者基于RWG函数的互阻抗矩阵,记为Z ij ;步骤44,分别利用分解部件b i 和分解部件b j 上定义的综合函数展开系数矩阵P i 、P j 对互阻抗矩阵进行压缩,压缩方式为z ij =[P i ] T Z ij [P j ],其中,z ij 为分解部件b i 和部件b j 压缩后的互阻抗矩阵;步骤45,按照传统MoM的激励矩阵计算方法,计算外部照射波与分解部件b i 上定义的RWG函数的内积,记为V i ;步骤46,利用分解部件b i 上定义的综合函数展开系数矩阵P i 对内积进行压缩,压缩方式为v i =[P i ] T V i ,其中,v i 为照射波对于分解部件b i 的压缩激励矩阵。
2.如权利要求1所述的曲面共形阵列目标电磁特性分析方法,其特征在于,所述步骤1包括:步骤11,将目标按照阵列结构和阵元分布特征分解为若干分解部件并分别编号;步骤12,对每个分解部件进行非结构化三角网格剖分;步骤13,在分解部件外部构建一个虚拟三维封闭表面,对构建的虚拟三维封闭表面进行非结构化三角网格剖分;步骤14,根据分解部件自身的三角剖分网格和构建的虚拟三维封闭表面的三角剖分网格分别定义RWG函数,其中,针对虚拟三维封闭曲面定义的RWG函数即是外部等效源。
3.如权利要求1所述的曲面共形阵列目标电磁特性分析方法,其特征在于,所述步骤3包括:步骤31,根据电场积分方程,利用子线程负责的分解部件自身的RWG函数和外部等效源构建反应方程,求解反应方程获得综合函数的解空间;步骤32,对综合函数的解空间进行奇异值分解,获得综合函数的展开系数矩阵;步骤33,将综合函数的展开系数矩阵与分解部件定义的RWG函数进行线性组合,获得分解部件上定义的综合函数,将每个分解部件构建的综合函数存储至本地内存并上报主线程。
4.如权利要求3所述的曲面共形阵列目标电磁特性分析方法,其特征在于,所述步骤31中的反应方程为:[Z][R]=V-V e (1)上式中,Z表示待分析目标的分解部件基于RWG函数的自阻抗矩阵,R表示综合函数解空间矩阵,V表示外部照射波对于待分析部件产生的激励矩阵,V e 表示待分析目标的分解部件与外部等效源之间的互阻抗矩阵;求解公式(1)得到综合函数的解空间:R=Z -1 (V-V e );步骤32中,对解空间R进行奇异值分解可得:R=UρV H其中,U和V为正交矩阵,即矩阵的每一列之间相互独立正交;ρ表示R的奇异值矩阵;步骤33中,综合函数的展开系数矩阵即为正交矩阵U的前K列,即P=[U 1 ,U 2 ,…,U K ]得到综合函数的展开系数矩阵后,综合函数的构成即为展开系数矩阵P与部件自身定义的RWG函数的线性组合:上式中,F k 表示第k个综合函数k=1,2,…,K,f n 表示定义在部件表面的第n个RWG函数,n=1,2,…,M,M表示部件自身定义的RWG函数总数量,P n,k 表示综合函数展开系数矩阵P的第n行第k列元素。
5.如权利要求1所述的曲面共形阵列目标电磁特性分析方法,其特征在于,所述步骤5包括:步骤51,根据压缩阻抗矩阵和压缩激励矩阵构建目标整体的压缩矩阵方程;步骤52,求解目标整体的压缩矩阵方程获得目标上定义的综合函数电流系数;步骤53,根据综合函数的定义获得待求解目标表面的感应电流源;步骤54,根据理想电流源在自由空间的辐射特性,利用目标表面的感应电流源获得目标在空间的电磁特性。
6.如权利要求5所述的曲面共形阵列目标电磁特性分析方法,其特征在于,所述步骤51中的压缩矩阵方程为:上式中,I i 表示分解部件b i 对应的综合函数电流系数矩阵。
7.一种曲面共形阵列目标电磁特性分析系统,包括存储器和处理器,所述存储器存储有曲面共形阵列目标电磁特性分析程序,所述处理器在运行所述曲面共形阵列目标电磁特性分析程序时执行权利要求1~6任一项方法的步骤。