CN107831480B - 弹载雷达和差通道稳健自适应杂波抑制方法 - Google Patents

弹载雷达和差通道稳健自适应杂波抑制方法 Download PDF

Info

Publication number
CN107831480B
CN107831480B CN201710954594.3A CN201710954594A CN107831480B CN 107831480 B CN107831480 B CN 107831480B CN 201710954594 A CN201710954594 A CN 201710954594A CN 107831480 B CN107831480 B CN 107831480B
Authority
CN
China
Prior art keywords
time
constraint
target
space
missile
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
CN201710954594.3A
Other languages
English (en)
Other versions
CN107831480A (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 CN201710954594.3A priority Critical patent/CN107831480B/zh
Publication of CN107831480A publication Critical patent/CN107831480A/zh
Application granted granted Critical
Publication of CN107831480B publication Critical patent/CN107831480B/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
    • G01S7/414Discriminating targets with respect to background clutter

Landscapes

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

Abstract

本发明公开了一种基于和差通道体制的稳健自适应杂波抑制方法,解决目标约束不确定和目标污染情况下传统弹载雷达和差体制面临的目标检测性能损失严重的问题。其实现步骤为:建立回波信号模型,得到空时快拍数据;用空时快拍数据估计协方差矩阵;设计幅相联合多点约束的空时二维和波束;设计联合零点和导数约束的空时差波束。本发明直接利用弹载雷达已有的和差通道输出信号,联合空域和时域信息并使用多种约束进行空时自适应处理。克服目标污染和目标约束的不确定性,具有较好的鲁棒性,做到了主瓣保形,提升了目标角度和多普勒参数估计的精度和稳健性,增强了动目标检测性能,用于弹载雷达领域。

Description

弹载雷达和差通道稳健自适应杂波抑制方法
技术领域
本发明属于雷达技术领域,主要涉及弹载雷达的杂波抑制,具体是一种弹载雷达和差通道稳健自适应杂波抑制方法,可用于动目标检测,提高角度和多普勒参数估计的性能。
背景技术
当弹载雷达处于下视工作状态时,杂波多普勒会扩散,导致微弱目标信号被杂波淹没。空时自适应处理(Space-Time Adaptive Processing,STAP)联合空间和时间二维信息能有效提高杂波抑制能力和动目标检测性能。传统STAP方法的角度估计取决于波束主瓣的指向,多普勒估计则由对应的检测多普勒通道给出,因此角度估计的精度受制于主瓣的宽度,多普勒估计的精度则受制于相干处理时间即频率分辨率。这些限制使得传统STAP方法的目标参数估计精度难以满足一些实际工程需要。
单脉冲测角技术已经广泛在跟踪雷达系统中用于实时角度参数估计,目前弹载雷达也多采用单脉冲天线。将经典的相控阵多波束STAP算法从阵元域推广到只有和差两个波束的波束域,即直接对弹载雷达的和差两个通道信号进行STAP处理的方法称之为和差STAP,利于工程实现,也能提高目标参数估计的精度。
R.D.Brown等在论文“STAP for clutter suppression with sum anddifference beams”(IEEE Trans.Aerosp.Electron.Syst.,vol.36,no.2,pp.634-646,Apr.2000.)中对和差STAP进行了介绍,其和差通道是在射频端形成的,即空间已经微波合成了和差方向图,弹载雷达通常都是这种情形。文中方法没有考虑相干干扰等非理想因素造成的自适应检测器性能损失,会引起主瓣畸变,目标临近自适应方向图凹口,导致目标检测性能恶化。
在STAP雷达中,存在协方差矩阵和目标导向矢量估计不准确的问题,当目标没有被准确约束保护时,会被当作干扰进行抑制,这会导致自适应空时响应主瓣的严重畸变。
实际中,STAP雷达还需要解决目标污染的问题,即当若干相距较近的目标出现在训练样本中时,真实目标外的相干干扰会引起目标污染,造成真实目标被抑制。
综上,传统和差STAP技术受相干干扰,协方差矩阵及导向矢量估计误差和目标污染问题的影响较大,会出现严重的目标检测性能下降。
发明内容
本发明的目的在于针对弹载雷达和差体制存在的目标检测性能损失严重的问题,提出一种受协方差矩阵及导向矢量估计误差影响较小的,能较好地克服目标污染并能精确估计目标参数的弹载雷达和差通道稳健自适应杂波抑制方法。
本发明是一种弹载雷达和差通道稳健自适应杂波抑制方法,其特征在于,包括有如下步骤:
(1)构造雷达回波数据:以前视弹载雷达为平台,基于微波和差通道的回波输出构造空(和差通道)时(相干脉冲)快拍数据,考虑发射K个相干脉冲,假定接收到L个距离单元的回波,则第l个距离单元的空时快拍数据表示为:
x(l)=xs(l)+xc(l)+xn(l)
其中xs(l)为第l个距离单元的运动目标信号,xc(l)为第l个距离单元的杂波回波信号,xn(l)为第l个距离单元的加性高斯白噪声;
(2)估计回波数据协方差矩阵:由L个距离单元的空时快拍数据估计协方差矩阵
Figure BDA0001433653650000024
Figure BDA0001433653650000021
(3)设计幅相联合多点约束的空时二维和波束:
3.1先在主瓣区域内约束目标(θ0,fd0)的临近角度和临近多普勒频率选取i个约束点(θi,fdi)。其中θ0是各约束点的方位角,fd0是各约束点的多普勒频率,θi是各约束点的方位角,fdi是各约束点的多普勒频率;
3.2然后对应地将约束目标和各约束点的空时导向矢量S(θi,fdi)作为列向量组成方向矩阵CΣ
3.3再利用与静态和波束幅度相位响应相匹配的复值响应矢量fΣ,得到幅相联合多点约束的和波束自适应处理器:
Figure BDA0001433653650000022
式中的上标H代表共扼转置运算,其中,和波束自适应权值wΣ表示为:
Figure BDA0001433653650000023
(4)设计零点和导数联合约束的时域差波束自适应处理器:
Figure BDA0001433653650000031
式中s(θ0,fd0)表示约束目标的空时导向矢量,
Figure BDA0001433653650000032
表示约束目标空时导向矢量的时域导数,时域差波束自适应权值wΔ表示为:
Figure BDA0001433653650000033
式中的上标-1代表求逆运算,其中CΔ为方向矩阵,fΔ为响应矢量。
本发明提供的稳健的自适应和差STAP方法,用以提高目标约束不确定和目标污染条件下的目标检测性能。
本发明与现有技术相比,具有以下优点:
(a)本发明提出的基于和差通道体制的稳健自适应杂波抑制技术,通过利用幅相联合多点约束、导数约束、零点约束使得自适应后的和差波束方向图在主瓣区域内逼近静态和差波束方向图的幅度和相位响应特性。
(b)相对于传统方法,本发明通过利用幅相联合多点约束和波束、导数和零点联合约束差波束克服了目标污染造成的性能损失,对于目标导向矢量估计误差和协方差矩阵估计误差具有一定的稳健性,改善了杂波抑制和动目标参数估计的性能。
(c)相对于传统方法,本发明通过将问题扩展到高维空间,增加了波束形成器的自由度,克服了多约束STAP方法自由度的损失问题,提高了杂波抑制和目标检测的性能。
附图说明
图1为本发明的实现总流程图;
图2为本发明的原理图;
图3为本发明所采用的前视弹载雷达的模型示意图;
图4为幅相联合多点约束方法在主瓣区域内约束目标的临近角度和临近多普勒选取的约束点分布示意图;
图5为传统MVDR-STAP、传统LCMV-STAP以及本发明所用的幅相联合多点约束STAP方法得到的自适应和波束方向图;
图6为传统MVDR-STAP、传统LCMV-STAP以及本发明所用的导数零点联合约束STAP方法得到的自适应时域差波束方向图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明。
实施例1
现有和差STAP方法在协方差矩阵和目标导向矢量估计不准确和目标污染的情况下,自适应空时响应主瓣会发生畸变,造成真实目标被抑制,目标检测性能受到损失。
针对这种现状,本发明展开了研究与创新,提出一种弹载雷达和差通道稳健自适应杂波抑制方法,该发明基于和差体制,因此能利用目前弹载雷达已有的和差通道数据直接进行后续空时处理,即联合空域信息和时域信息并利用多种约束进行杂波抑制和动目标检测。参照图1,弹载雷达和差通道稳健自适应杂波抑制方法包括有如下步骤:
(1)构造雷达回波数据:以前视弹载雷达为平台,基于微波和差通道的回波信号输出构造空时快拍数据,所述“空”指的是空间中的和差通道,“时”指的是时间上的K个相干脉冲,假定接收到仿真区域内L个距离单元的回波,则第l个距离单元的空时快拍数据表示为:
x(l)=xs(l)+xc(l)+xn(l)
其中xs(l)为第l个距离单元的运动目标信号,xc(l)为第l个距离单元的杂波回波信号,xn(l)为第l个距离单元的加性高斯白噪声。
(2)估计回波数据协方差矩阵:由L个距离单元的空时快拍数据估计回波协方差矩阵
Figure BDA0001433653650000043
Figure BDA0001433653650000041
式中L代表回波仿真区域内距离单元总数,l代表计算所用距离单元的序号,取值范围为[1,L]。
(3)设计幅相联合多点约束的空时二维和波束:
3.1先在主瓣区域内约束目标(θ0,fd0)的临近角度和临近多普勒频率选取i个约束点(θi,fdi)。其中θ0是各约束点的方位角,fd0是各约束点的多普勒频率,θi是各约束点的方位角,fdi是各约束点的多普勒频率;
3.2然后对应地将约束目标和各约束点的空时导向矢量S(θi,fdi)作为列向量组成方向矩阵CΣ
3.3再利用与静态和波束幅度相位响应相匹配的复值响应矢量fΣ,得到幅相联合多点约束的和波束自适应处理器:
Figure BDA0001433653650000042
本发明的和波束自适应处理器通过对步骤3.1中主瓣区域内的约束目标和所取的约束点进行幅度相位联合约束,使得自适应后的和波束方向图在主瓣区域内逼近静态和波束方向图的幅度和相位响应特性,因此克服了目标污染的情况下的主瓣畸变,做到了主瓣保形,提升了自适应和波束的稳健型。同时通过对自适应处理器的输出进行最小化,达到了杂波抑制的目的。
式中的上标H代表共扼转置运算,其中,和波束自适应权值wΣ表示为:
Figure BDA0001433653650000051
(4)设计零点和导数联合约束的时域差波束自适应处理器:
Figure BDA0001433653650000052
式中s(θ0,fd0)表示约束目标的空时导向矢量,
Figure BDA0001433653650000053
表示约束目标空时导向矢量的时域导数,时域差波束自适应权值wΔ表示为:
Figure BDA0001433653650000054
式中的上标-1代表求逆运算,其中CΔ为方向矩阵,fΔ为响应矢量。
相对于传统方法,本发明采用了幅相联合多点约束的空时二维和波束以及零点和导数联合约束的时域差波束,克服了目标污染造成的性能损失,对于目标导向矢量估计误差和协方差矩阵估计误差具有一定的稳健性,改善了动目标检测的性能。
本发明的技术原理是,直接对弹载雷达现有的微波和差通道的模拟回波信号输出进行采样,然后将得到的时域数字信号转换到频域上,最后通过空时自适应处理器得到目标信号输出,完成目标信号检测。参照图2,同时对图2的和通道输出模拟信号Σ和差通道输出模拟信号Δ进行A/D采样,完成模拟信号到数字信号的转换。之后通过对数字信号施加K个延时并通过FFT变化将信号转换到频域,最后在频域和空域上对回波信号进行空时自适应处理,完成杂波抑制和目标检测。
实施例2
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1,本发明步骤3.2中的方向矩阵CΣ是由步骤3.1中约束目标和各约束点的空时导向矢量S(θi,fdi)作为列向量组成的,表示为:
Figure BDA0001433653650000055
式中,方向矩阵CΣ的维度为2K*5,
参照图4,图4为幅相联合多点约束方法在主瓣区域内约束目标的临近角度和临近多普勒选取的约束点分布示意图。在主瓣区域内约束目标
Figure BDA0001433653650000061
的临近多普勒频率上选取了2个约束点,其中(θ1,fd1)见图中的①,(θ2,fd2)见图中的②,临近角度上也选取了2个约束点,其中(θ3,fd3)见图中的③,(θ4,fd4)见图中的④。
通过在多普勒维和角度维选取多个紧邻的约束点构造一个约束区域,并利用副相联合约束,使得即使约束区域内出现目标污染,主瓣也不会出现严重畸变。
实施例3
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1-2,本发明步骤3.3中所述的与静态和波束幅度相位响应相匹配的复值响应矢量fΣ表示为:
Figure BDA0001433653650000062
其中上标T代表转置运算,α代表复值响应矢量元素的幅度,β代表复值响应矢量元素的相位:
αi=|sHi,fdi)s(θ0,fd0)|=[aΣ Hi)aΣ0)]*[bH(fdi)b(fd0)]+[aΔ Hi)aΔ0)]*[bH(fdi)b(fd0)]
βi=∠{sHi,fdi)s(θ0,fd0)}
式中αi为第i个约束点对应的复值响应矢量元素的幅度,βi为第i个约束点对应的复值响应矢量元素的相位。aΣi)为第i个约束点对应的和方向图取值,aΔi)为第i个约束点对应的差方向图取值,b(fdi)为第i个约束点对应的时域导向矢量。
本发明通过利用幅相联合多点约束,使得自适应后的和波束方向图在主瓣区域内逼近静态和波束方向图的幅度和相位响应特性。
实施例4
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1-3,本发明步骤4中的时域差波束自适应处理器的自适应权值表达式中的方向矩阵CΔ和响应矢量fΔ分别表示为:
Figure BDA0001433653650000063
fΔ=[1,0]T
其中
Figure BDA0001433653650000071
是对估计的目标空时导向矢量s(θ0,fd0)在时域的3dB主瓣区域内进行求导得到的导数:
Figure BDA0001433653650000072
联合利用零点和导数约束的目的是形成多普勒维的方向图凹口,得到时域差波束。
下面给出一个更加详细完整的例子,对本发明进一步说明。
实施例5
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1-4,
参照图1,本发明的实现步骤包括有如下步骤:
步骤1:构造雷达回波数据
以前视弹载雷达为平台,基于微波和差通道的模拟回波输出构造空(和差通道)时(相干脉冲)快拍数据,考虑发射K个相干脉冲,假定接收到L个距离单元的回波,则第l个距离单元的空时快拍数据表示为:
x(l)=xs(l)+xc(l)+xn(l)
其中xs(l)为第l个距离单元的运动目标信号,xc(l)为第l个距离单元的杂波回波信号,xn(l)为第l个距离单元的加性高斯白噪声;
运动目标信号
Figure BDA0001433653650000073
可以表示为:
Figure BDA0001433653650000074
其中上标T代表转置运算,σs代表运动目标信号的复系数,θs和fds分别代表目标的空间锥角和多普勒频率,
Figure BDA0001433653650000075
代表克罗内克(kronecker)积,aΣs)和aΔs)分别代表和通道与差通道在目标处的输出,
Figure BDA0001433653650000076
则为时域导向矢量,具体可以表示为:
Figure BDA0001433653650000077
fr代表脉冲重复频率(PRF),fds代表动目标的多普勒频率:
Figure BDA0001433653650000078
υp代表导弹平台的飞行速度,vs代表动目标的速度,
Figure BDA0001433653650000081
为目标与导弹飞行方向的锥角。
地杂波回波
Figure BDA0001433653650000082
由等距离环内的许多散射单元的回波相干叠加而成。即表示为:
Figure BDA0001433653650000083
其中Na为距离模糊数,Nc为一个距离环内独立散射单元的个数,σm,q(l)为第m次模糊距离环的第q个散射单元的复系数,其幅度分布由地面特性决定,比如草地的散射复系数分布为瑞利分布,θm,q和fdm,q分别为相应散射单元的角度和多普勒频率,对于前视弹载雷达,地面散射单元的多普勒、角度的耦合关系满足:
Figure BDA0001433653650000084
Figure BDA0001433653650000085
为相应的地面散射单元与导弹飞行方向的锥角。
步骤2:估计回波数据协方差矩阵
由L个距离单元的空时快拍数据x(l)估计回波协方差矩阵
Figure BDA0001433653650000086
Figure BDA0001433653650000087
式中xH(l)的上标H代表共轭转置运算符。
步骤3:设计幅相联合多点约束的空时二维和波束
如图4所示,先在主瓣区域内约束目标(θ0,fd0)的临近角度和临近多普勒选取i个约束点(θi,fdi),然后对应地将约束目标和各约束点的空时导向矢量S(θi,fdi)作为列向量组成方向矩阵CΣ
Figure BDA0001433653650000088
第i个约束点的空时导向矢量
Figure BDA0001433653650000089
表示为:
Figure BDA00014336536500000810
约束目标的空时导向矢量
Figure BDA00014336536500000811
表示为:
Figure BDA00014336536500000812
θ0为估计的目标与线阵轴线的锥角,fd0是估计的动目标多普勒频率。
Figure BDA00014336536500000813
vs0为估计的目标径向速度,
Figure BDA0001433653650000091
为估计的目标与导弹飞行方向的锥角。
不同于传统LCMV方法采用主波束内全1约束,本发明所用和波束的响应矢量是利用约束目标的空时导向矢量s(θ0,fd0)得到的与静态和波束幅度相位响应相匹配的复值响应矢量:
Figure BDA0001433653650000092
式中:
αi=|sHi,fdi)s(θ0,fd0)|=[aΣ Hi)aΣ0)]*[bH(fdi)b(fd0)]+[aΔ Hi)aΔ0)]*[bH(fdi)b(fd0)]
βi=∠{sHi,fdi)s(θ0,fd0)}。
幅相联合多点约束的空时二维和波束对应的自适应处理器表示为:
Figure BDA0001433653650000093
自适应权值wΣ表示为
Figure BDA0001433653650000094
步骤4:设计零点和导数联合约束的时域差波束
对估计的目标空时导向矢量s(θ0,fd0)在时域的3dB主瓣区域内进行求导得到:
Figure BDA0001433653650000095
联合零点和导数约束的差波束自适应处理器为:
Figure BDA0001433653650000096
求得时域差波束的自适应权值
Figure BDA0001433653650000097
为:
Figure BDA0001433653650000098
式中
Figure BDA0001433653650000099
fΔ=[1,0]T
鉴于增加的约束条件会消耗波束形成器的自由度(DOF)进而会削减抑制杂波及干扰的能力,因此我们将以上差波束形成器扩展到高维来增加自由度,进而增强抑制杂波和干扰的能力:
Figure BDA0001433653650000101
其中
Figure BDA0001433653650000102
为权值矩阵,并且可以从代价函数看出上述求解最优权的问题对于权值矩阵WΔ是线性的。虽然上式增加了自由度但对权值矩阵WΔ秩为1的约束会导致非凸优化的问题。我们可以采用半正定松弛的优化方法去解决,也就是可以将秩为1的约束直接丢弃将非凸优化问题转化为凸优化问题:
Figure BDA0001433653650000103
相对于传统方法,本发明通过将问题扩展到高维空间,增加了波束形成器的自由度,克服了多约束STAP方法自由度的损失问题,提高了杂波抑制和目标检测的性能。
本发明的效果可以通过以下仿真实验做进一步说明。
实施例6
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1-5
一.实验环境
本发明所用的前视弹载雷达模型如图3所示,坐标系原点O为导弹平台在水平面的投影点,仿真中采用均匀线阵,天线阵元数N为8,d为阵元间距并等于波长的一半,y轴为导弹运动速度υp的方向,H为导弹的高度,αs为目标方位角,φS为目标的俯仰角,R为目标与导弹的距离,σs为目标信号的复系数,vs为目标的径向速度。
本发明的实例所用的各种参数如表1,另外假设仿真环境中空域和时域均为无方向性的高斯白噪声,两个角度多普勒域的临近目标被包含在训练样本中,两个目标的参数基本保持一致,分别位于第145号距离门和第150号距离门,其中待检测期望目标在第150号距离门,相应的保护单元为从148到152号距离门,145号距离门的目标则被视为相干干扰。
表1 STAP雷达参数
Figure BDA0001433653650000104
Figure BDA0001433653650000111
二.仿真内容
在所述仿真条件下,进行如下实验。
仿真实验1,比较传统MVDR-STAP、传统LCMV-STAP以及本发明所用的幅相联合多点约束STAP方法得到的自适应和波束方向图,仿真结果如图5。
仿真实验2,比较传统MVDR-STAP、传统LCMV-STAP以及本发明所用的导数零点联合约束STAP方法得到的自适应时域差波束方向图,仿真结果如图6。
三.结果分析
在实验1中,图5为空时自适应和波束方向图的结果比较。图5(a)是采用MVDR-STAP方法得到的和波束方向图,由图5(a)可以看出传统MVDR-STAP方法对目标污染较为敏感,方向图主瓣已经畸变,导致目标在角度多普勒域上的位置(十字所在位置)非常靠近方向图的凹口,这意味着MVDR-STAP方法误把真实目标当作干扰并对其进行了抑制。图5(b)是采用LCMV-STAP方法得到的和波束方向图,由图5(b)可以看出传统LCMV-STAP方法的方向图主瓣严重失真,目标依然临近凹口,目标被抑制。这是由于其采用的全1响应矢量与静态和波束存在失配问题,这种情况会导致目标检测性能的下降和参数估计精度的严重恶化。图5(c)是采用本发明所用方法得到的和波束方向图,由图5(c)则可以看出本发明所用的幅相联合多点约束STAP方法好于传统方法,由于采用与静态和波束相匹配的响应矢量,该方法能够实现方向图主瓣保形进而保证目标信号的无失真输出。
实施例7
弹载雷达和差通道稳健自适应杂波抑制方法同实施例1-5,仿真条件和仿真内容同实施例6
在实验2中,图6为空时自适应时域差波束方向图的结果比较。图6(a)是采用MVDR-STAP方法得到的时域差波束方向图,由图6(a)可以看出传统MVDR-STAP的时域差波束方向图凹口的位置被自适应地固定在了相干干扰的位置,这将导致严重的多普勒估计性能的下降。图6(b)是采用LCMV-STAP方法得到的时域差波束方向图,由图6(b)可以看出传统LCMV-STAP方法的方向图相对比较模糊,主瓣失真,差方向图的凹口也相对较宽,同样会引起多普勒估计性能的损失。图6(c)是采用导数和零点联合约束的方法得到的时域差波束方向图,由图6(c)可以看出本发明所用的导数和零点联合约束的方法好于传统方法,该方法能够在波束指向处形成清晰,较窄且较深的时域差方向图凹口,做到了主瓣保形并保证了较好的多普勒估计性能。
本发明利用幅相联合多点约束、导数约束、零点约束使得自适应和差波束方向图在主瓣区域内逼近静态和差波束方向图,实现了和差波束主瓣保形。而差波束设计中的导数约束因为不受阵列几何结构的制约,所以可以用于任何几何结构的阵列包括对称或者非对称阵列,另外本发明可以通过合理的零点约束实现角度多普勒解耦合。仿真实验也证明本发明能够较好地克服目标污染,对于目标导向矢量估计误差和协方差矩阵估计误差具有一定的稳健性,在目标参数估计精度和杂波抑制方面有着较大的提升。
简而言之,本发明公开的一种基于和差通道体制的稳健自适应杂波抑制方法,主要解决目标约束不确定和目标污染情况下传统弹载雷达和差体制面临的目标检测性能损失严重的问题。其实现步骤为:(1)建立回波信号模型,得到空时快拍数据;(2)用空时快拍数据估计协方差矩阵;(3)设计幅相联合多点约束的空时二维和波束;(4)设计零点和导数联合约束的空时差波束。
本发明用于弹载雷达领域,直接利用弹载雷达已有的和差通道输出信号,联合空域和时域信息并使用多种约束进行空时自适应处理。
本发明对于克服目标污染和目标约束不确定具有较好的鲁棒性,做到了主瓣保形,提升了目标角度和多普勒参数估计的精度和稳健性,增强了动目标检测性能。

Claims (4)

1.一种弹载雷达和差通道稳健自适应杂波抑制方法,其特征在于,包括有如下步骤:
(1)构造弹载雷达和差通道回波数据:以前视弹载雷达为平台,基于微波和差通道的回波输出构造空时快拍数据,考虑发射K个相干脉冲,假定接收到仿真区域内L个距离单元的回波,则第l个距离单元的空时快拍数据表示为:
x(l)=xs(l)+xc(l)+xn(l)
其中xs(l)为第l个距离单元的运动目标信号,xc(l)为第l个距离单元的杂波回波信号,xn(l)为第l个距离单元的加性高斯白噪声;
(2)估计和差通道回波数据协方差矩阵:由L个距离单元的空时快拍数据估计协方差矩阵
Figure FDA0002790298790000011
Figure FDA0002790298790000012
(3)设计幅相联合多点约束的空时二维和波束:
3.1先在主瓣区域内约束目标(θ0,fd0)的临近角度和临近多普勒频率选取i个约束点(θi,fdi),其中θ0是各约束点的方位角,fd0是各约束点的多普勒频率,θi是各约束点的方位角,fdi是各约束点的多普勒频率;
3.2然后对应地将约束目标和各约束点的和差体制空时导向矢量S(θi,fdi)作为列向量组成方向矩阵CΣ
3.3再利用与静态和波束幅度相位响应相匹配的复值响应矢量fΣ,得到幅相联合多点约束的和波束自适应处理器:
Figure FDA0002790298790000013
式中的上标H代表共扼转置运算,其中,和波束自适应权值wΣ表示为:
Figure FDA0002790298790000014
(4)设计零点和导数联合约束的时域差波束自适应处理器:
Figure FDA0002790298790000021
式中s(θ0,fd0)表示约束目标的空时导向矢量,
Figure FDA0002790298790000022
表示约束目标空时导向矢量的时域导数,时域差波束自适应权值wΔ表示为:
Figure FDA0002790298790000023
式中的上标-1代表求逆运算,其中CΔ为方向矩阵,fΔ为响应矢量。
2.根据权利要求1所述的弹载雷达和差通道稳健自适应杂波抑制方法,其特征在于,步骤3.2中的方向矩阵CΣ是由步骤3.1中约束目标和各约束点的和差体制空时导向矢量S(θi,fdi)作为列向量组成的,表示为:
Figure FDA0002790298790000024
3.根据权利要求1所述的弹载雷达和差通道稳健自适应杂波抑制方法,其特征在于,步骤3.3中所述的与静态和波束幅度相位响应相匹配的复值响应矢量fΣ表示为:
Figure FDA0002790298790000025
其中上标T代表转置运算,α代表复值响应矢量各元素的幅度,β代表复值响应矢量各元素的相位:
αi=|sHi,fdi)s(θ0,fd0)|=[aΣ Hi)aΣ0)]*[bH(fdi)b(fd0)]+[aΔ Hi)aΔ0)]*[bH(fdi)b(fd0)]
βi=∠{sHi,fdi)s(θ0,fd0)};
式中αi为第i个约束点对应的复值响应矢量元素的幅度,βi为第i个约束点对应的复值响应矢量元素的相位,aΣi)为第i个约束点对应的和方向图取值,aΔi)为第i个约束点对应的差方向图取值,b(fdi)为第i个约束点对应的时域导向矢量。
4.根据权利要求1所述的弹载雷达和差通道稳健自适应杂波抑制方法,其特征在于,步骤4中的时域差波束自适应处理器的自适应权值表达式中的方向矩阵CΔ和响应矢量fΔ分别表示为:
Figure FDA0002790298790000031
fΔ=[1,0]T
其中
Figure FDA0002790298790000033
是对估计的目标和差体制空时导向矢量s(θ0,fd0)在时域的3dB主瓣区域内进行求导得到的导数:
Figure FDA0002790298790000032
θ=θ0,fd=fd0
联合利用零点和导数约束的目的是形成多普勒维的方向图凹口,得到时域差波束。
CN201710954594.3A 2017-10-13 2017-10-13 弹载雷达和差通道稳健自适应杂波抑制方法 Active CN107831480B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710954594.3A CN107831480B (zh) 2017-10-13 2017-10-13 弹载雷达和差通道稳健自适应杂波抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710954594.3A CN107831480B (zh) 2017-10-13 2017-10-13 弹载雷达和差通道稳健自适应杂波抑制方法

Publications (2)

Publication Number Publication Date
CN107831480A CN107831480A (zh) 2018-03-23
CN107831480B true CN107831480B (zh) 2021-06-01

Family

ID=61648168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710954594.3A Active CN107831480B (zh) 2017-10-13 2017-10-13 弹载雷达和差通道稳健自适应杂波抑制方法

Country Status (1)

Country Link
CN (1) CN107831480B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108563144B (zh) * 2018-03-30 2021-06-29 西安电子科技大学 一种弹载雷达信号处理半实物仿真测试系统
CN108919205B (zh) * 2018-07-12 2022-07-08 中国船舶重工集团公司第七二四研究所 基于辅助阵元的频率分集阵列欺骗干扰抑制方法
CN109061594B (zh) * 2018-08-03 2022-10-28 中国航空工业集团公司雷华电子技术研究所 一种强杂波环境下基于子阵空时自适应处理的测角方法
CN109188389B (zh) * 2018-10-16 2023-03-28 哈尔滨工业大学 超视距多基被动雷达中的解时差测量模糊的方法
CN109361443B (zh) * 2018-11-21 2022-04-19 北京遥感设备研究所 一种自适应数字波束形成器及形成方法
CN111352101B (zh) * 2018-12-20 2023-08-18 中国人民解放军空军预警学院 一种相控阵机载雷达空时二维数字差通道形成方法
CN110596692B (zh) * 2019-08-19 2022-10-14 电子科技大学 基于联合约束的自适应单脉冲测向方法
CN110346764B (zh) * 2019-08-21 2022-03-29 上海无线电设备研究所 一种mimo雷达目标距离角度解耦方法
CN111257863B (zh) * 2019-12-26 2023-06-06 电子科技大学 一种高精度多点线性约束的自适应单脉冲测向方法
CN111856410B (zh) * 2020-07-21 2022-08-05 中国人民解放军国防科技大学 恒模约束下机载mimo雷达的稳健发射接收联合优化方法
CN111999727B (zh) * 2020-07-29 2022-07-05 中国人民解放军国防科技大学 一种基于主瓣幅度响应控制的机载频率分集阵列雷达快速运动目标检测方法
CN113885019B (zh) * 2021-09-24 2024-05-14 西安电子科技大学 基于和差天线体制雷达空时自适应处理的测角方法
CN114137483A (zh) * 2021-11-03 2022-03-04 广州辰创科技发展有限公司 一维相扫雷达自适应干扰抑制方法、介质及一维相扫雷达

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879767A (zh) * 2012-06-28 2013-01-16 北京理工大学 一种用于空时自适应处理的干扰目标检测方法
CN104808180A (zh) * 2015-04-09 2015-07-29 大连大学 杂波环境下mimo雷达稳健波形优化方法
CN105137409A (zh) * 2015-07-24 2015-12-09 西安电子科技大学 基于幅相约束的目标信号稳健空时自适应处理方法
CN105445703A (zh) * 2015-11-27 2016-03-30 西安电子科技大学 一种机载雷达空时回波数据的两级空时自适应处理方法
CN106383342A (zh) * 2016-09-09 2017-02-08 深圳大学 一种基于存在测量误差的阵列流形先验知识的稳健stap方法
CN106546964A (zh) * 2016-10-26 2017-03-29 深圳大学 一种机载雷达杂波秩估计方法及其系统
CN107229048A (zh) * 2017-06-06 2017-10-03 电子科技大学 一种高分宽幅sar动目标速度估计与成像方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9772402B2 (en) * 2014-06-09 2017-09-26 Src, Inc. Multiplatform GMTI radar with adaptive clutter suppression

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879767A (zh) * 2012-06-28 2013-01-16 北京理工大学 一种用于空时自适应处理的干扰目标检测方法
CN104808180A (zh) * 2015-04-09 2015-07-29 大连大学 杂波环境下mimo雷达稳健波形优化方法
CN105137409A (zh) * 2015-07-24 2015-12-09 西安电子科技大学 基于幅相约束的目标信号稳健空时自适应处理方法
CN105445703A (zh) * 2015-11-27 2016-03-30 西安电子科技大学 一种机载雷达空时回波数据的两级空时自适应处理方法
CN106383342A (zh) * 2016-09-09 2017-02-08 深圳大学 一种基于存在测量误差的阵列流形先验知识的稳健stap方法
CN106546964A (zh) * 2016-10-26 2017-03-29 深圳大学 一种机载雷达杂波秩估计方法及其系统
CN107229048A (zh) * 2017-06-06 2017-10-03 电子科技大学 一种高分宽幅sar动目标速度估计与成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
机载雷达超低空多径目标稳健STAP方法;冯阳 等;《系统工程与电子技术》;20170731;第39卷(第7期);全文 *
静态和差波束匹配的空时自适应处理方法;刘斌 等;《西安电子科技大学学报(自然科学版)》;20161204;第44卷(第3期);第1节信号模型与问题描述、第2节静态和差波束匹配的STAP方法 *

Also Published As

Publication number Publication date
CN107831480A (zh) 2018-03-23

Similar Documents

Publication Publication Date Title
CN107831480B (zh) 弹载雷达和差通道稳健自适应杂波抑制方法
CN109581352B (zh) 一种基于毫米波雷达的超分辨测角系统
CN111948619B (zh) 一种多类型外辐射源照射下目标被动协同探测方法及系统
CN110673086A (zh) 一种基于数字阵列雷达的二维角度超分辨方法
CN111007503B (zh) 基于频率谱精确定位的运动目标聚焦和定位方法及系统
Tan et al. An efficient range-Doppler domain ISAR imaging approach for rapidly spinning targets
CN112834980A (zh) 一种基于传播算子的涡旋电磁波超分辨测向方法
CN110133646B (zh) 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法
CN105824016B (zh) 运动平台雷达检测超低空目标的稳健空时自适应处理方法
Liu et al. Enhanced forward-looking missile-borne bistatic SAR imaging with electromagnetic vortex
Zhu et al. Robust moving targets detection and velocity estimation using multi-channel and multi-look SAR images
CN109597034B (zh) 一种基于欧几里得距离的空时自适应处理方法
Xiong et al. Nonstationary clutter suppression based on four dimensional clutter spectrum for airborne radar with conformal array
CN111812607A (zh) 基于波束空间的米波mimo雷达低仰角估计方法
Yili et al. Short-range clutter suppression method combining oblique projection and interpolation in airborne CFA radar
Li et al. Range-dependent clutter cancellation method in bistatic MIMO-STAP radars
De et al. Angle estimation using modified subarray level monopulse ratio algorithm and s-curve in digital phased array radar
CN114779199A (zh) 端射阵机载雷达自适应互耦补偿与杂波抑制方法
Garry et al. Framework and results on passive ISAR
Yuan et al. 3DRUDAT: 3D robust unambiguous Doppler beam sharpening using adaptive threshold for forward-looking region
Xiong et al. Space time adaptive processing for airborne MIMO radar based on space time sampling matrix
Zhu et al. Application of frequency diverse array to resolve range ambiguity for SAR imaging
Xue et al. An applied frequency scaling algorithm based on local stretch factor for near-field miniature millimeter-wave radar imaging
CN111505597A (zh) 一种消除雷达波束赋形栅瓣的方法及天线
Chen et al. Space-time adaptive monopulse based on space-time uniform constraint

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