弹载SAR回波模拟关键参数并行计算方法.pdf

第30卷 第11期2023年11月仪器仪表用户INSTRUMENTATIONVol.302023 No.11弹载SAR回波模拟关键参数并行计算方法孙厚鹏（北方工业大学 信息学院，北京 100144）摘 要：针对 SAR 雷达回波模拟，本文首先介绍回波模拟信号的产生原理，然后根据回波模拟信号的计算方法，针对其中的关键参数给出具体的计算流程。结合 TMS320C6678 多核 DSP 架构，对系统的任务进行合理划分并实现多核同步，提出了一种基于多核并行的 SAR 雷达回波模拟关键参数计算方案。最后，通过数据的对比，验证本文所提方案的有效性与实时性。关键词：多核并行；回波模拟；合成孔径雷达中图分类号：TN955+.2 文献标志码：AParallel Calculation Method for Key Parameters of Missile Borne SAR Echo SimulationSun Houpeng（School of Information,North China University of Technology,Beijing,100144,China）Abstract：For SAR radar echo simulation,this article first introduces the principle of generating echo simulation signals,and then provides a specific calculation process for the key parameters based on the calculation method of echo simulation signals.Com-bining the TMS320C6678 multi-core DSP architecture,the system tasks are reasonably divided and multi-core synchronization is achieved.A key parameter calculation scheme for SAR radar echo simulation based on multi-core parallel is proposed.Finally,the effectiveness and real-time performance of the proposed scheme are verified through data comparison.Key words：multi-core parallelism；echo simulation；synthetic aperture radar收稿日期：2023-08-25作者简介：孙厚鹏（1998-），男，天津人，硕士，学生，研究方向：雷达信号处理。DOI:10.3969/j.issn.1671-1041.2023.11.004文章编号：1671-1041(2023)11-0012-040 引言弹载 SAR 是一种特殊的，无法重复利用的产品，对其真实性能的检测和验证无法在其实际工作的环境下进行。部分功能性的验证可以在机载挂飞的实验条件下进行，但由于该方法无法模拟导弹真实的飞行条件，其核心性能与系统整体设计的完成度必须依赖地面上的仿真系统，即雷达回波模拟器。雷达回波模拟器可以对导弹真实工作条件下的环境进行模拟，包括导弹的速度、导弹的高度、SAR波束以及 SAR 回波参数等信息，实现对弹载 SAR 的成像功能、高度测量功能、波束指向功能及射频输出信号的功率和波形等指标的测试。本系统的应用场景是弹载 SAR 回波模拟，弹载 SAR雷达照射地面区域范围大，区域变化速度快，计算量很大，但同时又对实时性要求非常高，所以本系统采用了 FPGA与 DSP 相结合的处理方式，二者通过 SRIO 接口传输数据。其中，DSP 负责给 FPGA 提供当前时刻雷达波束照射地面区域的坐标信息，FPGA 负责给 DSP 提供弹道参数并通过DSP 计算的坐标信息读取基准图中的灰度信息（RCS 信息）和高程信息，本文主要讨论 DSP 的处理流程。1 关键参数计算系统流程介绍1.1 回波产生原理合成孔径雷达（SAR）的基本工作原理是通过装载基孙厚鹏弹载SAR回波模拟关键参数并行计算方法第11期13准图等数据包来实时计算回波特征传递函数，并将其与雷达激励信号实时卷积，从而产生合成回波信号。当模拟器接收到当前的弹道参数、雷达参数和基准图数据时，它会实时计算回波信号，以模拟雷达的实际工作情况。在这个过程中，模拟器会使用一系列的数学算法和模型来模拟信号传输和回波信号的形成，以产生准确的模拟结果。地面上某一时刻照射到的某点的回波信号公式为（1）公式（1）中，S 表示为经过时间延迟的雷达激励信号；表示地面上一个点的复反射系数。相对雷达距离为 R 的复反射系数 G 可以表示为（2）式（2）中，A 为复反射系数的幅度，表示该点后向散射系数值（复数 RCS）。反映了由于距离延迟引起回波相位延迟，该项取决于距离 R 和波长。地面上的每个散射点都有自己对应的幅度值和距离。现定义一个一维数组，该数组的序号表示与导弹的距离。将所有散射点累加到数组中，按照距离的不同将散射点分配到不同的数组单元中，可以得到回波特征传递函数，该传递函数按照距离一维排列。公式表现为（3）即雷达激励信号经延迟、调幅、调相后叠加。采用该种方式，实现公式的计算量会很大。式（3）可进一步写为（4）即回波信号为一个回波特征传递函数和雷达激励信号的卷积，回波特征传递函数定义为（5）式（5）中，An为模拟器从基准图中读取获得，为此需要计算所有参与回波计算的散射点的具体位置，获得基准图中散射点的坐标，即块号。1.2 块号计算方法首先，根据当前时刻成像系下雷达的位置、雷达的波束指向和波束宽度确定场景坐标系下照射区域的范围，再根据场景坐标系下照射区域范围内的坐标，经过坐标系变换，得到基准图坐标系下的坐标，即块号坐标。场景坐标系建立的是 LL 的坐标系，一般 L 取值为雷达最大作用距离的两倍，可以认为该坐标系是北天东坐标系，X 为东向，Y 为北向，Z 为天向，所有的辅助数据都是根据此坐标系生成的；基准图坐标系是以基准图的左上角作为（0,0）点的，因此在该坐标系中所有点的坐标值全为非负数；成像系将导弹在基准图上的投影点（弹下点）作为成像系的原点，以导弹速度方向在基准图上的投影方向作为 X 轴方向，Y 轴为天向，Z 轴与 X/Y 轴满足右手法则，且导弹的速度矢量在基准图上的投影以相对于正北方向的航向角来定义。模拟器使用的辅助数据是将场景坐标系中每个散射点的坐标进行处理，根据同心圆算法的要求，预制出的便于DSP 按照距离-角度方式读取坐标信息的数据。利用辅助数据，DSP 可以在确定了雷达波束照射区域距离弹下点最近和最远距离以及照射区域方位向边界与弹下点的夹角之后，直接读取照射区域的坐标值。为了得到基准图坐标系下的块号坐标，还需要进行坐标系变换，将场景坐标系下的坐标转换到基准图坐标系下。具体的坐标系转换规则如下：以成像系 Y 轴为轴旋转成像系，使得成像系的 X 轴和Z 轴分别与基准图的 X轴和 Y轴方向相同，建立新坐标系X1，Y1，Z1。n从 X 轴原点向正向看，左向偏转至正北向（180 范围）角度为正，右向偏转至正北向（180 范围）角度为负，则 XOZ 需要逆时针转 到 X1-O-Z1 系，根据右手系以 Y 轴旋转的坐标旋转公式，平台位置在 X1OZ1 坐标系下的坐标为（6）角从 Y 轴正向向负向看，若逆时针旋转则为正值，顺时针旋转则为负值，因此采用逆时针旋转，则=90-n 0或者顺时针选转，则=360-(n-90)。根据以上公式同样可以算出场景中心点 XC，ZC 在 X1-O-Z1 坐标系下坐标(X1C,Z1C)；由于场景中心点在基准图系 X-O-Y下的坐标为(L/4,L/4)，因而根据坐标系平移公式，算出平台位置在基准图系的坐标为（7）（8）搜索到的照射区域是在场景系的坐标(X2,Y2)，再根据平台位置在场景系的坐标将搜索到的区域平移到基准图系的相对位置(X3,Y3)。（9）（10）第30卷14 仪器仪表用户 INSTRUMENTATION根据上述坐标系变换原则，可以得到基准图坐标系下的块号坐标。2 多核并行搜索方案目前常见的 DSP 并行方法有：数据并行（Data Parallelism）可以通过主从方式实现数据并行，实现方法是通过主处理器进行任务划分，各处理器执行相同的操作。对于数据资源占用较少，相互依赖关系低的数据采用这种方式可以提高处理性能。主从并行方式就是数据并行的主要方法；任务并行（Task Parallelism）可以通过流水方式实现，流水模式是各个核处理同一块数据，每个核处理不同的算法，数据从一个核传给另外一个核用于更深一层次的处理。本系统采用主从方式作为并行设计方法，选取核0 作为主核，核 1 4 作为从核，5 个核进行并行处理。2.1 分段方法为了保证系统函数计算的正确性，搜索区域的长度不能超过系统函数长度与距离分辨率的乘积，但当搜索区域的长度小于系统函数长度与距离分辨率的乘积时，不影响计算的正确性，所以在确定实际参与计算的搜索区域时需要对其适当缩小，使其略小于系统函数长度与距离分辨率的乘积，以保证系统函数计算的正确性。当需要对照射区域进行删减时，以波束中心指向与地面的交点为起点，分别向距离增加和减小的方向移动系统函数长度与距离分辨率乘积的距离的一半，得到照射区域的最远距离和最近距离；根据照射区域距离向长度的缩小换算成距离向波束宽度角度值的缩小，同比例缩小方位向波束宽度的角度值，得到方位向的搜索范围。2.2 核间通信主核和从核可以通过 IPC 进行通信，IPC 中常用的通信模块为 Notify 模块和 MessageQ 模块。本系统采用MessageQ 模块进行多核通信，便于传递搜索参数。主核首先计算好搜索参数，然后向从核发送 MessageQ 消息，消息中包含开始搜索标志位与搜索参数。从核 A 接收到消息后，根据搜索参数确定 A 核的计算任务。当任务计算完成后，向主核发送结束搜索标志位与核心编号 A，用来告知主核已经完成计算任务的从核 ID。主核通过接收到的从核发送的 MessageQ 消息，判断出已经完成任务的核心数量与编号，当所有核心完成计算任务后，主核将所有核心的计算结果打包发送给 FPGA 并等待下一个计算任务。2.3 任务分配本系统中主核的任务是计算搜索参数，给从核分配计算任务，同时负责与 FPGA 进行数据的收发。本系统中从核根据主核分配的计算任务进行计算，任务分配的原则是每个从核平分计算任务，尽可能地缩短计算时间。分发的原则是根据等距离圆的序号进行分发，具体来说，如果需要计算的等距离圆总数为 X，核心数量为M，那么核 m 计算等距离圆序号（m-1）X/M+1）到（m-1）X/M+X/M）之间的数据，每个从核计算的任务量基本相同。3 结果分析3.1 块号搜索结果对比C6678 中搜索块号的结果如图 1 所示，相同的参数用MATLAB 搜索的结果如图 2 所示，将二者的搜索结果用MATLAB 绘制到一起进行对比的结果如图 3 所示。可以看出 C6678 的搜索结果与 MATLAB 的搜索结果完全一致。图1 C6678结果Fig.1 C6678 Results图2 MATLAB结果Fig.2 MATLAB Results图3 DSP与MATLAB结果对比Fig.3 Comparison of DSP and MATLAB results孙厚鹏弹载SAR回波模拟关键参数并行计算方法第11期153.2 单核计算耗时与多核计算耗时对比（时钟数）C6678 单核计算系统耗时见表 1，多核计算系统耗时（主核）见表 2 所示，多核计算系统耗时（从核）见表 3。通过单核计算与多核并行计算的系统耗时对比可以看出，C6678 通过并行计算能够很好地提高系统的计算速度，使系统达到实时性要求。4 结束语本文首先介绍了 SAR 回波模拟中回波信号的产生原理和块号的计算方法，然后介绍了计算块号任务的并行方法，核间通信和任务分配，给出了详细的 DSP 主核和从核处理流程图，最后通过仿真数据与实际结果的对比，以及 DSP运算速度的检测，验证了系统的正确性和实时性，证明了该系统能够满足 SAR 回波模拟关键参数实时计算的技术要求。参考文献：韦闽峰,石瑞军.合成孔径雷达测试方法的设计与实现J.遥测遥控,2013,34(05):54-57,76.表1 单核计算系统耗时Table 1 Time consumption of single core computing system表2 多核计算系统耗时（主核）Table 2 Time consumption of multi-core computing system(main core)表3 多核计算系统耗时（从核）Table 3 Time consumption of multi-core computing system(from core)傅东宁,鲁旭东,范金平.波束照射区域实时搜索方法、装置及电子设备:CN109856607AP.2019-06-07.朱火龙.基于多核DSP的弹载SAR成像信号处理系统设计D.西安:西安电子科技大学,2014.肖亮,张学健.基于多核DSP的超声成像处理算法的并行实现J.电子技术应用,2013,39(06):27-30.Cutrona L J,Hall G O.A comparison of techniques for achieving fine azimuth resolutionJ.IRE Transactions on Military Electronics,1962(02):119-121.23451石油化工自动化智库专家论坛在鄂尔多斯站召开7 月 29 日，由中国自动化产业链创新联合体、中自智库、中自联（北京)科技产业发展有限公司主办，中国神华煤制油化工有限公司鄂尔多斯煤制油分公司联合主办，北京力控元通科技有限公司等 6 家公司承办的以“推动高端化仪表装备应用 助力企业高质量快速发展”为主题的石油化工自动化智库专家论坛鄂尔多斯站在内蒙古鄂尔多斯市召开，150 多人与会交流。中国神华煤制油化工有限公司鄂尔多斯煤制油分公司机动部副经理张宝，内蒙古伊泰集团有限公司调度中心副主任贺澎分别主持会议。中国自动化产业链创新联合体发起人刘哲鸣、中国神华煤制油化工有限公司鄂尔多斯煤制油分公司执行董事胡庆斌致辞。原神华宁煤 400 万吨煤制油仪表运行部部长徐中海等多位专家做主题发言。会议还举行了内蒙古汇中源能源有限公司与北京斯克维思仪表有限公司等 4 家公司的签约仪式。