CN109581366A - 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 - Google Patents
一种基于目标导向矢量失配的离散旁瓣杂波识别方法 Download PDFInfo
- Publication number
- CN109581366A CN109581366A CN201811557571.XA CN201811557571A CN109581366A CN 109581366 A CN109581366 A CN 109581366A CN 201811557571 A CN201811557571 A CN 201811557571A CN 109581366 A CN109581366 A CN 109581366A
- Authority
- CN
- China
- Prior art keywords
- discrete
- target
- mismatch
- clutter recognition
- steering vector
- 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/886—Radar or analogous systems specially adapted for specific applications for alarm systems
-
- 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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- 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/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于目标导向矢量失配的离散旁瓣杂波识别方法,属于雷达滤波技术领域。本发明通过分析离散旁瓣杂波和真实目标的特点,采用与目标空域导向矢量失配的导向矢量进行杂波抑制,易知此时会造成目标信号的丢失,而对离散旁瓣杂波影响不大,以此为依据识别并剔除离散旁瓣杂波,该方法具有良好的有效性和实用性。本发明方法建立在前期STAP杂波抑制基础上,通过改变目标导向约束,比较不同权系数对应的杂波抑制结果,从而有效剔除离散旁瓣杂波点。算法实现简单,计算复杂度低;本发明方法对于从旁瓣进入的任何类型的离散旁瓣杂波均可有效识别剔除,因此适用范围广,稳健性高。
Description
技术领域
本发明涉及一种基于目标导向矢量失配的离散旁瓣杂波识别方法,属于雷达滤波技术领域。
背景技术
目前尚无在轨运行的天基预警雷达,公开文献资料中的天基预警雷达相关资料主要集中在天基预警雷达概念体制研究、杂波分析及抑制等几个方面,根据调研尚无天基预警雷达离散旁瓣杂波识别剔除的相关文献资料。对于机载雷达***,目前的离散杂波剔除方法主要包括两种:保护通道法和检波后时间灵敏度控制法。
保护通道的工作原理是通过比较两个并行接收通道的输出,一个与主天线相连接,另外一个与保护天线连接。其中要求保护天线的主瓣波束为宽波束,要求其主瓣方向图高于主天线的旁瓣。分别对两幅天线接收到的回波进行恒虚警检测,对二者中检测到的目标进行比对,如果保护通道的回波强则为离散杂波,反之为真实目标。该方法在机载预警中具有良好的试用价值,但是对于天基雷达并不适用,主要原因在于天基雷达波束完全指向地球,并且距离环也远大于机载雷达,因此其杂波能量远高于机载雷达,在进行杂波抑制前难以进行离散旁瓣杂波是否剔除的判断,甚至会经常造成目标信号被误认杂波而被剔除的情况(保护通道内的发射旁瓣杂波会高于杂波抑制后的目标信号)。
检波后STC处理的工作过程为:首先对所有检测到的目标进行距离解模糊处理,然后补偿距离引起的回波强度变化,再次设计检测门限,小于门限的为离散旁瓣杂波,反之为目标。该方法主要适用于特殊的机载雷达应用背景,其认为相同距离下,雷达的低旁瓣设计足以抵消目标RCS比离散旁瓣杂波弱的劣势,因此只要将距离引起的回波幅度变化补偿掉即可,该方法在机载预警雷达具有有效性。但是对于天基雷达并不适用,主要原因在于天基雷达的距离变化与机载雷达相比并不大(以轨道高度1000公里卫星为例,其星下点距离为1000公里,入射角为80度时,斜距不到3000公里,相差只有三倍;而机载雷达从十几公里至几百公里不等,相差数十倍)
总之,目前针对空中高速目标探测的天基预警雷达离散旁瓣杂波识别方法尚无相关公开文献。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种基于目标导向矢量失配的离散旁瓣杂波识别方法,通过分析离散旁瓣杂波和真实目标的特点,利用导向矢量约束对来自旁瓣的强杂波点进行抑制。
本发明的技术解决方案是:一种基于目标导向矢量失配的离散旁瓣杂波识别方法,包括如下步骤:
接收雷达天线阵的回波数据,并由目标未失配时的空时导向矢量对所述回波数据进行第一次杂波抑制;所述天线阵包括N个天线子阵;
对第一次杂波抑制后的回波数据进行恒虚警检测,将检测到的目标的距离和速度的集合作为第一集合;
由目标失配时的空域导向矢量计算目标失配时的空时导向矢量,并由所述目标失配时的空时导向矢量对原始接收的回波数据进行第二次杂波抑制;
对第二次杂波抑制后的回波数据进行恒虚警检测,将检测到的目标的距离和速度的集合作为第二集合;
逐一比对所述第一集合和第二集合中的目标的距离和速度,判定两个集合中距离和速度均相同的目标为离散旁瓣杂波。
进一步地,所述目标失配时的空时导向矢量为其中,St为经过降维处理后的时域导向矢量,Ss'为目标失配的空域导向矢量。
进一步地,所述第二次杂波抑制的方法为采用第二滤波器对原始接收的回波数据进行滤波;所述第二滤波器的权系数为其中,为估计的杂波协方差矩阵。
进一步地,所述目标失配时的空域导向矢量为其中,2≤k<N,k为整数,N为天线子阵的个数,d为每个天线子阵的长度,θ为天线阵3dB主瓣波束宽度,λ为雷达天线阵发射信号的波长。
进一步地,所述其中,L为所述天线阵的总长度。
进一步地,所述d=L/N;其中,L为所述天线阵的总长度。
进一步地,所述第一次杂波抑制的方法为采用第一滤波器对原始接收的回波数据进行滤波;所述第一滤波器的权系数为其中,为估计的杂波协方差矩阵,S为目标的空时导向矢量。
进一步地,所述目标未失配时的空时导向矢量为其中,St为经过降维处理后的时域导向矢量,Ss为目标空域导向矢量。
进一步地,所述目标空域导向矢量Ss=[1,1,…,1]T。
进一步地,所述经过降维处理后的时域导向矢量为St(fd0)=[g1,g2,…,gM]T;其中,fd_j为第j个多普勒通道,1≤j≤200,gn为时域导向矢量系数,1≤n≤M,M为降维后的时域导向矢量维数。
本发明与现有技术相比的优点在于:
(1)本发明方法建立在前期STAP杂波抑制基础上,通过改变目标导向约束,比较不同权系数对应的杂波抑制结果,从而有效剔除离散旁瓣杂波点。算法实现简单,计算复杂度低;
(2)本发明方法对于从旁瓣进入的任何类型的离散旁瓣杂波均可有效识别剔除,因此适用范围广,稳健性高。
附图说明
图1为本发明方法实施流程图;
图2为第一次杂波抑制后检测到的4个目标;
图3为目标导向失配后检测到的3个目标。
具体实施方式
本发明包括如下步骤:
步骤1假设雷达整个天线方位阵面总长度为L米,对应天线3分贝主瓣波束宽度其中,λ为雷达发射信号的波长,天线方位阵面安装方向与卫星平台速度方向相同;方位阵面共分为N个天线子阵,因此每个天线子阵长度为d=L/N米;雷达采用全天线发射,N个子阵接收;脉冲重复频率表示为PRF。
步骤2***采用空时自适应处理技术对雷达接收N个子阵的回波数据进行杂波抑制,假设目标的空时导向矢量为:其中Ss为N×1维的目标空域导向矢量,其与发射波束指向对应的导向矢量相同,在接收数据端Ss=[1,1,…,1]T;St假设为经过降维处理后的时域导向矢量,表示为St(fd0)=[g1,g2,…,gM]T(降维处理相关内容参见L.E.Brennan,D.J.Piwinski,and F.M.Staudaher.Comparison of space-time adaptive processingapproaches using experimental airborne radar data[C].The Record of NationalRadar Conference,IEEE,1993:176-181),M为降维后的时域导向矢量维数;为估计的杂波协方差矩阵,协方差矩阵估计方法属于STAP处理的一项基本原理。因此对应的空时自适应滤波器权系数为以此对雷达回波数据进行杂波抑制处理。
步骤3对杂波抑制后的雷达数据进行恒虚警检测,假设检测到P个候选目标,第i个目标对应的距离和速度为(Ri,Vi),i=1,2,…,P,其中检测到的候选目标中既包含真实目标,也包含离散旁瓣杂波。
步骤4计算目标失配的空域导向矢量其中k满足2≤k<N;从而得到新的失配空时导向矢量进一步计算目标空时导向失配的滤波器权系数目标失配是指目标空域导向约束偏差与真实的雷达发射波束指向失配。
步骤5采用步骤4中的对雷达回波数据重新进行杂波抑制处理,然后进行恒虚警检测,由于此时目标空域导向约束偏差与真实的雷达发射波束指向完全失配,因此真实目标信号被抑制,剩余的一定为离散旁瓣杂波。假设检测到Q个离散旁瓣杂波(Q必定小于等于P),第i(i=1,2,…,Q)个离散旁瓣杂波对应的距离和速度为(Ri,Vi)。
步骤6根据距离速度信息,将步骤3中目标与步骤5中检测的虚假目标逐一比对,如果存在距离速度相同的情况则予以剔除,剩余的即为真实目标。
下面对本发明实施例及效果作进一步的详细描述,具体实施流程图如图1所示。
本发明的使用场景为:假设雷达整个天线方位阵面总长度为40米,雷达发射信号的波长λ为0.2m,方位阵面共分为10个天线子阵,子阵长度d=4m,天线方位阵面安装方向与卫星平台速度方向相同,单个相参处理时间间隔(CPI)内脉冲数为200,雷达回波中包含3个离散旁瓣杂波和1个目标信号
离散旁瓣杂波识别具体实现步骤如下:
步骤1根据天基雷达***参数,天线3分贝主瓣波束宽度θ为0.28°。
步骤2***采用空时自适应处理技术对雷达接收N个子阵的回波数据进行杂波抑制,假设目标的空时导向矢量为:其中Ss为10×1维的目标空域导向矢量,其与发射波束指向对应的导向矢量相同,在接收数据端Ss=[1,1,…,1]T;降维处理后的时域导向矢量表示为St(fd_j),fd_j为第j(1≤j≤200)个多普勒通道。因此可以获得空时自适应滤波器权系数为以此对雷达回波数据进行杂波抑制处理。
步骤3对杂波抑制后的雷达数据进行恒虚警检测,假设对于检测到的某个候选目标,记录其对应的距离和速度为(Ri,Vi),如图2所示此时检测到4个可能目标,4个目标距离速度分别为(1200km,70m/s),(1210km,50m/s),(1230km,70m/s),(1235km,-80m/s)。
步骤4令目标失配的空域导向矢量为从而得到新的失配空时导向矢量进一步计算目标空时导向失配的滤波器权系数
步骤5采用步骤4中的对雷达回波数据重新进行杂波抑制处理,然后进行恒虚警检测,由于此时目标空域导向约束偏差与真实的雷达发射波束指向完全失配,因此真实目标信号被抑制,剩余的一定为离散旁瓣杂波,如图3所示此时检测到3个目标,对应的距离和速度分别为:(1200km,70m/s),(1230km,70m/s),(1235km,-80m/s)。
步骤6根据距离速度信息,将步骤3中目标与步骤5中检测的虚假目标逐一比对,如果存在距离速度相同的情况则予以剔除,剩余的即为真实目标,通过对比可知(1210km,50m/s)为真实目标,离散旁瓣杂波被剔除掉。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (10)
1.一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于,包括如下步骤:
接收雷达天线阵的回波数据,并由目标未失配时的空时导向矢量对所述回波数据进行第一次杂波抑制;所述天线阵包括N个天线子阵;
对第一次杂波抑制后的回波数据进行恒虚警检测,将检测到的目标的距离和速度的集合作为第一集合;
由目标失配时的空域导向矢量计算目标失配时的空时导向矢量,并由所述目标失配时的空时导向矢量对原始接收的回波数据进行第二次杂波抑制;
对第二次杂波抑制后的回波数据进行恒虚警检测,将检测到的目标的距离和速度的集合作为第二集合;
逐一比对所述第一集合和第二集合中的目标的距离和速度,判定两个集合中距离和速度均相同的目标为离散旁瓣杂波。
2.根据权利要求1所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述目标失配时的空时导向矢量为其中,St为经过降维处理后的时域导向矢量,S′s为目标失配的空域导向矢量。
3.根据权利要求2所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述第二次杂波抑制的方法为采用第二滤波器对原始接收的回波数据进行滤波;所述第二滤波器的权系数为其中,为估计的杂波协方差矩阵。
4.根据权利要求2所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述目标失配时的空域导向矢量为其中,2≤k<N,k为整数,N为天线子阵的个数,d为每个天线子阵的长度,θ为天线阵3dB主瓣波束宽度,λ为雷达天线阵发射信号的波长。
5.根据权利要求4所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述其中,L为所述天线阵的总长度。
6.根据权利要求4所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述d=L/N;其中,L为所述天线阵的总长度。
7.根据权利要求1所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述第一次杂波抑制的方法为采用第一滤波器对原始接收的回波数据进行滤波;所述第一滤波器的权系数为其中,为估计的杂波协方差矩阵,S为目标的空时导向矢量。
8.根据权利要求7所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述目标未失配时的空时导向矢量为其中,St为经过降维处理后的时域导向矢量,Ss为目标空域导向矢量。
9.根据权利要求8所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述目标空域导向矢量Ss=[1,1,…,1]T。
10.根据权利要求2或8所述的一种基于目标导向矢量失配的离散旁瓣杂波识别方法,其特征在于:所述经过降维处理后的时域导向矢量为St(fd0)=[g1,g2,…,gM]T;其中,fd_j为第j个多普勒通道,1≤j≤200,gn为时域导向矢量系数,1≤n≤M,M为降维后的时域导向矢量维数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811557571.XA CN109581366B (zh) | 2018-12-19 | 2018-12-19 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811557571.XA CN109581366B (zh) | 2018-12-19 | 2018-12-19 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109581366A true CN109581366A (zh) | 2019-04-05 |
CN109581366B CN109581366B (zh) | 2021-02-09 |
Family
ID=65931057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811557571.XA Active CN109581366B (zh) | 2018-12-19 | 2018-12-19 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109581366B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113567978A (zh) * | 2021-07-29 | 2021-10-29 | 电子科技大学 | 一种多基分布式雷达协同成像方法 |
CN116990773A (zh) * | 2023-09-27 | 2023-11-03 | 广州辰创科技发展有限公司 | 基于自适应门限的低慢小目标检测方法、装置和存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5499030A (en) * | 1994-03-18 | 1996-03-12 | The United States Of America As Represented By The Secretary Of The Air Force | Expert system constant false alarm rate (CFAR) processor |
CN101556328A (zh) * | 2009-05-08 | 2009-10-14 | 西安电子科技大学 | 基于杂波协方差矩阵的机载雷达空时二维滤波器构建方法 |
CN101710175A (zh) * | 2009-11-13 | 2010-05-19 | 北方工业大学 | 一种雷达信号处理方法 |
US8897571B1 (en) * | 2011-03-31 | 2014-11-25 | Raytheon Company | Detection of targets from hyperspectral imagery |
CN104931935A (zh) * | 2015-06-15 | 2015-09-23 | 电子科技大学 | 基于旁瓣对消的稳健波束形成方法 |
CN107167803A (zh) * | 2017-05-25 | 2017-09-15 | 河海大学 | 基于导向矢量失配估计的鲁棒波束域自适应波束形成方法 |
-
2018
- 2018-12-19 CN CN201811557571.XA patent/CN109581366B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5499030A (en) * | 1994-03-18 | 1996-03-12 | The United States Of America As Represented By The Secretary Of The Air Force | Expert system constant false alarm rate (CFAR) processor |
CN101556328A (zh) * | 2009-05-08 | 2009-10-14 | 西安电子科技大学 | 基于杂波协方差矩阵的机载雷达空时二维滤波器构建方法 |
CN101710175A (zh) * | 2009-11-13 | 2010-05-19 | 北方工业大学 | 一种雷达信号处理方法 |
US8897571B1 (en) * | 2011-03-31 | 2014-11-25 | Raytheon Company | Detection of targets from hyperspectral imagery |
CN104931935A (zh) * | 2015-06-15 | 2015-09-23 | 电子科技大学 | 基于旁瓣对消的稳健波束形成方法 |
CN107167803A (zh) * | 2017-05-25 | 2017-09-15 | 河海大学 | 基于导向矢量失配估计的鲁棒波束域自适应波束形成方法 |
Non-Patent Citations (2)
Title |
---|
YUZE SUN等: ""Robust sparse Bayesian learning STAP method for discrete interference suppression in nonhomogeneous clutter"", 《2017 IEEE RADAR CONFERENCE》 * |
袁涛等: ""失配滤波器在数字脉压旁瓣抑制技术中的应用"", 《火控雷达技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113567978A (zh) * | 2021-07-29 | 2021-10-29 | 电子科技大学 | 一种多基分布式雷达协同成像方法 |
CN113567978B (zh) * | 2021-07-29 | 2023-04-25 | 电子科技大学 | 一种多基分布式雷达协同成像方法 |
CN116990773A (zh) * | 2023-09-27 | 2023-11-03 | 广州辰创科技发展有限公司 | 基于自适应门限的低慢小目标检测方法、装置和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109581366B (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6809682B1 (en) | Method and device for the detection and track of targets in high clutter | |
CN102156279B (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
US7535408B2 (en) | Apparatus and methods for detection of multiple targets within radar resolution cell | |
Hu et al. | Optimal signal processing in ground-based forward scatter micro radars | |
US8884810B2 (en) | Compact beacon radar and full ATC services system | |
CN111736126B (zh) | 一种基于低成本安防雷达的信号处理方法 | |
KR102013205B1 (ko) | 레이더 신호 처리 모의 장치 및 방법 | |
Budillon et al. | Along-track interferometric SAR systems for ground-moving target indication: Achievements, potentials, and outlook | |
Aubry et al. | Radar detection and range estimation using oversampled data | |
EP3326000B1 (en) | Obstacle detection radar using a polarization test | |
WO2010039299A1 (en) | Counter target acquisition radar and acoustic adjunct for classification | |
Hao et al. | Advances in Adaptive Radar Detection and Range Estimation | |
CN109581366A (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 | |
EP0851238A2 (en) | Doppler radar | |
Abdalla et al. | Design and Implementation of Proposed Low-Cost Dual-Channel IF Receiver for SSR | |
Maresca et al. | Theoretical limitation of the sea on the detection of low Doppler targets by over-the-horizon radar | |
CN1321331C (zh) | 一种利用海洋回波进行阵列通道校正的方法 | |
Boyle et al. | Comparison of monostatic and bistatic bearing estimation performance for low RCS targets | |
Fabrizio | High frequency over-the-horizon radar | |
WO2005052630A2 (en) | Doppler polarimetric detection and classification methods | |
KR20190124487A (ko) | 레이다 시스템에서 DoA 추정을 위한 신호개수 판별 방법 및 장치 | |
Zhao et al. | Using alpha-phase modulation method to solve range ambiguity for high frequency surface wave radar | |
Sedivy et al. | Doppler frequency estimation using moving target detection filter bank | |
Fabrizio et al. | Passive radar in the high frequency band | |
Barès et al. | A multifrequency HF-VHF radar system for aircraft identification |
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 |