[发明专利]一种新型的空时联合抗干扰方法

说明书

技术领域

本发明涉及一种卫星导航设备上的空时阵列天线抗干扰方法，尤其是涉及基于先利用频域滤波方法滤除导航信号方向上的窄带干扰，再利用变步长算法求取空时二维权矢量的卫星导航信号抗干扰方法，具体为一种新型的空时联合抗干扰方法。

背景技术

随着卫星导航应用电磁空间环境的日益复杂，导航终端精密定位也变的日益困难，各种导航终端抑制干扰技术随之产生。传统的接收机抗干扰技术有空域滤波技术、时域滤波技术和频域滤波技术，这些抗干扰技术在性能上均存在一定的缺陷。空时自适应处理技术在不增加阵元的情况下，可以抑制的干扰数目大大增加，且具有在宽带内处理干扰的能力，从而使接收机的抗干扰能力有质的提高。

目前，空时抗干扰技术的关键点在于空时权矢量求解的实现方法。权矢量求解的一般流程图如图(1)所示。在阵列天线的各个阵元上增加相同数目的延迟抽头，从而形成空时二维联合处理结构，从单个阵元的角度来看，形成时域FIR结构；从同一个时间节点来看，形成空域滤波结构。通过空时联合处理可以抑制单频、窄带、宽带干扰信号。

空时二维自适应抗干扰处理技术通常使用线性约束最小方差(LCMV)准则，原理是根据卫星信号功率远远低于噪声和干扰功率，从而对空时联合处理结构中的空域和时域进行约束以调整权值，使输出信号方差最小，大大削弱干扰能量。最小方差无失真响应(MVDR)准则是LCMV准则的一种特例。当卫星信号空时导向矢量S己知时，约束卫星信号空时方向上的响应g＝1，输出信号功率最小即：基于线性约束最小方差准则的最优权值闭合解可以表达为:式中，W＝[ω₁₁,ω₁₂,…,ω_MP]^T为MP×1维向量，M为阵列天线的阵元数，P为延迟单元数目；R_X＝E(XX^H)为MP×MP维的接收信号协方差矩阵；X＝(x₁₁,x₁₂,…,x_MP)^T为MP×l信号向量；S为MP×l维的约束向量。直接利用(2)式进行权值计算时，会引入矩阵求逆公式，当M，P维数较大时，可获得较好的抗干扰效果，并且增加抗干扰数目。但其运算量增大，大的运算量会导致系统复杂性增加，导致系统的实时性变差。鉴于此，一般采用递推算法来求解权值。

设阵元1接收到的信号x₁(n)为期望信号，即d(n)＝x₁(n)。则滤波器的输出为：y(n)＝W‘^HX′，e(n)＝W₁^Hd(n)-y(n)(3)，其中W‘＝(ω₂₁,ω₂₂,…,ω_MP)^T,X′＝(x₂₁,x₂₂,…,x_MP)^T,W₁＝[1,1,…,1]^T，利用最陡下降法可以推导出W的递推公式为:W_n+1′＝W_n′+2μX_ne(n)^*(4)，

式中μ为一常数，称为收敛因子。根据式(4)不断迭代更新就可以得到最优权值W_opt的估计值。稳定后的权值输出为：W＝[W₁^T,W_opt^T]^T(5)，

由上述可知，步长固定的迭代算法能够避免矩阵的求逆过程，从而使计算量得到明显的降低。但是固定步长的迭代算法在求解阵列加权值会遇到收敛速度与稳态误差不兼容的问题，这会导致算法在抗干扰能力方面有所下降。

当干扰信号的方向与导航信号相同的时候，会在其来波方向产生很深的零陷,导致导航信号受到严重减弱。目前常用结合IIR陷波滤波器和空时FIR滤波器的STAP算法,可以解决这方面的问题。但是由于FIR滤波器抽头的限制,自适应处理器的频率分辨力很低,同时由于窄带干扰的频谱和统计特性是未知的,因此该方法不适合捕获快速变化的信号。

发明内容

本发明为了有效地抑制与导航信号同方向的窄带干扰，避免空时抗干扰算法计算过程复杂、难以工程实现的缺点，提供了一种新型的空时联合抗干扰方法。

本发明是采用如下的技术方案实现的：一种新型的空时联合抗干扰方法，包括以下步骤：