CN108562877B - 一种基于信号包络特征的欺骗式干扰抑制方法 - Google Patents

一种基于信号包络特征的欺骗式干扰抑制方法 Download PDF

Info

Publication number
CN108562877B
CN108562877B CN201810101471.XA CN201810101471A CN108562877B CN 108562877 B CN108562877 B CN 108562877B CN 201810101471 A CN201810101471 A CN 201810101471A CN 108562877 B CN108562877 B CN 108562877B
Authority
CN
China
Prior art keywords
target
threshold
envelope
value
time domain
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
Application number
CN201810101471.XA
Other languages
English (en)
Other versions
CN108562877A (zh
Inventor
郑坚
程增飞
赵春光
刘子威
童建文
翟海涛
白俊奇
王远程
龙超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Lesi Electronic Equipment Co., Ltd.
Original Assignee
Nanjing Laisi Electronic Equipment Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Laisi Electronic Equipment Co ltd filed Critical Nanjing Laisi Electronic Equipment Co ltd
Priority to CN201810101471.XA priority Critical patent/CN108562877B/zh
Publication of CN108562877A publication Critical patent/CN108562877A/zh
Application granted granted Critical
Publication of CN108562877B publication Critical patent/CN108562877B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/2813Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于信号包络特征的欺骗式干扰抑制方法,主要解决现有技术不能有效对欺骗式干扰进行抑制的问题,其实现步骤为:1)根据系统参数确定背景噪声功率,相干处理时间内的重复周期个数,真实目标回波时域和空域包络特征门限参数;2)采用慢门限检测方法检测目标,得到过门限检测点的距离单元;3)提取过门限检测点的时域包络并计算特征参数;4)根据时域特征参数对真假目标进行第一次判别;5)对第一次判别中为真目标的点迹,提取其空域包络并计算其特征参数。6)利用空域特征参数对一次判别结果进行二次判断。本发明利用真假目标时域和空域包络的特征差异,通过对其进行量化,实现了对欺骗式假目标干扰的有效抑制。

Description

一种基于信号包络特征的欺骗式干扰抑制方法
技术领域
本发明属于雷达干扰抑制领域,尤其涉及一种基于信号包络特征的欺骗式干扰抑制方法。
背景技术
随着电磁环境的日益复杂,雷达系统受到越来越多的干扰。为了保证系统的有效工作,要求其必须具备一定的抗干扰能力。广义上而言干扰可以分为无源干扰和有源干扰,无源干扰一般是由地物、箔条等无源反射体反射雷达发射信号造成的,而有源干扰是通过有源器件人为故意产生的。根据对雷达的干扰原理的不同,有源干扰可以分为欺骗式干扰和压制式干扰。欺骗式干扰通过应答或转发侦测到的雷达发射信号来产生类似真实目标回波的假目标回波信号,来达到干扰雷达系统的目的。压制式干扰通过发射大功率的类噪声信号来淹没真实目标回波,从而实现目标干扰。欺骗式干扰是未来灵巧式干扰的基础,其对雷达系统具有严重的威胁,因此研究有效的抗欺骗式干扰方法具有重要意义。
目前,抑制欺骗式干扰的技术手段主要有频率和波形捷变、超低旁瓣天线、旁瓣相消和匿影等。频率和波形捷变技术通过实时改变雷达系统的工作频率和发射波形使得干扰系统不能有效测量出系统的工作参数,从而抑制假目标的产生。但是随着储频转发技术的发展,产生假目标干扰信号不再需要对信号参数进行测量,降低了该方法的有效性。超低旁瓣天线技术通过降低天线旁瓣来减弱入射干扰信号的能量,从而对由旁瓣进入系统的干扰具有较好的抑制作用。旁瓣相消和匿影技术通过设置辅助天线,利用干扰信号在主辅天线内的空域相关性来对其进行抑制。另外,上述技术都仅对由天线旁瓣进入的干扰有效,而对由主瓣干扰完全失效。在无源侦测领域,Adam K.等人在文献“Radar emitter recognitionusing intrapulse data,IEEE Transactions Theory,2005,20(5):435–438”中提出了采用脉冲特征进行辐射源识别的思路。在上述思路的基础上,王宏伟等人在文献“基于脉冲包络前沿波形的雷达辐射源个体识别,航天电子对抗,2009,25(2):35–38”给出了一种基于脉冲包络前沿特征的辐射源识别方法。基于上述思想,由于假目标干扰和真实目标信号本质上是由不同的辐射源产生的,其信号特征间必然存在差别,因此可以采用信号包络特征对干扰信号进行抑制。
发明内容
本发明的目的在于针对上述已有技术的不足,提出一种基于单个回波信号的时域和多个回波信号的空域包络特征的欺骗式干扰抑制方法,以改善雷达系统抗欺骗式(特别是主瓣欺骗式干扰)干扰的能力。
实现本发明的技术方案是:通过提取检测到的目标信号的时域和空域包络,计算出信号时域和空域包络的特征参数,对假目标和真实目标信号进行判别,从而实现对假目标干扰的抑制,其包括如下步骤:
步骤1,根据雷达系统参数确定背景噪声功率、相干处理时间内的重复周期个数、真实目标回波时域和空域包络特征门限参数;
步骤2,采用慢门限检测方法检测目标,得到过门限检测点的距离单元;
步骤3,提取过门限检测点的时域包络并计算时域包络对应的特征参数;
步骤4,根据时域包络对应的特征参数对真假目标进行第一次判别;
步骤5,对第一次判别中为真目标的点迹,提取其空域包络并计算其空域包络特征参数;
步骤6,利用空域包络特征参数对一次判别结果进行二次判断。
步骤1包括如下步骤:
步骤1-1,若在雷达系统休止期内的第i个采样信号点为yi,i=1,2,…,I,I为采样信号点总数,则得到背景噪声功率为
Figure GDA0002519020340000021
步骤1-2,对于机扫雷达,若雷达系统发射脉冲周期为Z,雷达系统波束宽度为B,扫描速率为θ,则得到相干处理时间内的重复周期个数
Figure GDA0002519020340000022
其中
Figure GDA0002519020340000023
表示向下取整;对于相控阵扫描雷达,若雷达系统驻留时间为Ti,发射脉冲周期为Z,则重复周期个数
Figure GDA0002519020340000024
步骤1-3,设定真实目标回波时域包络特征门限参数αt和目标回波空域包络特征门限参数αs
步骤2包括如下步骤:
步骤2-1,初始化目标集合U为空集合,设定目标个数Q=0,初始化处理周期计数器m=1,并设定总处理周期数M;
步骤2-2,接收第m个周期的回波信号xm,根据背景噪声功率σn采用慢门限检测方法对xm内的目标进行检测,计算慢检测门限G=30σn,则由此对回波信号xm内的目标进行慢门限检测,即若|xmh|2≥G,则第m个距离单元内存在目标,否则不存在目标,其中xmh为xm的第h个元素;若xmh为在xm检测到的第k个过门限的检测点,则第m个周期内检测到的第k个检测点umk所在距离单元lmk=h,k=1,2,…,Km,其中Km为第m个周期内检测到的目标个数。
步骤3包括如下步骤:
步骤3-1,初始化目标个数计数器k=1;
步骤3-2,由回波信号xm,提取出第k个目标的时域包络信号
Figure GDA0002519020340000031
Figure GDA0002519020340000032
其中
Figure GDA0002519020340000033
为xm的第lmk个元素;
步骤3-3,找出ptmk中值最大的点,记其幅度和下标分别为Am1和Im1,在ptmk中寻找下标小于Im1且幅度最接近Am1/10的点的幅度和下标,并分别记为Am2和Im2,则由此计算第k个目标的时域包络特征参数αtmk=(Am1-Am2)/(Im1-Im2);
步骤3-4,将计数器k的值加1,并判断k是否大于Km,如果是则执行步骤4,否则返回步骤3-2。
步骤4包括如下步骤:
步骤4-1,根据真实目标回波时域包络特征门限参数αt,对过检测门限的目标进行判别,即若αtmk≤αt,则判别目标umk为假目标,否则判别目标umk为真实目标;
步骤4-2,记第m个周期内判定为假的目标个数为S,第s个判定为假的目标记为ums,其对应的距离单元为lms,s=1,2,…,S;若目标集合U为空,则执行步骤4-3,否则执行步骤4-4;
步骤4-3,将第m个周期内判定为假的目标加入集合U中,改变集合U中目标个数Q的值为S,集合U中第q个元素记为uq,其对应的距离单元为lq,q=1,2,…,Q,返回步骤2-2;
步骤4-4,初始化目标个数计数器s=1,建立临时目标集合UT,初始化临时目标个数T=0,建立标志数组F=[f1,f2,…,fq,…,fQ],标志数组的第q个元素fq=0,q=1,2,…,Q;
步骤4-5,查找集合U中是否存在元素uq满足|lq-lms|≤rg,其中rg为距离门限,其值一般取为2~5,若存在则改变lq的值为(lms+lq)/2,将标志数组的第q个元素fq置1,否则将ums加入临时目标集合UT中,并将临时目标个数T的值加1;
步骤4-6,将目标个数计数器s的值加1,并判断s的值是否大于S,如果是则执行步骤4-7,否则返回步骤4-5;
步骤4-7,提取出标志数组F中值为0的元素对应的目标,构成二次筛选目标集合UJ={uq|fq=0}={u1,…,uj,…,uJ},其中第j个目标uj对应的距离单元为lj,j取值为1~J,将{uq|fq=0}从集合U中移除,同时将临时目标集合UT中的元素并入集合U中,改变集合U中目标个数Q的值为Q-J+T;
步骤4-8,将检测周期计数器m的值加1,并判断检测周期计数器m是否大于N,如果是则执行步骤5,否则返回步骤2-2;
步骤5包括如下步骤:
步骤5-1,初始化二次筛选计数器j=1;
步骤5-2,根据第j个目标uj所在距离单元lj,从回波信号xm-N+1,xm-N+2,……,xm中提取出目标的空域包络数据
Figure GDA0002519020340000041
其中
Figure GDA0002519020340000042
为xm-N+1的第lj个元素;
步骤5-3,找出psmj中值最大的点,记其幅度和下标分别为Bm1和Dm1,在psmj中寻找下标小于Bm1且幅度最接近Bm1/10的点的幅度和下标,并分别记为Bm2和Dm2,则由此计算第j个目标的空域包络特征参数αsmj=(Bm1-Bm2)/(Dm1-Dm2);
步骤5-4,将计数器j的值加1,并判断j是否大于J,如果是则执行步骤6,否则返回步骤5-2。
步骤6包括如下步骤:
步骤6-1,根据目标回波空域包络特征门限参数αs,对过门限的检测目标进行判别,即若第j个目标的空域包络参数αsmj满足αsmj≤αs,则判别目标uj为假目标,否则判别目标uj为真实目标;
步骤6-2,将检测周期计数器m的值加1,并判断m的值是否大于M,若是结束流程,否则返回步骤2-2。
有益效果:1.本发明方法仅通过目标回波信号的时域和空域包络特征对欺骗式干扰目标进行识判别并抑制,对主瓣和旁瓣干扰都具有较好的抑制性能;
2.本发明方法运算量小,计算简单,具有良好的实时性。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述或其他方面的优点将会变得更加清楚。
图1为本发明方法的实施流程图;
图2a~图2d为实测目标回波时域包络和空域包络示意图;
图3a和图3b为目标时域和空域包络特征二维分布图;
图4a~图4d为采用本发明方法进行欺骗式干扰抑制前后的点迹分布图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
参照图1,本发明进行欺骗式干扰抑制的步骤如下:
步骤1,根据雷达系统参数确定背景噪声功率、相干处理时间内的重复周期个数、真实目标回波时域和空域包络特征门限参数;
步骤2,采用慢门限检测方法检测目标,得到过门限检测点的距离单元;
步骤3,提取过门限检测点的时域包络并计算时域包络对应的特征参数;
步骤4,根据时域包络对应的特征参数对真假目标进行第一次判别;
步骤5,对第一次判别中为真目标的点迹,提取其空域包络并计算其空域包络特征参数;
步骤6,利用空域包络特征参数对一次判别结果进行二次判断。
步骤1包括:
步骤1-1,若在系统休止期内的采样信号点为yi,i=1,2,…,I,则由此可以估计得到背景噪声功率为
Figure GDA0002519020340000061
步骤1-2,对于机扫雷达,若雷达系统发射脉冲周期为Z,雷达系统波束宽度为B,扫描速率为θ,则得到相干处理时间内的重复周期个数
Figure GDA0002519020340000062
其中
Figure GDA0002519020340000063
表示向下取整;对于相控阵扫描雷达,若雷达系统驻留时间为Ti,发射脉冲周期为Z,则重复周期个数
Figure GDA0002519020340000064
步骤1-3,设定真实目标回波时域包络特征门限参数αt和目标回波空域包络特征门限参数αs
步骤2包括如下步骤:
步骤2-1,初始化目标集合U为空集合,设定目标个数Q=0,初始化处理周期计数器m=1,并设定总处理周期数M;
步骤2-2,接收第m个周期的回波信号xm,根据背景噪声功率σn采用慢门限检测方法对xm内的目标进行检测,计算慢检测门限G=30σn,则由此对回波信号xm内的目标进行慢门限检测,即若|xmh|2≥G,则第m个距离单元内存在目标,否则不存在目标,其中xmh为xm的第h个元素;若xmh为在xm检测到的第k个过门限的检测点,则第m个周期内检测到的第k个检测点umk所在距离单元lmk=h,k=1,2,…,Km,其中Km为第m个周期内检测到的目标个数。
步骤3包括如下步骤:
步骤3-1,初始化目标个数计数器k=1;
步骤3-2,由回波信号xm,提取出第k个目标的时域包络信号
Figure GDA0002519020340000065
Figure GDA0002519020340000071
其中
Figure GDA0002519020340000072
为xm的第lmk个元素;
步骤3-3,找出ptmk中值最大的点,记其幅度和下标分别为Am1和Im1,在ptmk中寻找下标小于Im1且幅度最接近Am1/10的点的幅度和下标,并分别记为Am2和Im2,则由此计算第k个目标的时域包络特征参数αtmk=(Am1-Am2)/(Im1-Im2);
步骤3-4,将计数器k的值加1,并判断k是否大于Km,如果是则执行步骤4,否则返回步骤3-2。
步骤4包括如下步骤:
步骤4-1,根据真实目标回波时域包络特征门限参数αt,对过检测门限的目标进行判别,即若αtmk≤αt,则判别目标umk为假目标,否则判别目标umk为真实目标;
步骤4-2,记第m个周期内判定为假的目标个数为S,第s个判定为假的目标记为ums,其对应的距离单元为lms,s=1,2,…,S;若目标集合U为空,则执行步骤4-3,否则执行步骤4-4;
步骤4-3,将第m个周期内判定为假的目标加入集合U中,改变集合U中目标个数Q的值为S,集合U中第q个元素记为uq,其对应的距离单元为lq,q=1,2,…,Q,返回步骤2-2;
步骤4-4,初始化目标个数计数器s=1,建立临时目标集合UT,初始化临时目标个数T=0,建立标志数组F=[f1,f2,…,fq,…,fQ],标志数组的第q个元素fq=0,q=1,2,…,Q;
步骤4-5,查找集合U中是否存在元素uq满足|lq-lms|≤rg,其中rg为距离门限,其值一般取为2~5,若存在则改变lq的值为(lms+lq)/2,将标志数组的第q个元素fq置1,否则将ums加入临时目标集合UT中,并将临时目标个数T的值加1;
步骤4-6,将目标个数计数器s的值加1,并判断s的值是否大于S,如果是则执行步骤4-7,否则返回步骤4-5;
步骤4-7,提取出标志数组F中值为0的元素对应的目标,构成二次筛选目标集合UJ={uq|fq=0}={u1,…,uj,…,uJ},其中第j个目标uj对应的距离单元为lj,j取值为1~J,将{uq|fq=0}从集合U中移除,同时将临时目标集合UT中的元素并入集合U中,改变集合U中目标个数Q的值为Q-J+T;
步骤4-8,将检测周期计数器m的值加1,并判断检测周期计数器m是否大于N,如果是则执行步骤5,否则返回步骤2-2;
步骤5包括如下步骤:
步骤5-1,初始化二次筛选计数器j=1;
步骤5-2,根据第j个目标uj所在距离单元lj,从回波信号xm-N+1,xm-N+2,……,xm中提取出目标的空域包络数据
Figure GDA0002519020340000081
其中
Figure GDA0002519020340000082
为xm-N+1的第lj个元素;
步骤5-3,找出psmj中值最大的点,记其幅度和下标分别为Bm1和Dm1,在psmj中寻找下标小于Bm1且幅度最接近Bm1/10的点的幅度和下标,并分别记为Bm2和Dm2,则由此计算第j个目标的空域包络特征参数αsmj=(Bm1-Bm2)/(Dm1-Dm2);
步骤5-4,将计数器j的值加1,并判断j是否大于J,如果是则执行步骤6,否则返回步骤5-2。
步骤6包括如下步骤:
步骤6-1,根据目标回波空域包络特征门限参数αs,对过门限的检测目标进行判别,即若第j个目标的空域包络参数αsmj满足αsmj≤αs,则判别目标uj为假目标,否则判别目标uj为真实目标;
步骤6-2,将检测周期计数器m的值加1,并判断m的值是否大于M,若是结束流程,否则返回步骤2-2。
1)实验条件:
本发明数值实验中软件仿真平台为MATLAB R2012a,实验基于两型雷达的实测目标和干扰数据进行。
2)实验内容及结果:
仿真1,为了对本发明方法的合理性进行说明,首先从两型雷达的实测数据中提取出多批欺骗式干扰和真实目标回波的时域和空域包络。图2a~图2d给出了两型雷达实测数据的时域和空域包络特征,其中图2a为雷达1实测干扰和信号的归一化时域包络,图2b为雷达1实测干扰和信号的归一化空域包络,图2c为雷达2实测干扰和信号的归一化时域包络,图2d为雷达2实测干扰和信号的归一化空域包络。
由图2a~图2d可知,欺骗式干扰和真实目标回波的时域和空域包络存在一定差别,通过对其进行量化,可以对欺骗式干扰和真实目标进行判别。上述时域特征差异出现的原因是由于产生欺骗式干扰和真实目标回波的收发器件的特性不同,而在对信号进行匹配滤波时采用的是同真实目标信号匹配的滤波器系数,因此其对欺骗式干扰信号存在适配,使得欺骗式干扰信号的时域包络展宽。而导致空域特征不同的原因是发射欺骗式干扰的干扰机对雷达系统的脉冲重复周期进行估计时存在估计误差,使得不同周期内的干扰存在一定的异步特性,从而导致多周期内干扰信号的空域包络存在较大起伏。
实验2,对实验1中提取的多批欺骗式干扰和目标回波的时间和幅度包络,分别采用本发明方法对其特征进行量化,得到雷达1和雷达2的时域与空域包络特征参数二维分布图如图3a和图3b所示,其中图3a为雷达1采集干扰和目标信号的特征分布图,图3b为雷达2采集干扰和目标信号的特征分布图。
由图3a和图3b可知,采用本发明方法中的量化方法对信号的时域及空域包络特征进行量化,可以对欺骗式干扰和目标进行有效判别。
实验3采用本发明方法对两型雷达采集到的目标和密集欺骗式干扰的数据进行处理,试验中设定雷达1的时域包络特征门限参数αt=0.035,空域包络特征门限参数αs=0.085,雷达2的时域包络特征门限参数αt=0.065,空域包络特征门限参数αs=0.09,两雷达相干处理时间内的重复周期个数N分别为24和16,得到原始过检测门限的点迹及经过本发明方法处理后的剩余点迹如图4a~图4d所示。其中图4a和图4b为雷达1采集数据的原始点迹和采用本发明方法处理后的点迹,图4c和图4d为雷达2采集数据的原始点迹和采用本发明方法处理后的点迹。另外,为了对本发明方法的干扰抑制性能进行量化,定义抑制掉的干扰点迹数和总干扰点迹数的比值为干扰抑制率,定义抑制掉的真实目标点迹数与真实目标点迹数的比值为信号损失率,得到图4a~图4d中两型雷达的干扰抑制率和信号损失率如表1所示。
表1本发明方法的干扰抑制率和信号损失率
Figure GDA0002519020340000101
由图4a~图4d和表1可知,本发明方法能够有效对欺骗式干扰进行抑制,且对真实目标信号损失较小,是一种切实可行的欺骗式干扰抑制方法。综上可知,本发明方法能有效对欺骗式干扰信号进行抑制。
本发明提供了一种基于信号包络特征的欺骗式干扰抑制方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。

Claims (7)

1.一种基于信号包络特征的欺骗式干扰抑制方法,其特征在于,包括以下步骤:
步骤1,根据雷达系统参数确定背景噪声功率、相干处理时间内的重复周期个数、真实目标回波时域和空域包络特征门限参数;
步骤2,采用慢门限检测方法检测目标,得到过门限检测点的距离单元;
步骤3,提取过门限检测点的时域包络并计算时域包络对应的特征参数;
步骤4,根据时域包络对应的特征参数对真假目标进行第一次判别;
步骤5,对步骤4判别中为真目标的点迹,提取其空域包络并计算其空域包络特征参数;
步骤6,利用空域包络特征参数对一次判别结果进行二次判断。
2.根据权利要求1所述的方法,其特征在于,步骤1包括如下步骤:
步骤1-1,若在雷达系统休止期内的第i个采样信号点为yi,i=1,2,…,I,I为采样信号点总数,则得到背景噪声功率为
Figure FDA0002519020330000011
步骤1-2,对于机扫雷达,若雷达系统发射脉冲周期为Z,雷达系统波束宽度为B,扫描速率为θ,则得到相干处理时间内的重复周期个数
Figure FDA0002519020330000012
其中
Figure FDA0002519020330000013
表示向下取整;对于相控阵扫描雷达,若雷达系统驻留时间为Ti,发射脉冲周期为Z,则重复周期个数
Figure FDA0002519020330000014
步骤1-3,设定真实目标回波时域包络特征门限参数αt和目标回波空域包络特征门限参数αs
3.根据权利要求2所述的方法,其特征在于,步骤2包括如下步骤:
步骤2-1,初始化目标集合U为空集合,设定目标个数Q=0,初始化处理周期计数器m=1,并设定总处理周期数M;
步骤2-2,接收第m个周期的回波信号xm,根据背景噪声功率σn采用慢门限检测方法对xm内的目标进行检测,计算慢检测门限G=30σn,则由此对回波信号xm内的目标进行慢门限检测,即若|xmh|2≥G,则第m个距离单元内存在目标,否则不存在目标,其中xmh为xm的第h个元素;若xmh为在xm检测到的第k个过门限的检测点,则第m个周期内检测到的第k个检测点umk所在距离单元lmk=h,k=1,2,…,Km,其中Km为第m个周期内检测到的目标个数。
4.根据权利要求3所述的方法,其特征在于,步骤3包括如下步骤:
步骤3-1,初始化目标个数计数器k=1;
步骤3-2,由回波信号xm,提取出第k个目标的时域包络信号
Figure FDA0002519020330000021
Figure FDA0002519020330000022
其中
Figure FDA0002519020330000023
为xm的第lmk个元素;
步骤3-3,找出ptmk中值最大的点,记其幅度和下标分别为Am1和Im1,在ptmk中寻找下标小于Im1且幅度最接近Am1/10的点的幅度和下标,并分别记为Am2和Im2,则由此计算第k个目标的时域包络特征参数αtmk=(Am1-Am2)/(Im1-Im2);
步骤3-4,将计数器k的值加1,并判断k是否大于Km,如果是则执行步骤4,否则返回步骤3-2。
5.根据权利要求4所述的方法,其特征在于,步骤4包括如下步骤:
步骤4-1,根据真实目标回波时域包络特征门限参数αt,对过检测门限的目标进行判别,即若αtmk≤αt,则判别目标umk为假目标,否则判别目标umk为真实目标;
步骤4-2,记第m个周期内判定为假的目标个数为S,第s个判定为假的目标记为ums,其对应的距离单元为lms,s=1,2,…,S;若目标集合U为空,则执行步骤4-3,否则执行步骤4-4;
步骤4-3,将第m个周期内判定为假的目标加入集合U中,改变集合U中目标个数Q的值为S,集合U中第q个元素记为uq,其对应的距离单元为lq,q=1,2,…,Q,返回步骤2-2;
步骤4-4,初始化目标个数计数器s=1,建立临时目标集合UT,初始化临时目标个数T=0,建立标志数组F=[f1,f2,…,fq,…,fQ],标志数组的第q个元素fq=0,q=1,2,…,Q;
步骤4-5,查找集合U中是否存在元素uq满足|lq-lms|≤rg,其中rg为距离门限,若存在则改变lq的值为(lms+lq)/2,将标志数组的第q个元素fq置1,否则将ums加入临时目标集合UT中,并将临时目标个数T的值加1;
步骤4-6,将目标个数计数器s的值加1,并判断s的值是否大于S,如果是则执行步骤4-7,否则返回步骤4-5;
步骤4-7,提取出标志数组F中值为0的元素对应的目标,构成二次筛选目标集合UJ={uq|fq=0}={u1,…,uj,…,uJ},其中第j个目标uj对应的距离单元为lj,j取值为1~J,将{uq|fq=0}从集合U中移除,同时将临时目标集合UT中的元素并入集合U中,改变集合U中目标个数Q的值为Q-J+T;
步骤4-8,将检测周期计数器m的值加1,并判断检测周期计数器m是否大于N,如果是则执行步骤5,否则返回步骤2-2。
6.根据权利要求5所述的方法,其特征在于,步骤5包括如下步骤:
步骤5-1,初始化二次筛选计数器j=1;
步骤5-2,根据第j个目标uj所在距离单元lj,从回波信号xm-N+1,xm-N+2,……,xm中提取出目标的空域包络数据
Figure FDA0002519020330000031
其中
Figure FDA0002519020330000032
为xm-N+1的第lj个元素;
步骤5-3,找出psmj中值最大的点,记其幅度和下标分别为Bm1和Dm1,在psmj中寻找下标小于Bm1且幅度最接近Bm1/10的点的幅度和下标,并分别记为Bm2和Dm2,则由此计算第j个目标的空域包络特征参数αsmj=(Bm1-Bm2)/(Dm1-Dm2);
步骤5-4,将计数器j的值加1,并判断j是否大于J,如果是则执行步骤6,否则返回步骤5-2。
7.根据权利要求6所述的方法,其特征在于,步骤6包括如下步骤:
步骤6-1,根据目标回波空域包络特征门限参数αs,对过门限的检测目标进行判别,即若第j个目标的空域包络参数αsmj满足αsmj≤αs,则判别目标uj为假目标,否则判别目标uj为真实目标;
步骤6-2,将检测周期计数器m的值加1,并判断m的值是否大于M,若是结束流程,否则返回步骤2-2。
CN201810101471.XA 2018-02-01 2018-02-01 一种基于信号包络特征的欺骗式干扰抑制方法 Active CN108562877B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810101471.XA CN108562877B (zh) 2018-02-01 2018-02-01 一种基于信号包络特征的欺骗式干扰抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810101471.XA CN108562877B (zh) 2018-02-01 2018-02-01 一种基于信号包络特征的欺骗式干扰抑制方法

Publications (2)

Publication Number Publication Date
CN108562877A CN108562877A (zh) 2018-09-21
CN108562877B true CN108562877B (zh) 2020-07-24

Family

ID=63530158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810101471.XA Active CN108562877B (zh) 2018-02-01 2018-02-01 一种基于信号包络特征的欺骗式干扰抑制方法

Country Status (1)

Country Link
CN (1) CN108562877B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581304B (zh) * 2018-12-05 2020-08-07 航天南湖电子信息技术股份有限公司 一种反侦察抗干扰装置及其方法
CN109917367B (zh) * 2019-03-19 2023-11-14 安徽雷炎电子科技有限公司 一种雷达目标距离和方位两维包络信息的匹配方法及其应用
CN113037427B (zh) * 2021-03-03 2022-08-16 四川九洲空管科技有限责任公司 一种应用于敌我识别系统的抗欺骗应答方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5943001A (en) * 1974-03-08 1999-08-24 The United States Of America As Represented By The Secretary Of The Navy Swept comb method of jamming
CN103760530A (zh) * 2014-01-15 2014-04-30 西安电子科技大学 基于信号锥的雷达压制式干扰和欺骗式干扰识别方法
CN103954943A (zh) * 2014-05-13 2014-07-30 西安电子科技大学 网络化雷达系统抗欺骗式干扰方法
CN104991232A (zh) * 2015-06-26 2015-10-21 西安电子科技大学 目标信号相关下信号级融合的组网雷达抗欺骗式干扰方法
CN104991233A (zh) * 2015-06-29 2015-10-21 西安电子科技大学 基于信号级融合的组网雷达抗欺骗式干扰方法
CN105259541A (zh) * 2015-11-26 2016-01-20 西安电子科技大学 一种多站雷达抗有源欺骗式干扰的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5943001A (en) * 1974-03-08 1999-08-24 The United States Of America As Represented By The Secretary Of The Navy Swept comb method of jamming
CN103760530A (zh) * 2014-01-15 2014-04-30 西安电子科技大学 基于信号锥的雷达压制式干扰和欺骗式干扰识别方法
CN103954943A (zh) * 2014-05-13 2014-07-30 西安电子科技大学 网络化雷达系统抗欺骗式干扰方法
CN104991232A (zh) * 2015-06-26 2015-10-21 西安电子科技大学 目标信号相关下信号级融合的组网雷达抗欺骗式干扰方法
CN104991233A (zh) * 2015-06-29 2015-10-21 西安电子科技大学 基于信号级融合的组网雷达抗欺骗式干扰方法
CN105259541A (zh) * 2015-11-26 2016-01-20 西安电子科技大学 一种多站雷达抗有源欺骗式干扰的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Radar emitter recognition using intrapulse data;Adam KAWALEC et al.;《15th International Conference on Microwaves, Radar and Wireless Communications (IEEE Cat. No.04EX824)》;20041122;第435—438页 *
Signal Fusion-Based Algorithms to Discriminate Between Radar Targets and Deception Jamming in Distributed Multiple-Radar Architectures;Shanshan Zhao et al.;《IEEE Sensors Journal》;20150603;第15卷(第11期);第6697—6706页 *
基于信号包络的辐射源细微特征提取方法;张国柱等;《系统工程与电子技术》;20060630;第28卷(第6期);第795—797、936页 *
基于脉冲包络前沿波形的雷达辐射源个体识别;王宏伟等;《航天电子对抗》;20090609;第25卷(第2期);第35—38页 *

Also Published As

Publication number Publication date
CN108562877A (zh) 2018-09-21

Similar Documents

Publication Publication Date Title
CN101881822B (zh) 一种针对共享谱雷达同频干扰的抑制方法
Orlando et al. Track-before-detect strategies for STAP radars
CN108562877B (zh) 一种基于信号包络特征的欺骗式干扰抑制方法
US8138963B1 (en) Method for detecting targets using space-time adaptive processing and shared knowledge of the environment
CN104267379A (zh) 一种基于波形设计的主被动雷达协同抗干扰方法
CN105137423B (zh) 一种穿墙雷达对多运动目标实时探测、分离的方法
CN103616729A (zh) 一种基于uwb生物雷达的多个人体目标估算方法及系统
CN105445710B (zh) 隐形雷达低截获概率性能评估方法
CN106054153A (zh) 一种基于分数阶变换的海杂波区目标检测与自适应杂波抑制方法
CN105510887A (zh) 一种杂波背景下对机载雷达的有源压制式干扰抑制方法
CN112014806A (zh) 一种复杂干扰场景下的机载雷达无意干扰抑制方法
CN105093196A (zh) 基于逆伽马纹理复合高斯模型下的相干检测方法
Zhang et al. Weak target detection within the nonhomogeneous ionospheric clutter background of HFSWR based on STAP
Du et al. NLOS target localization with an L-band UWB radar via grid matching
CN111999715B (zh) 异质杂波下目标知识辅助自适应融合检测方法
Wu et al. NLOS positioning for building layout and target based on association and hypothesis method
CN106353743A (zh) 匹配于等效形状参数的近最优雷达目标检测方法
CN112612007B (zh) 基于近场效应的超稀疏阵列机载雷达动目标距离解模糊方法
CN103913724B (zh) 基于先验地形覆盖数据的杂波抑制方法
CN111157966B (zh) 一种假目标干扰识别方法
CN109655793B (zh) 高频雷达从瞬态干扰中提取回波的方法
CN112763984A (zh) 一种雷达自适应旁瓣相消算法
Yamaguchi et al. CFAR detection from noncoherent radar echoes using Bayesian Theory
CN115144830B (zh) 一种强地物杂波干扰下的小目标雷达检测方法
CN108983227A (zh) 一种基于白化滤波的极化mimo雷达检测方法

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
TA01 Transfer of patent application right

Effective date of registration: 20181024

Address after: 210007 the 5 building of Tianan Digital City, 36 Yongfeng Avenue, Qinhuai District, Nanjing, Jiangsu.

Applicant after: Nanjing Lesi Electronic Equipment Co., Ltd.

Address before: 210007 1 East Street, alfalfa garden, Nanjing, Jiangsu

Applicant before: China Electronic Technology Corporation (Group) 28 Institute

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant