一种基于压缩感知的窄带干扰估计和消除方法及装置制造方法
【专利摘要】本发明公开了一种基于压缩感知的窄带干扰估计和消除方法及装置，该方法包括：S1.对一组相邻信号帧的帧头中的L点无帧体干扰部分或L点准无帧体干扰部分进行逐点相减，得到一段时域差分序列，所述时域差分序列长度为L；S2.用时域差分序列构成窄带干扰信号的时域采样序列，根据所述时域采样序列，得到压缩感知算法模型；S3.根据压缩感知算法模型，采用压缩感知算法，估计窄带干扰信号的频率、幅度和相位，得到窄带干扰信号的估计；S4.将接收到的信号帧的帧体数据块进行窄带干扰消除，得到消除窄带干扰信号后的帧体数据。本发明方法能够在各种复杂信道中，不消耗冗余资源，较精确地估计并消除窄带干扰，提升数据传输性能。
【专利说明】一种基于压缩感知的窄带干扰估计和消除方法及装置
【技术领域】
[0001]本发明属于通信【技术领域】，特别涉及一种基于压缩感知的窄带干扰估计方法及装置。
【背景技术】
[0002]在宽带高速数据传输中，存在频率选择性衰落、时间选择性衰落和窄带噪声干扰(Narrowband Interference, NBI)等，影响数据传输质量。
[0003]目前，正交频分复用(Orthogonal Frequency Division Multiplex, 0FDM)技术在电力线通信系统、数字电视广播系统中都得到了广泛的应用。由于OFDM具有很好的频率选择性，它已经应用到各种数字信号传输系统中，如国际电联的电力线通信系统标准、无线局域网、欧洲数字视频地面广播以及中国地面数字电视传输标准(Digital TelevisionMultimedia Broadcast, DTMB)等。其中，DTMB系统采用了时域同步正交频分复用技术TDS-OFDM(Time Domain Synchronous OFDM),使用时域训练序列作为信号巾贞的保护间隔填充，时域训练序列还可用于同步和信道估计等。DTMB系统的多层复帧结构如图1所示。
[0004]在传输便捷、资源丰富的电力线信道中进行高效数字通信得到了广泛的研究和应用，然而电力线信道环境恶劣，尤其是有严重的窄带噪声干扰；地面数字电视广播系统中也存在窄带干扰等问题。同时，在TDS-OFDM系统中，窄带干扰会降低信道估计准确性，影响数据的正确解映射和解码，严重影响数据的正确传输，因此，需要设计对抗窄带干扰的方法。
[0005]现有的应对窄带干扰的方法性能不够理想，尤其在较高强度的窄带干扰、较严重的信道多径衰落或较强的信道噪声下，数据的传输性能会严重恶化。如传统的频域交织等手段，虽然可以一定程度上降低窄带干扰的影响，但是无法消除窄带干扰，在窄带干扰下效果不佳。现有的基于帧头结构优化设计对抗窄带干扰对同步性能的影响的方法，也不可以消除窄带干扰。目前常见的基于线性预测算法估计窄带干扰的方法，需要依赖于某种冗余资源，例如插入一定数量的虚拟子载波，或者必须进行过采样，或者在接收端需要使用多天线等；或者依赖于接收数据的二阶统计量进行窄带干扰估计。现有方法估计窄带干扰所要占用的冗余资源降低了资源利用效率，而依赖于接收数据的统计信息的方法则使窄带干扰估计性能受到接收数据统计信息准确性的影响，导致估计性能不佳。
【发明内容】
[0006](一)要解决的技术问题
[0007]本发明所要解决的技术问题是现有的窄带干扰信号估计方法在较高强度的窄带干扰、较严重的信道多径衰落或较强的信道噪声下，无法准确估计窄带干扰信号的频率和大小，以及现有窄带干扰估计方法必须占用冗余资源的问题。
[0008](二)技术方案
[0009]为解决上述问题，本发明提供一种基于压缩感知的窄带干扰估计和消除方法，该方法包括:[0010]S1.对一组相邻信号帧的帧头中的L点无帧体干扰部分或L点准无帧体干扰部分进行逐点相减，得到一段时域差分序列，所述时域差分序列长度为L ；
[0011]S2.用时域差分序列构成窄带干扰信号的时域采样序列，根据所述时域采样序列，得到压缩感知算法模型；
[0012]S3.根据压缩感知算法模型，采用压缩感知算法，估计窄带干扰信号的频率、幅度和相位，得到窄带干扰信号的估计；
[0013]S4.将接收到的信号帧的帧体数据块进行窄带干扰消除，得到消除窄带干扰信号后的帧体数据。
[0014]其中，在步骤SI中，所述信号帧包括帧头和帧体，其中，帧头为训练序列，帧体为单载波数据块或OFDM数据块，所述训练序列包括一段已知的训练序列、一段已知的训练序列及其循环前缀或两段相同的已知的训练序列；所述已知的训练序列包括频域二值伪随机序列的离散傅里叶逆变换或时域二值伪随机序列。
[0015]其中，在步骤SI中，所述无帧体干扰部分或准无帧体干扰部分包括整个帧头训练序列、包含部分信号帧间多径干扰区域的训练序列或两段相同已知训练序列构成的帧头中的第二段已知训练序列。
[0016]其中，在步骤S2中，所述窄带干扰信号的时域采样序列包括当前信号帧与下一信号帧进行逐点相减所得的时域差分序列。
[0017]其中，在步骤S2中，所述压缩感知算法模型为窄带干扰信号时频关系等式，即所述窄带干扰时域采样序列等于傅里叶逆变换矩阵乘以所述待估计的窄带干扰频域信号加上时域噪声信号。
[0018]其中，在步骤S2中，所述压`缩感知算法为凸优化算法或者贪婪算法，其中，所述凸优化算法包括内点法、一阶范数最小化算法；所述贪婪算法包括压缩采样匹配追踪法、正交匹配追踪法及其衍生算法、稀疏自`适应匹配追踪法以及子空间追踪法。
[0019]较佳的，在步骤S3和步骤S4之间，该方法还包括:
[0020]S31.根据得到的窄带干扰信号估计，将傅里叶逆变换矩阵乘以窄带干扰信号的估计，再计算其与所述窄带干扰信号的时域采样序列的残差平方，得到残差平方式，对残差平方式进行最小二乘计算，得到频域窄带干扰信号的估计；
[0021]较佳的，在步骤S31和步骤S4之间，该方法还包括:
[0022]S32.将所述频域窄带干扰信号的估计除以频域衰减因子，其中第k个子载波上的频域衰减因子为:
【权利要求】
1.一种基于压缩感知的窄带干扰估计和消除方法，其特征在于，该方法包括: 51.对一组相邻信号帧的帧头中的L点无帧体干扰部分或L点准无帧体干扰部分进行逐点相减，得到一段时域差分序列，所述时域差分序列长度为L ； 52.用时域差分序列构成窄带干扰信号的时域采样序列，根据所述时域采样序列，得到压缩感知算法模型； 53.根据压缩感知算法模型，采用压缩感知算法，估计窄带干扰信号的频率、幅度和相位，得到窄带干扰信号的估计； 54.将接收到的信号帧的帧体数据块进行窄带干扰消除，得到消除窄带干扰信号后的帧体数据。
2.根据权利要求1所述的方法，其特征在于，在步骤SI中，所述信号帧包括帧头和帧体，其中，帧头为训练序列，帧体为单载波数据块或OFDM数据块，所述训练序列包括一段已知的训练序列、一段已知的训练序列及其循环前缀或两段相同的已知的训练序列；所述已知的训练序列包括频域二值伪随机序列的离散傅里叶逆变换或时域二值伪随机序列。
3.根据权利要求1所述的方法，其特征在于，在步骤SI中，所述无帧体干扰部分或准无帧体干扰部分包括整个帧头训练序列、包含部分信号帧间多径干扰区域的训练序列或两段相同已知训练序列构成的帧头中的第二段已知训练序列。
4.根据权利要求1所述的方法，其特征在于，在步骤S2中，所述窄带干扰信号的时域采样序列包括当前信号帧 与下一信号帧进行逐点相减所得的时域差分序列。
5.根据权利要求1所述的方法，其特征在于，在步骤S2中，所述压缩感知算法模型为窄带干扰信号时频关系等式，即所述窄带干扰时域采样序列等于傅里叶逆变换矩阵乘以所述待估计的窄带干扰频域信号加上时域噪声信号。
6.根据权利要求1所述的方法，其特征在于，在步骤S2中，所述压缩感知算法为凸优化算法或者贪婪算法，其中，所述凸优化算法包括内点法、一阶范数最小化算法；所述贪婪算法包括压缩采样匹配追踪法、正交匹配追踪法及其衍生算法、稀疏自适应匹配追踪法以及子空间追踪法。
7.根据权利要求1所述的方法，其特征在于，在步骤S3和步骤S4之间，该方法还包括: 531.根据得到的窄带干扰信号估计，将傅里叶逆变换矩阵乘以窄带干扰信号的估计，再计算其与所述窄带干扰信号的时域采样序列的残差平方，得到残差平方式，对残差平方式进行最小二乘计算，得到频域窄带干扰信号的估计。
8.根据权利要求7所述的方法，其特征在于，在步骤S31和步骤S4之间，该方法还包括: 532.将所述频域窄带干扰信号的估计除以频域衰减因子，其中第k个子载波上的频域衰减因子为:
9.根据权利要求1所述的方法，其特征在于，在步骤S4中，所述窄带干扰消除包括:将接收到的信号帧的时域帧体数据块减去所述频域窄带干扰信号估计的离散傅里叶逆变换，得到消除窄带干扰信号后的帧体数据。
10.一种基于压缩感知的窄带干扰估计装置，其特征在于，该装置包括: 帧头差分模块，用于对一组相邻信号帧的帧头中的L点无帧体干扰部分或L点准无帧体干扰部分进行逐点相减，得到一段时域差分序列，所述序列长度为L ； 压缩感知估计模块，用于用时域 差分序列构成窄带干扰信号的时域采样序列，根据所述时域采样序列，得到压缩感知算法模型，并根据压缩感知算法模型，采用压缩感知算法，估计窄带干扰信号的频率、幅度和相位，得到窄带干扰信号的估计。 干扰消除模块，用于对帧体数据块进行窄带干扰消除，得到消除窄带干扰后的帧体数据。
【文档编号】H04L25/08GK103763227SQ201410006974
【公开日】2014年4月30日 申请日期:2014年1月7日 优先权日:2014年1月7日
【发明者】杨昉, 刘思聪, 宋健, 潘长勇 申请人:清华大学