基于多输入多输出合成孔径雷达的二维混合基线抗欺骗干扰方法
2019-04-04张福博陈龙永卜祥玺周思言
曾 铮 张福博 陈龙永* 卜祥玺 周思言
①(中国科学院电子学研究所微波成像技术重点实验室 北京 100190)
②(中国科学院大学 北京 100049)
1 引言
合成孔径雷达(Synthetic Aperture Radar,SAR)是主动式微波成像雷达,能全天候、全天时、实时获取大地域的地面图像,具有一定的植被和地面穿透能力,对资源勘察和军事侦察有重大意义,已广泛应用于民用和军事领域。现代战争中,电子干扰已大范围普及,对SAR的欺骗干扰可以在没有目标区域产生虚假目标,使SAR判断失误,从而掩护真实目标[1-3]。为了使SAR能够真正应用于军事领域,对SAR的抗欺骗干扰研究在当前是必要的和紧迫的。从干扰性质上来分,合成孔径雷达干扰主要分为压制性干扰和欺骗性干扰。当前阶段,压制性干扰已有较好的对抗措施[4,5],而欺骗干扰技术在信号处理增益上有着明显的优势,且针对SAR欺骗干扰的抗干扰措施还非常有限[6,7],因此研究SAR的抗欺骗性干扰具有重大意义。
目前,国内外的学者针对SAR的抗欺骗性干扰问题开展了许多研究工作。文献[8]提出根据欺骗干扰回波和真实目标回波方向不同的原理,利用距离向多孔径,对接收回波使用距离向数字波束形成(Digital Beam Forming, DBF)技术,通过谱估计得到欺骗干扰方向,然后计算加权系数,在干扰方向形成零陷,实现空域滤波,从而达到抗欺骗干扰的目的,但此方法要求系统有足够高的自由度和角分辨精度,天线设计难度较大,系统复杂度较高。文献[9]提出利用脉冲随机相位波形调制技术来对抗有源假目标干扰,并提出一种称为“惩罚函数”的方法实现对有源假目标干扰抑制。该方法实质上属于脉冲捷变体制的抗干扰方法,要求在发射信号带宽不变的前提下,对每个发射信号波形都要进行不同的调制操作,增加了信号调制和收发的难度,同时也大大增加了成像处理的复杂性,在现阶段实施起来有相当难度。文献[10]中提出一种新的脉内混沌调相-脉间随机编码的超混沌信号,增强抗欺骗干扰的能力。但实际工程中,混沌信号产生困难,对分布目标成像效果差,在一定程度上牺牲了SAR的成像性能。文献[11]中提出利用双孔径天线获取的两幅合成孔径雷达图像的相位差来检测欺骗干扰的新方法,但是该方法只利用了交轨干涉对欺骗干扰目标进行检测并消除,但消除效果有限。文献[12]中分析了双通道对消抑制合成孔径雷达欺骗干扰带来的真实目标的损失问题,并给出了目标损失函数和参数分析,这也正是双通道抗干扰所存在的缺陷。文献[13]中提出使用小卫星分布式合成孔径雷达系统的多通道干扰抑制技术实现抗干扰。但是该方法需要多颗卫星,保持星间同步,同时需要对编队构形优化,系统代价较大。文献[14]提出利用真假目标时频分布不同实现抗欺骗干扰目标重建,并能够实现对目标形状的提取,但当虚假目标与真实目标相距较近时,该方法不能有效抗干扰。
根据收发端的天线数目,可将雷达系统分成单发单收、单发多收及多发多收体制。随着天线数目的增加,系统自由度逐步提高,雷达性能稳健增长。然而,受平台尺寸与功率孔径积限制,不能单纯依靠提高接收天线数目来提升雷达的性能。在此背景之下,MIMO-SAR就是一种通过多发射天线多接收天线来获取丰富系统自由度的先进合成孔径雷达系统,它使用多个孔径同时发射信号,多个孔径接收回波,极大地提高了系统自由度[15],是当前SAR领域研究的热点,其在抗干扰问题上也具有相应的潜力。本文针对MIMO-SAR和欺骗干扰的特征,提出先分别获取不同干涉情况下的相位差并判别是否有欺骗干扰,然后进行相位补偿后相干积累,从而对干扰目标进行抑制,实现抗欺骗干扰的目的。
本文的结构安排如下:第2节主要介绍MIMOSAR抗欺骗干扰模型和算法原理,第3节为仿真实验验证,第4节对本文内容作总结。
2 基于MIMO-SAR的2维混合基线抗欺骗干扰算法
2.1 基于MIMO-SAR的编码方案
本文所提MIMO-SAR系统为4发4收系统,发射4路信号,先构建一个发射矩阵如式(1)所示
其中,f0为发射信号载频,Kr为信号的线性调频率,g(τ) 为门函数,τ为快时间。
则这4路发射信号分别为
其中,ω(η) 为 方位包络,fac为第2路发射信号的多普勒载频,η为慢时间,且有如式(3)所示关系
因此,第n个接收天线的回波可表示为
将式(4)变换至距离-多普勒域,可得
由式(5)可见,来自于不同发射天线的信号在距离多普勒域中处于不同的多普勒频率中心(见图1),因此可在距离多普勒域中采用带通滤波方法,有效分离多路发射天线的信号。
2.2 基于2维混合基线的欺骗干扰目标的判别
目前,SAR的2维假目标欺骗干扰方式主要是距离向延时和方位向多普勒频率调制联合形成[16-19],其中距离向延时欺骗干扰可在距离向形成假目标,方位向多普勒频率调制可在方位向形成假目标,两者联合可形成2维假目标。针对该类欺骗干扰方式,本文提出基于MIMO-SAR的2维混合基线抗欺骗干扰方法。本方法的主要原理如下:
(1) 多维信息识别欺骗干扰
真实信号的回波和干扰信号存在着某些特征上的区别,例如相位特征上的区别,这些区别因为延时或者多普勒频率调制被隐藏,但如果能够提取多维信息,如相位信息等,就可以发现真实信号和干扰信号的区别,来区别是否有干扰,这就有了SAR抗欺骗干扰的前提;
(2) 相位补偿抑制干扰目标
多通道相干积累的必要条件是通道间相位同步,否则为非相干积累,这一原理可以在SAR抗欺骗干扰上得到应用。本文可以通过相位补偿等手段,使得真实目标的相位差为0,满足相干积累,而干扰目标的相位差不为0,则不能相干积累。这样本文就能在判别干扰的前提下实现对干扰的抑制。
2.2.1 距离向欺骗干扰的判别
本节首先介绍距离向欺骗干扰产生的机理,然后介绍基于MIMO-SAR的距离向欺骗干扰的判别。首先,SAR几何模型图如图2所示,X轴为方位向,Y轴为距离向,Z轴为高度向。平台高度为h,飞行速度为Vr。 在t=0时刻,SAR经过坐标(0,0,h)处。A点为一真实点目标,以A点为干扰源,欲在B点位置处产生一个假目标。
产生机理如下,假设SAR只存在单一收发共用天线,t=0时,天线位于APC1点,坐标为(0, 0,h)。
设SAR发射为线性调频信号,同式(1),位于A点处的干扰源接收到的回波信号表达式为
式中,τ为快时间,η为慢时间,c 为光速,R(η)为瞬时斜距。在A点处的干扰源接收到上述回波信号后,干扰机将接收信号存储,并延时转发,转发后的回波信号表达式为
t1为延时转发的时间,如此,SAR接收到的回波信号比真实信号延迟了t1时间,B点处即为假目标出现的位置,设A点和B点之间距离为d,则
图1 多普勒频谱示意图Fig.1 Diagrammatic sketch of Doppler spectrum
针对上述欺骗干扰,基于MIMO-SAR的交轨干涉抗距离向延时欺骗干扰的机理如下:取4发4收的MIMO-SAR系统中的一对天线为例,如图2,t=0时,天线1位于APC1(0, 0,h),天线2位于APC2(0,s,h),交轨干涉的干涉基线长即为s。在B(B′)位置处,如果该点目标是真实点目标,本文称之为B点,如果该点目标是欺骗干扰所产生的假目标,本文称之为B′点,现分别比较A点和B(B′)点干涉相位差。
首先计算A点干涉相位,A点相对于通道1的瞬时斜距为
图2 距离向抗欺骗干扰几何模型Fig.2 Geometric model of range anti-deception jamming
其中,R0为最近斜距,z为雷达距离地面的高度。A点相对于通道2的瞬时斜距为
再计算真实目标B点的干涉相位,B点相对于通道1的瞬时斜距为
B点相对于通道2的瞬时斜距为
于是,B点的干涉相位为
所以,由式(11)和式(14)对比可以看出,如果是真实目标B点,则其交轨干涉相位和参考点A点是不一致的。而对于距离延时欺骗干扰目标B′点来说,式(12)和式(13)中,A点与B′点之间的距离d是通过快时间 τ的变化生成的,在计算瞬时斜距时,式(12)和式(13)应当与式(9)和式(10)相同。也就是说对A点做距离延时所生成的假目标B′点干涉相位和A点一致,它们的视角也是相同的。因此,对于距离向的某位置处未知真假点目标而言,本文可以根据干涉相位与参考点A点是否一致来判别,一致的是B′点假目标,不一致的是B点真实目标。
2.2.2 方位向欺骗干扰的判别
首先,SAR几何模型图如图3所示,X轴为方位向,Y轴为距离向,Z轴为高度向。平台高度为h,飞行速度为Vr。 在t=0时刻,SAR经过坐标(0,0,h)处。A点为一真实点目标,以A点为干扰源,欲在B点位置处产生一个假目标。
图3 方位向抗欺骗干扰几何模型Fig.3 Geometric model of azimuth anti-deception jamming
方位向欺骗干扰产生的机理与距离向有着明显不同,方位向假目标的产生是通过方位向多普勒频率调制所产生的。产生机理如下,假设SAR只存在单一收发共用天线,t=0时,天线也位于APC1点,坐标为(0, 0,h)。在A点处的干扰源接收到SAR回波信号后,采用方位向多普勒频率调制来产生方位向的假目标,其多普勒频率差值只与方位向上的间距有关。则多普勒频率差值表示为
其中,d表示A点和B点之间的距离。因此,为使得干扰信号与真实信号产生上述的多普勒频率差值,干扰机调制后的干扰回波信号为
针对上述欺骗干扰,基于MIMO-SAR的顺轨干涉抗方位向多普勒频率调制欺骗干扰的机理如下:取4发4收的MIMO-SAR系统中的一对天线为例,t=0时,天线1位于APC1(0, 0,h),天线2位于APC2(s,0,h),顺轨干涉的干涉基线长即为s。同样地,本文分别比较A点和B(B′)点干涉相位差。
首先计算A点干涉相位,A点相对于通道1的瞬时斜距为
A点相对于通道2的瞬时斜距为
于是,A点的干涉相位为
再计算真实目标B点的干涉相位,B点相对于通道1的瞬时斜距为
B点相对于通道2的瞬时斜距为
于是,B点的干涉相位为
又因为A, B两点的方位位置不同,导致多普勒中心到达A, B点的时刻不一样。假设飞机是正侧视的情况下,多普勒中心到达A点时刻为
对B点而言,多普勒中心到达B点时刻为
代入则有
将式(25)代入式(22),则式(22)退化为式(19)。所以,由此看出,如果是真实目标B点,则其顺轨干涉相位和参考点A点一致。而对于方位多普勒频率调制欺骗干扰目标B′点,B′点是由A点方位多普勒频率调制产生,式(25)不再成立,故而对B′点而言,干涉相位的表达式应该同式(22),此时不能退化为式(19)。因此,本文可以根据干涉相位与参考点A点是否一致来判别,一致的是B点真目标,不一致的是B′点假目标。
2.2.3 2 维混合基线判别欺骗干扰
实际战场环境中,仅仅实现距离向或者方位向的欺骗干扰是不够的,面目标才有更好的欺骗性和迷惑性,而实现面目标欺骗干扰的前提是实现2维点目标欺骗干扰。同时采用距离延时和方位向多普勒频率调制可以产生2维欺骗干扰。针对这种情况,则需要联合交轨干涉相位和顺轨干涉相位来做判别,因此至少需要2维混合基线,本文所提4发4收MIMO-SAR几何模型图如图4所示。
图4 2维抗欺骗干扰几何模型Fig.4 Geometric model for two-dimensional anti-deception jamming
图4 中,X轴为方位向,Y轴为距离向,Z轴为高度向。平台高度为h,飞行速度为Vr。 在t=0时刻,SAR经过坐标APC1(0, 0,h)处。APC2, APC3,APC4为另外3天线位置。取其中一对为例,APC1和APC2同时观测时,可以提取交轨干涉相位,APC1和APC3同时观测时,提取顺轨干涉相位。满足可以提取交轨干涉相位和顺轨干涉相位,那么判别上述欺骗干扰就有了依据,因此可以用2维混合基线的相位信息来判别目标。
2.3 基于MIMO-SAR的相位补偿抑制欺骗干扰原理
实现抗干扰意味着要对干扰目标有一定的抑制效果,而MIMO-SAR系统的相位补偿就是实现抑制的关键步骤。前面已经说到,对真目标和假目标而言,干涉相位之差一定不同,令任意位置处的点目标而言,若为真实目标则为B点,若为假目标则为B′点。B点相对通道1的相位为φ1,相对通道2的相位为 φ2, B点干涉相位为φ2-φ1=Δφ 。B′点相对通道1的相位为φ1′,相对通道2的相位为φ2′, B′点干涉相位为φ2′-φ′1=Δφ′。相位补偿就是统一将所有目标按真目标的相位做补偿,即将φ1通过补偿Δφ 的相位成为φ2。而对假目标而言,φ′1通过补偿Δφ 的相位不能变成φ2′,这是因为Δ φ 和Δ φ′不相等的缘故。这样一来,补偿后再多通道相干积累时,由于真目标相位差为0,可以实现相干积累,而假目标相位差不为0,只能实现非相干积累,相干积累的目标能量一定比非相干积累的目标能量强,从而实现了对假目标的抑制。
本文所提MIMO-SAR系统为4发4收系统,考虑等效相位中心的位置,以第1个通道发射第1个通道接收为基准,分别计算在方位频域上补偿相位因子。假设φ1为 交轨干涉的相位差,φ2为顺轨干涉的相位差,基线越长,需要补偿的相位差就越大。则通道1发1收需要补偿的相位为0,通道1发2收需补偿的相位为 φ1/2,通道1发3收需补偿的相位为φ2/2,通道1发4收需补偿的相位为( φ1+φ2)/2,依此类推,通道4发4收需补偿的相位为( φ1+φ2)。补偿后各个通道间累积,累积效果为
其中,s1为通道1发1收的信号能量,s∑为所有通道的信号能量累积和。
总结上述处理步骤,本文所提SAR抗欺骗干扰方法流程图如图5所示。
2.4 欺骗干扰抑制评估
图5 方法流程Fig.5 Method flow
为评估假目标的抑制效果,本文采用了雷达抗干扰改善因子(Electronic Counter Countermeasures Improvement Factor, EIF)作为评价因子,雷达抗干扰改善因子由Johnston提出,是目前为止唯一被IEEE采用的雷达抗干扰度量模型,其定义为:雷达未采用抗干扰措施时系统输出的干信比与采用抗干扰措施后系统输出的干信比的比值,即
其中,(J/S)为雷达未采用抗干扰措施时系统输出的干信比;(J/S)′为雷达采用抗干扰措施后系统输出的干信比。EIF值越大,表明雷达采取抗干扰措施后,要想有效地干扰雷达,必须付出更大的代价或更难以实现,因此,雷达的抗干扰性能就越好。
有了雷达抗干扰改善因子后,本文就能更好地说明4发4收的MIMO-SAR系统的好处。不失一般性地,本文用单发4收的SAR系统和4发4收的MIMO-SAR系统的EIF仿真,通过比较证明在相同天线数目和相同基线长度下,MIMO-SAR系统比多通道SAR系统抗干扰性能更好。
当系统为单发4收时,由上节可知,考虑等效相位中心的位置,以第1个通道发射第1个通道接收为基准,分别计算在方位频域上补偿相位因子,通道1发1收需要补偿的相位为0, 1发2收需补偿的相位为 φ1/2, 1发3收需补偿的相位为φ2/2, 1发4收需补偿的相位为 (φ1+φ2)/2,则单发4收的EIF值应为
其中,φ1为 交轨干涉的相位差,φ2为顺轨干涉的相位差,分子为通道的个数。同样地,当系统为4发4收时,参考上节的通道累积公式,则方型基线的EIF值应为
最后,为了进一步直观地说明4发4收的MIMO系统比单发4收的多通道系统所获得的雷达抗干扰改善因子(EIF)更好,设置了如下的仿真,X轴和Y轴分别和交轨干涉和顺轨干涉的相位差,取值介于0~2之间,Z轴表示EIF值,红色部分是4发4收的MIMO系统,绿色部分是单发4收的多通道系统。从图6中可以看出,红色部分一直处于绿色部分的上方,也就是说在相同天线数目和相同基线长度下,MIMO-SAR系统比单发多收SAR系统抗干扰性能更好。计算取平均值,MIMO-SAR系统比单发多收SAR系统雷达抗干扰改善因子提升3倍。
图6 4发4收与单发4收对比Fig.6 Comparison of 4-in-4-out and 4-in-1-out
3 仿真分析
仿真系统为机载C波段SAR,仿真参数如表1所示。
表1 仿真参数Tab.1 Simulation parameters
图7 点目标成像Fig.7 Point target imaging
图8 真实目标干涉相位Fig.8 Interferometric phase of real target
图9 干扰目标干涉相位Fig.9 Interferometric phase of jamming target
当欺骗干扰只采用了距离延时欺骗干扰时,利用2维混合基线的交轨干涉相位判别欺骗干扰。图7-图9中左边的是A点,也就是干扰源所处的位置,图7-图9右边的是B点,是真实目标,对比的是处于相同位置的B′点,是距离延时欺骗干扰目标,图8,图9方框中的数值表示交轨干涉相位。从图8和图9中对比可以看出,B点是真实目标时,干涉相位和A点处的干涉相位并不相同,而且相差较大。而B′点是欺骗干扰目标时,干涉相位和A点处的干涉相位基本相同,而且相差不大,与理论分析相符。
同样地,当欺骗干扰只采用了方位多普勒频率调制欺骗干扰时,利用2维混合基线的顺轨干涉相位判别欺骗干扰。图10-图12中上边的是A点,也就是干扰源所处的位置,图10-图12下边的是B点,是真实目标,对比的处于相同位置的B′点是多普勒频率调制欺骗干扰目标,图11,图12方框中的数值表示顺轨干涉相位。从图11和图12中对比可以看出,B点是真实目标时,干涉相位和A点处的干涉相位基本相同,而且相差不大。而B′点是欺骗干扰目标时,干涉相位和A点处的干涉相位并不相同,而且相差较大,与理论分析相符。
当欺骗干扰同时采用了距离延时和方位向多普勒频率调制欺骗干扰时,利用2维混合基线的交轨干涉和顺轨干涉相位联合判别。本文设计了如下的仿真方案:图13中,中心点A点有一真实点目标,同时也是干扰源所处的位置,其余周围8个点都是由A点处采用距离延时、方位多普勒频率调制或两者结合的方式产生的假目标,左右两端是两个真实点目标。提取交轨干涉和顺轨干涉相位分别如表2和表3所示。表中各处的干涉相位与上述对比,均符合距离向和方位向的假目标判决依据。而图14是抗干扰后的目标成像图,抗干扰采用的是4发4收的MIMO-SAR系统,处理流程为本文所介绍的算法处理流程。
对比图13和14可以发现:第1,图14中假目标得到了很大程度的抑制,雷达抗干扰改善因子(EIF)约为29.4 dB。第2,左右两端和中心点的真目标幅值变化0.17 dB,说明该抗干扰算法只针对了假目标进行抑制,而不会影响到真实点目标成像。
图10 点目标成像Fig.10 Point target imaging
图12 干扰目标干涉相位Fig.12 Interferometric phase of jamming target
4 结束语
图13 干扰目标成像Fig.13 Jamming target imaging
表2 9点目标交轨干涉相位Tab.2 Cross-track interference phase of nine-point target
表3 9点目标顺轨干涉相位Tab.3 Along-track interference phase of nine-point target
本文针对的是距离延时和方位多普勒频率调制欺骗干扰所带来的SAR抗欺骗干扰问题,解决问题的突破点是利用多维对抗1维具有优势,实现途径是通过相位补偿和相干积累来实现假目标的抑制,并通过对比说明了4发4收的MIMO-SAR系统所带来的好处,但不局限于4发4收。在相干积累过程中,如果能使用超分辨算法将真假目标回波的主瓣变窄,这将有利于假目标的抑制效果,这也是本文的进一步研究方向。
图14 抗干扰后目标成像Fig.14 Target imaging after anti-jamming