一种同时SAR/GMT回波模拟方法及其系统转让专利
申请号 : CN202210190103.3
文献号 : CN114252855B
文献日 : 2022-05-06
发明人 : 王娟 , 曹庆刚 , 王九灵 , 哈敏 , 占庆祥
申请人 : 南京雷电信息技术有限公司
摘要 :
权利要求 :
1.一种同时SAR/GMT回波模拟方法,其特征在于,包括:步骤1:在上位机中设置试验参数,试验参数包括SAR场景参数、GMT个数、GMT与SAR场景的初始几何关系、雷达位置初始信息;
步骤2:导入试验参数,运行GPU架构下对SAR/GMT分布式散射面源并行算法,得到独立的以及合并后的SAR/GMT一维距离像数据和回波数据;
所述步骤2中运行GPU架构下对SAR/GMT分布式散射面源并行算法包括:步骤A:读取上位机设置的试验参数,导入雷达GPS航迹解析文件、导入SAR场景后向散射图像,导入GMT的后向散射图像,设计GMT航迹;
步骤B、根据预设的SAR场景、GMT轨迹范围以及雷达信号参数,对GMT后向散射图像进行插值处理,使其与SAR场景的像素比例相融合;
步骤C、调节SAR场景与GMT后向散射系数强度以及相位,使目标的信杂比满足试验需求;
步骤D、将雷达航迹文件、SAR/GMT后向散射数据、GMT航迹文件从CPU搬移到GPU显存中;
步骤E、将SAR场景与GMT按照各自的散射面源个数设置GRID的大小,并将它们独立、离散分布到GPU的计算内核上:
步骤F、利用GPU的并行计算性能,分别对SAR和每个GMT的散射面源同时进行矢量并行计算,计算当前脉冲所有散射面源的幅度相位信息,并判断其所在的距离门;
步骤G、对当前脉冲的SAR场景的所有散射面源的幅度相位信息进行累积得到SAR场景一维距离像信息,同时对当前脉冲的每个GMT分别进行所有散射面源的幅度相位信息累积得到每个GMT的一维距离像信息,合并当前脉冲的SAR/GMT的一维距离像信息;
步骤H:根据导入的试验设置对步骤G产生的一维距离像数据进行卷积处理,生成并保存当前脉冲的SAR场景回波数据、每个GMT的回波数据以及合并后的SAR/GMT回波数据;
步骤I、重复步骤F H,直到所有脉冲的一维距离像数据和回波数据计算完成;
~
步骤3:将试验参数和对应的一维距离像数据和回波数据打包发送至存储板中;
步骤4:上位机输出包含回波类型的试验参数至实时控制系统,实时控制系统根据试验参数实时控制信号处理板卡;
步骤5:信号处理板卡根据回波类型选择调用存储板中的数据,并将选择的数据打包后通过光纤网络输出;
步骤6:若试验参数的回波类型为直接回放,实时控制系统控制信号处理板卡接收雷达参考信号,对光纤网络传输的一维距离像数据进行特征调制,然后进行数字上变频、数字模拟变换,再通过射频上变分系统处理后输出;
步骤7:若试验参数的回波类型为卷积回放,实时控制系统控制信号处理板卡接收光纤网络传输的回波数据,然后对该回波数据进行数字上变频、数字模拟变换,再通过射频上变分系统处理后输出。
2.根据权利要求1所述的同时SAR/GMT回波模拟方法,其特征在于:所述雷达参考信号由外部产生,经射频下变分系统处理后输入至所述信号处理板卡;或由所述实时控制系统根据试验参数输出波形信息至信号处理板卡,信号处理板卡对波形信息实时计算基带波形生成雷达参考信号。
3.根据权利要求2所述的同时SAR/GMT回波模拟方法,其特征在于,所述步骤E中包括:步骤E‑a:据雷达航迹起始位置、雷达波束照射参数、雷达信号参数计算两维分辨率,计算波束照射区域的大小,根据分辨率的一半计算初始波束照射区的距离向 和方位向的散射面源个数 ;
步骤E‑b:根据 、 、 的大小以及GPU计算资源的限制,设置并行计算的线程个数,并将散射面源的计算分布在尽可能多的GPU的计算核心上, 表示合成孔径雷达的单个孔径累计脉冲个数。
4.根据权利要求3所述的同时SAR/GMT回波模拟方法,其特征在于,所述 、 的计算过程如下:
;
其中, 表示雷达到波束照射中心的作用距离,H表示雷达距离地面的高度, 、分别表示方位向、距离向波束宽度,为波束下俯角; 、 分别表示方位向、距离向散射面源的尺寸,计算过程如下
;
其中, 和 分别表示方位向和距离向的分辨率,计算过程如下;
其中, 表示雷达平台飞行速度, 表示雷达发射信号工作频率, 表示光速单位m/s, 表示雷达发射信号的带宽;
选取 、 的最小值作为散射面源的尺寸:所述步骤E‑b中合成孔径雷达的单个孔径累计脉冲个数的 计算过程如下:其中, 表示脉冲重复周期。
5.根据权利要求4所述的同时SAR/GMT回波模拟方法,其特征在于,所述步骤F包括:对于SAR散射面源,SAR场景的计算输入参数有场景后向散射系数、雷达航迹数据、雷达波束照射参数、以及场景散射面源的位置坐标,其计算公式如下:;
以雷达第一个航迹位置时的波束照射中心建立北天东直角坐标系,其中, 表示第p个脉冲北向第n个天向第a个东向第e个散射面源的SAR的幅度相位信息; 表示第p个脉冲北向第n个天向第a个东向第e个散射面源的后向散射系数; 表示雷达发射信号的波长; 、 、 分别表示雷达平台第p个脉冲时刻的北、天、东位置; 、、 分别表示SAR场景中北向第n个天向、第a个东向、第e个散射面源的位置坐标,j是指复数运算中取虚数部分;
对于GMT散射面源,GMT的计算输入参数有场景后向散射系数、雷达航迹数据、雷达波束照射参数、GMT插值后的散射面源的位置坐标,每个GMT的航迹,其计算公式如下:;
其中, 表示第p个脉冲第i个GMT的北向第n个天向第a个东向第e个散射面源幅度相位信息; 表示第p个脉冲第i个GMT的北向第n个天向第a个东向第e个散射面源的后向散射系数; 、 、 分别表示第i个GMT的中北向第n个天向、第a个东向、第e个散射面源的位置坐标,这三个位置信息根据目标运动轨迹插值得到。
6.根据权利要求5所述的同时SAR/GMT回波模拟方法,其特征在于,所述步骤G中合并前脉冲的SAR/GMT的一维距离像信息过程如下:。
7.根据权利要求6所述的同时SAR/GMT回波模拟方法,其特征在于,所述步骤H中当前脉冲的SAR场景回波数据 、每个GMT的回波数据 以及合并后的SAR/GMT回波数据 分别如下:
;
;
其中,LFM为雷达参考信号。
8.根据权利要求7所述的同时SAR/GMT回波模拟方法,其特征在于:所述步骤I的计算过程如下:
。
9.一种同时SAR/GMT回波模拟系统,其特征在于:包括上位机、实时控制系统、GPU加速模拟单元,存储板、射频上变分系统、射频下变分系统、信号处理板卡,信号处理板卡与实时控制系统、射频上变分系统、射频下变频分系统、存储板分别相连;所述上位机用于设置试验参数,试验参数包括SAR场景参数、GMT个数、GMT与SAR场景的初始几何关系、雷达位置初始信息;所述GPU加速模拟单元用于接收试验参数、调用试验参数,并通过以下步骤A至步骤I计算得到独立的以及合并后的SAR/GMT一维距离像数据和回波数据并存储至所述存储板;
步骤A:读取上位机设置的试验参数,导入雷达GPS航迹解析文件、导入SAR场景后向散射图像,导入GMT的后向散射图像,设计GMT航迹;
步骤B、根据预设的SAR场景、GMT轨迹范围以及雷达信号参数,对GMT后向散射图像进行插值处理,使其与SAR场景的像素比例相融合;
步骤C、调节SAR场景与GMT后向散射系数强度以及相位,使目标的信杂比满足试验需求;
步骤D、将雷达航迹文件、SAR/GMT后向散射数据、GMT航迹文件从CPU搬移到GPU显存中;
步骤E、将SAR场景与GMT按照各自的散射面源个数设置GRID的大小,并将它们独立、离散分布到GPU的计算内核上:
步骤F、利用GPU的并行计算性能,分别对SAR和每个GMT的散射面源同时进行矢量并行计算,计算当前脉冲所有散射面源的幅度相位信息,并判断其所在的距离门;
步骤G、对当前脉冲的SAR场景的所有散射面源的幅度相位信息进行累积得到SAR场景一维距离像信息,同时对当前脉冲的每个GMT分别进行所有散射面源的幅度相位信息累积得到每个GMT的一维距离像信息,合并当前脉冲的SAR/GMT的一维距离像信息;
步骤H:根据导入的试验设置对步骤G产生的一维距离像数据进行卷积处理,生成并保存当前脉冲的SAR场景回波数据、每个GMT的回波数据以及合并后的SAR/GMT回波数据;
步骤I、重复步骤F H,直到所有脉冲的一维距离像数据和回波数据计算完成;
~
所述实时控制系统用于接受上位机输出的试验参数并输出控制信息至所述信号处理板卡;所述信号处理板卡接收射频下变分系统输入的外部雷达参考信号或根据实时控制系统输出的控制信息产生对应的雷达参考信号,读取存储卡的一维距离像数据和回波数据,按照回波类型进行输出选择:将一维距离像数据和回波数据进行数据打包后通过光纤网络发送到指定终端;直接回放,实时读取回波数据并进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出;卷积回放,对一维距离像数据进行特征调制后进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出。
说明书 :
一种同时SAR/GMT回波模拟方法及其系统
技术领域
背景技术
气和遮盖物的影响,可以穿透云层和植被探测地物信息,是地形地物探测装备的最重要和
主要的手段。几十年来越来越广泛地应用到民用、军事领域,是各个国家雷达技术和信号处
理技术的重点研究方向。
SAR/GMT回波模拟系统采用GPU架构,将SAR场景和GMT进行分离处理,利用散射面元计算法
同时计算SAR场景和GMT的回波,最后对SAR场景和分布式GMT回波进行合成得到总的SAR回
波。SAR/GMT回波模拟系统相比于常规的SAR回波模拟系统应用场景更加广泛,相比对多通
道GMT回波模拟系统它减少了硬件资源,极大地节约了设备成本。
回波模拟系统主要用于在SAR/GMTI雷达上,用于检测该类雷达在SAR图像上显示运动目标
的位置,检测运动目标个数、对运动目标进行测速、对运动目标进行定位等的性能测试。但
SAR回波信号仿真的计算量非常大,加上多个运动目标的回波信号仿真,计算量倍增,无法
满足SAR/GMTI雷达性能测试需求。
发明内容
波信号,应用场景更加广泛,极大地节约了设备成本。
字模拟变换,再通过射频上变分系统处理后输出;
上变分系统处理后输出。
处理板卡,信号处理板卡对波形信息实时计算基带波形生成雷达参考信号。
累积得到每个GMT的一维距离像信息,合并当前脉冲的SAR/GMT的一维距离像信息;
向的散射面源个数 ;
单个孔径累计脉冲个数。
距离向散射面源的尺寸,计算过程如下
第p个脉冲北向第n个天向第a个东向第e个散射面源的后向散射系数; 表示雷达发射
信号的波长; 、 、 分别表示雷达平台第p个脉冲时刻的北、天、东位置; 、
、 分别表示SAR场景中北向第n个天向、第a个东向、第e个散射面源的位置坐
标,j是指复数运算中取虚数部分;
的后向散射系数; 、 、 分别表示第i个GMT的中北向第n个天向、第a个东
向、第e个散射面源的位置坐标,这三个位置信息根据目标运动轨迹插值得到。
板卡与实时控制系统、射频上变分系统、射频下变频分系统、存储板分别相连;所述上位机
用于设置试验参数;所述GPU加速模拟单元用于接收试验参数、调用试验参数,计算得到独
立的以及合并后的SAR/GMT一维距离像数据和回波数据并存储至所述存储板;所述实时控
制系统用于接受上位机输出的试验参数并输出控制信息至所述信号处理板卡;
据,按照回波类型进行输出选择:将一维距离像数据和回波数据进行数据打包后通过光纤
网络发送到指定终端;直接回放,实时读取回波数据并进行数字上变频、数字模拟变换,再
经所述射频上变分系统处理后输出;卷积回放,对一维距离像数据进行特征调制后进行数
字上变频、数字模拟变换,再经所述射频上变分系统处理后输出。
到独立的以及合并后SAR/GMT一维距离像数据和回波数据,这几组数据通过光纤接口/万兆
以太网传输到存储板卡中,待试验开始时,通过PCIe总线将数据传输到信号处理板卡中,根
据试验设置选择回波数据对其进行DUC和DA处理或者选择一维距离像数据进行目标特征调
制再对回波做DUC和DA,最后进行上变频生成SAR/GMT回波信号。
SAR/GMT合成的SAR雷达回波信号,应用场景更加广泛,相比对多通道GMT回波模拟系统减少
了硬件资源,极大地节约了设备成本;本系统在GPU平台上采用分布式散射面源并行算法大
大提高了同时SAR/GMT回波模拟的速度,将运动目标模拟与SAR场景模拟分离开来,采用了
场景目标分离式的处理方法离线计算数据然后实时回放的方式模拟回波信号,降低了并行
计算的计算量且有利于系统扩展。
据。本系统采用GPU并行计算架构计算SAR场景和GMT的一维距离像数据/回波数据,采用
sysBios的实时操作系统与上位机进行光纤通讯,采用PCIe高速传输接口与信号处理板卡
和存储板卡进行数据传输。这种方法有利于控制运动目标个数,实时切换雷达工作模式、适
应任意雷达发射信号波形,使得回波模拟能够适应SAR/GMTI雷达性能测试需求,其一维距
离像的计算加速比可达10的三次方量级。
附图说明
具体实施方式
件数据对散射元进行并行加速计算,得到独立的以及合并后的SAR/GMT一维距离像数据和
回波数据,这几组数据通过光纤接口/万兆以太网传输到存储板卡中。待试验开始时,通过
PCIe总线将数据传输到FPGA信号处理板卡中。信号处理板卡功能包括:根据雷达参考信号
对一维距离像数据进行特征调制,实时回放目标回波数据,对回波数据进行DUC及DA变换。
回波类型输出选择由上位机试验设置,发送到实时控制系统接收,然后控制转发至信号处
理板卡进行实现。这里的雷达参考信号是外部的输入,它不由本系统产生,当没有外部输入
的这个雷达参考信号时,则由上位机设置雷达波形参数,并发送给实时控制系统交由信号
处理板卡模拟产生一个理想的雷达参考信号替代外部雷达参考信号。
信号;2‑卷积回放,这种设置条件下选择一维距离像数据进行目标特征调制然后再进行DUC
+DA处理即可输出中频回波信号,最后进行上变频生成SAR/GMT回波信号。
NVIDIA 公司的 RTX A2000系列的显卡,其CUDA核心数为3328个,每个SM上运行的线程个数
最大为256。
GPU有多少可控并发执行的线程就能对算法进行相应数量的加速,但由于SAR场景的计算量
级是4096×4096×4096,可并发执行的线程显然是不够用的,下面按照分布式散射面源并
行算法将散射面源的计算自动分布到并发的线程中,其算法计算流程如下:
面源个数:
表示雷达到波束照射中心的作用距离,表示雷达距离地面的高度,为波束下俯角。
的最小值作为散射面源的尺寸:
计算分布在尽可能多的GPU的计算核心上;
SAR的幅度相位信息; 表示第p个脉冲北向第n个天向第a个东向第e个散射面源的后向
散射系数; 表示雷达发射信号的波长; 、 、 分别表示雷达平台第p个脉冲时刻
的北、天、东位置; 、 、 分别表示SAR场景中北向第n个天向、第a个东
向、第e个散射面源的位置坐标;
向散射系数; 、 、 分别表示第i个GMT的中北向第n个天向、第a个东向、
第e个散射面源的位置坐标,这个三个位置信息根据目标运动轨迹插值得到;
积得到每个GMT的一维距离像信息;合并当前脉冲的SAR/GMT的一维距离像信息:
过光纤网络发送到指定终端;
数据进行流水卷积调制;或对回波数据的逐脉冲读取,最后将回波数据进行DUC及DA变换得
到中频回波信号;
统扩展;能够在单通道中模拟同时SAR/GMT合成的SAR雷达回波信号,应用场景更加广泛,相
比对多通道GMT回波模拟系统减少了硬件资源,极大地节约了设备成本。
其在形式上和细节上作出各种变化。