IJAP 国际期刊的天线和传播 1687 - 5877 1687 - 5869 Hindawi出版公司 10.1155 / 2016/6367205 6367205 研究文章 高效的宽带RCS计算使用奇异值分解增强改善超宽频特征基函数方法 http://orcid.org/0000 - 0002 - 7407 - 5537 Wen-yan 1 http://orcid.org/0000 - 0002 - 5709 - 0727 Zhong-gen 2 Ding-Bing 1 机电工程学院 淮南师范大学 淮南 安徽232001年 中国 hnnu.edu.cn 2 电气与信息工程学院 安徽科技大学 淮南 安徽232001年 中国 aust.edu.cn 2016年 14 12 2016年 2016年 26 08年 2016年 17 11 2016年 2016年 版权©2016 Wen-yan聂和Zhong-gen王。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

奇异值分解)增强改进超宽频特征基函数方法(IUCBFM)提出了有效分析宽带散射问题。传统IUCBFM,圣言只应用于减少之间的线性相关性特征基函数(cbf)由于过高的入射面波。然而,增加目标的大小根据分析需要大量的入射面波,它会变得非常耗时的解决这些数字的矩阵方程。摘要激励矩阵压缩通过圣言会为了减少许多矩阵方程的解决方案和cbf的数量与传统IUCBFM相比。此外,减少矩阵的尺寸和减少矩阵填充时间明显减少了。数值结果表明,该方法是准确和有效的。

中国国家自然科学基金 61401003 安徽省教育部门的自然科学基金 KJ2016A669 淮南师范大学自然科学基金 2015年xj09zd
1。介绍

系统响应在一个宽的频带要求在许多应用程序中,如现代雷达目标识别、微波遥感、微波成像。矩量法(MoM) [ 1)是最流行的一种数值方法用于雷达截面(RCS)的预测,但这地方一个沉重的负担对记忆和解决时间处理电大尺寸问题。许多不同的方法提出了规避这个问题快速解决方案,如快速多极子方法(FMM) [ 2),多层快速多极方法(MLFMM) ( 3, 4),自适应积分法(AIM) ( 5),自适应交叉近似(ACA)算法( 6),特征基函数方法(CBFM) [ 7, 8]。许多电磁应用程序需要的解决方案来自天线的辐射或散射问题在一个宽的频带,而不是单一频率点。然而,问题的解决方案使用传统CBFM或其他上述快速解决方法需要计算在每个频率点执行。CPU时间这将创建一个沉重的负担,尤其是在分析电大对象。现有计算宽带RCS的方法是基于插值妈妈矩阵( 9)或利用频率和频率导数数据( 10)以降低频率样本。然而,这些方法可能会过高的CPU时间,当身体的电尺寸很大。因此,在[ 11),一个超宽频特征基函数方法(UCBFM)提出了分析宽带电磁散射问题。超宽频特征基函数(UCBFs)建造在最高频率点感兴趣的范围可以重用在较低频率点没有重复的耗时的过程生成特征基函数(cbf)。然而,随着表示在 11],UCBFs通常是大型的RCS计算错误以较低的频率点。程序的原因是在较低的频率点使用离散化进行最高频率导致越来越多的条件,同时计算阻抗矩阵。此外,这些UCBFs一直在建造最高频率点,和他们的数量是高于必要水平低的频率点。这将导致减少矩阵的尺寸比要求。不必要的UCBFs数量的增加与减少的频率。结果,计算时间减少矩阵的计算和解决方案将远高于低频点。为了减轻这些问题,已经提出了一些改进的方法。在[ 12],UCBFs的建设已被充分考虑改进的子群之间的相互耦合影响cbf (SCBFs),获得二级水平,改善UCBFs (IUCBFs)含有更多的当前信息,提高计算精度较低频率点。在[ 13),提出了一种适应性IUCBFs施工方法和自适应IUCBFs取得了最高的频率点在每一次能带导致少数IUCBFs和显著减少在低频段的解决时间。在[ 14),去掉了不必要的UCBFs设置为基础的频率降低了。然而,在生成步骤的IUCBFs [ 12- - - - - - 14),入射波数的选择(PWs)通常是由经验和冗余浆通常选择增加cbf的解决时间。为了缓解这个问题,奇异值分解)适应加强cbf生成的效率。在生成cbf的新方案,减少矩阵方程的解决方案需要和cbf减少的数量。此外,减少矩阵的尺寸和减少矩阵填充时间明显减少了。

本文的其余部分组织如下。节 2,传统的IUCBFM简要描述。部分 3描述了圣言的使用有效地加快cbf代减少PWs的数量。节 4给出了一些数值结果,结论部分 5

2。改善超宽频特征基函数方法

IUCBFM [ 12首先将物体分成<我nline-formula> 块。然后,它建立了一个模型在最高频率点。Multiangle PWs将照射每一块。假设<我nline-formula> N θ 和<我nline-formula> N ϕ 代表的数量规模的方向<我nline-formula> θ 和<我nline-formula> ϕ ,分别<我nline-formula> N p w 年代 = 2 N θ N ϕ (两个偏振模式被认为是),指出<我nline-formula> E N p w 年代 。获得每个块的PCBFs,必须解决以下系统: (1) Z · J C B F = E N p w 年代 , 在哪里<我nline-formula> Z 是一个<我nline-formula> N × N 固有阻抗矩阵的块<我nline-formula> ,因为<我nline-formula> = 1、2 , , ;<我nline-formula> N 代表未知的数量扩展块<我nline-formula> ;<我nline-formula> E N p w 年代 是一个<我nline-formula> N × N p w 年代 激励矩阵;和<我nline-formula> J C B F PCBFs矩阵的维度<我nline-formula> N × N p w 年代 。块的PCBFs<我nline-formula> 可以直接得到解决( 1)。为了提高精度较低的频率点,构造SCBFs利用Foldy-Lax方程理论( 12, 15]。的SCBFs<我nline-formula> J 年代 1 和<我nline-formula> J 年代 2 块<我nline-formula> 计算了 (2) Z J 年代 1 = - - - - - - j = 1 j Z j J j P , Z J 年代 2 = - - - - - - j = 1 j Z j J j 年代 1 , 在哪里<我nline-formula> Z j 的阻抗矩阵块吗<我nline-formula> 并阻止<我nline-formula> j 。上面描述的过程后,<我nline-formula> 2 N θ N ϕ J P ,<我nline-formula> 2 N θ N ϕ J 年代 1 ,<我nline-formula> 2 N θ N ϕ J 年代 2 可以获得。通常,PWs的数量被用来生成cbf的数量将超过自由度与块相关联的,因此,需要使用奇异值分解过程去除这些cbf的冗余。只有相对的奇异值超过一定阈值,例如,<我nline-formula> 1.0 E - - - - - - 3 IUCBFs,保留。假设有<我nline-formula> K IUCBFs保留在每一块圣言,表面电流可以表示为一个线性组合的IUCBFs如下: (3) J = = 1 k = 1 K α k J C B F k , 在哪里<我nline-formula> J C B F k 代表了<我nline-formula> k th IUCBFs块<我nline-formula> 和<我nline-formula> α k 代表未知的重量系数。伽辽金方法用于传统的妈妈方程转化为线性方程系数矩阵<我nline-formula> α 。一个<我nline-formula> K × K 降低矩阵可以得到: (4) Z R · α = V R , 在哪里<我nline-formula> V R = J T · E ,<我nline-formula> T 代表换位,<我nline-formula> Z R 代表的阻抗矩阵维度<我nline-formula> K × K 。详细的计算表达式可以表示如下: (5) Z j R = J T · Z j · J , j , 在哪里<我nline-formula> Z j 代表之间的阻抗矩阵<我nline-formula> 和<我nline-formula> j 块。通常情况下,减少矩阵的维度<我nline-formula> Z R 小于生成通过传统的妈妈,然后呢<我nline-formula> α 可以直接得到解决( 4)。通过这种方式,可以获得任意频率点的表面电流。虽然IUCBFs可以提高计算精度较低频率点,应该注意的是,增加目标的大小根据分析将增加的数量规模和生成过程的PCBFs IUCBFM仍将需要昂贵的时间。需要加快PCBFs代减少的数量规模。

3所示。圣言会增强cbf代

为了获得足够的cbf对于一个给定的计算阈值和任意几何形状的块,冗余浆通常选择( 11- - - - - - 14]。然而,规模增加的冗余cbf的解决时间。为了缓解这个问题,计算过程适应加强cbf生成的效率。

在(做 12),每一块辐照multiangle PWs和一个激励矩阵<我nline-formula> E N p w 年代 获得与大小<我nline-formula> N × N p w 年代 。通常,PWs的数量被用来生成cbf的数量将超过自由度与块相关联的,所以激励矩阵将包含线性依赖关系,可以消除冗余进行压缩的事件。本文应用奇异值分解过程去除冗余的激励矩阵。这是通过表达这后者 (6) E N p w 年代 = U D V T , 在哪里<我nline-formula> U 是一个<我nline-formula> N × N 正交矩阵,<我nline-formula> V T 是一个<我nline-formula> N p w 年代 × N p w 年代 正交矩阵,<我nline-formula> D 是一个<我nline-formula> N × N p w 年代 对角矩阵。的上标<我nline-formula> T 表示转置操作。作者保留从左边的列矩阵奇异值<我nline-formula> U 超过一个阈值是谁的奇异值<我nline-formula> ε ,通常选择<我nline-formula> 1.0 E - - - - - - 3 。因此,一个新的激励矩阵<我nline-formula> E N e w 获得和事件的数量规模却降低了。为简单起见,假设<我nline-formula> l 浆为每个事件块奇异值分解后的尺寸<我nline-formula> E N e w 是<我nline-formula> N × l ,在那里<我nline-formula> l 总是小于<我nline-formula> N p w 年代 。后<我nline-formula> E N p w 年代 取而代之的是<我nline-formula> E N e w 在( 1),PCBFs<我nline-formula> J P 块<我nline-formula> 可以通过求解以下线性方程组: (7) Z · J C B F = E N e w

一旦解决,每一块的PCBFs SCBFs的生成(描述的是一样的 12]。矩阵方程的总数产生cbf所需的解决方案<我nline-formula> 3 · l ,这是小于<我nline-formula> 3 · N p w 年代 在IUCBFM<我nline-formula> l N p w 年代 。上面描述的过程后,<我nline-formula> 3 l 为每个块(包括cbf会存在<我nline-formula> l PCBFs和<我nline-formula> 2 l SCBFs)。再一次,计算过程是用来进一步减少这些cbf之间的线性相关性。为简单起见,假定所有的块包含相同的号码<我nline-formula> l N e w IUCBFs。的尺寸减少矩阵将减少<我nline-formula> l N e w × l N e w 相比,这是一个显著减少<我nline-formula> K × K 在IUCBFM<我nline-formula> l N e w K

4所示。数值结果

为了验证该方法的准确性和效率,给出了几个数值例子从压电陶瓷散射物体。所有的计算都是在个人电脑进行3.0 GHz Intel (R)奔腾(R) G2030 CPU和4 GB RAM(只使用一个核心)。第二层次的SCBFs计算。在传统IUCBFM阈值生成cbf的圣言<我nline-formula> 1.0 E - - - - - - 3 。在该方法(SVD-IUCBFM),圣言的阈值压缩规模和cbf将生成<我nline-formula> 1.0 E - - - - - - 3 和<我nline-formula> 1.0 E - - - - - - 2 ,分别。为了降低计算复杂度较低的频率点,使用自适应IUCBFs施工方法( 13]。

首先,散射问题的压电陶瓷球体半径为0.1米是0.2到2 GHz的频率范围。对象被分为2200个三角形补丁平均长度<我nline-formula> λ / 10 在2 GHz和未知的总数是5790。指的是( 12),每一块被multiangle浆<我nline-formula> 0 ° θ < 180年 ° 和<我nline-formula> 0 ° φ < 360年 ° 与<我nline-formula> N θ = 8 和<我nline-formula> N ϕ = 8 ,导致384规模。在SVD-IUCBFM 8 PWs的方向<我nline-formula> θ 和<我nline-formula> ϕ 设置。圣言后,只有37个新激励向量(平均值)相对于最大奇异值保留在每个块,结果在一个新的激励矩阵<我nline-formula> E N e w 。与IUCBFM相比,激励矩阵的维度都显著降低。通过激动人心的每一块新的激励矩阵和求解线性方程组,37岁<我nline-formula> J P 37岁的<我nline-formula> J 年代 1 ,37<我nline-formula> J 年代 2 可以在每个块。圣言会再次用于进一步减少这些cbf之间的线性相关性。IUCBFs保留在每个块的数量和两种方法的计算时间在每个子带表 1在细节。它可以很容易地观察到的数量IUCBFs获得使用获得的SVD-IUCBFM显然小于使用IUCBFM和奇异值分解的时间减少了,因为cbf减少矩阵的维度。收发分置的RCS的<我nline-formula> θ θ 极化计算IUCBFM,妈妈(FEKO), SVD-IUCBFM 1100 MHz如图 1。可以看出SVD-IUCBFM显示与妈妈一个很好的协议。宽带RCS(10在每个部分波段频率采样点)获得通过使用IUCBFM和SVD-IUCBFM在0.2到2 GHz的频率范围,如图 2。结果计算通过使用SVD-IUCBFM IUCBFM获得的结果吻合较好。表 1浆和IUCBFs列表的数量和两种方法的计算时间。可以看出SVD-IUCBFM降低了浆通过圣言与IUCBFM相比。此外,生成cbf的效率和减少矩阵构造都是改善。

浆和IUCBFs数量和计算时间的球体。

IUCBFM SVD-IUCBFM
频率子带(GHz) 0.2 - -0.65 0.65 - -1.1 1.1 - -2.0 0.2 - -0.65 0.65 - -1.1 1.1 - -2.0
数量的规模 128年 128年 128年 23 33 54
计算时间(年代) 9.55 9.03 8.25 5.61 6.09 8.29
IUCBFs数量 344年 551年 809年 316年 517年 706年
IUCBFs施工时间(年代) 216.54 213.27 207.21 98.94 111.94 152.79
降低矩阵填充(s) 108.27 279.17 564.01 87.56 207.21 418.26
总时间(s) 2135.25 4159.47 7016.52 1919.54 3336.54 5365.49

收发分置的RCS的球体在1100 MHz。

宽带RCS的球体。

第二,压电陶瓷缸0.1到2 GHz的频率范围。圆柱体的半径和高度是0.1米和0.4米,分别。对象被分为5674个三角形补丁和未知数的数目是12363。气缸分为6块。宽带RCS(10频率采样点在每个部分波段)计算传统IUCBFM和SVD-IUCBFM比较图 3。表 2总结了数字的PWs IUCBFs和圆柱体的计算时间。可以看出SVD-IUCBFM优于传统IUCBFM, cbf代和RCS计算时间。特别是PWs的数量和降低矩阵填充时间显著减少,收益约58%和25%,分别。

浆和IUCBFs和油缸的计算时间。

IUCBFM SVD-IUCBFM
频率子带(GHz) 0.1 - -0.575 0.575 - -1.05 1.05 - -2.0 0.1 - -0.575 0.575 - -1.05 1.05 - -2.0
数量的规模 128年 128年 128年 31日 48 81年
计算时间(年代) 82.70 79.27 79.29 65.66 77.32 81.86
IUCBFs数量 402年 649年 1044年 334年 549年 951年
IUCBFs施工时间(年代) 1318.23 1315.72 1311.26 875.08 1012.41 1252.8
降低矩阵填充(s) 5684.77 9976.09 15501.19 4293.39 7150.81 12070.43
RCS计算时间(年代) 63203.93 102016.86 161399.16 48847.90 77564.81 125723.17

宽带RCS的缸。

最后,派克立方体边长0.1米的0.1到5 GHz的频率范围。三角片离散化是在13607年5 GHz导致许多未知数。几何分为8块。宽带RCS计算IUCBFM和SVD-IUCBFM图所示 4。从图可以看出一个好的协议。浆和IUCBFs的数目和计算时间如表所示 3。可以看出SVD-IUCBFM需要一个较小的规模利用奇异值分解过程与传统IUCBFM相比。此外,减少矩阵维度和减少矩阵填充时间都减少了。

浆和IUCBFs和油缸的计算时间。

IUCBFM SVD-IUCBFM
频率子带(GHz) 0.1 - -0.7125 0.7125 - -1.325 1.325 - -2.55 2.55 - -5.0 0.1 - -0.7125 0.7125 - -1.325 1.325 - -2.55 2.55 - -5.0
数量的规模 128年 128年 128年 128年 20. 24 41 73年
计算时间(年代) 50.37 50.55 51.09 51.69 30.91 32.73 43.21 73.27
IUCBFs数量 365年 476年 728年 1341年 247年 355年 617年 1191年
IUCBFs施工时间(年代) 1421.12 1430.03 1428.43 1436.81 555.94 612.81 820.01 1072.36
降低矩阵填充(s) 3416.11 5895.31 9772.54 19744.71 2324.82 5747.99 8591.35 14608.87
RCS计算时间(年代) 41810.12 66620.11 105411.94 205182.12 29994.712 48222.06 92996.69 153532.42

宽带RCS的多维数据集。

5。结论

本文提出了一种新方法有效地计算出压电陶瓷宽带RCS对象。建议的方法,所需的PWs的数量已经明显减少进一步压缩事件PWs使用奇异值分解,导致更少的矩阵方程的解决方案,减少时间的cbf的一代。此外,减少矩阵的尺寸和RCS计算时间都减少了与传统IUCBFM相比。的结果表明,该SVD-IUCBFM能够更有效地计算宽带RCS与传统IUCBFM不影响精度。

相互竞争的利益

作者宣称没有利益冲突。

确认

这项工作是由中国国家自然科学基金支持下批准号61401003,安徽省自然科学基金的教育部门批准号KJ2016A669,淮南师范大学自然科学基金批准号下2015 xj09zd。

哈林顿 r F。 字段由矩量法计算 1992年 纽约,纽约,美国 IEEE出版社 Coifman R。 Rokhlin V。 Wandzura 年代。 波动方程的快速多极方法:一个行人的处方 IEEE天线和传播杂志 1993年 35 3 7 12 10.1109/74.250128 2 - s2.0 - 0027609499 首歌 J。 c c。 咀嚼 w . C。 多层快速多极算法电磁散射的大型复杂的对象 IEEE天线和传播 1997年 45 10 1488年 1493年 10.1109/8.633855 2 - s2.0 - 0031258501 M。 r S。 x Q。 增强MLFMM从压电陶瓷散射分析对象与精细结构 应用计算电磁学社会杂志上 2011年 26 5 418年 428年 2 - s2.0 - 80052467530 Bleszynski E。 Bleszynski M。 Jaroszewicz T。 目的:自适应积分方法求解大规模的电磁散射和辐射问题 无线电科学 1996年 31日 5 1225年 1251年 10.1029/96 rs02504 2 - s2.0 - 0030242018 K。 Vouvakis m . N。 肯尼迪。 自适应交叉近似算法加速矩量法计算电磁兼容问题 IEEE电磁兼容性 2005年 47 4 763年 773年 10.1109 / temc.2005.857898 2 - s2.0 - 31744437470 普拉卡什 诉诉。 Mittra R。 特征基函数方法:有效解的新技术时刻矩阵方程的方法 微波和光学技术的信件 2003年 36 2 95年 One hundred. 10.1002 / mop.10685 2 - s2.0 - 0037454813 Lucente E。 Monorchio 一个。 Mittra R。 iteration-free妈妈方法基于激发独立特征基函数求解大型多尺度的电磁散射问题 IEEE天线和传播 2008年 56 4 999年 1007年 10.1109 / TAP.2008.919166 2 - s2.0 - 42549130424 纽曼 e . H。 代的矩量法的宽带数据通过插值阻抗矩阵(EM问题) IEEE天线和传播 1988年 36 12 1820年 1824年 10.1109/8.14404 2 - s2.0 - 0024171880 伯克 g . J。 米勒 e·K。 Chakrabarti 年代。 Demarest K。 使用基于模型的参数估计效率的增加计算电磁传输功能 IEEE磁学 1989年 25 4 2807年 2809年 10.1109/20.34291 2 - s2.0 - 0024704069 德格雷戈里奥 M。 Tiberi G。 Monorchio 一个。 Mittra R。 宽带散射问题的解决方案使用特征基函数方法 专业微波、天线和传播 2012年 6 1 60 66年 10.1049 / iet-map.2011.0309 2 - s2.0 - 84857621813 W.-Y。 Z.-G。 宽带散射问题解决方案通过使用改进的超宽波段特征基函数方法 电磁学研究快报的进展 2016年 58 37 43 10.2528 / pierl15080801 2 - s2.0 - 84951263852 w . Y。 z G。 使用自适应宽带散射的分析对象改善超宽波段特征基函数 电磁学研究快报的进展 2016年 60 45 51 Koc s . N。 Koksal 一个。 宽带分析使用CBFM平面可伸缩的天线和压电陶瓷的身体 土耳其的电气工程和计算机科学杂志》上 2016年 24 1652年 1662年 10.3906 /麋鹿- 1401 - 48 z G。 太阳 y F。 g . H。 从完美的导电目标电磁散射分析使用改进的特征基函数法和快速偶极子法 杂志的电磁波和应用程序 2014年 28 7 893年 902年 10.1080 / 09205071.2014.895425 2 - s2.0 - 84896284724