一种基于自适应加窗的分段线性调频干扰消除方法
技术领域
本发明属于无线电接收机干扰抑制领域，涉及一种基于自适应加窗的分段线性调频干扰消除方法。
背景技术
线性调频干扰信号是卫星导航接收机常见的干扰类型，该类干扰可能产生于雷达、或者恶意干扰机，不易被时域或频域干扰抑制技术消除。为了保障GNSS服务的连续性和可靠性，研究人员从信号与干扰在时频域、空域以及空时域的差异着手，一些基于多天线和单天线的干扰抑制方法被提出。其中单天线接收机接收天线的相位中心几乎不变、不存在多阵元通道的幅相特性不一致等误差，定位精度较高，且硬件可空间成本较低，适用于高精度和小型平台。
适用于单天线接收机的变换域干扰检测与抑制方法仍是研究热点。常规单天线方法为将接收信号转换到时频域，检测干扰参数，再利用滤波器或者消隐技术消除干扰成分。典型的时频变换方法有：短时傅里叶变换(Short-Time Fourier Transform,STFT)、小波变换(Wavelet Transform,WT)、wigner-ville分布(Wigner-Ville Distribution,WVD),分数阶傅里叶变换(Fractional Fourier transform,FrFT)。其中FrFT的方法估计精度高，但窗口内出现不同长度不同调频斜率干扰时，会影响分数域聚集程度，同时多次分数域处理会增加有用信号的损失。
因此，如何提升干扰的聚集程度并增加干扰抑制的精确度，进一步增加对于分段线性调频干扰信号的适应能力是亟待解决的技术问题。
发明内容
针对上述现有技术，本发明要解决的技术问题是提供一种适用于单天线接收机的基于自适应加窗的分段线性调频干扰消除方法，在对线性调频干扰抑制过程中，提高干扰在分数阶域的聚集程度，减少干扰与期望信号的交叠程度，在消除干扰成分时对期望信号损伤较小。
为解决上述技术问题，本发明的一种基于自适应加窗的分段线性调频干扰消除方法，包括以下步骤：
步骤1：对单天线接收机接收到的射频信号进行射频处理和数字化处理，获得N个采样点的数字接收信号矢量X；X＝[X(1) X(2),…,X(n),…,X(N)]T，其中n＝1,2,…N，n表示采样时刻，[·]T表示转置；
步骤2：获得信号X的差分图，记为D(n)＝X(n)-X(n-1)，对所得D(n)取其共轭并与自身相乘，获得|D(n)|2＝D(n)·D*(n),(·)*表示共轭，|·|表示模值；
步骤3：|D(n)|2乘一个以q为窗中心的滑动窗口，窗长度为2Q+1，对每个时刻的窗口中的值取中值，即Dm(q)2＝Med[|D(n)|2Wq(n)]，其中
步骤4：对Dm(q)2二次差分获得ΔD(q)2：ΔD(q)2＝Dm(q)2-Dm(q-1)2，检测ΔD(q)2中的序列的周期性；
步骤5：根据获得的干扰估计信息对接收信号进行分段，再消除各段数据中的干扰成分；然后将处理后的各段信号去除残余干扰。
本发明还包括：
1.步骤4中检测D(q)2中的序列的周期性具体为：
步骤4.1：初始化阈值β为设定的门限权值，表示均值，将超过阈值的局部峰值识别为有效峰值，并将这些峰的位置记录为Pi，如果存在Pi-Pi-1＜Hi，Hi为给定的最小周期长度门限，则删除Pi-1；计算获得周期起止位置估计值下角标d＝1,2,···v，d为周期成分标号，v为周期总个数；
步骤4.2：获取线性调频干扰周期长度估计值 并对信号按的位置分段，获得分段后的信号Xd(n)；
步骤4.3：通过离散多项式的方法将局部的Xd(n)信号转换为Md(n)，Md(n)＝Xd(n)Xd *(n-a)，a为设定的时延长度，L为设定的中心提取长度，计算获取线性调频干扰调频率粗估计值
2.步骤5中根据获得的干扰估计信息对接收信号进行分段，再消除各段数据中的干扰成分；然后将处理后的各段信号去除残余干扰具体为：
步骤5.1：在Xd(n)前后添加w/4个零，其中k为可调整数值，当w值超过给定的单次最大处理点数时，调节k使w小于给定值；
步骤5.2：对于补零后的分段信号，加长度为w的布莱克曼窗，窗位置为获得加窗后的信号X'd(n)；
步骤5.3：对于X'd(n)，若d＝1，转到步骤5.4，否则，将调频率估计值与内的估计值依次比较，下角标i为用作对比的成分标号，如果有令pd＝pi，pi为[p1 p2 ...,pi ...,pd-1]中的对应值，转到步骤5.5，否则，则转到步骤5.4；
步骤5.4：通过粗估计最优阶数 根据在给定范围内使用二分法扫描获取最优阶数pd；
步骤5.5：对X'd(n)进行pd阶分数阶傅里叶变换，通过均值统计提取分数域门限值P为给定的门限权值，使用谱线切除的方法去除干扰能量，获得处理后的信号Fd(n)；
步骤5.6：对Fd(n)进行pd阶分数阶傅里叶逆变换，将前后各w/4的样点抛弃获得处理后的信号Yd(n)，并在时域使用谱线切除的方法去除残余脉冲。
本发明的有益效果：本申请发明针对现有适用于单天线卫星导航接收机的线性调频干扰检测与抑制算法对期望卫星信号损伤较大，对于变调频率分段线性调频干扰信号抑制效果不佳的问题，提出一种适用于单天线接收机的基于自适应加窗的分段线性调频干扰消除方法。首先，利用周期线性调频信号的差分图仍是周期函数且在周期内具有单调性的特点，提出一种基于差分图周期参数估计方法，使用局部窗中值处理的方法降低噪声和频率突变对估计结果的影响，估计出接收信号中周期调频分量的调频周期长度、调频起止位置、调频率变化情况；对于分段的接收信号，在前后补零进而根据数据长度自适应加窗处理，保证干扰在分数域的聚集程度，最后通过分数域联合时域完成干扰抑制。本发明能够将不同调频率的干扰分段，防止同一窗口中出现多种调频率的连续调频干扰，降低期望信号与干扰信号的交叠程度，减小周期间相位跳变与多次分数域处理对期望信号的影响。本发明适用单一调频率与变调频率分段线性调频干扰场景。
附图说明
图1基于自适应加窗的变调频率线性调频干扰消除方法原理框图。
具体实施方式
下面结合说明书附图和具体实施方式对本发明做进一步说明。
本发明对接收数据进行差分图干扰参数估计，获取周期调频成分的调频周期长度、调频起止位置与调频率变化情况；根据调频干扰的周期信息将接收的信号进行分段，对分段信号进行补零并自适应地调整分数阶傅里叶变换所需的窗长度以提高干扰在分数阶域的聚集程度，最后通过分数域干扰抑制获得消除干扰后的信号并在时域去除残余干扰。具体实施步骤为：
(1)对单天线接收机接收到的射频信号，经过射频处理和数字化处理，获得N个采样点的数字接收信号矢量X；X＝[X(1) X(2),…,X(n),…,X(N)]T，其中n＝1,2,…N，n表示采样时刻，[·]T表示转置；
(2)获得信号X的差分图，记为D(n)＝X(n)-X(n-1)，其中n为采样点数，对所得D(n)取其共轭并与自身相乘，获得|D(n)|2＝D(n)D*(n)，(·)*表示共轭，|·|表示模值；
(3)|D(n)|2乘一个滑动窗口，窗中心为q，窗长度为2Q+1，对每个时刻的窗口中的值取中值，即Dm(q)2＝Med[|D(n)|2Wq(n)]，其中
(4)检测ΔD(q)2中的序列的周期性，具体步骤如下：
①对Dm(q)2二次差分获得ΔD(q)2：ΔD(q)2＝Dm(q)2-Dm(q-1)2；
②初始化阈值β(12＜β＜16)为门限权值，表示均值，将超过阈值的局部峰值识别为有效峰值，并将这些峰的位置记录为Pi，如果存在Pi-Pi-1＜Hi，Hi为给定的最小周期长度门限，则删除Pi-1；计算获得周期起止位置估计值下角标d＝1,2,···v，d为周期成分标号，v为周期总个数。
③获取线性调频干扰周期长度估计值 并对信号按的位置分段，获得分段后的信号Xd(n)；
④通过离散多项式的方法将局部的Xd(n)信号转换为Md(n)，Md(n)＝Xd(n)Xd *(n-a)，这里a为设定的时延长度，，L为设定的中心提取长度，计算获取线性调频干扰调频率粗估计值
(5)根据获得的干扰估计信息对接收信号进行分段，再消除各段数据中的干扰成分；然后将处理后的各段信号去除残余干扰，具体步骤如下：
①在Xd(n)前后添加w/4个零，其中k为可调整数值；当w值超过给定的单次最大处理点数时，调节k使w小于给定值。
②对于补零后的分段信号，加长度为w的布莱克曼窗，窗位置为获得加窗后的信号X'd(n)；
③对于X'd(n)，若d＝1，转到步骤④，否则，将调频率估计值与内的估计值依次比较，下角标i为用作对比的成分标号，如果有令pd＝pi，pi为[p1 p2 ...,pi ...,pd-1]中的对应值，转到步骤⑤，否则，则转到步骤④；
④通过粗估计最优阶数 根据在给定范围内使用二分法扫描获取最优阶数pd。
⑤对X'd(n)进行pd阶分数阶傅里叶变换，通过均值统计提取分数域门限值P(2<P<4)为门限权值，使用谱线切除的方法去除干扰能量，获得处理后的信号Fd(n)。
⑥对Fd(n)进行pd阶分数阶傅里叶逆变换，将前后各w/4的样点抛弃获得处理后的信号Yd(n)，并在时域使用谱线切除的方法去除残余脉冲。
下面结合具体参数给出具体实施例：
本申请实施例根据周期线性调频干扰的差分图单调性，提出一种基于自适应加窗的分段线性调频干扰消除方法。所申请方法针对分段线性调频干扰在分数域会出现多种聚集峰，影响干扰抑制的问题，用基于差分图的调频周期估计方法对干扰参数进行估计并对信号分组处理。所申请方法可以对分段线性调频干扰分段补零，减少了分数域干扰与期望信号的交叠程度，在消除干扰成分时对期望信号损伤较小。
为更加清晰的说明所申请方法，本申请实施例通过仿真实验进行流程说明以及效果展示，但不限制本申请实施例的范围。实验条件为：有20个周期连续线性调频干扰信号、1个卫星导航信号，经过射频、下变频及数字化处理，卫星导航信号信噪比(SNR)为-20dB，五个连续干扰信号的干噪比(INR)为40dB，干扰信号其他参数如表1所示；接收机的下变频后的中心频率为1.25MHz，数字化采样频率为10.24MHz。周期的最大值设为70μs，最长数据处理长度为26个采样点。
图1是本发明所述方法的原理框图，该方法包括：
S110.数字化接收信号：
x＝[x(1) x(2),…,x(n),…,x(N)]T
其中n＝1,2,…N，表示采样时刻，[·]T表示转置，x(n)表示接收信号模型：
其中jk(n)为第k个周期调频信号，例如d＝1,2,...,v为周期调频信号数，η(n)表示噪声及其他信号。
S120.通过连续差分处理获取差分图估计干扰信号周期信息，具体步骤如下：
①计算D(n)＝X(n)-X(n-1)，乘自身共轭获得|D(n)|2＝D(n)·D*(n)；
②使用滑动窗口取中值Dm(q)2＝Med[|D(n)|2Wq(n)]，其中并对D(q)2二次差分：ΔD(q)2＝Dm(q)2-Dm(q-1)2；
③提取门限(12<β<16)，将超过阈值的局部峰值识别为有效峰值，并将这些峰的位置记录为Pi，如果存在Pi-Pi-1＜50，删除Pi-1。以此获得周期起止位置估计值
④计算周期长度估计值 并对信号按的位置分段，获得分段后的信号Xd(n)；
⑤过离散多项式的方法将局部的Xd(n)信号转换为Md(n)，Md(n)＝Xd(n)Xd *(n-8)，这里计算
S130.对分段信号自适应加窗并去除干扰，具体步骤如下：
①在每段信号前后添加w/4个零，其中
②加长度为w的布莱克曼窗，窗位置为获得加窗后的信号X'd(n)；
③对于X'd(n)，若d＝1，转到步骤④，否则，将调频率估计值与内的估计值依次比较，下角标i为用作对比的成分标号，如果有令pd＝pi，pi为[p1 p2 ...,pi ...,pd-1]中的对应值，转到步骤⑤，否则，则转到步骤④；
④通过粗估计最优阶数 根据在给定范围内使用二分法扫描获取最优阶数pd。
⑤对X'd(n)进行pd阶分数阶傅里叶变换，提取门限P＝2.97，使用谱线切除的方法去除干扰能量，获得Fd(n)；
⑥对Fd进行pd阶分数阶傅里叶逆变换，将前后各w/4的样点抛弃获得处理后的信号Yd(n)，并在时域使用谱线切除的方法去除残余脉冲干扰。
经过干扰消除后，干噪比40dB下输出信号的输出信干噪比(SINRout)、卫星导航信号与原始卫星导航信号的归一化均方误差(NMSE)、卫星信号的相关捕获因子(AF)结果如表2所示，对比方法为传统重叠加窗分数阶傅里叶干扰抑制方法。可以发现本发明所提出方法可以在有效降低干扰抑制处理对期望卫星信号的损伤，并能够有效消除干扰，保证卫星导航接收机在多干扰情况下的工作效能。
综上，本实施例的方法能够根据线性调频信号的差分图特性，使用9个采样点的滑动窗口减少噪声等因素对判别结果的影响，并对含分段线性调频干扰的导航信号进行分段处理。所申请方法可以将一个片段内的干扰能量集中在同一段，减少了分数域最优阶数的搜索次数，同时减少了干扰与期望信号在分数域的交叠程度，在消除干扰成分时对期望信号损伤较小。
表1干扰信号参数
名称
类型
起始频率
调频率
调频周期
1
线性调频
0.6MHz
20GHZ/S
70.7μs
2
线性调频
0.4MHz
30GHZ/S
60.16μs
3
线性调频
0.4MHz
40GHZ/S
60.16μs
4
线性调频
0.1MHz
60GHZ/S
50μs
5
线性调频
0.1MHz
80GHZ/S
50μs
表2本发明所述方法在干噪比40dB抗干扰效果
本领域技术人员可以理解，在本申请具体实施方式的上述方法中，各步骤的序号大小并不意味着执行顺序的先后，各步骤的执行顺序应以其功能和内在逻辑确定，而不应对本申请具体实施方式的实施过程构成任何限定。
最后应说明的是，以上实施例仅用以描述本发明的技术方案而不是对本技术方法进行限制，本发明在应用上可以延伸为其他的修改、变化、应用和实施例，并且因此认为所有这样的修改、变化、应用、实施例都在本发明的精神和教导范围内。
