[发明专利]抗干扰方法和装置

说明书

技术领域

本发明涉及抗干扰滤波领域，尤其涉及一种应用于各种来向干扰环境下的卫星导航接收机以及其他抗干扰扩频通信系统的抗干扰方法和装置。

背景技术

全球卫星导航系统目前有美国全球定位系统(GPS)、俄罗斯全球卫星导航系统(GLONASS)、欧洲伽利略卫星定位系统(Galileo)、以及发展中的中国北斗定位系统，它们能提供全天候、实时、连续的高精度位置信息，已经广泛应用于各类军用和民用目标的定位、导航、授时和精密测量。但是卫星信号在传输和接收过程中会受到有意和无意的干扰。自适应阵列天线通过自适应调节各阵元信号的加权幅度和相位来改变阵列的方向图形状，使波束指向期望方向，而零点指向干扰方向，从而提高天线的增益和信干噪比，也称为空域滤波技术，是目前最有效的干扰抑制技术之一。

然而，单纯的空域滤波受到阵列自由度个数的限制，不能满足复杂的干扰以及多径环境下的应用，而且空域滤波对空间不同于信号来向的干扰可以有效地抑制，但对于与信号同向的窄带干扰抑制程度不够，同时滤除部分有用信号。空域频域联合滤波对点频干扰处理时有效位、系数精度要求较高，同样有效位精度下抑制窄带干扰效果比抑制宽带干扰效果差，通过频域预处理窄带干扰的方法解决了该问题。

发明内容

本发明提供了一种抗干扰方法和装置，解决了单一滤波方式局限性的问题。

一种抗干扰方法，包括：

对中频信号进行频域窄带干扰预处理，得到频域数组；

对所述频域数组按照多个频点分别进行空域滤波；

对空域滤波后的所述频域数组进行多点IFFT变换，得到抗干扰处理后的时域中频数据。

优选的，对中频信号进行频域窄带干扰预处理，得到频域数组包括：

积累一定长度的中频AD采样信号，所述中频AD采样信号中包含卫星信号、窄带、宽带干扰以及背景噪声，将所述中频AD采样信号进行FFT变换，得到幅度谱ξ_n；

根据以下表达式统计带内幅度谱模值的平均值：

μ=1KΣn=1K|ξn|;]]>

根据二阶矩算法计算干扰门限T＝μ+βσ，其中β为门限优化因子，初始值设为1，σ为幅度标准方差，根据表达式确定；

对所述中频AD采样信号进行FFT变换后得到的频域数组进行窄带干扰预处理，超过干扰门限T的频域数据置为平均值μ，否则不予处理，并记录超过门限值的次数Q；

当Q>10％带内频点数时，β+0.1；当Q<2％带内频点数时，β-0.1。

优选的，对所述频域数组按照多个频点分别进行空域滤波包括：

每个阵元接收中频AD采样信号经过相同的频域窄带干扰预处理后，经过多点FFT变换，得到频域数组X：

X＝[X₁₁,X₁₂,…X_1K；X₂₁,X₂₂,…X_2K；…X_M1,X_M2,…X_MK]；

根据所述频域数组，计算协方差矩阵R_k＝E[X_kX_k^H](M×M维)，其中，X_k＝[X_1k,X_2k,…X_Mk]^T k＝1,2,...,K；

求解线性约束的最小方差方程，由线性约束最小方差准则，确定空频滤波器的最优化问题如下：