CN106295120B - 一种目标电磁散射特性快速计算方法 - Google Patents
一种目标电磁散射特性快速计算方法 Download PDFInfo
- Publication number
- CN106295120B CN106295120B CN201610594570.7A CN201610594570A CN106295120B CN 106295120 B CN106295120 B CN 106295120B CN 201610594570 A CN201610594570 A CN 201610594570A CN 106295120 B CN106295120 B CN 106295120B
- Authority
- CN
- China
- Prior art keywords
- target
- electromagnetic
- matrix
- subdivision
- boundary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 230000006870 function Effects 0.000 claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 39
- 239000013598 vector Substances 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 230000006698 induction Effects 0.000 claims abstract description 12
- 230000005855 radiation Effects 0.000 claims abstract description 10
- 230000003044 adaptive effect Effects 0.000 claims abstract description 9
- 230000005284 excitation Effects 0.000 claims description 11
- 230000005672 electromagnetic field Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 238000005290 field theory Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 239000002184 metal Substances 0.000 claims description 3
- 238000006880 cross-coupling reaction Methods 0.000 claims description 2
- 229940050561 matrix product Drugs 0.000 claims 1
- 230000005611 electricity Effects 0.000 abstract 2
- 230000006835 compression Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000790 scattering method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种目标电磁散射特性快速计算方法，对建立的几何外形非共形网格模型进行整体剖分，在目标剖分网格上选用电磁基函数模拟目标表面的感应电磁流，利用不连续伽辽金积分方程和边界条件构建阻抗矩阵元素，采用树形分块分组策略对目标进行分组，利用自适应交叉近似算法加速矩阵‑矢量相乘速度，采用广义最小残量法迭代计算照射源下的目标表面电磁流感应系数，开展目标表面的感应电磁流的辐射场计算，获取二次散射场特性数据。本发明适应性强、计算精度高，极大的拓展了积分方程的计算效率和计算能力，减少了剖分难度，易于工程实现。
Description
技术领域
本发明涉及一种目标电磁散射特性快速计算方法，尤其涉及一种基于自适应交叉近似加速的不连续伽辽金积分方程计算方法。
背景技术
目标电磁散射特性快速精确计算是分析目标隐身性能及反隐身设计的重要手段之一，随着进气道、机翼前沿、雷达舱、频率选择表面等复杂部件在军事目标上的大量应用，急需开展快速精确计算电大复杂目标的电磁散射特性，为武器系统参数设计和电磁波感应机理统计提供依据。如包含耶路撒冷十字频率选择(FSS)表面的某导弹雷达罩，若根据基函数尺寸网格剖分，则可得到96327个三角形面元，但基函数整体网格不足以描述频率选择表面，必须进一步细分剖分网格尺度，一个FSS单元至少需12个三角形贴片，此时雷达罩整体产生835748个三角形贴片，极大的增加了计算量。另外一种可取的剖分策略为：FSS上选用小贴片单元，雷达罩主体采用大贴片单元，此时三角形贴片个数为261853个，可在保证计算精度的情况下极大的提高计算效率。
吕志清，安翔等在2009年发表在全国微波毫米波会议优秀论文集上的非专利文献“基于自适应交叉近似的非匹配网格区域分解快速算法”给出了一种基于自适应交叉近似加速的非匹配网格的区域分解算法，并利用超松弛迭代法加速收敛过程，在贴片天线阵列、频率选择表面等周期结构电磁计算方面展示了优越的性能。但该方法是基于有限元方法给出的，计算精度不足且计算规模不大。
Zhen Peng，Kheng-Hwee Lim等在2013年发表在IEEE Transactions on Antennasand Propagation上的非专利文献“A Discontinuous Galerkin Surface IntegralEquation Method for Electromagnetic Wave Scattering From NonpenetrableTargets”公开了一种不连续伽辽金表面积分方程，详细给出了计算非共形网格电大尺寸目标的迭代求解方法，相比于传统的电磁散射积分方程求解方程，该方法利用L2型基函数处理不同部件交界处的连续性问题，保证了目标表面电磁流的连续性。但L2型基函数增加了阻抗矩阵规模，限制了目标的求解能力。
Yulei Zhao，CHuangming Tong等在2011年发表在Cross Strait Quad-RegionalRadio Science and Wireless Technology Conference上的非专利文献“The ACAAlgorithm in the Electromagnetic Scatter of Multiple Target”提出了一种计算多尺度目标的快速电磁散射方法，该方法利用纯代数压缩手段极大的提升了积分方程求解目标电磁散射问题的能力，具有很好的扩展性和实用性。但该方法针对目标整体计算，尚未考虑目标局部细节与目标主体间的对应关系。
发明内容
本发明提供一种目标电磁散射特性快速计算方法，解决了现有技术中多尺度电大复杂目标的电子散射快速仿真计算的难题，为多尺度电大尺寸目标电磁散射特性获取提供技术手段，适应性强、计算精度高，极大的拓展了积分方程的计算效率和计算能力，减少了剖分难度，易于工程实现。
为了达到上述目的，本发明提供一种目标电磁散射特性快速计算方法，包含以下步骤：
步骤S1、通过测绘获取多尺度目标的几何外形数据及部件间的相互位置关系，建立包含目标精细部件的几何外形非共形网格模型；
步骤S2、根据目标几何外形及局部细节，不同部件采用不同尺寸的剖分网格，对目标进行整体剖分；
步骤S3、在目标剖分网格上选用电磁基函数模拟目标表面的感应电磁流，利用不连续伽辽金积分方程和边界条件构建阻抗矩阵元素，建立待测目标感应电磁流与照射场间的散射关系；
电磁基函数为RWG基函数，在部件连接处选用半个基函数；
步骤S4、采用树形分块分组策略对目标进行分组，形成便于计算的阻抗矩阵方程；
步骤S5、利用自适应交叉近似算法加速矩阵-矢量相乘速度；
对树形结构，非相邻组包含的基函数的互耦矩阵Zm×n进行压缩，形成两个子矩阵Um×r,Vr×n，其中m,n为两个盒子包含的基函数个数，r为压缩后得到的矩阵维数；
步骤S6、采用广义最小残量法迭代计算照射源下的目标表面电磁流感应系数；
S7、根据目标表面的感应电磁流系数，开展目标表面的感应电磁流的辐射场计算，获取二次散射场特性数据。
所述的步骤S2中，电磁波频率范围为3.0GHz，照射波激励方向为(俯仰，方位)＝(0°,0°)，散射场方向为(俯仰，方位)＝(0～180°，180°)，间隔1°，大部件剖分网格为1/5波长，小部件细节剖分网格为1/8波长，大小部件连接处的网格是非共形的。
所述的步骤S3中，在部件内部的RWG基函数定义在两个三角形贴片的公共边上，在部件边界的贴片上定义的RWG基函数仅在一个三角形贴片上，RWG的定义形式为：
其中，r为三角形贴片内部的位置矢量，ρ为三角形的边对应的顶点到r的位置矢量，表示与基函数相关联的两个三角形贴片，±为RWG基函数对应的两个三角形贴片，如果RWG基函数对应的三角形在部件的边界上，则没有对应的三角形对形成互耦，此时RWG仅定义在一个贴片上；
利用构造的RWG基函数组构造阻抗矩阵元素，根据电磁场理论，任意金属目标的边界积分方程为：
其中，为目标表面的激励电场切向分量，k0为波数，A(r)为自由空间的矢量位函数；
根据步骤S2中的剖分策略，目标整体分为多个部件，目标边界可描述为： 为目标原始边界，Sn为每个剖分部件的边界，利用不连续伽辽金方法，每个剖分部件的感应电磁流的残差可描述为：
对上式两边利用测试函数组tn做内积运算，可得到参数的弱形式为：
同样的，可以得到每个部件边界上的磁场残差余量的弱形式为：
其中，为第m个剖分部件的切向磁场激励分量，J(r)为目标感应电流，Fk(r)为自由空间的磁矢量位函数；
为保证互相连接的部件间的电磁流连续性，在两个相连部件边界轮廓上强加边界条件：
整个矩阵运算的误差为式(4)、(5)、(6)三个部分的线性累积，通过迭代运算求解使得整个残差最小的电磁流系数即可相应的感应电磁流值。
所述的步骤S4中，采用树形结构对目标进行分组，首先找到能包围整个目标的最小长方体盒子，然后在x、y、z三个方向对等分组成八个子盒子，依次将每个子盒子细分，直至最细层盒子尺寸接近半个波长停止，分组过程中，不包含目标网格的盒子舍弃，仅保留有效盒子组。
所述的步骤S5具体包含以下步骤：
步骤S5.1、初始化第一个行索引I1，此时近似矩阵
步骤S5.2、初始化近似误差矩阵的第一行：
步骤S5.3、在第一行中找寻最大值从而确定第一个列索引J1:
步骤S5.4、得到V矩阵的第一行：
步骤S5.5、初始化近似误差矩阵的第一列：
步骤S5.6、计算U矩阵的第一列：
步骤S5.7、计算：
其中，||·||为矩阵的二范数；
步骤S5.8、判断收敛误差：如果迭代终止，否则继续；
步骤S5.9、找到第一列中的最大值作为第二个索引I2，且I2≠I1；
重复上述步骤，至到循环终止跳出，最终两个非近邻盒子形成的组行元素压缩为：
Zm×n≈Um×r×Vr×n (7)
其中，r为迭代步数；
在迭代求解运算中，矩阵和矢量的乘积表示为：
Zm×n·I＝Um×r·(Vr×n·I) (8)。
自适应交叉近似算法将矩阵计算的存储量和复杂度由O(m×n)降到O(r×(m+n))。
所述的步骤S7中，根据感应电磁流系数求解辐射场信息，计算远区电磁场值及近区相位变化，在保证求解精度的前提下极大的提升求解效率。
本发明具有以下优点：采用非共形网格对目标几何模型进行剖分，极大的简化了剖分模拟难度，对于复杂的目标可采用更加灵活的方式处理，采用RWG基函数模拟目标表面的感应电磁流，在不同部件连接处采用半个RWG基函数的形式极大的节省了计算机资源，采用自适应交叉近似算法加速迭代求解过程，不需通过复杂的场源关系及微分变化过程即可改善求解能力，计算简单，代码移植性高。
附图说明
图1为本发明提供的一种目标电磁散射特性快速计算方法的流程图。
图2为实施例中扫描得到某导弹几何模型。
图3为实施例中剖分网格与现有网格技术分布对比。
图4为实施例采用的RWG基函数结构。
图5为实施例采用的树形分组结构。
图6为实施例中目标感应电流分布与现有方法计算结果对比。
图7为实施例中目标远场RCS分布与商业软件计算对比结果。
具体实施方式
以下根据图1～图7，具体说明本发明的较佳实施例。
如图1所示，本发明提供一种目标电磁散射特性快速计算方法，包含以下步骤：
步骤S1、通过测绘获取多尺度目标的几何外形数据及部件间的相互位置关系，建立包含目标精细部件的几何外形非共形网格模型；
如图2所示，本实施例中，测绘目标为某导弹模型，目标原始尺寸(m)为：5.58×2.49×1.06，测绘仪器选用激光扫描仪，扫描精度为1mm，构建某导弹目标模型，目标几何外形建模精度为99％；
步骤S2、根据目标几何外形及局部细节，不同部件采用不同尺寸的剖分网格，对目标进行整体剖分；
将步骤S1中得到的多尺度目标几何外形进行分类区分，依据照射电磁波频率，选择足以描述目标几何、电磁模型的剖分网格尺寸对目标进行整体拟合，为降低对目标整体剖分带来的压力；
采用局部细分的方法减轻计算量，即针对不同的部件，采用不同的三角形网格进行拟合，且不同部件间连接处不必是连续的，即非共形的；
电磁波频率范围为3.0GHz，照射波激励方向为(俯仰，方位)＝(0°,0°)，散射场方向为(俯仰，方位)＝(0～180°，180°)，间隔1°，大部件剖分网格为1/5波长，小部件细节剖分网格为1/8波长，大小部件连接处的网格是非共形的，非共形网格处采用特殊的边界条件保证电磁流连续性，进而保证目标整体感应电磁流模拟的准确性和稳定性；
如图3所示，本实施例中，原始剖分网格得到261853个三角形贴片，采用本发明的方法，照射波频率为3.0GHz，目标主体剖分网格尺寸为20mm，机翼部分剖分网格为12.8mm，进气道剖分网格为10mm，机翼、进气道与主体连接处的网格不必是连续的，得到109638个三角形贴片，大大减少了贴片数量；
步骤S3、在目标剖分网格上选用合适的电磁基函数模拟目标表面的感应电磁流，利用不连续伽辽金积分方程和边界条件构建阻抗矩阵元素，建立待测目标感应电磁流与照射场间的散射关系；
本实施例中，电磁基函数为RWG基函数，在部件连接处选用半个基函数；
如图4所示，在部件内部的RWG基函数定义在两个三角形贴片的公共边上，在部件边界的贴片上定义的RWG基函数仅在一个三角形贴片上，RWG的定义形式为：
其中，r为三角形贴片内部的位置矢量，ρ为三角形的边对应的顶点到r的位置矢量，表示与基函数相关联的两个三角形贴片，±为RWG基函数对应的两个三角形贴片，如果RWG基函数对应的三角形在部件的边界上，则没有对应的三角形对形成互耦，此时RWG仅定义在一个贴片上；
利用构造的RWG基函数组构造阻抗矩阵元素，根据电磁场理论，任意金属目标的边界积分方程为：
其中，为目标表面的激励电场切向分量，k0为波数，A(r)为自由空间的矢量位函数；
根据步骤S2中的剖分策略，目标整体分为多个部件，目标边界可描述为： 为目标原始边界，Sn为每个剖分部件的边界，利用不连续伽辽金方法，每个剖分部件的感应电磁流的残差可描述为：
对上式两边利用测试函数组tn(本发明中测试函数与基函数表达形式一致)做内积运算，可得到参数的弱形式为：
同样的，可以得到每个部件边界上的磁场残差余量的弱形式为：
其中，为第m个剖分部件的切向磁场激励分量，J(r)为目标感应电流，Fk(r)为自由空间的磁矢量位函数；
为保证互相连接的部件间的电磁流连续性，在两个相连部件边界轮廓上强加边界条件：
整个矩阵运算的误差为式(4)、(5)、(6)三个部分的线性累积，通过迭代运算求解使得整个残差最小的电磁流系数即可相应的感应电磁流值；
步骤S4、采用树形分块分组策略对目标进行分组，形成便于计算的阻抗矩阵方程；
在目标剖分过程中，模拟目标几何模型的网格数是巨大的，即便当前的超级计算机也无法在短时间内完成对式(4)、(5)、(6)的迭代运算，例如对一个网格贴片数为100000的目标，其所需的存储空间为200GHz的内存空间，且随着贴片网格的增加，其所需内存空间按照平方量级增加，在分析电大尺寸目标时是非常耗时的；
如图5所示，本发明采用树形结构对目标进行分组，首先找到能包围整个目标的最小长方体盒子，然后在x、y、z三个方向对等分组成八个子盒子，依次将每个子盒子细分，直至最细层盒子尺寸接近半个波长停止，分组过程中，不包含目标网格的盒子舍弃，仅保留有效盒子组；
本实施例中，采用的树形结构为八叉树形式，最底层八叉树的边长为40mm，通过迭代形式最终形成矩阵块；
步骤S5、利用自适应交叉近似算法加速矩阵-矢量相乘速度；
自适应交叉近似算法是一种高效的纯代数矩阵压缩方法，利用压缩计算减少阻抗矩阵的存储量，本实施例中，自适应交叉近似算法采用的压缩门限为0.01，迭代收敛门限为0.003，在目标电磁场计算中，目标不同部件间的电磁波互耦随距离变化剧烈，对步骤S4中的树形结构，非相邻组包含的基函数的互耦矩阵Zm×n进行压缩，形成两个子矩阵Um×r,Vr×n，其中m,n为两个盒子包含的基函数个数，r为压缩后得到的矩阵维数；
所述的自适应交叉近似算法具体包含以下步骤：
步骤S5.1、初始化第一个行索引I1，此时近似矩阵
步骤S5.2、初始化近似误差矩阵的第一行：
步骤S5.3、在第一行中找寻最大值从而确定第一个列索引J1:
步骤S5.4、得到V矩阵的第一行：
步骤S5.5、初始化近似误差矩阵的第一列：
步骤S5.6、计算U矩阵的第一列：
步骤S5.7、计算：
其中，||·||为矩阵的二范数；
步骤S5.8、判断收敛误差：如果迭代终止，否则继续
步骤S5.9、找到第一列中的最大值作为第二个索引I2，且I2≠I1；
重复上述步骤，至到循环终止跳出，最终两个非近邻盒子形成的组行元素压缩为：
Zm×n≈Um×r×Vr×n (7)
其中，r为迭代步数；
在迭代求解运算中，矩阵和矢量的乘积表示为：
Zm×n·I＝Um×r·(Vr×n·I) (8)
显然，自适应交叉近似算法将矩阵计算的存储量和复杂度由O(m×n)降到O(r×(m+n))，特别是r往往远低于m,n，从而达到节省计算量的目的。
步骤S6、采用广义最小残量法(Gmres)迭代计算照射源下的目标表面电磁流感应系数；
采用Gmres迭代算法仅利用一次矩阵矢量相乘即可完成迭代运算，利用矩阵求解空间向残量小的方向开展，从而逐次逼近真实解；
S7、根据目标表面的感应电磁流系数，开展目标表面的感应电磁流的辐射场计算，获取二次散射场特性数据；
得到诸如RCS等近远场特性数据；
根据感应电磁流系数求解辐射场信息，计算远区电磁场值及近区相位变化，在保证求解精度的前提下极大的提升求解效率。
图6和图7给出了利用本发明的计算结果与商业软件计算结果的对比，图6中左侧为商业软件电流结果，右侧为本发明的电流计算结果，表明本发明计算的感应电流分布与商业软件一致，验证了本发明的正确性，图7中给出了通过感应电流计算的远场散射RCS值，验证了计算远区辐射场的正确性。
本发明具有以下优点：采用非共形网格对目标几何模型进行剖分，极大的简化了剖分模拟难度，对于复杂的目标可采用更加灵活的方式处理，采用RWG基函数模拟目标表面的感应电磁流，在不同部件连接处采用半个RWG基函数的形式极大的节省了计算机资源，采用自适应交叉近似算法加速迭代求解过程，不需通过复杂的场源关系及微分变化过程即可改善求解能力，计算简单，代码移植性高。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
Claims (7)
1.一种目标电磁散射特性快速计算方法，其特征在于，包含以下步骤：
步骤S1、通过测绘获取多尺度目标的几何外形数据及部件间的相互位置关系，建立包含目标精细部件的几何外形非共形网格模型；
步骤S2、根据目标几何外形及局部细节，不同部件采用不同尺寸的剖分网格，对目标进行整体剖分；
步骤S3、在目标剖分网格上选用电磁基函数模拟目标表面的感应电磁流，利用不连续伽辽金积分方程和边界条件构建阻抗矩阵元素，形成阻抗矩阵Z，建立待测目标感应电磁流与照射场间的散射关系；
电磁基函数为RWG基函数，在部件连接处选用半个基函数；
步骤S4、采用树形分块分组策略对目标进行分组，然后对树形结构中非相邻组包含的基函数的互耦矩阵Z进行压缩，形成两个子矩阵U,V的矩阵乘积，将矩阵Z与矢量I的相乘替换为U×(V×I)，提升相乘速度；
步骤S5、采用广义最小残量法迭代求解电磁场激励下的激励向量Z×I＝U×(V×I)＝b，得到未知感应电磁流的展开系数I，然后根据电磁场辐射方程计算照射源下的目标表面电磁流感应系数；
步骤S6、根据目标表面的电磁流感应系数，开展目标表面的感应电磁流的辐射场计算，获取二次散射场特性数据。
2.如权利要求1所述的目标电磁散射特性快速计算方法，其特征在于，所述的步骤S2中，网格尺寸包含：电磁波频率范围为3.0GHz，照射波激励方向为(俯仰，方位)＝(0°,0°)，散射场方向为(俯仰，方位)＝(0～180°，180°)，间隔1°，大部件剖分网格为1/5波长，小部件细节剖分网格为1/8波长，大小部件连接处的网格是非共形的。
3.如权利要求1所述的目标电磁散射特性快速计算方法，其特征在于，所述的步骤S3中，在部件内部的RWG基函数定义在两个三角形贴片的公共边上，在部件边界的贴片上定义的RWG基函数仅在一个三角形贴片上，RWG的定义形式为：
其中，r为三角形贴片内部的位置矢量，ρ为三角形的边对应的顶点到r的位置矢量，表示与基函数相关联的两个三角形贴片，±为RWG基函数对应的两个三角形贴片，如果RWG基函数对应的三角形在部件的边界上，则没有对应的三角形对形成互耦，此时RWG仅定义在一个贴片上；
利用构造的RWG基函数组构造阻抗矩阵元素，根据电磁场理论，任意金属目标的边界积分方程为：
其中，为目标表面的激励电场切向分量，k0为波数，A(r)为自由空间的矢量位函数；
根据步骤S2中的剖分策略，目标整体分为多个部件，目标边界可描述为： 为目标原始边界，Sn为每个剖分部件的边界，利用不连续伽辽金方法，每个剖分部件的感应电磁流的残差可描述为：
对上式两边利用测试函数组tn做内积运算，可得到参数的弱形式为：
同样的，可以得到每个部件边界上的磁场残差余量的弱形式为：
其中，为第m个剖分部件的切向磁场激励分量，J(r)为目标感应电流，Fk(r)为自由空间中对应基函数序列号k的磁矢量位函数；
为保证互相连接的部件间的电磁流连续性，在两个相连部件边界轮廓上强加边界条件：
整个矩阵运算的误差为式(4)、(5)、(6)三个部分的线性累积，通过迭代运算求解使得整个残差最小的电磁流感应系数即相应的感应电磁流值。
4.如权利要求1所述的目标电磁散射特性快速计算方法，其特征在于，所述的步骤S4中，采用树形结构对目标进行分组，首先找到能包围整个目标的最小长方体盒子，然后在x、y、z三个方向对等分组成八个子盒子，依次将每个子盒子细分，直至最细层盒子尺寸接近半个波长停止，分组过程中，不包含目标网格的盒子舍弃，仅保留有效盒子组。
5.如权利要求1所述的目标电磁散射特性快速计算方法，其特征在于，所述的步骤S4具体包含以下步骤：
步骤S4.1、初始化第一个行索引I1，此时近似矩阵
步骤S4.2、初始化近似误差矩阵的第一行：
步骤S4.3、在第一行中找寻最大值从而确定第一个列索引J1:
步骤S4.4、得到V矩阵的第一行：
步骤S4.5、初始化近似误差矩阵的第一列：
步骤S4.6、计算U矩阵的第一列：
步骤S4.7、计算：
其中，||·||为矩阵的二范数；
步骤S4.8、判断收敛误差：如果迭代终止，否则继续；
步骤S4.9、找到第一列中的最大值作为第二个索引I2，且I2≠I1；
重复上述步骤，至到循环终止跳出，最终两个非近邻盒子形成的组行元素压缩为：
Zm×n≈Um×k×Vk×n (7)
其中，k为迭代步数；
在迭代求解运算中，矩阵和矢量的乘积表示为：
Zm×n·I＝Um×k·(Vk×n·I) (8)。
6.如权利要求5所述的目标电磁散射特性快速计算方法，其特征在于，自适应交叉近似算法将矩阵计算的存储量和复杂度由O(m×n)降到O(k ×(m+n))。
7.如权利要求1所述的目标电磁散射特性快速计算方法，其特征在于，所述的步骤S6中，根据感应电磁流感应系数求解辐射场信息，计算远区电磁场值及近区相位变化，在保证求解精度的前提下极大的提升求解效率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594570.7A CN106295120B (zh) | 2016-07-26 | 2016-07-26 | 一种目标电磁散射特性快速计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594570.7A CN106295120B (zh) | 2016-07-26 | 2016-07-26 | 一种目标电磁散射特性快速计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106295120A CN106295120A (zh) | 2017-01-04 |
| CN106295120B true CN106295120B (zh) | 2019-10-18 |
Family
ID=57652846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610594570.7A Active CN106295120B (zh) | 2016-07-26 | 2016-07-26 | 一种目标电磁散射特性快速计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106295120B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090255A (zh) * | 2017-11-29 | 2018-05-29 | 上海无线电设备研究所 | 一种计算等离子体覆盖目标电磁散射的高低频混合算法 |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107526855B (zh) * | 2016-06-22 | 2021-10-08 | 南京理工大学 | 分析不确定性等离子体特性的不连续伽辽金时域有限元法 |
| CN106844919B (zh) * | 2017-01-11 | 2020-09-01 | 三峡大学 | 一种输电线路散射场的求解方法 |
| CN107330162A (zh) * | 2017-06-13 | 2017-11-07 | 电子科技大学 | 针对非曲线边界二维模型的网格划分及信息获取方法 |
| CN107967377B (zh) * | 2017-11-02 | 2021-01-19 | 上海无线电设备研究所 | 基于蒙特卡罗方法的非均匀等离子体电磁特性仿真方法 |
| CN108051792B (zh) * | 2017-12-05 | 2021-07-20 | 上海无线电设备研究所 | 目标与粗糙面耦合散射源的分布描述方法 |
| CN108052738B (zh) * | 2017-12-13 | 2021-10-15 | 电子科技大学 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
| CN108170647B (zh) * | 2017-12-13 | 2019-12-27 | 南京理工大学 | 一种低频电磁特性分析的自适应嵌套交叉近似方法 |
| CN109376485A (zh) * | 2018-12-03 | 2019-02-22 | 上海无线电设备研究所 | 基于aca-mlfma加速的区域分解非共形网格的快速仿真建模方法 |
| CN109765538B (zh) * | 2019-01-11 | 2022-12-02 | 西安电子科技大学 | 非均匀介质目标体的散射场确定方法 |
| CN109992819A (zh) * | 2019-01-25 | 2019-07-09 | 南京师范大学镇江创新发展研究院 | 一种基于自适应交叉近似算法的电缆串扰处理方法 |
| CN110059422B (zh) * | 2019-04-23 | 2023-01-20 | 北京环境特性研究所 | 频率选择表面曲面天线罩的电磁散射特性仿真方法 |
| CN111274703B (zh) * | 2020-01-20 | 2021-10-01 | 中国人民解放军国防科技大学 | 一种复杂结构电大尺寸目标电磁特性分析方法及系统 |
| CN112364362B (zh) * | 2020-11-16 | 2023-12-29 | 宁波九寰适创科技有限公司 | 一种面向流体仿真方向的并行多层自适应局部加密方法 |
| CN112733347B (zh) * | 2020-12-31 | 2022-09-02 | 中国人民解放军空军工程大学 | 一种导体-介质复合目标与环境电磁散射快速计算方法 |
| CN112949079B (zh) * | 2021-03-19 | 2023-07-18 | 西安电子科技大学 | 金属目标散射的快速仿真方法 |
| CN113792257B (zh) * | 2021-09-18 | 2024-01-12 | 南京航空航天大学 | 一种基于mbrwg和网格自适应加密的电磁散射求解方法 |
| CN117150824B (zh) * | 2023-10-31 | 2024-03-19 | 南京理工大学 | 基于部分重叠rwg基函数的复杂目标电磁散射分析方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
| CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
| CN105224780A (zh) * | 2014-06-28 | 2016-01-06 | 南京理工大学 | 分析导体瞬态电磁散射特性的时域高阶Nystrom方法 |
| CN105303022A (zh) * | 2014-07-28 | 2016-02-03 | 南京辉宏电子信息有限公司 | 快速获取目标电磁散射特性的高斯波束方法 |
| CN105608311A (zh) * | 2015-12-17 | 2016-05-25 | 三峡大学 | 基于高阶Radau方法的高压输电线路电磁暂态快速计算方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6904374B2 (en) * | 2003-09-29 | 2005-06-07 | The Boeing Company | Methods and systems for predicting electromagnetic scattering |
-
2016
- 2016-07-26 CN CN201610594570.7A patent/CN106295120B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
| CN105224780A (zh) * | 2014-06-28 | 2016-01-06 | 南京理工大学 | 分析导体瞬态电磁散射特性的时域高阶Nystrom方法 |
| CN105303022A (zh) * | 2014-07-28 | 2016-02-03 | 南京辉宏电子信息有限公司 | 快速获取目标电磁散射特性的高斯波束方法 |
| CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
| CN105608311A (zh) * | 2015-12-17 | 2016-05-25 | 三峡大学 | 基于高阶Radau方法的高压输电线路电磁暂态快速计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 隐身目标电磁散射特性的图形电磁计算方法研究;钟全华 等;《宇航计测技术》;20070831;第27卷(第4期);第54-57、64页 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090255A (zh) * | 2017-11-29 | 2018-05-29 | 上海无线电设备研究所 | 一种计算等离子体覆盖目标电磁散射的高低频混合算法 |
| CN108090255B (zh) * | 2017-11-29 | 2021-08-31 | 上海无线电设备研究所 | 一种计算等离子体覆盖目标电磁散射的高低频混合方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106295120A (zh) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106295120B (zh) | 一种目标电磁散射特性快速计算方法 | |
| Pan et al. | A fast algorithm for multiscale electromagnetic problems using interpolative decomposition and multilevel fast multipole algorithm | |
| CN105930567A (zh) | 一种基于子区域自适应积分的电磁散射特性获取方法 | |
| CN104112051B (zh) | 基于单元形函数的反射面天线机电集成设计方法 | |
| CN106529082A (zh) | 一种快速计算电大尺寸目标电磁散射特征的方法 | |
| CN113158527B (zh) | 一种基于隐式fvfd计算频域电磁场的方法 | |
| CN102708229A (zh) | 复杂分层媒质结构的矩阵分解结合新奇异值分解方法 | |
| CN112733364A (zh) | 一种基于阻抗矩阵分块的箔条云散射快速计算方法 | |
| Taygur et al. | A bidirectional ray-tracing method for antenna coupling evaluation based on the reciprocity theorem | |
| CN112329285A (zh) | 泰勒级数展开的混合算法分析多尺度结构瞬态响应方法 | |
| CN112329204A (zh) | 考虑载体平台耦合的重复性结构电磁特性特征模快速分析方法 | |
| CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
| Na et al. | A Multiple Huygens Surface Based Ray Tracing Framework with GPU Acceleration | |
| Jing et al. | Experiments with Lanczos biconjugate A-orthonormalization methods for MoM discretizations of Maxwell's equations | |
| CN103246827B (zh) | 复杂外形金属目标电磁散射的无网格仿真方法 | |
| Taygur et al. | Computation of antenna transfer functions with a bidirectional ray-tracing algorithm utilizing antenna reciprocity | |
| Lai et al. | Analysis of nurbs surfaces modeled geometries with higher-order MoM based aim | |
| CN111368398B (zh) | 一种不确定结构电大目标电磁散射特性分析方法及装置 | |
| Ozgun et al. | Iterative leap-field domain decomposition method: A domain decomposition finite element algorithm for 3D electromagnetic boundary value problems | |
| Wang et al. | An effective MoM solution with nested complex source beam method for electromagnetic scattering problems | |
| Cáez | Fast and Efficient Electromagnetic Simulation of Electrically Extra-Large Problems Using Phase Information and Mesh Automation | |
| Ding et al. | Multiresolution preconditioner combined with matrix decomposition–singular value decomposition algorithm for fast analysis of electromagnetic scatterers with dense discretisations | |
| Gu et al. | Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts | |
| Zhang et al. | Comparison of Satellite Target Recognition Methods | |
| Garcıa et al. | Analysis of the parameters of an approach that combines the characteristic basis function method and the multilevel fast multipole |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |