CN102129523B - 基于mda和mlssm的分析复杂目标电磁散射的方法 - Google Patents
基于mda和mlssm的分析复杂目标电磁散射的方法 Download PDFInfo
- Publication number
- CN102129523B CN102129523B CN 201110070211 CN201110070211A CN102129523B CN 102129523 B CN102129523 B CN 102129523B CN 201110070211 CN201110070211 CN 201110070211 CN 201110070211 A CN201110070211 A CN 201110070211A CN 102129523 B CN102129523 B CN 102129523B
- Authority
- CN
- China
- Prior art keywords
- mlssm
- cube
- mda
- matrix
- target
- 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
Images
Landscapes
- Aerials With Secondary Devices (AREA)
Abstract
本发明公开了一种基于MDA和MLSSM的分析复杂目标电磁散射的方法。该方法首先对待求物体建立八叉树结构,通过MDA方法来获得离散电场积分方程形成的稠密阻抗矩阵中对应于远场作用部分的稀疏低秩表示,然后将阻抗矩阵中低秩表示的部分用改进的MLSSM表示出来。最后基于改进的MLSSM表示,构造一种适用于迭代求解的快速计算矩阵矢量乘的算法。本发明能够大大减少用电场积分方程方法分析复杂目标电磁散射的计算时间和内存需求。
Description
技术领域
本发明涉及一种积分方程方法求解电磁问题里加速矩阵矢量乘计算的快速技术,特别是一种基于MDA和MLSSM的分析复杂目标电磁散射的方法。
背景技术
近年来,分析任意形状三维目标的电磁散射特性,已经受到越来越多国内外研究学者的关注。其应用主要表现在:雷达系统的设计与目标识别、军用武器的隐身与反隐身、复杂环境中的电磁兼容问题等等。由于雷达工作在微波频段,常见军用目标,如战车、导弹、飞机等等不但具有复杂的几何结构,而且具有超大的电尺寸,增加了散射特性分析的复杂度。因此,如何高效计算这类大型复杂目标的电磁散射特性,如目标的雷达截面积RCS 参数的提取,一直是从事雷达总体设计以及隐身与反隐身研究的学者、工程师们所共同关心的问题。
对于目标电磁散射特性的分析,传统的解析方法仅局限于解决几何形状为球体或无限长圆柱体等简单目标。随着计算机技术与计算电磁学理论的迅猛发展,出现了许多可用于分析任意形状复杂目标电磁散射特性的数值方法。这些方法大致可划分为基于高频近似的高频RCS预估方法和基于微分方法或积分方程的严格数值方法两大类。
基于高频近似的高频方法,如几何光学法(GO)、物理光学法(PO)、物理绕射理论(PTD)、弹跳射线方法(SBR)等等,由于具有计算速度快,消耗机器内存少的优点被广泛应用于各类复杂目标电磁散射特性的分析。然而,由于高频近似条件的引入,使得高频方法的计算结果精度低,这是制约高频方法发展及其应用的一个重要的因素。用于散射特性分析的微分方程法,包括时域有限差分法(FDTD)和有限元法(FEM),虽然得到稀疏矩阵,但是难以精确拟合目标的曲面边界。并且需要在传播空间内进行网格划分,存在网格数值色散。对于开域问题,微分方程的求解必须施加吸收边界条件,常用于分析目标的电磁散射特性。积分方程法又可分为体积分方程法和表面积分方程法两种。体积分方程法适合于非均匀介质材料的电磁散射特性分析,而表面积分方程法适合于理想导体或均匀介质材料的电磁散射特性分析。
尽管基于积分方程的矩量法具有严格的理论模型,但是其生成的矩阵为满阵。假定N为未知量个数,则存储该稠密矩阵将要耗费O(N2)的内存量。同时,如果利用直接法来求解矩量法的阻抗矩阵方程,其计算复杂度为O(N2)。对于电大尺寸的目标, 问题变得越发严重。正是由于受到计算机存储量的限制,长期以来矩量法仅仅局限于分析低频区到谐振区范围内目标的电磁散射特性。近几十年来,由于计算机技术与计算电磁学理论的迅猛发展,矩量法(MoM)作为一种严格的数值方法,重新引起了人们极大的关注。特别是各种高效的快速积分方法的提出,大幅度地降低了矩量法分析所需的计算量以及内存量。其中自由空间里多层快速多极子(MLFMA)是其中最有代表性的快速算法。MLFMA是依赖于Green函数,因此应用到多层媒质的情形,MLFMA就需要做修改。而基于矩阵低秩分解的代数类快速算法是不依赖于Green函数,因此很容易直接应用到多层媒质的情形。低秩类方法主要有ACA,UV,MDA。其中MDA主要利用等效源的思想,最初是应用在二维的情况(E.Michielssen and A. Boag, “A multilevel matrix decomposition algorithm for analyzing scattering from largestructures”, IEEE Trans. Antennas Propag., vol. 44,no. 8, pp. 1086-1093, August 1996.),后来由J.M. Rius推广到三维的情况(J.M. Rius, J. Parr′on, E. beda, J.R.Mosig, “Multilevel Matrix Decomposition Algorithm forAnalysis of Electrically Large Electromagnetic Problemsin 3-D”, Microwave and Optical Technology Letters,Vol. 22, No.3, pp. 177-182, 5th August 1999.)。相对于其它方法而言,MDA在计算平面结构的目标在计算时间和内存需求上有优势。
由于积分方程方法得到的阻抗矩阵是稠密的,F. X. Canning提出Multilevel Simple Sparse Method (MLSSM) 矩阵稀疏化表示的方法。最初是F. X. Canning 在利用直接方法来求解电场积分方程时存储阻抗矩阵 (1. F. X. Canning and K. Rogovin, "Auniversal matrix solver for integral-equation based problems," IEEE Antennas and Propagation Magazine, vol.45, February, 2003; 2. A. Zhu, R. J. Adams, and F.X. Canning, "Modified simply sparse method for electromagnetic scattering by PEC," IEEE Antennas and Propagation Society International Symposium, Washington, DC, 2005, pp. 427-430.) 。随着目标电尺寸增大,待求的未知量增多,迭代求解过程中矩阵矢量乘的计算复杂度是求解的关键技术。基于MLSSM稀疏表示的迭代解法还没有报道。
发明内容
本发明所解决的技术问题在于提供一种基于MDA和MLSSM的分析复杂目标电磁散射的方法。
实现本发明目的的技术解决方案为:一种基于MDA和MLSSM的分析复杂目标电磁散射的方法,包括以下步骤:
步骤1、对入射平面波参数进行设置,根据实际物理背景,对入射平面波的频率和入射角度进行设置,入射平面波的频率为微波频段,入射角度在球坐标下为(θ,φ),其中θ∈[0°,180°], φ∈[0°,360°];
步骤2、利用Ansys软件建立待求目标的几何模型,之后利用平面三角形进行离散剖分;
步骤3、采用八叉树结构对剖分后的目标模型进行分组,用一个立方体将目标体包围住,该立方体定义为第零层的第一个且是最后一个组结点,把该立方体等分为八个子立方体结点形成第一层组结点,然后再对每个子立方体进行与上一步相同的细分,并以此类推直到最底层立方体的电尺寸达到所需合适的大小为止;
步骤4、根据步骤3得到的八叉树结构,首先进行Morton编码,然后将每个立方体相邻的立方体组设定为近场区组,之后设定远场区,所述远场区为包含该立方体的父层立方体区域的近场区域中除掉本层区域的近场区组;
步骤5、根据立方体电尺寸的大小,在其表面设置等效源,利用等效源对阻抗矩阵的远作用子块进行低秩表示,并直接计算阻抗矩阵的近作用子块;
步骤6、利用QR和SVD对步骤5得到的低秩子块表示进一步压缩处理,得到阻抗矩阵改进的MLSSM表示;
步骤7、根据步骤6得到改进的MLSSM表示,构造迭代求解的快速矩阵矢量乘算法;
步骤8、求解步骤7中的阻抗矩阵方程,得到模型表面电流系数,并且根据电流系数得到目标表面的电流分布,从而得到目标的电磁散射特性参数RCS。
本发明与现有技术相比,其显著优点为:1)利用本方法最终形成的MLSSM表示降低了计算机内存的需求,从而在有限的计算机资源下处理更大点尺寸问题;2)在迭代法求解的时候,能实现快速的矩阵矢量乘,从而减少迭代求解时间。因此,本发明可以在有限计算资源下更快处理更大的问题。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明基于MDA和MLSSM的快速矩阵矢量乘技术的流程图。
图2是sphere-cone在Ansys下的剖分示意图。
图3是三维的八叉树结构的建立过程示意图,从左至右分别是第零层,第一层,第二层。
图4 阻抗矩阵的示意图。
图5 根据等效原理建立远场矩阵子块的低秩表示的过程示意图。
图7是sphere-cone模型随电尺寸变化的矩阵矢量乘时间曲线示意图。
图8是sphere-cone模型随电尺寸变化的内存需求曲线示意图。
图9是入射波频率为8 GHz ,入射角是θ=0°,φ=0°时, sphere-cone模型的双站RCS曲线图。
具体实施方式
结合图1,本发明的一种基于MDA和MLSSM的分析复杂目标电磁散射的方法,包括以下步骤:
步骤1、对入射平面波参数进行设置,根据实际物理背景,对入射平面波的频率和入射角度进行设置;入射平面波的频率为微波频段,入射角度在球坐标下为(θ,φ),其中θ∈[0°,180°], φ∈[0°,360°]。
步骤2、利用Ansys软件建立待求目标的几何模型,之后利用平面三角形进行离散剖分;利用平面三角形进行离散剖分时,采用基于RWG基函数的三角形网格对复杂目标表面进行剖分,每平方电波长内剖分的三角形数目大于120。
步骤3、采用八叉树结构对剖分后的目标模型进行分组,用一个立方体将目标体包围住,该立方体定义为第零层的第一个且是最后一个组结点,把该立方体等分为八个子立方体结点形成第一层组结点,然后再对每个子立方体进行与上一步相同的细分,并以此类推直到最底层立方体的电尺寸达到所需合适的大小为止;最底层立方体的电尺寸为0.2~0.4个入射波波长。
步骤4、根据步骤3得到的八叉树结构,首先进行Morton编码,然后将每个立方体相邻的立方体组设定为近场区组,之后设定远场区,所述远场区为包含该立方体的父层立方体区域的近场区域中除掉本层区域的近场区组;
步骤5、根据立方体电尺寸的大小,在其表面设置等效源,利用等效源对阻抗矩阵的远作用子块进行低秩表示,并直接计算阻抗矩阵的近作用子块;所述等效源的个数Q按照以下标准确定:
其中Ss,So,d分别是远作用的时候场所在的立方体和源所在的立方体的电尺寸,以及它们之间的距离,k是入射波的波数,每一个等效源就是一个RWG基函数,而且这些等效源均匀分布在立方体的表面。利用等效源对阻抗矩阵的远作用子块进行低秩表示具体为:
其中[Zmn]∈Cm×n是阻抗矩阵中对应于远场作用的子块,[U]∈Cm×r,[V]∈Cn×r,m、n、r为正整数,其中r小于min(m,n)。
步骤6、利用QR和SVD对步骤5得到的低秩子块表示进一步压缩处理,得到阻抗矩阵改进的MLSSM表示;阻抗矩阵改进的MLSSM表示为:
步骤7、根据步骤6得到改进的MLSSM表示,构造迭代求解的快速矩阵矢量乘算法;构造迭代求解的快速矩阵矢量乘算法按照以下函数执行:
步骤8、求解步骤7中的阻抗矩阵方程,得到模型表面电流系数,并且根据电流系数得到目标表面的电流分布,从而得到目标的电磁散射特性参数RCS。
下面以图2所示sphere-cone结构为例,对本发明的具体步骤作进一步详细描述。
根据本发明分析图2所示sphere-cone 结构的基于Matrix Decomposition Algorithm和Multilevel Simple Sparse Method的分析复杂目标电磁散射的快速方法,其具体操作步骤如下:
第一步,入射波参数的设置。根据实际物理背景,设置入射波为平面波。这里为了验证我们快速方法的有效性, 对于相同的sphere-cone结构, 入射波的频率设置为2 GHz, 4 GHz, 8 GHz, 得到电尺寸逐渐增大的结构。入射角度是θ=0°,φ=0°
第二步,目标的几何建模。根据目标模型的几何参数,用Ansys建立目标的几何模型。这里建立的模型的时候按十分之一波长来剖分sphere-cone。
第三步,采用八叉树结构对剖分后的sphere-cone结构进行分组,先用一个立方体将第一层sphere-cone包围住,该立方体定义为第零层的第一个且是最后一个组结点,把该立方体等分为八个子立方体结点形成第一层组结点,然后再对每个子立方体进行与上一步相同的细分,并以此类推直到最底层立方体的尺寸达到所需要求为止,一般为0.2~0.4个电波长。
第四步,对第三步建立的树形结构采用Morton编码。首先查找每一层非空的小正方体,其次通过树形结构得到每一个非空小正方体近场作用和远场作用小正方体的信息。这里我们采用Morton编码可以快速进行索引,从而提高计算效率。建立树形结构的目的是为了确定阻抗矩阵中需要满秩填充(近场作用)和低秩填充(远场作用)的子块。
第五步,利用矩阵分解算法(Matrix Decomposition Algorithm)对阻抗矩阵的远场作用子块进行填充。根据等效原理(如图5所示),可以对远场作用子块表示为以下形式:
其中Zmn是待填充的子块,维数是m×n;Zmp,Zqp,Zqn维数分别是m×p,p×q,q×n;p,q是等效源的数目,m,n>>p,q。也就是将待填充的子块Zmn用Zmp,Zqp,Zqn的乘积来低秩表示。
第六步,利用QR分解(QR Decomposition)和奇异值分解(Singular Value Decomposition —— SVD分解)对上面得到的Zmn低秩表示进一步压缩。具体过程如下:
做SVD分解有
做QR分解有
从而有
因此可以得到
这样就得到压缩后的低秩表示矩阵。
第七步,利用上面得到远作用子块的低秩表示,得到阻抗矩阵Z改进的Multilevel Simple Sparse Method表示如下:
构造改进的Multilevel Simple Sparse Method具体过程如下:
(1) 先对L层每一个非空场立方体Li,根据八叉树结构找到l(3≤l≤L)层的远场作用源立方体lj。而每个源立方体对应阻抗矩阵的一个低秩表示子块U(Li,lj)ω(Li,lj)V(Li,lj),然后将U(Li,lj)ω(Li,lj)排成一行形成矩阵Z(Li,l)。
(2) 将(1)得到的矩阵Z(Li,l)排成一行形成矩阵Z(Li)。对Z(Li)做SVD操作得到Z(Li)=U(Li)S(Li)VH(Li)。因此,对L层每一个非空场立方体Li对应了一个矩阵U(Li)。将这些子块按对角线排列就形成了块对角的酉矩阵UL。令VL=(UL)T就得到了VL。
(3) 令ZL-1(Li,lj)表示ZL-1中与场立方体Li和源立方体lj相联系的子块,其中 l(3≤l≤L)。按照以下方式得到 。
(4) 对于l(3≤l≤L-1), 按照上述(1)-(3)同样的步骤就可以得到, Ul , Vl。在表示完成之后,只需要存储Ul, (2≤l≤L)。图5显示第四层和第三层改进的MLSSM表示中的U4,,V4,U3,,V3示意图。根据上面得到改进的MLSSM表示,构造以下快速矩阵矢量乘算法。
其中x,y分别为输入向量和输出向量。如图7-8所示,相比MDA-SVD,现有的改进的MLSSSM表示在内存需求和矩阵矢量乘的时间有很大的优势。
第八步,利用广义最小余量法求解线性方程组的解得到电流系数。根据电流系数得到目标的散射参数RCS。
Claims (7)
1.一种基于MDA和MLSSM的分析复杂目标电磁散射的方法,其特征在于,包括以下步骤:
步骤1、对入射平面波参数进行设置,根据实际物理背景,对入射平面波的频率和入射角度进行设置,入射平面波的频率为微波频段,入射角度在球坐标下为(θ,φ),其中θ∈[0°,180°], φ∈[0°,360°];
步骤2、利用Ansys软件建立待求目标的几何模型,之后利用平面三角形进行离散剖分;
步骤3、采用八叉树结构对剖分后的目标模型进行分组,用一个立方体将目标体包围住,该立方体定义为第零层的第一个且是最后一个组结点,把该立方体等分为八个子立方体结点形成第一层组结点,然后再对每个子立方体进行与上一步相同的细分,并以此类推直到最底层立方体的电尺寸达到所需合适的大小为止;
步骤4、根据步骤3得到的八叉树结构,首先进行Morton编码,然后将每个立方体相邻的立方体组设定为近场区组,之后设定远场区,所述远场区为包含该立方体的父层立方体区域的近场区域中除掉本层区域的近场区组;
步骤5、根据立方体电尺寸的大小,在其表面设置等效源,利用等效源对阻抗矩阵的远作用子块进行低秩表示,并直接计算阻抗矩阵的近作用子块;
步骤6、利用QR和SVD对步骤5得到的低秩子块表示进一步压缩处理,得到阻抗矩阵改进的MLSSM表示;
步骤7、根据步骤6得到改进的MLSSM表示,构造迭代求解的快速矩阵矢量乘算法;
步骤8、求解步骤7中的阻抗矩阵方程,得到模型表面电流系数,并且根据电流系数得到目标表面的电流分布,从而得到目标的电磁散射特性参数RCS。
2.根据权利要求1所述的基于MDA和MLSSM的分析复杂目标电磁散射的方法,其特征在于,步骤2中利用平面三角形进行离散剖分时,采用基于RWG基函数的三角形网格对复杂目标表面进行剖分,每平方电波长内剖分的三角形数目大于120。
3.根据权利要求1所述的基于MDA和MLSSM的分析复杂目标电磁散射的方法,其特征在于,步骤3中最底层立方体的电尺寸为0.2~0.4个入射波波长。
4.根据权利要求1所述的基于MDA和MLSSM的分析复杂目标电磁散射的方法,其特征在于,步骤5中等效源的个数Q按照以下标准确定:
其中Ss,So,d分别是远作用的时候场所在的立方体和源所在的立方体的电尺寸,以及它们之间的距离,k是入射波的波数,每一个等效源就是一个RWG基函数,而且这些等效源均匀分布在立方体的表面。
5.根据权利要求1所述的基于MDA和MLSSM的分析复杂目标电磁散射的方法,其特征在于,步骤5中利用等效源对阻抗矩阵的远作用子块进行低秩表示具体为:
其中[Zmn] ∈Cm×n是阻抗矩阵中对应于远场作用的子块,[U] ∈Cm×r ,[V] ∈Cn×r,m、n、r为正整数,其中r小于min(m,n)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110070211 CN102129523B (zh) | 2011-03-23 | 2011-03-23 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110070211 CN102129523B (zh) | 2011-03-23 | 2011-03-23 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102129523A CN102129523A (zh) | 2011-07-20 |
CN102129523B true CN102129523B (zh) | 2013-04-24 |
Family
ID=44267606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110070211 Active CN102129523B (zh) | 2011-03-23 | 2011-03-23 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102129523B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150415B (zh) * | 2012-12-14 | 2016-02-24 | 中国科学院深圳先进技术研究院 | 微带电路全波分析方法 |
CN103176082B (zh) * | 2013-03-01 | 2015-06-24 | 中国舰船研究设计中心 | 基于高低频算法校模的电特大尺寸目标散射预估方法 |
CN103218487B (zh) * | 2013-04-09 | 2016-08-31 | 南京理工大学 | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 |
CN103246827B (zh) * | 2013-05-31 | 2016-09-21 | 南京理工大学 | 复杂外形金属目标电磁散射的无网格仿真方法 |
CN103439698B (zh) * | 2013-08-27 | 2015-06-10 | 西北工业大学 | 获取雷达散射面积的方法 |
CN103530469B (zh) * | 2013-10-23 | 2016-03-30 | 武汉大学 | 一种角反射器雷达截面积的计算方法 |
CN105277927B (zh) * | 2014-06-19 | 2018-09-04 | 南京理工大学 | 飞行器编队瞬态电磁特性时域阶数步进分析方法 |
CN106156394B (zh) * | 2015-04-22 | 2019-06-25 | 南京理工大学 | 基于显式差分格式的电磁特性提取方法 |
CN106202599B (zh) * | 2015-05-07 | 2019-12-27 | 南京理工大学 | 加速求解混合目标电磁散射特性的多层复点源方法 |
CN106202595B (zh) * | 2015-05-07 | 2019-12-27 | 南京理工大学 | 分析介质目标瞬态电磁散射特性的时域非共形网格方法 |
CN105184010A (zh) * | 2015-09-25 | 2015-12-23 | 天津城建大学 | 基于快速多极间接边界元法的高频地震波散射模拟方法 |
CN106066941B (zh) * | 2016-06-08 | 2018-10-19 | 南京航空航天大学 | 一种基于cbfm和smw算法的电磁散射快速分析方法 |
CN106295120B (zh) * | 2016-07-26 | 2019-10-18 | 上海无线电设备研究所 | 一种目标电磁散射特性快速计算方法 |
CN106503349B (zh) * | 2016-10-25 | 2019-08-13 | 中国人民解放军国防科学技术大学 | 一种类周期结构目标电磁散射特性快速计算方法 |
CN107884754A (zh) * | 2017-11-06 | 2018-04-06 | 电子科技大学 | 雷达散射面有效面积检测装置与方法 |
JP6956958B2 (ja) * | 2017-11-29 | 2021-11-02 | 富士通株式会社 | 推定プログラム、推定装置および推定方法 |
CN107943755B (zh) * | 2017-12-12 | 2020-12-04 | 电子科技大学 | 复杂色散媒质中指数时间积分法的矩阵指数降维方法 |
CN108230414A (zh) * | 2017-12-13 | 2018-06-29 | 南京理工大学 | 一种低高频混合多尺度目标低秩压缩分析方法 |
CN112329203B (zh) * | 2020-10-08 | 2022-09-06 | 南京理工大学 | 宽带宽角扫描相控阵天线阵元失效后的方向图重构及修复方法 |
CN112329204B (zh) * | 2020-10-11 | 2022-08-16 | 南京理工大学 | 考虑载体平台耦合的重复性结构电磁特性特征模快速分析方法 |
CN112949079B (zh) * | 2021-03-19 | 2023-07-18 | 西安电子科技大学 | 金属目标散射的快速仿真方法 |
CN115935802B (zh) * | 2022-11-23 | 2023-08-29 | 中国人民解放军军事科学院国防科技创新研究院 | 电磁散射边界元计算方法、装置、电子设备及存储介质 |
CN116430350B (zh) * | 2023-06-14 | 2023-09-26 | 南京理工大学 | 一种空间动态群目标属性散射中心的建模方法及系统 |
CN116956472B (zh) * | 2023-09-20 | 2024-01-05 | 中国空气动力研究与发展中心空天技术研究所 | 一种mlfma伴随求解的rcs表面敏感度计算方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6947869B2 (en) * | 2002-03-29 | 2005-09-20 | The United States Of America As Represented By The Secretary Of The Navy | Efficient near neighbor search (ENN-search) method for high dimensional data sets with noise |
FR2914750B1 (fr) * | 2007-04-03 | 2009-07-03 | Commissariat Energie Atomique | Methode d'estimation de surface equivalente radar |
CN101923166B (zh) * | 2009-06-10 | 2012-11-14 | 中国科学院电子学研究所 | 一种三维复杂目标的合成孔径雷达图像仿真方法 |
-
2011
- 2011-03-23 CN CN 201110070211 patent/CN102129523B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102129523A (zh) | 2011-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102129523B (zh) | 基于mda和mlssm的分析复杂目标电磁散射的方法 | |
CN102081690B (zh) | 复杂电路的矩阵分解结合新奇异值分解方法 | |
CN102033985A (zh) | 基于*-矩阵算法的高效时域电磁仿真方法 | |
CN102708229A (zh) | 复杂分层媒质结构的矩阵分解结合新奇异值分解方法 | |
CN103425864B (zh) | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 | |
CN102156764A (zh) | 一种分析天线辐射和电磁散射的多分辨预条件方法 | |
CN106529082A (zh) | 一种快速计算电大尺寸目标电磁散射特征的方法 | |
CN113158527B (zh) | 一种基于隐式fvfd计算频域电磁场的方法 | |
CN106294894B (zh) | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 | |
CN103425816A (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
CN102054094A (zh) | 一种平面微带电路的快速多层方向性仿真方法 | |
Burkholder et al. | Fast dual-MGS block-factorization algorithm for dense MoM matrices | |
CN108038313A (zh) | 一种剖分不均匀的目标电磁散射特性的分析方法 | |
Pan et al. | Hierarchical interpolative decomposition multilevel fast multipole algorithm for dynamic electromagnetic simulations | |
CN103235193A (zh) | 毫米波段内卫星电磁散射特性的数值方法 | |
CN106066941A (zh) | 一种基于cbfm和smw算法的电磁散射快速分析方法 | |
CN113919128A (zh) | 一种适用于隐身灵敏度计算的电磁变分方法 | |
Shao et al. | Fast simulation of array structures using T-EPA with hierarchical LU decomposition | |
CN104346488B (zh) | 电大复杂外形金属目标混合建模及电磁散射快速仿真方法 | |
CN112784459A (zh) | 一种基于压缩型有限元撕裂对接法的电磁仿真方法 | |
CN103279612A (zh) | 复杂目标雷达回波快速获取的多重网格预条件方法 | |
CN106294283A (zh) | 基于泰勒级数展开的时域积分方程快速算法 | |
CN105550451A (zh) | 基于辅助微分方程的一维左手材料Crank-Nicolson完全匹配层实现算法 | |
Wan et al. | Fast analysis of three‐dimensional electromagnetic problems using dual‐primal finite‐element tearing and interconnecting method combined with ℋ‐matrix technique | |
CN106156479A (zh) | 快速分析金属目标电磁散射特性的嵌套多层复点源方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |