[发明专利]一种曲面共形阵列目标电磁特性分析方法及系统

权利要求书

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]^TZ_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]^TZ_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]^TV_i，其中，v_i为照射波对于分解部件b_i的压缩激励矩阵。

2.如权利要求1所述的曲面共形阵列目标电磁特性分析方法，其特征在于，所述步骤1包括：

步骤11，将目标按照阵列结构和阵元分布特征分解为若干分解部件并分别编号；

步骤12，对每个分解部件进行非结构化三角网格剖分；

步骤13，在分解部件外部构建一个虚拟三维封闭表面，对构建的虚拟三维封闭表面进行非结构化三角网格剖分；

步骤14，根据分解部件自身的三角剖分网格和构建的虚拟三维封闭表面的三角剖分网格分别定义RWG函数，其中，针对虚拟三维封闭曲面定义的RWG函数即是外部等效源。