CN103986562B - 抑制同频干扰的方法和装置 - Google Patents

抑制同频干扰的方法和装置 Download PDF

Info

Publication number
CN103986562B
CN103986562B CN201410245825.XA CN201410245825A CN103986562B CN 103986562 B CN103986562 B CN 103986562B CN 201410245825 A CN201410245825 A CN 201410245825A CN 103986562 B CN103986562 B CN 103986562B
Authority
CN
China
Prior art keywords
msub
mrow
data point
predeterminable area
msup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410245825.XA
Other languages
English (en)
Other versions
CN103986562A (zh
Inventor
汪浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Honor Device Co Ltd
Original Assignee
Huawei Technologies Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Priority to CN201410245825.XA priority Critical patent/CN103986562B/zh
Publication of CN103986562A publication Critical patent/CN103986562A/zh
Priority to PCT/CN2015/080498 priority patent/WO2015184967A1/zh
Application granted granted Critical
Publication of CN103986562B publication Critical patent/CN103986562B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/08—Modifications for reducing interference; Modifications for reducing effects due to line faults ; Receiver end arrangements for detecting or overcoming line faults

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Noise Elimination (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明实施例提供一种抑制同频干扰的方法和装置，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。具体方案为：首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对接收信号y进行干扰抑制合并。本发明的实施例用于抑制同频干扰。

Description

抑制同频干扰的方法和装置
技术领域
本发明实施例涉及通信技术领域，尤其涉及一种抑制同频干扰的方法和装置。
背景技术
随着无线宽带通信技术的发展,用户对通信系统的性能提出更高要求，比如更高的峰值吞吐量和平均吞吐量，使得位于小区边缘的用户也能获得质量高的通信话音和数据业务服务。第3代合作伙伴计划(3rd Generation Partnership Project,3GPP)开始了通用移动通信系统(Universal Mobile Telecommunications System,UMTS)技术的长期演进(Long Term Evolution,LTE)系统。多输入多输出(Multiple-Input Multiple-Output，MIMO)和正交频分复用(Orthogonal Frequency Division Multiplex,OFDM)是LTE最关键的两个技术。
在LTE实际场景中，用户设备(User Equipment,UE)可能受到邻小区被调度UE相同时间、频率资源上的干扰，会严重降低UE解调数据信道的性能。基站高密度和异构是LTE网络结构的演进大方向，这显然会带来更加严重的小区间同频干扰，一般采用干扰抑制合并(Interference Rejection Combining,IRC)等算法来抑制干扰，IRC算法利用导频点或数据点的相关数据估计干扰信号的统计特性并基于该统计特性进行干扰白化以减少干扰信号对UE解调数据信道性能的影响。
IRC算法的关键在于准确估计干扰信号的统计特性也即准确计算干扰自相关矩阵，但是现有技术中，计算干扰自相关矩阵时由于要利用解调译码后再进行反编码得到的数据，在中低信噪比(Signal-to-Noise Ratio,SNR)时反编码得到的数据准确度很差，这就使得计算出的干扰自相关矩阵会出现很大偏差，从而无法有效抑制干扰。
发明内容
本发明实施例提供一种抑制同频干扰的方法和装置，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
第一方面，本发明的实施例提供一种抑制同频干扰的方法，所述方法包括：
利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对所述预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵其中，所述预设区域包括预设个数的资源块，所述资源块中包含有所述数据点和所述导频点；
对所述和进行多入多出MIMO检测解调，获取所述数据点的比特对数似然比；
基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值；
基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'；
利用所述Ruu'对所述接收信号y进行干扰抑制合并。
结合第一方面，在第一种可能的实现方式中，所述获取所述数据点的比特对数似然比还包括：
对经过所述检测解调后获取的所述数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将所述新的比特对数似然比作为所述数据点的比特对数似然比。
结合第一方面或第一方面的第一种可能的实现方式，在第二种可能的实现方式中，所述基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值包括：
根据均值计算公式计算所述数据点承载的数据的均值；
所述均值计算公式包括：
其中，E(xj)为所述数据点承载的数据的均值，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)。
结合第一方面或第一方面的第一种或第二种可能的实现方式，在第三种可能的实现方式中，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'包括：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
结合第一方面或第一方面的第一种或第二种可能的实现方式，在第四种可能的实现方式中，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'还包括：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
结合第一方面，在第五种可能的实现方式中，在所述获取所述数据点的比特对数似然比之后还包括：
计算所述第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当所述比值大于第一预设阈值，则需要挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当所述比值小于或等于所述第一预设阈值，则任意选取所述预设区域内的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
结合第一方面的第五种可能的实现方式，在第六种可能的实现方式中，所述挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点包括：
将获取的所述预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，所述每个数据点承载的数据的方差包括预设层数数据的方差；
若所述数据点任一层数据的方差大于所述第二预设阈值，则删除所述数据点；
若所述数据点每层数据的方差小于或等于所述第二预设阈值，则确定所述数据点满足预设条件作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
结合第一方面，在第七种可能的实现方式中，在所述获取所述数据点的比特对数似然比之前还包括：
将用户设备的秩与接收天线数进行比较；
若所述秩小于所述接收天线数，则执行所述获取所述数据点的比特对数似然比；
若所述秩大于或等于所述接收天线数，则利用所述Ruu对所述接收信号y进行干扰抑制合并。
第二方面，本发明的实施例提供一种抑制同频干扰的装置，所述装置包括：
第一抑制单元，用于利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对所述预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵其中，所述预设区域包括预设个数的资源块，所述资源块中包含有所述数据点和所述导频点；
第一获取单元，用于对所述和进行多入多出MIMO检测解调，获取所述数据点的比特对数似然比；
第二获取单元，用于基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值；
第三获取单元，用于基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'；
第二抑制单元，用于利用所述Ruu'对所述接收信号y进行干扰抑制合并。
结合第二方面，在第一种可能的实现方式中，所述第一获取单元还用于：
对经过所述检测解调后获取的所述数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将所述新的比特对数似然比作为所述数据点的比特对数似然比。
结合第二方面或第二方面的第一种可能的实现方式，在第二种可能的实现方式中，所述第二获取单元具体用于：
根据均值计算公式计算所述数据点承载的数据的均值；
所述均值计算公式包括：
其中，E(xj)为所述数据点承载的数据的均值，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)。
结合第二方面或第二方面的第一种或第二种可能的实现方式，在第三种可能的实现方式中，所述第三获取单元具体用于：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
结合第二方面或第二方面的第一种或第二种可能的实现方式，在第四种可能的实现方式中，所述第三获取单元还具体用于：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
结合第二方面的第四种可能的实现方式，在第五种可能的实现方式中，所述装置还包括：
第一判断单元，用于计算所述第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当所述比值大于第一预设阈值，则需要挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当所述比值小于或等于所述第一预设阈值，则任意选取所述预设区域内的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
结合第二方面的第五种可能的实现方式，在第六种可能的实现方式中，所述装置还包括：
挑选单元，用于挑选所述预设区域内满足预设条件的数据点作为所述第三获取单元使用的数据点，所述挑选所述预设区域内满足预设条件的数据点作为所述第三获取单元使用的数据点包括：
将获取的所述预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，所述每个数据点承载的数据的方差包括预设层数数据的方差；
若所述数据点任一层数据的方差大于所述第二预设阈值，则删除所述数据点；
若所述数据点每层数据的方差小于或等于所述第二预设阈值，则确定所述数据点满足预设条件作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
结合第二方面，在第七种可能的实现方式中，所述装置还包括：
第二判断单元，用于在获取所述数据点的比特对数似然比之前，将用户设备的秩与接收天线数进行比较；
若所述秩小于所述接收天线数，则执行所述获取所述数据点的比特对数似然比；
若所述秩大于或等于所述接收天线数，则利用所述Ruu对所述接收信号y进行干扰抑制合并。
本发明实施例提供的抑制同频干扰的方法和装置，首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对接收信号y进行干扰抑制合并。这样，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的实施例提供的抑制同频干扰的方法的流程示意图一；
图2为本发明的实施例提供的抑制同频干扰的方法的流程示意图二；
图3为本发明的实施例中获取预设区域内的数据点的比特对数似然比的流程示意图；
图4为本发明的实施例中获取预设区域内的数据点的比特对数似然比的效果示意图；
图5为本发明的实施例中获取预设区域内的第二干扰自相关矩阵Ruu'的效果示意图；
图6为本发明的实施例提供的抑制同频干扰的装置的结构示意图一；
图7为本发明的实施例提供的抑制同频干扰的装置的结构示意图二；
图8为本发明的实施例提供的抑制同频干扰的装置的结构示意图三。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供一种抑制同频干扰的方法，如图1所示，该方法包括：
步骤101、利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵
其中，预设区域包括预设个数的资源块，资源块中包含有数据点和导频点。
步骤102、对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比。
步骤103、基于均值算法，根据比特对数似然比计算数据点承载的数据的均值。
步骤104、基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'。
步骤105、利用Ruu'对接收信号y进行干扰抑制合并。
其中，利用干扰自相关矩阵对数据点的接收信号、等效信道矩阵进行干扰抑制合并可以是对数据点的接收信号、等效信道矩阵进行干扰白化，干扰白化只是干扰抑制合并的其中一种实现方式，本发明实施例对干扰抑制合并的具体实现方式不做限定，本领域技术人员可以根据实际选择具体的实现方式。
本发明实施例提供的抑制同频干扰的方法，首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对所述接收信号y进行干扰抑制合并。这样，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
为了使本领域技术人员能够更清楚地理解本发明实施例提供的技术方案，下面通过具体的实施例，对本发明的实施例提供的抑制同频干扰的方法进行详细说明，如图2所示，该方法包括：
步骤201、将用户设备的秩与接收天线数进行比较。
若用户设备的秩大于或等于接收天线数，则执行步骤202；
若用户设备的秩小于接收到天线数，则执行步骤203及后续步骤。
示例性的，UE的秩(Rank Index，RI)，即数据层数。若UE的秩小于UE的接收天线数，则执行步骤202，步骤202实际上就是现有技术中IRC算法使用预设区域内导频点数据估算预设区域内干扰自相关矩阵并用得到的干扰自相关矩阵抑制干扰。当RI大于或等于接收天线数时，利用本发明实施例提供的计算干扰自相关矩阵相比现有技术所获得的准确度的增益是比较有限的。在UE复杂度受限的时候，可以在RI大于或等于接收天线数的情况下选择使用现有技术的方案。
步骤202、利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y进行干扰抑制合并。
其中，预设区域包括预设个数的资源块(Resource Block，RB)，该资源块RB在时域上包含连续的若干个OFDM符号，在频域上包含连续的若干个子载波，每个RB中包括若干个资源单元(Resource Element，RE)，RE是导频点或数据点，RB的个数以及其中导频点和数据点的个数根据实际需要选取，本发明对此不作限定。
具体的，步骤202实际上就是现有技术中抑制同频干扰的方案，下面进行示例性的说明：
给定一个LTE系统，假设UE侧接收天线个数为N，用x、y和H分别表示某个RE的发射符号向量、接收向量和等效频域信道矩阵，那么信道模型可以表示为：y＝Hx+u，其中u是N维向量，表示UE受到的同频干扰和加性高斯白噪声。
计算预设区域的第一干扰自相关矩阵Ruu的公式如下：
其中，yk和Hk分别为该区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，Lp是导频点个数。需要说明的是，上述公式中，计算该预设区域内所有导频点的干扰自相关矩阵，然后取平均值作为该区域内导频点的第一干扰自相关矩阵。
对上述Ruu进行Cholesky分解，即Ruu＝LLH，其中L是下三角阵。
然后利用L-1对数据点的接收信号y进行干扰抑制合并，例如，可以对接收信号y(叠加有干扰信号)滤波也即对y中的干扰信号进行干扰白化，得到：
即接收信号y中干扰信号的相关性被去除，其中为干扰白化后的接收信号，上述干扰白化只是干扰抑制合并的其中一种实现方式，本发明实施例对干扰抑制合并的具体实现方式不做限定，本领域技术人员可以根据实际选择具体的实现方式。
步骤203、获取数据点的比特对数似然比。
具体的，如图3所示，获取数据点的比特对数似然比的步骤包括：
步骤203a、获取预设区域内的导频点的第一干扰自相关矩阵Ruu。
示例性的，按照步骤202中给出的方法计算预设区域内的导频点的第一干扰自相关矩阵Ruu。
步骤203b、利用Ruu对数据点的接收信号y和等效信道矩阵H进行干扰抑制合并得到干扰抑制合并后的接收信号和等效频域信道矩阵
示例性的，按照步骤202给出的方法对预设区域内数据点的接收信号y和等效信道矩阵H做干扰白化。
对Ruu进行Cholesky分解，即Ruu＝LLH，利用L-1对数据点的接收信号y干扰白化，得到 为干扰白化后的接收信号；利用L-1对数据点的等效信道矩阵H干扰白化，得到 为干扰白化后的等效信道矩阵。
步骤203c、对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比。
示例性的，根据步骤203b得到干扰白化后的数据点的接收信号和等效信道矩阵然后经过MIMO检测解调后得到比特对数似然比(Byte Log Likelihood Ratio，BLLR)。
其中，MIMO检测解调获取比特对数似然比可以使用多种方法，以下为一种实现算法：
例如，基于最大似然准则，xj的第i个比特bj,i的比特对数似然比γj,i为：
其中，Pr(bj,i)为对应xj的第i个比特bj,i的概率，xj为数据点的发射符号，为干扰抑制合并后的接收信号，为干扰抑制合并后的等效频域信道矩阵。
经过步骤203a、203b以及步骤203c已经得到比特对数似然比，为了获得更准确的比特对数似然比还可以对步骤203c得到的比特对数似然比进行进一步处理。
步骤203d、对经过检测解调后获取的数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比。
示例性的，利用检测解调后得到的比特对数似然比，经过软值速率匹配、软值编码块(Coded Block，CB)级联和软值加扰(和传统编码中的相关操作基本相同，区别仅在输入从0/1比特变成了比特对数似然比)操作后得到长度/顺序和解调后得到的比特对数似然比完全匹配的新的比特对数似然比。
步骤203e、将新的比特对数似然比作为数据点的比特对数似然比。
具体的，将步骤203d获取的新的比特对数似然比作为步骤203要获取的数据点的比特对数似然比。
示意性的，图4所示为上述两种获取比特对数似然比的方式。
需要说明的是，经过步骤203a、步骤203b以及步骤203c后得到的比特对数似然比与经过步骤203a至步骤203e后得到的比特对数似然比相比，采用后一种方式得到的比特对数似然比准确度更高些。
步骤204、计算导频点的第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值。
示例性的，对于根据步骤203a已经得到的矩阵Ruu，计算Ruu对角线元素模的平均值与非对角线元素模的平均值的比值的公式如下：
其中，ci和di分别是Ruu中第i个对角线和非对角线元素，Nd和Nnd分别是对角线和非对角线元素的个数。
将Ratio和第一预设阈值Ts比较，若Ratio＞Ts，则需要挑选预设区域内满足预设条件的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点，执行步骤205、步骤206、步骤207或步骤208、步骤209，其中步骤206为挑选数据点的步骤；若Ratio≤Ts，则任意选取预设区域内的数据点数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点，执行步骤205、步骤207或步骤208、步骤209。
步骤205、基于均值算法，根据比特对数似然比计算数据点承载的数据的均值以及根据该均值计算数据点承载的数据的方差。
示例性的，根据比特对数似然比以及均值计算公式和方差计算公式计算数据点承载的数据的均值与方差；
均值计算公式包括：
方差计算公式包括：
其中，E(xj)为数据点承载的数据的均值，Var(xj)为数据点承载的数据的方差，xj为数据点的发射符号，Pr(bj,i)为对应xj的第i个比特bj,i的概率，根据比特对数似然比计算Pr(bj,i)。
其中，计算Pr(bj,i)的公式为：和γj,i是比特bj,i的比特对数似然比。
步骤206、挑选预设区域内满足预设条件的数据点作为后续计算预设区域内第二干扰自相关矩阵Ruu'时使用的数据点。
具体的，将获取的预设区域的每个数据点承载的数据的方差与第二预设阈值比较，其中，每个数据点承载的数据的方差包括预设层数数据的方差，预设层数即UE的数据层数；若该数据点任一层数据的方差大于第二预设阈值，则确定该数据点不满足预设条件；若该数据点每层数据的方差小于或等于第二预设阈值，则确定该数据点满足预设条件；若该数据点满足预设条件则作为后续步骤207或步骤208计算预设区域内第二干扰自相关矩阵Ruu'时使用的数据点；若该数据点不满足预设条件则在后续步骤207或步骤208计算预设区域内第二干扰自相关矩阵Ruu'时不使用该数据点。
示例性的，将根据步骤205计算得到的数据点的每层数据的方差Var(xj)与第二预设阈值Tvar比较，若任一层数据的方差Var(xj)＞Tvar，则该数据点不满足预设条件，若每层数据的方差Var(xj)≤Tvar，则认为该数据点满足预设条件。
步骤207、将均值作为第一协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'。
具体的，其中，第一协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
步骤208、将均值与方差作为第二协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'。
具体的，其中，第二协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
示意性的，图5所示为上述获取预设区域内的第二干扰自相关矩阵Ruu'的效果示意图。
需要说明的是，现有技术中当译码质量很差，计算干扰自相关矩阵的时候，反编码错误的符号与真实传输符号的误差一般较大，而且计算干扰自相关矩阵未考虑误差的存在。在本发明的实施例中，用均值替代真实传输信号，因为是平均值，误差会小很多；另外，在计算干扰自相关矩阵时还可以通过引入方差进一步减小误差。这样，通过本发明实施例提供的方法计算的干扰自相关矩阵相对现有技术计算的干扰自相关矩阵更加准确。
步骤209、利用Ruu'对接收信号y进行干扰抑制合并。
其中，对接收信号y进行干扰抑制合并可以是对接收信号y进行干扰白化，步骤209可以采用本领域通用的利用干扰自相关矩阵进行干扰抑制合并，干扰白化只是给出的其中一种实现方式，本发明对干扰抑制合并的具体实现方式不做限定。
示例性的，对根据步骤207或步骤208得到的Ruu'进行Cholesky分解，即Ruu'＝LLH，其中L是下三角阵。
然后可以利用L-1对预设区域内的干扰信号进行干扰抑制合并，例如，可以对导频点的接收信号y(叠加有干扰信号)滤波也即对y中的干扰信号进行干扰白化，得到：
即接收信号y中干扰信号的相关性被去除，其中为干扰白化后的接收信号。
需要说明的是，本发明的技术方案，可以应用于各种通信系统，例如：全球移动通信系统(Global System for Mobile Communications，GSM)、码分多址(Code DivisionMultiple Access，CDMA)系统、宽带码分多址(Wideband Code Division MultipleAccess，WCDMA)系统、时分同步码分多址(Time Division-Synchronous Code DivisionMultiple Access，TD-SCDMA)系统、LTE系统等。本发明对此不做限定。
本发明实施例提供的抑制同频干扰的方法，首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对接收信号y进行干扰抑制合并。这样，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
本发明的实施例提供一种抑制同频干扰的装置00，如图6所示，该装置00包括：第一抑制单元10、第一获取单元20、第二获取单元30、第三获取单元40以及第二抑制单元50。
第一抑制单元10，用于利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵
其中，预设区域包括预设个数的资源块，该资源块中包含有数据点和导频点。
第一获取单元20，用于对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比。
示例性的，预设区域包括预设个数的资源块，该资源块RB在(时域上包含连续的若干个OFDM符号，在频域上包含连续的若干个子载波)，每个RB中包括若干个资源单元(Resource Element，RE)，RE是导频点或数据点，RB的个数以及其中导频点和数据点的个数根据实际需要选取，本发明对此不作限定。
给定一个LTE系统，假设UE侧接收天线个数为N，用x、y和H分别表示某个RE的发射符号向量、接收向量和等效频域信道矩阵，那么信道模型可以表示为：y＝Hx+u，其中u是N维向量，表示UE受到的同频干扰和加性高斯白噪声。
计算预设区域的导频点的第一干扰自相关矩阵Ruu的公式如下：
其中，yk和Hk分别为该区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，Lp是导频点个数。需要说明的是，上述公式中，计算该预设区域内所有导频点的干扰自相关矩阵，然后取平均值作为该区域内导频点的第一干扰自相关矩阵。
对上述Ruu进行Cholesky分解，即Ruu＝LLH，利用L-1对数据点的接收信号y干扰白化，得到 为干扰白化后的接收信号；利用L-1对数据点的等效信道矩阵H干扰白化，得到 为干扰白化后的等效信道矩阵。
根据上述得到的干扰白化后的数据点的接收信号和等效信道矩阵然后经过MIMO检测解调后得到比特对数似然比。
其中，MIMO检测解调获取比特对数似然比可以使用多种方法，以下为一种实现算法：
例如，基于最大似然准则，xj的第i个比特bj,i的比特对数似然比γj,i为：
其中，Pr(bj,i)为对应xj的第i个比特bj,i的概率，xj为数据点的发射符号，为干扰抑制合并后的接收信号，为干扰抑制合并后的等效频域信道矩阵。
第二获取单元30，用于基于均值算法，根据比特对数似然比计算数据点承载的数据的均值。
第三获取单元40，用于基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'。
第二抑制单元50，用于利用Ruu'对接收信号y进行干扰抑制合并。
可选的，第一获取单元20还可以具体用于：
对经过检测解调后获取的数据点的比特对数似然比进行译码处理、软值速率匹配处理、软值编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将新的比特对数似然比作为数据点的比特对数似然比。
示例性的，利用检测解调后得到的比特对数似然比，经过软值速率匹配、软值编码块级联和软值加扰(和传统编码中的相关操作基本相同，区别仅在输入从0/1比特变成了比特对数似然比)操作后得到长度/顺序和解调后得到的比特对数似然比完全匹配的新的比特对数似然比，将新的比特对数似然比作为第一获取单元20要获取的预设区域内数据点的比特对数似然比。
可选的，第二获取单元30可以具体用于：
根据均值计算公式计算数据点承载的数据的均值；
均值计算公式包括：
其中，E(xj)为数据点承载的数据的均值，xj为数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据比特对数似然比计算Pr(bj,i)。
示例性的，其中，计算Pr(bj,i)的公式为：和γj,i是比特bj,i的比特对数似然比。
可选的，第三获取单元40可以具体用于：
将均值作为第一协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'；
第一协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的第二干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
可选的，第三获取单元40还可以具体用于：
根据均值以及方差计算公式计算数据点承载的数据的方差；
方差计算公式包括：
其中E(xj)为数据点承载的数据的均值，Var(xj)为数据点承载的数据的方差，xj为数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据比特对数似然比计算Pr(bj,i)；
将均值与方差作为第二协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'；
第二协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的第二干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
可选的，如图7所示，该装置00还包括：
第一判断单元60，用于计算第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当比值大于第一预设阈值，则需要挑选预设区域内满足预设条件的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当比值小于或等于第一预设阈值，则任意选取预设区域内的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点。
示例性的，根据第一干扰自相关矩阵Ruu计算Ruu对角线元素模的平均值与非对角线元素模的平均值的比值，公式如下：
其中，ci和di分别是Ruu中第i个对角线和非对角线元素，Nd和Nnd分别是对角线和非对角线元素的个数。
将Ratio和第一预设阈值Ts比较，若Ratio＞Ts，则需要挑选预设区域内满足预设条件的数据点作为第三获取单元40中使用的数据点；若Ratio≤Ts，则任意选取预设区域内的数据点作为第三获取单元40中使用的数据点。
挑选单元70，用于挑选预设区域内满足预设条件的数据点作为第三获取单元40中使用的数据点。
具体的，将获取的预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，每个数据点承载的数据的方差包括预设层数数据的方差；
若该数据点任一层数据的方差大于第二预设阈值，则确定该数据点不满足预设条件剔除该数据点；
若该数据点每层数据的方差小于或等于第二预设阈值，则确定该数据点满足预设条件作为第三获取单元40计算第二干扰自相关矩阵Ruu'时使用的数据点。
示例性的，将计算得到的数据点的每层数据的方差Var(xj)与第二预设阈值Tvar比较，若任一层数据的方差Var(xj)＞Tvar，则该数据点不满足预设条件，若每层数据的方差Var(xj)≤Tvar，则认为该数据点满足预设条件。
第二判断单元80，用于在获取数据点的比特对数似然比之前，将用户设备的秩与接收天线数进行比较；
若秩小于接收天线数，则执行获取数据点的比特对数似然比；
若秩大于或等于接收天线数，则利用Ruu对接收信号y进行干扰抑制合并。
示例性的，UE的秩，即数据层数。若UE的秩小于UE的接收天线数，则利用现有技术中IRC算法使用预设区域内导频点数据估算预设区域内干扰信号协方差矩阵并用得到的协方差矩阵抑制干扰。这么做的原因是：当RI大于或等于接收天线数时，利用本发明实施例提供的计算干扰信号协方差矩阵相比现有技术所获得的准确度的增益是比较有限的。在UE复杂度受限的时候，可以在RI大于或等于接收天线数的情况下选择使用现有技术的方案。
本实施例用于实现上述各方法实施例，本实施例中各个单元的工作流程和工作原理参见上述各方法实施例中的描述，在此不再赘述。
另外，上述实施例中的装置可以集成在基于各种通信系统的所有同频干扰下的多天线接收机中，可以在基站的接收机中也可以在UE的接收机中，本发明对此不作限制。
本发明实施例提供的抑制同频干扰的装置，首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对接收信号y进行干扰抑制合并。这样，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
本发明实施例还提供了一种抑制同频干扰的装置90，如图8所示，该装置90包括：总线94；以及连接到总线94的处理器91、存储器92和接口93，其中该接口93用于通信；该存储器92用于存储指令，处理器91用于执行该指令用于：
利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵其中，预设区域包括预设个数的资源块，资源块中包含有数据点和导频点；
对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；
基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；
基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；
利用Ruu'对接收信号y进行干扰抑制合并。
可选地，处理器91执行该指令用于获取数据点的比特对数似然比，还可以包括：
对经过检测解调后获取的数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将新的比特对数似然比作为数据点的比特对数似然比。
可选的，处理器91执行该指令用于基于均值算法，根据比特对数似然比计算数据点承载的数据的均值，具体可以包括：
根据均值计算公式计算数据点承载的数据的均值；
均值计算公式包括：
其中，E(xj)为数据点承载的数据的均值，xj为数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据比特对数似然比计算Pr(bj,i)。
可选的，处理器91执行该指令用于基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'，具体可以包括：
将均值作为第一协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'；
第一协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的第二干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi，M))T，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
可选的，处理器91执行该指令用于基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'，具体还可以包括：
根据均值以及方差计算公式计算数据点承载的数据的方差；
方差计算公式包括：
其中E(xj)为数据点承载的数据的均值，Var(xj)为数据点承载的数据的方差，xj为数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据比特对数似然比计算Pr(bj,i)；
将均值与方差作为第二协方差矩阵计算公式的输入计算预设区域内的第二干扰自相关矩阵Ruu'；
第二协方差矩阵计算公式包括：
其中，Ruu'为预设区域内的第二干扰自相关矩阵，yk和Hk分别为预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
可选的，在获取数据点的比特对数似然比之后，处理器91执行该指令还用于：
计算第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当比值大于第一预设阈值，则需要挑选预设区域内满足预设条件的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当比值小于或等于第一预设阈值，则任意选取预设区域内的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点。
可选的，处理器91执行该指令用于挑选预设区域内满足预设条件的数据点作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点，具体可以包括：
将获取的预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，每个数据点承载的数据的方差包括预设层数数据的方差；
若数据点任一层数据的方差大于第二预设阈值，则删除数据点；
若数据点每层数据的方差小于或等于第二预设阈值，则确定数据点满足预设条件作为计算第二干扰自相关矩阵Ruu'步骤中使用的数据点。
可选的，在获取数据点的比特对数似然比之前，处理器91执行该指令还用于：
将用户设备的秩与接收天线数进行比较；
若秩小于接收天线数，则执行获取数据点的比特对数似然比；
若秩大于或等于接收天线数，则利用Ruu对接收信号y进行干扰抑制合并。
本发明实施例提供的抑制同频干扰的装置，首先利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵然后对和进行多入多出MIMO检测解调，获取数据点的比特对数似然比；基于均值算法，根据比特对数似然比计算数据点承载的数据的均值；基于协方差矩阵算法，根据均值，计算预设区域内的第二干扰自相关矩阵Ruu'；最后利用Ruu'对接收信号y进行干扰抑制合并。这样，能够提高计算干扰自相关矩阵的准确度，并且利用该矩阵可以有效抑制干扰。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (20)

1.一种抑制同频干扰的方法，其特征在于，包括：
利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对所述预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵其中，所述预设区域包括预设个数的资源块，所述资源块中包含有所述数据点和所述导频点；
对所述和进行多入多出MIMO检测解调，获取所述数据点的比特对数似然比；
基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值；
基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'；
利用所述Ruu'对所述接收信号y进行干扰抑制合并。
2.根据权利要求1所述的方法，其特征在于，所述获取所述数据点的比特对数似然比还包括：
对经过所述检测解调后获取的所述数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将所述新的比特对数似然比作为所述数据点的比特对数似然比。
3.根据权利要求1或2所述的方法，其特征在于，所述基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值包括：
根据均值计算公式计算所述数据点承载的数据的均值；
所述均值计算公式包括：
<mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mo>{</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>}</mo> </mrow> </munder> <msub> <mi>x</mi> <mi>j</mi> </msub> <munder> <mo>&amp;Pi;</mo> <mi>i</mi> </munder> <mi>Pr</mi> <mrow> <mo>(</mo> <msub> <mi>b</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
其中，E(xj)为所述数据点承载的数据的均值，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)。
4.根据权利要求1或2所述的方法，其特征在于，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'包括：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mo>&amp;lsqb;</mo> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>&amp;rsqb;</mo> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
5.根据权利要求3所述的方法，其特征在于，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'包括：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mo>&amp;lsqb;</mo> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>&amp;rsqb;</mo> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
6.根据权利要求1或2所述的方法，其特征在于，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'还包括：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mover> <mrow> <msup> <mi>HVH</mi> <mi>H</mi> </msup> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
7.根据权利要求3所述的方法，其特征在于，所述基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'还包括：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mover> <mrow> <msup> <mi>HVH</mi> <mi>H</mi> </msup> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
8.根据权利要求1所述的方法，其特征在于，在所述获取所述数据点的比特对数似然比之后还包括：
计算所述第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当所述比值大于第一预设阈值，则需要挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当所述比值小于或等于所述第一预设阈值，则任意选取所述预设区域内的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
9.根据权利要求8所述的方法，其特征在于，所述挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点包括：
将获取的所述预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，所述每个数据点承载的数据的方差包括预设层数数据的方差；
若所述数据点任一层数据的方差大于所述第二预设阈值，则删除所述数据点；
若所述数据点每层数据的方差小于或等于所述第二预设阈值，则确定所述数据点满足预设条件作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
10.根据权利要求1所述的方法，其特征在于，在所述获取所述数据点的比特对数似然比之前还包括：
将用户设备的秩与接收天线数进行比较；
若所述秩小于所述接收天线数，则执行所述获取所述数据点的比特对数似然比；
若所述秩大于或等于所述接收天线数，则利用所述Ruu对所述接收信号y进行干扰抑制合并。
11.一种抑制同频干扰的装置，其特征在于，包括：
第一抑制单元，用于利用时频域上预设区域内的导频点的第一干扰自相关矩阵Ruu对所述预设区域内的数据点的接收信号y和等效信道矩阵H进行干扰抑制合并，得到干扰抑制合并后的接收信号和等效频域信道矩阵其中，所述预设区域包括预设个数的资源块，所述资源块中包含有所述数据点和所述导频点；
第一获取单元，用于对所述和进行多入多出MIMO检测解调，获取所述数据点的比特对数似然比；
第二获取单元，用于基于均值算法，根据所述比特对数似然比计算所述数据点承载的数据的均值；
第三获取单元，用于基于协方差矩阵算法，根据所述均值，计算所述预设区域内的第二干扰自相关矩阵Ruu'；
第二抑制单元，用于利用所述Ruu'对所述接收信号y进行干扰抑制合并。
12.根据权利要求11所述的装置，其特征在于，所述第一获取单元还用于：
对经过所述检测解调后获取的所述数据点的比特对数似然比进行译码处理、软值速率匹配处理、编码块级联处理以及软值加扰处理后获取新的比特对数似然比；
将所述新的比特对数似然比作为所述数据点的比特对数似然比。
13.根据权利要求11或12所述的装置，其特征在于，所述第二获取单元具体用于：
根据均值计算公式计算所述数据点承载的数据的均值；
所述均值计算公式包括：
<mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mo>{</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>}</mo> </mrow> </munder> <msub> <mi>x</mi> <mi>j</mi> </msub> <munder> <mo>&amp;Pi;</mo> <mi>r</mi> </munder> <mi>Pr</mi> <mrow> <mo>(</mo> <msub> <mi>b</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
其中，E(xj)为所述数据点承载的数据的均值，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)。
14.根据权利要求11或12所述的装置，其特征在于，所述第三获取单元具体用于：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mo>&amp;lsqb;</mo> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>&amp;rsqb;</mo> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
15.根据权利要求13所述的装置，其特征在于，所述第三获取单元具体用于：
将所述均值作为第一协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第一协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mo>&amp;lsqb;</mo> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> 4
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值为：E(xi)＝(E(xi，1),…,E(xi,M))T，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，上标H表示对矩阵共轭转置。
16.根据权利要求11或12所述的装置，其特征在于，所述第三获取单元还具体用于：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mover> <mrow> <msup> <mi>HVH</mi> <mi>H</mi> </msup> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
17.根据权利要求13所述的装置，其特征在于，所述第三获取单元还具体用于：
根据所述均值以及方差计算公式计算所述数据点承载的数据的方差；
所述方差计算公式包括：
其中E(xj)为所述数据点承载的数据的均值，Var(xj)为所述数据点承载的数据的方差，xj为所述数据点的发射符号，Pr(bj,i)为对应xj的i个比特bj,i的概率，根据所述比特对数似然比计算Pr(bj,i)；
将所述均值与方差作为第二协方差矩阵计算公式的输入计算所述预设区域内的第二干扰自相关矩阵Ruu'；
所述第二协方差矩阵计算公式包括：
<mrow> <msup> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mi>u</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mi>i</mi> </msub> <mi>E</mi> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mover> <mrow> <msup> <mi>HVH</mi> <mi>H</mi> </msup> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow>
其中，Ruu'为所述预设区域内的第二干扰自相关矩阵，yk和Hk分别为所述预设区域内对应第k个导频点的接收信号和等效信道矩阵，pk为第k个导频符号，yi和Hi分别为所述预设区域内对应第i个数据点的接收信号和等效信道矩阵，xi为第i个数据点的发射符号，当xi为M维向量即xi为M层数据时，其均值和方差分别为：E(xi)＝(E(xi，1),…,E(xi,M))T和Vi＝diag{Var(xi,1),…,Var(xi,M)}，L为参与求和的所述预设区域内的导频点个数k与数据点个数i的和，用于补偿误差，上标H表示对矩阵共轭转置。
18.根据权利要求11所述的装置，其特征在于，所述装置还包括：
第一判断单元，用于计算所述第一干扰自相关矩阵Ruu中对角线元素模的平均值与非对角线元素模的平均值的比值；
当所述比值大于第一预设阈值，则需要挑选所述预设区域内满足预设条件的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点；
当所述比值小于或等于所述第一预设阈值，则任意选取所述预设区域内的数据点作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
19.根据权利要求18所述的装置，其特征在于，所述装置还包括：
挑选单元，用于挑选所述预设区域内满足预设条件的数据点作为所述第三获取单元使用的数据点，所述挑选所述预设区域内满足预设条件的数据点作为所述第三获取单元使用的数据点包括：
将获取的所述预设区域内的每个数据点承载的数据的方差与第二预设阈值比较，其中，所述每个数据点承载的数据的方差包括预设层数数据的方差；
若所述数据点任一层数据的方差大于所述第二预设阈值，则删除所述数据点；
若所述数据点每层数据的方差小于或等于所述第二预设阈值，则确定所述数据点满足预设条件作为计算所述第二干扰自相关矩阵Ruu'步骤中使用的数据点。
20.根据权利要求11所述的装置，其特征在于，所述装置还包括：
第二判断单元，用于在获取所述数据点的比特对数似然比之前，将用户设备的秩与接收天线数进行比较；
若所述秩小于所述接收天线数，则执行所述获取所述数据点的比特对数似然比；
若所述秩大于或等于所述接收天线数，则利用所述Ruu对所述接收信号y进行干扰抑制合并。
CN201410245825.XA 2014-06-04 2014-06-04 抑制同频干扰的方法和装置 Active CN103986562B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201410245825.XA CN103986562B (zh) 2014-06-04 2014-06-04 抑制同频干扰的方法和装置
PCT/CN2015/080498 WO2015184967A1 (zh) 2014-06-04 2015-06-01 抑制同频干扰的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410245825.XA CN103986562B (zh) 2014-06-04 2014-06-04 抑制同频干扰的方法和装置

Publications (2)

Publication Number Publication Date
CN103986562A CN103986562A (zh) 2014-08-13
CN103986562B true CN103986562B (zh) 2017-09-12

Family

ID=51278390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410245825.XA Active CN103986562B (zh) 2014-06-04 2014-06-04 抑制同频干扰的方法和装置

Country Status (2)

Country Link
CN (1) CN103986562B (zh)
WO (1) WO2015184967A1 (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103986562B (zh) * 2014-06-04 2017-09-12 华为技术有限公司 抑制同频干扰的方法和装置
CN104243377B (zh) * 2014-09-01 2017-08-29 华为技术有限公司 一种干扰抑制方法和装置
CN104683282B (zh) 2015-02-16 2018-03-30 深圳市中兴微电子技术有限公司 一种支持发射分集的干扰抑制合并方法和装置
CN107204817B (zh) * 2016-03-16 2021-08-17 大唐移动通信设备有限公司 一种干扰检测方法及终端
CN109150364B (zh) * 2017-06-28 2020-01-31 展讯通信（上海）有限公司 移动终端及其干扰噪声估计方法、计算机可读存储介质
CN111385229B (zh) * 2018-12-29 2022-10-11 中兴通讯股份有限公司 多普勒频移的确定方法及装置
CN116260547B (zh) * 2023-05-11 2023-09-22 武汉能钠智能装备技术股份有限公司四川省成都市分公司 一种抑制同频干扰的系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102647387A (zh) * 2012-04-09 2012-08-22 华为技术有限公司 同频干扰的消除方法及装置
CN103731384A (zh) * 2013-12-19 2014-04-16 华为技术有限公司 一种抑制干扰的方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5180224B2 (ja) * 2006-11-06 2013-04-10 クゥアルコム・インコーポレイテッド オンタイムシンボル干渉除去を備えるｍｉｍｏ検波
CN103986562B (zh) * 2014-06-04 2017-09-12 华为技术有限公司 抑制同频干扰的方法和装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102647387A (zh) * 2012-04-09 2012-08-22 华为技术有限公司 同频干扰的消除方法及装置
CN103731384A (zh) * 2013-12-19 2014-04-16 华为技术有限公司 一种抑制干扰的方法和装置

Also Published As

Publication number Publication date
CN103986562A (zh) 2014-08-13
WO2015184967A1 (zh) 2015-12-10

Similar Documents

Publication Publication Date Title
CN103986562B (zh) 抑制同频干扰的方法和装置
CN101039137B (zh) Mimo-ofdm系统基于码本搜索减少预编码反馈比特数的方法及装置
CN104618061B (zh) 一种大规模多天线系统中多用户信号的检测方法
CN102201847B (zh) 接收装置和接收方法
CN102571659B (zh) 一种干扰噪声估计和干扰抑制方法及相应系统
CN100395975C (zh) 多载波扩频信号的接收
CN102480444B (zh) 一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统
CN104243377B (zh) 一种干扰抑制方法和装置
Aref et al. Deep learning-aided successive interference cancellation for MIMO-NOMA
US20140334561A1 (en) Method and System for Symbol Detection Using Matrix Decomposition
CN109981151A (zh) 大规模mimo系统中改进的高斯树近似消息传递检测算法
CN106330284A (zh) 一种低复杂度大规模mimo信道估计方法
CN101312359B (zh) 多小区联合信道估计、多小区联合检测装置及方法
Yang et al. Deep learning based OFDM channel estimation using frequency-time division and attention mechanism
CN102546483B (zh) 一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统
CN103384180B (zh) 用于多输入多输出系统的空时译码方法和装置
CN102480453B (zh) 一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统
EP2804332A1 (en) Method and system for symbol detection using matrix decomposition
Shieh et al. Low-complexity soft-output sphere decoding with modified repeated tree search strategy
CN101997657A (zh) Mimo系统中宽度优先球形译码检测方法
CN105471778B (zh) 一种信号检测方法及装置
CN105337685B (zh) 一种图拨检测方法及其装置
CN103138820A (zh) Vamos模式下的信号检测方法及装置
CN102739282B (zh) Td-scdma系统的多天线联合检测系统及其方法
CN102480440B (zh) 一种宽带同频干扰环境下的信道估计方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210422

Address after: Unit 3401, unit a, building 6, Shenye Zhongcheng, No. 8089, Hongli West Road, Donghai community, Xiangmihu street, Futian District, Shenzhen, Guangdong 518040

Patentee after: Honor Device Co.,Ltd.

Address before: 518129 Bantian HUAWEI headquarters office building, Longgang District, Guangdong, Shenzhen

Patentee before: HUAWEI TECHNOLOGIES Co.,Ltd.

TR01 Transfer of patent right