利用频域检索的单音干扰消除算法

2018-03-07陆炫宇程郁凡唐万斌

信号处理 2018年12期

陆炫宇柏果程郁凡唐万斌

(1. 四川九洲电器集团有限责任公司，四川绵阳 621000；2. 电子科技大学通信抗干扰技术国家级重点实验室，四川成都 611731)

1 引言

在现代无线通信中，正交频分复用(Orthogonal Frequency Division Multiplexing， OFDM)系统以其高频谱效率、灵活的带宽与资源分配等特点，得到了广泛的应用，如第四代移动通信、数字电视地面广播传输系统等[1]。作为一种宽带通信系统，OFDM很容易受到干扰[2]。其中，单音干扰作为一种典型的窄带干扰，可以对OFDM系统的信道误码率造成很恶劣的影响[3]。因此，如何有效抑制在OFDM系统中的单音干扰显得尤为迫切。

单音干扰抑制技术包括干扰频点幅度抑制技术和干扰重构消除技术。其中干扰幅度抑制技术包括干扰频点置零算法[4]，干扰频点限幅算法[5]等。这类技术可以降低单音干扰对系统的影响，却不能改善信道误码率，因此其干扰抑制性能有限。干扰重构消除技术的核心在于干扰参数的估计，其精度直接决定干扰消除后系统的误码性能。通过精确的干扰参数估计与重构，单音干扰可以被基本消除。因此，国内外学者对单音干扰参数估计进行了大量的研究。

单音干扰参数包括幅度、初相和频率。其中频率估计是其他参数估计的基础。单音干扰频率估计算法主要包括相位差分算法和基于离散傅里叶变换(Discrete Fourier Transform， DFT)的算法。其中相位差分算法实现简单，但其估计精度容易受噪声影响，特别是在低干噪比时[6]。基于DFT的算法一般将频率估计分为两步，首先通过DFT得到频率粗估计值，然后通过插值公式得到频率精估计值。文献[7]提出了一种三点插值公式，提高了频率估计精度。文献[8]通过对接收数据补零和一个两点插值公式，进一步提高了频谱分辨率。文献[9]在文献[8]的基础上提出了一种迭代校正算法，提高了频率估计的稳定性。文献[10]结合了加权相位平均(Weighted Phase Average，WPA)算法[11]和DFT，提出了一种DFT-WPA算法，但WPA需要较多的OFDM符号。文献[2]提出了一种DFT-MLE算法，该算法直接对所有接收数据进行DFT，然后进行三点插值，得到精确频率估计值，当OFDM符号较多时，该算法具有较高估计精度，但因其DFT变换点数过多，其计算复杂度较高。

上述频率估计算法只考虑了噪声对频率估计的影响，而没有考虑OFDM调制符号对频率估计的影响。当一次处理多个OFDM符号时，一个数据子载波将对应多个不同的调制符号，根据大数定律，该子载波对应的调制信息近似为一个零均值的噪声。而当一次处理的OFDM符号较少时，特别是只有一个OFDM符号时，大数定律将不能被满足，此时调制符号将对单音干扰参数估计造成较大偏差，从而使干扰消除基本失效，甚至引入新的干扰。为此，本文提出一种利用频域检索的单音干扰消除算法。该算法首先对补零后的接收数据进行DFT，得到干扰频点的粗估计值，然后检索干扰频点附近子载波的调制符号，消除其对单音干扰频域主瓣的影响，然后进行干扰参数估计和干扰消除，根据干扰消除后对应子载波剩余信号与其假设的调制符号组合的欧氏距离，得到最终的调制符号组合，并以此实现精确的单音干扰参数估计和消除。仿真结果表明，新算法可以有效提高单音干扰参数估计精度，实现精确的干扰重构和消除，并可以降低系统误码率。

本文内容安排如下：第2节给出系统模型；第3节介绍现有的一些单音干扰参数估计算法；第4节详细阐述了本文提出的基于频域检索的干扰消除算法；第5节进行了仿真分析；最后对全文进行了总结。

2 信号模型

OFDM信号经过加性高斯白噪声(AWGN，Additive White Gaussian Noise)信道和单音干扰后，基带接收复信号可以表示为

r[n]=s[n]+j[n]+w[n]

(1)

其中s[n]为OFDM信号,j[n]为单音干扰，w[n]为加性高斯白噪声，且w[n]～CN(0，2σ2)。j[n]可以表示为

j[n]=Aej (2πfTsn+φ)

(2)

其中f为干扰频率，A为幅度，φ为初相，Ts为采样间隔且Ts=1/fs,fs为采样率。

假设OFDM系统具有N个子载波，循环前缀长度为C，一次处理L个OFDM符号。为提高频谱分辨率，对去除循环前缀后的每个OFDM符号补N个0，然后分别进行2N点DFT变换，则此时第l个OFDM符号的频域信号可以表示为

Rl[k]=Sl[k]+Jl[k]+Wl[k]

(3)

其中k=0，1，…，2N-1，l=0，1，…，L-1，

(4)

(5)

(6)

其中NC=N+C，φl为第l个OFDM符号的初相。

(7)

(8)

因此，根据式(4)和式(8)，可得

(9)

3 现有单音干扰参数估计算法

本节主要介绍现有的一些单音干扰参数估计算法，包括文献[2]中的DFT-MLE算法，文献[8]中的Fang算法以及文献[10]中的DFT-WPA算法。

3.1 DFT-MLE算法

(10)

(11)

其中

(12)

(13)

其中

(14)

3.2 Fang算法

(15)

(16)

3.3 DFT-WPA算法

(17)

其中

(18)

(19)

4 本文提出的单音干扰消除算法

4.1 新算法原理

(20)

其中n=0，1，…，N-1。故此时可根据式(21)得到干扰消除后的N点频域信号。

(21)

(22)

(23)

(24)

(25)

(26)

(27)

然后根据式(28)，得到调制符号在非子载波频点km-1和km+1产生的信息。

x=-1，1

(28)

(29)

(30)

(31)

(32)

(33)

(34)

其中n=0，1，…，NC-1。最后根据式(35)，得到干扰消除后的接收信号。

(35)

4.2 计算复杂度

最后，对本算法的计算复杂度进行分析。计算复杂度主要包括复数乘法次数。为简化计算量，可提前存储ej 2πγ的值，其中γ=0，1/Λ，2/Λ，…，(Λ-1)/Λ，Λ为一个较大的整数，此时利用查表法，即可得某个角度θ下ej θ的值。

对每一个OFDM符号，计算点2N点频谱需要2N2次复数乘法，由于N一般较大，可忽略式(16)和式(17)的乘法次数，根据式(20)重构干扰需N次复数乘法，根据式(21)得到N点频谱需N2次复数乘法，然后，针对每一种调制符号组合，当km为奇数时，式(23)需N次复数乘法，由于p一般远小于N，可忽略式(25)和式(26)的乘法次数，同理，当km为偶数时，式(28)需2N次复数乘法，综上所述，当km为奇数时，本算法共需约L[3N2+(1+pM)N]次复数乘法，当km为偶数时，本算法共需约L[3N2+(1+2qM)N]次复数乘法。

DFT-MLE算法对所有接收数据进行DFT，故计算频谱需(NCL)2复数乘法，重构干扰需NCL次复数乘法，再次对每个符号进行N点DFT，得到N点频谱共需LN2次复数乘法次数，故其共需约(NCL)2+NCL+LN2次复数乘法。

DFT-WPA算法需L≥2，对每个OFDM符号进行2N点DFT，共需2LN2次复数乘法，WPA算法共需约L(log2L+2)次复数乘法，重构干扰需约LN次复数乘法，再次对每个符号进行N点DFT，得到N点频谱共需LN2次复数乘法次数，故共需3LN2+L(log2L+2+N)。

表1 算法复杂度对比

5 性能仿真分析

本节仿真分析了当OFDM符号数L=1和L=4时，本文所提算法与第3节所介绍算法的干扰抑制度和系统误码率。其中干扰抑制度GJ定义为

(36)

其中PJ为单音干扰时域功率。干信比定义为PJ/Ps，其中Ps为OFDM符号时域功率。

仿真参数如下：N=128，C=0, fs=2000Hz, f～U[0, fs)，φ～U[0，2π)，调制方式为BPSK，OFDM符号时域功率Ps=1，干信比PJ/Ps=1，信道为AWGN信道。

图1为L=1时，不同单音干扰算法干扰抑制度对比。由于L=1时DFT-WPA算法无法计算频率估计值，故此处没有对比该算法。从图1可以看出，本文所提算法干扰抑制度远优于DFT-MLE算法，且随着p与q的增大，性能逐渐提升，当p≥4，q≥3时，则几乎不再有性能提升，这是因为离干扰主瓣越远，干扰的影响越小，检索离干扰主瓣越远的子载波的性能增益也越少。

图1 L=1时，不同单音干扰消除算法干扰抑制度对比

图2为L=1时，不同单音干扰算法误码率对比。从图2中可以看出，本文所提算法和DFT-MLE算法都可以降低误码率，然而DFT-MLE算法在SNR>9 dB之后即出现误码平层，这是因为OFDM符号较少时，其调制符号将对DFT-MLE算法的参数估计会产生较大偏差，且该算法不能消除这种偏差，而本文所提算法没有出现误码平层，且误码率低于DFT-MLE算法，随着p与q的增大，本文算法性能逐渐提升，当p≥4，q≥3时，则几乎不再有性能提升，这与干扰抑制度的分析一致。

图2 L=1时，不同单音干扰消除算法误码率对比

图3为L=4时，不同单音干扰算法干扰抑制度对比。从图3可以看出，当SNR<3 dB时，DFT-MLE算法优于本文算法和DFT-WPA算法，这是因为随着DFT点数增多，DFT-MLE算法性能提升。当SNR>3 dB时，本文所提算法优于DFT-MLE算法和DFT-WPA算法，且随着p与q的增大，性能逐渐提升，当p≥4，q≥3时，性能提升将会变得很小。

图4为L=4时，不同单音干扰算法误码率对比。从图2中可以看出，本文所提算法、DFT-MLE和DFT-WPA算法都可以降低误码率，然而DFT-MLE和DFT-WPA算法在SNR>9 dB之后即出现误码平层，而本文所提算法没有出现误码平层，且误码率低于DFT-MLE和DFT-WPA算法，随着p与q的增大，本文算法性能逐渐提升，当p≥4，q≥3时，性能提升将会变得很小，这与干扰抑制度的分析一致。