杜 盈, 张劲东, 尹明月, 蒋宜林
(南京航空航天大学电子信息工程学院, 江苏南京 211106)
随着电子对抗技术的不断成熟和发展,雷达波形设计成为了雷达研究领域的热点。学者们通过设计具有低距离旁瓣的雷达发射波形来提升雷达的检测性能;而干扰信号经过滤波器后的输出水平可用滤波器与干扰信号协方差矩阵的二次型表示,通过最小化干扰输出功率可提升雷达的抗干扰性能。
很多学者基于模糊函数来设计发射波形。文献[1]采用半正定规划(Semi-Definite Programming, SDP)的方法最小化旁瓣电平,设计在强杂波环境下可被探测的宽带恒模信号,文献[2]提出一种FADPM方法(Fractional-Alternating Direction Penalty Method),通过对模糊函数的旁瓣控制来保证雷达的探测性能,文献[3-5]分别采用极大极小算法(Majorization-Minimization, MM)、梯度下降算法(Gradient Descent, GD)和循环算法最小化积分旁瓣电平(Integrated Sidelobe Level, ISL),但这些研究都只涉及到距离旁瓣的抑制。考虑到雷达对抗干扰的能力,文献[6-9]对稀疏信号的功率谱密度和自相关函数性能同时进行优化,设计稀疏频谱波形,可使雷达避免被其他辐射源信号干扰,且具有较强的探测弱小目标的能力,但这些文献的研究都是针对同频干扰的情况。
本文设计可对抗压制干扰的相位编码信号作为发射信号,将干扰信号经匹配滤波器的输出功率写进目标函数,联合最小化发射信号距离旁瓣和干扰信号输出功率,采用ADMM算法求解,使得雷达兼具良好的目标探测性能和抗干扰性能。同时,考虑到ADMM算法运算量大、计算复杂的问题,本文结合PMLI算法提出一种新的CA法求解问题模型。
压制式干扰就是采用干扰发射机发射大功率的噪声干扰信号或密集假目标信号,使干扰信号进入雷达接收机,使有用的回波信号被噪声或假目标覆盖,造成雷达接收机信噪比(Signal-to-Noise Ratio, SNR)大幅下降,从而使雷达检测或跟踪不到目标。压制干扰可分为射频噪声干扰、噪声调幅干扰、噪声调频干扰、噪声调相干扰等形式,本文仅以噪声调频信号为例。
一般噪声调频干扰的时域表达式可表示如下:
(1)
式中,为噪声调频信号的幅度,为噪声调频信号的中心频率,为调制过程中的调频斜率,其被积函数为调制噪声(′),一般地,它是均值为零、且广义平稳的随机过程,为随机变量,在[0,2π]内均匀分布,且与()相互独立。
一个点的相位编码信号序列可表示为
=[,,…,-1]=[ej,ej,…,ej-1],∈
(2)
发射信号旁瓣协方差矩阵为
(3)
所以发射信号距离旁瓣可表示为
=
(4)
设干扰信号协方差矩阵为,∈×,干扰信号经匹配滤波器的输出功率为
=
(5)
折中考虑干扰信号功率和发射信号旁瓣,引入权系数,那么离散相位形式约束下的联合最小化干扰信号输出功率和发射信号距离旁瓣问题模型可表示为
(6)
式中为离散相位个数。
问题的目标函数可改写为如下形式:
(7)
在该问题中引入辅助变量和约束=,问题即可表示为
(8)
根据问题写出增广拉格朗日方程:
(9)
式中和为拉格朗日乘子,为二次惩罚项系数。记=(+j),则式(9)可重写为
(10)
记(),(),()为第次迭代后的值,各值更新如下:
1) 更新,此时将(),()视为已知量:
(11)
这是一个无约束优化问题,使用基于内部映射牛顿法的子空间置信域法求解。令
(12)
(13)
(14)
利用式(12)~ (14)将目标函数式(11)由复数形式转化为实数形式:
(15)
对目标函数求梯度:
(16)
2) 更新,此时将(+1),()视为已知量:
(17)
的更新公式为
(18)
3) 更新:
(+1)=()+(+1)-(+1)
(19)
ADMM求解优化问题模型的步骤如表1所示。
2.1节中的ADMM算法在更新时需要进行
表1 ADMM求解优化问题P2步骤
数据形式的转化和梯度的计算,运算量大,计算时间长。针对这一问题,本小节将在ADMM方法的基础上作改进。
对于发射信号的更新,
‖-()+()‖
(20)
令=(()-()),忽略常数部分,式(20)可表示为
(21)
那么求解发射信号可表示为
(22)
(23)
这样可通过在每次迭代时求解如下最近向量问题来更新发射信号:
(24)
显然其最优解为
(+1)=()
(25)
在原ADMM算法的基础上引入PMLI方法,新的复合算法克服了ADMM算法在更新时计算复杂的问题,减少算法计算时间。
复合算法求解DPCQP问题的步骤如表2所示。
表2 CA求解优化问题P2步骤
从图1的脉压结果可以看到,初始的发射信号被压制干扰淹没,目标无法被检测,而优化后干扰被抑制,信号能够被检测到,两种优化算法得到的发射信号波形均可提升雷达对抗压制干扰的能力。
(a) 优化前
(b) ADMM优化后
(c) CA优化后图1 优化前后信号脉压结果
雷达的探测性能由距离旁瓣水平衡量,优化前后发射信号的自相关处理结果如图2所示。从图2可以看出,两种算法优化后的信号旁瓣均降低,雷达探测弱小目标的能力得以提升。
(a) 优化前
(b) ADMM优化后
(c) CA优化后图2 优化前后发射波形自相关处理结果
下面对比两种优化算法的性能。在Intel i5-9400HQ CPU、16G内存、Matlab2016a仿真平台下,进行100次蒙特卡洛实验取平均值,两种算法的运算时间、干扰输出功率减少量Δ、信号旁瓣减少量Δ如表3所示,ADMM和CA优化的目标函数迭代曲线如图3所示。
表3 两种算法性能对比
图3 两种算法迭代曲线
两种算法均能收敛,CA抑制旁瓣的能力不及ADMM算法,但对抗干扰的能力强于ADMM,这是由于CA中发射信号的更新取决于矩阵,在强干扰环境下,干扰信号对的影响要远大于发射信号,而ADMM算法中的更新是采用拟牛顿法,每一次迭代时只要求知道目标函数的梯度,受到梯度变化量的影响而非梯度本身数值大小。本文提出的ADMM和CA两种算法均可提升雷达的抗干扰性能和探测性能,但CA相对ADMM算法计算速度有了30多倍的提升。
本文以联合最小化干扰信号输出功率和发射信号距离旁瓣为准则建立目标函数,采用ADMM算法求解,设计具有较低距离旁瓣并能对抗压制干扰的相位编码波形。针对ADMM算法中更新发射信号时计算量大、运算复杂的问题,提出一种新的CA算法。仿真结果表明,本文提出的两种算法均能提升雷达的探测性能和抗干扰性能,CA算法在旁瓣抑制方面弱于ADMM算法,但抗干扰能力强于ADMM算法,且计算速度相对ADMM算法有很大提升。