靶标雷达特性快速预测与分析.pdf

第4期2023年8月Vol.21 No.4August 2023雷达科学与技术Radar Science and TechnologyDOI:10.3969/j.issn.16722337.2023.04.003靶标雷达特性快速预测与分析黄巍浩1，何思远1，张云华1，李婉聪1，杨泽望2，刘建2（1.武汉大学电子信息学院，湖北武汉 430072；2.中国人民解放军63921部队，北京 100094）摘要：含多次散射复杂靶标雷达特性的快速预估和分析是靶标设计优化的重要基础，同时也是电磁散射建模领域的重要研究课题之一。多次散射结构及其复杂散射机理导致靶标雷达特性仿真效率低下、特性分析困难。本文提出了一种CPU多核并行技术、GPU硬件加速技术和KDTree遍历技术相结合的靶标高频电磁散射加速方法，建立了带腔舰船、角反射器阵列与舰船等含多次散射复杂靶标的高频电磁散射模型，可以满足大型靶标进行快速评估的需求。本文进一步提出了多次散射距离像位置快速预测方法，通过多次散射射线分集结合射线路径相位理论预估，快速推导了多次散射等效视在中心的位置，揭示了含腔靶标的多路径散射机理与作用过程；并通过舰船目标的仿真与分析，建立了大型靶标多次散射形成的强散射贡献与其具体几何结构之间的映射关系。关键词：靶标雷达特性；高频算法；多次散射；并行计算；快速预测中图分类号：TN957.52文献标志码：A文章编号：16722337（2023）04037509Fast Prediction and Analysis of Radar Characteristics of TargetHUANG Weihao1，HE Siyuan1，ZHANG Yunhua1，LI Wancong1，YANG Zewang2，LIU Jian2（1.School of Electronic Information,Wuhan University,Wuhan 430072,China；2.Unit 63921 of PLA,Beijing 100094,China）Abstract:Fast prediction and analysis of radar characteristics of complex targets with multiple scattering is an important basis for target design optimization and is also one of the important research topics in the field of electromagneticscattering modeling.Multiple scattering structure and complex scattering mechanism lead to low simulation efficiencyand difficulty in characteristic analysis of radar target.In this paper,a high frequency electromagnetic scattering acceleration method for targets is proposed,which combines CPU multicore parallel technology,GPU hardware accelerationtechnology and KDTree traversal technology.The high frequency electromagnetic scattering model of complex targetsincluding ships with cavity,angular reflector arrays and ships is established,which can satisfy the needs of rapid evaluation of large targets.Further,a fast prediction method for multiplescattering range image position is proposed.By combining multiplescattering ray diversity with ray path phase prediction theory,the multiplescattering equivalent apparent position of the center is deduced rapidly,and the multipath scattering mechanism and action process of the targetcontaining cavity are revealed.Through the simulation and analysis of ship target,the mapping relationship between thestrong scattering contribution formed by multiple scattering of large target and its specific geometric structure is established.Key words:target radar characteristics;highfrequency algorithm;multiplescattering;parallel computing;fastprediction0引言靶标设计和优化对提升武器系统性能有着重要的支撑作用，靶标雷达特性预估和分析则是靶标设计优化的重要基础。多次散射结构广泛存在于各类靶标外形结构和靶场应用之中，如靶场角反射器及其阵列的使用、飞机进气道结构、舰船复杂上层建筑形成的多次散射结构等12。如何准确、高效地实现对含多次散射复杂靶标雷达特性的快速预估和分析，仍然是电磁散射建模领域的重要研究课题34。高频电磁散射建模方法具有计算效率高、物理机理清晰、工程应用性强的优点，非常适合于大型靶标宽带雷达特性的快速计算和分析5。收稿日期：20230204；修回日期：20230329基金项目：国家自然科学基金（No.61571333）雷达科学与技术第 21 卷 第 4 期然而，在靶场尤其是大型靶标设计的实际工程应用中，为了获取合理的设计结果，需要大量重复计算和优化，现有高频计算方法依然难以满足靶标设计高效性的需求，靶标雷达特性快速预估和分析面临诸多困难和挑战。这些挑战和困难主要来自以下几个方面。首先，对电大尺寸靶标的高分辨一维距离像（HRRP），二维雷达图像（SAR）等目标宽带雷达特性的需求加剧了电磁计算的困难6。靶标宽带电磁散射特性仿真必须克服雷达宽带信号多个频点和目标多姿态采样所带来计算量大的挑战7。其次，为了满足靶标特性计算精度的需求，含多次散射结构（靶场角反射器及其阵列、腔体等）在电磁建模过程中，需要计入三次及三次以上的电磁高阶耦合作用，开展三次以上的高频射线追踪。这类多次的射线追踪过程极其复杂，计算非常耗时，导致靶标雷达特性预估面临计算效率低下的问题8。第三，含多次散射结构靶标涉及复杂多次电磁散射机理，对其宽带雷达特性形成机制的分析同样面临巨大挑战。多次散射形成的等效视在中心（距离像位置或者散射中心位置）往往偏离目标区域9，对多次散射机理和等效视在中心位置的分析，目前还缺乏有效的分析手段。为了提升电大尺寸复杂目标高频方法计算效率，研究者们采用加速算法和并行计算技术改进高频仿真过程中最耗时的射线追踪模块。Suk等人提出利用多分辨率网格算法减少射线管的数量从而减少计算量10；Glassner等人提出空间划分算法，利用面元之间的空间关系，通过减少射线追踪过程中求交判断的次数来减少计算量11；随后KDTree和八叉树等算法在此基础上不断发展12；国内浙江大学团队利用GPU平台同时计算众多射线管，通过并行计算来提高计算效率1314。除了上述加速算法和GPU并行技术，CPU多核技术近年来迅猛发展，为并行电磁计算的发展创造了发展空间和条件1516。本文针对含多次散射复杂靶标雷达特性的快速计算难题，提出了一种 CPU 多核并行技术、GPU 硬件加速技术和 KDTree遍历技术相结合的靶标电磁散射加速计算方法。在此基础上，建立了带腔舰船、角反射器阵列与舰船等含多次散射复杂靶标的电磁散射（高频）模型，可以满足大型靶标进行快速评估的需求。针对多次散射特性分析困难，本文进一步提出了多次散射距离像位置快速预测方法。通过多次散射射线分集结合射线路径相位理论预估，预测了多次散射等效视在中心的位置。将多次散射等效视在中心理论预估位置与高频仿真雷达特性进行比较，揭示了腔体等含多次散射结构复杂的多路径散射机理与作用过程。通过具体数值算例，建立了靶标多次散射雷达特性与其多次散射几何结构之间的映射关系。本文结构安排如下：首先针对靶标雷达特性仿真效率低下的问题，提出了基于三种计算机加速技术的靶标电磁散射高频建模加速方法。然后，针对多次散射特性分析困难，基于电磁场射线理论相位预估，发展了一种多次散射距离像位置快速预测方法。最后，通过数值仿真和分析，表明了本文方法的有效性。1靶标电磁散射加速计算方法图 1 展示了传统的靶标雷达特性高频建模流程。,P1CAD4PM*-)?D44*3EC30FLLDFM(LD,4(SAR)4C.(HRRP)图1雷达目标特性获取流程图针对靶标建模对象，首先需要建立对应的高精度CAD模型，并将目标表面剖分为若干面元，后续根据目标的点面信息进行电磁计算。在对模型预处理之后，将目标散射贡献按照高频散射机理进行分解，可将复杂目标电磁散射贡献分解为边376黄巍浩：靶标雷达特性快速预测与分析2023 年第 4 期缘绕射场、面元散射及多次射线场等贡献之和。在给定计算姿态下，通过射线追踪确定电磁场的传播路径，并结合物理光学法计算该姿态下的总射线场。进一步可结合时频变换获取宽带雷达特性数据，包括高分辨一维距离像数据以及合成孔径雷达二维图像数据。上述建模流程中，复杂多次射线追踪和大量的目标宽姿态宽频带采样是电磁计算耗时的两个主要原因。在射线追踪过程中，需要所有射线和所有面元两两进行照亮判断，此时的时间复杂度为O(N2)，N为面元数量。针对电大尺寸复杂靶标高频区散射，N往往数以万计，此时要获取雷达宽带特性数据的计算代价是巨大的。尤其是当目标存在多次散射结构时，每一次反射点记录需要重复射线与面元两两判断的过程，成倍地增加了计算时间。考虑到射线追踪的过程中，各个射线之间的追踪互不相关，非常适合大规模并行计算。本文通过 GPU 平台来对多次射线追踪进行加速。GPU上有大量的计算单元，在每个计算单元上创建线程进行射线的追踪运算，实现所有射线并行追踪的功能，此时射线追踪的时间复杂度降为O(N)。同时，利用目标各面元的空间位置关系，进一步引入KDTree数据结构，仅对射线与发射方向附近的面元进行求交测试，此时射线追踪的时间复杂度进一步下降到O(log2N)与O(N)之间。传统仿真串行地依次计算每一个姿态和频点下的回波数据，使得目标宽频带和宽姿态采样的计算非常耗时。考虑到各个姿态下的计算流程相同，且计算过程中的中间数据互不相关，本文采用功能更为强大的CPU来创建指定数量的多线程实现并行计算。由CPU创建的每一个子线程在射线追踪环节都将访问GPU，之后继续在CPU端完成后续的射线场计算。而GPU作为计算过程中的共享资源，在被其中一个子线程访问时，需要阻止其他线程的访问。本文利用OpenMP编程模型来协调线程和管理共享资源，使得雷达目标多姿态并行计算顺利进行。如图2所示，通过OpenMP中的critical指导语句设立临界区，临界区内的内容同一时刻只能有一个线程访问，将访问GPU的代码块放置在临界区内，可以保证同一时刻只有一个线程在访问GPU，从而实现了共享资源的保护。41234GPU3EC3C13C23C33C4LD,RCS(*)(=)30303030图2多姿态并行计算模型在上述过程中，子线程数量并非越多越好，需要确定合理的子线程数量。由于不同线程之间对共享资源GPU存在竞争关系，当创建一定数量的线程之后，GPU将达到满负载工作状态，整体效率不会再因为线程数量增多而增加。本文根据线程在GPU和CPU上运行的时间长短来创建合理数量的线程。如果把GPU端进行射线追踪的时间记为A，把CPU端后续射线场计算部分的时间记为B，那么最合理的线程数N可以由式（1）得到：N=1+BA（1）图3展示了两种不同情况的用时占比。第一种情况下，GPU上用时与CPU上用时相等，此时创建的最佳线程数为2。第二种情况下CPU上用时情况为GPU上的2倍，最佳线程数为3。两种情况下，GPU上线程运行状况如图4所示。可以看出，此时GPU已经达到满负载工作，无法通过创建更多的线程来提高系统效率。GPU*CPU*3/13/23/33/13/212图3GPU用时不同占比的两种情况3/1GPUD=1GPUD=23/13/13/13/23/23/3图4两种情况下GPU运行情况2多次散射距离像位置预测方法本文提出了多次散射距离像位置快速预测方377雷达科学与技术第 21 卷 第 4 期法。通过多次散射射线分集结合射线路径相位理论预估，预测了多次散射等效视在中心的位置。针对复杂目标多次散射，本文通过射线追踪记录电磁波与目标相互作用的反射点，获取电磁波传播过程中的射线路径集合，计算每一条射线的光程。根据射线光程可以快速推导多次散射等效视在中心的位置，即距离像位置或者散射中心位置。一维距离像峰值的射线理论预测值X可以按式（2）计算：X=L2-R（2）式中L为光程，R为目标参考系原点到雷达的距离。为了说明光程的计算过程及其物理含义，本文选取多次散射典型结构，以矩形腔体为例，进行具体的分析。当雷达入射方向与腔体横截面平行时，射线在一个平面内传播，此时可以清晰简洁地分析横截面上射线的传播路径，并给出峰值预测结果。以腔体口径面为界，将射线分为腔体内部射线以及腔体外部射线，然后分别计算光程。图5给出了腔体内部一条射线在口径面某一点入射，经四次反射后返回雷达的传播路径。DL3L1L4L3L4L2L5L1图5腔体内部射线追踪情况按反射情况可以将腔体内部的射线分为5段，分别标记为L1,L2,L3,L4,L5。将L1,L3,L4段以腔内壁为轴翻转，根据几何反射原理，翻转后的L1,L2,L3在同一条直线上，L4,L5在同一条直线上，由此可以得到该射线长度为L=L1+L2+L3+L4+L5=(L1+L2+L3)+(L4+L5)=2D/cos（3）由式（3）可以看出，射线光程在腔体内的长度仅由入射角度和腔体深度D决定，与入射位置无关。所以对于在口径面其他点入射的射线，光程均相等。而当反射次数增加时，同样可以通过上述方式计算光程。图6用于说明如何计算腔体内外总的射线长度。为方便分析，如图建立二维坐标系，并将入射射线与口径面的交点和出射射线与口径面的交点记为y1,y2，与雷达的距离记为a,b，则腔体外射线长度可记为a+b。根据几何关系，可以得到入射方向和反射方向平行，此时腔体外光程可以为y1y2中点到雷达距离的两倍，所以需要先求得y1y2中点。根据几何关系，有如下关系式：y1+m+H=Dtany2+n=Dtan（4）式中m,n为腔体底部被接触点分割的两部分，且二者之和为腔体底部长度H，可以求得(y1+y2)/2=Dtan-H（5）式（5）说明入射点与出射点的中点位置坐标仅受入射角度和腔体尺寸影响，从不同点入射的射线在腔体外的长度也相等。所以，针对矩形腔体，在同一入射姿态下形成的所有4次反射的射线路径光程相同，一维距离像位置相同。LDDmHnxyy1y2ba图6自由空间中射线路径图图 7 是预测腔体目标一维距离像的示意图。A,B分别为射线入射时与口径面的交点以及反射378黄巍浩：靶标雷达特性快速预测与分析2023 年第 4 期时与口径面的交点在雷达径向方向的投影。O点为目标坐标系原点在雷达径向方向上的投影，记坐标原点与雷达之间的距离为R，OA，OB分别为O点和A，B点之间的距离，有R=a+OA=b+OB。一维距离像的峰值坐标可以由式（6）求出，其中OA+OB的值可以通过上文求出的y1y2中点和雷达间的距离求出。X=(a+b+L-2R)/2=(L-OA-OB)/2=(2D/cos-OA-OB)/2（6）LDOABxyab图7自由空间中腔体目标成像示意图3靶标雷达特性仿真与分析三面角和腔体是靶标设计中常见的具有多次散射的典型结构，本节首先通过本文加速方法，获取三面角和矩形腔体两种结构的雷达特性。并将矩形腔体距离像射线理论预估位置与高频仿真距离像特性进行比较，两者吻合一致。进一步，针对带腔舰船、角反射器阵列与舰船复合目标开展了仿真与特性分析，仿真结果表明了本文方法的有效性。本文仿真环境如下：CPU型号为 Intel Corei79700，内存大小为16G，内含8核。GPU型号为NVIDIA GeForce GTX 1660Ti，显存大小为8 192 MB。3.1三面角三面角的几何模型如图 8 所示，面元数为34 292，点数量为17 361。入射射线与z轴夹角记为俯仰角，入射射线在XOY平面上的投影与x轴夹角记为方位角。三面角雷达参数设置如下：角范围为090，角为45，频率为3 GHz，极化方向为垂直极化。xyz图8三面角模型图9是本文计算方法得到的RCS曲线和商用软件FEKO计算得到的RCS曲线对比图，商用软件FEKO 所 用 计 算 方 法 为 快 速 多 极 子 算 法（MLFMM），两种计算方法结果基本吻合。0102030405060708090/()RCS/dBsm-5051015202530MLFMM图9三面角RCS结果对比分别在无加速、利用GPU和KDTree数据结构加速射线追踪过程以及CPU和GPU混合编程三种情况下记录所需用时，加速情况如表1所示。此时线程调用 GPU 的时间与在 CPU 端运行的时间相近，所以通过OpenMP的协调，与不调用OpenMP相比提高了约两倍的计算效率。表1不同方法计算三面角RCS曲线耗时对比加速情况无加速GPU+KDTreeCPU+GPU+KDTree总时间/s4 282.8834.3018.113加速比1124.86236.453.2腔体目标矩形腔体算例尺寸如图10所示。379雷达科学与技术第 21 卷 第 4 期400 mm300 mm300 mm200 mmyzx600 mm800 mm(a)(b)(c)(d)图10矩形腔体模型尺寸腔体目标雷达参数设置如下：角为90，角范围为090，频率为10 GHz，极化方向为垂直极化，商用软件FEKO所用方法为高频方法，计算结果相比如图11所示。0102030405060708090-40-30-20-1001020304FEKO/()RCS/dBsm图11腔体目标RCS结果对比选取RCS曲线中的极大值点，计算当=37时的一维距离像。不模糊距离为3 m，成像分辨率为0.01 m。图12给出了腔体目标在该入射角度下一维距离像的数值结果。仿真结果表明，腔体目标在该入射角度下形成了一个强散射中心，距离像峰值坐标为0.43 m。用本文射线路径相位理论预估的方法计算得到的结果为0.433 1 m，符合仿真结果，说明通过目标的几何结构预测目标的一维距离像峰值坐标方法的可靠性。-2-10120510E#C.N图12腔体目标一维距离像（=90,=37）选取RCS曲线中的其他极值点，分别计算角度为 51.3，60.2，66.1和 73情况下的一维距离像，不模糊距离为 5 m，成像分辨率为 0.01 m。结果如图13所示。-2-1012051051.360.266.173E#C.N图13不同入射角度下的腔体目标一维距离像不同姿态下的射线长度不同，所以得到的距离像峰值位置也不同。由图可见，随着入射角增大，即入射射线与腔体口径法向的夹角增大，射线经由腔体反射形成的强散射中心逐渐远离参考原点。同时，口径面在雷达方向的投影面积减小导致进入腔体的射线减少，一维距离像的峰值强度也逐渐减小。在实际工程中，射线不一定在一个平面内传播，多次反射的结构也可能是不规则的，此时同样可以根据本文射线追踪方法得到射线光程，进而预测目标的一维距离像。选取俯仰角=82，方位角=37的姿态进行射线追踪和分析，目标姿态如图14所示。图14腔体目标姿态图（=82,=37）通过记录所有射线路径，可以计算得到所有射线长度均为1.824 0 m。坐标原点到雷达接收点之间的距离为0.471 7 m，根据射线路径相位理论预估方法可计算目标一维距离像中峰值的理论坐标X=0.440 3 m。该姿态下仿真得到的一维距离像结果如图15所示，图中峰值坐标为0.44 m，与通过射线追踪预测的坐标一致，说明了通过射线长380黄巍浩：靶标雷达特性快速预测与分析2023 年第 4 期度预测目标一维距离像峰值的普适性。-1.5-1.0-0.500.51.01.50510E#C.N图15腔体目标一维距离像（=82,=37）3.3带腔舰船目标图16所示为某一驱逐舰船模型，舰船长166.7 m、宽17.2 m、高33 m，船舱部件之间存在多种多次散射结构。在舰船的船头、中部和船尾存在明显的腔体部件，腔体部件长8 m、宽5 m、高3 m。后文分别以前腔、中腔和后腔指代。舰船除3个腔体部件以外关于x轴对称。XYZ0.0035.0070.00m17.5052.50图16带腔舰船几何模型雷达参数设置如下：角为 90，角范围为0360，频率为10 GHz，极化方向为垂直极化，计算结果如图17所示。角为0时探测方向正对船头，散射贡献主要来自船体前部的平板结构。在90和270处RCS达到峰值，此时主要散射贡献来源于船侧的平板结构。在俯仰角为90的情况下，关于x轴对称目标的全姿态RCS曲线也具有对称性，关于方位角为0和180对称。后腔是唯一不关于x轴对称的部件，在方位角为20110的范围内后腔散射贡献主要来自于腔体内部的多次散射作用，在方位角为 110、200和 290时，后腔散射贡献主要来自于腔体外壁的一次作用，这些方位角对应的RCS值略高于对称姿态下的RCS值。对含腔目标而言，如果射线追踪的次数过少，则无法计入部分角度下腔体的散射贡献，会影响对应角度下RCS的计算精度。050100150200250300350020406080100/()RCS/dBsm图17带腔舰船全姿态RCS结果图图18为带腔舰船目标在俯仰角为65，方位角为37姿态下的一维距离像，简易舰船为不包含3个腔体部件的舰船模型。对带腔舰船目标的几个较强的峰值进行分析，形成峰值的等效散射源与目标几何结构的对应关系如图19所示，射线理论预估和仿真结果对比如表2所示。其中1，3，4号峰值点分别由3个腔体部件形成，前腔和中腔开口方向与船头一致，然而前腔截面偏离XOY面角度约3，导致其距离像幅度较低。2号峰值点为船体中部船舱各平板之间的三次作用形成。-100-80-60-40-2002040608010051015202530788-188-E#C.N图18带腔舰船与简易舰船一维距离像对比图1234图19一维距离像峰值与目标几何结构对应关系示意图表2舰船目标一维距离像位置验证序号1234射线理论预估值/m-30.348 330.346 542.288 669.154 2距离像位置仿真值/m-30.429 730.482 642.927 569.174 7腔体结构不仅在某些角度下能形成强散射中心，腔体内部的多次散射还能影响其等效视在中381雷达科学与技术第 21 卷 第 4 期心与部件几何中心投影之间的位置关系。如图20所示，前腔与中腔的开口方向一致，在该入射姿态下具有相同的光程，视在中心与几何中心投影均相距5.2 m。后腔在该入射姿态下的散射贡献来源于腔体内壁与底部的二次作用，视在中心与几何中心投影距离为3.4 m。可以看出，多次散射使得目标的等效散射源偏离了原有目标空间范围，偏离程度随着射线光程的增加而增大。70 m3.4 m5.2 m5.2 m-51 m64646KLN6KLN60607774LD6KLN60777464图20腔体等效视在中心与腔体几何结构关系示意图3.4角反阵列与舰船复合目标本文以图21所示的角反阵列与舰船复合模型为例，快速评估复合模型的散射特性及角反阵列回波干扰信号的有效性。舰船长 119 m、宽 10.8 m、高12.88 m，角反射器阵列分别包含边长为0.5 m、1 m、2 m的角反射器各4个，这些角反射器开口朝向各不相同，复合模型的网格数量为278 000。xyz图21角反射器阵列与舰船目标复合模型当雷达参数为俯仰角60，中心方位角45，中心频率 10 GHz，极化方式为垂直极化，分辨率为1 m时，获取图 23所示的复合模型 SAR图像数据需对22 500个不同参数（150个频点、150个方位角点）的目标RCS进行采样计算，不同方法的耗时如表3所示。本文提出的CPU多核并行技术、GPU硬件加速技术和KDTree遍历技术相结合的靶标高频电磁散射加速方法与传统高频仿真方法相比，将普通性能计算机获取单幅大型舰船与角反射器阵列复合散射SAR图像的耗时由186 h降至0.35 h，可实现大型舰船靶标散射特性的单机快速评估。此时线程调用GPU的时间远小于在CPU端运行的时间，所以通过OpenMP的协调，与不调用OpenMP相比提高了约6倍的计算效率。表3不同方法计算单幅大型舰船与角反射器阵列复合散射SAR图像的耗时对比加速情况无加速GPU+KDTreeCPU+GPU+KDTree总时间/h1861.990.35加速比193531图22为舰船与角反阵列复合散射的一维距离像特性，雷达参数设置如下：俯仰角为60，方位角为45，极化方向为垂直极化，不模糊距离为150 m，分辨率为1 m。图中蓝色点状区域为目标在雷达视线方向的投影示意，依据距离像峰值与目标几何结构的对应关系可知，6个峰值点全部由开口朝向雷达的9个角反射结构三次散射形成。其中，1，4，5号峰值点分别由3对位于同一雷达分辨单元的角反射结构形成，而2，3，6号峰值点则分别由单个角反射器形成。图23为舰船与角反阵列复合模型的SAR图像，图中的9个强散射中心同样全部由这9个角反射结构的三次散射机理形成。通过上述分析可以发现，在雷达观测的一维距离像和SAR图像特性中，边长0.5 m的角反射器产生的强后向散射回波足以掩盖舰船目标自身的回波信号，且由于雷达对角反射器三次散射机理的等效视在位置位于三面角顶点，因此，在舰船周围放置合适尺寸及方向的角反阵列，可以完全改变舰船的雷达特性。0-20123456-40-60050100150200250204060E#C.N,图22三面角阵列与舰船复合模型距离像382黄巍浩：靶标雷达特性快速预测与分析2023 年第 4 期-500500.51.01.52.02.53.03.54.04.55.05.5C.NN图23三面角阵列与舰船复合模型二维图像4结束语本文针对靶标特性建模的需求，提出了结合OpenMP和GPU混合编程的加速算法，提高了获取雷达目标回波数据的效率。对于存在多次作用结构体的目标，本文计算了雷达目标包含高阶散射场的总体回波信息，通过与FEKO对比验证了计算结果的准确性。通过对不同方法计算的用时说明了本文方法的高效性，并利用雷达目标特性数据获得目标的一维距离像和SAR图像，根据射线理论预测一维距离像峰值坐标，解释了距离像峰值的物理意义，雷达图像所展示出的特征与目标几何形状有较好的关联性，说明了计算结果的可靠性。参考文献：1赵京城，娄长玉，李家碧，等.一种方形腔体的目标散射特性及测量方法 J.北京航空航天大学学报，2022，48（12）:24152424.2郭立新，魏仪文.复杂动态海面与目标电磁散射及回波仿真研究现状与展望 J.雷达学报，2023，12（1）:76109.3朱雪慧，余泽太，覃潇潇，等.基于动态RCS的舰船雷达回波仿真与分析 J.雷达科学与技术，2019，17（5）:543549.4孙佳兴，姚慧婧.基于联合处理的复杂目标RCS估计方法 J.雷达科学与技术，2022，20（2）:181186.5王超.高频电磁散射建模方法及工程应用 D.北京：中国传媒大学，2009.6殷红成，郭琨毅.目标电磁散射特性研究的若干热点和难点问题 J.电波科学学报，2020，35（1）:128134.7郝晓军，杨晓帆，赵宏宇，等.基于RCS静态宽带数据的雷达目标成像算法 J.太赫兹科学与电子信息学报,2022，20（3）:208212.8殷红成，朱国庆，董纯柱，等.基于自适应射线管分裂的多次反射计算方法 J.系统工程与电子技术，2013，35（4）:700706.9XIAO Guangliang,GUO Kunyi,SHENG Xinqing.Parametric Scattering Center Modeling for a Conducting Deep Cavity J.IEEE Antennas and Wireless Propagation Letters,2021,20（8）:14191423.10SUK S,SEO T I,PARK H S,et al.Multiresolution GridAlgorithm in the SBR and Its Application to the RCS CalculationJ.Microwave&Optical Technology Letters,2001,29（6）:394397.11GLASSNER A S.Space Subdivision for Fast Ray TracingJ.IEEE Computer Graphics and Applications,1984,4（10）:1524.12MACDONALD J D,BOOTH K S.Heuristics for RayTracing Using Space SubdivisionJ.Visual Computer,1990,6（3）:153166.13SONG Tao,LIAO Huanyu,SUBBARAYAN G.EfficientLocal Refinement near Parametric Boundaries Using KDTree Data Structure and Algebraic Level Sets J.Algorithms,2022,15（7）:245.14TAO Yubo,LIN Hai.GPUBased Shooting and BouncingRay Method for Fast RCS Prediction J.IEEE Trans onAntennas&Propagation,2010,58（2）:494502.15PAJOT A,BARTHE L,PAULIN M,et al.CombinatorialBidirectional PathTracing for Efficient Hybrid CPU/GPURendering J.Computer Graphics Forum,2011,30（2）:315324.16GAO Pengcheng,TAO Yubo,BAI Zhihui,et al.Mappingthe SBR and TWILDCs to Heterogeneous CPUGPU Architecture for Fast Computation of Electromagnetics Scattering J.Progress In Electromagnetics Research,2012,122（1）:137154.作者简介：黄巍浩男，1998年生，在读硕士研究生，主要研究方向为雷达目标特性、高性能计算。（下转第390页）383雷达科学与技术第 21 卷 第 4 期7BUCRIS Y,COHEN I,DORON M A.Bayesian Focusingfor Coherent Wideband BeamformingJ.IEEE Trans onAudio Speech&Language Processing,2012,20（4）:12821296.8CHEN J C,HUDSON R E,KUNG Y.MaximumLikelihood Source Localization and Unknown Sensor LocationEstimation for Wideband Signals in the NearFieldJ.IEEE Trans on Signal Processing,2002,50（8）:1843 1854.9WAX M,SHAN Tiejun,KAILATH T.SpatioTemporalSpectral Analysis by Eigenstructure MethodsJ.IEEETrans on Acoustics,Speech,and Signal Processing,1984,32（4）:817827.10VALAEE S,KABAL P.Wideband Array Processing Using a TwoSided Correlation TransformationJ.IEEETrans on Signal Processing,1995,43（1）:160172.11周豪，胡国平，汪云.基于自适应步长萤火虫多重信号分类算法的低空目标波达方向估计 J.雷达学报，2015，4（3）:309316.12杜永兴，李晨璐，秦岭，等.多径环境下阵列波达方向估计及互耦校正 J.微波学报，2019，35（3）:2433.13李磊，李国林，刘润杰.基于相干积累矩阵重构的波达方向估计新方法 J.雷达学报，2015，4（2）:178184.14刁鸣，安春莲.独立信号和相干信号DOA估计新方法J.系统工程与电子技术，2011，33（12）:25822586.15CHEN Xin,LIU Zhen,WEI Xizhang.Fast FRFTBasedAlgorithm for 3D LFM Source Localization with UniformCircular Array J.IEEE Access,2018,6:21302135.16杨利民，苏卫民，顾红，等.基于二维映射的LFM UWB雷达动目标速度估计 J.雷达学报，2012，1（3）:232237.17JUNG T J,LEE K K.ClosedForm Algorithm for 3D SingleSource Localization with Uniform Circular Array J.IEEE Antennas and Wireless Propagation Letters,2014,13:10961099.18王永良，陈辉，彭应宁，等.空间谱估计理论与算法M.北京：清华大学出版社，2004.作者简介：陈鑫男，1992年生，博士，航天飞行器生存技术与效能评估实验室助理研究员，主要研究方向为阵列信号处理和无源定位。户盼鹤男，1991年生，博士，国防科技大学电子科学学院副教授，主要研究方向为雷达系统与技术、阵列信号处理和无源定位。龚政辉男，1988 年生，国防科技大学电子科学学院助理研究员，主要研究方向为压缩感知和雷达信号处理。苏晓龙男，1994年生，博士，国防科技大学电子科学学院在读博士研究生，主要研究方向为阵列信号处理和深度学习。刘振男，1983年生，博士，国防科技大学电子科学学院教授，主要研究方向为雷达波形设计、雷达抗干扰和无源定位。何思远女，1982年生，博士研究生，博士生导师，主要研究方向为计算电磁学、电磁散射建模。张云华男，1981年生，博士研