一种运动干扰稳健滤波方法与流程

本发明属于雷达技术领域，更进一步涉及运动干扰条件下自适应波束形成信号处理技术，可用于先验信息不足条件下的运动干扰稳健滤波。

背景技术：

传统的自适应信号处理技术利用协方差求逆得到自适应权值，在干扰处形成凹口。这种方法适用于天线阵列与干扰相对静止的条件下，然而实际环境中往往存在干扰源快速变化的情况，这使得协方差矩阵估计快拍数难以满足，干扰不能被方向图凹口完全抑制，导致系统滤波性能降低。因此阵列运动干扰稳健滤波问题亟需解决。

目前，针对运动干扰抑制问题，已有多种自适应抗干扰方法被提出。实时rls估计方法通过不断更新自适应权值达到抗运动干扰的目的，但是考虑到干扰运动是随机的，且很难进行实时预测，因此这种方法形成的阵列凹口无法精确对准干扰源。导数约束法将方向图在凹口处高阶导数置零，从而使干扰处零陷展宽。mailloux-zatman算法(m-z算法)将接收协方差矩阵与零陷展宽矩阵相乘，从而得到凹口展宽的方向图，相比导数约束法，计算量有所减小。为了进一步提高自适应波束形成算法的稳健性，有些学者将m-z方法与对角加载方法结合，减小了噪声和误差对算法性能的影响。制约上述运动干扰滤波性能最关键的因素是凹口展宽大小的确定，过大的凹口会引起周围目标信号能量损失，太小的凹口无法完成运动干扰信号抑制。

目前国内外对于抗运动干扰中凹口展宽量的研究未见公开，因此亟需一种能不依赖凹口展宽大小的稳健运动干扰抑制方法，这也是阵列抗干扰工程化需考虑的重要问题之一。

技术实现要素：

本发明解决的技术问题是：克服现有技术的不足，提供公开一种运动干扰稳健滤波方法，解决在实际工程中先验信息不足条件下的自适应阵列抗运动干扰问题

本发明的技术方案是：一种运动干扰稳健滤波方法，步骤如下：

步骤1：获取接收数据的协方差矩阵r；

步骤2：根据协方差矩阵获得新的零陷展宽协方差矩阵r1：

步骤3：对步骤2中零陷展宽协方差矩阵引入对角加载技术，得到新的协方差矩阵r2；

步骤4：根据步骤3得到的新的协方差矩阵，求得自适应权值进而得到能够在干扰位置处形成零陷的自适应方向图。其中θ0表示感兴趣的目标方向，d代表天线间距，

所述步骤1中获取接收数据的协方差矩阵r的具体表达形式为：

其中k是快拍数；x(i)是目标所在距离的第i个脉冲的多通道数据，x(i)^h表示对x(i)进行共轭转置处理。

所述步骤2中新的零陷展宽协方差矩阵r1的具体表达形式为：

r1＝cr+cbrb^h；

矩阵c中元素定义如下：

ck,l＝sinc((k-l)δ/π)；

其中1≤k,l≤n，n为天线阵元数，δ为展宽系数；b＝diag(0,1,…n-1)，对b其归一化后得

步骤3所述的对角加载技术，对角加载量大小选取方法如下：

1)对r进行特征值分解；

2)将特征值按绝对值大小排序，当接收数据中存在干扰信号时，特征值曲线必定存在拐点，求得拐点以及拐点两侧特征值v1、v2和v3，则对角加载量大小为σ＝v1+v2+v3/3。

所述步骤3中计算得到新的协方差矩阵r2的具体方法为：r2＝r1+σi；其中i是维度为n×n的单位阵，σ为对角加载量大小。

所述步骤4中自适应权值其中θ0表示感兴趣的目标方向，d代表天线间距，

本发明与现有技术相比有益效果为：

本发明方法属于自适应阵列抗运动干扰算法，相比经典的几种运动干扰抑制算法，所提算法将零陷展宽技术与导数约束方法相结合，用得到的新协方差矩阵求自适应权值。在此基础上，为了提高自适应波束形成对噪声和误差的稳健性，引入了对角加载技术，加载量通过寻找协方差矩阵的特征值曲线拐点得到，有效避免了加载量过大导致的干扰信号被淹没以及加载量过小导致噪声分量引起性能损失的问题。该发明对零陷展宽系数选取依赖性较低，且对角加载量不需要设置固定的经验值，有利于稳健运动干扰滤波，适用于实际工程中先验信息不足条件下的自适应阵列抗运动干扰情况。

附图说明

图1为本发明方法处理流程。

图2(a)是在展宽系数δ＝0.05的情况下采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法形成的自适应阵列方向图，图2(b)是在展宽系数δ＝0.00005的情况下采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法形成的自适应阵列方向图。

图3(a)是在展宽系数δ＝0.05的情况下采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的输出信干噪比曲线，图3(b)是在展宽系数δ＝0.00005的情况下采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的输出信干噪比曲线。

图4(a)是在展宽系数δ＝0.05的情况下得到的协方差矩阵特征值分布。图4(b)是在展宽系数δ＝0.00005的情况下得到的协方差矩阵特征值分布。

图5(a)给出了展宽系数δ＝0.05、目标与干扰方位角较近时(干扰方位角18°)采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的自适应方向图，图5(b)给出了展宽系数δ＝0.00005、目标与干扰方位角较近时(干扰方位角18°)采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的自适应方向图。

图6(a)给出了展宽系数δ＝0.05、目标与干扰方位角较近时(干扰方位角18°)采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的输出信干噪比曲线，图6(b)给出了展宽系数δ＝0.00005、目标与干扰方位角较近时(干扰方位角18°)采用协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法得到的输出信干噪比曲线。

具体实施方式

下面对本发明实施及效果作进一步的详细描述。

本发明的使用场景为：本发明可应用于先验信息不足条件下的自适应阵列抗运动干扰情况中。将零陷展宽技术与导数约束方法相结合，用得到的新协方差矩阵求自适应权值。同时引入了自适应对角加载技术提高自适应波束形成对噪声和误差的稳健性。其实现步骤如下：

步骤1获取接收协方差矩阵：

其中k是快拍数；x(i)是目标所在距离的第i个脉冲的多通道数据，x(i)^h表示对x(i)进行共轭转置处理。

步骤2根据协方差矩阵获得新的零陷展宽协方差矩阵r1：

矩阵c中元素定义如下：

ck,l＝sinc((k-l)δ/π)

其中1≤k,l≤n，n为天线阵元数，δ为展宽系数。b＝diag(0,1,…n-1)，对其归一化后

步骤3对步骤2中零陷展宽协方差矩阵引入对角加载技术，对角加载量大小选取方法如下：

1)对r进行特征值分解。

2)将特征值按绝对值大小排序，当接收数据中存在干扰信号时，特征值曲线必定存在拐点，求得拐点以及拐点两侧特征值v1、v2和v3，对角加载量大小定义为σ＝v1+v2+v3/3。

3)新的协方差矩阵定义为：r2＝r1+σi。

步骤4根据最终的协方差矩阵求得自适应权值得到归一化方向图。其中θ0表示感兴趣的目标方向，d代表天线间距。

下面通过仿真实验对本发明的效果做进一步说明。

图2对比了协方差矩阵求逆算法、本发明算法、m-z算法以及对角加载算法形成的阵列方向图，可以看出相比其他三种算法，专利所提算法对零陷展宽系数的选取依赖性较小，可以在不同展宽系数条件下获得较好的运动干扰抑制效果。另一方面，专利中对角量选取方法可以保证阵列抗运动干扰的同时对噪声等具有一定的稳健性。图3中不同方法的输出信干噪比进一步验证了上述结论。图4为不同展宽系数对应的协方差特征值分布，当存在干扰时，从特征值曲线中很容易找到拐点以及拐点左右的特征值得到对角加载量。其中图2、图3和图4仿真参数设置如下：目标信噪比为5db，干扰对应干噪比为40db，噪声功率0db，目标和干扰对应方位角分别为10°和35°，阵元数和快拍数设置为16和300。

为了进一步验证本发明抗运动干扰的稳健性，图5和图6分别给出了目标干扰方位角较近时(干扰方位角18°)不同方法对应的方向图和输出信干比曲线，当展宽系数较小时，m-z方法方向图凹口展宽有限，因此滤波性能不如专利方法；当展宽系数较大时，m-z方法方向图凹口展宽程度较大，但由于凹口与目标信号较近，导致方向图出现畸变，所提方法通过合理选取对角加载量，保证了抑制运动干扰的同时不影响目标方向图的指向。

仿真结论：仿真结果表明，相比经典的几种自适应阵列抗干扰方法，本发明方法提升了在没有运动干扰先验信息条件下的滤波性能，在不同零陷展宽系数以及临近干扰等情况下，均有较稳健的抗运动干扰能力。