基于STFRFT的脉冲干扰抑制方法研究

来源：核心期刊咨询网时间：12

摘要：摘 要：为了解决航空无线电导航服务频段附近的脉冲干扰影响接收机工作的问题，并进一步提高GNSS系统的抗干扰能力，提出了基于STFRFT(short time fractional Fourier transform)的脉冲干扰抑制方法。首先确定接收信号的最佳旋转阶次，再对信号进行STFRFT，通

　　摘 要：为了解决航空无线电导航服务频段附近的脉冲干扰影响接收机工作的问题，并进一步提高GNSS系统的抗干扰能力，提出了基于STFRFT(short time fractional Fourier transform)的脉冲干扰抑制方法。首先确定接收信号的最佳旋转阶次，再对信号进行STFRFT，通过坐标旋转得到时频面的二维分布，在最佳阶次下采用自适应时变滤波器实现对干扰和信号的分离，进而抑制干扰。仿真结果表明，和几种传统方法相比，基于STFRFT的自适应滤波算法剔除干扰效果最好，在抑制脉冲干扰的同时保留了更多有用信号。结合STFRFT技术和自适应时变滤波技术提高了抑制脉冲干扰的能力，可为航空无线电导航服务的脉冲抗干扰技术提供参考依据。

　　关键词：信号检测;STFRFT;脉冲干扰抑制;最佳阶次;坐标旋转;时变滤波器

　　作者:王晓君; 薛琳博; 王彦朋

　　在GNSS的建设和进行现代化改造的过程中会面临与航空无线电设备的电磁兼容的问题。伽利略E5、北斗B2和GPS L5等这些不同频点的卫星信号都是GNSS的重要组成部分，它们的中心频率和现有的航空无线电导航服务(aeronautical radio navigation services，ARNS)频段会发生重叠。ARNS频段附近工作的设备主要有测距仪(distance measuring equipment，DME)系统和塔康(tactical air navigation system，TACAN)系统，这些设备工作时都不可避免地产生脉冲干扰信号。研究调查表明，测距仪系统发射的大功率脉冲信号是影响GPS L5等卫星信号的主要因素[1]，所以必须采取有效的干扰抑制手段，避免脉冲干扰影响卫星导航系统的正常工作。

　　目前抑制脉冲干扰常用的有时域脉冲消隐法、频域陷波法和时域频域结合等方法。时域消隐方法首先检测输入信号的幅度，然后把超过所设门限的信号置零，以此来抑制干扰信号，其缺点是会滤除在脉冲干扰较为集中区域的有用卫星信号，损失较为严重。文献[2]提出了频域陷波法，在频域里对脉冲干扰进行抑制，然后通过陷波滤波法剔除干扰，干扰抑制效果好于时域消隐法，但会滤除掉脉冲同频率处的有用信号。时频域结合法是上述2种方法的综合，干扰抑制效果优于上述2种方法。近年来，学者针对脉冲干扰又提出了很多有效的干扰抑制方法。文献[3]提出了小波包抑制脉冲干扰的方法，将接收信号通过小波包变换转变到小波包系数域中，通过比较系数的绝对值和阈值的大小界定干扰，但如果阈值选取不当，该方法在一定程度上也会丢失有用卫星信号。文献[4]提出了自适应空域脉冲抗干扰算法以剔除脉冲干扰，该方法涉及到阵列信号结构且需要计算求得信号的协方差矩阵，抑制效果较好但运算复杂。

　　由于来自不同DME站台的脉冲干扰信号是非平稳的，与噪声干扰和卫星信号之间存在较强的时频耦合，使得常见的干扰抑制方法难以有效地实现信号和干扰的分离。随着科学技术的快速发展，分数阶傅里叶变换(fractional Fourier Transform， FRFT)的研究逐渐成为热点。其对信号有良好的时频聚焦性质，是分析非平稳信号的新工具，在水波信号、地震信号、通信信号等领域应用广泛，且在处理复杂的时频重叠信号领域拥有广阔的应用前景。

　　本文从脉冲干扰的非平稳特性入手，研究了分数阶傅里叶变换抑制脉冲干扰的方法，在此基础上对信号进行加窗处理，探究短时分数阶傅里叶变换下脉冲的干扰抑制方法和性能，并与其他干扰抑制方法进行了分析比对。

　　1 系统模型

　　本文以GPS接收机为例，把包含噪声和脉冲干扰的卫星信号采样到中频进行处理，此时接收信号的模型为

　　式中fJ和θJ为载波频率和载波相位。基带和调制后的DME脉冲信号时域图如图1所示。

　　2 分数阶傅里叶变换

　　分数阶傅里叶变换是现代信号处理的重要工具，是当下使用广泛的一种信号分析方法。分数阶傅里叶变换通过信号在时间轴上旋转合适的角度，不仅能反映时域信息，还能反映频域信息。利用分数阶傅里叶变换将信号和干扰分离，可以极大程度地滤除干扰，保留有用的GNSS信号。接收信號的分数阶傅里叶变换定义如下[6-7]：

　　3 基于短时分数阶傅里叶域的干扰抑制方法

　　3.1 短时分数阶变换的干扰抑制原理

　　利用STFRFT抑制脉冲干扰的基本思想[10-11]是：首先对接收的GNSS信号做FRFT，计算出(p，u)平面下的分布，在此平面搜索峰值，确定它的最优阶次p0，然后对接收信号做p0阶次的STFRFT，由于脉冲干扰和GNSS信号在分数阶域的能量聚集度不同，在u0域采用自适应时变滤波器剔除脉冲干扰。基于STFRFT抑制脉冲干扰的基本原理框图，如图2所示。

　　由于对干扰抑制都是基于最优阶次下进行的，为了得到更好的干扰抑制效果，确定最优阶次尤为重要。传统的二维搜索方法在搜索过程中必须选择较小的搜索步长才会满足精度要求，这便增加了计算的复杂度。本文采用牛顿迭代方法[12]，其是一种用目标函数在点x0处的二阶泰勒展开式代替目标函数的近似方法，迭代公式为

　　3.2 自适应时变滤波器

　　自适应时变滤波器是对传统带通滤波器的一个改进，能更好地适应脉冲实时变化的特性，在每个时刻时变滤波器的中心频率会随着接收信号的瞬时频率自适应发生改变，且对于处理多分量信号有着良好的效果。设接收信号为xi(t)，自适应时变滤波器的基本思想如下[13-15]。

　　在STFRFT的基础上对信号进行时变滤波的方法可以概括为3个步骤：分析、加重和综合。其中加重函数M(t，u)在通域内的值为1，在通域之外的值为0。Mi(t，u)代表第i个接收信号的自适应时变滤波器，时变滤波器可以理解为在GNSS有用信号的时频域能量聚集处M(t，u)=1，在干扰区域M(t，u)=0。

　　推荐阅读：无线电工程期刊水平怎么样

转载请注明来自：http://www.qikan2017.com/lunwen/dzi/18505.html