CN114252855B - 一种同时sar/gmt回波模拟方法及其系统 - Google Patents
一种同时sar/gmt回波模拟方法及其系统 Download PDFInfo
- Publication number
- CN114252855B CN114252855B CN202210190103.3A CN202210190103A CN114252855B CN 114252855 B CN114252855 B CN 114252855B CN 202210190103 A CN202210190103 A CN 202210190103A CN 114252855 B CN114252855 B CN 114252855B
- Authority
- CN
- China
- Prior art keywords
- gmt
- sar
- radar
- echo
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012360 testing method Methods 0.000 claims abstract description 59
- 238000012545 processing Methods 0.000 claims abstract description 47
- 238000006243 chemical reaction Methods 0.000 claims abstract description 37
- 239000013307 optical fiber Substances 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 31
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000012856 packing Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000036278 prepulse Effects 0.000 claims description 2
- 238000011056 performance test Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical group OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开同时SAR/GMT回波模拟方法及其系统,包括:在上位机中设置试验参数;导入试验参数,运行GPU架构下对SAR/GMT分布式散射面源并行算法,得到独立的以及合并后SAR/GMT一维距离像数据和回波数据并发送至存储板中;上位机输出包含回波类型的试验参数至实时控制系统;信号处理板卡根据回波类型选择调用存储板中的数据后通过光纤网络输出;根据试验参数的回波类型,实时控制系统控制信号处理板卡对一维距离像数据或回波数据进行处理,然后进行DUC、DA变换,再通过射频上变分系统处理后输出。本发明采用GPU并行计算架构,能够在单通道中模拟同时SAR/GMT合成的SAR雷达回波信号。
Description
技术领域
本发明涉及雷达性能测试技术领域,具体涉及一种同时SAR/GMT回波模拟方法及其系统。
背景技术
合成孔径雷达(Synthtic Aperture Radar,SAR)是一种利用脉冲压缩技术和合成孔径原理获取两维(距离维和方位维)高分辨率的雷达系统。相比于光学雷达系统它不受天气和遮盖物的影响,可以穿透云层和植被探测地物信息,是地形地物探测装备的最重要和主要的手段。几十年来越来越广泛地应用到民用、军事领域,是各个国家雷达技术和信号处理技术的重点研究方向。
在常规的SAR图像中,地面运动目标(Ground Moving Target,GMT)会出现距离走动、散焦、方位位置偏移等现象,GMT往往只能以模糊的形式叠加在静止场景SAR图像上。SAR/GMT回波模拟系统采用GPU架构,将SAR场景和GMT进行分离处理,利用散射面元计算法同时计算SAR场景和GMT的回波,最后对SAR场景和分布式GMT回波进行合成得到总的SAR回波。SAR/GMT回波模拟系统相比于常规的SAR回波模拟系统应用场景更加广泛,相比对多通道GMT回波模拟系统它减少了硬件资源,极大地节约了设备成本。
同时SAR/GMT回波模拟系统是检验雷达对环境感知性能的重要装置,特别对具有地面运动目标检测(Ground Moving Target Indication,GMTI)功能的雷达具有重要意义。回波模拟系统主要用于在SAR/GMTI雷达上,用于检测该类雷达在SAR图像上显示运动目标的位置,检测运动目标个数、对运动目标进行测速、对运动目标进行定位等的性能测试。但SAR回波信号仿真的计算量非常大,加上多个运动目标的回波信号仿真,计算量倍增,无法满足SAR/GMTI雷达性能测试需求。
发明内容
发明目的:本发明目的在于针对现有技术的不足,提供一种同时SAR/GMT回波模拟方法及其系统,采用GPU并行计算架构,能够在单通道中模拟同时SAR/GMT合成的SAR雷达回波信号,应用场景更加广泛,极大地节约了设备成本。
技术方案:本发明提供的同时SAR/GMT回波模拟方法,包括:
步骤1:在上位机中设置试验参数;
步骤2:导入试验参数,运行GPU架构下对SAR/GMT分布式散射面源并行算法,得到独立的以及合并后的SAR/GMT一维距离像数据和回波数据;
步骤3:将试验参数和对应的一维距离像数据和回波数据打包发送至存储板中;
步骤4:上位机输出包含回波类型的试验参数至实时控制系统,实时控制系统根据试验参数实时控制信号处理板卡;
步骤5:信号处理板卡根据回波类型选择调用存储板中的数据,并将选择的数据打包后通过光纤网络输出;
步骤6:若试验参数的回波类型为直接回放,实时控制系统控制信号处理板卡接收雷达参考信号,对光纤网络传输的一维距离像数据进行特征调制,然后进行数字上变频、数字模拟变换,再通过射频上变分系统处理后输出;
步骤7:若试验参数的回波类型为卷积回放,实时控制系统控制信号处理板卡接收光纤网络传输的回波数据,然后对该回波数据进行数字上变频、数字模拟变换,再通过射频上变分系统处理后输出。
进一步完善上述技术方案,所述雷达参考信号由外部产生,经射频下变分系统处理后输入至所述信号处理板卡;或由所述实时控制系统根据试验参数输出波形信息至信号处理板卡,信号处理板卡对波形信息实时计算基带波形生成雷达参考信号。
进一步地,所述步骤1中的试验参数包括SAR场景参数、GMT个数、GMT与SAR场景的初始几何关系、雷达位置初始信息;
所述步骤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,直到所有脉冲的一维距离像数据和回波数据计算完成。
进一步地,所述步骤E中包括:
进一步地,所述步骤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个散射面源的位置坐标,这三个位置信息根据目标运动轨迹插值得到。
进一步地,所述步骤G中合并前脉冲的SAR/GMT的一维距离像信息过程如下:
其中,LFM为雷达参考信号。
所述步骤I的计算过程如下:
用于实现上述同时SAR/GMT回波模拟方法的模拟系统,包括上位机、实时控制系统、GPU加速模拟单元,存储板、射频上变分系统、射频下变分系统、信号处理板卡,信号处理板卡与实时控制系统、射频上变分系统、射频下变频分系统、存储板分别相连;所述上位机用于设置试验参数;所述GPU加速模拟单元用于接收试验参数、调用试验参数,计算得到独立的以及合并后的SAR/GMT一维距离像数据和回波数据并存储至所述存储板;所述实时控制系统用于接受上位机输出的试验参数并输出控制信息至所述信号处理板卡;
所述信号处理板卡接收射频下变分系统输入的外部雷达参考信号或根据实时控制系统输出的控制信息产生对应的雷达参考信号,读取存储卡的一维距离像数据和回波数据,按照回波类型进行输出选择:将一维距离像数据和回波数据进行数据打包后通过光纤网络发送到指定终端;直接回放,实时读取回波数据并进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出;卷积回放,对一维距离像数据进行特征调制后进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出。
在本发明中,先将预设的试验参数放在上位机的数据管理中心,GPU矢量并行计算程序调用该接口文件对SAR场景和GMT数据、航迹文件数据对散射元进行并行加速计算,得到独立的以及合并后SAR/GMT一维距离像数据和回波数据,这几组数据通过光纤接口/万兆以太网传输到存储板卡中,待试验开始时,通过PCIe总线将数据传输到信号处理板卡中,根据试验设置选择回波数据对其进行DUC和DA处理或者选择一维距离像数据进行目标特征调制再对回波做DUC和DA,最后进行上变频生成SAR/GMT回波信号。
有益效果:与现有技术相比,本发明的优点在于:本发明针对SAR/GMTI雷达性能测试需求而提供的一种回波仿真模拟系统,采用GPU并行计算架构,能够在单通道中模拟同时SAR/GMT合成的SAR雷达回波信号,应用场景更加广泛,相比对多通道GMT回波模拟系统减少了硬件资源,极大地节约了设备成本;本系统在GPU平台上采用分布式散射面源并行算法大大提高了同时SAR/GMT回波模拟的速度,将运动目标模拟与SAR场景模拟分离开来,采用了场景目标分离式的处理方法离线计算数据然后实时回放的方式模拟回波信号,降低了并行计算的计算量且有利于系统扩展。
考虑到试验过程中在轨迹不同时刻存在工作模式切换的可能,因而需要在离线计算数据时,把不同时刻不同模式的数据都准备完成,以便试验中根据指令来切换回放的数据。本系统采用GPU并行计算架构计算SAR场景和GMT的一维距离像数据/回波数据,采用sysBios的实时操作系统与上位机进行光纤通讯,采用PCIe高速传输接口与信号处理板卡和存储板卡进行数据传输。这种方法有利于控制运动目标个数,实时切换雷达工作模式、适应任意雷达发射信号波形,使得回波模拟能够适应SAR/GMTI雷达性能测试需求,其一维距离像的计算加速比可达10的三次方量级。
附图说明
图1是同时SAR/GMT回波模拟系统框图;
图2是基于GPU的SAR/GMT分布式散射面源并行算法计算流程图;
图3是SAR场景图像;
图4是GMT图像及运动轨迹;
图5是某时刻同时SAR/GMT合成图像。
具体实施方式
下面通过附图对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
如图1所示,本发明提供的同时SAR/GMT回波模拟系统,先将预设的试验参数放在上位机的数据管理中心,GPU矢量并行计算程序调用接口文件对SAR场景、GMT数据、航迹文件数据对散射元进行并行加速计算,得到独立的以及合并后的SAR/GMT一维距离像数据和回波数据,这几组数据通过光纤接口/万兆以太网传输到存储板卡中。待试验开始时,通过PCIe总线将数据传输到FPGA信号处理板卡中。信号处理板卡功能包括:根据雷达参考信号对一维距离像数据进行特征调制,实时回放目标回波数据,对回波数据进行DUC及DA变换。回波类型输出选择由上位机试验设置,发送到实时控制系统接收,然后控制转发至信号处理板卡进行实现。这里的雷达参考信号是外部的输入,它不由本系统产生,当没有外部输入的这个雷达参考信号时,则由上位机设置雷达波形参数,并发送给实时控制系统交由信号处理板卡模拟产生一个理想的雷达参考信号替代外部雷达参考信号。
信号处理板卡根据试验设置选择回波数据,包括两种类型:1-回波数据直接回放,这种设置条件下对该数据进行数字上变频+数字模拟转换(DUC+DA)处理即可输出中频回波信号;2-卷积回放,这种设置条件下选择一维距离像数据进行目标特征调制然后再进行DUC+DA处理即可输出中频回波信号,最后进行上变频生成SAR/GMT回波信号。
基于同时SAR/GMT回波模拟系统的模拟方法,包括:
步骤1:在上位机中设置试验初始状态:SAR场景参数,GMT个数,GMT与SAR场景的初始几何关系,雷达位置初始信息,
步骤2:导入试验参数,运行GPU架构的下对SAR/GMT分布式散射面源并行算法,生成多组待试验的SAR/GMT合路一维距离像数据/回波数据,主要步骤如下:
A、读取上位机试验参数,导入雷达GPS航迹解析文件、导入SAR场景后向散射图像,导入GMT的后向散射图像,设计GMT航迹;
B、根据预设的SAR场景参数、GMT航迹范围,以及雷达信号参数,对GMT后向散射图像进行插值处理,使其与SAR场景的像素比例相融合;
C、调节SAR场景与GMT的后向散射系数强度以及相位,使GMT的信杂比满足测试需求;
D、将雷达航迹文件、SAR/GMT后向散射数据、GMT的航迹文件从CPU搬移到GPU显存中;
E、将SAR场景与GMT按照各自的散射面源个数设置GRID的大小,并将它们独立、离散分布到GPU的计算内核上:
F、假设SAR场景的散射面源个数为4096×4096,GMT数目为2个,GMT插值后的散射面源个数均为32*32,累计方位向脉冲个数为4096个,脉冲重复周期为1000us。这里选用NVIDIA 公司的 RTX A2000系列的显卡,其CUDA核心数为3328个,每个SM上运行的线程个数最大为256。
将SAR场景的散射面源与每个GMT的散射面源分别独立地加载到不同的GPU计算核心上,同时并发地所有面源的一维距离像进行计算。理论上,在不受共享内存限制的条件下GPU有多少可控并发执行的线程就能对算法进行相应数量的加速,但由于SAR场景的计算量级是4096×4096×4096,可并发执行的线程显然是不够用的,下面按照分布式散射面源并行算法将散射面源的计算自动分布到并发的线程中,其算法计算流程如下:
a、根据雷达航迹起始位置、雷达波束照射参数、雷达信号参数计算两维分辨率,计算波束照射区域的大小,根据分辨率的一半计算初始波束照射区的距离向和方位向的散射面源个数:
其中、分别表示方位向和距离向散射面源的个数,、分别表示方位
向、距离向波束宽度,、分别表示方位向、距离向散射面源的尺寸,单位m,
表示雷达到波束照射中心的作用距离,表示雷达距离地面的高度,为波束下俯角。
理论上,
b、根据、、(表示合成孔径雷达的单个孔径累计脉冲个数)的大小、GPU计
算资源的限制,设置并行计算的线程个数(须为2的幂指数个,且大于128),并将散射面源的
计算分布在尽可能多的GPU的计算核心上;
c、利用GPU的并行计算性能,分别对SAR和各个GMT的散射面源同时进行矢量并行计算,计算当前脉冲所有散射面源的幅度相位特征信息,并判断其所在的距离门:
d、如图2所示,对于SAR散射面源,SAR场景的计算输入参数有场景后向散射系数、雷达航迹数据、雷达波束照射参数、以及场景散射面源的位置坐标,其计算公式如下:
公式中j表示复数中的虚部,以雷达第一个航迹位置时的波束照射中心建立北天
东直角坐标系,其中,表示第p个脉冲北向第n个天向第a个东向第e个散射面源的
SAR的幅度相位信息;表示第p个脉冲北向第n个天向第a个东向第e个散射面源的后向
散射系数;表示雷达发射信号的波长;、、分别表示雷达平台第p个脉冲时刻
的北、天、东位置;、、分别表示SAR场景中北向第n个天向、第a个东
向、第e个散射面源的位置坐标;
对于GMT散射面源,GMT的计算输入参数有场景后向散射系数、雷达航迹数据、雷达波束照射参数、目标插值后的散射面源的位置坐标、每个GMT的航迹,其计算公式如下:
表示第p个脉冲第i个GMT的北向第n个天向第a个东向第e个散射面源幅度
相位信息;表示第p个脉冲第i个GMT的北向第n个天向第a个东向第e个散射面源的后
向散射系数;、、分别表示第i个GMT的中北向第n个天向、第a个东向、
第e个散射面源的位置坐标,这个三个位置信息根据目标运动轨迹插值得到;
G、对当前脉冲的SAR场景的所有散射面源的幅度相位信息进行累积得到SAR场景一维距离像信息,同时对当前脉冲的所有GMT的所有散射面源的幅度相位信息分别进行累积得到每个GMT的一维距离像信息;合并当前脉冲的SAR/GMT的一维距离像信息:
H、根据导入的试验设置对上一步骤产生的一维距离像数据进行卷积处理,生成并保存当前脉冲的SAR场景回波数据、每个GMT的回波数据以及合并后的SAR/GMT回波数据:
其中,LFM为雷达参考信号。
I、重复步骤F~H,直到所有脉冲的一维距离像数据和回波数据计算完成;
步骤3:将试验参数和对应的一维距离像数据/回波数据打包发送到存储板中;
步骤4:上位机输出试验参数至实时控制系统,实时控制系统根据试验参数实时控制信号处理板卡;
步骤5:在sysBios实时控制系统中,根据试验参数的回波类型输出选择调用存储板中的一维距离像数据/回波模拟数据发送到信号处理板卡,信号处理板卡将数据打包通过光纤网络发送到指定终端;
步骤6:信号处理板卡根据实时控制系统的实时控制命令,接收外部输入的雷达参考信号或根据实时控制系统输出波形信息生成基带波形,然后按照工作模式对一维距离像数据进行流水卷积调制;或对回波数据的逐脉冲读取,最后将回波数据进行DUC及DA变换得到中频回波信号;
步骤7:中频回波信号注入到射频上变分系统中,经上变频和功率放大等处理后形成同时SAR/GMT射频回波信号在射频端输出,如图5所示。
本发明采用了分布式散射面源并行算法,大大提高了同时SAR/GMT回波模拟的速度的同时,将运动目标模拟与SAR场景模拟分离开来,降低了并行计算的计算量且有利于系统扩展;能够在单通道中模拟同时SAR/GMT合成的SAR雷达回波信号,应用场景更加广泛,相比对多通道GMT回波模拟系统减少了硬件资源,极大地节约了设备成本。
如上所述,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (9)
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回波模拟方法,其特征在于:所述雷达参考信号由外部产生,经射频下变分系统处理后输入至所述信号处理板卡;或由所述实时控制系统根据试验参数输出波形信息至信号处理板卡,信号处理板卡对波形信息实时计算基带波形生成雷达参考信号。
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的航迹,其计算公式如下:
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,直到所有脉冲的一维距离像数据和回波数据计算完成;
所述实时控制系统用于接受上位机输出的试验参数并输出控制信息至所述信号处理板卡;所述信号处理板卡接收射频下变分系统输入的外部雷达参考信号或根据实时控制系统输出的控制信息产生对应的雷达参考信号,读取存储卡的一维距离像数据和回波数据,按照回波类型进行输出选择:将一维距离像数据和回波数据进行数据打包后通过光纤网络发送到指定终端;直接回放,实时读取回波数据并进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出;卷积回放,对一维距离像数据进行特征调制后进行数字上变频、数字模拟变换,再经所述射频上变分系统处理后输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210190103.3A CN114252855B (zh) | 2022-03-01 | 2022-03-01 | 一种同时sar/gmt回波模拟方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210190103.3A CN114252855B (zh) | 2022-03-01 | 2022-03-01 | 一种同时sar/gmt回波模拟方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114252855A CN114252855A (zh) | 2022-03-29 |
CN114252855B true CN114252855B (zh) | 2022-05-06 |
Family
ID=80800153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210190103.3A Active CN114252855B (zh) | 2022-03-01 | 2022-03-01 | 一种同时sar/gmt回波模拟方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114252855B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116609742B (zh) * | 2023-06-05 | 2023-09-22 | 广州博远装备科技有限公司 | 基于实时地图反演sar回波模拟与干扰评估方法及系统 |
CN117908026B (zh) * | 2024-02-27 | 2025-01-03 | 南京雷电信息技术股份有限公司 | 一种用于地面雷达威力测试的外场试验方法 |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0042785A1 (fr) * | 1980-06-20 | 1981-12-30 | ETAT-FRANCAIS représenté par le Délégué Général pour l' Armement | Procédé et dispositif pour la simulation d'échos radar |
DE3248879A1 (de) * | 1982-06-18 | 1984-04-12 | Dornier System Gmbh, 7990 Friedrichshafen | Verfahren und vorrichtung zur erzeugung kuenstlicher zielmarken in der abbildung eines radars mit synthetischer apertur (sar) |
US4523196A (en) * | 1981-03-27 | 1985-06-11 | Dornier System Gmbh | Test equipment for a synthetic aperture radar system |
WO2008041252A1 (en) * | 2006-10-06 | 2008-04-10 | Space Engineering S.P.A. | Active transponder, particularly for synthetic aperture radar, or sar, systems |
CN101295019A (zh) * | 2008-05-23 | 2008-10-29 | 北京航空航天大学 | 一种线性调频连续波合成孔径雷达视频信号的仿真方法 |
CN103176170A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种基于gpu并行计算的sar回波模拟方法 |
EP2762918A2 (en) * | 2014-01-16 | 2014-08-06 | Institute of Electronics, Chinese Academy of Sciences | Method for analyzing effect of sub-band interference on imaging performance in synthetic aperture radar |
CN104076341A (zh) * | 2014-07-04 | 2014-10-01 | 西安电子科技大学 | 基于fpga和dsp的雷达回波实时模拟方法 |
CN104569938A (zh) * | 2015-02-09 | 2015-04-29 | 中国科学院电子学研究所 | 一种合成孔径雷达回波仿真器 |
CN104698445A (zh) * | 2015-03-30 | 2015-06-10 | 北京润科通用技术有限公司 | 一种雷达回波模拟方法及雷达回波信号模拟器 |
CN105911532A (zh) * | 2016-06-29 | 2016-08-31 | 北京化工大学 | 基于深度协同的合成孔径雷达回波并行模拟方法 |
CN105954730A (zh) * | 2016-06-21 | 2016-09-21 | 电子科技大学 | 一种sar回波快速时域生成方法 |
CN107229051A (zh) * | 2017-05-26 | 2017-10-03 | 西安电子科技大学 | 基于gpu的视频sar回波仿真并行实现方法 |
WO2018004503A1 (en) * | 2016-07-01 | 2018-01-04 | Havelsan Hava Elektronik Sanayi Ve Ticaret Anonim Sirketi | Radar video simulator |
CN107765226A (zh) * | 2017-09-18 | 2018-03-06 | 北京空间飞行器总体设计部 | 一种sar卫星雷达回波模拟方法、系统和介质 |
CN108519590A (zh) * | 2018-03-26 | 2018-09-11 | 北京理工大学 | 激光成像雷达目标回波信号模拟方法及模拟器 |
CN109901165A (zh) * | 2019-03-28 | 2019-06-18 | 河南九乾电子科技有限公司 | 星载sar回波的模拟装置及模拟方法 |
CN110515050A (zh) * | 2019-08-30 | 2019-11-29 | 北京空间飞行器总体设计部 | 一种基于gpu的星载sar实时回波模拟器 |
CN209821374U (zh) * | 2019-04-02 | 2019-12-20 | 南京雷电信息技术有限公司 | 一种雷达目标模拟器 |
CN112698280A (zh) * | 2020-12-09 | 2021-04-23 | 南京长峰航天电子科技有限公司 | 一种基于dsp和fpga架构的双基sar实时回波模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002334708A1 (en) * | 2001-10-01 | 2003-04-14 | Kline And Walker, Llc | Pfn/trac system faa upgrades for accountable remote and robotics control |
CN111025246B (zh) * | 2019-11-28 | 2021-09-07 | 北京遥测技术研究所 | 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法 |
-
2022
- 2022-03-01 CN CN202210190103.3A patent/CN114252855B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0042785A1 (fr) * | 1980-06-20 | 1981-12-30 | ETAT-FRANCAIS représenté par le Délégué Général pour l' Armement | Procédé et dispositif pour la simulation d'échos radar |
US4523196A (en) * | 1981-03-27 | 1985-06-11 | Dornier System Gmbh | Test equipment for a synthetic aperture radar system |
DE3248879A1 (de) * | 1982-06-18 | 1984-04-12 | Dornier System Gmbh, 7990 Friedrichshafen | Verfahren und vorrichtung zur erzeugung kuenstlicher zielmarken in der abbildung eines radars mit synthetischer apertur (sar) |
WO2008041252A1 (en) * | 2006-10-06 | 2008-04-10 | Space Engineering S.P.A. | Active transponder, particularly for synthetic aperture radar, or sar, systems |
CN101295019A (zh) * | 2008-05-23 | 2008-10-29 | 北京航空航天大学 | 一种线性调频连续波合成孔径雷达视频信号的仿真方法 |
CN103176170A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种基于gpu并行计算的sar回波模拟方法 |
EP2762918A2 (en) * | 2014-01-16 | 2014-08-06 | Institute of Electronics, Chinese Academy of Sciences | Method for analyzing effect of sub-band interference on imaging performance in synthetic aperture radar |
CN104076341A (zh) * | 2014-07-04 | 2014-10-01 | 西安电子科技大学 | 基于fpga和dsp的雷达回波实时模拟方法 |
CN104569938A (zh) * | 2015-02-09 | 2015-04-29 | 中国科学院电子学研究所 | 一种合成孔径雷达回波仿真器 |
CN104698445A (zh) * | 2015-03-30 | 2015-06-10 | 北京润科通用技术有限公司 | 一种雷达回波模拟方法及雷达回波信号模拟器 |
CN105954730A (zh) * | 2016-06-21 | 2016-09-21 | 电子科技大学 | 一种sar回波快速时域生成方法 |
CN105911532A (zh) * | 2016-06-29 | 2016-08-31 | 北京化工大学 | 基于深度协同的合成孔径雷达回波并行模拟方法 |
WO2018004503A1 (en) * | 2016-07-01 | 2018-01-04 | Havelsan Hava Elektronik Sanayi Ve Ticaret Anonim Sirketi | Radar video simulator |
CN107229051A (zh) * | 2017-05-26 | 2017-10-03 | 西安电子科技大学 | 基于gpu的视频sar回波仿真并行实现方法 |
CN107765226A (zh) * | 2017-09-18 | 2018-03-06 | 北京空间飞行器总体设计部 | 一种sar卫星雷达回波模拟方法、系统和介质 |
CN108519590A (zh) * | 2018-03-26 | 2018-09-11 | 北京理工大学 | 激光成像雷达目标回波信号模拟方法及模拟器 |
CN109901165A (zh) * | 2019-03-28 | 2019-06-18 | 河南九乾电子科技有限公司 | 星载sar回波的模拟装置及模拟方法 |
CN209821374U (zh) * | 2019-04-02 | 2019-12-20 | 南京雷电信息技术有限公司 | 一种雷达目标模拟器 |
CN110515050A (zh) * | 2019-08-30 | 2019-11-29 | 北京空间飞行器总体设计部 | 一种基于gpu的星载sar实时回波模拟器 |
CN112698280A (zh) * | 2020-12-09 | 2021-04-23 | 南京长峰航天电子科技有限公司 | 一种基于dsp和fpga架构的双基sar实时回波模拟方法 |
Non-Patent Citations (3)
Title |
---|
基于GPU的视频SAR回波仿真研究与实现;陈文俊;《中国优秀硕士学位论文全文数据库信息科技辑》;20190215(第2期);全文 * |
宽带高速目标回波模拟系统设计与实现;麻哲;《中国优秀硕士学位论文全文数据库信息科技辑》;20220115(第1期);全文 * |
雷达目标回波模拟器射频部件研究;曹庆刚;《中国优秀硕士学位论文全文数据库信息科技辑》;20040415(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114252855A (zh) | 2022-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114252855B (zh) | 一种同时sar/gmt回波模拟方法及其系统 | |
CN112415482B (zh) | 一种射频注入式数字信号合成仿真测试系统 | |
CN104614713B (zh) | 一种适合于艇载雷达系统的雷达回波信号模拟器 | |
CN112505643B (zh) | 一种雷达与红外复合导引头开环半实物仿真方法及系统 | |
CN109001697B (zh) | 一种多目标雷达回波模拟器 | |
CN104050716B (zh) | 一种海上多目标sar图像可视化建模方法 | |
CN109901165B (zh) | 星载sar回波的模拟装置及模拟方法 | |
CA1210117A (en) | Algorithm for radar coordinate conversion in digital scan converters | |
CN105334498B (zh) | 一种适用于多平台雷达的地物回波快速生成算法 | |
CN112286079B (zh) | 一种高拟真度无人机航电半实物实景仿真系统 | |
CN114442051B (zh) | 一种高逼真度弹载雷达回波模拟方法 | |
CN109085552A (zh) | 一种基于试飞数据的杂波半实物化仿真测试方法及系统 | |
CN111025246A (zh) | 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法 | |
CN108983240A (zh) | 基于正交调制体制的防撞毫米波雷达目标信号模拟系统及方法 | |
CN115128559B (zh) | 机载相控阵雷达实时信号级回波模拟方法及装置 | |
CN106980117A (zh) | 基于重叠孔径的毫米波雷达视频成像方法 | |
Lebarbu et al. | Complete radar simulation chain: Application to maritime patrol surveillance using SAR/ISAR modes | |
RU75058U1 (ru) | Цифровой имитатор бортовых радиолокационных систем | |
CN113391275A (zh) | 一种采用电子角反的大型复杂目标等效模拟方法和装置 | |
CN114791592B (zh) | 一种合成孔径雷达对抗模拟方法、装置、电子设备及介质 | |
CN112698280B (zh) | 一种基于dsp和fpga架构的双基sar实时回波模拟方法 | |
RU186130U1 (ru) | Многофункциональный имитатор радиолокационных целей | |
CN106154277A (zh) | 无人船监控区域相控阵超声波低频成像装置 | |
Karboski et al. | Airborne Polarimetric Doppler Phased Array Weather Radar: Digital Twin of the Active Electronically Scanned Array | |
CN101509973B (zh) | 基于孔径合成技术的三维高分辨目标散射系数测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 210000 Room 301, floor 3, building 75, zone B, entrepreneurship and innovation city, No. 15, Fengji Avenue, Yuhuatai District, Nanjing, Jiangsu Province Patentee after: Nanjing Thunderbolt Information Technology Co.,Ltd. Address before: 210000 15 Fengji Avenue, Yuhuatai District, Nanjing, Jiangsu Province Patentee before: NANJING LEADING INFORMATION TECHNOLOGY Co.,Ltd. |
|
CP03 | Change of name, title or address |