CN112444810B - 一种雷达对空多目标超分辨方法 - Google Patents
一种雷达对空多目标超分辨方法 Download PDFInfo
- Publication number
- CN112444810B CN112444810B CN202011161591.2A CN202011161591A CN112444810B CN 112444810 B CN112444810 B CN 112444810B CN 202011161591 A CN202011161591 A CN 202011161591A CN 112444810 B CN112444810 B CN 112444810B
- Authority
- CN
- China
- Prior art keywords
- target
- resolution
- grower
- radar
- virtual array
- 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
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9094—Theoretical aspects
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/426—Scanning radar, e.g. 3D radar
-
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种雷达对空多目标超分辨方法,应用于雷达探测与成像技术领域,针对传统实孔径雷达对空多目标探测时存在的角分辨率低的问题;本发明首先通过扫描雷达方位序列回波的虚拟阵列信号表征,将扫描雷达波束域的角超分辨问题转化为虚拟阵元域的到达角估计问题;然后针对虚拟阵列到达角估计的非凸优化求解问题,采用一种基于虚拟阵列孔径生长的多维迭代寻优超分辨方法进行求解,提升了雷达对空多目标角度分辨能力;相比传统技术,本发明方法具有更高的角分辨率,可用于空中编队目标的超分辨率探测。
Description
技术领域
本发明属于雷达探测与成像技术领域,特别涉及一种雷达对空多目标探测与分辨技术。
背景技术
雷达对空目标分辨,对于战略预警、目标探测与追踪、飞行管制等领域具有重要意义。然而,现有的实孔径雷达探测方法分辨率受孔径限制,无法实现远距离条件下多个目标的分辨;逆合成空孔径雷达(ISAR)成像方法虽然可以获得较高的分辨率,但是对成像构型要求较高,且信号处理过程复杂;动目标检测(MTD,Microwave Traffic Detector)技术受到机理限制,无法分辨具有相同运动速度的多个目标。
针对雷达对空多目标分辨难题,文献:Blair W D,Brandt-Pearce M.MonopulseDOA estimation of two unresolved Rayleigh targets[J].Aerospace and ElectronicSystems,IEEE Transactions on,2001,37(2):452-469.,提出一种采用单脉冲技术方法,主要适用于单个强点目标的定位,但对于存在多散射中心的复杂目标环境下,将会出现严重的角闪烁现象;文献Guan J,Huang Y,Yang J,et al.Improving angular resolutionbased on maximum a posteriori criterion for scanning radar[C]//RadarConference(RADAR),2012IEEE.IEEE,2012: 0451-0454.,该文章提出一种贝叶斯框架下的最大后验解卷积方法。虽然该方法能够实现多个目标的分辨,但依赖于很高的信噪比,实际应用局限性较大。
发明内容
为解决传统雷达对空多目标探测时存在的角分辨率低的问题,本发明提出一种雷达对空多目标超分辨方法,相比传统技术,本发明方法具有更高的角分辨率,可用于空中编队目标的超分辨率探测。
本发明采用的技术方案为:一种雷达对空多目标超分辨方法,将扫描雷达波束域的角超分辨问题转化为虚拟阵元域的到达角估计问题;然后采用基于虚拟阵列孔径生长的多维迭代寻优超分辨方法进行求解,得到目标的散射系数分布的估计值。
将扫描雷达波束域的角超分辨问题转化为虚拟阵元域的到达角估计问题,具体为:
A1、根据同一距离单元扫描雷达方位向回波,计算目标散射系数的空间频谱;
A2、根据目标散射系数的空间频谱,得到目标散射系数的空间时域和空间频域的映射关系;
A3、采用虚拟阵列表征目标散射系数的空间时域和空间频域的映射关系。
虚拟阵列孔径生长的过程为:
B1、初始时,设置阵列孔径一端的生长长度、生长器的阶数,并定义起始时刻的前向和后向生长误差;
B2、计算一阶生长器的权系数,并定义反射系数;
B3、计算多阶生长器的权系数与反射系数;
B4、最后各阶生长器的权系数与反射系数,得到虚拟孔径生长后的目标散射系数空间频谱数据。
所述多维迭代过程为:
C1、对虚拟阵列数据进行预处理,计算扫描区域中Q个点的目标散射系数的初始估计值;
C2、对Q个点的目标散射系数的初始估计值进行对角排列形成对角矩阵;
C3、根据步骤C2的对角矩阵定义自相关矩阵;
C4、根据步骤C3的自相关矩阵得到代价函数;
C5、依次遍历扫描区域中的点,根据每次迭代得到的目标散射系数更新自相关矩阵的逆矩阵;
C6、重复步骤C5,直至迭代结束,得到目标的散射系数分布的估计值。
本发明的有益效果:本发明首先通过扫描雷达方位序列回波的虚拟阵列信号表征,将扫描雷达波束域的角超分辨问题转化为数学上可解的虚拟阵元域的到达角估计问题,从机理上突破雷达对空目标分辨瓶颈;然后针对虚拟阵列到达角估计的非凸优化求解难题,采用一种多维迭代超分辨方法进行求解。相比传统技术,本发明方法具有更高的角分辨率,大大推远了空中目标的可分辨探测距离。
附图说明
图1为本发明的方法流程图;
图2为扫描雷达回波模型示意图;
图3为本发明实施例提供的天线方向图;
图4为本发明实施例提供的实波束扫描的分辨结果;
图5为本发明实施例提供的本发明方法的超分辨结果。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
本发明的扫描雷达对空目标超分辨成像方法流程图如图1所示,包括以下步骤:
步骤一:扫描雷达方位序贯回波的虚拟阵列信号表征
如图2所示,雷达发射线性调频信号以提高距离向分辨率。为获得二维微波图像,以雷达位置为圆心,在发射脉冲信号的同时以扫描速度w在方位向上进行波束扫描。在相邻发射脉冲间隔时间段内,该雷达接收回波信号。针对同一距离单元,扫描雷达方位向回波可以表示为:
sm=h*σm+nm (1)
其中和分别表示第m个距离单元处的方位回波数据、目标散射系数和加性噪声。h=[h(θ1),h(θ2),…,h(θL)]T表示天线方向图的采样,符号*表示卷积关系,Q、K分别表示方位向和目标区域的采样点数。从公式(1)计算得出目标散射系数的空间频谱可表示为:
利用公式(2)的目标散射系数的空间频谱,目标散射系数的空间时域和空间频域的映射关系如下:
其中Ω代表波束扫描区域的范围,θi(i=1,2,…Q)代表扫描范围内的Q个采样角。
根据傅里叶转向矩阵F的特性,可以用虚拟阵列模型来描述公式(3)中的模型。假设均匀线阵的天线数目为B′,相邻阵元之间的间距为d。对于位于θi方向的目标,第b根天线与参考阵元之间的相位差为:
其中λ是发射信号的波长,比较公式(4)和公式(5)中傅里叶转向矩阵的相位分布,得出扫描雷达方位序贯回波的虚拟阵列可表征为:
该虚拟阵列表征了实孔径扫描雷达与等效虚拟阵列之间的映射关系。虚拟阵列模型可以表示为:
步骤二:虚拟阵列孔径生长
其中:阵列孔径一端的生长长度为q,生长器的阶数为p。扫描雷达虚拟阵列孔径生长的过程如下:
(1)令m=0时,定义起始时刻的前向和后向生长误差:
计算:
(3)当m≥2时,计算m阶生长器的权系数和反射系数:
步骤三:虚拟阵列到达角估计问题的非凸优化求解
多维迭代超分辨方法的过程如下:
1)计算扫描区域中Q个点的散射系数初始估计值:
3)迭代计算
对扫描区域所有的点进行遍历(k=1,…,Q),每次遍历时做如下操作:
②对上次迭代结果进行修正,并实现计算的非负性:
③更新自相关矩阵的逆矩阵:
以下结合具体数据对本发明的实现过程进行进一步地阐述:
步骤一:扫描雷达方位序贯回波的虚拟阵列信号表征
表1列出了扫描雷达平台的仿真参数,例如载波频率fc、发射信号的带宽和时宽、距离向和方位向的采样频率等等,采样速率满足奈奎斯特采样定律。
表1扫描雷达系统参数
假设在扫描区域中,每个方位采样点上都有目标存在,令这些目标的位置参数为θ=(θ1,θ2,...θK),幅度参数为σ=(σ1,σ2,...,σK),则这些目标回波信号经相干解调后可表示为:
其中:σk(k=1,2,…,K)表示位置在θk处的目标的散射系数,t表示距离快时间,变化范围由发射机到目标的双程距离决定,τ为方位慢时间,rect(·)为矩形窗函数,exp(·)为指数函数,a(θk,τ)表示τ时刻指向角度为θk的天线方向图函数值,K表示目标的个数;Kr表示发射信号的调频斜率;Rk表示位置在θk处的目标相对于雷达的距离(这里把公式的R0改成Rk);λ表示发射信号的波长;c代表光速,数值为3×108m/s。
对于接收到的回波信号进行预处理,即进行距离向脉冲压缩。首先对回波进行距离向 FFT得到公式:
其中,B表示发射信号的带宽。
然后通过距离向乘以匹配滤波器频域匹配函数H(fr),得到距离压缩频域数据 SS(t,fr)。
Rk(t)表示位置在θk处的目标相对于雷达的距离函数;
最后对公式(23)的距离压缩频域数据进行IFFT变换,得到距离时域数据:
对于各个距离单元,方位扫描成像的回波模型及处理方式是一致的,因此任意选取一个距离单元的回波数据S。这里设置的信噪比为20dB。根据公式(6)对应的虚拟阵列映射关系,构造阵列信号的导向矩阵A。
步骤二:虚拟阵列孔径生长
选取一个距离单元的回波数据S,根据公式(2)计算出目标散射系数的估计值,令表示目标散射系数的空间频谱的估计值;表示维度为B×1的复数空间,设置阵列孔径一端的生长长度q,生长器的阶数为p。前向生长误差和后向生长误差可表示为:
对于公式(25),设计如下所示的格型滤波器结构:
步骤三:虚拟阵列到达角估计问题的非凸优化求解
经过步骤二的孔径生长过程,虚拟阵元的数量得到增加。各距离单元回波此时可以统一采用虚拟阵列模型进行表征:
fp=ln|R|+S(n)HR-1S(n) (32)
图3为天线方向图,图4为传统实波束扫描分辨的结果,可以看出分辨不出两个预先设置好的点目标;图5为采用本发明方法的分辨结果,可以看出很好的分辨出了设置的两个点目标;分辨结果如表2所示,分辨出的两个点目标位置是-1度和1度,和表2设置的吻合,并且散射系数也是吻合的:
表5具体实施方式采用的目标场景参数
目标编号 | 方位向角度 | 散射系数 |
1 | -1° | 0.5 |
2 | 1° | 1 |
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (2)
1.一种雷达对空多目标超分辨方法,其特征在于,将扫描雷达波束域的角超分辨问题转化为虚拟阵元域的到达角估计问题;然后采用基于虚拟阵列孔径生长的多维迭代寻优超分辨方法进行求解,得到目标的散射系数分布的估计值;
将扫描雷达波束域的角超分辨问题转化为虚拟阵元域的到达角估计问题,具体为:
A1、根据同一距离单元扫描雷达方位向回波,计算目标散射系数的空间频谱;
A2、根据目标散射系数的空间频谱,得到目标散射系数的空间时域和空间频域的映射关系;
A3、采用虚拟阵列表征目标散射系数的空间时域和空间频域的映射关系;
虚拟阵列孔径生长的过程为:
B1、初始时,设置阵列孔径一端的生长长度、生长器的阶数,并定义起始时刻的前向和后向生长误差;
B2、计算一阶生长器的权系数,并定义反射系数;所述生长器的权系数计算式为:
B3、计算多阶生长器的权系数与反射系数;
B4、最后各阶生长器的权系数与反射系数,得到虚拟孔径生长后的目标散射系数空间频谱数据;步骤B4计算得到虚拟孔径生长后的目标散射系数空间频谱数据,表达式为:
2.根据权利要求1所述的一种雷达对空多目标超分辨方法,其特征在于,所述多维迭代过程为:
C1、对虚拟阵列数据进行预处理,计算扫描区域中Q个点的目标散射系数的初始估计值;
C2、对Q个点的目标散射系数的初始估计值进行对角排列形成对角矩阵;
C3、根据步骤C2的对角矩阵定义自相关矩阵;
C4、根据步骤C3的自相关矩阵得到代价函数;
C5、依次遍历扫描区域中的点,根据每次迭代得到的目标散射系数更新自相关矩阵的逆矩阵;
C6、重复步骤C5,直至迭代结束,得到目标的散射系数分布的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011161591.2A CN112444810B (zh) | 2020-10-27 | 2020-10-27 | 一种雷达对空多目标超分辨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011161591.2A CN112444810B (zh) | 2020-10-27 | 2020-10-27 | 一种雷达对空多目标超分辨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112444810A CN112444810A (zh) | 2021-03-05 |
CN112444810B true CN112444810B (zh) | 2022-07-01 |
Family
ID=74736194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011161591.2A Active CN112444810B (zh) | 2020-10-27 | 2020-10-27 | 一种雷达对空多目标超分辨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112444810B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640762B (zh) * | 2021-07-23 | 2024-02-09 | 北京理工大学 | 基于多维参数空间自聚焦的雷达目标方位超分辨估计方法 |
CN113960530A (zh) * | 2021-10-11 | 2022-01-21 | 中国科学院声学研究所 | 一种基于直达波和海面反射波到达角的声源被动定位方法 |
CN114384516A (zh) * | 2022-01-12 | 2022-04-22 | 电子科技大学 | 一种基于重建前检测的实孔径雷达实时角超分辨方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323583A (zh) * | 2011-09-13 | 2012-01-18 | 电子科技大学 | 一种超分辨线阵三维合成孔径雷达成像方法 |
CN111081266A (zh) * | 2019-12-18 | 2020-04-28 | 暗物智能科技(广州)有限公司 | 一种训练生成对抗网络、语音增强方法及系统 |
CN111781575A (zh) * | 2020-07-13 | 2020-10-16 | 深圳大学 | 一比特信号单快拍波达方向估计方法及相关组件 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8184043B2 (en) * | 2010-03-12 | 2012-05-22 | The Boeing Company | Super-resolution imaging radar |
CN103955701B (zh) * | 2014-04-15 | 2017-02-22 | 浙江工业大学 | 多层次结合的多视合成孔径雷达图像目标识别方法 |
CN104020469B (zh) * | 2014-05-30 | 2016-08-17 | 哈尔滨工程大学 | 一种mimo雷达距离-角度二维超分辨率成像算法 |
CN104698430B (zh) * | 2015-03-27 | 2017-06-06 | 北京邮电大学 | 一种用于提高精度的基于虚拟天线阵列的角度估计方法 |
CN105260222B (zh) * | 2015-10-13 | 2018-10-26 | 哈尔滨工程大学 | 一种可重构编译器中循环流水迭代间启动间距优化方法 |
WO2017187243A1 (en) * | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | Vehicular radar sensing system utilizing high rate true random number generator |
CN106405487B (zh) * | 2016-10-17 | 2018-08-21 | 重庆邮电大学 | 一种通用基于扩展esprit技术的空间谱估计方法 |
CN108933745B (zh) * | 2018-07-16 | 2020-07-10 | 北京理工大学 | 一种基于超分辨率角度和时延估计的宽带信道估计方法 |
-
2020
- 2020-10-27 CN CN202011161591.2A patent/CN112444810B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323583A (zh) * | 2011-09-13 | 2012-01-18 | 电子科技大学 | 一种超分辨线阵三维合成孔径雷达成像方法 |
CN111081266A (zh) * | 2019-12-18 | 2020-04-28 | 暗物智能科技(广州)有限公司 | 一种训练生成对抗网络、语音增强方法及系统 |
CN111781575A (zh) * | 2020-07-13 | 2020-10-16 | 深圳大学 | 一比特信号单快拍波达方向估计方法及相关组件 |
Also Published As
Publication number | Publication date |
---|---|
CN112444810A (zh) | 2021-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112444810B (zh) | 一种雷达对空多目标超分辨方法 | |
CN110275166B (zh) | 基于admm的快速稀疏孔径isar自聚焦与成像方法 | |
CN109212526B (zh) | 用于高频地波雷达的分布式阵列目标角度测量方法 | |
CN105445701B (zh) | Ddma‑mimo雷达目标的单脉冲角度估计方法 | |
CN109613532B (zh) | 一种机载雷达实时多普勒波束锐化超分辨成像方法 | |
CN106353744A (zh) | 基于双基地fda‑mimo雷达的多参数联合估计方法 | |
CN106772253B (zh) | 一种非均匀杂波环境下的雷达杂波抑制方法 | |
CN108693511B (zh) | 时分复用mimo雷达的运动目标角度计算方法 | |
CN101881822B (zh) | 一种针对共享谱雷达同频干扰的抑制方法 | |
CN104502898B (zh) | 将修正rft和修正mdcft相结合的机动目标参数估计方法 | |
CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN110174650B (zh) | 基于两维联合插值的气象雷达风电场杂波抑制方法 | |
CN108226928B (zh) | 基于期望传播算法的逆合成孔径雷达成像方法 | |
CN113777599B (zh) | 一种基于fda-mimo雷达的速度解模糊方法 | |
CN109471083A (zh) | 基于空时级联的机载外辐射源雷达杂波抑制方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN111060879A (zh) | 一种基于二维匹配滤波结果的联合副瓣抑制方法 | |
CN109613507A (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
CN111060909B (zh) | 一种机载雷达斜前视超分辨成像方法 | |
CN115291207A (zh) | 基于mimo雷达的小型旋翼无人机多目标检测方法 | |
CN111965612A (zh) | 一种基于子空间投影的杂波抑制方法 | |
CN113376607B (zh) | 机载分布式雷达小样本空时自适应处理方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN110579737A (zh) | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 | |
CN113900099A (zh) | 稀疏孔径isar机动目标成像与定标方法 |
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 |