CN112731399A - 一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 - Google Patents
一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 Download PDFInfo
- Publication number
- CN112731399A CN112731399A CN202011534816.4A CN202011534816A CN112731399A CN 112731399 A CN112731399 A CN 112731399A CN 202011534816 A CN202011534816 A CN 202011534816A CN 112731399 A CN112731399 A CN 112731399A
- Authority
- CN
- China
- Prior art keywords
- matrix
- data matrix
- obtaining
- target
- sparse
- 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.)
- Granted
Links
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/91—Radar or analogous systems specially adapted for specific applications for traffic control
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明提供了一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,包括得到回波数据矩阵、得到地杂波对消后的数据矩阵、得到二维稀疏结构数据矩阵、得到二值数据矩阵、完成所述二值图聚类和得到架次信息;该方法基于二维稀疏结构模型,同时利用了飞机目标回波的多普勒信息和距离信息,大大降低了杂波和噪声的干扰,提高了架次识别的准确性;同时该方法可实现回波数据丢失情况下的目标架次识别。该方法弥补了窄带识别中对相参积累时间的要求和宽带识别中对信噪比要求较高的不足。
Description
技术领域
本发明涉及雷达架次识别领域,尤其涉及一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法。
背景技术
随着现代科技的快速发展,各类人造飞行器逐渐普及,特别地,得益于飞行器性能的提升和无线电网络的互联发展,无论是军事战场还是民事应用,空中编队目标越来越多。常规雷达作为能够提供目标位置和运动参数的传感器,已远远不能满足现代军民各领域的需求,要想获得编队目标的详细架次信息,需要基于现有信息资源提出准确的架次识别方法。雷达架次识别技术是以雷达作为传感器,对探测到的目标的数量进行自动识别,对空域监视、目标规模判定乃至国土防空预警具有十分重要的作用。空中编队目标架次识别主要是利用雷达回波信号提取编队目标特征从而对架次进行识别,按照雷达体制,编队目标架次识别技术可分为窄带识别方法和宽带识别方法。
窄带识别方法是根据编队目标之间的多普勒频率差异,通过时频分析类的方法实现目标架次的识别。这类方法对雷达信号的相干性和相参积累时间有着一定的要求,否则无法达到要求的多普勒频率分辨率;此外,该类方法无法应用于飞机相对雷达径向飞行的情况。宽带识别方法是利用目标的各散射中心在径向距离上的高分辨力,得到目标的一维高分辨距离像,通过对距离像中峰值数量的提取来实现目标架次的识别。但是仅仅依靠飞机目标的高分辨距离像也很难确定其架次,主要原因在于回波信号中含有噪声和杂波,特别是当目标距离雷达较远,信噪比较低时,难以确定检测阈值。
发明内容
本发明的目的在于针对上述现有技术的不足,提供了一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,基于二维稀疏结构模型,同时利用了飞机目标回波的多普勒信息和距离信息,大大降低了杂波和噪声的干扰,提高了架次识别的准确性。
为实现上述目的,本发明采用了如下技术方案:
本发明提供了一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,其特征在于,包括以下步骤:
S1、得到回波数据矩阵SM×N:对采集的雷达回波数据进行混频、相位检波和A/D采样,得到零中频的数字信号;
S2、得到地杂波对消后的数据矩阵SL×N:采用三脉冲MTI技术对所述矩阵SM×N的每一列数据做运算:si,j-2×si-1,j+si+1,j;
其中,i=2,3L M-1,j=1,L N,得到数据矩阵SL×N;其中,L=M-2;
S3、得到二维稀疏结构数据矩阵Psparse:利用同步正交匹配追踪算法,优化求解下式:
其中,ST为雷达回波所述数据矩阵SL×N的转置矩阵;
P为待求解的二维稀疏结构图矩阵;Pi,j表述矩阵第i行第j列的元素;Φ为稀疏观测矩阵;
所述稀疏观测矩阵的构造过程如下:
Φ=(F-1·H·F)-1
其中,F为尺寸为N×N的傅里叶矩阵;(·)-1矩阵的取逆运算;
其中,j为虚数单位;γ为发射信号的调频率;diag(·)表示将向量转换为对角矩阵;
其中,fs是采样频率,Δf=fs/N;
S4、得到二值数据矩阵P2,所述二值数据矩阵P2分两步获得:
单次回波检测数据矩阵Λ,其中:所述Psparse为所述二维稀疏结构数据矩阵;
其中,xi和yi是待检测单元左、右的参考单元;
S4.2、在得到所述单次回波检测数据矩阵Λ的基础上,先根据m次回波的二分层检测结果积累,再与给定的门限K进行比较以得到最终的二值数据矩阵P2;
S5、完成所述二值数据矩阵P2聚类:循环搜索,将所述二值图P2中相邻的1值点融合为一类,统计所述二值数据矩阵P2中类的数量;
S6、得到架次信息。
进一步,所述同步正交匹配追踪算法的求解步骤为:
进一步,所述S4.2包括以下步骤:
S4.2.3、确定所述Psparse中检测单元的开始点数START=M+1+N和检测单元的结束点数STOP=C-M–N;
S4.2.6、确定所述单次回波检测数据矩阵Λ的行数R和列数C;
进一步,所述S5的包括以下步骤:
S5.1、所述二值数据矩阵P2的元素个数为M×N,即矩阵P2有M行N列的元素。逐个像素搜索,相邻的值为1的元素归为同一类;
S5.2、记所述二值数据矩阵P2中一共只有Q类,每一类记为一个单独的集合,Cq(q=1,2,L,Q)。
进一步,所述S6包括以下步骤:
S6.1、统计所述S5中每一类所包含的像素个数,记为Yq(q=1,2,L,Q);
S6.4、统计所剩类的数量Qr,即为架次数。
本发明的有益效果为:该方法基于二维稀疏结构模型,同时利用了飞机目标回波的多普勒信息和距离信息,大大降低了杂波和噪声的干扰,提高了架次识别的准确性;同时该方法可实现回波数据丢失情况下的目标架次识别。该方法弥补了窄带识别中对相参积累时间的要求和宽带识别中对信噪比要求较高的不足。
1)利用了多测量矢量模型,对雷达在不同慢时间接收到的回波信号进行了建模,有助于提高信号恢复的精度。
2)在获得距离多普勒二值图前,通过对待检测单元左右参考单元数据的估计和合理选择的二次检测门限最大程度的保留了目标的各散射点,为后续聚类图像的获得奠定了基础。
3)通过对雷达回波信号的稀疏处理,可实现回波数据缺失或污染条件下的目标识别。
4)本算法同时利用了飞机目标回波的多普勒信息和距离信息弥补了单纯窄带识别中对相参积累时间的要求和单纯宽带识别中对信噪比要求较高的不足。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,其特征在于,包括以下步骤:
S1、得到回波数据矩阵SM×N:对采集的雷达回波数据进行混频、相位检波和A/D采样,得到零中频的数字信号;
S2、得到地杂波对消后的数据矩阵SL×N:采用三脉冲MTI技术对所述矩阵SM×N的每一列数据做运算:si,j-2×si-1,j+si+1,j;
其中,i=2,3L M-1,j=1,L N,得到数据矩阵SL×N;其中,L=M-2;
S3、得到二维稀疏结构数据矩阵Psparse:利用同步正交匹配追踪算法,优化求解下式:
其中,ST为雷达回波所述数据矩阵SL×N的转置矩阵;
P为待求解的二维稀疏结构图矩阵;Pi,j表述矩阵第i行第j列的元素;Φ为稀疏观测矩阵;
所述稀疏观测矩阵的构造过程如下:
Φ=(F-1·H·F)-1
其中,F为尺寸为N×N的傅里叶矩阵;(·)-1矩阵的取逆运算;
其中,j为虚数单位;γ为发射信号的调频率;diag(·)表示将向量转换为对角矩阵;
其中,fs是采样频率,Δf=fs/N;
所述同步正交匹配追踪算法的求解步骤为:
S4、得到二值数据矩阵P2,所述二值数据矩阵P2分两步获得:
单次回波检测数据矩阵Λ,其中:所述Psparse为所述二维稀疏结构数据矩阵;
其中,xi和yi是待检测单元左、右的参考单元;
S4.2、在得到所述单次回波检测数据矩阵Λ的基础上,先根据m次回波的二分层检测结果积累,再与给定的门限K进行比较以得到最终的二值数据矩阵P2;其具体步骤如下:
S4.2.3、确定所述Psparse中检测单元的开始点数START=M+1+N和检测单元的结束点数STOP=C-M–N;
S4.2.6、确定所述单次回波检测数据矩阵Λ的行数R和列数C;
S5、完成所述二值数据矩阵P2聚类:循环搜索,将所述二值数据矩阵P2中相邻的1值点融合为一类,统计所述二值数据矩阵P2中类的数量,其具体步骤如下:
S5.1、所述二值数据矩阵P2的元素个数为M×N,即矩阵P2有M行N列的元素。逐个像素搜索,相邻的值为1的元素归为同一类;
其中,相邻是指如下表的情况:
1 | 1 | 1 |
1 | 1 | 1 |
1 | 1 | 1 |
其中,对角线的中心处表示当前搜索的像素点,其周围的八个像素中,值为1则被归为当前搜索像素点的类。
S5.2、记所述二值数据矩阵P2中一共只有Q类,每一类记为一个单独的集合,Cq(q=1,2,L,Q)。
S6、得到架次信息,其具体步骤如下:
S6.1、统计所述S5中每一类所包含的像素个数,记为Yq(q=1,2,L,Q);
S6.4、统计所剩类的数量Qr,即为架次数。
以上所述实施例仅表达了本发明的实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (5)
1.一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,其特征在于,包括以下步骤:
S1、得到回波数据矩阵SM×N:对采集的雷达回波数据进行混频、相位检波和A/D采样,得到零中频的数字信号;
S2、得到地杂波对消后的数据矩阵SL×N:采用三脉冲MTI技术对所述矩阵SM×N的每一列数据做运算:si,j-2×si-1,j+si+1,j;
其中,i=2,3L M-1,j=1,L N,得到数据矩阵SL×N;其中,L=M-2;
S3、得到二维稀疏结构数据矩阵Psparse:利用同步正交匹配追踪算法,优化求解下式:
其中,ST为雷达回波所述数据矩阵SL×N的转置矩阵;
P为待求解的二维稀疏结构图矩阵;Pi,j表述矩阵第i行第j列的元素;Φ为稀疏观测矩阵;
所述稀疏观测矩阵的构造过程如下:
Φ=(F-1.H·F)-1
其中,F为尺寸为N×N的傅里叶矩阵;(·)-1矩阵的取逆运算;
H=diag(-jπfr 2/γ)
其中,j为虚数单位;γ为发射信号的调频率;diag(·)表示将向量转换为对角矩阵;
其中,fs是采样频率,Δf=fs/N;
S4、得到二值数据矩阵P2,所述二值数据矩阵P2分两步获得:
单次回波检测数据矩阵Λ,其中:所述Psparse为所述二维稀疏结构数据矩阵;
其中,xi和yi是待检测单元左、右的参考单元;
S4.2、在得到所述单次回波检测数据矩阵Λ的基础上,先根据m次回波的二分层检测结果积累,再与给定的门限K进行比较以得到最终的二值数据矩阵P2;
S5、完成所述二值数据矩阵P2聚类:循环搜索,将所述二值数据矩阵P2中相邻的1值点融合为一类,统计所述二值数据矩阵P2中类的数量;
S6、得到架次信息。
3.根据权利要求1所述的一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,其特征在于,所述S4.2包括以下步骤:
S4.2.3、确定所述Psparse中检测单元的开始点数START=M+1+N和检测单元的结束点数STOP=C-M–N;
S4.2.6、确定所述单次回波检测数据矩阵Λ的行数R和列数C;
4.根据权利要求1所述的一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法,其特征在于,所述S5的包括以下步骤:
S5.1、所述二值数据矩阵P2的元素个数为M×N,即矩阵P2有M行N列的元素。逐个像素搜索,相邻的值为1的元素归为同一类;
S5.2、记所述二值数据矩阵P2中一共只有Q类,每一类记为一个单独的集合,Cq(q=1,2,L,Q)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011534816.4A CN112731399B (zh) | 2020-12-22 | 2020-12-22 | 一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011534816.4A CN112731399B (zh) | 2020-12-22 | 2020-12-22 | 一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112731399A true CN112731399A (zh) | 2021-04-30 |
CN112731399B CN112731399B (zh) | 2023-06-02 |
Family
ID=75604304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011534816.4A Active CN112731399B (zh) | 2020-12-22 | 2020-12-22 | 一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112731399B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217676A (zh) * | 2013-05-06 | 2013-07-24 | 西安电子科技大学 | 基于双谱去噪的噪声背景下雷达目标识别方法 |
CN104122540A (zh) * | 2014-07-31 | 2014-10-29 | 电子科技大学 | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 |
CN104865568A (zh) * | 2015-06-02 | 2015-08-26 | 西安电子科技大学 | 基于稀疏重构的宽带雷达高速群目标分辨方法 |
CN105572649A (zh) * | 2015-12-11 | 2016-05-11 | 中北大学 | 基于稀疏傅里叶变换的雷达目标检测方法 |
CN107132535A (zh) * | 2017-04-07 | 2017-09-05 | 西安电子科技大学 | 基于变分贝叶斯学习算法的isar稀疏频带成像方法 |
CN108008385A (zh) * | 2017-11-20 | 2018-05-08 | 西安电子科技大学 | 基于稀疏贝叶斯学习的干扰环境isar高分辨成像方法 |
-
2020
- 2020-12-22 CN CN202011534816.4A patent/CN112731399B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217676A (zh) * | 2013-05-06 | 2013-07-24 | 西安电子科技大学 | 基于双谱去噪的噪声背景下雷达目标识别方法 |
CN104122540A (zh) * | 2014-07-31 | 2014-10-29 | 电子科技大学 | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 |
CN104865568A (zh) * | 2015-06-02 | 2015-08-26 | 西安电子科技大学 | 基于稀疏重构的宽带雷达高速群目标分辨方法 |
CN105572649A (zh) * | 2015-12-11 | 2016-05-11 | 中北大学 | 基于稀疏傅里叶变换的雷达目标检测方法 |
CN107132535A (zh) * | 2017-04-07 | 2017-09-05 | 西安电子科技大学 | 基于变分贝叶斯学习算法的isar稀疏频带成像方法 |
CN108008385A (zh) * | 2017-11-20 | 2018-05-08 | 西安电子科技大学 | 基于稀疏贝叶斯学习的干扰环境isar高分辨成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112731399B (zh) | 2023-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108802722B (zh) | 一种基于虚拟谱的弱目标检测前跟踪方法 | |
CN106226751B (zh) | 基于dp-tbd的机动目标检测与跟踪方法 | |
CN103439697B (zh) | 基于动态规划的目标检测方法 | |
CN113808174B (zh) | 基于全卷积网络和卡尔曼滤波的雷达小目标跟踪方法 | |
CN104714225B (zh) | 一种基于广义似然比的动态规划检测前跟踪方法 | |
CN107861107A (zh) | 一种适用于连续波雷达的双门限cfar与点迹凝聚方法 | |
CN107576959B (zh) | 一种基于区域映射解模糊的高重频雷达目标检测前跟踪方法 | |
CN110018461B (zh) | 基于高分辨距离像和单脉冲测角的群目标识别方法 | |
CN104459661B (zh) | 检测快速火炮类微弱目标的方法 | |
CN110058222B (zh) | 一种基于传感器选择的双层粒子滤波检测前跟踪方法 | |
CN110146873B (zh) | 分布式非相参雷达的目标位置和速度的估计方法 | |
CN105447867B (zh) | 基于isar图像的空间目标姿态估计方法 | |
CN109521420B (zh) | 基于多特征匹配的多目标跟踪方法 | |
CN107356922A (zh) | 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 | |
CN109100696A (zh) | 基于点状杂波图的慢速运动目标消除方法 | |
CN108414992B (zh) | 一种基于相位信息杂波图的目标检测方法 | |
CN105445710A (zh) | 隐形雷达低截获概率性能评估方法 | |
CN111208479A (zh) | 一种降低深度网络检测中虚警概率的方法 | |
CN113221631B (zh) | 基于卷积神经网络的序列脉冲抗干扰目标检测方法 | |
CN113960587A (zh) | 基于类别信息反馈的毫米波雷达多目标跟踪方法 | |
CN112731399A (zh) | 一种基于二维稀疏结构的宽带雷达飞机目标架次识别方法 | |
CN110221289A (zh) | 用于三坐标有源相控阵雷达的目标检测方法 | |
CN115657005A (zh) | 一种基于多普勒雷达多信息融合的无人机识别方法 | |
CN108983192B (zh) | 基于gps辐射源的雷达运动目标参数估计方法 | |
CN113009450A (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 |