本发明公开一种快速自适应地生成激励无关特征基函数的方法，首先对已分块的电磁目标做分块扩展，针对各个分块，利用快速自适应交叉近似算法自适应地选择一定数量的具有不同极化方式的入射平面波，生成激励矩阵并表示成三个矩阵相乘的形式；接着由扩展块上的自阻抗矩阵和激励矩阵的左矩阵计算出感应电流，并利用该感应电流矩阵生成初始特征基函数；然后，利用正交分解算法和截断奇异值分解算法去掉特征基函数中的冗余部分，得到初始特征基函数的截断奇异值分解形式；最后由该截断奇异值分解形式而得到最终的特征基函数。本发明针对不同电磁目标和分块情况可以自适应地选取不同极化方式不同入射角度的平面波数量，有效地提升生成特征基函数的效率。
1.一种快速自适应地生成激励无关特征基函数的方法，其特征在于，步骤如下：步骤1.1针对导体目标的表面用三角形面片进行离散，在每个相邻的三角形面片对上定义RWG基函数；根据导体目标表面边界条件建立用于散射计算的表面积分方程，利用所定义的RWG基函数和矩量法对表面积分方程进行离散；
步骤1.2对所有RWG基函数进行分块，并对每一个分块做块扩展；
步骤1.3针对每个扩展块，利用快速自适应交叉近似算法，设置入射平面波数量及极化方式，自适应地生成激励矩阵，并利用矩阵分解算法将其表示成三个矩阵相乘的形式；接着，依据该激励矩阵的分解形式生成初始特征基函数；然后，利用正交分解算法和截断奇异值分解算法对初始特征基函数中的矩阵形式进行分解和压缩，计算初始特征基函数的截断奇异值分解形式，并得到最终的特征基函数；
以第i个扩展块为例，步骤1.3中，利用快速自适应交叉近似算法，自适应地选择重要的元素，为生成激励矩阵提供条件；该算法的迭代过程中的第n次迭代，n＝1,2,3,…，步骤如下：步骤4.1在单位球面内均匀设置 个入射平面波，由对应的扩展块中选择 个RWG基函数，并选择 个入射平面波，其取值定义为：式中， 是第一次迭代时选取的入射平面波个数的初始值，i+指第i个扩展块， 是第i个扩展块包含的RWG基函数总数；其中，入射波在每个入射角下包含θ极化和 极化这两种极化方式，这样就可以将任意极化方向的入射波表示出来；
步骤4.2根据选择的RWG基函数，利用 个入射波的激励向量，组成激励矩阵步骤4.3利用自适应交叉近似算法对激励矩阵 压缩、分解，其分解形式为其中矩阵U(n)的大小为 矩阵V(n)的大小为 这里， 是
步骤4.4针对步骤4.3得到的结果，判断其是否满足收敛判据：
若满足，则停止迭代，不满足则继续迭代；其中，β∈[0.01,0.1]； 是第n-1次迭代得到的入射平面波数量； 是第n-1次迭代由自适应交叉近似算法得到的激励矩阵的有效秩；
步骤4.5对于满足收敛判据的情况，记录下迭代收敛对应迭代次数下步骤4.3中选择的个行、列的序号；其中，nt指满足收敛判据的迭代次数， 即为第nt次迭代由自适应交叉近似算法得到的激励矩阵的有效秩；
步骤1.3中自适应地生成激励矩阵，并利用矩阵分解算法将其表示成三个矩阵相乘的形式；接着，依据该激励矩阵的分解形式生成初始特征基函数，步骤如下：步骤5.1根据矩阵分解算法，生成激励矩阵并表示成三个矩阵相乘的形式式中，nt指快速自适应交叉近似算法的满足收敛判据的迭代次数； 是该扩展块第nt次迭代得到的激励矩阵的形式，其矩阵大小为 C是由 中被选中的 个列按序组成的，其矩阵大小是 RT是由 中被选中的 个行按序组成的，其矩阵大小是D是由 中被选中的 行、 列所共有的元素，按序组成的矩阵，D-1为其逆矩阵，其矩阵大小为 其中， 是第nt次迭代时选取的入射平面波个数， 是第nt次迭代由自适应交叉近似算法得到的激励矩阵的有效秩；
步骤5.2由式(4)分解形式的左边矩阵C，求解感应电流矩阵：
式中，Ii即是第i个扩展块中所有基函数上由入射波引起的感应电流矩阵， 是第i个扩展块的自阻抗矩阵的逆， 是 的子矩阵，它是通过将 中的第i个分块对应的扩展块扩展部分的基函数对应的行去掉而得到的；
步骤1.3中利用正交分解算法和截断奇异值分解算法对初始特征基函数中的矩阵形式进行分解和压缩，计算初始特征基函数的截断奇异值分解形式，并得到最终的特征基函数，步骤如下：步骤6.1利用正交分解算法，使得感应电流矩阵Ii和矩阵R有如下分解形式：Ii＝Qi，1Ri，1 (7)
式中，Qi，1、Qi，2是正交矩阵，Ri，1、Ri,2是上三角矩阵；
步骤6.3利用截断奇异值分解算法，将矩阵内 矩阵表示成
式中，UR、VR分别是 的左、右奇异矩阵，SR是对角阵且其每个元素都是的奇异值；
步骤6.4将式(10)代入式(9)，得到初始特征基函数的截断奇异值分解形式为式中，矩阵向量积Qi,1UR的每一列都是Ji的左奇异向量，且这些列即是最终的特征基函数。
2.根据权利要求1所述的一种快速自适应地生成激励无关特征基函数的方法，其特征在于，步骤1.1中，针对导体目标的表面用三角形面片进行离散，该三角形面片的边长设置为0.05λ～0.1λ，其中λ为计算散射特性时入射平面波的频率所对应的波长。
3.根据权利要求1所述的一种快速自适应地生成激励无关特征基函数的方法，其特征在于，步骤1.2中，对每一个分块做块扩展，在原有分块基础上，向外扩展0.1λ的大小，其中λ为计算散射特性时入射平面波的频率所对应的波长。
一种快速自适应地生成激励无关特征基函数的方法
技术领域：
[0001] 本发明涉及一种快速分析电大尺寸导体目标电磁散射的方法，尤其涉及一种快速自适应地生成激励无关特征基函数的方法。背景技术：
[0002] 电大目标的电磁散射问题一直受到国内外学者的广泛关注。矩量法(Method of Moments,MoM)将电磁积分方程转化成矩阵方程，是计算目标散射特性的有效途径。但是传统矩量法的迭代求解的复杂度为O(N2)，这里N是未知量的数目，如此高的复杂度限制着传统矩量法在计算电大目标的应用。
[0003] 激励无关(Excitation Independent,EI)的特征基函数法(Characteristic Basis Function Method,CBFM)用不同入射方向和不同极化方式的平面波照射每块所产生的响应作为初始特征基函数。由于初始特征基函数对应于各个方向和各种极化的平面波在该块的响应，所以初始特征基函数能够表征出该块在平面波照射下的感应电流的所有特征。然后通过奇异值分解(Singular Value Decomposition,SVD)去除这些初始特征基函数之间的相关性，可以减少初始特征基函数中的冗余。由于该CBFM产生的特征基函数适用于任意激励，对不同的激励不需要重新生成特征基函数和缩减矩阵，所以在求解多激励电磁问题时具有明显优势，例如单站雷达散射截面(Radar Cross Section，RCS)的求解。
[0004] 近年来，国内外学者基于激励无关的初始特征基函数，提出了一些提升其计算和奇异值分解的效率的方法。通常情况下，每种极化方式的入射平面波都是在单位球面内均匀地选取的，对于不同电磁目标选取的入射平面波的数量，是完全依靠经验去设置的。显然，仅仅凭借经验去选取合适的入射平面波数量，是不可靠的。为了解决这个问题，Raj Mittra的课题组于2009年给出了这样一个建议数值，每种极化方向的入射平面波可选取的最大数量为8π2(rλ+1)2，其中rλ是能够完全包含该扩展块的最小球面的半径。然而，这种方法仅仅考虑了电磁目标的分块及扩展块的大小，忽略了电磁目标的结构复杂度性的影响，以及截断奇异值分解的误差的影响。对于简单结构的电磁目标，这种方法会使其有很大的冗余量，严重影响特征基函数的生成效率；对于复杂结构的电磁目标，平面波数量易于选择不足够，导致生成的特征基函数不能完整地体现出其特性。发明内容：
[0005] 发明目的：为了解决生成激励无关特征基函数效率不足的问题，本发明提出了一种快速自适应地生成激励无关特征基函数的方法。该方法基于快速自适应交叉近似算法，在生成特征基函数时，能够根据不同的电磁目标和分块自适应地选择每种极化方式的入射平面波数目，避免了选择数量不足导致其计算数值不准确，以及选择数量过多导致计算效率降低的问题。
[0006] 为了达到上述目的，本发明的技术方案实现的基本步骤如下：
[0007] 第一步：针对导体目标的表面用三角形面片进行离散，在每个相邻的三角形面片对上定义RWG基函数；根据导体目标表面边界条件建立用于散射计算的表面积分方程，利用所定义的RWG基函数和矩量法对表面积分方程进行离散；
[0008] 第二步：对所有RWG基函数进行分块，并对每一个分块做块扩展；
[0009] 第三步：针对每个扩展块，利用快速自适应交叉近似算法，设置入射平面波数量及极化方式，自适应地生成激励矩阵，并利用矩阵分解算法将其表示成三个矩阵相乘的形式；接着，依据该激励矩阵的分解形式生成初始特征基函数；然后，利用正交分解算法和截断奇异值分解算法对初始特征基函数中的矩阵形式进行分解和压缩，计算初始特征基函数的截断奇异值分解形式，并得到最终的特征基函数。
[0010] 与现有技术相比，本发明的优势在于：对于不同电磁目标的不同分块的结构，能够自适应的给出合适的不同极化方式的入射平面波的数目，使其不再需要人为通过经验来设置；在通过截断奇异值分解来获得最终特征基函数的过程中，利用正交三角分解显著降低了传统奇异值分解的计算复杂度。总体而言，相对于现有生成特征基函数的方法，本发明能够自适应地完成特征基函数的生成，有效地提升了生成效率，且该方法具有普适性。附图说明：
[0011] 图1是本发明树形结构分块示意图。
[0012] 图2是本发明块扩展方法示意图。
[0013] 图3是本发明快速自适应交叉近似算法的结果示意图。
[0014] 图4是本发明初始特征基函数表示成三个矩阵的形式。
[0015] 图5是本发明初始特征基函数的截断奇异值分解形式。
[0016] 图6是本发明理想电导体球在不同分块方法下 面的结果对比。
[0017] 图7是本发明理想电导体球在不同分块方法下 面的结果对比。具体实施方案：
[0018] 下面结合附图对技术方案的实施作进一步的详细描述：
[0019] 第一步：针对导体目标的表面用三角形面片进行离散，在每个相邻的三角形面片对上定义RWG(Rao-Wilton-Glisson)基函数。通常情况下，对目标离散的三角形面片的边长，依据需要计算的散射所在的频率下的相应波长λ，设置为0.1λ。
[0020] 接着，根据导体目标表面边界条件建立用于散射计算的表面积分方程(Surface Integral Equation,SIE)
[0021]
[0022] 并利用所定义的基函数和矩量法对表面积分方程进行离散；其中，是场点位置矢量，在这里被认为是满足远场条件下的任意一点； 是源点位置矢量，在这里被认为是导体目标表面的离散点； 是入射平面波电场，j是虚数单位， 是自由空间格林函数， 是源点 处的感应电流；∫s·ds'是针对源表面的面积分；是针对场点的梯度算子； 是针对源场的散度算子。
[0023] 第二步：对所有基函数进行分块(如附图1)，并对每一块做适当的扩展；通常情况下，扩展的尺寸为Δ＝0.1λ(如附图2)。该步骤的主要目的在于降低由分块带来的电流奇异性误差；
[0024] 第三步：生成激励无关的特征基函数，以第i个扩展块为例，其基本步骤是：
[0025] (一)利用快速自适应交叉近似(Fast Adaptive Cross Approximation,FACA)算法，自适应地选择重要的元素，为生成激励矩阵提供条件。这是一个迭代的过程，主要包括以下5步：
[0026] ①在单位球面内均匀设置 个入射平面波，由对应的扩展块中选择 个RWG基函数，并选择 个入射波，其取值定义为：
[0027]
[0028]
[0029] 式中， 是第一次迭代时选取的入射平面波个数的初始值，i+指代第i个扩展块，是第i个扩展块包含的RWG基函数总数。其中，第i个扩展块是由第i个分块扩展得到的；入射波在每个入射角下包含θ极化和 极化这两种极化方式，这样就可以将任意极化方向的入射波表示出来。
[0030] ②根据选择的RWG基函数，利用SIE计算得到 个入射波的激励向量 的子阵并组成激励矩阵：
[0031]
[0032] 式中，Ei,i＝1,2,3,...,M是对应于不同入射波在 个RWG基函数上的激励向量，其中 矩阵 的大小为 矩阵 的大小为 相当于从中选择出 行 列。
[0033] ③利用自适应交叉近似算法(Adaptive Cross Approximation,ACA)对 压缩、分解，其分解形式为
[0034]
[0035] 其中矩阵U(n)的大小为 矩阵V(n)的大小为 且二者为 的ACA分解矩阵。这里， 是由ACA分解得到的矩阵 的有效秩，这一步快速而高效地选取了激励矩阵 当中的较为重要的部分。这里， 中元素较为重要的这 行 列的序号是需要被记录下的。
[0036] ④对上述步骤的结果判断，若收敛判据
[0037]
[0038] 是满足的，则停止迭代，不满足则继续迭代；其中， 是第n-1次迭代得到的入射平面波数量； 是第n-1次迭代由自适应交叉近似算法得到的激励矩阵的有效秩；β∈[0.01,0.1]。这里，当式(6)已满足后，其不等号左边形式的数值基本趋于稳定；若仅仅在设置β为更小的值，以减少其误差，其效果并不明显，且其误差的减少是需要更多的计算时间来获得的。
[0039] ⑤对于满足收敛判据的情况，假设第nt步迭代收敛，记录下该次迭代下的步骤③中选择的 个行、列的序号。
[0040] (二)根据矩阵分解算法(Matrix Decomposition Algorithm,MDA)，生成激励矩阵并表示成三个矩阵相乘的形式
[0041]
[0042] 式中， 是第nt代收敛得到的激励矩阵的形式，其矩阵大小为 C是由中被选中的 个列按序组成的，其矩阵大小是 RT是由 中被选中的 个
行按序组成的，其矩阵大小是 D是由 中被选中的 行、 列所共有的元
-1
素，按序组成的矩阵，D 为其逆矩阵，其矩阵大小为 其示意图如附图3，图中假设为18×12矩阵，并选中了其中3行3列。
[0043] (三)由式(7)分解形式的左边矩阵C，以第i个扩展块为例，求解感应电流矩阵：
[0044]
[0045] 式中，Ii即是第i个扩展块中所有基函数上由入射波引起的感应电流矩阵， 是第i个扩展块的自阻抗矩阵的逆， 是 的子矩阵，它是通过将 中的第i个扩展块扩展部分的基函数对应的行去掉而得到的。
[0046] (四)生成初始特征基函数(如附图4)：
[0047]
[0048] 通过该方法生成的初始特征基函数，而不是直接求解 目的在于减少矩阵求解的复杂度。
[0049] (五)利用正交分解算法，即QR分解，使得感应电流矩阵Ii和矩阵R有如下分解形式：
[0050] Ii＝Qi，1Ri，1 (10)
[0051] R＝Qi,2Ri,2 (11)
[0052] 式中，Qi，1(Qi，2)是正交矩阵，Ri，1(Ri,2)是上三角矩阵。那么，初始特征基函数可表示为：
[0053]
[0054] (六)利用截断SVD算法，将矩阵内 矩阵表示成
[0055]
[0056] 式中，UR(VR)是 的左(右)奇异矩阵，SR是对角阵且其每个元素都是的奇异值。该步骤通过截断的方式去除了那些奇异值为零或者很小的部分，以此来实现了对矩阵的压缩。
[0057] (七)将式(13)代入式(12)，得到初始特征基函数的截断SVD形式(如附图5)为[0058]
[0059] 式中，矩阵向量积Qi,1UR的每一列都是Ji的左奇异向量，且这些列即是最终的特征基函数。
[0060] 第四步：按第三步相同的原理，计算每一个分块的激励无关的特征基函数，最终得到所有的激励无关的特征基函数；接着，在基于特征基函数的算法中，利用这些已生成的特征基函数进行求解计算。
[0061] 下面以一具体实例对本发明方法作进一步说明：
[0062] 本发明以一个半径2.5米的理想电导体(Perfect Electric Conductor,PEC)球为例，按照0.1λ的间距离散，得到28890个RWG基函数。这里，λ是入射平面波的波长，该入射波是由 入射的300MHz的 方向极化的平面波。该算例中，生成特征基函数过程中的SVD的截断阈值ε＝10-3，FACA收敛判据中的参数β＝0.1。
[0063] 该算例中，对目标的分块采用了两种取值，分别是：1.0m和2.0m。表1给出了两种取值情况下的分块的数量、特征基函数的数量和 面的计算结果与解析解之间的误差。分块为2.0m时，传统方法设置入射平面波个数为1044个时，其生成特征基函数的时间为806s，使用本发明的方法，其生成特征基函数的时间为369s。本发明和解析解的结果对比，如附图6和附图7所示，其中附图6是 面的双站RCS值的对比，附图6是 面的双站RCS值的对比；从图中可以看出，两种分块方法的计算结果都吻合良好。该算例验证了本发明的准确性与适用性。在本次实验中，我们使用提出的自适应方法生成特征基函数，不需要人为的设置平面波的数目。
[0064] 表1
[0065]分块大小 分块数量 特征基函数数量 结果误差
1.0m 128 5767 0.08dBsm
2.0m 26 3053 0.01dBsm
