一种新型的空时联合抗干扰方法
【技术领域】
[0001] 本发明涉及一种卫星导航设备上的空时阵列天线抗干扰方法，尤其是涉及基于先 利用频域滤波方法滤除导航信号方向上的窄带干扰，再利用变步长算法求取空时二维权矢 量的卫星导航信号抗干扰方法，具体为一种新型的空时联合抗干扰方法。
【背景技术】
[0002] 随着卫星导航应用电磁空间环境的日益复杂，导航终端精密定位也变的日益困 难，各种导航终端抑制干扰技术随之产生。传统的接收机抗干扰技术有空域滤波技术、时域 滤波技术和频域滤波技术，这些抗干扰技术在性能上均存在一定的缺陷。空时自适应处理 技术在不增加阵元的情况下，可以抑制的干扰数目大大增加，且具有在宽带内处理干扰的 能力，从而使接收机的抗干扰能力有质的提高。
[0003] 目前，空时抗干扰技术的关键点在于空时权矢量求解的实现方法。权矢量求解的 一般流程图如图（1)所示。在阵列天线的各个阵元上增加相同数目的延迟抽头，从而形成 空时二维联合处理结构，从单个阵元的角度来看，形成时域FIR结构；从同一个时间节点来 看，形成空域滤波结构。通过空时联合处理可以抑制单频、窄带、宽带干扰信号。
[0004] 空时二维自适应抗干扰处理技术通常使用线性约束最小方差（LCMV)准 贝1J，原理是根据卫星信号功率远远低于噪声和干扰功率，从而对空时联合处理结 构中的空域和时域进行约束以调整权值，使输出信号方差最小，大大削弱干扰能 量。最小方差无失真响应（MVDR)准则是LCMV准则的一种特例。当卫星信号空时 导向矢量S己知时，约束卫星信号空时方向上的响应g= 1，输出信号功率最小即： Wopt= arg minVVHR-^ s.t WhS = I (1)，基于线性约束最小方差准则的最优 权值闭合解可以表达为：W = RfS (2)，式中，W= [ωη，ω12，…，ωΜΡ]τ为MPXl维向 量，M为阵列天线的阵元数，P为延迟单元数目；Rx= E(XXh)为MPXMP维的接收信号协方 差矩阵；X= (X11, X12，…，xmp)tSMPX1信号向量；S为MPXl维的约束向量。直接利用（2) 式进行权值计算时，会引入矩阵求逆公式，当M，P维数较大时，可获得较好的抗干扰效果， 并且增加抗干扰数目。但其运算量增大，大的运算量会导致系统复杂性增加，导致系统的实 时性变差。鉴于此，一般采用递推算法来求解权值。
[0005] 设阵元1接收到的信号X1 (η)为期望信号，即d(n) =X1 (η)。则滤波器的输出 为：y (n) = W 'HX '，e (n) = Wjd (n) -y (η) (3)，其中 W ' = (ω 21，ω 22，…，ω ΜΡ)τ，X ' = (χ21，χ22, 一,Xm^W1= [1，1，…，1] τ，利用最陡下降法可以推导出W的递推公式为：Wn+1'= Wn' +2μΧηθ(ηΓ (4)，式中μ为一常数，称为收敛因子。根据式（4)不断迭代更新就可以 得到最优权值Wtjpt的估计值。稳定后的权值输出为：W = [W Λ W。/]1 (5)，由上述可知，步 长固定的迭代算法能够避免矩阵的求逆过程，从而使计算量得到明显的降低。但是固定步 长的迭代算法在求解阵列加权值会遇到收敛速度与稳态误差不兼容的问题，这会导致算法 在抗干扰能力方面有所下降。
[0006] 当干扰信号的方向与导航信号相同的时候，会在其来波方向产生很深的零陷，导 致导航信号受到严重减弱。目前常用结合IIR陷波滤波器和空时FIR滤波器的STAP算法， 可以解决这方面的问题。但是由于FIR滤波器抽头的限制，自适应处理器的频率分辨力很 低，同时由于窄带干扰的频谱和统计特性是未知的，因此该方法不适合捕获快速变化的信 号。
【发明内容】
[0007] 本发明为了有效地抑制与导航信号同方向的窄带干扰，避免空时抗干扰算法计算 过程复杂、难以工程实现的缺点，提供了一种新型的空时联合抗干扰方法。
[0008] 本发明是采用如下的技术方案实现的：一种新型的空时联合抗干扰方法，包括以 下步骤：
[0009] 每个阵元首先对阵列接收信号中的每一路进行采样，然后进行射频下变频转换并 数子化，得到数子化中频?曰号 Xi，Xi= (X ii，Xi2，· · ·，xip)，i = 1，2,3，· · · M ;
[0010] 将数字化中频信号Xi输入到LMS频域滤波器中，通过滤波器中的频域LMS 迭代算法将导航信号方向的窄带干扰滤除，得到过滤信号2"1〇 = (zn，zi2，…，ziP) (i = 1，2,…，M)，频域LMS迭代算法的频域权值系数迭代公式为：U(k+1)= a U (k) +2uE (k) *F (k)，其中，U (k+1)，U (k)分别为k+1和k时刻的滤波系数向量；α为权值 因子；μ为收敛因子；F(k)为数字化中频信号傅里叶变换；E(k) = F(k)-U(k)*F(k)， 为误差信号；*为向量点乘，则ZiGO = Γ1 (U(K) *F (K))，其中Γ1为傅里叶逆变换；
[0011] 将过滤信号ZiGO送入空时阵列信号处理滤波器中，进一步消除窄带宽带干扰，该 滤波器算法的步长更新公式为：
【主权项】
1. 一种新型的空时联合抗干扰方法，其特征在于包括以下步骤： 每个阵元首先对阵列接收信号中的每一路进行采样，然后进行射频下变频转换并数字 化，得到数字化中频信号Xi， Xi(xn，Xu，? ? ?，Xip)，i1，2，3，...M; 将数字化中频信号\输入到LMS频域滤波器中，通过滤波器中的频域LMS迭代算法将 导航信号方向的窄带干扰滤除，得到过滤信号Zi(k) = (zn,zi2，…，ziP)(i= 1，2, "^M)， 频域LMS迭代算法的频域权值系数迭代公式为：U(k+1) =aU(k)+2yE(k)*F(k)，其中， U(k+1),U(k)分别为k+1和k时刻的滤波系数向量；a为权值泄漏因子；y为收敛因子； F(k)为数字化中频信号Xi的傅里叶变换为；E(k) =F(k)-U(k)F(k)为误差信号；*为矩阵 点乘，则Zi(k) =r1(U(K)*F(K))，其中r1为傅里叶逆变换； 将过滤信号Zjk)送入空时阵列信号处理滤波器中，进一步消除窄带宽带干扰，该滤波 器算法的步长更新公式为
，变步长 最小均方误差算法迭代出最优权值的步骤为：阵列信号处理滤波器的第2个到第M个阵元 接收的信号可以表示为：Z(k) = (Z2，Z3，…ZM)H，经空时阵列信号处理后可以表示为：Y(k) =W'HZ(k)，将第一个阵元接收到的信号与接收信号Y(k)作差得到阵元的误差信号：e(k) zW^GO-YGO;利用最陡下降法可以推导出W的递推公式为：Wk+1' =Wk' +2y(k)Z(k) e(k)%则稳定后的最优加权值为W。#，此时空时阵列信号处理滤波器输出纯净的导航信号
【专利摘要】本发明涉及一种卫星导航设备上的空时阵列天线抗干扰方法，具体为一种新型的空时联合抗干扰方法。本发明提出了一种新型的空时联合抗干扰方法，该方法利用改进的频域窄带干扰滤出技术来提前抑制掉导航信号方向上的窄带干扰，再结合变步长最小均方误差实现空时阵列权值的迭代求解，以降低该方法的运算复杂度。本发明由于能够抑制导航信号方向的干扰，算法复杂度低，实现简单，特别适合未来小型化高动态抗干扰终端的应用。
【IPC分类】G01S19-21
【公开号】CN104808219
【申请号】CN201510201630
【发明人】李灯熬, 刘金强, 赵菊敏, 牛文慧, 李帅
【申请人】太原理工大学
【公开日】2015年7月29日
【申请日】2015年4月24日