CN107153178B - 外辐射源雷达参考信号含有多径干扰时的目标检测方法 - Google Patents

外辐射源雷达参考信号含有多径干扰时的目标检测方法 Download PDF

Info

Publication number
CN107153178B
CN107153178B CN201710352530.6A CN201710352530A CN107153178B CN 107153178 B CN107153178 B CN 107153178B CN 201710352530 A CN201710352530 A CN 201710352530A CN 107153178 B CN107153178 B CN 107153178B
Authority
CN
China
Prior art keywords
signal
reference signal
time domain
echo
signals
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
CN201710352530.6A
Other languages
English (en)
Other versions
CN107153178A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201710352530.6A priority Critical patent/CN107153178B/zh
Publication of CN107153178A publication Critical patent/CN107153178A/zh
Application granted granted Critical
Publication of CN107153178B publication Critical patent/CN107153178B/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/41Details 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于雷达技术领域,公开了一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,包括:获取外辐射源雷达接收到的信号,包含参考信号和回波信号,且参考信号中含有多径干扰信号,回波信号中含有直达波信号和多径干扰信号;对参考信号和回波信号依次进行放大、混频、A/D采样和数字下变频处理,得到数字参考信号和数字回波信号;对数字参考信号进行盲均衡滤波,得到多径干扰抑制后的参考信号;对数字回波信号进行时域干扰相消,得到时域干扰相消后的回波信号;对时域干扰相消后的回波信号进行脉冲压缩处理,从而得到目标检测结果;可以实现抑制参考信号中的多径干扰成分,消除多径干扰造成的虚假目标,得到较好的目标检测结果。

Description

外辐射源雷达参考信号含有多径干扰时的目标检测方法
技术领域
本发明属于雷达技术领域,尤其涉及一种外辐射源雷达参考信号含有多径干扰时的目标检测方法。
背景技术
外辐射源雷达探测目标属于被动探测,依靠第三方辐射源信号进行目标探测。外辐射雷达接收站直接接收直达波信号作为参考信号,当参考信号中含有多径干扰成分时,直接利用含有多径干扰信号的参考信号做匹配滤波会导致虚假目标的出现。所以对多径干扰的抑制程度直接影响着外辐射雷达探测目标的性能。
外辐射源雷达系统采用双通道信号处理结构,即参考通道与回波通道,分别接收参考信号和回波信号。外辐射源雷达参考通道利用指向照射源的参考天线接收直达波信号作为参考信号,参考信号用于后续回波通道的时域干扰相消和匹配滤波处理,所以获取纯净的参考信号尤为重要。
由于近地建筑或山丘的折射效应,受到折射的直达波信号相对于直线传播的直达波信号就有了时间延迟,这些延迟后的直达波信号就形成多径干扰信号。参考通道和回波通道都会受到多径干扰的影响,严重影响外辐射雷达系统的目标探测性能。
针对回波通道中的直达波和多径干扰,可采用时域干扰相消的方法来消除回波通道中的直达波和多径干扰。目前国内外主要的时域干扰相消算法主要有:SMI(采样矩阵求逆,Sample Matrix Inversion)算法、 ECA-B(分段扩展相消,Batch Version ofExtensive Cancellation)算法、 LMS(最小均方,Least Mean Square)算法及其改进的算法等。
参考通道中的多径干扰是较难消除的,一般都采用空域滤波的方法,通过自适应的波束形成,在多径干扰来波方向形成方向图零点。但是,当多径信号与直达波信号来波方向比较接近时,多径信号与直达波信号在时域、空域和频域都是不可分的,目前还并没有行之有效的方法来消除参考信号中的多径干扰信号。
发明内容
针对上述现有技术的不足,本发明的目的在于提供一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,在外辐射源雷达参考信号中含有多径干扰时,可以实现抑制参考信号中的多径干扰成分,消除多径干扰造成的虚假目标,得到较好的目标检测结果。
实现本发明目的的技术思路是:来自第三方照射源的直达波信号和多径传播信道的都是未知的,在无法获取训练信号的情况下,通过传统滤波手段无法消除多径干扰。针对照射源信号具有恒模的特点,利用改进的恒模算法(CMA,Constant Modulus Algorithm)来对含有多径干扰信号的参考信号进行盲均衡处理,抑制参考信号中的多径干扰成分,改善外辐射源雷达系统的目标探测性能。
为达到上述目的,本发明采用如下技术方案予以实现。
一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,所述方法包括如下步骤:
步骤1,获取外辐射源雷达接收到的信号,所述外辐射源雷达接收到的信号包含参考信号和回波信号,且所述参考信号中含有多径干扰信号,所述回波信号中含有直达波信号和多径干扰信号;分别对所述参考信号和所述回波信号依次进行放大、混频、A/D采样和数字下变频处理,得到数字参考信号和数字回波信号;
步骤2,对所述数字参考信号进行盲均衡滤波,得到多径干扰抑制后的参考信号;
步骤3,根据所述多径干扰抑制后的参考信号对所述数字回波信号进行时域干扰相消,得到时域干扰相消后的回波信号;
步骤4,根据所述多径干扰抑制后的参考信号,对所述时域干扰相消后的回波信号进行脉冲压缩处理,从而得到目标检测结果。
本发明与现有技术相比具有以下优点:本发明将CMA算法改进并应用到外辐射源雷达抑制参考信号中多径干扰的背景下。传统滤波方法无法抑制参考信号中多径干扰,并且传统CMA算法对较小时延的多径干扰抑制效果较差,本发明改进盲均衡算法中的CMA算法,并对(含有多径干扰的)参考信号进行盲均衡处理,该改进算法具有收敛速度快、稳定性强的优点,可以更好抑制参考信号中的多径干扰。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种外辐射源雷达参考信号含有多径干扰时的目标检测方法的流程示意图;
图2为本发明实施例提供的外辐射源雷达双基地场景配置示意图;
图3为本发明仿真实验中外辐射源雷达参考信号含有多径干扰信号时的距离-多普勒示意图
图4为本发明仿真实验中外辐射源雷达利用传统CMA抑制参考信号中多径干扰后的距离-多普勒示意图
图5为本发明仿真实验中外辐射源雷达利用改进CMA抑制参考信号中多径干扰后的距离-多普勒示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,如图1所示,所述方法包括如下步骤:
步骤1,获取外辐射源雷达接收到的信号,所述外辐射源雷达接收到的信号包含参考信号和回波信号,且所述参考信号中含有多径干扰信号,所述回波信号中含有直达波信号和多径干扰信号;分别对所述参考信号和所述回波信号依次进行放大、混频、A/D采样和数字下变频处理,得到数字参考信号和数字回波信号。
在步骤1之前,需要进行场景设置,具体的,如图2所示,外辐射源雷达上设置有一副参考天线和一副回波接收天线,所述参考天线指向辐射源,所述回波接收天线指向观测区,且所述参考天线接收的信号为参考信号,所述回波接收天线接收的信号为回波信号,同时,所述场景中还包括第三方辐射源(即有源脉冲雷达),所述第三方辐射源设置于所述外辐射源雷达接收站的远场作为发射站,所述第三方辐射源用于发射信号。
步骤2,对所述数字参考信号进行盲均衡滤波,得到多径干扰抑制后的参考信号。
基本CMA算法的权值迭代公式中,FIR滤波器的权值向量w(n)依赖于n,并采用下面的最小均方(LMS)算法来更新:
w(n+1)=w(n)+μx(n)e*(n)
其中,μ为迭代更新步长。误差信号为e(n):e(n)=y(n)[|y(n)|22]
其中:
Figure BDA0001298206230000051
式中,s(n)表示源信号(即不含多径干扰的参考信号)。
本发明实施例中,步骤2具体包括如下子步骤:
(2a)设置具有单位范数的向量为盲均衡滤波器的初始权值向量 w(n)=[w-L(n),w-L+1(n),…,wL(n)]T;并构造盲均衡滤波器的输入信号 x(n)=sref(n),其中,sref(n)为数字参考信号,n表示离散时间序列;
令n的初值为L+1,且L+1≤n≤N-L,其中,N为数字参考信号sref(n) 的数据长度,盲均衡滤波器的阶数为2L+1;且N≥2L+1
(2b)计算盲均衡滤波器的输出信号y(n):
Figure BDA0001298206230000052
其中,k表示盲均衡滤波器的抽头序号,(·)*表示求共轭操作,(·)H表示求共轭转置操作;
(2c)计算盲均衡滤波器的误差信号e(n):
e(n)=y(n)[|y(n)|22]
其中,|·|2表示求模值平方操作,且γ2=1;
(2d)令n的值加1,采用下式对盲均衡滤波器的权值向量进行更新:
Figure BDA0001298206230000053
其中,α、β为正实数;取α=1.0,β=2.0。
(2e)若n<N-L,则返回子步骤(2b)继续执行;否则,将最后得到的盲均衡滤波器的输出信号y(n)作为多径干扰抑制后的参考信号
Figure BDA0001298206230000061
需要说明的是,改进的CMA算法,根据最小化干扰原理,权值向量应该以最小的增量来更新,这可以看作是一种约束优化。权值向量的自适应增量为:
△w(n+1)=w(n+1)-w(n)
利用拉格朗日乘子法来解这个约束优化问题,其代价函数为:
J(n)=||△w(n+1)||2+Re[λ*(d(n)-wH(n+1)x(n))]
Re表示取实部运算,λ为拉格朗日乘子,d(n)为期望信号,即不含多径干扰的参考信号;
对代价函数求梯度得到:
Figure BDA0001298206230000062
当梯度等于0时,得到下列等式:
Figure BDA0001298206230000063
进一步得到:
Figure BDA0001298206230000064
从而得到改进后的CMA算法的权值向量更新公式:
Figure BDA0001298206230000065
其中,α为正实数,来控制权值向量的增量;β为正常数,用来避免输入信号||x(n)||2值较小。
步骤3,根据所述多径干扰抑制后的参考信号对所述数字回波信号进行时域干扰相消,得到时域干扰相消后的回波信号。
步骤3具体包括如下子步骤:
(3a)记所述多径干扰抑制后的参考信号
Figure BDA0001298206230000071
和数字回波信号的信号长度为M,且M<N,将所述多径干扰抑制后的参考信号
Figure BDA0001298206230000072
和所述数字回波信号ssur(n)分别分成B段,每段多径干扰抑制后的参考信号长度为M/B,且每段数字回波信号的长度为M/B;
令b的初值为1,且b=1,...,B;
(3b)构造第b段多径干扰抑制后的参考信号及其时延张成的子空间矩阵
Figure BDA0001298206230000073
Figure BDA0001298206230000074
其中,p为数字回波信号中需要消除的多径干扰的时延长度,T表示转置;
(3c)第b段数字回波信号
Figure BDA0001298206230000075
为:
Figure BDA0001298206230000076
(3d)根据所述第b段多径干扰抑制后的参考信号及其时延张成的子空间矩阵
Figure BDA0001298206230000077
构造第b段数字回波信号
Figure BDA0001298206230000078
的时域相消子空间投影矩阵 wb,并根据所述时域相消子空间投影矩阵wb对第b段数字回波信号
Figure BDA0001298206230000079
进行时域干扰相消,得到时域干扰相消后的第b段数字回波信号
Figure BDA00012982062300000710
具体的,求解如下最小优化问题:
Figure BDA00012982062300000711
目标函数梯度为0的地方即最小值所在的位置:
Figure BDA0001298206230000081
整理得到第b段回波信号
Figure BDA0001298206230000082
的时域相消子空间投影矩阵:
Figure BDA0001298206230000083
经过ECA-B时域干扰相消之后的第b段的回波信号为:
Figure BDA0001298206230000084
(3e)令b的值加1,并重复执行子步骤(3b)至(3d),得到时域干扰相消后的B段数字回波信号,将所述时域干扰相消后的B段数字回波信号组合成一路信号,并作为时域干扰相消后的回波信号
Figure BDA0001298206230000085
即:
Figure BDA0001298206230000086
步骤4,根据所述多径干扰抑制后的参考信号,对所述时域干扰相消后的回波信号进行脉冲压缩处理,从而得到目标检测结果。
步骤4具体包括如下子步骤:
(4a)分别对所述多径干扰抑制后的参考信号
Figure BDA0001298206230000087
所述时域干扰相消后的回波信号
Figure BDA0001298206230000088
进行快速傅里叶变换,分别得到多径干扰抑制后的频域参考信号
Figure BDA0001298206230000089
和时域干扰相消后的频域回波信号
Figure BDA00012982062300000810
其中,FFT{·}表示快速傅里叶变换操作;
(4b)对所述多径干扰抑制后的频域参考信号Sref(f)和时域干扰相消后的频域回波信号Ssur(f)进行脉冲压缩处理,得到脉冲压缩处理的频域输出S0(f):
Figure BDA0001298206230000091
其中,(·)*表示取共轭操作,sd(n)表示多径干扰抑制后的参考信号中的时域直达波信号,n(n)表示参考通道时域噪声,sechoj(n)表示时域干扰相消后的回波信号中的第j个时域目标回波信号,
Figure BDA0001298206230000092
表示回波通道时域噪声,Sd(f)表示多径干扰抑制后的参考信号中的频域直达波信号, n(f)表示参考通道频域噪声,Sechoj(f)表示时域干扰相消后的回波信号中的第j个频域目标回波信号,
Figure BDA0001298206230000093
表示回波通道频域噪声,j=1,2,...,m, m表示时域干扰相消后的回波信号中包含的目标回波信号的个数,表示脉冲压缩处理后总的频域噪声信号;
(4c)将所述脉冲压缩处理的频域输出S0(f)变换到时域,得到时域输出s0(t):
Figure BDA0001298206230000094
其中,IFFT{·}表示快速逆傅里叶变换操作,sRj(n)表示脉冲压缩处理得到的第j个时域目标回波信号,sn(n)表示脉冲压缩处理后的时域噪声信号平台。
1、本发明实施例的仿真条件:
本发明实验中信号源为LFM脉冲信号,脉冲宽度为300μs,频率为 88MHz,带宽为2.5MHz,采样频率为5MHz,脉冲相干积累数64;参考信号中的直达波信噪比为35dB,参考信号中多径干扰信号个数为3,其相对直达波时延分别为24μs,30μs和44μs,其干噪比分别为30dB、29dB 和28dB。仿真目标1信噪比为-10dB,距离为40km,多普勒频移50Hz;仿真目标2信噪比为-8dB,距离为50km,多普勒频移100Hz。发射站高度分别为200m,初始基线距离为13Km。
2、本发明实验的仿真结果分析:
图3为在参考信号含有多径干扰情况下,外辐射源雷达目标检测的结果。图3(a)为距离-多普勒示意图,图3(b)为距离-多普勒示意图的时延维。从图中明显看出,距离位于40km,50km,多普勒为50Hz,100Hz的真实目标1和2可以被检测到,但是在真实目标所在的多普勒单元内,出现三个虚假目标,其与真实目标在时延上相差24μs,30μs和44μs,这些虚假目标正是参考信号中多径干扰与目标回波脉冲压缩处理形成的。在参考信号含有多径干扰的情况下,在真实目标所在的多普勒单元出现虚假目标,这些虚假目标超过检测门限,使系统的虚假概率升高,严重影响了雷达系统的探测性能。
图4为外辐射源雷达采用传统CMA算法抑制参考信号中多径干扰后目标检测结果。图4(a)为距离-多普勒示意图,图4(b)为距离-多普勒示意图的时延维。可以看出,真实目标1和2可以被检测到,说明CMA算法对真实目标的检测并没有产生影响。同时,多径干扰造成的虚假目标幅度降低,表明CMA算法对抑制多径干扰起到一定效果。但是距离真实目标最近的虚假目标峰值依然较高,说明传统CMA算法对时延较小的多径干扰的抑制效果比较差。
图5为外辐射源雷达采用改进CMA算法抑制参考信号中多径干扰后目标检测结果。图5(a)为距离-多普勒示意图,图5(b)为距离-多普勒示意图的时延维。可以看出,图中只有真实目标1和2的峰值,多径干扰造成的虚假目标峰值淹没在噪声平台之下,无法被检测到。表明改进的CMA算法对多径干扰具有更好的抑制作用,消除了多径干扰造成的虚假目标,改善了外辐射源雷达系统的探测性能。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (4)

1.一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,其特征在于,所述方法包括如下步骤:
步骤1,获取外辐射源雷达接收到的信号,所述外辐射源雷达接收到的信号包含参考信号和回波信号,且所述参考信号中含有多径干扰信号,所述回波信号中含有直达波信号和多径干扰信号;分别对所述参考信号和所述回波信号依次进行放大、混频、A/D采样和数字下变频处理,得到数字参考信号和数字回波信号;
步骤2,对所述数字参考信号进行盲均衡滤波,得到多径干扰抑制后的参考信号;
步骤3,根据所述多径干扰抑制后的参考信号对所述数字回波信号进行时域干扰相消,得到时域干扰相消后的回波信号;
步骤4,根据所述多径干扰抑制后的参考信号,对所述时域干扰相消后的回波信号进行脉冲压缩处理,从而得到目标检测结果;
其中,步骤2具体包括如下子步骤:
(2a)设置具有单位范数的向量为盲均衡滤波器的初始权值向量w(n)=[w-L(n),w-L+1(n),...,wL(n)]T,(·)T表示求转置操作;并构造盲均衡滤波器的输入信号x(n)=sref(n),其中,sref(n)为数字参考信号,n表示离散时间序列;
令n的初值为L+1,且L+1≤n≤N-L,其中,N为数字参考信号sref(n)的数据长度,盲均衡滤波器的阶数为2L+1;且N≥2L+1
(2b)计算盲均衡滤波器的输出信号y(n):
Figure FDA0002276854510000011
其中,wk(n)表示权值向量中的元素,k表示盲均衡滤波器的抽头序号,(·)*表示求共轭操作,(·)H表示求共轭转置操作,x(n-k)表示将盲均衡滤波器的输入信号x(n)延迟k个抽头序号的后的信号;
(2c)计算盲均衡滤波器的误差信号e(n):
e(n)=y(n)[|y(n)|22]
其中,|·|2表示求模值平方操作,且γ2=1;
(2d)令n的值加1,采用下式对盲均衡滤波器的权值向量进行更新:
Figure FDA0002276854510000021
其中,α、β为正实数;
(2e)若n<N-L,则返回子步骤(2b)继续执行;否则,将迭代停止时得到的盲均衡滤波器的输出信号y(n)作为多径干扰抑制后的参考信号
Figure FDA0002276854510000022
2.根据权利要求1所述的一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,其特征在于,步骤3具体包括如下子步骤:
(3a)记所述多径干扰抑制后的参考信号
Figure FDA0002276854510000023
和数字回波信号的信号长度分别为M,且M<N,将所述多径干扰抑制后的参考信号
Figure FDA0002276854510000024
和所述数字回波信号ssur(n)分别分成B段,每段多径干扰抑制后的参考信号长度为M/B,且每段数字回波信号的长度为M/B;
令b的初值为1,且b=1,...,B;
(3b)构造第b段多径干扰抑制后的参考信号及其时延张成的子空间矩阵
Figure FDA0002276854510000025
Figure FDA0002276854510000031
其中,p为数字回波信号中需要消除的多径干扰的时延长度,T表示转置;
(3c)第b段数字回波信号
Figure FDA0002276854510000032
为:
Figure FDA0002276854510000033
(3d)根据所述第b段多径干扰抑制后的参考信号及其时延张成的子空间矩阵
Figure FDA0002276854510000034
构造第b段数字回波信号
Figure FDA0002276854510000035
的时域相消子空间投影矩阵wb,并根据所述时域相消子空间投影矩阵wb对第b段数字回波信号
Figure FDA0002276854510000036
进行时域干扰相消,得到时域干扰相消后的第b段数字回波信号
Figure FDA0002276854510000037
(3e)令b的值加1,并重复执行子步骤(3b)至(3d),得到时域干扰相消后的B段数字回波信号,将所述时域干扰相消后的B段数字回波信号组合成一路信号,并作为时域干扰相消后的回波信号
Figure FDA0002276854510000038
3.根据权利要求2所述的一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,其特征在于,子步骤(3d)具体包括:
构建最小优化目标函数:
Figure FDA0002276854510000039
其中,
Figure FDA00022768545100000310
表示求使得表达式·最小时的wb值;||·||2表示求二范数;
计算所述最小化目标函数的梯度:
Figure FDA00022768545100000311
其中,
Figure FDA00022768545100000312
表示求导操作;
当所述最小化目标函数的梯度为零时,求解得到第b段数字回波信号
Figure FDA0002276854510000041
的时域相消子空间投影矩阵
Figure FDA0002276854510000042
根据所述时域相消子空间投影矩阵wb对第b段数字回波信号
Figure FDA0002276854510000043
进行时域干扰相消,得到时域干扰相消后的第b段数字回波信号
Figure FDA0002276854510000044
4.根据权利要求1所述的一种外辐射源雷达参考信号含有多径干扰时的目标检测方法,其特征在于,步骤4具体包括如下子步骤:
(4a)分别对所述多径干扰抑制后的参考信号
Figure FDA0002276854510000045
所述时域干扰相消后的回波信号
Figure FDA0002276854510000046
进行快速傅里叶变换,分别得到多径干扰抑制后的频域参考信号
Figure FDA0002276854510000047
和时域干扰相消后的频域回波信号
Figure FDA0002276854510000048
其中,FFT{·}表示快速傅里叶变换操作;
(4b)对所述多径干扰抑制后的频域参考信号Sref(f)和时域干扰相消后的频域回波信号Ssur(f)进行脉冲压缩处理,得到脉冲压缩处理的频域输出S0(f):
Figure FDA0002276854510000049
其中,(·)*表示取共轭操作,sd(n)表示多径干扰抑制后的参考信号中的时域直达波信号,n(n)表示参考通道时域噪声,sechoj(n)表示时域干扰相消后的回波信号中的第j个时域目标回波信号,
Figure FDA00022768545100000410
表示回波通道时域噪声,Sd(f)表示多径干扰抑制后的参考信号中的频域直达波信号,n(f)表示参考通道频域噪声,Sechoj(f)表示时域干扰相消后的回波信号中的第j个频域目标回波信号,
Figure FDA0002276854510000051
表示回波通道频域噪声,j=1,2,...,m,m表示时域干扰相消后的回波信号中包含的目标回波信号的个数,Sn(f)表示脉冲压缩处理后总的频域噪声信号;
(4c)将所述脉冲压缩处理的频域输出S0(f)变换到时域,得到时域输出s0(t),并将所述时域输出s0(t)作为目标检测结果;所述时域输出s0(t)表示为:
Figure FDA0002276854510000052
其中,IFFT{·}表示快速逆傅里叶变换操作,sRj(n)表示脉冲压缩处理得到的第j个时域目标回波信号,sn(n)表示脉冲压缩处理后的时域噪声信号平台。
CN201710352530.6A 2017-05-18 2017-05-18 外辐射源雷达参考信号含有多径干扰时的目标检测方法 Active CN107153178B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710352530.6A CN107153178B (zh) 2017-05-18 2017-05-18 外辐射源雷达参考信号含有多径干扰时的目标检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710352530.6A CN107153178B (zh) 2017-05-18 2017-05-18 外辐射源雷达参考信号含有多径干扰时的目标检测方法

Publications (2)

Publication Number Publication Date
CN107153178A CN107153178A (zh) 2017-09-12
CN107153178B true CN107153178B (zh) 2020-04-10

Family

ID=59792868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710352530.6A Active CN107153178B (zh) 2017-05-18 2017-05-18 外辐射源雷达参考信号含有多径干扰时的目标检测方法

Country Status (1)

Country Link
CN (1) CN107153178B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108919206B (zh) * 2018-07-18 2022-03-04 武汉大学 一种基于子载波处理的外辐射源雷达极化滤波方法
CN109143184B (zh) * 2018-10-29 2022-05-24 北京理工大学 一种扫描雷达的双门限检测方法
CN109358318B (zh) * 2018-11-20 2023-05-05 南京理工大学 一种外辐射源雷达盲源分离提取目标回波与直达波的方法
CN110646769B (zh) * 2019-09-03 2021-07-20 武汉大学深圳研究院 一种适用于lte外辐射源雷达的时域杂波抑制方法
CN111650563B (zh) * 2020-06-15 2022-05-31 桂林电子科技大学 一种外辐射源雷达同频干扰时延和能量快速估计系统及方法
CN113433523B (zh) * 2021-06-17 2022-09-16 同济大学 一种基于双框架循环rls的汽车雷达干扰抑制方法
CN113419219B (zh) * 2021-06-18 2022-04-22 桂林电子科技大学 基于空域特征认知的外辐射源雷达同频干扰级联相消方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014022399A2 (en) * 2012-07-30 2014-02-06 Telecommunication Systems Inc. Navigation redirect using cmas emergency alerts
CN103401824B (zh) * 2013-07-22 2016-03-02 西安电子科技大学 基于修正牛顿法的频率选择性mimo系统空时盲均衡器方法
CN103630891B (zh) * 2013-12-03 2016-04-27 西安电子科技大学 利用gpu实现外辐射源雷达中估计目标来波方向的方法
CN104777466B (zh) * 2015-03-23 2017-05-03 河海大学 基于降维恒模盲均衡的无源双基地雷达信号处理方法
CN104749562B (zh) * 2015-03-23 2017-03-15 河海大学 基于循环常数模盲均衡的无源双基地雷达信号处理方法
CN106054143B (zh) * 2016-07-12 2019-02-19 西安电子科技大学 一种外辐射源雷达同频干扰的消除方法

Also Published As

Publication number Publication date
CN107153178A (zh) 2017-09-12

Similar Documents

Publication Publication Date Title
CN107153178B (zh) 外辐射源雷达参考信号含有多径干扰时的目标检测方法
Guan et al. Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter
CN103728594B (zh) 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法
CN103984676A (zh) 一种基于协方差矩阵重构的正交投影自适应波束形成方法
CN105527610B (zh) 基于分数阶时延估计的多天线联合优化杂波抑制方法
CN102955150B (zh) 基于天线主瓣占优强度约束的sar方位模糊抑制方法
CN106772253B (zh) 一种非均匀杂波环境下的雷达杂波抑制方法
CN105137437B (zh) 一种基于空域相位方差加权的目标检测方法
CN107561507B (zh) 一种外辐射源雷达的杂波对消方法
CN110646769B (zh) 一种适用于lte外辐射源雷达的时域杂波抑制方法
CN105302936A (zh) 基于相关计算和协方差矩阵重构的自适应波束形成方法
CN109143190B (zh) 一种零陷展宽的宽带稳健自适应波束形成方法
CN111693971B (zh) 一种用于弱目标检测的宽波束干扰抑制方法
CN103760529A (zh) 一种无源探测高效级联空时自适应处理方法
CN111045002A (zh) 基于TRT和SNuFFT的机动目标相参积累方法
Rabideau Clutter and jammer multipath cancellation in airborne adaptive radar
CN105044688B (zh) 基于迭代子空间跟踪算法的雷达稳健空时自适应处理方法
CN115372925A (zh) 一种基于深度学习的阵列稳健自适应波束形成方法
CN109799484B (zh) 一种外辐射源雷达系统多径抑制方法、系统及存储介质
CN103926570B (zh) 基于多通道rls的外辐射源雷达海杂波干扰抑制方法
CN110850421A (zh) 基于混响对称谱的空时自适应处理的水下目标检测方法
CN111044996A (zh) 一种基于降维近似消息传递的lfmcw雷达目标检测方法
CN105334502A (zh) 基于fir滤波器的雷达空时二维自适应处理方法
CN114152918A (zh) 基于压缩感知的抗间歇式主瓣干扰方法
Tian et al. Main-Lobe Interrupted Sampling Repeater Jamming Suppression Method based on Signal Reconstruction and Mismatched Filtering

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